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

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.

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.

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.

Conformational transitions in the γ subunit of the archaeal translation initiation factor 2

In eukaryotes and archaea, the heterotrimeric translation initiation factor 2 (e/aIF2) is pivotal for the delivery of methionylated initiator tRNA (Met-tRNA(i)) to the ribosome. It acts as a molecular switch that cycles between inactive (GDP-bound) and active (GTP-bound) states. Recent studies show that eIF2 can also exist in a long-lived eIF2γ-GDP-P(i) (inorganic phosphate) active state. Here, four high-resolution crystal structures of aIF2γ from Sulfolobus solfataricus are reported: aIF2γ-GDPCP (a nonhydrolyzable GTP analogue), aIF2γ-GDP-formate (in which a formate ion possibly mimics P(i)), aIF2γ-GDP and nucleotide-free aIF2γ. The structures describe the different states of aIF2γ and demonstrate the conformational transitions that take place in the aIF2γ `life cycle'.

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.