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

FT-IR Method Limitations for β-Glucan Analysis

β-glucans are known as biological response modifiers. However, different sources can result in structural differences and as a result differences in their biological activity. The hot water extraction method allows to obtain, high molecular weight β-glucans without altering their structure by using strong chemicals, such as alkalis or acids. Analysis of β-glucans by FT-IR and NMR spectroscopy in solid state is superior to analysis in solution as it allows researchers to study the preserved structure of the extracted polysaccharides. FT-IR spectroscopy was used in this study to make side-by-side comparison analysis of hot water extracted β-glucans from different yeast sources. NMR spectroscopy was used to confirm findings made by FT-IR spectroscopy. Extracted β-glucans exhibit characteristic structure of β-1,3/1,6-linked glucans with noticeable levels of proteins, possibly in a form of oligopeptides, chitin and other impurities. β-glucans obtained from C. guilliermondii, P. pastoris and S. pastorianus exhibited higher protein content. Differences in mannan, chitin and α-glucan content were also observed; however, the species-specific structure of obtained β-glucans could not be confirmed without additional studies. Structural analysis of high molecular weight β-glucans in solid state by FT-IR spectroscopy is difficult or limited due to band intensity changes and overlapping originating from different molecules.

Differences in the Nature of the Phosphoryl Transfer Transition State in Protein Phosphatase 1 and Alkaline Phosphatase: Insights from QM Cluster Models

Quantum mechanical (QM) cluster models are used to probe effects on the catalytic properties of protein phosphatase 1 (PP1) and alkaline phosphatase (AP) due to metal ions and active site residues. The calculations suggest that the phosphoryl transfer transition states in PP1 are synchronous in nature with a significant degree of P-Olg cleavage, while those in AP are tighter with a modest degree of P-Olg cleavage and a range of P-Onuc formation. Similar to observations made in our recent work, a significant degree of cross talk between the forming and breaking P-O bonds complicates the interpretation of the Brønsted relation, especially in regard to AP for which the computed βlg/βEQ,lg value does not correlate with the degree of P-Olg cleavage regardless of the metal ions in the active site. By comparison, the correlation between βlg/βEQ,lg and the P-Olg bond order is more applicable to PP1, which generally exhibits less variation in the transition state than AP. Results for computational models with swapped metal ions between PP1 and AP suggest that the metal ions modulate both the nature of the transition state and the degrees of sensitivity of the transition state to the leaving group. In the reactant state, the degree of the scissile bond polarization is also different in the two enzymes, although this difference appears to be largely determined by the active site residues rather than the metal ions. Therefore, both the identity of the metal ion and the positioning of polar or charged residues in the active site contribute to the distinct catalytic characteristics of these enzymes. Several discrepancies observed between the QM cluster results and the available experimental data highlight the need for further QM/MM method developments for the quantitative analysis of metalloenzymes that contain open-shell transition metal ions.

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.

Titration of ionizable groups in proteins using multiple neutron data sets from a single crystal: application to the small GTPase Ras

Neutron protein crystallography (NPC) reveals the three-dimensional structures of proteins, including the positions of H atoms. The technique is particularly suited to elucidate ambiguous catalytic steps in complex biochemical reactions. While NPC uniquely complements biochemical assays and X-ray structural analyses by revealing the protonation states of ionizable groups at and around the active site of enzymes, the technique suffers from a major drawback: large single crystals must be grown to compensate for the relatively low flux of neutron beams. However, in addition to revealing the positions of hydrogens involved in enzyme catalysis, NPC has the advantage over X-ray crystallography that the crystals do not suffer radiation damage. The lack of radiation damage can be exploited to conduct in crystallo parametric studies. Here, the use of a single crystal of the small GTPase Ras to collect three neutron data sets at pD 8.4, 9.0 and 9.4 is reported, enabling an in crystallo titration study using NPC. In addition to revealing the behavior of titratable groups in the active site, the data sets will allow the analysis of allosteric water-mediated communication networks across the molecule, particularly regarding Cys118 and three tyrosine residues central to these networks, Tyr32, Tyr96 and Tyr137, with pK_a values expected to be in the range sampled in our experiments.

Comparing the catalytic strategy of ATP hydrolysis in biomolecular motors

ATP-driven biomolecular motors utilize the chemical energy obtained from the ATP hydrolysis to perform vital tasks in living cells. Understanding the mechanism of enzyme-catalyzed ATP hydrolysis reaction has substantially progressed lately thanks to combined quantum/classical molecular mechanics (QM/MM) simulations. Here, we present a comparative summary of the most recent QM/MM results for myosin, kinesin and F1-ATPase motors. These completely different motors achieve the acceleration of ATP hydrolysis through a very similar catalytic mechanism. ATP hydrolysis has high activation energy because it involves the breaking of two strong bonds, namely the Pγ-Oβγ bond of ATP and the H-O bond of lytic water. The key to the four-fold decrease in the activation barrier by the three enzymes is that the breaking of the Pγ-Oβγ bond precedes the deprotonation of the lytic water molecule, generating a metaphosphate hydrate complex. The resulting singly charged trigonal planar PγO3(-) metaphosphate is a better electrophilic target for attack by an OaH(-) hydroxyl group. The formation of this OaH(-) is promoted by a strong polarization of the lytic water: in all three proteins, this water is forming a hydrogen-bond with a backbone carbonyl group and interacts with the carboxylate group of glutamate (either directly or via an intercalated water molecule). This favors the shedding of one proton by the attacking water. The abstracted proton is transferred to the γ-phosphate via various proton wires, resulting in a H2PγO4(-)/ADP(3-) product state. This catalytic strategy is so effective that most other nucleotide hydrolyzing enzymes adopt a similar approach, as suggested by their very similar triphosphate binding sites.

Kinesin motility is driven by subdomain dynamics

The microtubule (MT)-associated motor protein kinesin utilizes its conserved ATPase head to achieve diverse motility characteristics. Despite considerable knowledge about how its ATPase activity and MT binding are coupled to the motility cycle, the atomic mechanism of the core events remain to be found. To obtain insights into the mechanism, we performed 38.5 microseconds of all-atom molecular dynamics simulations of kinesin-MT complexes in different nucleotide states. Local subdomain dynamics were found to be essential for nucleotide processing. Catalytic water molecules are dynamically organized by the switch domains of the nucleotide binding pocket while ATP is torsionally strained. Hydrolysis products are 'pulled' by switch-I, and a new ATP is 'captured' by a concerted motion of the α0/L5/switch-I trio. The dynamic and wet kinesin-MT interface is tuned for rapid interactions while maintaining specificity. The proposed mechanism provides the flexibility necessary for walking in the crowded cellular environment.

Motor proteins called kinesins perform a number of different roles inside cells, including transporting cargo and organizing filaments called microtubules to generate the force needed for a cell to divide. Kinesins move along the microtubules, with different kinesins moving in different ways: some ‘walk’, some jump, and some destroy the microtubule as they travel along it. All kinesins power their movements using the same molecule as fuel – adenosine triphosphate, known as ATP for short.

Energy stored in ATP is released by a chemical reaction known as hydrolysis, which uses water to break off specific parts of the ATP molecule. The site to which ATP binds in a kinesin has a similar structure to the ATP binding site of many other proteins that use ATP. However, little was known about the way in which kinesin uses ATP as a fuel, including how ATP binds to kinesin and is hydrolyzed, and how the products of hydrolysis are released. These events are used to power the motor protein.

Hwang et al. have used powerful computer simulation methods to examine in detail how ATP interacts with kinesin whilst moving across a microtubule. The simulations suggest that regions (or 'domains') of kinesin near the ATP binding site move around to help in processing ATP. These kinesin domains trap a nearby ATP molecule from the environment and help to deliver water molecules to ATP for hydrolysis. Hwang et al. also found that the domain motion subsequently helps in the release of the hydrolysis products by kinesin.

The domains around the ATP pocket vary among the kinesins and these differences may enable kinesins to fine-tune how they use ATP to move. Further investigations will help us understand why different kinesin families behave differently. They will also contribute to exploring how kinesin inhibitors might be used as anti-cancer drugs.

Elucidation of Single Hydrogen Bonds in GTPases via Experimental and Theoretical Infrared Spectroscopy

Time-resolved Fourier transform infrared (FTIR) spectroscopy is a powerful tool to elucidate label-free the reaction mechanisms of proteins. After assignment of the absorption bands to individual groups of the protein, the order of events during the reaction mechanism can be monitored and rate constants can be obtained. Additionally, structural information is encoded into infrared spectra and can be decoded by combining the experimental data with biomolecular simulations. We have determined recently the infrared vibrations of GTP and guanosine diphosphate (GDP) bound to Gα_i1, a ubiquitous GTPase. These vibrations are highly sensitive for the environment of the phosphate groups and thereby for the binding mode the GTPase adopts to enable fast hydrolysis of GTP. In this study we calculated these infrared vibrations from biomolecular simulations to transfer the spectral information into a computational model that provides structural information far beyond crystal structure resolution. Conformational ensembles were generated using 15 snapshots of several 100 ns molecular-mechanics/molecular-dynamics (MM-MD) simulations, followed by quantum-mechanics/molecular-mechanics (QM/MM) minimization and normal mode analysis. In comparison with other approaches, no time-consuming QM/MM-MD simulation was necessary. We carefully benchmarked the simulation systems by deletion of single hydrogen bonds between the GTPase and GTP through several Gα_i1 point mutants. The missing hydrogen bonds lead to blue-shifts of the corresponding absorption bands. These band shifts for α-GTP (Gα_i1-T48A), γ-GTP (Gα_i1-R178S), and for both β-GTP/γ-GTP (Gα_i1-K46A, Gα_i1-D200E) were found in agreement in the experimental and the theoretical spectra. We applied our approach to open questions regarding Gα_i1: we show that the GDP state of Gα_i1 carries a Mg²⁺, which is not found in x-ray structures. Further, the catalytic role of K46, a central residue of the P-loop, and the protonation state of the GTP are elucidated.

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

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

Local Mode Analysis: Decoding IR Spectra by Visualizing Molecular Details

Integration of experimental and computational approaches to investigate chemical reactions in proteins has proven to be very successful. Experimentally, time-resolved FTIR difference-spectroscopy monitors chemical reactions at atomic detail. To decode detailed structural information encoded in IR spectra, QM/MM calculations are performed. Here, we present a novel method which we call local mode analysis (LMA) for calculating IR spectra and assigning spectral IR-bands on the basis of movements of nuclei and partial charges from just a single QM/MM trajectory. Through LMA the decoding of IR spectra no longer requires several simulations or optimizations. The novel approach correlates the motions of atoms of a single simulation with the corresponding IR bands and provides direct access to the structural information encoded in IR spectra. Either the contributions of a particular atom or atom group to the complete IR spectrum of the molecule are visualized, or an IR-band is selected to visualize the corresponding structural motions. Thus, LMA decodes the detailed information contained in IR spectra and provides an intuitive approach for structural biologists and biochemists. The unique feature of LMA is the bidirectional analysis connecting structural details to spectral features and vice versa spectral features to molecular motions.

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.

Structural role of two histidines in the (6-4) photolyase reaction

Photolyases (PHRs) are DNA repair enzymes that revert UV-induced photoproducts, either cyclobutane pyrimidine dimers (CPD) or (6-4) photoproducts (PPs), into normal bases to maintain genetic integrity. (6-4) PHR must catalyze not only covalent bond cleavage, but also hydroxyl or amino group transfer, yielding a more complex mechanism than that postulated for CPD PHR. Previous mutation analysis revealed the importance of two histidines in the active center, H354 and H358 for Xenopus (6-4) PHR, whose mutations significantly lowered the enzymatic activity. Based upon highly sensitive FTIR analysis of the repair function, here we report that both H354A and H358A mutants of Xenopus (6-4) PHR still maintain their repair activity, although the efficiency is much lower than that of the wild type. Similar difference FTIR spectra between the wild type and mutant proteins suggest a common mechanism of repair in which (6-4) PP binds to the active center of each mutant, and is released after repair, as occurs in the wild type. Similar FTIR spectra also suggest that a decrease in volume by the H-to-A mutation is possibly compensated by the addition of water molecule( s). Such a modified environment is sufficient for the repair function that is probably controlled by proton-coupled electron transfer between the enzyme and substrate. On the other hand, two histidines must work in a concerted manner in the active center of the wild-type enzyme, which significantly raises the repair efficiency.

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.

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.

Quantifying bond distortions in transient enzyme species by a combination of density functional theory calculations and time-resolved infrared difference spectroscopy. Implications for the mechanism of dephosphorylation of the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA1a)

The sarcoplasmic Ca(2+)-ATPase (SERCA1a) forms two phosphoenzyme intermediates during Ca(2+) pumping. The second intermediate E2P hydrolyzes rapidly, which is essential for the rapid removal of Ca(2+) from the cytosol of muscle cells. The present work studies whether a weakening of the scissile PO bond in the E2P ground state facilitates dephosphorylation. To this end, the experimentally known vibrational spectrum of the E2P phosphate group was calculated with density functional theory (DFT) using structural models at two levels of structural complexity: (i) Models of acetyl phosphate in simple environments and (ii) ~150 atom models of the catalytic site. It was found that the enzyme environment distorts the structure of the phosphate group: one of the terminal PO bonds is shorter in the catalytic site indicating weaker interactions than in water. However, the bond that bridges phosphate and Asp351 is unaffected. This indicates that the scissile PO bond is not weakened by the enzyme environment of E2P. A second finding was that the catalytic site of the E2P state in aqueous solution appears to adopt a structure as in the crystals with BeF3(-), where the ATPase is in a non-reactive conformation. The reactant state of the dephosphorylation reaction differs from the E2P ground state: Glu183 faces Asp351 and positions the attacking water molecule. This state has a 0.04Å longer, and thus weaker, bridging PO bond. The reactant state is not detected in our experiments, indicating that its energy is at least 1kcal/mol higher than that of the E2P ground state.

Asymmetric environmental effects on the structure and vibrations of cis-[Pt(NH3)2Cl2] in condensed phases

We report the structural and vibrational properties of anticancer drug cisplatin (cis-[Pt(NH3)2Cl2]) in gas phase, in solid phase, and in aqueous solution using density functional theory (DFT) calculations, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics, and effective normal modes analysis. In contrast with the gas-phase case, asymmetric hydrogen bonding environments are found in both solid phase and aqueous solution. It is shown that the discrepancy of the molecular geometry between previous gas phase calculations and the X-ray crystal structure can be resolved by considering intermolecular hydrogen bonds in the calculations of solid phase. In addition, our simulations in solid phase and aqueous solution reveal that asymmetric environmental effects lead to several spectral features observed in experiments, such as the blue-shift in the N-H stretching region and the frequency splitting of NH3 symmetric deformation modes. Furthermore, a similar decoupling and localization of several vibrational modes of cisplatin is found in solid phase and aqueous solution, in comparison to those of O-H stretching modes of water molecules in liquid water [ J. Phys. Chem. Lett. 2013 , 4 ( 19 ), 3245 - 3250 ].

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.

Lessons from computer simulations of Ras proteins in solution and in membrane

BACKGROUND

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

SCOPE OF REVIEW

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

MAJOR CONCLUSIONS

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

GENERAL SIGNIFICANCE

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.

Ras and GTPase-activating protein (GAP) drive GTP into a precatalytic state as revealed by combining FTIR and biomolecular simulations

Members of the Ras superfamily regulate many cellular processes. They are down-regulated by a GTPase reaction in which GTP is cleaved into GDP and P(i) by nucleophilic attack of a water molecule. Ras proteins accelerate GTP hydrolysis by a factor of 10(5) compared to GTP in water. GTPase-activating proteins (GAPs) accelerate hydrolysis by another factor of 10(5) compared to Ras alone. Oncogenic mutations in Ras and GAPs slow GTP hydrolysis and are a factor in many cancers. Here, we elucidate in detail how this remarkable catalysis is brought about. We refined the protein-bound GTP structure and protein-induced charge shifts within GTP beyond the current resolution of X-ray structural models by combining quantum mechanics and molecular mechanics simulations with time-resolved Fourier-transform infrared spectroscopy. The simulations were validated by comparing experimental and theoretical IR difference spectra. The reactant structure of GTP is destabilized by Ras via a conformational change from a staggered to an eclipsed position of the nonbridging oxygen atoms of the γ- relative to the β-phosphates and the further rotation of the nonbridging oxygen atoms of α- relative to the β- and γ-phosphates by GAP. Further, the γ-phosphate becomes more positive although two of its oxygen atoms remain negative. This facilitates the nucleophilic attack by the water oxygen at the phosphate and proton transfer to the oxygen. Detailed changes in geometry and charge distribution in the ligand below the resolution of X-ray structure analysis are important for catalysis. Such high resolution appears crucial for the understanding of enzyme catalysis.

The role of magnesium for geometry and charge in GTP hydrolysis, revealed by quantum mechanics/molecular mechanics simulations

The coordination of the magnesium ion in proteins by triphosphates plays an important role in catalytic hydrolysis of GTP or ATP, either in signal transduction or energy conversion. For example, in Ras the magnesium ion contributes to the catalysis of GTP hydrolysis. The cleavage of GTP to GDP and P(i) in Ras switches off cellular signaling. We analyzed GTP hydrolysis in water, Ras, and Ras·Ras-GTPase-activating protein using quantum mechanics/molecular mechanics simulations. By comparison of the theoretical IR-difference spectra for magnesium ion coordinated triphosphate to experimental ones, the simulations are validated. We elucidated thereby how the magnesium ion contributes to catalysis. It provides a temporary storage for the electrons taken from the triphosphate and it returns them after bond cleavage and P(i) release back to the diphosphate. Furthermore, the Ras·Mg(2+) complex forces the triphosphate into a stretched conformation in which the β- and γ-phosphates are coordinated in a bidentate manner. In this conformation, the triphosphate elongates the bond, which has to be cleaved during hydrolysis. Furthermore, the γ-phosphate adopts a more planar structure, driving the conformation of the molecule closer to the hydrolysis transition state. GTPase-activating protein enhances these changes in GTP conformation and charge distribution via the intruding arginine finger.

Influence of the molecular environment on phosphorylated amino acid models: a density functional theory study

A protein environment can affect the structure and charge distribution of substrate molecules. Here, the structure and partial charges were studied for different phosphorylated amino acid models in varying environments using density functional theory. The three systems investigated, acetyl phosphate, methyl phosphate, and p-tolyl phosphate are representative models for aspartyl phosphate, serine or threonine phosphate, and tyrosine phosphate, respectively. Combined with the CPCM continuum model, explicit HF and H(2)O molecules were added in order to model environmental effects and interactions that may occur in a protein matrix. We show how the different interactions affect the scissile P-O(R) bond and that the elongation can be explained by an anomeric effect. An increasing scissile bond length will result in transfer of negative charge to the leaving group and in a widening of the angle between the terminal oxygens of the phosphate molecule, features that can expose the phosphate group to attacking nucleophiles. Lastly, calculations were performed on the active site of the Ca(2+)-ATPase E2P intermediate, which provide an example of how a protein environment facilitates the formation of a destabilized ground state.

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.

FTIR study of light-dependent activation and DNA repair processes of (6-4) photolyase

The UV component of sunlight threatens all life on the earth by damaging DNA. The photolyase (PHR) DNA repair proteins maintain genetic integrity by harnessing blue light to restore intact bases from the major UV-induced photoproducts, cyclobutane pyrimidine dimers (CPD), and (6-4) photoproducts ((6-4) PPs). The (6-4) PHR must catalyze not only covalent bond cleavage between two pyrmidine bases but also hydroxyl or amino group transfer from the 5'- to 3'-pyrimidine base, requiring a more complex mechanism than that postulated for CPD PHR. In this paper, we apply Fourier transform infrared (FTIR) spectroscopy to (6-4) PHR and report difference FTIR spectra that correspond to its photoactivation, substrate binding, and light-dependent DNA repair processes. The presence of DNA carrying a single (6-4) PP uniquely influences vibrations of the protein backbone and a protonated carboxylic acid, whereas photoactivation produces IR spectral changes for the FAD cofactor and the surrounding protein. Difference FTIR spectra for the light-dependent DNA damage repair reaction directly show significant DNA structural changes in the (6-4) lesion and the neighboring phosphate group. Time-dependent illumination of samples with different enzyme:substrate stoichiometries successfully distinguished signals characteristic of structural changes in the protein and the DNA resulting from binding and catalysis.

New insight into the dynamic properties and the active site architecture of H-Ras p21 revealed by X-ray crystallography at very high resolution

Background

In kinetic crystallography, the usually static method of X-ray diffraction is expanded to allow time-resolved analysis of conformational rearrangements in protein structures. To achieve this, reactions have to be triggered within the protein crystals of interest, and optical spectroscopy can be used to monitor the reaction state. For this approach, a modified form of H-Ras p21 was designed which allows reaction initiation and fluorescence readout of the initiated GTPase reaction within the crystalline state. Rearrangements within the crystallized protein due to the progressing reaction and associated heterogeneity in the protein conformations have to be considered in the subsequent refinement processes.

Results

X-ray diffraction experiments on H-Ras p21 in different states along the reaction pathway provide detailed information about the kinetics and mechanism of the GTPase reaction. In addition, a very high data quality of up to 1.0 Å resolution allowed distinguishing two discrete subconformations of H-Ras p21, expanding the knowledge about the intrinsic flexibility of Ras-like proteins, which is important for their function. In a complex of H-Ras•GppNHp (guanosine-5'-(β,γ-imido)-triphosphate), a second Mg²⁺ ion was found to be coordinated to the γ-phosphate group of GppNHp, which positions the hydrolytically active water molecule very close to the attacked γ-phosphorous atom.

Conclusion

For the structural analysis of very high-resolution data we have used a new 'two-chain-isotropic-refinement' strategy. This refinement provides an alternative and easy to interpret strategy to reflect the conformational variability within crystal structures of biological macromolecules. The presented fluorescent form of H-Ras p21 will be advantageous for fluorescence studies on H-Ras p21 in which the use of fluorescent nucleotides is not feasible.

Unravelling the mechanism of dual-specificity GAPs

The molecular mechanism by which dual-specificity RasGAPs of the Gap1 subfamily activate the GTP hydrolysis of both Rap and Ras is an unresolved phenomenon. RasGAPs and RapGAPs use different strategies to stimulate the GTPase reaction of their cognate G-proteins. RasGAPs contribute an arginine finger to orient through the Gln61 of Ras the nucleophilic water molecule. RapGAP contributes an asparagine (Asn thumb) into the active site to substitute for the missing Gln61. Here, by using steady-state kinetic assays and time-resolved Fourier-transform infrared spectroscopy (FTIR) experiments with wild type and mutant proteins, we unravel the remarkable mechanism for the specificity switch. The plasticity of GAP1(IP4BP) and RASAL is mediated by the extra GTPase-activating protein (GAP) domains, which promote a different orientation of Ras and Rap's switch-II and catalytic residues in the active site. Thereby, Gln63 in Rap adopts the catalytic role normally taken by Gln61 of Ras. This re-orientation requires specific interactions between switch-II of Rap and helix-alpha6 of GAPs. This supports the notion that the specificities of fl proteins versus GAP domains are potentially different.

The hydrolysis activity of adenosine triphosphate in myosin: a theoretical analysis of anomeric effects and the nature of the transition state

Combined quantum mechanical/molecular mechanical (QM/MM) calculations with density functional theory are employed to analyze two issues related to the hydrolysis activity of adenosine triphosphate (ATP) in myosin. First, we compare the geometrical properties and electronic structure of ATP in the open (post-rigor) and closed (pre-powerstroke) active sites of the myosin motor domain. Compared to both solution and the open active site cases, the scissile P(gamma)-O(3beta) bond of ATP in the closed active site is shown to be substantially elongated. Natural bond orbital (NBO) analysis clearly shows that this structural feature is correlated with the stronger anomeric effects in the closed active site, which involve charge transfers from the lone pairs in the nonbridging oxygen in the gamma-phosphate to the antibonding orbital of the scissile bond. However, an energetic analysis finds that the ATP molecule is not significantly destabilized by the P(gamma)-O(3beta) bond elongation. Therefore, despite the notable perturbations in the geometry and electronic structure of ATP as its environment changes from solution to the hydrolysis-competent active site, ground-state destabilization is unlikely to play a major role in enhancing the hydrolysis activity in myosin. Second, two-dimensional potential energy maps are used to better characterize the energetic landscape near the hydrolysis transition state. The results indicate that the transition-state region is energetically flat and a range of structures representative of different mechanisms according to the classical nomenclature (e.g., "associative", "dissociative", and "concerted") are very close in energy. Therefore, at least in the case of ATP hydrolysis in myosin, the energetic distinction between different reaction mechanisms following the conventional nomenclature is likely small. This study highlights the importance of (i) explicitly evaluating the relevant energetic properties for determining whether a factor is essential to catalysis and (ii) broader explorations of the energy landscape beyond saddle points (even on free-energy surface) for characterizing the molecular mechanism of catalysis.

Role of the arginine finger in Ras.RasGAP revealed by QM/MM calculations

In the Ras.Ras.GAP complex, hydrolysis of guanosine triphosphate is strongly accelerated GAP as compared to Ras alone. This is largely attributed to the arginine finger R789(GAP) pointing to AlF(x) in the transition state analogue. We performed QM/MM simulations where triphosphate was treated using the quantum mechanical method of density functional theory, while the protein complex and water environment were described classically using MD. Compared to Ras, the crucial electron shift, bond stretching and distortion towards an eclipsed gamma-to-beta orientation are much more pronounced. The arginine finger is shown to act by displacing water out of the binding niche. The resulting enhanced electrostatic field catalyses the cleavage step.

Infrared spectroscopy of proteins

This review discusses the application of infrared spectroscopy to the study of proteins. The focus is on the mid-infrared spectral region and the study of protein reactions by reaction-induced infrared difference spectroscopy.

Development of effective quantum mechanical/molecular mechanical (QM/MM) methods for complex biological processes

Motivated by the long-term goal of understanding vectorial biological processes such as proton transport (PT) in biomolecular ion pumps, a number of developments were made to establish combined quantum mechanical/molecular mechanical (QM/MM) methods suitable for studying chemical reactions involving significant charge separation in the condensed phase. These developments were summarized and discussed with representative problems. Specifically, free energy perturbation and boundary potential methods for treating long-range electrostatics were implemented to test the robustness of QM/MM results for protein systems. It was shown that consistent models with sufficient sampling were able to produce quantitatively satisfactory results, such as pKa for titritable groups in the interior of T4-lysozyme, while an inconsistent treatment of electrostatics or lack of sufficient sampling may produce incorrect results. Modifications were made to an approximate density functional theory (SCC-DFTB) to improve the description of proton affinity and hydrogen-bonding, which are crucial for the treatment of PT in polar systems. Test calculations on water autoionization showed clearly that both improvements are necessary for quantitatively reliable results. Finally, the newly established SCC-DFTB/MM-GSBP protocol was used to explore mechanistic issues in carbonic anhydrase (CA). Preliminary results suggest that PT in CA occurs mainly through short water wires containing two water molecules in a thermally activated fashion. Although longer water wires occur with similar frequencies, PT along those pathways, on average, has substantially higher barriers, a result not expected based on previous studies. The fluctuations of water molecules peripheral to the water wire were found to make a larger impact on the PT energetics compared to polar protein residues in the active site, which are largely pre-organized and therefore have less tendency to reorganize during the reaction.

Surface change of Ras enabling effector binding monitored in real time at atomic resolution

Ras, the prototype of the Ras superfamily, acts as a molecular switch for cell growth. External growth signals induce a GDP-to-GTP exchange. This modifies the Ras surface (Ras(on)GTP) and enables effector binding, which then activates signal-transduction pathways. GTP hydrolysis, catalysed by Ras and GAP, returns the signal to "off" (Ras(off)GDP). Oncogenic mutations in Ras prevent this hydrolysis, and thereby cause uncontrolled cell growth. In the Ras(off)-to-Ras(on) transition, the Ras surface is changed by a movement of the switch I loop that controls effector binding. We monitored this surface change at atomic resolution in real time by time-resolved FTIR (trFTIR) spectroscopy. In the transition from Ras(off) to Ras(on) a GTP-bound intermediate is now identified, in which effector binding is still prevented (Ras(off)GTP). The loop movement from Ras(off)GTP to Ras(on)GTP was directly monitored by the C=O vibration of Thr35. The structural change creates a binding site with a rate constant of 5 s(-1) at 260 K. A small molecule that shifted the equilibrium from the Ras(on)GTP state towards the Ras(off)GTP state would prevent effector binding, even if hydrolysis were blocked by oncogenic mutations. We present a spectroscopic fingerprint of both states that can be used as an assay in drug screening for such small molecules.

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.

Structure and Raman Spectrum of Clavulanic Acid in Aqueous Solution

The calculation of the vibrational Raman spectrum of enzyme-bound beta-lactamase inhibitors may be of help to understand the mechanisms responsible for bacterial drug resistance. Here, we present a study of the solvation structure and the vibrational properties of clavulanate, an important beta-lactamase inhibitor, in aqueous solution as obtained from full quantum and hybrid empirical/quantum molecular dynamics simulations at ambient conditions. The analysis of the vibrational density of states indicates that hybrid empirical/quantum mechanical simulations are able to properly describe the vibrational levels of clavulanate in solution. In addition, we propose a computationally efficient protocol to calculate the vibrational Raman effect for large solute molecules in water, which is able to faithfully reproduce the experimentally recorded clavulanate Raman spectrum and discloses the possibility to employ hybrid simulations to assign the experimental Raman spectra of inhibitors bound to beta-lactamases.

On the mechanism of hydrolysis of phosphate monoesters dianions in solutions and proteins

The nature of the hydrolysis of phosphate monoester dianions in solutions and in proteins is a problem of significant current interest. The present work explores this problem by systematic calculations of the potential surfaces of the reactions of a series of phosphate monoesters with different leaving groups. These calculations involve computational studies ranging from ab initio calculations with implicit solvent models to ab initio QM/MM free energy calculations. The calculations reproduce the observed linear free energy relationship (LFER) for the solution reaction and thus are consistent with the overall experimental trend and can be used to explore the nature of the transition state (TS) region, which is not accessible to direct experimental studies. It is found that the potential surface for the associative and dissociative paths is very flat and that the relative height of the associative and dissociative TS is different in different systems. In general, the character of the TS changes from associative to dissociative upon decrease in the pKa of the leaving group. It is also demonstrated that traditional experimental markers such as isotope effects and the LFER slope cannot be used in a conclusive way to distinguish between the two classes of transition states. In addition it is found that the effective charges of the TS do not follow the previously assumed simple rule. Armed with that experience we explore the free energy surface for the GTPase reaction of the RasGap system. In this case it is found that the surface is flat but that the lowest TS is associative. The present study indicates that the nature of the potential surfaces for the phosphoryl transfer reactions in solution and proteins is quite complicated and cannot be determined in a conclusive way without the use of careful theoretical studies that should, of course, reproduce the available experimental information.

Structural basis of effector regulation and signal termination in heterotrimeric Galpha proteins

This chapter addresses, from a molecular structural perspective gained from examination of x-ray crystallographic and biochemical data, the mechanisms by which GTP-bound Galpha subunits of heterotrimeric G proteins recognize and regulate effectors. The mechanism of GTP hydrolysis by Galpha and rate acceleration by GAPs are also considered. The effector recognition site in all Galpha homologues is formed almost entirely of the residues extending from the C-terminal half of alpha2 (Switch II) together with the alpha3 helix and its junction with the beta5 strand. Effector binding does not induce substantial changes in the structure of Galpha*GTP. Effectors are structurally diverse. Different effectors may recognize distinct subsets of effector-binding residues of the same Galpha protein. Specificity may also be conferred by differences in the main chain conformation of effector-binding regions of Galpha subunits. Several Galpha regulatory mechanisms are operative. In the regulation of GMP phospodiesterase, Galphat sequesters an inhibitory subunit. Galphas is an allosteric activator and inhibitor of adenylyl cyclase, and Galphai is an allosteric inhibitor. Galphaq does not appear to regulate GRK, but is rather sequestered by it. GTP hydrolysis terminates the signaling state of Galpha. The binding energy of GTP that is used to stabilize the Galpha:effector complex is dissipated in this reaction. Chemical steps of GTP hydrolysis, specifically, formation of a dissociative transition state, is rate limiting in Ras, a model G protein GTPase, even in the presence of a GAP; however, the energy of enzyme reorganization to produce a catalytically active conformation appears to be substantial. It is possible that the collapse of the switch regions, associated with Galpha deactivation, also encounters a kinetic barrier, and is coupled to product (Pi) release or an event preceding formation of the GDP*Pi complex. Evidence for a catalytic intermediate, possibly metaphosphate, is discussed. Galpha GAPs, whether exogenous proteins or effector-linked domains, bind to a discrete locus of Galpha that is composed of Switch I and the N-terminus of Switch II. This site is immediately adjacent to, but does not substantially overlap, the Galpha effector binding site. Interactions of effectors and exogenous GAPs with Galpha proteins can be synergistic or antagonistic, mediated by allosteric interactions among the three molecules. Unlike GAPs for small GTPases, Galpha GAPs supply no catalytic residues, but rather appear to reduce the activation energy for catalytic activation of the Galpha catalytic site.

A phosphoryl transfer intermediate in the GTPase reaction of Ras in complex with its GTPase-activating protein

The hydrolysis of nucleoside triphosphates by enzymes is used as a regulation mechanism in key biological processes. Here, the GTP hydrolysis of the protein complex of Ras with its GTPase-activating protein is monitored at atomic resolution in a noncrystalline state by time-resolved FTIR spectroscopy. At 900 ms, after the attack of water at the gamma-phosphate, there appears a H2PO4- intermediate that is shown to be hydrogen-bonded in an eclipsed conformation to the beta-phosphate of GDP. The H2PO4- intermediate is in a position where it can either reform GTP or be released from the protein in 7 s in the rate-limiting step of the GTPase reaction. We propose that such an intermediate also occurs in other GTPases and ATPases.

A polarizable continuum approach for the study of heterogeneous dielectric environments

We present a computational method, exploiting some features of the polarizable continuum model (PCM) to describe heterogeneous media; it belongs to the family of electrostatic embedding mixed methods, such as the more common quantum-mechanical (QM)/molecular mechanics approaches, with the electrostatic long range effects accounted for by a polarized continuum instead of atomic point charges. Provided effective dielectric constants are determined for the various parts of the system, the method is much faster than its atomistic counterpart, and allows for high-level QM calculations on the fragment of interest, using all the highly efficient computational tools developed for homogeneous PCM. Two case studies (the calculation of the pKa of solvent exposed acidic residues in a model protein, and the calculation of the electron spin resonance spectrum of a typical spin probe partially embedded in a membrane) are analyzed in some detail, to illustrate the application of the method to complex systems.

Empirical rules facilitate the search for binding sites on protein surfaces

Computational surface screening of 3D protein structures is a valuable means of finding possible docking sites for substrates, effectors and similar molecules. It can be improved by considering properties of molecules which are known to bind to protein surfaces, and thus reflect the required properties of binding sites. In-depth studies are available on drugs and lead compounds as binding partners with statistically assured properties. Here we present a simple strategy for finding binding sites, which is based on the empirical rule-of-five by Lipinski et al. for oral drugs and the rule-of-three by Congreve et al. for leads. The fast automated search with the new C-code TRIDOCK yields a preliminary set of sites, thus facilitating further investigation by visual, comparative and quantitative work. Possible binding sites are tagged by pseudo-atoms added to the structure file for molecular graphical evaluation. Usually, the strategy yields not just a few single sites, but an accumulation of several sites in known substrate binding pockets. Clusters are also found at known or putative protein-protein docking interfaces. A comparison of the activated and inactivated form of the GTPase Ras reveals clear differences and identifies a niche, which is possibly a suitable new target for compounds that bind specifically to activated Ras.