Converting structural information into an allosteric-energy-based picture for elongation factor Tu activation by the ribosome

Molecular dynamics simulations for the structure-based drug design: targeting small-GTPases proteins

INTRODUCTION

Molecular Dynamics (MD) simulations can support mechanism-based drug design. Indeed, MD simulations by capturing biomolecule motions at finite temperatures can reveal hidden binding sites, accurately predict drug-binding poses, and estimate the thermodynamics and kinetics, crucial information for drug discovery campaigns. Small-Guanosine Triphosphate Phosphohydrolases (GTPases) regulate a cascade of signaling events, that affect most cellular processes. Their deregulation is linked to several diseases, making them appealing drug targets. The broad roles of small-GTPases in cellular processes and the recent approval of a covalent KRas inhibitor as an anticancer agent renewed the interest in targeting small-GTPase with small molecules.

AREA COVERED

This review emphasizes the role of MD simulations in elucidating small-GTPase mechanisms, assessing the impact of cancer-related variants, and discovering novel inhibitors.

EXPERT OPINION

The application of MD simulations to small-GTPases exemplifies the role of MD simulations in the structure-based drug design process for challenging biomolecular targets. Furthermore, AI and machine learning-enhanced MD simulations, coupled with the upcoming power of quantum computing, are promising instruments to target elusive small-GTPases mutations and splice variants. This powerful synergy will aid in developing innovative therapeutic strategies associated to small-GTPases deregulation, which could potentially be used for personalized therapies and in a tissue-agnostic manner to treat tumors with mutations in small-GTPases.

Catalytic mechanism study of ATP-citrate lyase during citryl-CoA synthesis process

ATP-citrate lyase (ACLY) is a critical metabolic enzyme and promising target for drug development. The structure determinations of ACLY have revealed its homotetramer states with various subunit symmetries, but catalytic mechanism of ACLY tetramer and the importance of subunit symmetry have not been clarified. Here, we constructed the free energy landscape of ACLY tetramer with arbitrary subunit symmetries and investigated energetic and conformational coupling of subunits during citryl-CoA synthesis process. The optimal conformational pathway indicates that ACLY tetramer encounters three critical conformational barriers and undergoes a loss of rigid-D2 symmetry to gain an energetic advantage. Energetic coupling of conformational changes and biochemical reactions suggests that these biological events are not independent but rather coupled with each other, showing a comparable energy barrier to the experimental data for the rate-limiting step. These findings could contribute to further research on catalytic mechanism, functional modulation, and inhibitor design of ACLY.

Application of Coarse-Grained (CG) Models to Explore Conformational Pathway of Large-Scale Protein Machines

Protein machines are clusters of protein assemblies that function in order to control the transfer of matter and energy in cells. For a specific protein machine, its working mechanisms are not only determined by the static crystal structures, but also related to the conformational transition dynamics and the corresponding energy profiles. With the rapid development of crystallographic techniques, the spatial scale of resolved structures is reaching up to thousands of residues, and the concomitant conformational changes become more and more complicated, posing a great challenge for computational biology research. Previously, a coarse-grained (CG) model aiming at conformational free energy evaluation was developed and showed excellent ability to reproduce the energy profiles by accurate electrostatic interaction calculations. In this study, we extended the application of the CG model to a series of large-scale protein machine systems. The spike protein trimer of SARS-CoV-2, ATP citrate lyase (ACLY) tetramer, and P4-ATPases systems were carefully studied and discussed as examples. It is indicated that the CG model is effective to depict the energy profiles of the conformational pathway between two endpoint structures, especially for large-scale systems. Both the energy change and energy barrier between endpoint structures provide reasonable mechanism explanations for the associated biological processes, including the opening of receptor binding domain (RBD) of spike protein, the phospholipid transportation of P4-ATPase, and the loop translocation of ACLY. Taken together, the CG model provides a suitable alternative in mechanistic studies related to conformational change in large-scale protein machines.

The energetics of subunit rotation in the ribosome

Protein synthesis in the cell is controlled by an elaborate sequence of conformational rearrangements in the ribosome. The composition of a ribosome varies by species, though they typically contain ∼ 50-100 RNA and protein molecules. While advances in structural techniques have revolutionized our understanding of long-lived conformational states, a vast range of transiently visited configurations can not be directly observed. In these cases, computational/simulation methods can be used to understand the mechanical properties of the ribosome. Insights from these approaches can then help guide next-generation experimental measurements. In this short review, we discuss theoretical strategies that have been deployed to quantitatively describe the energetics of collective rearrangements in the ribosome. We focus on efforts to probe large-scale subunit rotation events, which involve the coordinated displacement of large numbers of atoms (tens of thousands). These investigations are revealing how the molecular structure of the ribosome encodes the mechanical properties that control large-scale dynamics.

Peculiarities in Activation of Hydrolytic Activity of Elongation Factors

Translational GTPases (trGTPases) belong to the family of G proteins and play key roles at all stages of protein biosynthesis on the ribosome. Unidirectional and cyclic functioning of G proteins is ensured by their ability to switch between the active and inactive states due to GTP hydrolysis accelerated by the auxiliary GTPase-activating proteins. Although trGTPases interact with the ribosomes in different conformational states, they bind to the same conserved region, which, unlike in classical GTPase-activating proteins, is represented by ribosomal RNA. The resulting catalytic sites have almost identical structure in all elongation factors suggesting a common mechanism of GTP hydrolysis. However, fine details of the activated state formation and significantly different rates of GTP hydrolysis indicate the existence of distinctive features upon GTP hydrolysis catalyzed by the different factors. Here, we present a contemporary view on the mechanism of GTPase activation and GTP hydrolysis by the elongation factors EF-Tu, EF-G, and SelB based on the analysis of structural, biochemical, and bioinformatics data.

A steric gate controls P/E hybrid-state formation of tRNA on the ribosome

The ribosome is a biomolecular machine that undergoes multiple large-scale structural rearrangements during protein elongation. Here, we focus on a conformational rearrangement during translocation, known as P/E hybrid-state formation. Using a model that explicitly represents all non-hydrogen atoms, we simulated more than 120 spontaneous transitions, where the tRNA molecule is displaced between the P and E sites of the large subunit. In addition to predicting a free-energy landscape that is consistent with previous experimental observations, the simulations reveal how a six-residue gate-like region can limit P/E formation, where sub-angstrom structural perturbations lead to an order-of-magnitude change in kinetics. Thus, this precisely defined set of residues represents a novel target that may be used to control functional dynamics in bacterial ribosomes. This theoretical analysis establishes a direct relationship between ribosome structure and large-scale dynamics, and it suggests how next-generation experiments may precisely dissect the energetics of hybrid formation on the ribosome.

The ribosome undergoes multiple large-scale structural rearrangements during protein elongation. Here the authors present an all-atom model of the ribosome to study the energetics of P/E hybrid-state formation, an early conformational rearrangement occurring during translocation.

Structural basis of elongation factor 2 switching

Archaebacterial and eukaryotic elongation factor 2 (EF-2) and bacterial elongation factor G (EF-G) are five domain GTPases that catalyze the ribosomal translocation of tRNA and mRNA. In the classical mechanism of activation, GTPases are switched on through GDP/GTP exchange, which is accompanied by the ordering of two flexible segments called switch I and II. However, crystal structures of EF-2 and EF-G have thus far not revealed the conformations required by the classical mechanism. Here, we describe crystal structures of Methanoperedens nitroreducens EF-2 (MnEF-2) and MnEF-2-H595N bound to GMPPCP (GppCp) and magnesium displaying previously unreported compact conformations. Domain III forms interfaces with the other four domains and the overall conformations resemble that of SNU114, the eukaryotic spliceosomal GTPase. The gamma phosphate of GMPPCP is detected through interactions with switch I and a P-loop structural element. Switch II is highly ordered whereas switch I shows a variable degree of ordering. The ordered state results in a tight interdomain arrangement of domains I-III and the formation of a portion of a predicted monovalent cation site involving the P-loop and switch I. The side chain of an essential histidine residue in switch II is placed in the inactive conformation observed for the “on” state of elongation factor EF-Tu. The compact conformations of MnEF-2 and MnEF-2-H595N suggest an “on” ribosome-free conformational state.

•

Crystal structures of ribosome-free elongation factor 2 (EF-2) bound to GTP analog and magnesium.

•

Compact conformation and P-loop, switch I, and switch II structures suggest “on” state.

•

Arrangement of domains I-III similar to that of ribosome-bound EF-2/EF-G complexed with GTP analog.

•

Switch II histidine shows inactive conformation observed for “on” state of ribosome-free EF-Tu.

tRNA Dissociation from EF-Tu after GTP Hydrolysis: Primary Steps and Antibiotic Inhibition

In each round of ribosomal translation, the translational GTPase elongation factor Tu (EF-Tu) delivers a transfer RNA (tRNA) to the ribosome. After successful decoding, EF-Tu hydrolyzes GTP, which triggers a conformational change that ultimately results in the release of the tRNA from EF-Tu. To identify the primary steps of these conformational changes and how they are prevented by the antibiotic kirromycin, we employed all-atom explicit-solvent molecular dynamics simulations of the full ribosome-EF-Tu complex. Our results suggest that after GTP hydrolysis and Pi release, the loss of interactions between the nucleotide and the switch 1 loop of EF-Tu allows domain D1 of EF-Tu to rotate relative to domains D2 and D3 and leads to an increased flexibility of the switch 1 loop. This rotation induces a closing of the D1-D3 interface and an opening of the D1-D2 interface. We propose that the opening of the D1-D2 interface, which binds the CCA tail of the tRNA, weakens the crucial EF-Tu-tRNA interactions, which lowers tRNA binding affinity, representing the first step of tRNA release. Kirromycin binds within the D1-D3 interface, sterically blocking its closure, but does not prevent hydrolysis. The resulting increased flexibility of switch 1 explains why it is not resolved in kirromycin-bound structures.

Cryo-EM shows stages of initial codon selection on the ribosome by aa-tRNA in ternary complex with GTP and the GTPase-deficient EF-TuH84A

The GTPase EF-Tu in ternary complex with GTP and aminoacyl-tRNA (aa-tRNA) promotes rapid and accurate delivery of cognate aa-tRNAs to the ribosomal A site. Here we used cryo-EM to study the molecular origins of the accuracy of ribosome-aided recognition of a cognate ternary complex and the accuracy-amplifying role of the monitoring bases A1492, A1493 and G530 of the 16S rRNA. We used the GTPase-deficient EF-Tu variant H84A with native GTP, rather than non-cleavable GTP analogues, to trap a near-cognate ternary complex in high-resolution ribosomal complexes of varying codon-recognition accuracy. We found that ribosome complexes trapped by GTPase-deficicent ternary complex due to the presence of EF-Tu^H84A or non-cleavable GTP analogues have very similar structures. We further discuss speed and accuracy of initial aa-tRNA selection in terms of conformational changes of aa-tRNA and stepwise activation of the monitoring bases at the decoding center of the ribosome.

GTP Hydrolysis Without an Active Site Base: A Unifying Mechanism for Ras and Related GTPases

GTP hydrolysis is a biologically crucial reaction, being involved in regulating almost all cellular processes. As a result, the enzymes that catalyze this reaction are among the most important drug targets. Despite their vital importance and decades of substantial research effort, the fundamental mechanism of enzyme-catalyzed GTP hydrolysis by GTPases remains highly controversial. Specifically, how do these regulatory proteins hydrolyze GTP without an obvious general base in the active site to activate the water molecule for nucleophilic attack? To answer this question, we perform empirical valence bond simulations of GTPase-catalyzed GTP hydrolysis, comparing solvent- and substrate-assisted pathways in three distinct GTPases, Ras, Rab, and the G_αi subunit of a heterotrimeric G-protein, both in the presence and in the absence of the corresponding GTPase activating proteins. Our results demonstrate that a general base is not needed in the active site, as the preferred mechanism for GTP hydrolysis is a conserved solvent-assisted pathway. This pathway involves the rate-limiting nucleophilic attack of a water molecule, leading to a short-lived intermediate that tautomerizes to form H₂PO₄^- and GDP as the final products. Our fundamental biochemical insight into the enzymatic regulation of GTP hydrolysis not only resolves a decades-old mechanistic controversy but also has high relevance for drug discovery efforts. That is, revisiting the role of oncogenic mutants with respect to our mechanistic findings would pave the way for a new starting point to discover drugs for (so far) "undruggable" GTPases like Ras.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

How Ricin Damages the Ribosome

Ricin belongs to the group of ribosome-inactivating proteins (RIPs), i.e., toxins that have evolved to provide particular species with an advantage over other competitors in nature. Ricin possesses RNA N-glycosidase activity enabling the toxin to eliminate a single adenine base from the sarcin-ricin RNA loop (SRL), which is a highly conserved structure present on the large ribosomal subunit in all species from the three domains of life. The SRL belongs to the GTPase associated center (GAC), i.e., a ribosomal element involved in conferring unidirectional trajectory for the translational apparatus at the expense of GTP hydrolysis by translational GTPases (trGTPases). The SRL represents a critical element in the GAC, being the main triggering factor of GTP hydrolysis by trGTPases. Enzymatic removal of a single adenine base at the tip of SRL by ricin blocks GTP hydrolysis and, at the same time, impedes functioning of the translational machinery. Here, we discuss the consequences of SRL depurination by ricin for ribosomal performance, with emphasis on the mechanistic model overview of the SRL modus operandi.

Exploring allosteric communication in multiple states of the bacterial ribosome using residue network analysis

Antibiotic resistance is one of the most important problems of our era and hence the discovery of new effective therapeutics is urgent. At this point, studying the allosteric communication pathways in the bacterial ribosome and revealing allosteric sites/residues is critical for designing new inhibitors or repurposing readily approved drugs for this enormous machine. To shed light onto molecular details of the allosteric mechanisms, here we construct residue networks of the bacterial ribosomal complex at four different states of translation by using an effective description of the intermolecular interactions. Centrality analysis of these networks highlights the functional roles of structural components and critical residues on the ribosomal complex. High betweenness scores reveal pathways of residues connecting numerous sites on the structure. Interestingly, these pathways assemble highly conserved residues, drug binding sites, and known allosterically linked regions on the same structure. This study proposes a new residue-level model to test how distant sites on the molecular machine may be linked through hub residues that are critically located on the contact topology to inherently form communication pathways. Findings also indicate intersubunit bridges B1b, B3, B5, B7, and B8 as critical targets to design novel antibiotics.

Translation in Prokaryotes

This review summarizes our current understanding of translation in prokaryotes, focusing on the mechanistic and structural aspects of each phase of translation: initiation, elongation, termination, and ribosome recycling. The assembly of the initiation complex provides multiple checkpoints for messenger RNA (mRNA) and start-site selection. Correct codon-anticodon interaction during the decoding phase of elongation results in major conformational changes of the small ribosomal subunit and shapes the reaction pathway of guanosine triphosphate (GTP) hydrolysis. The ribosome orchestrates proton transfer during peptide bond formation, but requires the help of elongation factor P (EF-P) when two or more consecutive Pro residues are to be incorporated. Understanding the choreography of transfer RNA (tRNA) and mRNA movements during translocation helps to place the available structures of translocation intermediates onto the time axis of the reaction pathway. The nascent protein begins to fold cotranslationally, in the constrained space of the polypeptide exit tunnel of the ribosome. When a stop codon is reached at the end of the coding sequence, the ribosome, assisted by termination factors, hydrolyzes the ester bond of the peptidyl-tRNA, thereby releasing the nascent protein. Following termination, the ribosome is dissociated into subunits and recycled into another round of initiation. At each step of translation, the ribosome undergoes dynamic fluctuations between different conformation states. The aim of this article is to show the link between ribosome structure, dynamics, and function.

EF-Tu and EF-G are activated by allosteric effects

Many cellular processes are controlled by GTPases, and gaining quantitative understanding of the activation of such processes has been a major challenge. In particular, it is crucial to obtain reliable free-energy surfaces for the relevant reaction paths both in solution and in GTPases active sites. Here, we revisit the energetics of the activation of EF-G and EF-Tu by the ribosome and explore the nature of the catalysis of the GTPase reaction. The comparison of EF-Tu to EF-G allows us to explore the impact of possible problems with the available structure of EF-Tu. Additionally, mutational effects are used for a careful validation of the emerging conclusions. It is found that the reaction may proceed by both a two-water mechanism and a one-water (GTP as a base) mechanism. However, in both cases, the activation involves a structural allosteric effect, which is likely to be a general-activation mechanism for all GTPases.

Computer simulations of the catalytic mechanism of wild-type and mutant β-phosphoglucomutase

β-Phosphoglucomutase (β-PGM) has served as an important model system for understanding biological phosphoryl transfer.

Ribosome dynamics during decoding

Elongation factors Tu (EF-Tu) and SelB are translational GTPases that deliver aminoacyl-tRNAs (aa-tRNAs) to the ribosome. In each canonical round of translation elongation, aa-tRNAs, assisted by EF-Tu, decode mRNA codons and insert the respective amino acid into the growing peptide chain. Stop codons usually lead to translation termination; however, in special cases UGA codons are recoded to selenocysteine (Sec) with the help of SelB. Recruitment of EF-Tu and SelB together with their respective aa-tRNAs to the ribosome is a multistep process. In this review, we summarize recent progress in understanding the role of ribosome dynamics in aa-tRNA selection. We describe the path to correct codon recognition by canonical elongator aa-tRNA and Sec-tRNA^Sec and discuss the local and global rearrangements of the ribosome in response to correct and incorrect aa-tRNAs. We present the mechanisms of GTPase activation and GTP hydrolysis of EF-Tu and SelB and summarize what is known about the accommodation of aa-tRNA on the ribosome after its release from the elongation factor. We show how ribosome dynamics ensures high selectivity for the cognate aa-tRNA and suggest that conformational fluctuations, induced fit and kinetic discrimination play major roles in maintaining the speed and fidelity of translation.

This article is part of the themed issue ‘Perspectives on the ribosome’.

The Conformational Change in Elongation Factor Tu Involves Separation of Its Domains

Elongation factor Tu (EF-Tu) is a highly conserved GTPase that is responsible for supplying the aminoacylated tRNA to the ribosome. Upon binding to the ribosome, EF-Tu undergoes GTP hydrolysis, which drives a major conformational change, triggering the release of aminoacylated tRNA to the ribosome. Using a combination of molecular simulation techniques, we studied the transition between the pre- and post-hydrolysis structures through two distinct pathways. We show that the transition free energy is minimal along a non-intuitive pathway that involves "separation" of the GTP binding domain (domain 1) from the OB folds (domains 2 and 3), followed by domain 1 rotation, and, eventually, locking the EF-Tu conformation in the post-hydrolysis state. The domain separation also leads to a slight extension of the linker connecting domain 1 to domain 2. Using docking tools and correlation-based analysis, we identified and characterized the EF-Tu conformations that release the tRNA. These calculations suggest that EF-Tu can release the tRNA before the domains separate and after domain 1 rotates by 25°. We also examined the EF-Tu conformations in the context of the ribosome. Given the high degrees of sequence similarity with other translational GTPases, we predict a similar separation mechanism is followed.

The effect of magnesium ions on triphosphate hydrolysis

The role of metal ions in catalyzing phosphate ester hydrolysis has been the subject of much debate, both in terms of whether they change the transition state structure or mechanistic pathway. Understanding the impact of metal ions on these biologically critical reactions is central to improving our understanding of the role of metal ions in the numerous enzymes that facilitate them. In the present study, we have performed density functional theory studies of the mechanisms of methyl triphosphate and acetyl phosphate hydrolysis in aqueous solution to explore the competition between solvent- and substrate-assisted pathways, and examined the impact of Mg2+ on the energetics and transition state geometries. In both cases, we observe a clear preference for a more dissociative solvent-assisted transition state, which is not significantly changed by coordination of Mg2+. The effect of Mg2+ on the transition state geometries for the two pathways is minimal. While our calculations cannot rule out a substrate-assisted pathway as a possible solution for biological phosphate hydrolysis, they demonstrate that a significantly higher energy barrier needs to be overcome in the enzymatic reaction for this to be an energetically viable reaction pathway.

High-Dimensional Mutant and Modular Thermodynamic Cycles, Molecular Switching, and Free Energy Transduction

Understanding how distinct parts of proteins produce coordinated behavior has driven and continues to drive advances in protein science and enzymology. However, despite consensus about the conceptual basis for allostery, the idiosyncratic nature of allosteric mechanisms resists general approaches. Computational methods can identify conformational transition states from structural changes, revealing common switching mechanisms that impose multistate behavior. Thermodynamic cycles use factorial perturbations to measure coupling energies between side chains in molecular switches that mediate shear during domain motion. Such cycles have now been complemented by modular cycles that measure energetic coupling between separable domains. For one model system, energetic coupling between domains has been shown to be quantitatively equivalent to that between dynamic side chains. Linkages between domain motion, switching residues, and catalysis make nucleoside triphosphate hydrolysis conditional on domain movement, confirming an essential yet neglected aspect of free energy transduction and suggesting the potential generality of these studies.

Regulation and Plasticity of Catalysis in Enzymes: Insights from Analysis of Mechanochemical Coupling in Myosin

The mechanism of ATP hydrolysis in the myosin motor domain is analyzed using a combination of DFTB3/CHARMM simulations and enhanced sampling techniques. The motor domain is modeled in the pre-powerstroke state, in the post-rigor state, and as a hybrid based on the post-rigor state with a closed nucleotide-binding pocket. The ATP hydrolysis activity is found to depend on the positioning of nearby water molecules, and a network of polar residues facilitates proton transfer and charge redistribution during hydrolysis. Comparison of the observed hydrolysis pathways and the corresponding free energy profiles leads to detailed models for the mechanism of ATP hydrolysis in the pre-powerstroke state and proposes factors that regulate the hydrolysis activity in different conformational states. In the pre-powerstroke state, the scissile P_γ-O_3β bond breaks early in the reaction. Proton transfer from the lytic water to the γ-phosphate through active site residues is an important part of the kinetic bottleneck; several hydrolysis pathways that feature distinct proton transfer routes are found to have similar free energy barriers, suggesting a significant degree of plasticity in the hydrolysis mechanism. Comparison of hydrolysis in the pre-powerstroke state and the closed post-rigor model suggests that optimization of residues beyond the active site for electrostatic stabilization and preorganization is likely important to enzyme design.

Review: Translational GTPases

Translational GTPases (trGTPases) play key roles in facilitating protein synthesis on the ribosome. Despite the high degree of evolutionary conservation in the sequences of their GTP-binding domains, the rates of GTP hydrolysis and nucleotide exchange vary broadly between different trGTPases. EF-Tu, one of the best-characterized model G proteins, evolved an exceptionally rapid and tightly regulated GTPase activity, which ensures rapid and accurate incorporation of amino acids into the nascent chain. Other trGTPases instead use the energy of GTP hydrolysis to promote movement or to ensure the forward commitment of translation reactions. Recent data suggest the GTPase mechanism of EF-Tu and provide an insight in the catalysis of GTP hydrolysis by its unusual activator, the ribosome. Here we summarize these advances in understanding the functional cycle and the regulation of trGTPases, stimulated by the elucidation of their structures on the ribosome and the progress in dissecting the reaction mechanism of GTPases. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 463-475, 2016.

The GTPase hGBP1 converts GTP to GMP in two steps via proton shuttle mechanisms

GTPases play a crucial role in the regulation of many biological processes by catalyzing the hydrolysis of GTP into GDP. The focus of this work is on the dynamin-related large GTPase human guanine nucleotide binding protein-1 (hGBP1) which is able to hydrolyze GTP even to GMP. Here, we studied the largely unknown mechanisms of both GTP and GDP hydrolysis steps utilizing accelerated ab initio QM/MM metadynamics simulations to compute multi-dimensional free energy landscapes. We find an indirect substrate-assisted catalysis (SAC) mechanism for GTP hydrolysis involving transfer of a proton from the water nucleophile to a nonbridging phosphoryl oxygen via a proton relay pathway where the rate-determining first step is concerted-dissociative nature. A "composite base" consisting of Ser73, Glu99, a bridging water molecule, and GTP was found to activate the nucleophilic water, thus disclosing the complex nature of the general base in hGBP1. A nearly two-fold reduction in the free energy barrier was obtained for GTP hydrolysis in the enzyme in comparison to bulk solvent. The subsequent GDP hydrolysis in hGBP1 was also found to follow a water-mediated proton shuttle mechanism. It is expected that the proton shuttle mechanisms unravelled for hGBP1 apply to many classes of GTPases/ATPases that possess an optimally-arranged hydrogen bonding network, which connects the catalytic water to a proton acceptor.

The pathway to GTPase activation of elongation factor SelB on the ribosome

In all domains of life, selenocysteine (Sec) is delivered to the ribosome by selenocysteine-specific tRNA (tRNA^Sec) with the help of a specialized translation factor, SelB in bacteria. Sec-tRNA^Sec recodes a UGA stop codon next to a downstream mRNA stem-loop. Here we present the structures of six intermediates on the pathway of UGA recoding in Escherichia coli by single-particle cryo-electron microscopy. The structures explain the specificity of Sec-tRNA^Sec binding by SelB and show large-scale rearrangements of Sec-tRNA^Sec. Upon initial binding of SelB-Sec-tRNA^Sec to the ribosome and codon reading, the 30S subunit adopts an open conformation with Sec-tRNA^Sec covering the sarcin-ricin loop (SRL) on the 50S subunit. Subsequent codon recognition results in a local closure of the decoding site, which moves Sec-tRNA^Sec away from the SRL and triggers a global closure of the 30S subunit shoulder domain. As a consequence, SelB docks on the SRL, activating the GTPase of SelB. These results reveal how codon recognition triggers GTPase activation in translational GTPases.

Structural characterization of ribosome recruitment and translocation by type IV IRES

Viral mRNA sequences with a type IV IRES are able to initiate translation without any host initiation factors. Initial recruitment of the small ribosomal subunit as well as two translocation steps before the first peptidyl transfer are essential for the initiation of translation by these mRNAs. Using electron cryomicroscopy (cryo-EM) we have structurally characterized at high resolution how the Cricket Paralysis Virus Internal Ribosomal Entry Site (CrPV-IRES) binds the small ribosomal subunit (40S) and the translocation intermediate stabilized by elongation factor 2 (eEF2). The CrPV-IRES restricts the otherwise flexible 40S head to a conformation compatible with binding the large ribosomal subunit (60S). Once the 60S is recruited, the binary CrPV-IRES/80S complex oscillates between canonical and rotated states (Fernández et al., 2014; Koh et al., 2014), as seen for pre-translocation complexes with tRNAs. Elongation factor eEF2 with a GTP analog stabilizes the ribosome-IRES complex in a rotated state with an extra ~3 degrees of rotation. Key residues in domain IV of eEF2 interact with pseudoknot I (PKI) of the CrPV-IRES stabilizing it in a conformation reminiscent of a hybrid tRNA state. The structure explains how diphthamide, a eukaryotic and archaeal specific post-translational modification of a histidine residue of eEF2, is involved in translocation.

DOI:http://dx.doi.org/10.7554/eLife.13567.001

Steric interactions lead to collective tilting motion in the ribosome during mRNA-tRNA translocation

Translocation of mRNA and tRNA through the ribosome is associated with large-scale rearrangements of the head domain in the 30S ribosomal subunit. To elucidate the relationship between 30S head dynamics and mRNA–tRNA displacement, we apply molecular dynamics simulations using an all-atom structure-based model. Here we provide a statistical analysis of 250 spontaneous transitions between the A/P–P/E and P/P–E/E ensembles. Consistent with structural studies, the ribosome samples a chimeric ap/P–pe/E intermediate, where the 30S head is rotated ∼18°. It then transiently populates a previously unreported intermediate ensemble, which is characterized by a ∼10° tilt of the head. To identify the origins of head tilting, we analyse 781 additional simulations in which specific steric features are perturbed. These calculations show that head tilting may be attributed to specific steric interactions between tRNA and the 30S subunit (PE loop and protein S13). Taken together, this study demonstrates how molecular structure can give rise to large-scale collective rearrangements.

During protein elongation, the translocation of mRNA and tRNA molecules across the 30S ribosomal subunit is associated with large-scale motions of the 30S head domain. Here the authors carry out MD simulations to probe the associated steric interactions and identify novel tilting motions during the late stages of translocation.

Theoretical Insights on the Mechanism of the GTP Hydrolysis Catalyzed by the Elongation Factor Tu (EF-Tu)

The purpose of this work is to have a better understanding of the mechanism of GTP hydrolysis catalyzed by the elongation factor Tu. Two main aspects are being discussed in the literature: the associative or dissociative character of the process and the nature of nucleophile activation. The calculations of the QM subsystem have been done by means of the M06-2X density functional and the split valence triple-ζ 6-311+G(d,p) basis set. The environmental effect has been introduced through the continuum SMD method. We have studied three models of increasing complexity in order to analyze the different factors that intervene in the catalytic action. The results obtained in this paper confirm that the protonated His84 plays a fundamental role in the catalytic mechanism, but we have also found that the crystallographic sodium ion has a notable effect in the catalysis. So, our work has permitted a new insight, complementary to those obtained with QM/MM calculations, into this very complex process.

Exceptionally large entropy contributions enable the high rates of GTP hydrolysis on the ribosome

Protein synthesis on the ribosome involves hydrolysis of GTP in several key steps of the mRNA translation cycle. These steps are catalyzed by the translational GTPases of which elongation factor Tu (EF-Tu) is the fastest GTPase known. Here, we use extensive computer simulations to explore the origin of its remarkably high catalytic rate on the ribosome and show that it is made possible by a very large positive activation entropy. This entropy term (TΔS(‡)) amounts to more than 7 kcal/mol at 25 °C. It is further found to be characteristic of the reaction mechanism utilized by the translational, but not other, GTPases and it enables these enzymes to attain hydrolysis rates exceeding 500 s(-1). This entropy driven mechanism likely reflects the very high selection pressure on the speed of protein synthesis, which drives the rate of each individual GTPase towards maximal turnover rate of the whole translation cycle.

Cooperative Electrostatic Interactions Drive Functional Evolution in the Alkaline Phosphatase Superfamily

It is becoming widely accepted that catalytic promiscuity, i.e., the ability of a single enzyme to catalyze the turnover of multiple, chemically distinct substrates, plays a key role in the evolution of new enzyme functions. In this context, the members of the alkaline phosphatase superfamily have been extensively studied as model systems in order to understand the phenomenon of enzyme multifunctionality. In the present work, we model the selectivity of two multiply promiscuous members of this superfamily, namely the phosphonate monoester hydrolases from Burkholderia caryophylli and Rhizobium leguminosarum. We have performed extensive simulations of the enzymatic reaction of both wild-type enzymes and several experimentally characterized mutants. Our computational models are in agreement with key experimental observables, such as the observed activities of the wild-type enzymes, qualitative interpretations of experimental pH-rate profiles, and activity trends among several active site mutants. In all cases the substrates of interest bind to the enzyme in similar conformations, with largely unperturbed transition states from their corresponding analogues in aqueous solution. Examination of transition-state geometries and the contribution of individual residues to the calculated activation barriers suggest that the broad promiscuity of these enzymes arises from cooperative electrostatic interactions in the active site, allowing each enzyme to adapt to the electrostatic needs of different substrates. By comparing the structural and electrostatic features of several alkaline phosphatases, we suggest that this phenomenon is a generalized feature driving selectivity and promiscuity within this superfamily and can be in turn used for artificial enzyme design.

Using Molecular Simulation to Model High-Resolution Cryo-EM Reconstructions

An explosion of new data from high-resolution cryo-electron microscopy (cryo-EM) studies has produced a large number of data sets for many species of ribosomes in various functional states over the past few years. While many methods exist to produce structural models for lower resolution cryo-EM reconstructions, high-resolution reconstructions are often modeled using crystallographic techniques and extensive manual intervention. Here, we present an automated fitting technique for high-resolution cryo-EM data sets that produces all-atom models highly consistent with the EM density. Using a molecular dynamics approach, atomic positions are optimized with a potential that includes the cross-correlation coefficient between the structural model and the cryo-EM electron density, as well as a biasing potential preserving the stereochemistry and secondary structure of the biomolecule. Specifically, we use a hybrid structure-based/ab initio molecular dynamics potential to extend molecular dynamics fitting. In addition, we find that simulated annealing integration, as opposed to straightforward molecular dynamics integration, significantly improves performance. We obtain atomistic models of the human ribosome consistent with high-resolution cryo-EM reconstructions of the human ribosome. Automated methods such as these have the potential to produce atomistic models for a large number of ribosome complexes simultaneously that can be subsequently refined manually.

Role of a ribosomal RNA phosphate oxygen during the EF-G-triggered GTP hydrolysis

Elongation factor-catalyzed GTP hydrolysis is a key reaction during the ribosomal elongation cycle. Recent crystal structures of G proteins, such as elongation factor G (EF-G) bound to the ribosome, as well as many biochemical studies, provide evidence that the direct interaction of translational GTPases (trGTPases) with the sarcin-ricin loop (SRL) of ribosomal RNA (rRNA) is pivotal for hydrolysis. However, the precise mechanism remains elusive and is intensively debated. Based on the close proximity of the phosphate oxygen of A2662 of the SRL to the supposedly catalytic histidine of EF-G (His87), we probed this interaction by an atomic mutagenesis approach. We individually replaced either of the two nonbridging phosphate oxygens at A2662 with a methyl group by the introduction of a methylphosphonate instead of the natural phosphate in fully functional, reconstituted bacterial ribosomes. Our major finding was that only one of the two resulting diastereomers, the SP methylphosphonate, was compatible with efficient GTPase activation on EF-G. The same trend was observed for a second trGTPase, namely EF4 (LepA). In addition, we provide evidence that the negative charge of the A2662 phosphate group must be retained for uncompromised activity in GTP hydrolysis. In summary, our data strongly corroborate that the nonbridging proSP phosphate oxygen at the A2662 of the SRL is critically involved in the activation of GTP hydrolysis. A mechanistic scenario is supported in which positioning of the catalytically active, protonated His87 through electrostatic interactions with the A2662 phosphate group and H-bond networks are key features of ribosome-triggered activation of trGTPases.

Structural Insights into tRNA Dynamics on the Ribosome

High-resolution structures at different stages, as well as biochemical, single molecule and computational approaches have highlighted the elasticity of tRNA molecules when bound to the ribosome. It is well acknowledged that the inherent structural flexibility of the tRNA lies at the heart of the protein synthesis process. Here, we review the recent advances and describe considerations that the conformational changes of the tRNA molecules offer about the mechanisms grounded in translation.

Modeling the mechanisms of biological GTP hydrolysis

Enzymes that hydrolyze GTP are currently in the spotlight, due to their molecular switch mechanism that controls many cellular processes. One of the best-known classes of these enzymes are small GTPases such as members of the Ras superfamily, which catalyze the hydrolysis of the γ-phosphate bond in GTP. In addition, the availability of an increasing number of crystal structures of translational GTPases such as EF-Tu and EF-G have made it possible to probe the molecular details of GTP hydrolysis on the ribosome. However, despite a wealth of biochemical, structural and computational data, the way in which GTP hydrolysis is activated and regulated is still a controversial topic and well-designed simulations can play an important role in resolving and rationalizing the experimental data. In this review, we discuss the contributions of computational biology to our understanding of GTP hydrolysis on the ribosome and in small GTPases.

Identification of a second GTP-bound magnesium ion in archaeal initiation factor 2

Eukaryotic and archaeal translation initiation processes involve a heterotrimeric GTPase e/aIF2 crucial for accuracy of start codon selection. In eukaryotes, the GTPase activity of eIF2 is assisted by a GTPase-activating protein (GAP), eIF5. In archaea, orthologs of eIF5 are not found and aIF2 GTPase activity is thought to be non-assisted. However, no in vitro GTPase activity of the archaeal factor has been reported to date. Here, we show that aIF2 significantly hydrolyses GTP in vitro. Within aIF2γ, H97, corresponding to the catalytic histidine found in other translational GTPases, and D19, from the GKT loop, both participate in this activity. Several high-resolution crystal structures were determined to get insight into GTP hydrolysis by aIF2γ. In particular, a crystal structure of the H97A mutant was obtained in the presence of non-hydrolyzed GTP. This structure reveals the presence of a second magnesium ion bound to GTP and D19. Quantum chemical/molecular mechanical simulations support the idea that the second magnesium ion may assist GTP hydrolysis by helping to neutralize the developing negative charge in the transition state. These results are discussed in light of the absence of an identified GAP in archaea to assist GTP hydrolysis on aIF2.

The conformation of a catalytic loop is central to GTPase activity on the ribosome

The translational GTPases hydrolyze GTP on the ribosome at several stages of the protein synthesis cycle. Because of the strong conservation of their catalytic center, these enzymes are expected to operate through a universal hydrolysis mechanism, in which a critical histidine residue together with the sarcin-ricin loop of the large ribosomal subunit is necessary for GTPase activation. Here we examine different possible pathways for GTP hydrolysis by EF-Tu through extensive computer simulations. We show that a conformational change of the peptide plane preceding this histidine has a decisive effect on the energetics of the reaction. This transition was predicted earlier by us and has recently been confirmed experimentally. It is found to promote early proton transfer from water to the γ-phosphate group of GTP, followed by nucleophilic attack by hydroxide ion. The calculated reaction energetics is in good agreement with available kinetic data, for both wild-type and mutant versions of EF-Tu, and indicates that the latter may enforce a change in mechanism toward more concerted pathways.

Resolving apparent conflicts between theoretical and experimental models of phosphate monoester hydrolysis

Understanding phosphoryl and sulfuryl transfer is central to many biochemical processes. However, despite decades of experimental and computational studies, a consensus concerning the precise mechanistic details of these reactions has yet to be reached. In this work we perform a detailed comparative theoretical study of the hydrolysis of p-nitrophenyl phosphate, methyl phosphate and p-nitrophenyl sulfate, all of which have served as key model systems for understanding phosphoryl and sulfuryl transfer reactions, respectively. We demonstrate the existence of energetically similar but mechanistically distinct possibilities for phosphate monoester hydrolysis. The calculated kinetic isotope effects for p-nitrophenyl phosphate provide a means to discriminate between substrate- and solvent-assisted pathways of phosphate monoester hydrolysis, and show that the solvent-assisted pathway dominates in solution. This preferred mechanism for p-nitrophenyl phosphate hydrolysis is difficult to find computationally due to the limitations of compressing multiple bonding changes onto a 2-dimensional energy surface. This problem is compounded by the need to include implicit solvation to at least microsolvate the system and stabilize the highly charged species. In contrast, methyl phosphate hydrolysis shows a preference for a substrate-assisted mechanism. For p-nitrophenyl sulfate hydrolysis there is only one viable reaction pathway, which is similar to the solvent-assisted pathway for phosphate hydrolysis, and the substrate-assisted pathway is not accessible. Overall, our results provide a unifying mechanistic framework that is consistent with the experimentally measured kinetic isotope effects and reconciles the discrepancies between theoretical and experimental models for these biochemically ubiquitous classes of reaction.

Ribosome-induced tuning of GTP hydrolysis by a translational GTPase

GTP hydrolysis by elongation factor Tu (EF-Tu), a translational GTPase that delivers aminoacyl-tRNAs to the ribosome, plays a crucial role in decoding and translational fidelity. The basic reaction mechanism and the way the ribosome contributes to catalysis are a matter of debate. Here we use mutational analysis in combination with measurements of rate/pH profiles, kinetic solvent isotope effects, and ion dependence of GTP hydrolysis by EF-Tu off and on the ribosome to dissect the reaction mechanism. Our data suggest that--contrary to current models--the reaction in free EF-Tu follows a pathway that does not involve the critical residue H84 in the switch II region. Binding to the ribosome without a cognate codon in the A site has little effect on the GTPase mechanism. In contrast, upon cognate codon recognition, the ribosome induces a rearrangement of EF-Tu that renders GTP hydrolysis sensitive to mutations of Asp21 and His84 and insensitive to K(+) ions. We suggest that Asp21 and His84 provide a network of interactions that stabilize the positions of the γ-phosphate and the nucleophilic water, respectively, and thus play an indirect catalytic role in the GTPase mechanism on the ribosome.

A monovalent cation acts as structural and catalytic cofactor in translational GTPases

Translational GTPases are universally conserved GTP hydrolyzing enzymes, critical for fidelity and speed of ribosomal protein biosynthesis. Despite their central roles, the mechanisms of GTP-dependent conformational switching and GTP hydrolysis that govern the function of trGTPases remain poorly understood. Here, we provide biochemical and high-resolution structural evidence that eIF5B and aEF1A/EF-Tu bound to GTP or GTPγS coordinate a monovalent cation (M(+)) in their active site. Our data reveal that M(+) ions form constitutive components of the catalytic machinery in trGTPases acting as structural cofactor to stabilize the GTP-bound "on" state. Additionally, the M(+) ion provides a positive charge into the active site analogous to the arginine-finger in the Ras-RasGAP system indicating a similar role as catalytic element that stabilizes the transition state of the hydrolysis reaction. In sequence and structure, the coordination shell for the M(+) ion is, with exception of eIF2γ, highly conserved among trGTPases from bacteria to human. We therefore propose a universal mechanism of M(+)-dependent conformational switching and GTP hydrolysis among trGTPases with important consequences for the interpretation of available biochemical and structural data.

Quantitative exploration of the molecular origin of the activation of GTPase

GTPases play a major role in cellular processes, and gaining quantitative understanding of their activation demands reliable free energy surfaces of the relevant mechanistic paths in solution, as well as the interpolation of this information to GTPases. Recently, we generated ab initio quantum mechanical/molecular mechanical free energy surfaces for the hydrolysis of phosphate monoesters in solution, establishing quantitatively that the barrier for the reactions with a proton transfer (PT) step from a single attacking water (1 W) is higher than the one where the PT is assisted by a second water (2 W). The implication of this finding on the activation of GTPases is quantified here, by using the ab initio solution surfaces to calibrate empirical valence bond surfaces and then exploring the origin of the activation effect. It is found that, although the 2 W PT path is a new element, this step is not rate determining, and the catalytic effect is actually due to the electrostatic stabilization of the pre-PT transition state and the subsequent plateau. Thus, the electrostatic catalytic effect found in our previous studies of the Ras GTPase activating protein (RasGAP) and the elongation factor-Tu (EF-Tu) with a 1 W mechanism is still valid for the 2 W path. Furthermore, as found before, the corresponding activation appears to involve a major allosteric effect. Overall, we believe that our finding is general to both GTPases and ATPases. In addition to the biologically relevant finding, we also provide a critical discussion of the requirements from reliable surfaces for enzymatic reactions.

How conformational flexibility stabilizes the hyperthermophilic elongation factor G-domain

Proteins from thermophilic organisms are stable and functional well above ambient temperature. Understanding the molecular mechanism underlying such a resistance is of crucial interest for many technological applications. For some time, thermal stability has been assumed to correlate with high mechanical rigidity of the protein matrix. In this work we address this common belief by carefully studying a pair of homologous G-domain proteins, with their melting temperatures differing by 40 K. To probe the thermal-stability content of the two proteins we use extensive simulations covering the microsecond time range and employ several different indicators to assess the salient features of the conformational landscape and the role of internal fluctuations at ambient condition. At the atomistic level, while the magnitude of fluctuations is comparable, the distribution of flexible and rigid stretches of amino-acids is more regular in the thermophilic protein causing a cage-like correlation of amplitudes along the sequence. This caging effect is suggested to favor stability at high T by confining the mechanical excitations. Moreover, it is found that the thermophilic protein, when folded, visits a higher number of conformational substates than the mesophilic homologue. The entropy associated with the occupation of the different substates and the thermal resilience of the protein intrinsic compressibility provide a qualitative insight on the thermal stability of the thermophilic protein as compared to its mesophilic homologue. Our findings potentially open the route to new strategies in the design of thermostable proteins.

Energetics of activation of GTP hydrolysis on the ribosome

Several of the steps in protein synthesis on the ribosome utilize hydrolysis of guanosine triphosphate (GTP) as the driving force. This reaction is catalyzed by translation factors that become activated upon binding to the ribosome. The recently determined crystal structure of an elongation factor-Tu ternary complex bound to the ribosome allows the energetics of GTP activation to be explored by computer simulations. A central problem regards the role of the universally conserved histidine, which has been proposed to act as a general base for guanosine triphosphate hydrolysis. Here we report a detailed energetic and structural analysis of different possible protonation states that could be involved in activation of the reaction. We show that the histidine cannot act as a general base, but must be protonated and in its active conformation to promote GTP hydrolysis. We further show that the sarcin-ricin loop of the ribosome spontaneously drives the histidine into the correct conformation for GTP activation.

Mechanism of activation of elongation factor Tu by ribosome: catalytic histidine activates GTP by protonation

Elongation factor Tu (EF-Tu) is central to prokaryotic protein synthesis as it has the role of delivering amino-acylated tRNAs to the ribosome. Release of EF-Tu, after correct binding of the EF-Tu:aa-tRNA complex to the ribosome, is initiated by GTP hydrolysis. This reaction, whose mechanism is uncertain, is catalyzed by EF-Tu, but requires activation by the ribosome. There have been a number of mechanistic proposals, including those spurred by a recent X-ray crystallographic analysis of a ribosome:EF-Tu:aa-tRNA:GTP-analog complex. In this work, we have investigated these and alternative hypotheses, using high-level quantum chemical/molecular mechanical simulations for the wild-type protein and its His85Gln mutant. For both proteins, we find previously unsuggested mechanisms as being preferred, in which residue 85, either His or Gln, directly assists in the reaction. Analysis shows that the RNA has a minor catalytic effect in the wild-type reaction, but plays a significant role in the mutant by greatly stabilizing the reaction's transition state. Given the similarity between EF-Tu and other members of the translational G-protein family, it is likely that these mechanisms of ribosome-activated GTP hydrolysis are pertinent to all of these proteins.

Elongation factor G bound to the ribosome in an intermediate state of translocation

A key step of translation by the ribosome is translocation, which involves the movement of messenger RNA (mRNA) and transfer RNA (tRNA) with respect to the ribosome. This allows a new round of protein chain elongation by placing the next mRNA codon in the A site of the 30S subunit. Translocation proceeds through an intermediate state in which the acceptor ends of the tRNAs have moved with respect to the 50S subunit but not the 30S subunit, to form hybrid states. The guanosine triphosphatase (GTPase) elongation factor G (EF-G) catalyzes the subsequent movement of mRNA and tRNA with respect to the 30S subunit. Here, we present a crystal structure at 3 angstrom resolution of the Thermus thermophilus ribosome with a tRNA in the hybrid P/E state bound to EF-G with a GTP analog. The structure provides insights into structural changes that facilitate translocation and suggests a common GTPase mechanism for EF-G and elongation factor Tu.

Quantifying the mechanism of phosphate monoester hydrolysis in aqueous solution by evaluating the relevant ab initio QM/MM free-energy surfaces

Understanding the nature of the free-energy surfaces for phosphate hydrolysis is a prerequisite for understanding the corresponding key chemical reactions in biology. Here, the challenge has been to move to careful ab initio QM/MM (QM(ai)/MM) free-energy calculations, where obtaining converging results is very demanding and computationally expensive. This work describes such calculations, focusing on the free-energy surface for the hydrolysis of phosphate monoesters, paying special attention to the comparison between the one water (1W) and two water (2W) paths for the proton-transfer (PT) step. This issue has been explored before by energy minimization with implicit solvent models and by nonsystematic QM/MM energy minimization, as well as by nonsystematic free-energy mapping. However, no study has provided the needed reliable 2D (3D) surfaces that are necessary for reaching concrete conclusions. Here we report a systematic evaluation of the 2D (3D) free-energy maps for several relevant systems, comparing the results of QM(ai)/MM and QM(ai)/implicit solvent surfaces, and provide an advanced description of the relevant energetics. It is found that the 1W path for the hydrolysis of the methyl diphosphate (MDP) trianion is 6-9 kcal/mol higher than that the 2W path. This difference becomes slightly larger in the presence of the Mg(2+) ion because this ion reduces the pKa of the conjugated acid form of the phosphate oxygen that accepts the proton. Interestingly, the BLYP approach (which has been used extensively in some studies) gives a much smaller difference between the 1W and 2W activation barriers. At any rate, it is worth pointing out that the 2W transition state for the PT is not much higher that the common plateau that serves as the starting point of both the 1W and 2W PT paths. Thus, the calculated catalytic effects of proteins based on the 2W PT mechanistic model are not expected to be different from the catalytic effects predicted using the 1W PT mechanistic model, which was calibrated on the observed barrier in solution and in which the TS charge distribution was similar to the that of the plateau (as was done in all of our previous EVB studies).