Mechanistic analysis of the observed linear free energy relationships in p21ras and related systems

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

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.

Water Distribution within Wild-Type NRas Protein and Q61 Mutants during Unrestrained QM/MM Dynamics

Point mutations in p21^ras are associated with ∼30% of human tumors by disrupting its GTP hydrolysis cycle, which is critical to its molecular switch function in cellular signaling pathways. In this work, we investigate the impact of Gln 61 substitutions in the structure of the p21^N-ras active site and particularly focus on water reorganization around GTP, which appears to be crucial to evaluate favorable and unfavorable hydration sites for hydrolysis. The NRas-GTP complex is analyzed using a hybrid quantum mechanics/molecular mechanics approach, treating for the first time to our knowledge transient water molecules at the ab initio level and leading to results that account for the electrostatic coupling between the protein complex and the solvent. We show that for the wild-type protein, water molecules are found around the GTP γ-phosphate group, forming an arch extended from residues 12 to 35. Two density peaks are observed, supporting previous results that suggest the presence of two water molecules in the active site, one in the vicinity of residue 35 and a second one stabilized by hydrogen bonds formed with nitrogen backbone atoms of residues 12 and 60. The structural changes observed in NRas Gln 61 mutants result in the drastic delocalization of water molecules that we discuss. In mutants Q61H and Q61K, for which water distribution is overlocalized next to residue 60, the second density peak supports the hypothesis of a second water molecule. We also conclude that Gly 60 indirectly participates in GTP hydrolysis by correctly positioning transient water molecules in the protein complex and that Gln 61 has an indirect steric effect in stabilizing the preorganized catalytic site.

The protonation states of GTP and GppNHp in Ras proteins

The small GTPase Ras transmits signals in a variety of cellular signaling pathways, most prominently in cell proliferation. GTP hydrolysis in the active center of Ras acts as a prototype for many GTPases and is the key to the understanding of several diseases, including cancer. Therefore, Ras has been the focus of intense research over the last decades. A recent neutron diffraction crystal structure of Ras indicated a protonated γ-guanylyl imidodiphosphate (γ-GppNHp) group, which has put the protonation state of GTP in question. A possible protonation of GTP was not considered in previously published mechanistic studies. To determine the detailed prehydrolysis state of Ras, we calculated infrared and NMR spectra from quantum mechanics/molecular mechanics (QM/MM) simulations and compared them with those from previous studies. Furthermore, we measured infrared spectra of GTP and several GTP analogs bound to lipidated Ras on a membrane system under near-native conditions. Our findings unify results from previous studies and indicate a structural model confirming the hypothesis that γ-GTP is fully deprotonated in the prehydrolysis state of Ras.

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.

Theory of single-molecule controlled rotation experiments, predictions, tests, and comparison with stalling experiments in F1-ATPase

A recently proposed chemomechanical group transfer theory of rotary biomolecular motors is applied to treat single-molecule controlled rotation experiments. In these experiments, single-molecule fluorescence is used to measure the binding and release rate constants of nucleotides by monitoring the occupancy of binding sites. It is shown how missed events of nucleotide binding and release in these experiments can be corrected using theory, with F₁-ATP synthase as an example. The missed events are significant when the reverse rate is very fast. Using the theory the actual rate constants in the controlled rotation experiments and the corrections are predicted from independent data, including other single-molecule rotation and ensemble biochemical experiments. The effective torsional elastic constant is found to depend on the binding/releasing nucleotide, and it is smaller for ADP than for ATP. There is a good agreement, with no adjustable parameters, between the theoretical and experimental results of controlled rotation experiments and stalling experiments, for the range of angles where the data overlap. This agreement is perhaps all the more surprising because it occurs even though the binding and release of fluorescent nucleotides is monitored at single-site occupancy concentrations, whereas the stalling and free rotation experiments have multiple-site occupancy.

A Microscopic Capacitor Model of Voltage Coupling in Membrane Proteins: Gating Charge Fluctuations in Ci-VSD

The voltage sensitivity of membrane proteins is reflected in the response of the voltage sensing domains (VSDs) to the changes in membrane potential. This response is implicated in the displacement of positively charged residues, associated with the conformational changes of VSDs. The displaced charges generate nonlinear (i.e., voltage-dependent) capacitance current called the gating current (and its corresponding gating charge), which is a key experimental quantity that characterizes voltage activation in VSMP. However, the relevant theoretical/computational approaches, aimed to correlate the structural information on VSMP to electrophysiological measurements, have been rather limited, posing a broad challenge in computer simulations of VSMP. Concomitant with the development of our coarse-graining (CG) model of voltage coupling, we apply our theoretical framework for the treatments of voltage effects in membrane proteins to modeling the VSMP activation, taking the VSDs (Ci-VSD) derived from the Ciona intestinalis voltage sensitive phosphatase (Ci-VSP) as a model system. Our CG model reproduces the observed gating charge of Ci-VSD activation in several different perspectives. In particular, a new closed-form expression of the gating charge is evaluated in both nonequilibrium and equilibrium ways, while considering the fluctuation-dissipation relation that connects a nonequilibrium measurement of the gating charge to an equilibrium measurement of charge fluctuations (i.e., the voltage-independent linear component of membrane capacitance). In turn, the expression uncovers a novel link that connects an equilibrium measurement of the voltage-independent linear capacitance to a nonequilibrium measurement of the voltage-dependent nonlinear capacitance (whose integral over voltage is equal to the gating charge). In addition, our CG model yields capacitor-like voltage dependent free energy parabolas, resulting in the free energy difference and the free energy barrier for the Ci-VSD activation at "zero" (depolarization) membrane potential. Significantly, the resultant voltage dependent energetics enables a direct evaluation of capacitance-voltage relationship (C-V curve) as well as charge-voltage relationship (Q-V curve) that is in a good agreement with the observed measurement of Ci-VSD voltage activation. Importantly, an extension of our kinetic/thermodynamic model of voltage dependent activation in VSMP allows for novel derivations of voltage-dependent rate constants, whose parameters are expressed by the intrinsic properties of VSMP. These novel closed-form expressions offer a physicochemical foundation for the semiempirical Eyring-type voltage dependent rate equations that have been the cornerstone for the phenomenological (kinetic) descriptions of gating and membrane currents in the mechanistic study of ion channels and transporters. Our extended theoretical framework developed in the present study has potential implications on the roles played by gating charge fluctuations for the spike generations in nerve cells within the framework of the Hodgkin-Huxley-type model.

Neutron Crystal Structure of RAS GTPase Puts in Question the Protonation State of the GTP γ-Phosphate

RAS GTPase is a prototype for nucleotide-binding proteins that function by cycling between GTP and GDP, with hydrogen atoms playing an important role in the GTP hydrolysis mechanism. It is one of the most well studied proteins in the superfamily of small GTPases, which has representatives in a wide range of cellular functions. These proteins share a GTP-binding pocket with highly conserved motifs that promote hydrolysis to GDP. The neutron crystal structure of RAS presented here strongly supports a protonated γ-phosphate at physiological pH. This counters the notion that the phosphate groups of GTP are fully deprotonated at the start of the hydrolysis reaction, which has colored the interpretation of experimental and computational data in studies of the hydrolysis mechanism. The neutron crystal structure presented here puts in question our understanding of the pre-catalytic state associated with the hydrolysis reaction central to the function of RAS and other GTPases.

Theory for rates, equilibrium constants, and Brønsted slopes in F1-ATPase single molecule imaging experiments

A theoretical model of elastically coupled reactions is proposed for single molecule imaging and rotor manipulation experiments on F1-ATPase. Stalling experiments are considered in which rates of individual ligand binding, ligand release, and chemical reaction steps have an exponential dependence on rotor angle. These data are treated in terms of the effect of thermodynamic driving forces on reaction rates, and lead to equations relating rate constants and free energies to the stalling angle. These relations, in turn, are modeled using a formalism originally developed to treat electron and other transfer reactions. During stalling the free energy profile of the enzymatic steps is altered by a work term due to elastic structural twisting. Using biochemical and single molecule data, the dependence of the rate constant and equilibrium constant on the stall angle, as well as the Børnsted slope are predicted and compared with experiment. Reasonable agreement is found with stalling experiments for ATP and GTP binding. The model can be applied to other torque-generating steps of reversible ligand binding, such as ADP and Pi release, when sufficient data become available.

Overview of simulation studies on the enzymatic activity and conformational dynamics of the GTPase Ras

Over the last 40 years, we have learnt a great deal about the Ras onco-proteins. These intracellular molecular switches are essential for the function of a variety of physiological processes, including signal transduction cascades responsible for cell growth and proliferation. Molecular simulations and free energy calculations have played an essential role in elucidating the conformational dynamics and energetics underlying the GTP hydrolysis reaction catalysed by Ras. Here we present an overview of the main lessons from molecular simulations on the GTPase reaction and conformational dynamics of this important anti-cancer drug target. In the first part, we summarise insights from quantum mechanical and combined quantum mechanical/molecular mechanical simulations as well as other free energy methods and highlight consensus viewpoints as well as remaining controversies. The second part provides a very brief overview of new insights emerging from large-scale molecular dynamics simulations. We conclude with a perspective regarding future studies of Ras where computational approaches will likely play an active role.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.

On the Origins of the Linear Free Energy Relationships: Exploring the Nature of the Off-Diagonal Coupling Elements in S(N)2 Reactions

Understanding the relationship between the adiabatic free energy profiles of chemical reactions and the underlining diabatic states is central to the description of chemical reactivity. The diabatic states form the theoretical basis of Linear Free Energy Relationships (LFERs) and thus play a major role in physical organic chemistry and related fields. However, the theoretical justification for some of the implicit LFER assumptions has not been fully established by quantum mechanical studies. This study follows our earlier works(1,2) and uses the ab initio frozen density functional theory (FDFT) method(3) to evaluate both the diabatic and adiabatic free energy surfaces and to determine the corresponding off-diagonal coupling matrix elements for a series of S(N)2 reactions. It is found that the off-diagonal coupling matrix elements are almost the same regardless of the nucleophile and the leaving group but change upon changing the central group. Furthermore, it is also found that the off diagonal elements are basically the same in gas phase and in solution, even when the solvent is explicitly included in the ab initio calculations. Furthermore, our study establishes that the FDFT diabatic profiles are parabolic to a good approximation thus providing a first principle support to the origin of LFER. These findings further support the basic approximation of the EVB treatment.

The specific vibrational modes of GTP in solution and bound to Ras: a detailed theoretical analysis by QM/MM simulations

The hydrolysis of guanosine triphosphate (GTP) in general, and especially by GTPases like the Ras protein, is in the focus of biological investigations. A huge amount of experimental data from Fourier-transformed infrared studies is currently available, and many vibrational bands of free GTP, GTP·Mg(2+), and Ras·GTP·Mg(2+) in solution have been assigned by isotopic labeling. In the Ras environment, bands between 800 cm(-1) and 1300 cm(-1) have already been assigned, but not those below 800 cm(-1). The combination of quantum and molecular mechanics (QM/MM) methods takes the quantum effects for selected relevant atoms into account. This provides structural details, vibrational frequencies and electron distributions of the region of interest. We therefore used MM and QM/MM simulations to investigate the normal vibrational modes of GTP, GTP·Mg(2+), and Ras·GTP·Mg(2+) in solution, and assigned the vibrational frequencies for each normal vibration mode. In this study, the quantum box contains the nucleoside and the Mg(2+). The comparison of calculated and experimental vibrational spectra provides a very good control for the quality of the calculations. Structurally, MM and QM/MM simulations reveal a stable tridentate coordination of the Mg(2+) by GTP in water, and a stable bidentate coordination by GTP in complex with Ras. For validation, we compare the calculated frequencies and isotopic shifts with the experimental results available in the range of 800 cm(-1) to 1300 cm(-1). For the first time we suggest band assignments of the vibrational modes below 800 cm(-1) by comparison of calculated and experimental spectra.

ATP hydrolysis in Eg5 kinesin involves a catalytic two-water mechanism

Motor proteins couple steps in ATP binding and hydrolysis to conformational switching both in and remote from the active site. In our kinesin.AMPPPNP crystal structure, closure of the active site results in structural transformations appropriate for microtubule binding and organizes an orthosteric two-water cluster. We conclude that a proton is shared between the lytic water, positioned for gamma-phosphate attack, and a second water that serves as a general base. To our knowledge, this is the first experimental detection of the catalytic base for any ATPase. Deprotonation of the second water by switch residues likely triggers subsequent large scale structural rearrangements. Therefore, the catalytic base is responsible for initiating nucleophilic attack of ATP and for relaying the positive charge over long distances to initiate mechanotransduction. Coordination of switch movements via sequential proton transfer along paired water clusters may be universal for nucleotide triphosphatases with conserved active sites, such as myosins and G-proteins.

The EVB as a quantitative tool for formulating simulations and analyzing biological and chemical reactions

Recent years have seen dramatic improvements in computer power, allowing ever more challenging problems to be approached. In light of this, it is imperative to have a quantitative model for examining chemical reactivity, both in the condensed phase and in solution, as well as to accurately quantify physical organic chemistry (particularly as experimental approaches can often be inconclusive). Similarly, computational approaches allow for great progress in studying enzyme catalysis, as they allow for the separation of the relevant energy contributions to catalysis. Due to the complexity of the problems that need addressing, there is a need for an approach that can combine reliability with an ability to capture complex systems in order to resolve long-standing controversies in a unique way. Herein, we will demonstrate that the empirical valence bond (EVB) approach provides a powerful way to connect the classical concepts of physical organic chemistry to the actual energies of enzymatic reactions by means of computation. Additionally, we will discuss the proliferation of this approach, as well as attempts to capture its basic chemistry and repackage it under different names. We believe that the EVB approach is the most powerful tool that is currently available for studies of chemical processes in the condensed phase in general and enzymes in particular, particularly when trying to explore the different proposals about the origin of the catalytic power of enzymes.

Dineopentyl phosphate hydrolysis: evidence for stepwise water attack

Phosphate ester hydrolysis is ubiquitous in biology, playing a central role in energy production, signaling, biosynthesis, and the regulation of protein function among other things. Although the mechanism of action of the enzymes regulating this reaction has been the focus of intensive research in the past few decades, the correct description of this apparently simple reaction remains controversial. A clear understanding of the mechanism that takes place in solution is crucial to be able to evaluate whether proposals for the enzyme-catalyzed mechanisms are reasonable. For the pH-independent hydrolysis of phosphate diesters, several kinetically equivalent mechanisms are plausible, including hydroxide attack on the neutral phosphate. However, it is very difficult to measure the rate of this reaction directly by experimental methods, so it has been evaluated by examining the rate of hydrolysis of neutral phosphate triesters, where a methyl group has replaced a proton. This may not be an accurate model of the neutral phosphate diester and does not provide information about a reaction pathway that is concerted with nucleophilic attack to generate a similar phosphorane. We have carefully mapped out free energy surfaces for both hydroxide and water attack on the dineopentyl phosphate anion and for water attack on the neutral diester. In doing so, we have accurately reproduced existing experimental data and demonstrate that water attack proceeds through an associative mechanism with proton transfer to the phosphate to generate a phosphorane intermediate. Our data show that the substrate-as-base mechanism is viable for phosphate ester hydrolysis, which may have important implications for the studies of phosphate ester hydrolysis by enzymes.

Semisynthesis of H-Ras with a glutamic acid methylester at position 61

The mechanism of the guanosine triphosphate (GTP) hydrolysis reaction of small G-proteins such as Ras is generally understood; however, some important molecular details are still missing. One example concerns the role of Gln61 in the catalysis of the GTP hydrolysis reaction. This amino acid is frequently mutated in oncogenic Ras leading to constitutively active variants of the protein. To elucidate the role of Gln61, subtle structural changes were introduced at this position by exchanging the natural occurring glutamine against a glutamic acid methyl ester (GluOme). Thereby the H-bond donor properties of this residue are changed and analysis of the GTP hydrolysis reaction can provide information on the function of the native carboxamide moiety. Using a semisynthetic approach, Ras(1-166)Gln61GluOMe was synthesized by sequential native chemical ligation of three unprotected peptide segments. Peptides Ras(1-50) and Ras(51-79)Gln61GluOMe were synthesized using Boc chemistry. The C-terminal peptide Ras(80-166) was expressed in E. coli. Initial tests of this semisynthetic strategy were performed by synthesis of the N- and C-terminally truncated protein variant Ras(39-101)Gln61GluOMe. The identified optimal reaction conditions were then applied to the synthesis of Ras(1-166)Gln61GluOMe. Refolding of the semisynthetic product in the presence of GTP was successful and revealed intrinsic GTPase activity of Ras(1-166)Gln61GluOMe.

On the interpretation of the observed linear free energy relationship in phosphate hydrolysis: a thorough computational study of phosphate diester hydrolysis in solution

The hydrolysis of phosphate esters is crucially important to biological systems, being involved in, among other things, signaling, energy transduction, biosynthesis, and the regulation of protein function. Despite this, there are many questions that remain unanswered in this important field, particularly with regard to the preferred mechanism of hydrolysis of phosphate esters, which can proceed through any of multiple pathways that are either associative or dissociative in nature. Previous comparisons of calculated and observed linear free energy relationships (LFERs) for phosphate monoester dianions with different leaving groups showed that the TS character gradually changes from associative to dissociative with the increasing acidity of the leaving group, while reproducing the experimental LFER. Here, we have generated ab initio potential energy surfaces for the hydrolysis of phosphate diesters in solution, with a variety of leaving groups. Once again, the reaction changes from a compact concerted pathway to one that is more expansive in character when the acidity of the leaving group increases. When such systems are examined in solution, it is essential to take into consideration the contribution of solute to the overall activation entropy, which remains a major computational challenge. The popular method of calculating the entropy using a quasi-harmonic approximation appears to markedly overestimate the configurational entropy for systems with multiple occupied energy wells. We introduce an improved restraint release approach for evaluating configurational entropies and apply this approach to our systems. We demonstrate that when this factor is taken into account, it is possible to reproduce the experimental LFER for this system with reasonable accuracy.

H and other transfers in enzymes and in solution: theory and computations, a unified view. 2. Applications to experiment and computations

Equations obtained in part I for the free-energy barrier to one-step enzymatic reactions between bound reactants are discussed. The rate is expressed in terms of lambdao (protein reorganization energy), DeltaG(o) (standard free energy of reaction of the H-transfer step), bond breaking/bond forming term, w (work terms), and H-transmission property. Two alternative approximations for the coupling of the bond breaking/bond forming and protein are distinguished experimentally in favorable cases by the DeltaG(o) where the maximum deuterium kinetic isotope effect occurs. Plots of log rate versus DeltaG(o) and properties such as DeltaS* and DeltaS(o) are discussed. The weak or zero T-dependence of the kinetic isotope effect for wild-type enzymes operating under physiological conditions is interpreted in terms of vanishing (or isotopically insensitive) w plus transfer from the lowest H-state. Static and dynamic protein flexibility is discussed. While the many correlations accessible for electron transfers are not available for H-transfers in enzymes, a combination of experiment, computation, and analytical approaches can assist in evaluating the utility of the present equations and in suggesting further experiments and computations. A protein reorganization energy lambdao is obtained in the literature from the extended valence bond formalism where diabatic electronic states are used. A method is suggested for extracting it when instead a bond distance difference coordinate is used. The results may provide a bridge between the two approaches.

Mechanisms of guanosine triphosphate hydrolysis by Ras and Ras-GAP proteins as rationalized by ab initio QM/MM simulations

The hydrolysis reaction of guanosine triphosphate (GTP) by p21(ras) (Ras) has been modeled by using the ab initio type quantum mechanical-molecular mechanical simulations. Initial geometry configurations have been prompted by atomic coordinates of the crystal structure (PDBID: 1QRA) corresponding to the prehydrolysis state of Ras in complex with GTP. Multiple searches of minimum energy geometry configurations consistent with the hydrogen bond networks have been performed, resulting in a series of stationary points on the potential energy surface for reaction intermediates and transition states. It is shown that the minimum energy reaction path is consistent with an assumption of a two-step mechanism of GTP hydrolysis. At the first stage, a unified action of the nearest residues of Ras and the nearest water molecules results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low barrier (16.7 kcal/mol) transition state TS1. At the second stage, the inorganic phosphate is formed in consequence of proton transfers mediated by two water molecules and assisted by the Gln61 residue from Ras. The highest transition state at this segment, TS3, is estimated to have an energy 7.5 kcal/mol above the enzyme-substrate complex. The results of simulations are compared to the previous findings for the GTP hydrolysis in the Ras-GAP (p21(ras)-p120(GAP)) protein complex. Conclusions of the modeling lead to a better understanding of the anticatalytic effect of cancer causing mutation of Gln61 from Ras, which has been debated in recent years.

Thermodynamics of CFTR channel gating: a spreading conformational change initiates an irreversible gating cycle

CFTR is the only ABC (ATP-binding cassette) ATPase known to be an ion channel. Studies of CFTR channel function, feasible with single-molecule resolution, therefore provide a unique glimpse of ABC transporter mechanism. CFTR channel opening and closing (after regulatory-domain phosphorylation) follows an irreversible cycle, driven by ATP binding/hydrolysis at the nucleotide-binding domains (NBD1, NBD2). Recent work suggests that formation of an NBD1/NBD2 dimer drives channel opening, and disruption of the dimer after ATP hydrolysis drives closure, but how NBD events are translated into gate movements is unclear. To elucidate conformational properties of channels on their way to opening or closing, we performed non-equilibrium thermodynamic analysis. Human CFTR channel currents were recorded at temperatures from 15 to 35°C in inside-out patches excised from Xenopus oocytes. Activation enthalpies(ΔH^‡) were determined from Eyring plots. ΔH^‡ was 117 ± 6 and 69 ± 4 kJ/mol, respectively, for opening and closure of partially phosphorylated, and 96 ± 6 and 73 ± 5 kJ/mol for opening and closure of highly phosphorylated wild-type (WT) channels. ΔH^‡ for reversal of the channel opening step, estimated from closure of ATP hydrolysis–deficient NBD2 mutant K1250R and K1250A channels, and from unlocking of WT channels locked open with ATP+AMPPNP, was 43 ± 2, 39 ± 4, and 37 ± 6 kJ/mol, respectively. Calculated upper estimates of activation free energies yielded minimum estimates of activation entropies (ΔS^‡), allowing reconstruction of the thermodynamic profile of gating, which was qualitatively similar for partially and highly phosphorylated CFTR. ΔS^‡ appears large for opening but small for normal closure. The large ΔH^‡ and ΔS^‡ (TΔS^‡ ≥ 41 kJ/mol) for opening suggest that the transition state is a strained channel molecule in which the NBDs have already dimerized, while the pore is still closed. The small ΔS^‡ for normal closure is appropriate for cleavage of a single bond (ATP's beta-gamma phosphate bond), and suggests that this transition state does not require large-scale protein motion and hence precedes rehydration (disruption) of the dimer interface.

QM/MM modeling the Ras-GAP catalyzed hydrolysis of guanosine triphosphate

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras-GAP (p21(ras) - p120(GAP)) has been modeled by the quantum mechanical-molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two-step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low-barrier transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the gamma-phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H(2)O --> GDP + H(2)PO(4) (-) can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer-causing mutations of Ras, which has been debated in recent years.

Dissection of the GTPase mechanism of Ras protein by MD analysis of Ras mutants

Controlling the hydrolysis rate of GTP bound to the p21ras protein is crucial for the delicate timing of many biological processes. A few mechanisms were suggested for the hydrolysis of GTP. To gain more insight into the individual elementary events of GTP hydrolysis, we carried out molecular dynamic analysis of wild-type p21ras and some of its mutants. It was recently shown that Ras-related proteins and mutants generally follow a linear free energy relationship (LFER) relating the rate of reaction to the pK(a) of the gamma-phosphate group of the bound GTP, indicating that proton transfer from the attacking water to the GTP is the first elementary event in the GTPase mechanism. However, some exceptions were observed. Thus, the Gly12 --> Aspartic p21ras (G12D) mutant had a very low GTPase activity although its pK(a) was very close to that of the wild-type ras. Here we compared the molecular dynamics (MD) of wild-type Ras and G12D, showing that in the mutant the catalytic water molecule is displaced to a position where proton transfer to GTP is unfavorable. These results suggest that the mechanism of GTPase is indeed composed of an initial proton abstraction from water by the GTP, followed by a nucleophilic attack of the hydroxide ion on the gamma-phosphorus of GTP.

H662 is the linchpin of ATP hydrolysis in the nucleotide-binding domain of the ABC transporter HlyB

The ABC transporter HlyB is a central element of the HlyA secretion machinery, a paradigm of Type I secretion. Here, we describe the crystal structure of the HlyB-NBD (nucleotide-binding domain) with H662 replaced by Ala in complex with ATP/Mg2+. The dimer shows a composite architecture, in which two intact ATP molecules are bound at the interface of the Walker A motif and the C-loop, provided by the two monomers. ATPase measurements confirm that H662 is essential for activity. Based on these data, we propose a model in which E631 and H662, highly conserved among ABC transporters, form a catalytic dyad. Here, H662 acts as a 'linchpin', holding together all required parts of a complicated network of interactions between ATP, water molecules, Mg2+, and amino acids both in cis and trans, necessary for intermonomer communication. Based on biochemical experiments, we discuss the hypothesis that substrate-assisted catalysis, rather than general base catalysis might operate in ABC-ATPases.

Theoretical IR spectroscopy based on QM/MM calculations provides changes in charge distribution, bond lengths, and bond angles of the GTP ligand induced by the Ras-protein

The GTPase Ras p21 is a crucial switch in cellular signal transduction. Fourier transform infrared (FTIR) spectra of the substrate guanosine triphosphate (GTP) show remarkable changes when it binds to the enzyme. The reduced band widths indicate that the flexible GTP molecule is guided by the protein into a preferred conformation. The delocalized phosphate vibrations of unbound GTP become localized. The frequency shifts show an electron movement toward beta-phosphate, which probably contributes to catalysis by reducing the free activation energy. To quantify these qualitative observations we performed QM/MM molecular dynamics simulations of Ras.GTP and GTP in water. The triphosphate part of GTP was treated quantum mechanically using density functional theory (DFT). Vibrational spectra were calculated in harmonic approximation with an average deviation of 3% from the experimental frequencies. This provides a high confidence in the computational results as vibrational spectra are highly sensitive to conformation and charge distribution. As compared to GTP in water, Ras-bound GTP shows a shift of negative charge of approximately 0.2 e toward the beta-phosphate from gamma-phosphate and from alpha-phosphate due to the positive charge of the magnesium ion, to a lesser extent of Lys-16, and surprisingly without any effect of the P-loop backbone. Magnesium and Gly-13 twist and bend the gamma-O-beta bonds such that the crucial bond is stretched before cleaving.

GTP hydrolysis mechanism of Ras-like GTPases

The Ras-like GTPases regulate diverse cellular functions via the chemical cycle of binding and hydrolyzing GTP molecules. They alternate between GTP- and GDP-bound conformations. The GTP-bound conformation is biologically active and promotes a cellular function, such as signal transduction, cytoskeleton organization, protein synthesis/translocation, or a membrane budding/fusion event. GTP hydrolysis turns off the GTPase switch by converting it to the inactive GDP-bound conformation. The fundamental GTP hydrolysis mechanism by these GTPases has generated considerable interest over the last two decades but remained to be firmly established. This review provides an update on the catalytic mechanism with discussions on recent developments from kinetic, structural, and model studies in the context of the various GTP hydrolysis models proposed over the years.

Quantum chemical modeling of the GTP hydrolysis by the RAS-GAP protein complex

We present results of the modeling for the hydrolysis reaction of guanosine triphosphate (GTP) in the RAS-GAP protein complex using essentially ab initio quantum chemistry methods. One of the approaches considers a supermolecular cluster composed of 150 atoms at a consistent quantum level. Another is a hybrid QM/MM method based on the effective fragment potential technique, which describes interactions between quantum and molecular mechanical subsystems at the ab initio level of the theory. Our results show that the GTP hydrolysis in the RAS-GAP protein complex can be modeled by a substrate-assisted catalytic mechanism. We can locate a configuration on the top of the barrier corresponding to the transition state of the hydrolysis reaction such that the straightforward descents from this point lead either to reactants GTP+H(2)O or to products guanosine diphosphate (GDP)+H(2)PO(4)(-). However, in all calculations such a single-step process is characterized by an activation barrier that is too high. Another possibility is a two-step reaction consistent with formation of an intermediate. Here the Pgamma-O(Pbeta) bond is already broken, but the lytic water molecule is still in the pre-reactive state. We present arguments favoring the assumption that the first step of the GTP hydrolysis reaction in the RAS-GAP protein complex may be assigned to the breaking of the Pgamma-O(Pbeta) bond prior to the creation of the inorganic phosphate.

What really prevents proton transport through aquaporin? Charge self-energy versus proton wire proposals

The nature of the control of water/proton selectivity in biological channels is a problem of a fundamental importance. Most studies of this issue have proposed that an interference with the orientational requirements of the so-called proton wire is the source of selectivity. The elucidation of the structures of aquaporins, which have evolved to prevent proton transfer (PT), provided a clear benchmark for exploring the selectivity problem. Previous simulations of this system have not examined, however, the actual issue of PT, but only considered the much simpler task of the transfer of water molecules. Here we take aquaporin as a benchmark and quantify the origin of the water/proton selectivity in this and related systems. This is done by evaluating in a consistent way the free energy profile for transferring a proton along the channel and relating this profile to the relevant PT rate constants. It is found that the water/proton selectivity is controlled by the change in solvation free energy upon moving the charged proton from water to the channel. The reason for the focus on the elegant concept of the proton wire and the related Grotthuss-type mechanism is also considered. It is concluded that these mechanisms are clearly important in cases with flat free energy surfaces (e.g., in bulk water, in gas phase water chains, and in infinitely long channels). However, in cases of biological channels, the actual PT mechanism is much less important than the energetics of transferring the proton charge from water to different regions in the channels.

Kinetic isotope effects in Ras-catalyzed GTP hydrolysis: evidence for a loose transition state

A remote labeling method has been developed to determine (18)O kinetic isotope effects (KIEs) in Ras-catalyzed GTP hydrolysis. Substrate mixtures consist of (13)C-depleted GTP and [(18)O,(13)C]GTP that contains (18)O at phosphoryl positions of mechanistic interest and (13)C at all carbon positions of the guanosine moiety. Isotope ratios of the nonvolatile substrates and products are measured by using a chemical reaction interface/isotope ratio mass spectrometer. The isotope effects are 1.0012 (0.0026) in the gamma nonbridge oxygens, 1.0194 (0.0025) in the leaving group oxygens (the beta-gamma oxygen and the two beta nonbridge oxygens), and 1.0105 (0.0016) in the two beta nonbridge oxygens. The KIE in the beta-gamma bridge oxygen was computed to be 1.0116 or 1.0088 by two different methods. The significant KIE in the leaving group reveals that chemistry is largely rate-limiting whereas the KIEs in the gamma nonbridge oxygens and the leaving group indicate a loose transition state that approaches a metaphosphate. The KIE in the two beta nonbridge oxygens is roughly equal to that in the beta-gamma bridge oxygen. This indicates that, in the transition state, Ras shifts one-half of the negative charge that arises from P(gamma)-O(beta-gamma) fission from the beta-gamma bridge oxygen to the two beta nonbridge oxygens. The KIE effects, interpreted in light of structural and spectroscopic data, suggest that Ras promotes a loose transition state by stabilizing negative charge in the beta-gamma bridge and beta nonbridge oxygens of GTP.

Converting conformational changes to electrostatic energy in molecular motors: The energetics of ATP synthase

F1-ATPase is the catalytic component of the ATP synthase molecular machine responsible for most of the uphill synthesis of ATP in living systems. The enormous advances in biochemical and structural studies of this machine provide an opportunity for detailed understanding of the nature of its rotary mechanism. However, further quantitative progress in this direction requires development of reliable ways of translating the observed structural changes to the corresponding energies. This requirement is particularly challenging because we are dealing with a large system that couples major structural changes with a chemical process. The present work provides such a structure-function correlation by using the linear response approximation to describe the rotary mechanism. This approach allows one to evaluate the energy of transitions between different conformational states by considering only the changes in the corresponding electrostatic energies of the ligands. The relevant energetics are also obtained by calculating the linear response approximation-based free energies of transferring the ligands from water to the different sites of F1-ATPase in their different conformational states. We also use the empirical valence bond approach to evaluate the actual free-energy profile for the ATP synthesis in the different conformational states of the system. Integrating the information from the different approaches provides a semiquantitative structure-function correlation for F1-ATPase. It is found that the conformational changes are converted to changes in the electrostatic interaction between the protein and its ligands, which drives the ATP synthesis.

Essential role of histidine 84 in elongation factor Tu for the chemical step of GTP hydrolysis on the ribosome

Elongation factor Tu (EF-Tu) is a GTP-binding protein that delivers aminoacyl-tRNA to the A site of the ribosome during protein synthesis. The mechanism of GTP hydrolysis in EF-Tu on the ribosome is poorly understood. It is known that mutations of a conserved histidine residue in the switch II region of the factor, His84 in Escherichia coli EF-Tu, impair GTP hydrolysis. However, the partial reaction which is directly affected by mutations of His84 was not identified and the effect on GTP hydrolysis was not quantified. Here, we show that the replacement of His84 with Ala reduces the rate constant of GTP hydrolysis more than 10(6)-fold, whereas the preceding steps of ternary complex binding to the ribosome, codon recognition and, most importantly, the GTPase activation step are affected only slightly. These results show that His84 plays a key role in the chemical step of GTP hydrolysis. Rate constants of GTP hydrolysis by wild-type EF-Tu, measured using the slowly hydrolyzable GTP analog, GTPgammaS, showed no dependence on pH, indicating that His84 does not act as a general base. We propose that the catalytic role of His84 is to stabilize the transition state of GTP hydrolysis by hydrogen bonding to the attacking water molecule or, possibly, the gamma-phosphate group of GTP.

GTPase catalysis by Ras and other G-proteins: insights from Substrate Directed SuperImposition

Comparisons of different protein structures are commonly carried out by superimposing the coordinates of the protein backbones or selected parts of the proteins. When the objective is analysis of similarities and differences in the enzyme's active site, there is an inherent problem in using the same domains for the superimposition. In this work we use a comparative approach termed here "Substrate Directed SuperImposition" (SDSI). It entails the superimposition of multiple protein-substrate structures using exclusively the coordinates of the comparable substrates. SDSI has the advantage of unbiased comparison of the active-site environment from the substrate's point of view. Our analysis extends previous usage of similar approaches to comparison of enzyme catalytic machineries. We applied SDSI to various G-protein structures for dissecting the mechanism of the GTPase reaction that controls the signaling activity of this important family. SDSI indicates that dissimilar G-proteins stabilize the transition state of the GTPase reaction similarly and supports the commonality of the critical step in this reaction, the reorientation of the critical arginine and glutamine. Additionally, we ascribe the catalytic inefficiency of the small G-protein Ras to the great flexibility of its active site and downplay the possible catalytic roles of the Lys16 residue in Ras GTPase. SDSI demonstrated that in contrast to all other Gly12 Ras mutants, which are oncogenic, the Gly12-->Pro mutant does not interfere with the catalytic orientation of the critical glutamine. This suggests why this mutant has a higher rate of GTP hydrolysis and is non-transforming. Remarkably, SDSI also revealed similarities in the divergent catalytic machineries of G-proteins and UMP/CMP kinase. Taken together, our results promote the use of SDSI to compare the catalytic machineries of both similar and different classes of enzymes.

Oncogenic mutations and packing defects in protein structure

Oncogenic mutations in expressed proteins are of primary interest to understand tumor formation but their structural consequences bearing on protein function are not clearly understood. In this contribution I report on two illustrative examples, p21ras and p57, revealing that such mutations have an effect on specific structural deficiencies in the packing of the protein structure, i. e., on backbone hydrogen bonds insufficiently shielded from water attack. These structural deficiencies in the wild type are typically "corrected intermolecularly" by protein complexation or protein-ligand association. However, in the oncogenic mutants, these binding signals are partially or completely suppressed: the mutated residues properly wrap or desolvate the hydrogen bonds intramolecularly. Thus, the interactivity of the proteins becomes impaired: their binding affinity decreases sharply, as there is no thermodynamic benefit from removing water surrounding properly desolvated hydrogen bonds. The results, specialized for p21ras and p53, reveal how oncogenic mutations determine a hindrance to GAP-induced hydrolysis (p21) and decrease binding affinity for DNA (p53). Furthermore, the oncogenic potential of mutations in residues not directly engaged in the interface electrostatics is assessed. The results suggest that a high sensitivity of structural defects to genetic accident might be a necessary condition to establish the existence of a proto-oncogene, an angle that merits a systematic study.

Revisiting the structural flexibility of the complex p21(ras)-GTP: the catalytic conformation of the molecular switch II

The hydrolysis of GTP in p21(ras) triggers conformational changes that regulate the ras/ERK signaling pathway. An important active site residue is Gln61, which has been found to be mutated in 30% of human tumors. The dynamics of the active site conformation is studied by using molecular dynamics simulation of two independent structures of the GTP-bound uncomplexed enzyme. Two distinct conformations of the enzyme are observed, in which the side-chain residue Gln61 is in different orientations. Essential dynamics analysis is used to describe the essential motions in the transition between the two conformations. Results are compared with earlier simulations of p21(ras) and its complex with GTPase activating protein p21-GAP.

gamma-phosphate protonation and pH-dependent unfolding of the Ras.GTP.Mg2+ complex: a vibrational spectroscopy study

The interdependence of GTP hydrolysis and the second messenger functions of virtually all GTPases has stimulated intensive study of the chemical mechanism of the hydrolysis. Despite numerous mutagenesis studies, the presumed general base, whose role is to activate hydrolysis by abstracting a proton from the nucleophilic water, has not been identified. Recent theoretical and experimental work suggest that the gamma-phosphate of GTP could be the general base. The current study investigates this possibility by studying the pH dependence of the vibrational spectrum of the Ras.GTP.Mg(2+) and Ras.GDP.Mg(2+) complexes. Isotope-edited IR studies of the Ras.GTP.Mg(2+) complex show that GTP remains bound to Ras at pH as low as 2.0 and that the gamma-phosphate is not protonated at pH > or = 3.3, indicating that the active site decreases the gamma-phosphate pK(a) by at least 1.1 pK(a) units compared with solution. Amide I studies show that the Ras.GTP.Mg(2+) and Ras.GDP.Mg(2+) complexes partially unfold in what appear to be two transitions. The first occurs in the pH range 5.4-2.6 and is readily reversible. Differences in the pH-unfolding midpoints for the Ras.GTP.Mg(2+) and Ras.GDP.Mg(2+) complexes (3.7 and 4.8, respectively) reveal that the enzyme-gamma-phosphoryl interactions stabilize the structure. The second transition, pH 2.6-1.7, is not readily reversed. The pH-dependent unfolding of the Ras.GTP.Mg(2+) complex provides an alternative interpretation of the data that had been used to support the gamma-phosphate mechanism, thereby raising the issue of whether this mechanism is operative in GTPase-catalyzed GTP hydrolysis reactions.

Substrate assisted catalysis -- application to G proteins

The idea that both the substrate and the enzyme contribute to catalysis (substrate assisted catalysis; SAC) is applicable to guanine nucleotide-binding proteins (G proteins). Naturally occurring SAC uses GTP as a general base in the GTPase reaction catalyzed by G proteins. Engineered SAC has identified a putative rate-limiting step for the GTPase reaction and shown that GTPase-deficient oncogenic Ras mutants are not irreversibly impaired. Thus, anti-cancer drugs could potentially be designed to restore the blocked GTPase reaction.

Interaction kinetics of reversible inhibitors and substrates with acetylcholinesterase and its fasciculin 2 complex

Fasciculin 2 (Fas2), a three-fingered peptide of 61 amino acids, binds tightly to the peripheral site of acetylcholinesterases (AChE; EC ), occluding the entry portal into the active center gorge of the enzyme and inhibiting its catalytic activity. We investigated the mechanism of Fas2 inhibition by studying hydrolysis of cationic and neutral substrates and by determining the kinetics of interaction for fast equilibrating cationic and neutral reversible inhibitors with the AChE.Fas2 complex and free AChE. Catalytic parameters, derived by eliminating residual Fas2-resistant activity, reveal that Fas2 reduces k(cat)/K(m) up to 10(6)-fold for cationic substrates and less than 10(3)-fold for neutral substrates. Rate constants for association of reversible inhibitors with the active center of the AChE.Fas2 complex were reduced about 10(4)-fold for both cationic and neutral inhibitors, while dissociation rate constants were reduced 10(2)-to 10(3)-fold, compared with AChE alone. Rates of ligand association with both AChE and AChE.Fas2 complex were dependent on the protonation state of ionizable ligands but were also markedly reduced by protonation of enzyme residue(s) with pK(a) of 6.1-6.2. Linear free energy relationships between the equilibrium constant and the kinetic constants show that Fas2, presumably through an allosteric influence, markedly alters the position of the transition state in the reaction pathway. Since Fas2 complexation introduces an energetic barrier for hydrolysis of substrates that exceeds that found for association of reversible ligands, Fas2 influences catalytic parameters by a more complex mechanism than simple restriction of diffusional entry and exit from the active center. Conformational flexibility appears critical for facilitating ligand passage in the narrow active center gorge for both AChE and the AChE.Fas2 complex.

Application of linear free energy relationships to the serpin-proteinase inhibition mechanism

Linear free energy relationships can be used to link the changes in rate constant for a reaction to changes in the equilibrium caused by alterations in structure. While they have most often been used in the analysis of chemical reactions, they have also been employed to resolve questions in enzymology and protein folding. Here we analyze the reaction of a serpin with a panel of six serine proteinases, and observe that a linear free energy relationship exists between the true second-order rate constant for reaction, k(inh), and the inhibition constant, K(I), indicating that formation of the covalent serpin-enzyme complex may be reversible.