Free energy simulations of a GTPase: GTP and GDP binding to archaeal initiation factor 2

Divalent-Metal-Ion Selectivity of the CRISPR-Cas System-Associated Cas1 Protein: Insights from Classical Molecular Dynamics Simulations and Electronic Structure Calculations

CRISPR-associated protein 1 (Cas1) is a universally conserved essential metalloenzyme of the clustered regularly interspaced short palindromic repeat (CRISPR) immune system of prokaryotes (bacteria, archaea) that can cut and integrate a part of viral DNA to its host genome with the help of other proteins. The integrated DNA acts as a memory of viral infection, which can be transcribed to RNA and stop future infection by recognition (based on the RNA/DNA complementarity principle) followed by protein-mediated degradation of the viral DNA. It has been proposed that the presence of a single manganese (Mn²⁺) ion in a conserved divalent-metal-ion binding pocket (key residues: E190, H254, D265, D268) of Cas1 is crucial for its function. Cas1-mediated DNA degradation was proposed to be hindered by metal substitution, metal chelation, or mutation of the binding pocket residues. Cas1 is active toward dsDNA degradation with both Mn²⁺ and Mg²⁺. X-ray structures of Cas1 revealed an intricate atomic interaction network of the divalent-metal-ion binding pocket and opened up the possibility of modeling related metal ions (viz., Mg²⁺, Ca²⁺) in the binding pocket of wild-type (WT) and mutated Cas1 proteins for computational analysis, which includes (1) quantitative estimation of the energetics of the divalent-metal-ion preference and (2) exploring the structural and dynamical aspects of the protein in response to divalent-metal-ion substitution or amino acid mutation. Using the X-ray structure of the Cas1 protein from Pseudomonas aeruginosa as a template (PDB 3GOD), we performed (∼2.23 μs) classical molecular dynamics (MD) simulations to compare structural and dynamical differences between Mg²⁺- and Ca²⁺-bound binding pockets of wild-type (WT) and mutant (E190A, H254A, D265A, D268A) Cas1. Furthermore, reduced binding pocket models were generated from X-ray and molecular dynamics (MD) trajectories, and the resulting structures were subjected to quantum chemical calculations. Results suggest that Cas1 prefers Mg²⁺ binding relative to Ca²⁺ and the preference is the strongest for WT and the weakest for the D268A mutant. Quantum chemical calculations indicate that Mn²⁺ is the most preferred relative to both Mg²⁺ and Ca²⁺ in the wild-type and mutant Cas1. Substitution of Mg²⁺ by Ca²⁺ does not alter the interaction network between Cas1 and the divalent metal ion but increases the wetness of the binding pocket by introducing a single water molecule in the first coordination shell of the latter. The strength of metal-ion preference (Mg²⁺ versus Ca²⁺) seems to be dependent on the solvent accessibility of the divalent-metal-ion binding pocket, strongest for wild-type Cas1 (in which the metal-ion binding pocket is dry, which includes two water molecules) and the weakest for the D268A mutant (in which the metal-ion binding pocket is wet, which includes four water molecules).

Molecular dynamics free energy simulations of ATP:Mg2+ and ADP:Mg2+ using the polarizable force field AMOEBA

ATPases and GTPases are two important classes of protein that play critical roles in energy transduction, cellular signaling, gene regulation and catalysis. These proteins use cofactors such as nucleoside di and tri-phosphates (NTP, NDP) and can detect the difference between NDP and NTP which then induce different protein conformations. Mechanisms that drive proteins into the NTP or NDP conformation may depend on factors such as ligand structure and how Mg²⁺ coordinates with the ligand, amino acids in the pocket and water molecules. Here, we have used the advanced electrostatic and polarizable force field AMOEBA and molecular dynamics free energy simulations (MDFE) to examine the various binding mechanisms of ATP:Mg²⁺ and ADP:Mg²⁺.We compared the ATP:Mg²⁺ binding with previous studies using non-polarizable force fields and experimental data on the binding affinity. It was found that the total free energy of binding for ATP:Mg²⁺ (-7.00 ± 2.13 kcal/mol) is in good agreement with experimental values (-8.6 ± .2 kcal/mol)¹. In addition, parameters for relevant protonation states of ATP, ADP, GTP and GDP have been derived. These parameters will allow for researchers to investigate biochemical phenomena involving NTP's and NDP's with greater accuracy than previous studies involving non-polarizable force fields.

Why double-stranded RNA with 5'-monophosphate is a poor binder to retinoic acid-inducible gene-I with respect to 5'-hydroxyl analogue?

AbbreviationsAActiveABNRAdopted basis Newton-RaphsonCARDsCaspase activation and recruitment domainsCTDC-terminal domaindsRNADouble-stranded RNAHEL1Helicase1HEL2Helicase2HEL2iHelicase2 insertionMDMolecular dynamicsNAInactivePPincerPDBprotein data bankPMEParticle mesh EwaldRIG-IRetinoic acid-inducible gene-IRMSDRoot mean square deviationRMSFRoot mean square fluctuationSDSteepest descentTIThermodynamic IntegrationCommunicated by Ramaswamy H. Sarma.

Mg2+ vs Ca2+ bound active site of group II intron- A MD study

Group II introns are enzymes which undergo self-splicing and remove itself from pre-messenger RNA. X-ray structures of group II intron of Oceanobacillus iheyensis at various stages of the self-splicing pathway (pre-hydrolytic, post-hydrolytic, and ligand-free state) revealed intricate atomic interaction network in the active site of the intron. It has been confirmed that a heteronuclear metal ion cluster consisting of four metal ions (K1, K2 sites with K⁺ and M1, M2 sites with Mg²⁺) are crucial for function. Substitution of Mg²⁺ by Ca²⁺ results in loss of enzymatic activity. The X-ray structures not only opens up the possibility of modelling Mg²⁺ and Ca²⁺ bound active site of group II intron and quantitatively estimate the energetics of Mg²⁺ vs Ca²⁺ preference but also explore the relative structural and dynamical differences in response to divalent metal ion substitution. Thus, using X-ray structures as a template we performed molecular dynamics simulations to compare structural and dynamical differences between Mg²⁺ and Ca²⁺ bound active site of group II intron at various stages of the splicing pathway (i.e., pre-hydrolytic, post-hydrolytic, and ligand-free state). Quantitative estimation of Mg²⁺ vs Ca²⁺ selectivity at the M1, M2 sites confirmed Mg²⁺ preference at intron active sites relative to Ca²⁺. Ca²⁺ is relatively more hydrated in the intron active site relative to Mg²⁺. The local environment (bound nucleophilic water, interaction with scissile phosphate) around M2 is strikingly different between Mg²⁺ and Ca²⁺ bound pre-hydrolytic state. In the post-hydrolytic state, the exon part of the hydrolysis product is involved in direct interaction with the M1, whereas the intron part is highly flexible in our MD trajectories. Solvent exposure of M1, M2 sites are least in the pre-hydrolytic state, highest in the ligand-free state, and intermediate in the post-hydrolytic state.

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

Classical molecular dynamics simulations reveal size-dependent trends of alkaline earth metal ions binding to DNA are due to ion size and hydration behavior.

Principles of tRNAAla Selection by Alanyl-tRNA Synthetase Based on the Critical G3·U70 Base Pair

Throughout evolution, the presence of a single G3·U70 mismatch in the acceptor stem of tRNA^Ala is the major determinant for aminoacylation with alanine by alanyl-tRNA synthetase (AlaRS). Recently reported crystal structures of the complexes AlaRS-tRNA^Ala/G3·U70 and AlaRS-tRNA^Ala/A3·U70 suggest two very different conformations, representing a reactive and a nonreactive state, respectively. On the basis of these structures, it has been proposed that the G3·U70 base pair guides the -CCA end of the tRNA acceptor stem into the active site of AlaRS, thereby enabling aminoacylation. The crystal structures open up the possibility of directly computing the energetics of tRNA specificity by AlaRS. We have carried out molecular dynamics free-energy simulations to quantitatively estimate tRNA discrimination by AlaRS, focusing on the mutations of the single critical base pair G3·U70 to uncover the energetics underlying the accuracy of tRNA selection. The calculations show that the reactive complex is highly selective in favor of the cognate tRNA^Ala/G3·U70 over its noncognate analogues (A3·U70/G3·C70/A3·C70). In contrast, the nonreactive complex is predicted to be unselective between tRNA^Ala/G3·U70 and tRNA^Ala/A3·U70. Utilizing our calculated relative binding free energies, we show how a simple three-step kinetic scheme for aminoacylation, involving both an initial nonspecific binding step and a subsequent transition to a selective reactive complex, accounts for the observed kinetics of the process.

Polarizable Force Fields for Biomolecular Simulations: Recent Advances and Applications

Realistic modeling of biomolecular systems requires an accurate treatment of electrostatics, including electronic polarization. Due to recent advances in physical models, simulation algorithms, and computing hardware, biomolecular simulations with advanced force fields at biologically relevant timescales are becoming increasingly promising. These advancements have not only led to new biophysical insights but also afforded opportunities to advance our understanding of fundamental intermolecular forces. This article describes the recent advances and applications, as well as future directions, of polarizable force fields in biomolecular simulations.

Mutations in Parkinson's Disease Associated Protein DJ-1 Alter the Energetics of DJ-1 Dimerization

Patients suffering from familial Parkinson's disease are linked to mutated DJ-1 protein. Wild-type DJ-1 occurs as a homodimer, which appears to be crucial for its function. It has been established that mutation (L166P) in DJ-1 protein could destabilize the DJ-1 homodimer. Hence, dimerization aspect of DJ-1 is fundamentally important for understanding its link to the disease. X-ray structures of wild-type DJ-1 dimer have given an atomic insight into the interaction network at the dimer interface. However, the energetics of dimerization in the wild-type and its mutant protein is unknown. Using the X-ray structure of wild-type DJ-1 as the template, we report ∼1.55 μs of molecular dynamics simulations to quantitatively estimate the relative free energy of DJ-1 dimerization in the disease linked variant (L166P, A104T, and M26I) with respect to its wild-type analogue. The results suggest that dimerization is disfavored for L166P and A104T mutations, severely for the former. Notably, the M26I mutation does not alter the energetics of DJ-1 dimerization. The dynamics of the DJ-1 dimer is significantly altered in response to the L166P and A104T mutations, resulting in the significant loss of interactions at the dimer interface. L166P mutant showed the structural difference and increased flexibility in α6, α7, α8 regions with respect to the WT. A structural difference in the α6 region was noticeable between WT and A104T mutant of DJ-1. The interaction network in the dimer interface is identical for the wild-type protein and the M26I mutant. No significant change in secondary structural content was observed for DJ-1 mutants (L166P, A104T, M26I) with respect to its WT analogue.

Accurate PDZ/Peptide Binding Specificity with Additive and Polarizable Free Energy Simulations

PDZ domains contain 80-100 amino acids and bind short C-terminal sequences of target proteins. Their specificity is essential for cellular signaling pathways. We studied the binding of the Tiam1 PDZ domain to peptides derived from the C-termini of its Syndecan-1 and Caspr4 targets. We used free energy perturbation (FEP) to characterize the binding energetics of one wild-type and 17 mutant complexes by simulating 21 alchemical transformations between pairs of complexes. Thirteen complexes had known experimental affinities. FEP is a powerful tool to understand protein/ligand binding. It depends, however, on the accuracy of molecular dynamics force fields and conformational sampling. Both aspects require continued testing, especially for ionic mutations. For six mutations that did not modify the net charge, we obtained excellent agreement with experiment using the additive, AMBER ff99SB force field, with a root mean square deviation (RMSD) of 0.37 kcal/mol. For six ionic mutations that modified the net charge, agreement was also good, with one large error (3 kcal/mol) and an RMSD of 0.9 kcal/mol for the other five. The large error arose from the overstabilization of a protein/peptide salt bridge by the additive force field. Four of the ionic mutations were also simulated with the polarizable Drude force field, which represents the first test of this force field for protein/ligand binding free energy changes. The large error was eliminated and the RMS error for the four mutations was reduced from 1.8 to 1.2 kcal/mol. The overall accuracy of FEP indicates it can be used to understand PDZ/peptide binding. Importantly, our results show that for ionic mutations in buried regions, electronic polarization plays a significant role.

Structural characterization of nucleotide 5'-triphosphates by infrared ion spectroscopy and theoretical studies

The molecular family of nucleotide triphosphates (NTPs), with adenosine 5'-triphosphate (ATP) as its best-known member, is of high biochemical importance as their phosphodiester bonds form Nature's main means to store and transport energy. Here, gas-phase IR spectroscopic studies and supporting theoretical studies have been performed on adenosine 5'-triphosphate, cytosine 5'-triphosphate and guanosine 5'-triphosphate to elucidate the intrinsic structural properties of NTPs, focusing on the influence of the nucleobase and the extent of deprotonation. Mass spectrometric studies involving collision induced dissociation showed similar fragmentation channels for the three studied NTPs within a selected charge state. The doubly charged anions exhibit fragmentation similar to the energy-releasing hydrolysis reaction in nature, while the singly charged anions show different dominant fragmentation channels, suggesting that the charge state plays a significant role in the favorability of the hydrolysis reaction. A combination of infrared ion spectroscopy and quantum-chemical computations indicates that the singly charged anions of all NTPs are preferentially deprotonated at their β-phosphates, while the doubly-charged anions are dominantly αβ-deprotonated. The assigned three-dimensional structure differs for ATP and CTP on the one hand and GTP on the other, in the sense that ATP and CTP show no interaction between nucleobase and phosphate tail, while in GTP they are hydrogen bonded. This can be rationalized by considering the structure and geometry of the NTPs where the final three dimensional structure depends on a subtle balance between hydrogen bond strength, flexibility and steric hindrance.

Classical Drude Polarizable Force Field Model for Methyl Phosphate and Its Interactions with Mg2

Phosphate groups are essential components of nucleic acids and proteins, whose interactions with solvent, metal ions, and ionic side chains help control folding and binding. Methyl phosphate (MP) represents a simple analog of phosphate moieties that are post-translation modifications in proteins and present at the termini of nucleic acids, among other environments. In the present study, we optimized parameters for use in polarizable molecular dynamics simulations of MP in its mono- and dianionic forms, MP- ≡ CH3HPO4- and MP2- ≡ CH3PO42-, along with P i2- ≡ HPO42-, in the context of the classical Drude oscillator model. Parameter optimization was done in a manner consistent with the remainder of the Drude molecular mechanics force field, choosing atomic charges and polarizabilities to reproduce molecular properties from quantum mechanics as well as experimental hydration free energies. Optimized parameters were similar to existing dimethyl phosphate parameters, with a few significant differences. The developed parameters were then used to compute magnesium binding affinities in aqueous solution, using alchemical molecular dynamics free energy simulations. Good agreement with experiment was obtained, and outer sphere binding was shown to be predominant for MP- and MP2-.

Energetics of Preferential Binding of Retinoic Acid-Inducible Gene-I to Double-Stranded Viral RNAs with 5' Tri-/Diphosphate over 5' Monophosphate

Retinoic acid-inducible gene-I (RIG-I) is a cytosolic sensor protein that recognizes viral RNAs and triggers an innate immune response in cells. Panhandle-like base-paired blunt-ended 5' ppp/pp-dsRNA is a characteristic feature of viral RNAs. Structural studies of RIG-I C-terminal domain bound 5' ppp/pp-dsRNA complexes show the direct interaction between all the 5' terminal phosphates (α, β, and γ) and protein, suggesting γ phosphate might be a major recognition determinant for RIG-I binding. Biochemical studies, however, suggest that 5' pp-dsRNA is the minimal determinant for RIG-I binding and antiviral response. Despite biochemical and structural studies, the origin of viral RNA recognition by RIG-I is an unsolved problem. X-ray structures of RIG-I bound dsRNA not only provide atomic insight into the interaction network but also provide sufficiently good models for computational studies. We report structure-based molecular dynamics (MD) free energy calculations to quantitatively estimate the energetics of RIG-I binding to dsRNA containing 5' ppp, 5' pp, and 5' p. The results suggest that RIG-I weakly discriminates between 5' ppp-dsRNA and 5' pp-dsRNA (favoring former) and strongly disfavors 5' p-dsRNA with respect to the rest. Interestingly, direct interaction between γ phosphate of 5' ppp-dsRNA and RIG-I is a robust feature of the MD simulations. dsRNA binding to RIG-I is associated with Mg²⁺ dissociation from the 5' phosphate/s of dsRNA. The higher Mg²⁺ dissociation penalty from 5' ppp-dsRNA with respect to 5' pp-dsRNA offsets most of the favorable interaction between RIG-I and γ phosphate of 5' ppp-dsRNA. This leads to weak discrimination between 5' ppp-dsRNA and 5' pp-dsRNA. 5' p-dsRNA is discriminated strongly because of the loss of interaction with RIG-I.

AMOEBA Polarizable Force Field Parameters of the Heme Cofactor in Its Ferrous and Ferric Forms

We report the first parameters of the heme redox cofactors for the polarizable AMOEBA force field in both the ferric and ferrous forms. We consider two types of complexes, one with two histidine side chains as axial ligands and one with a histidine and a methionine side chain as ligands. We have derived permanent multipoles from second-order Møller-Plesset perturbation theory (MP2). The sets of parameters have been validated in a first step by comparison of AMOEBA interaction energies of heme and a collection of biologically relevant molecules with MP2 and Density Functional Theory (DFT) calculations. In a second validation step, we consider interaction energies with large aggregates comprising around 80 H₂O molecules. These calculations are repeated for 30 structures extracted from semiempirical PM7 DM simulations. Very encouraging agreement is found between DFT and the AMOEBA force field, which results from an accurate treatment of electrostatic interactions. We finally report long (10 ns) MD simulations of cytochromes in two redox states with AMOEBA testing both the 2003 and 2014 AMOEBA water models. These simulations have been carried out with the TINKER-HP (High Performance) program. In conclusion, owing to their ubiquity in biology, we think the present work opens a wide array of applications of the polarizable AMOEBA force field on hemeproteins.

QM/MM investigation of ATP hydrolysis in aqueous solution

Adenosine-5'-triphosphate (ATP) hydrolysis represents a most important reaction in biology. Despite extensive research efforts, the mechanism for ATP hydrolysis in aqueous solution still remains under debate. Previous theoretical studies often predefined reaction coordinates to characterize the mechanism for ATP hydrolysis in water with Mg(2+) by evaluating free energy profiles through these preassumed reaction paths. In the present work, a nudged elastic band method is applied to identify the minimum energy path calculated with a hybrid quantum mechanics and molecular mechanics approach. Along the reaction path, the free energy profile was obtained to have a single transition state and the activation energy of 32.5 kcal/mol. This transition state bears a four-centered structure that describes a concerted nature of the reaction. In the More-O'Ferrall-Jencks diagram, the results show that the reaction proceeds through a concerted path before the system reaches the transition state and along an associative path after the transition state. In addition, the calculated reaction free energy is -7.0 kcal/mol, in good agreement with experiment, capturing the exothermic feature of MgATP(2-) hydrolysis in aqueous solution, whereas the reaction was often shown to be endothermic in the previous theoretical studies. As Mg(2+) is required for ATP hydrolysis in cells, its role in the reaction is also elucidated.

Allosteric response and substrate sensitivity in peptide binding of the signal recognition particle

We characterize the conformational dynamics and substrate selectivity of the signal recognition particle (SRP) using a thermodynamic free energy cycle approach and microsecond timescale molecular dynamics simulations. The SRP is a central component of the co-translational protein targeting machinery that binds to the N-terminal signal peptide (SP) of nascent proteins. We determined the shift in relative conformational stability of the SRP upon substrate binding to quantify allosteric coupling between SRP domains. In particular, for dipeptidyl aminopeptidase, an SP that is recognized by the SRP for co-translational targeting, it is found that substrate binding induces substantial changes in the SRP toward configurations associated with targeting of the nascent protein, and it is found that the changes are modestly enhanced by a mutation that increases the hydrophobicity of the SP. However, for alkaline phosphatase, an SP that is recognized for post-translational targeting, substrate binding induces the reverse change in the SRP conformational distribution away from targeting configurations. Microsecond timescale trajectories reveal the intrinsic flexibility of the SRP conformational landscape and provide insight into recent single molecule studies by illustrating that 10-nm lengthscale changes between FRET pairs occur via the rigid-body movement of SRP domains connected by the flexible linker region. In combination, these results provide direct evidence for the hypothesis that substrate-controlled conformational switching in the SRP provides a mechanism for discriminating between different SPs and for connecting substrate binding to downstream steps in the protein targeting pathway.

Electrostatic free energies in translational GTPases: Classic allostery and the rest

GTPases typically switch between an inactive, OFF conformation and an active, ON conformation when a GDP ligand is replaced by GTP. Their ON/OFF populations and activity thus depend on the stabilities of four protein complexes, two apo-protein forms, and GTP/GDP in solution. A complete characterization is usually not possible experimentally and poses major challenges for simulations. We review the most important methodological challenges and we review thermodynamic data for two GTPases involved in translation of the genetic code: archaeal Initiation Factors 2 and 5B (aIF2, aIF5B). One main challenge is the multiplicity of states and conformations, including those of GTP/GDP in solution. Another is force field accuracy, especially for interactions of GTP/GDP with co-bound divalent Mg(2+) ions. The calculation of electrostatic free energies also poses specific challenges, and requires careful protocols. For aIF2, experiments and earlier simulations showed that it is a "classic" GTPase, with distinct ON/OFF conformations that prefer to bind GTP and GDP, respectively. For aIF5B, we recently proposed a non-classic mechanism, where the ON/OFF states differ only in the protonation state of Glu81 in the nucleotide binding pocket. This model is characterized here using free energy simulations. The methodological analysis should help future studies, while the aIF2, aIF5B examples illustrate the diversity of ATPase/GTPase mechanisms. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

Simulating GTP:Mg and GDP:Mg with a simple force field: a structural and thermodynamic analysis

Di- and tri-phosphate nucleotides are essential cofactors for many proteins, usually in an Mg(2+) -bound form. Proteins like GTPases often detect the difference between NDP and NTP and respond by changing conformations. To study such complexes, simple, fixed charge force fields have been used, which allow long simulations and precise free energy calculations. The preference for NTP or NDP binding depends on many factors, including ligand structure and Mg(2+) coordination and the changes they undergo upon binding. Here, we use a simple force field to examine two Mg(2+) coordination modes for the unbound GDP and GTP: direct, or "Inner Sphere" (IS) coordination by one or more phosphate oxygens and indirect, "Outer Sphere" (OS) coordination involving one or more bridging waters. We compare GTP: and GDP:Mg binding with OS and IS coordination; combining the results with experimental data then indicates that GTP prefers the latter. We also examine different kinds of IS coordination and their sensitivity to a key force field parameter: the optimal Mg:oxygen van der Waals distance Rmin . Increasing Rmin improves the Mg:oxygen distances, the GTP: and GDP:Mg binding affinities, and the fraction of GTP:Mg with β + γ phosphate coordination, but does not improve or change the GTP/GDP affinity difference, which remains much larger than experiment. It has no effect on the free energy of GDP binding to a GTPase.

Nucleotide recognition by the initiation factor aIF5B: free energy simulations of a neoclassical GTPase

The GTPase aIF5B is a universally conserved initiation factor that assists ribosome assembly. Crystal structures of its nucleotide complexes, X-ray(GTP) and X-ray(GDP), are similar in the nucleotide vicinity, but differ in the orientation of a distant domain IV. This has led to two, contradictory, mechanistic models. One postulates that X-ray(GTP) and X-ray(GDP) are, respectively, the active, "ON" and the inactive, "OFF" states; the other postulates that both structures are OFF, whereas the ON state is still uncharacterized. We study GTP/GDP binding using molecular dynamics and a continuum electrostatic free energy method. We predict that X-ray(GTP) has a ≈ 3 kcal/mol preference to bind GDP, apparently contradicting its assignment as ON. However, the preference arises mainly from a single, nearby residue from the switch 2 motif: Glu81, which becomes protonated upon GTP binding, with a free energy cost of about 4 kcal/mol. We then propose a different model, where Glu81 protonation/deprotonation defines the ON/OFF states. With this model, the X-ray(GTP):GTP complex, with its protonated Glu81, is ON, whereas X-ray(GTP):GDP is OFF. The model postulates that distant conformational changes such as domain IV rotation are "uncoupled" from GTP/GDP exchange and do not affect the relative GTP/GDP binding affinities. We analyze the model using a general thermodynamic framework for GTPases. It yields rather precise predictions for the nucleotide specificities of each state, and the state specificities of each nucleotide, which are roughly comparable to the homologues IF2 and aIF2, despite the lack of any conformational switching in the model.

Conformational selection through electrostatics: Free energy simulations of GTP and GDP binding to archaeal initiation factor 2

Archaeal Initiation Factor 2 is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and their dependence on the conformational state, molecular dynamics free energy simulations were performed. The ON state specificity was predicted to be weak, with a GTP/GDP binding free energy difference of -1 kcal/mol, favoring GTP. The OFF state specificity is larger, 4 kcal/mol, favoring GDP. The overall effects result from a competition among many interactions in several complexes. To interpret them, we use a simpler, dielectric continuum model. Several effects are robust with respect to the model details. Both nucleotides have a net negative charge, so that removing them from solvent into the binding pocket carries a desolvation penalty, which is large for the ON state, and strongly disfavors GTP binding compared to GDP. Short-range interactions between the additional GTP phosphate group and ionized sidechains in the binding pocket offset most, but not all of the desolvation penalty; more distant groups also contribute significantly, and the switch 1 loop only slightly. The desolvation penalty is lower for the more open, wetter OFF state, and the GTP/GDP difference much smaller. Short-range interactions in the binding pocket and with more distant groups again make a significant contribution. Overall, the simulations help explain how conformational selection is achieved with a single phosphate group.

Conformational selection by the aIF2 GTPase: a molecular dynamics study of functional pathways

Archaeal initiation factor 2 (aIF2) is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To improve our understanding of the role of each conformational state in the aIF2 "life cycle", we start from the state immediately after GTP hydrolysis, ON:GDP:P(i) (where P(i) is inorganic phosphate), and consider the possible next steps on the pathway to the OFF:GDP product. The first possibility is P(i) dissociation, leading to ON:GDP, which could then relax into OFF:GDP. We use molecular dynamics simulations to compute the P(i) dissociation free energy and show that dissociation is highly favorable. The second possibility is conformational relaxation into the OFF state before P(i) dissociation, to form OFF:GDP:P(i). We estimate the corresponding free energy approximately, 2 ± 3.5 kcal/mol, so that this is an uphill or weakly downhill process. A third possibility is relaxation into another conformation, neither ON nor OFF. Indeed, a third, "MIXED" conformation was seen recently in a crystal structure of the aIF2:GDP:P(i) complex. For this conformational state, P(i) dissociation is weakly unfavorable, in contrast to the ON and OFF states. From this, we will deduce that if the MIXED:GDP complex is not too unstable, the ON:GDP:P(i) → MIXED:GDP:P(i) transformation is a downhill process, which can occur spontaneously. This suggests that the MIXED state could be a functional intermediate.