The Iron Maiden. Cytosolic Aconitase/IRP1 Conformational Transition in the Regulation of Ferritin Translation and Iron Hemostasis

Maintaining iron homeostasis is fundamental for almost all living beings, and its deregulation correlates with severe and debilitating pathologies. The process is made more complicated by the omnipresence of iron and by its role as a fundamental component of a number of crucial metallo proteins. The response to modifications in the amount of the free-iron pool is performed via the inhibition of ferritin translation by sequestering consensus messenger RNA (mRNA) sequences. In turn, this is regulated by the iron-sensitive conformational equilibrium between cytosolic aconitase and IRP1, mediated by the presence of an iron–sulfur cluster. In this contribution, we analyze by full-atom molecular dynamics simulation, the factors leading to both the interaction with mRNA and the conformational transition. Furthermore, the role of the iron–sulfur cluster in driving the conformational transition is assessed by obtaining the related free energy profile via enhanced sampling molecular dynamics simulations.


Introduction
Iron is a fundamental element, which is involved in key biological processes [1][2][3][4], such as oxygen transport in hemoglobin, electron transfer, muscle contraction, and serves as a co-factor for many enzymes. As a matter of fact, a number of iron-containing proteins were identified, and they are involved in extremely diverse functions [5]. For instance, Fe-S clusters have been assigned both a structural and redox role in DNA repair proteins and other enzymes [6][7][8]. In addition to protein-embedded iron, a limited quantity of free iron is also available in cells, in what is called the free-iron pool (FIP) [9][10][11]. FIP is constituted by iron in both II and III oxidation states and the metal is usually complexed with different, and not very well characterized, ligands and chelatants [12]. Despite its limited characterization, FIP is essential for cells viability, and its dysregulation may be correlated with cellular instability, and at a systemic level, with the development of debilitating pathologies, including cancer [13][14][15][16] and mitochondria failure [17,18]. On the other hand, FIP iron accumulation, and especially the accumulation of Fe 2+ , also induces considerable oxidative stress, most notably due to the production of peroxyl radicals via the Fenton reaction [19]. As a logical consequence, cells have developed complex regulatory machinery to stabilize and regulate the amount of FIP [20,21].
It is recognized that the regulation of iron homeostasis involves tuning the expression of several genes involved in storage, outtake and intake of iron such as ferritin H and L (involved in iron storage), ferroportin (involved in iron outtake) and transferrin receptor It is recognized that the regulation of iron homeostasis involves tuning the expression of several genes involved in storage, outtake and intake of iron such as ferritin H and L (involved in iron storage), ferroportin (involved in iron outtake) and transferrin receptor (TFR1 involved in iron uptake) [22,23]. These regulations have to be monitored rapidly to react to all changes in iron availability. Thus, cells possess a post-transcriptional regulatory mechanism that involves the binding of messenger RNA (mRNA) consensus sites present in untranslated sequences (UTR) of the iron regulatory genes by specific proteins [24][25][26][27][28].
Indeed, in presence of a high iron level, cytosolic aconitase (cAco) embeds a Fe-S cluster which confers an enzymatic activity catalyzing the isomerization of citrate in ciscitrate via cis-aconitate. Conversely, when the FIP level decreases, cAco loses the Fe-S cluster triggering a conformational transition that is accompanied by the loss of the catalytic activity [29,30]. Furthermore, in its iron-free conformation, the protein, known as iron regulating protein 1 (IRP1), is able to bind a hairpin mRNA sequence, the iron regulating element (IRE), which is highly conserved and present in the ferritin, ferroportin and transferrin receptor transcripts of vertebrates, while it is absent in non-vertebrate despite the conservation of the cAco/IRP1 equilibrium [31][32][33]. The binding of IRP1 on IRE, located in 5′UTR, induces a blockage of ferritin and ferroportin translation, leading to their downexpression and ultimately to the increase in FIP level (Scheme 1) [34,35]. On the contrary, the binding of IRP1 on the IRE present in the 3′UTR of the transferrin receptor leads to mRNA stabilization and thus to an elevation of its translation allowing an increase in iron uptake. Once the cellular iron concentration rises again sufficiently, the equilibrium between cAco and IRP1 is altered releasing IRE and reinstating the translation of these genes. From a structural point of view, different crystal structures of cAco and IRP1, also complexed with the IRE hairpin sequence, were resolved [7,36]. Evidencing that the conformational translation responsible for homeostasis regulation can be summarized as the passage from a closed (cAco) to an open (IRP1) structure which exposes the RNA binding motifs. However, some important questions still hold concerning the molecular bases From a structural point of view, different crystal structures of cAco and IRP1, also complexed with the IRE hairpin sequence, were resolved [7,36]. Evidencing that the conformational translation responsible for homeostasis regulation can be summarized as the passage from a closed (cAco) to an open (IRP1) structure which exposes the RNA binding motifs. However, some important questions still hold concerning the molecular bases driving the mechanisms. In particular, the role of the Fe-S cluster in stabilizing the cAco form has not been clearly evaluated, as well as the reasons leading to the selective interaction with the IRE consensus sequence. Furthermore, the flexibility and intrinsic dynamic of IRP1 bound and unbound to RNA is also lacking and could provide important insights to clarify its mechanism of action.
In this contribution, we resort to all atom classical molecular dynamic (MD) simulations to explore the behavior of cAco, IRP1, and IRE and their mutual interactions using the consensus IRE of the ferritin transcript as a model. Furthermore, we also provide enhanced sampling MD simulation to assess the free energy profile for the cAco → IRP1 transition in the presence and in the absence of a Fe-S cluster, hence assessing the thermodynamics factors at the base of the first step of iron hemostasis.
The deregulation of the cellular iron level is associated with a series of debilitating pathologies, including cancer, yet its current comprehension, from a structural point of view, is still far from being ideal. Unraveling the key biophysical effects behind the complex structural transitions taking place may lead, in the future, to the development of novel therapeutic strategies, also involving therapeutic mRNA protocols.

Materials and Methods
The initial structure of the IRP1/IRE complex was obtained from pdb code 3snp [7], and the initial conditions for the MD simulations of the isolated IRP1 and IRE moieties were obtained from the same crystal structure. The structure of the human cAco enzyme was instead retrieved from the pdb code 2b3y [36]. Even if the two structures refer to different organisms, the sequence homology exceeds 95%, hence assuring meaningful comparison.
For all the MD simulations, the amber2014 force field was considered [37], while in the case of nucleic acid, the bsc1 corrections [38,39] to enhance the description of the backbone parameters were also taken into account. The Fe-S (4Fe-4S) cluster was modeled using a specifically tailored set of parameters designed by Carvalho and Swart [40] and already used by some of us in the modeling of endonucleases [41]. All the systems were then solvated in a cubic TIP3P [42] water box with a 9 Å buffer and K + cations added to ensure electroneutrality. MD simulations were performed using the NAMD code [43,44], and a constant temperature (300 K) and pressure (1 atm) were enforced using the Langevin thermostat [45] and barometer [46], hence yielding an NPT simulation. After 5000 minimization steps, the system was equilibrated for 9 ns during which constraints on the protein backbone and on the nucleic acid heavy atoms were progressively released, and finally, the production runs were performed. The equilibrium MD simulation reached 200 ns for the IRP1/IRE, IRP1, and IRE. In the case of cAco including the Fe-S cluster, we also performed 200 ns, while the MD simulation of Fe-S-missing cAco was propagated for 500 ns. A time step of 4 fs was consistently used thanks to the hydrogen mass reparation (HMR) strategy [47] in combination with the Rattle and Shake procedure [48]. Particle mesh Ewald (PME) [49] was used to treat long-range interactions with a cut-off of 9 Å.
For the enhanced sampling procedure, we chose to sample a collective variable based on the RMSD difference (∆RMSD) between two representative states, namely the open and closed conformations. While the closed conformation can be easily obtained from the human cAco equilibrium simulation, the open conformation was built starting from the equilibrium MD of the chicken IRP1 after manually mutating the residue to obtain the same sequence of the human protein. Note that the stability of the mutated system was checked and validated by an equilibrium MD. The enhanced sampling MD was performed on two situations, i.e., including or excluding the Fe-S cluster. Since the latter is not covalently bound to the cAco enzyme and to avoid its release in the medium, a constraint to the distance with the Ser 435 and Cys 436 residues was imposed. The potential of mean force (PMF) was obtained using a combination of extended adaptative biasing force (eABF) [50] and metadynamics [51] in what is known as meta-eABF [52,53]. Most notably, this choice allowed us to treat the system of interest without its separation in windows. Enhanced sampling MD simulations were performed using NAMD and its Colvar interface and the total sampling for the two systems reached up to 300 ns. Walls were placed at DRMSD of −9.7 and 9.7 Å, respectively, and Gaussians with a width of 0.1 were deposited with a height of 5 Å and a weight of 0.1 kcal/mol. All the simulations were visualized and analyzed using VMD [54].

Structure of IRE mRNA
As stated in the introduction the post-transcriptional regulation of ferritin expression is regulated by a well-conserved mRNA motif, IRE, which assumes a hairpin sequence [55,56]. Notably, while some sequence variations in IRE exist, both between species and in the same organism, the hairpin motif appears as fundamental to assure recognition by IRP1, with a particular role played by the loop region, which is particularly well conserved [35]. In addition, some hanging uncoupled nucleobases are also present and may lead to supplementary interactions with the regulatory protein. Representative snapshots extracted from the MD simulation of a solvated IRE strand are reported in Figure 1 together with the time series of the root mean square deviation (RMSD). Two contrasting pieces of evidence may be underlined from its inspection: while globally the strand appears quite flexible, leading to a relatively high value of the RMSD (~11 Å), the loop region is well conserved, as well as the global hairpin arrangement. Interestingly, the most important movement observed can be described as a global bending of the RNA strand. This in turn leads to a coexistence of straight and bent sequences, in which the loop area is approaching the double-strand arrangement and the terminal bases. Notably, two overhanging, solvent-exposed bases are also observed at the edges of the loop and in the double-strand region. If these extrahelical bases experience, as expected, enhanced flexibility, they do not induce any particular structural deformation. The global flexibility of the RNA hairpin is also confirmed by the fact that the RMSD spans a relatively large interval and experiences an oscillatory behavior. Coherently, the bending movement can also be evidenced by monitoring the distance between the hanging nucleobases U 6 (in the loop region) and U 19 (in the double-helical region), which is reported in Supplementary Information and which shows a behavior strongly correlated with the RMSD confirming the large-scale breathing mode.
windows. Enhanced sampling MD simulations were performed using NAMD and its Colvar interface and the total sampling for the two systems reached up to 300 ns. Walls were placed at DRMSD of −9.7 and 9.7 Å, respectively, and Gaussians with a width of 0.1 were deposited with a height of 5 Å and a weight of 0.1 kcal/mol.
All the simulations were visualized and analyzed using VMD [55].

Structure of IRE mRNA
As stated in the introduction the post-transcriptional regulation of ferritin expression is regulated by a well-conserved mRNA motif, IRE, which assumes a hairpin sequence [56,57]. Notably, while some sequence variations in IRE exist, both between species and in the same organism, the hairpin motif appears as fundamental to assure recognition by IRP1, with a particular role played by the loop region, which is particularly well conserved [35]. In addition, some hanging uncoupled nucleobases are also present and may lead to supplementary interactions with the regulatory protein. Representative snapshots extracted from the MD simulation of a solvated IRE strand are reported in Figure 1 together with the time series of the root mean square deviation (RMSD). Two contrasting pieces of evidence may be underlined from its inspection: while globally the strand appears quite flexible, leading to a relatively high value of the RMSD (~11 Å), the loop region is well conserved, as well as the global hairpin arrangement. Interestingly, the most important movement observed can be described as a global bending of the RNA strand. This in turn leads to a coexistence of straight and bent sequences, in which the loop area is approaching the double-strand arrangement and the terminal bases. Notably, two overhanging, solvent-exposed bases are also observed at the edges of the loop and in the double-strand region. If these extrahelical bases experience, as expected, enhanced flexibility, they do not induce any particular structural deformation. The global flexibility of the RNA hairpin is also confirmed by the fact that the RMSD spans a relatively large interval and experiences an oscillatory behavior. Coherently, the bending movement can also be evidenced by monitoring the distance between the hanging nucleobases U 6 (in the loop region) and U 19 (in the double-helical region), which is reported in Supplementary Information and which shows a behavior strongly correlated with the RMSD confirming the large-scale breathing mode.

Structure of the IRP/IRE Complex
As stated in the introduction the regulation of ferritin translation is modulated by the complexation of the IRE unit performed by the IRP1 protein, i.e., the same protein sequence forming cAco but lacking the Fe-S cluster. The results of the MD simulation for the IRP1/IRE complex are reported in Figure 2, which shows that the interaction between the RNA sequence and the protein leads to the formation of a persistent and rather rigid aggregate. Interestingly, while the central double-strand helical region of the nucleic acid stays rather solvent-exposed, the main contact points involve the loop and the edge bases, as well as the extrahelical extruded nucleotides, mainly uracil. A striking reduction in the flexibility of the mRNA is observed upon complexation, since its RMSD now barely reaches 3 Å. the RNA sequence and the protein leads to the formation of a persistent and rather rigid aggregate. Interestingly, while the central double-strand helical region of the nucleic acid stays rather solvent-exposed, the main contact points involve the loop and the edge bases, as well as the extrahelical extruded nucleotides, mainly uracil. A striking reduction in the flexibility of the mRNA is observed upon complexation, since its RMSD now barely reaches 3 Å. The same rigidity is also observed for the protein moiety, whose RMSD is also plateauing at around 3 Å with rare peaks at 4 Å. As mentioned, two main contact points can be observed and are highlighted in Figure 2A,B. They involve either the loop region (1) or the edge bases of the mRNA sequence (2) and appear as extremely persistent all along the trajectory. Anchoring the RNA at the two extremities may also justify the increased rigidity of the strand. The involvement of the loop region in the direct contact with IRP1 protein is also coherent with the experimental observations pointing to the fact that this motif is, indeed, well conserved in IRE regulatory area, and is necessary to regulate ferritin translation. In both regions, we may identify a complex network of amino acids interacting both with the RNA backbone and with specific and exposed nucleobases motifs as reported in Figure 3. Interestingly enough, while the prominent amino acids in the contact area are polar or basic, such as Asn, Asp, or Thr, no salt bridges involving the RNA backbone can be evidenced. This observation is in striking contrast with some of the most common protein/nucleic acid interaction patterns that involve the rather unspecific electrostatic binding of the negatively charged backbone of the DNA/RNA fragment. The interaction involving hydrogen bonds with specific nucleobases is also especially important in the loop region, most probably adding up to the global selectivity of the interaction. The same rigidity is also observed for the protein moiety, whose RMSD is also plateauing at around 3 Å with rare peaks at 4 Å. As mentioned, two main contact points can be observed and are highlighted in Figure 2A,B. They involve either the loop region (1) or the edge bases of the mRNA sequence (2) and appear as extremely persistent all along the trajectory. Anchoring the RNA at the two extremities may also justify the increased rigidity of the strand. The involvement of the loop region in the direct contact with IRP1 protein is also coherent with the experimental observations pointing to the fact that this motif is, indeed, well conserved in IRE regulatory area, and is necessary to regulate ferritin translation. In both regions, we may identify a complex network of amino acids interacting both with the RNA backbone and with specific and exposed nucleobases motifs as reported in Figure 3. Interestingly enough, while the prominent amino acids in the contact area are polar or basic, such as Asn, Asp, or Thr, no salt bridges involving the RNA backbone can be evidenced. This observation is in striking contrast with some of the most common protein/nucleic acid interaction patterns that involve the rather unspecific electrostatic binding of the negatively charged backbone of the DNA/RNA fragment. The interaction involving hydrogen bonds with specific nucleobases is also especially important in the loop region, most probably adding up to the global selectivity of the interaction.
Interestingly, and as reported in Supplementary Information, we may also identify some interactions which, while not present in the crystal structure, develop in the course of the MD simulation and are quite persistent, see for instance Arg 353 and Arg 148 developing salt bridges with the phosphate of guanine and uracil, respectively. This reorganization however is local and is not accompanied by any major structural deformation of either the protein or the nucleic acid; on the contrary, it confirms the high stability and affinity of the IRP1/IRE complex. In Figure 4 we also report the time series of some distances characterizing chosen paradigmatic interactions involving either the loop ( Figure 4A   Interestingly, and as reported in Supplementary Information, we may also identify some interactions which, while not present in the crystal structure, develop in the course   Interestingly, and as reported in Supplementary Information, we may also identify some interactions which, while not present in the crystal structure, develop in the course  Indeed, we may also underline that the negative RNA phosphates in the backbone are mainly stabilized via the formation of hydrogen bonds, as shown in Figure 4A involving the OH group of a Thr, or the amine group of Asn in Figure 4B. Note, however, that while the interaction in Figure 4A was present already in the crystal structure and leads to an extremely stable time series, the second one appears looser, showing some evidence of local and momentarily loss of the interaction. This was particularly true before equilibration, in the crystal structure the characteristic distance exceeds 6 Å and it populates values compatible with the formation of stable hydrogen bonds only after equilibration.
Conversely, in Figure 4C,D we evidence interaction taking place directly with the nucleobase elements. While Figure 4C shows the interaction between the lateral chain of Asn and the electron-rich oxygen atom of uracil, Figure 4B deals with the interaction between the carbonyl group of cytosine and the peptidic NH domain of an Asp residue. This latter evidence is important since it also shows that the protein backbone is directly participating in the stabilization, which also constitutes a most paradigmatic example of the prominence of Hydrogen bond patterns over salt bridges. Furthermore, while the interaction represented in Figure 4D is extremely stable, as evidenced by the time series, the more classic one in Figure 4C happens to be quite dynamic and flexible, and hence less important to assure the global stability of the complex.
As final notes, one should also point out that the hanging nucleobase of IRE participates in the recognition, and that because the contact points are located at the extremity of the RNA strand, IRP1 should assume a rather extended and open conformation to accommodate and efficiently sequestrate the mRNA strand by assuming a pincer conformation.

Structural Features of Unbound IRP1 and cAco
The analysis of the equilibrium MD simulation as seen in Figure 5 confirms that the protein remains globally open and maintains the pincer-like shape. However, an increase in flexibility is evidenced upon the removal of IRE. Indeed, in the apo form IRP1 exhibits a sort of oscillatory or breathing movement characterized by the partial closing and opening of the pincer, as can be observed in Figure 5A,B. If globally, the distance between the extremities of the pincers spans an interval of more than 25 Å, it is also evident that values below 50 Å are not observed in the course of the MD simulation. Note, as a comparison that a distance between the pincers edges of about 60 Å is observed both in the crystal structure and in the IRP1/IRE simulation. Hence, while the absence of IRE induces a global flexibilization of IRP1 significant energetic barriers prevent the protein from assuming more closed conformations, which would be incompatible with mRNA recognition and sequestration. In the case of cAco, which has the same sequence of IRP1 with the supplementary inclusion of the Fe-S cluster, the MD simulation shows different behavior. Indeed, cAco presents a very rigid and compact structure in which the pincers are totally closed, burying the IRE recognition area in the solvent inaccessible regions (Figure 6). Note also that In the case of cAco, which has the same sequence of IRP1 with the supplementary inclusion of the Fe-S cluster, the MD simulation shows different behavior. Indeed, cAco presents a very rigid and compact structure in which the pincers are totally closed, burying the IRE recognition area in the solvent inaccessible regions ( Figure 6). Note also that the Fe-S cluster is buried deep in cAco and is located at the basis of the pincers, where it interacts particularly with the residues Hys 206, Asn 297, Ser 435, Cys 436, and Asn 534, exerting a force sufficient to lock the pincers in the closed conformation. The crucial role of the Fe-S cluster can also be evidenced by the results of the MD simulation in which the former was removed from the cAco closed conformation. In this case ( Figure 6B) even though we do not recover the equilibrium open conformation of IRP1 we can note more pronounced flexibility of the protein and more importantly the appearance of a first partial opening of the pincers region.

Thermodynamics of the IRP/Aco Transition
To better assess the thermodynamic and the structural mechanisms involved in the cAco/IRP1 conformational transition we resorted to enhanced sampling techniques providing the free energy profile for the interconversion between closed and open conformations in presence or in absence of the Fe-S cluster ligand. To this aim and due to a rather complex character of the transition, we used the RMSD difference between the reference conformations as a collective variable. Indeed, even if RMSD can be degenerate and hence could lead to the exploration of unphysical states, it allows capturing global deformations like those that are relevant here.
As seen from Figure 7A in the case of the apo system, i.e., the Fe-S-missing equivalent of IRP1, we can identify a rather large and flat-free energy area corresponding to the absolute minimum extending roughly to the −2 to +2 Å ΔRMSD region. As also seen from the representative snapshots reported in the inlay this area corresponds to a closed conformation of the protein, which is equivalent to an IRP1 functional conformation. This

Thermodynamics of the IRP/Aco Transition
To better assess the thermodynamic and the structural mechanisms involved in the cAco/IRP1 conformational transition we resorted to enhanced sampling techniques providing the free energy profile for the interconversion between closed and open conformations in presence or in absence of the Fe-S cluster ligand. To this aim and due to a rather complex character of the transition, we used the RMSD difference between the reference conformations as a collective variable. Indeed, even if RMSD can be degenerate and hence could lead to the exploration of unphysical states, it allows capturing global deformations like those that are relevant here.
As seen from Figure 7A in the case of the apo system, i.e., the Fe-S-missing equivalent of IRP1, we can identify a rather large and flat-free energy area corresponding to the absolute minimum extending roughly to the −2 to +2 Å ∆RMSD region. As also seen from the representative snapshots reported in the inlay this area corresponds to a closed conformation of the protein, which is equivalent to an IRP1 functional conformation. This result is also coherent with the observations of the equilibrium MD which point to an increased flexibility of the IRE-unbound IRP1 and a slight deviation from the reference crystal structure which is obtained from the IRP1/IRE complex. On the other hand, closed conformations can be observed at ∆RMSD of about −8 Å, which are fully consistent with the crystallographic structure of the functional cAco. However, the closed conformation for the apo system clearly does not correspond to a free energy minimum, since only a flat surface can be observed and is, hence, unstable. Furthermore, a free energy penalty of about 10 kcal/mol is associated with this conformation as compared to the open conformation. Note that some deep minima can be observed, whose representative structures are given in  The free energy profile is drastically changed when considering the inclusion of the Fe-S cluster ( Figure 7B). Indeed, while we may still recognize the flat area close to a ΔRMSD of 0.0 Å and indicative of open conformations, this area is no more associated with the global free energy minimum. The closed conformation (ΔRMSD ~−8.0 Å) is now highly stabilized and presents a rather steep free energy well leading to a stabilization of around 15 kcal/mol compared to the open conformation. Interestingly, only a rather small barrier of about 5 kcal/mol should be overcome for passing from the open to the closed state, while the opposite transition would require bypassing 20 kcal/mol. Once again, these results are coherent with the observation raised from the equilibrium MD concerning the flexibility and the breathing mode of the apo system. The large stabilization of the closed conformation by the Fe-S cluster is mainly due to specific interactions taking place with amino acids belonging to the different protein regions which help to maintain the compact structure.

Discussion and Conclusions
By using extended MD simulations, including enhanced sampling and free energy methods, we provided a unified vision of the molecular basis involved in the complex regulation of iron hemostasis and in the regulation of the translation of ferritin. Indeed, from one side, we examined the specific interactions exerted between IRP1 protein and the IRE region of ferritin mRNA. Differently from other nucleic acid interacting proteins [42,[58][59][60][61][62][63], IRP1 is able to develop strong and specific hydrogen bonds not only with the backbone but also with specific IRE nucleobases; in particular, the role of the mRNA hairpin region in favoring the interaction is clearly evidenced. Since the interaction with IRP1 involves different IRE structural domains, the protein should maintain an open conformation to allow for the efficient sequestration of the mRNA. Interestingly, we evidenced that the interaction network between IRE and IRP1 is quite complex and involves also The free energy profile is drastically changed when considering the inclusion of the Fe-S cluster ( Figure 7B). Indeed, while we may still recognize the flat area close to a ∆RMSD of 0.0 Å and indicative of open conformations, this area is no more associated with the global free energy minimum. The closed conformation (∆RMSD~−8.0 Å) is now highly stabilized and presents a rather steep free energy well leading to a stabilization of around 15 kcal/mol compared to the open conformation. Interestingly, only a rather small barrier of about 5 kcal/mol should be overcome for passing from the open to the closed state, while the opposite transition would require bypassing 20 kcal/mol. Once again, these results are coherent with the observation raised from the equilibrium MD concerning the flexibility and the breathing mode of the apo system. The large stabilization of the closed conformation by the Fe-S cluster is mainly due to specific interactions taking place with amino acids belonging to the different protein regions which help to maintain the compact structure.

Discussion and Conclusions
By using extended MD simulations, including enhanced sampling and free energy methods, we provided a unified vision of the molecular basis involved in the complex regulation of iron hemostasis and in the regulation of the translation of ferritin. Indeed, from one side, we examined the specific interactions exerted between IRP1 protein and the IRE region of ferritin mRNA. Differently from other nucleic acid interacting proteins [41,[57][58][59][60][61][62], IRP1 is able to develop strong and specific hydrogen bonds not only with the backbone but also with specific IRE nucleobases; in particular, the role of the mRNA hairpin region in favoring the interaction is clearly evidenced. Since the interaction with IRP1 involves different IRE structural domains, the protein should maintain an open conformation to allow for the efficient sequestration of the mRNA. Interestingly, we evidenced that the interaction network between IRE and IRP1 is quite complex and involves also transient interactions which are not always captured by the crystallographic structure.
The role of the Fe-S cluster in driving the conformational transitions between the open (IRP1) form of the protein and its closed state (equivalent to cAco) was rationalized. Indeed, the specific interactions developed by the Fe-S cluster, which is buried inside the protein, favors the closed conformation, leading to a free energy gain of about 15 kcal/mol, only relatively small barriers (about 5 kcal/mol) should be overcome for achieving the open/closed transition. On the contrary, the open conformation is more stable in the absence of the stabilization network induced by Fe-S leading to a 10 kcal/mol preference for the open conformation. Note also that in the IRE/unbound state, the open conformation is characterized by a large and rather flat free energy surface, which is compatible with the observed breathing oscillation of the pincers.
Our results provide for the first time a clear atomistic view of the complex mechanisms at the base of FIP regulation. They are also instrumental to identify crucial regions of the IRP1/cAco protein which could be potentially druggable in the perspective of controlling iron hemostasis. The involvement of mRNA sequence in the process and the precise determination of its structure and role can also be potentially important in the framework of therapeutic RNA strategies, which are gaining significant momentum in cancer therapy.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/biom11091329/s1, Figure S1: Time series of the RMSD for the isolated IRP1, Figure S2: Time series of the Solvent Accessible Solvent Area (SASA) for isolated IRP1 and the IRE/IRP1 complex, Figure S3: Time series of the Solvent Accessible Solvent Area (SASA) for cAco with and without the inclusion of the Fe-S cluster, Figure S4: Time series of the distances between the interfacial residues of IRP1 and IRE, Figure