Modeling the mechanisms of biological GTP hydrolysis

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.

Unveiling the State Transition Mechanisms of Ras Proteins through Enhanced Sampling and QM/MM Simulations

In cells, wild-type RasGTP complexes exist in two distinct states: active State 2 and inactive State 1. These complexes regulate their functions by transitioning between the two states. However, the mechanisms underlying this state transition have not been clearly elucidated. To address this, we conducted a detailed simulation study to characterize the energetics of the stable states involved in the state transitions of the HRasGTP complex, specifically from State 2 to State 1. This was achieved by employing multiscale quantum mechanics/molecular mechanics and enhanced sampling molecular dynamics methods. Based on the simulation results, we constructed the two-dimensional free energy landscapes that provide crucial information about the conformational changes of the HRasGTP complex from State 2 to State 1. Furthermore, we also explored the conformational changes from the intermediate state to the product state during guanosine triphosphate hydrolysis. This study on the conformational changes involved in the HRas state transitions serves as a valuable reference for understanding the corresponding events of both KRas and NRas as well.

Far-reaching effects of tyrosine64 phosphorylation on Ras revealed with BeF3- complexes

Tyrosine phosphorylation on Ras by Src kinase is known to uncouple Ras from upstream regulation and downstream communication. However, the mechanisms by which phosphorylation modulates these interactions have not been detailed. Here, the major mono-phosphorylation level on tyrosine64 is quantified by 31P NMR and mutagenesis. Crystal structures of unphosphorylated and tyrosine64-phosphorylated Ras in complex with a BeF3- ground state analogue reveal "closed" Ras conformations very different from those of the "open" conformations previously observed for non-hydrolysable GTP analogue structures of Ras. They deliver new mechanistic and conformational insights into intrinsic GTP hydrolysis. Phosphorylation of tyrosine64 delivers conformational changes distant from the active site, showing why phosphorylated Ras has reduced affinity to its downstream effector Raf. 19F NMR provides evidence for changes in the intrinsic GTPase and nucleotide exchange rate and identifies the concurrent presence of a major "closed" conformation alongside a minor yet functionally important "open" conformation at the ground state of Ras. This study expands the application of metal fluoride complexes in revealing major and minor conformational changes of dynamic and modified Ras proteins.

Comparative genomic analysis of the RabGAP gene family in seven Rosaceae species, and functional identification of PbrRabGAP10 in controlling pollen tube growth by mediating cellulose deposition in pear

Rab GTPase-activating proteins (RabGAPs), serving as crucial signaling switches, play essential roles in several physiological processes related to plant growth and development. However, despite their importance, information regarding the RabGAP gene family and their biological functions remains unknown in the Rosaceae. In this study, we identified a total of 127 RabGAP genes in seven Rosaceae species, which were divided into five subfamilies. Our findings indicate that whole genome duplication (WGD) events or dispersed duplication events largely contributed to the expansion of RabGAP family members within Rosaceae species. Through tissue-specific expression analyses, we revealed that the PbrRabGAP genes exhibited distinct expression patterns in different pear tissues. Furthermore, by examining the expression pattern during pollen development and employing an antisense oligonucleotide approach, we demonstrated that PbrRabGAP10, located in the cytoplasm, mediates the imbalance of cellulose distribution, thus regulating pollen tube elongation. In conclusion, the present study offers an overview of the RabGAP family in Rosaceae genomes and serves as the basis for further functional studies.

Mechanism-Based Redesign of GAP to Activate Oncogenic Ras

Ras GTPases play a crucial role in cell signaling pathways. Mutations of the Ras gene occur in about one third of cancerous cell lines and are often associated with detrimental clinical prognosis. Hot spot residues Gly12, Gly13, and Gln61 cover 97% of oncogenic mutations, which impair the enzymatic activity in Ras. Using QM/MM free energy calculations, we present a two-step mechanism for the GTP hydrolysis catalyzed by the wild-type Ras.GAP complex. We found that the deprotonation of the catalytic water takes place via the Gln61 as a transient Brønsted base. We also determined the reaction profiles for key oncogenic Ras mutants G12D and G12C using QM/MM minimizations, matching the experimentally observed loss of catalytic activity, thereby validating our reaction mechanism. Using the optimized reaction paths, we devised a fast and accurate procedure to design GAP mutants that activate G12D Ras. We replaced GAP residues near the active site and determined the activation barrier for 190 single mutants. We furthermore built a machine learning for ultrafast screening, by fast prediction of the barrier heights, tested both on the single and double mutations. This work demonstrates that fast and accurate screening can be accomplished via QM/MM reaction path optimizations to design protein sequences with increased catalytic activity. Several GAP mutations are predicted to re-enable catalysis in oncogenic G12D, offering a promising avenue to overcome aberrant Ras-driven signal transduction by activating enzymatic activity instead of inhibition. The outlined computational screening protocol is readily applicable for designing ligands and cofactors analogously.

Bacterial GTPases as druggable targets to tackle antimicrobial resistance

Small molecules as antibacterial agents have contributed immensely to the growth of modern medicine over the last several decades. However, the emergence of drug resistance among bacterial pathogens has undermined the effectiveness of the existing antibiotics. Thus, there is an exigency to address the antibiotic crisis by developing new antibacterial agents and identifying novel drug targets in bacteria. In this review, we summarize the importance of guanosine triphosphate hydrolyzing proteins (GTPases) as key agents for bacterial survival. We also discuss representative examples of small molecules that target bacterial GTPases as novel antibacterial agents, and highlight areas that are ripe for exploration. Given their vital roles in cell viability, virulence, and antibiotic resistance, bacterial GTPases are highly attractive antibacterial targets that will likely play a vital role in the fight against antimicrobial resistance.

Transcriptional and Epigenetic Alterations in the Progression of Non-Alcoholic Fatty Liver Disease and Biomarkers Helping to Diagnose Non-Alcoholic Steatohepatitis

Non-alcoholic fatty liver disease (NAFLD) encompasses a broad spectrum of conditions from simple steatosis (non-alcoholic fatty liver (NAFL)) to non-alcoholic steatohepatitis (NASH), and its global prevalence continues to rise. NASH, the progressive form of NAFLD, has higher risks of liver and non-liver related adverse outcomes compared with those patients with NAFL alone. Therefore, the present study aimed to explore the mechanisms in the progression of NAFLD and to develop a model to diagnose NASH based on the transcriptome and epigenome. Differentially expressed genes (DEGs) and differentially methylated genes (DMGs) among the three groups (normal, NAFL, and NASH) were identified, and the functional analysis revealed that the development of NAFLD was primarily related to the oxidoreductase-related activity, PPAR signaling pathway, tight junction, and pathogenic Escherichia coli infection. The logistic regression (LR) model, consisting of ApoF, THOP1, and BICC1, outperformed the other five models. With the highest AUC (0.8819, 95%CI: 0.8128-0.9511) and a sensitivity of 97.87%, as well as a specificity of 64.71%, the LR model was determined as the diagnostic model, which can differentiate NASH from NAFL. In conclusion, several potential mechanisms were screened out based on the transcriptome and epigenome, and a diagnostic model was built to help patient stratification for NAFLD populations.

A biochemical view on the septins, a less known component of the cytoskeleton

The septins are a conserved family of guanine nucleotide binding proteins, often named the fourth component of the cytoskeleton. They self-assemble into non-polar filaments and further into higher ordered structures. Properly assembled septin structures are required for a wide range of indispensable intracellular processes such as cytokinesis, vesicular transport, polarity establishment and cellular adhesion. Septins belong structurally to the P-Loop NTPases. However, unlike the small GTPases like Ras, septins do not mediate signals to effectors through GTP binding and hydrolysis. The role of nucleotide binding and subsequent GTP hydrolysis by the septins is rather controversially debated. We compile here the structural features from the existing septin crystal- and cryo-EM structures regarding protofilament formation, inter-subunit interface architecture and nucleotide binding and hydrolysis. These findings are supplemented with a summary of available biochemical studies providing information regarding nucleotide binding and hydrolysis of fungal and mammalian septins.

The role of nucleoside triphosphate hydrolase metallochaperones in making metalloenzymes

Metalloenzymes catalyze a diverse set of challenging chemical reactions that are essential for life. These metalloenzymes rely on a wide range of metallocofactors, from single metal ions to complicated metallic clusters. Incorporation of metal ions and metallocofactors into apo-proteins often requires the assistance of proteins known as metallochaperones. Nucleoside triphosphate hydrolases (NTPases) are one important class of metallochaperones and are found widely distributed throughout the domains of life. These proteins use the binding and hydrolysis of nucleoside triphosphates, either adenosine triphosphate (ATP) or guanosine triphosphate (GTP), to carry out highly specific and regulated roles in the process of metalloenzyme maturation. Here, we review recent literature on NTPase metallochaperones and describe the current mechanistic proposals and available structural data. By using representative examples from each type of NTPase, we also illustrate the challenges in studying these complicated systems. We highlight open questions in the field and suggest future directions. This minireview is part of a special collection of articles in memory of Professor Deborah Zamble, a leader in the field of nickel biochemistry.

Mechanism of Guanosine Triphosphate Hydrolysis by the Visual Proteins Arl3-RP2: Free Energy Reaction Profiles Computed with Ab Initio Type QM/MM Potentials

We report the results of calculations of the Gibbs energy profiles of the guanosine triphosphate (GTP) hydrolysis by the Arl3-RP2 protein complex using molecular dynamics (MD) simulations with ab initio type QM/MM potentials. The chemical reaction of GTP hydrolysis to guanosine diphosphate (GDP) and inorganic phosphate (Pi) is catalyzed by GTPases, the enzymes, which are responsible for signal transduction in live cells. A small GTPase Arl3, catalyzing the GTP → GDP reaction in complex with the activating protein RP2, constitute an essential part of the human vision cycle. To simulate the reaction mechanism, a model system is constructed by motifs of the crystal structure of the Arl3-RP2 complexed with a substrate analog. After selection of reaction coordinates, energy profiles for elementary steps along the reaction pathway GTP + H2O → GDP + Pi are computed using the umbrella sampling and umbrella integration procedures. QM/MM MD calculations are carried out, interfacing the molecular dynamics program NAMD and the quantum chemistry program TeraChem. Ab initio type QM(DFT)/MM potentials are computed with atom-centered basis sets 6-31G** and two hybrid functionals (PBE0-D3 and ωB97x-D3) of the density functional theory, describing a large QM subsystem. Results of these simulations of the reaction mechanism are compared to those obtained with QM/MM calculations on the potential energy surface using a similar description of the QM part. We find that both approaches, QM/MM and QM/MM MD, support the mechanism of GTP hydrolysis by GTPases, according to which the catalytic glutamine side chain (Gln71, in this system) actively participates in the reaction. Both approaches distinguish two parts of the reaction: the cleavage of the phosphorus-oxygen bond in GTP coupled with the formation of Pi, and the enzyme regeneration. Newly performed QM/MM MD simulations confirmed the profile predicted in the QM/MM minimum energy calculations, called here the pathway-I, and corrected its relief at the first elementary step from the enzyme-substrate complex. The QM/MM MD simulations also revealed another mechanism at the part of enzyme regeneration leading to pathway-II. Pathway-II is more consistent with the experimental kinetic data of the wild-type complex Arl3-RP2, whereas pathway-I explains the role of the mutation Glu138Gly in RP2 slowing down the hydrolysis rate.

Unveiling the Dynamics of KRAS4b on Lipid Model Membranes

Small GTPase proteins are ubiquitous and responsible for regulating several processes related to cell growth and differentiation. Mutations that stabilize their active state can lead to uncontrolled cell proliferation and cancer. Although these proteins are well characterized at the cellular scale, the molecular mechanisms governing their functions are still poorly understood. In addition, there is limited information about the regulatory function of the cell membrane which supports their activity. Thus, we have studied the dynamics and conformations of the farnesylated KRAS4b in various membrane model systems, ranging from binary fluid mixtures to heterogeneous raft mimics. Our approach combines long time-scale coarse-grained (CG) simulations and Markov state models to dissect the membrane-supported dynamics of KRAS4b. Our simulations reveal that protein dynamics is mainly modulated by the presence of anionic lipids and to some extent by the nucleotide state (activation) of the protein. In addition, our results suggest that both the farnesyl and the polybasic hypervariable region (HVR) are responsible for its preferential partitioning within the liquid-disordered (Ld) domains in membranes, potentially enhancing the formation of membrane-driven signaling platforms.

Recent Advances in Understanding Biological GTP Hydrolysis through Molecular Simulation

GTP hydrolysis is central to biology, being involved in regulating a wide range of cellular processes. However, the mechanisms by which GTPases hydrolyze this critical reaction remain controversial, with multiple mechanistic possibilities having been proposed based on analysis of experimental and computational data. In this mini-review, we discuss advances in our understanding of biological GTP hydrolysis based on recent computational studies and argue in favor of solvent-assisted hydrolysis as a conserved mechanism among GTPases. A concrete understanding of the fundamental mechanisms by which these enzymes facilitate GTP hydrolysis will have significant impact both for drug discovery efforts and for unraveling the role of oncogenic mutations.

Diversity of mechanisms in Ras-GAP catalysis of guanosine triphosphate hydrolysis revealed by molecular modeling

The mechanism of the deceptively simple reaction of guanosine triphosphate (GTP) hydrolysis catalyzed by the cellular protein Ras in complex with the activating protein GAP is an important issue because of the significance of this reaction in cancer research. We show that molecular modeling of GTP hydrolysis in the Ras-GAP active site reveals a diversity of mechanisms of the intrinsic chemical reaction depending on molecular groups at position 61 in Ras occupied by glutamine in the wild-type enzyme. First, a comparison of reaction energy profiles computed at the quantum mechanics/molecular mechanics (QM/MM) level shows that an assignment of the Gln61 side chain in the wild-type Ras either to QM or to MM parts leads to different scenarios corresponding to the glutamine-assisted or the substrate-assisted mechanisms. Second, replacement of Gln61 by the nitro-analog of glutamine (NGln) or by Glu, applied in experimental studies, results in two more scenarios featuring the so-called two-water and the concerted-type mechanisms. The glutamine-assisted mechanism in the wild-type Ras-GAP, in which the conserved Gln61 plays a decisive role, switching between the amide and imide tautomer forms, is consistent with the known experimental results of structural, kinetic and spectroscopy studies. The results emphasize the role of the Ras residue Gln61 in Ras-GAP catalysis and explain the retained catalytic activity of the Ras-GAP complex towards GTP hydrolysis in the Gln61NGln and Gln61Glu mutants of Ras.

The Impact of the Stringent Response on TRAFAC GTPases and Prokaryotic Ribosome Assembly

Many facets of ribosome biogenesis and function, including ribosomal RNA (rRNA) transcription, 70S assembly and protein translation, are negatively impacted upon induction of a nutrient stress-sensing signalling pathway termed the stringent response. This stress response is mediated by the alarmones guanosine tetra- and penta-phosphate ((p)ppGpp), the accumulation of which leads to a massive cellular response that slows growth and aids survival. The 70S bacterial ribosome is an intricate structure, with assembly both complex and highly modular. Presiding over the assembly process is a group of P-loop GTPases within the TRAFAC (Translation Factor Association) superclass that are crucial for correct positioning of both early and late stage ribosomal proteins (r-proteins) onto the rRNA. Often described as 'molecular switches', members of this GTPase superfamily readily bind and hydrolyse GTP to GDP in a cyclic manner that alters the propensity of the GTPase to carry out a function. TRAFAC GTPases are considered to act as checkpoints to ribosome assembly, involved in binding to immature sections in the GTP-bound state, preventing further r-protein association until maturation is complete. Here we review our current understanding of the impact of the stringent response and (p)ppGpp production on ribosome maturation in prokaryotic cells, focusing on the inhibition of (p)ppGpp on GTPase-mediated subunit assembly, but also touching upon the inhibition of rRNA transcription and protein translation.

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.

Hybrid QM/MM vs Pure MM Molecular Dynamics for Evaluating Water Distribution within p21N-ras and the Resulting GTP Electronic Density

p21^ras protein activity, regulated by GTP hydrolysis, constitutes an active field of research for the development of cancer targeted therapies that would concern ∼30% of human tumors to which specific mutations have been associated. Indeed, the catalyzing mechanisms provided by the protein environment during GTP hydrolysis and how they are impaired by specific mutations remain to be fully elucidated. In this article, we present results from molecular mechanics (MM) molecular dynamics (MD) simulations and density functional theory (DFT) calculations carried out for wild-type p21 ^N-ras and six Gln 61 mutants. In the first part, we present the water distribution within the active site of the wild-type protein according to MM MD. Significant differences are observed when comparing the results to the previous distribution assessed through quantum mechanics/molecular mechanics (QM/MM) MD. Such method-dependent results highlight the importance of accounting for the electrostatic coupling between the protein complex and the solvent molecules in identifying hydration sites. In the second part, we present the results from DFT calculations performed to determine the electronic distribution of the GTP ligand, considering the wild-type active site arrangement according to both classical and hybrid approaches. Only in the QM/MM-based configuration is the ligand electronic density similar to that of a GDP-like state observed experimentally. For this reason, in the last set of calculations carried out for p21 ^N-ras Gln 61 mutants, only the active site structural conformations obtained through hybrid MD are considered. Through the analysis of the GTP electronic density, we conclude that the wild-type active site arrangement according to QM/MM MD is closer to a catalytically efficient conformation of the protein than the arrangement according to MM MD. Hence, water distribution according to the hybrid approach must correspond to the optimal placement of solvent in the active site. Within all of the studied Gln 61 substituted proteins, p21^ras major catalyzing effect, which consists of stabilizing a more GDP-like state, is lost.

A unique dynamin-related protein is essential for mitochondrial fission in Toxoplasma gondii

The single mitochondrion of apicomplexan protozoa is thought to be critical for all stages of the life cycle, and is a validated drug target against these important human and veterinary parasites. In contrast to other eukaryotes, replication of the mitochondrion is tightly linked to the cell cycle. A key step in mitochondrial segregation is the fission event, which in many eukaryotes occurs by the action of dynamins constricting the outer membrane of the mitochondria from the cytosolic face. To date, none of the components of the apicomplexan fission machinery have been identified and validated. We identify here a highly divergent, dynamin-related protein (TgDrpC), conserved in apicomplexans as essential for mitochondrial biogenesis and potentially for fission in Toxoplasma gondii. We show that TgDrpC is found adjacent to the mitochondrion, and is localised both at its periphery and at its basal part, where fission is expected to occur. We demonstrate that depletion or dominant negative expression of TgDrpC results in interconnected mitochondria and ultimately in drastic changes in mitochondrial morphology, as well as in parasite death. Intriguingly, we find that the canonical adaptor TgFis1 is not required for mitochondrial fission. The identification of an Apicomplexa-specific enzyme required for mitochondrial biogenesis and essential for parasite growth highlights parasite adaptation. This work paves the way for future drug development targeting TgDrpC, and for the analysis of additional partners involved in this crucial step of apicomplexan multiplication.

Exposing catalytic versatility of GTPases: taking reaction detours in mutants of hGBP1 enzyme without additional energetic cost

Our study reveals that the replacement of catalytically competent residues by the inert amino acid alanine, S73A and E99A, in hGBP1 opens a plethora of molecularly different reaction pathways featuring very similar energy barriers as the wild type.

Quantitative Measurement of Intrinsic GTP Hydrolysis for Carcinogenic Glutamine 61 Mutants in H-Ras

Mutations of human oncoprotein p²¹H-Ras (hereafter "Ras") at glutamine 61 are known to slow the rate of guanosine triphosphate (GTP) hydrolysis and transform healthy cells into malignant cells. It has been hypothesized that this glutamine plays a role in the intrinsic mechanism of GTP hydrolysis by interacting with an active site water molecule that stabilizes the formation of the charged transition state at the γ-phosphate during hydrolysis. However, there is no comprehensive data set of the effects of mutations to Q61 on the protein's intrinsic catalytic rate, structure, or interactions with water at the active site. Here, we present the first comprehensive and quantitative set of initial rates of intrinsic hydrolysis for all stable variants of RasQ61X. We further conducted enhanced molecular dynamics (MD) simulations of each construct to determine the solvent accessible surface area (SASA) of the side chain at position 61 and compared these results to previously measured changes in electric fields caused by RasQ61X mutations. For polar and negatively charged residues, we found that the rates are normally distributed about an optimal electrostatic contribution, close to that of the native Q61 residue, and the rates are strongly correlated to the number of waters in the active site. Together, these results support a mechanism of GTP hydrolysis in which Q61 stabilizes a transient hydronium ion, which then stabilizes the transition state while the γ-phosphate is undergoing nucleophilic attack by a second, catalytically active water molecule. We discuss the implications of such a mechanism on future strategies for combating Ras-based cancers.

Molecular Mechanism of ATP Hydrolysis in an ABC Transporter

Hydrolysis of nucleoside triphosphate (NTP) plays a key role for the function of many biomolecular systems. However, the chemistry of the catalytic reaction in terms of an atomic-level understanding of the structural, dynamic, and free energy changes associated with it often remains unknown. Here, we report the molecular mechanism of adenosine triphosphate (ATP) hydrolysis in the ATP-binding cassette (ABC) transporter BtuCD-F. Free energy profiles obtained from hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations show that the hydrolysis reaction proceeds in a stepwise manner. First, nucleophilic attack of an activated lytic water molecule at the ATP γ-phosphate yields ADP + HPO₄ ^2- as intermediate product. A conserved glutamate that is located very close to the γ-phosphate transiently accepts a proton and thus acts as catalytic base. In the second step, the proton is transferred back from the catalytic base to the γ-phosphate, yielding ADP + H₂PO₄ ^-. These two chemical reaction steps are followed by rearrangements of the hydrogen bond network and the coordination of the Mg²⁺ ion. The rate constant estimated from the computed free energy barriers is in very good agreement with experiments. The overall free energy change of the reaction is close to zero, suggesting that phosphate bond cleavage itself does not provide a power stroke for conformational changes. Instead, ATP binding is essential for tight dimerization of the nucleotide-binding domains and the transition of the transmembrane domains from inward- to outward-facing, whereas ATP hydrolysis resets the conformational cycle. The mechanism is likely relevant for all ABC transporters and might have implications also for other NTPases, as many residues involved in nucleotide binding and hydrolysis are strictly conserved.

Amide-imide tautomerization in the glutamine side chain in enzymatic and photochemical reactions in proteins

Amide-imide tautomerization presents a pervasive class of chemical transformations in organic chemistry of natural compounds. In this Perspective, we describe two distinctively different protein systems, in which the amide-imide tautomerization in the glutamine side chain takes place in enzymatic or photochemical reactions. First, hydrolysis of guanosine triphosphate (GTP) catalyzed by the Ras-GAP protein complex suggests the occurrence of the imide tautomer of glutamine in reaction intermediates. Second, photoexcitation of flavin-binding protein domains (BLUFs) initiates a chain of reactions in the chromophore-binding pocket, including amide-imide tautomerization of glutamine. Mechanisms of these reactions at the atomic level have been revealed in quantum mechanics/molecular mechanics (QM/MM) simulations. To reinforce conclusions on the critical role of amide-imide tautomerization of glutamine in these reactions we describe results of new quantum chemistry and QM/MM calculations for relevant molecular model systems. We reexamine results of the recent IR spectroscopy studies of BLUF domains, which provide experimental evidences of Gln tautomerization in proteins. We also propose to validate the glutamine-assisted mechanism of enzymatic GTP hydrolysis by using IR spectroscopy in a proper range of wavenumbers.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Challenges and advances in the computational modeling of biological phosphate hydrolysis

Phosphate ester hydrolysis is fundamental to many life processes, and has been the topic of substantial experimental and computational research effort. However, even the simplest of phosphate esters can be hydrolyzed through multiple possible pathways that can be difficult to distinguish between, either experimentally, or computationally. Therefore, the mechanisms of both the enzymatic and non-enzymatic reactions have been historically controversial. In the present contribution, we highlight a number of technical issues involved in reliably modeling these computationally challenging reactions, as well as proposing potential solutions. We also showcase examples of our own work in this area, discussing both the non-enzymatic reaction in aqueous solution, as well as insights obtained from the computational modeling of organophosphate hydrolysis and catalytic promiscuity amongst enzymes that catalyze phosphoryl transfer.

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.

Assessing the Influence of Mutation on GTPase Transition States by Using X-ray Crystallography, 19 F NMR, and DFT Approaches

We report X-ray crystallographic and ¹⁹ F NMR studies of the G-protein RhoA complexed with MgF₃^- , GDP, and RhoGAP, which has the mutation Arg85'Ala. When combined with DFT calculations, these data permit the identification of changes in transition state (TS) properties. The X-ray data show how Tyr34 maintains solvent exclusion and the core H-bond network in the active site by relocating to replace the missing Arg85' sidechain. The ¹⁹ F NMR data show deshielding effects that indicate the main function of Arg85' is electronic polarization of the transferring phosphoryl group, primarily mediated by H-bonding to O^3G and thence to P^G . DFT calculations identify electron-density redistribution and pinpoint why the TS for guanosine 5'-triphosphate (GTP) hydrolysis is higher in energy when RhoA is complexed with RhoGAP_Arg85'Ala relative to wild-type (WT) RhoGAP. This study demonstrates that ¹⁹ F NMR measurements, in combination with X-ray crystallography and DFT calculations, can reliably dissect the response of small GTPases to site-specific modifications.

Common mechanisms of catalysis in small and heterotrimeric GTPases and their respective GAPs

GTPases are central switches in cells. Their dysfunctions are involved in severe diseases. The small GTPase Ras regulates cell growth, differentiation and apoptosis by transmitting external signals to the nucleus. In one group of oncogenic mutations, the 'switch-off' reaction is inhibited, leading to persistent activation of the signaling pathway. The switch reaction is regulated by GTPase-activating proteins (GAPs), which catalyze GTP hydrolysis in Ras, and by guanine nucleotide exchange factors, which catalyze the exchange of GDP for GTP. Heterotrimeric G-proteins are activated by G-protein coupled receptors and are inactivated by GTP hydrolysis in the Gα subunit. Their GAPs are called regulators of G-protein signaling. In the same way that Ras serves as a prototype for small GTPases, Gαi1 is the most well-studied Gα subunit. By utilizing X-ray structural models, time-resolved infrared-difference spectroscopy, and biomolecular simulations, we elucidated the detailed molecular reaction mechanism of the GTP hydrolysis in Ras and Gαi1. In both proteins, the charge distribution of GTP is driven towards the transition state, and an arginine is precisely positioned to facilitate nucleophilic attack of water. In addition to these mechanistic details of GTP hydrolysis, Ras dimerization as an emerging factor in signal transduction is discussed in this review.

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.

The advantage of channeling nucleotides for very processive functions

Nucleoside triphosphate (NTP)s, like ATP (adenosine 5'-triphosphate) and GTP (guanosine 5'-triphosphate), have long been considered sufficiently concentrated and diffusible to fuel all cellular ATPases (adenosine triphosphatases) and GTPases (guanosine triphosphatases) in an energetically healthy cell without becoming limiting for function. However, increasing evidence for the importance of local ATP and GTP pools, synthesised in close proximity to ATP- or GTP-consuming reactions, has fundamentally challenged our view of energy metabolism. It has become evident that cellular energy metabolism occurs in many specialised 'microcompartments', where energy in the form of NTPs is transferred preferentially from NTP-generating modules directly to NTP-consuming modules. Such energy channeling occurs when diffusion through the cytosol is limited, where these modules are physically close and, in particular, if the NTP-consuming reaction has a very high turnover, i.e. is very processive. Here, we summarise the evidence for these conclusions and describe new insights into the physiological importance and molecular mechanisms of energy channeling gained from recent studies. In particular, we describe the role of glycolytic enzymes for axonal vesicle transport and nucleoside diphosphate kinases for the functions of dynamins and dynamin-related GTPases.

The advantage of channeling nucleotides for very processive functions

Nucleoside triphosphate (NTP)s, like ATP (adenosine 5'-triphosphate) and GTP (guanosine 5'-triphosphate), have long been considered sufficiently concentrated and diffusible to fuel all cellular ATPases (adenosine triphosphatases) and GTPases (guanosine triphosphatases) in an energetically healthy cell without becoming limiting for function. However, increasing evidence for the importance of local ATP and GTP pools, synthesised in close proximity to ATP- or GTP-consuming reactions, has fundamentally challenged our view of energy metabolism. It has become evident that cellular energy metabolism occurs in many specialised 'microcompartments', where energy in the form of NTPs is transferred preferentially from NTP-generating modules directly to NTP-consuming modules. Such energy channeling occurs when diffusion through the cytosol is limited, where these modules are physically close and, in particular, if the NTP-consuming reaction has a very high turnover, i.e. is very processive. Here, we summarise the evidence for these conclusions and describe new insights into the physiological importance and molecular mechanisms of energy channeling gained from recent studies. In particular, we describe the role of glycolytic enzymes for axonal vesicle transport and nucleoside diphosphate kinases for the functions of dynamins and dynamin-related GTPases.

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.

Mechanism of the intrinsic arginine finger in heterotrimeric G proteins

Heterotrimeric G proteins are crucial molecular switches that maintain a large number of physiological processes in cells. The signal is encoded into surface alterations of the Gα subunit that carries GTP in its active state and GDP in its inactive state. The ability of the Gα subunit to hydrolyze GTP is essential for signal termination. Regulator of G protein signaling (RGS) proteins accelerates this process. A key player in this catalyzed reaction is an arginine residue, Arg178 in Gα_i1, which is already an intrinsic part of the catalytic center in Gα in contrast to small GTPases, at which the corresponding GTPase-activating protein (GAP) provides the arginine "finger." We applied time-resolved FTIR spectroscopy in combination with isotopic labeling and site-directed mutagenesis to reveal the molecular mechanism, especially of the role of Arg178 in the intrinsic Gα_i1 mechanism and the RGS4-catalyzed mechanism. Complementary biomolecular simulations (molecular mechanics with molecular dynamics and coupled quantum mechanics/molecular mechanics) were performed. Our findings show that Arg178 is bound to γ-GTP for the intrinsic Gα_i1 mechanism and pushed toward a bidentate α-γ-GTP coordination for the Gα_i1·RGS4 mechanism. This movement induces a charge shift toward β-GTP, increases the planarity of γ-GTP, and thereby catalyzes the hydrolysis.

Theoretical vibrational spectroscopy of intermediates and the reaction mechanism of the guanosine triphosphate hydrolysis by the protein complex Ras-GAP

The structures and vibrational spectra of the reacting species upon guanosine triphosphate (GTP) hydrolysis to guanosine diphosphate and inorganic phosphate (Pi) trapped inside the protein complex Ras-GAP were analyzed following the results of QM/MM simulations. The frequencies of the phosphate vibrations referring to the reactants and to Pi were compared to those observed in the experimental FTIR studies. A good correlation between the theoretical and experimental vibrational data provides a strong support to the reaction mechanism of GTP hydrolysis by the Ras-GAP enzyme system revealed by the recent QM/MM modeling. Evolution of the vibrational bands associated with the inorganic phosphate Pi during the elementary stages of GTP hydrolysis is predicted.

Invited review: Small GTPases and their GAPs

Widespread utilization of small GTPases as major regulatory hubs in many different biological systems derives from a conserved conformational switch mechanism that facilitates cycling between GTP-bound active and GDP-bound inactive states under control of guanine nucleotide exchange factors (GEFs) and GTPase activating proteins (GAPs), which accelerate slow intrinsic rates of activation by nucleotide exchange and deactivation by GTP hydrolysis, respectively. Here we review developments leading to current understanding of intrinsic and GAP catalyzed GTP hydrolytic reactions in small GTPases from structural, molecular and chemical mechanistic perspectives. Despite the apparent simplicity of the GTPase cycle, the structural bases underlying the hallmark hydrolytic reaction and catalytic acceleration by GAPs are considerably more diverse than originally anticipated. Even the most fundamental aspects of the reaction mechanism have been challenging to decipher. Through a combination of experimental and in silico approaches, the outlines of a consensus view have begun to emerge for the best studied paradigms. Nevertheless, recent observations indicate that there is still much to be learned. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 431-448, 2016.

Reaction Mechanism of Guanosine Triphosphate Hydrolysis by the Vision-Related Protein Complex Arl3-RP2

Complexes of small GTPases with GTPase-activating proteins have been intensively studied with the main focus on the complex of H-Ras with p120GAP (Ras-GAP). The detailed mechanism of GTP hydrolysis is still unresolved. To clarify it, we calculated the energy profile of GTP hydrolysis in the active site of a recently characterized vision-related member of this family, the Arl3-RP2 complex. The mechanism suggested in this study retains the main features of GTP hydrolysis by the Ras-GAP complex, but the relative energies of the corresponding intermediates are different and an additional intermediate exists in the Arl3-RP2 complex compared with the Ras-GAP. These differences arise from small deviations in the catalytic arginine conformation of the active site. In the Arl3-RP2 complex, the first two intermediates, corresponding to the Pγ-Oβγ bond cleavage and the glutamine-assisted proton transfer, are almost isoenergetic with the ES complex. Numerical simulations of the kinetic curves demonstrate that the concentrations of these intermediates are comparable with that of ES during the reaction. The calculated IR spectra reveal specific vibrational bands, corresponding to these intermediates. These specific features of the Arl3-RP2 complex open the opportunity to identify spectroscopically two more reaction intermediates in GTP hydrolysis in addition to the ES and EP complexes.

(19)F NMR and DFT Analysis Reveal Structural and Electronic Transition State Features for RhoA-Catalyzed GTP Hydrolysis

Molecular details for RhoA/GAP catalysis of the hydrolysis of GTP to GDP are poorly understood. We use (19)F NMR chemical shifts in the MgF3(-) transition state analogue (TSA) complex as a spectroscopic reporter to indicate electron distribution for the γ-PO3(-) oxygens in the corresponding TS, implying that oxygen coordinated to Mg has the greatest electron density. This was validated by QM calculations giving a picture of the electronic properties of the transition state (TS) for nucleophilic attack of water on the γ-PO3(-) group based on the structure of a RhoA/GAP-GDP-MgF3(-) TSA complex. The TS model displays a network of 20 hydrogen bonds, including the GAP Arg85' side chain, but neither phosphate torsional strain nor general base catalysis is evident. The nucleophilic water occupies a reactive location different from that in multiple ground state complexes, arising from reorientation of the Gln-63 carboxamide by Arg85' to preclude direct hydrogen bonding from water to the target γ-PO3(-) group.

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.

Why does mutation of Gln61 in Ras by the nitro analog NGln maintain activity of Ras-GAP in hydrolysis of guanosine triphosphate?

Interpretation of the experiments showing that the Ras-GAP protein complex maintains activity in guanosine triphosphate (GTP) hydrolysis upon replacement of Glu61 in Ras with its unnatural nitro analog, NGln, is an important issue for understanding details of chemical transformations at the enzyme active site. By using molecular modeling we demonstrate that both glutamine and its nitro analog in the aci-nitro form participate in the reaction of GTP hydrolysis at the stages of proton transfer and formation of inorganic phosphate. The computed structures and the energy profiles for the complete pathway from the enzyme-substrate to enzyme-product complexes for the wild-type and mutated Ras suggest that the reaction mechanism is not affected by this mutation.

Hydrolysis of Guanosine Triphosphate (GTP) by the Ras·GAP Protein Complex: Reaction Mechanism and Kinetic Scheme

Molecular mechanisms of the hydrolysis of guanosine triphosphate (GTP) to guanosine diphosphate (GDP) and inorganic phosphate (Pi) by the Ras·GAP protein complex are fully investigated by using modern modeling tools. The previously hypothesized stages of the cleavage of the phosphorus-oxygen bond in GTP and the formation of the imide form of catalytic Gln61 from Ras upon creation of Pi are confirmed by using the higher-level quantum-based calculations. The steps of the enzyme regeneration are modeled for the first time, providing a comprehensive description of the catalytic cycle. It is found that for the reaction Ras·GAP·GTP·H2O → Ras·GAP·GDP·Pi, the highest barriers correspond to the process of regeneration of the active site but not to the process of substrate cleavage. The specific shape of the energy profile is responsible for an interesting kinetic mechanism of the GTP hydrolysis. The analysis of the process using the first-passage approach and consideration of kinetic equations suggest that the overall reaction rate is a result of the balance between relatively fast transitions and low probability of states from which these transitions are taking place. Our theoretical predictions are in excellent agreement with available experimental observations on GTP hydrolysis rates.

Protein synthesis during cellular quiescence is inhibited by phosphorylation of a translational elongation factor

In nature, most organisms experience conditions that are suboptimal for growth. To survive, cells must fine-tune energy-demanding metabolic processes in response to nutrient availability. Here, we describe a novel mechanism by which protein synthesis in starved cells is down-regulated by phosphorylation of the universally conserved elongation factor Tu (EF-Tu). Phosphorylation impairs the essential GTPase activity of EF-Tu, thereby preventing its release from the ribosome. As a consequence, phosphorylated EF-Tu has a dominant-negative effect in elongation, resulting in the overall inhibition of protein synthesis. Importantly, this mechanism allows a quick and robust regulation of one of the most abundant cellular proteins. Given that the threonine that serves as the primary site of phosphorylation is conserved in all translational GTPases from bacteria to humans, this mechanism may have important implications for growth-rate control in phylogenetically diverse organisms.