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.

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.

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.

Exploration of the zinc finger motif in controlling activity of matrix metalloproteinases.

Discovering ways to control the activity of matrix metalloproteinases (MMPs), zinc-dependent enzymes capable of degrading extracellular matrix proteins, is an important field of cancer research. We report here a novel strategy for assembling MMP inhibitors on the basis of oligopeptide ligands by exploring the pattern known as the zinc finger motif. Advanced molecular modeling tools were used to characterize the structural binding motifs of experimentally tested MMP inhibitors, as well as those of newly proposed peptidomimetics, in their zinc-containing active sites. The results of simulations based on the quantum mechanics/molecular mechanics (QM/MM) approach and Car-Parrinello molecular dynamics with QM/MM potentials demonstrate that, upon binding of Regasepin1, a known MMP-9 inhibitor, the Zn(2+)(His3) structural element is rearranged to the Zn(2+)(Cys2His2) zinc finger motif, in which two Cys residues are borrowed from the ligand. Following consideration of the crystal structure of MMP-2 with its inhibitor, the oligopeptide APP-IP, we proposed a new peptidomimetic with two replacements in the substrate, Tyr3Cys and Asp6Cys. Simulations show that this peptide variant blocks an enzyme active site by the Zn(2+)(Cys2His2) zinc finger construct. Similarly, a natural substrate of MMP-2, Ace-Gln-Gly ∼ Ile-Ala-Gly-Nme, can be converted to an inhibiting compound by two replacements, Ile by Cys and Gly by the d isomer of Cys, favoring formation of the zinc finger motif.

Metal ions-stimulated iron oxidation in hydroxylases facilitates stabilization of HIF-1 alpha protein.

The exposure of cells to several metal ions stabilizes HIF-1 alpha protein. However, the molecular mechanisms are not completely understood. They may involve inhibition of hydroxylation by either substitution of iron by metal ions or by iron oxidation in the hydroxylases. Here we provide evidence supporting the latter mechanism. We show that HIF-1 alpha stabilization in human lung epithelial cells occurred following exposure to various metal and metalloid ions, including those that cannot substitute for iron in the hydroxylases. In each case addition of the reducing agent ascorbic acid (AA)* abolished HIF-1 alpha protein stabilization. To better understand the role of iron oxidation in hydroxylase inhibition and to define the role of AA in the enzyme recovery we applied molecular modeling techniques. Our results indicate that the energy required for iron substitution by Ni(II) in the enzyme is high and unlikely to be achieved in a biological system. Additionally, computer modeling allowed us to identify a tridentate coordination of AA with the enzyme-bound iron, which explains the specific demand for AA as the iron reductant. Thus, the stabilization of HIF-1 alpha by numerous metal ions that cannot substitute for iron in the enzyme, the alleviation of this effect by AA, and our computer modeling data support the hypothesis of iron oxidation in the hydroxylases following exposure to metal ions.

Mechanism of the myosin catalyzed hydrolysis of ATP as rationalized by molecular modeling.

The intrinsic chemical reaction of adenosine triphosphate (ATP) hydrolysis catalyzed by myosin is modeled by using a combined quantum mechanics and molecular mechanics (QM/MM) methodology that achieves a near ab initio representation of the entire model. Starting with coordinates derived from the heavy atoms of the crystal structure (Protein Data Bank ID code 1VOM) in which myosin is bound to the ATP analog ADP.VO(4)(-), a minimum-energy path is found for the transformation ATP + H(2)O --> ADP + P(i) that is characterized by two distinct events: (i) a low activation-energy cleavage of the P(gamma) O(betagamma) bond and separation of the gamma-phosphate from ADP and (ii) the formation of the inorganic phosphate as a consequence of proton transfers mediated by two water molecules and assisted by the Glu-459-Arg-238 salt bridge of the protein. The minimum-energy model of the enzyme-substrate complex features a stable hydrogen-bonding network in which the lytic water is positioned favorably for a nucleophilic attack of the ATP gamma-phosphate and for the transfer of a proton to stably bound second water. In addition, the P(gamma) O(betagamma) bond has become significantly longer than in the unbound state of the ATP and thus is predisposed to cleavage. The modeled transformation is viewed as the part of the overall hydrolysis reaction occurring in the closed enzyme pocket after ATP is bound tightly to myosin and before conformational changes preceding release of inorganic phosphate.

Electronic Excitations of the Chromophore from the Fluorescent Protein asFP595 in Solutions.

We present the results of modeling spectral properties of the chromophore, 2-acetyl-4-(p-hydroxybenzylidene)-1-methyl-5-imidazolone (AHBMI), from the newly discovered fluorescent protein asFP595 in different solvents and compare computational and recent experimental data. The time-dependent density functional theory (TDDFT) method is used to estimate positions of spectral bands with large oscillator strengths for vertical transitions to excited states following geometry optimizations of chromophore coordinates in vacuo and in solutions. The performance of different TDDFT functionals in computing excitations for a simpler chromophore from the green fluorescent protein was tested at the preliminary stage. Properties of various protonation states (neutral, anionic, zwitterionic) for the cis and trans conformations of AHBMI are compared. By using the polarizable continuum model, the following solvents have been considered for AHBMI: water, ethanol, acetonitrile, and dimethyl sulfoxide. It is shown that the bands found experimentally in aqueous solution refer to the cis neutral and cis anionic (or trans zwitterionic) conformations. The computed band positions deviate from experimental ones in water by no more than 35 nm (0.23 eV). In accord with experimental studies, the band shifts in different solvents do not show correlation with the dielectric constant or dipole moment; however, the computed values of the shifts are much smaller than those measured experimentally for the ionic species.

Quantum chemical modeling of reaction mechanism for 2-oxoglutarate dependent enzymes: effect of substitution of iron by nickel and cobalt.

Enzymatic hydroxylation reactions carried out by 2-oxoglutarate (2OG) dependent iron-containing oxygenases were recently implicated in oxygen sensing. In addition to oxygen depletion, two metals, cobalt and nickel, are capable of inducing hypoxic stress in cells by inhibiting oxygenase activity. Two possible scenarios have been proposed for the explanation of the hypoxic effects of cobalt and nickel: oxidation of enzyme-bound iron following cobalt or nickel exposure, and substitution of iron by cobalt or nickel. Here, by using density functional theory calculations, we modeled the reaction route from the reaction components to the high-spin metal-oxide intermediate in the activation of oxygen molecule by 2OG-dependent enzymes for three metal ions Fe(II), Ni(II), and Co(II) in the active site. An initial molecular model was constructed based on the crystal structure of iron-containing asparaginyl hydroxylase (FIH-1). Nickel- and cobalt-containing enzymes were modeled by a consequent replacement of the iron in the active center. The energy profiles connecting stationary points on the potential surfaces were computed by using the intrinsic reaction coordinate (IRC) technique from the located transition states. The results of calculations show that the substitution of iron by nickel or cobalt modifies the reaction energy profile; however, qualitatively, the reaction mechanism remains essentially the same. Thus, we would postulate that if the iron ion in the active site were substitutable by nickel and/or cobalt ions enzyme activity would be considerably altered due to high activation barriers.

Computational study of a transition state analog of phosphoryl transfer in the Ras-RasGAP complex: AlF(x) versus MgF3-.

The structures of the complexes between Ras*GDP bound to RasGAP in the presence of three probable gamma-phosphate analogs (AlF3, AlF4- and MgF3-) for the transition state (TS) of the hydrolysis of guanosine triphosphate (GTP) by the Ras-RasGAP enzymes have been modeled by quantum mechanical-molecular mechanical (QM/MM) calculations. These simulations contribute to the dispute on the nature of the TS in the hydrolysis reaction, since medium resolution X-ray crystallography cannot discern among stereochemically similar isoelectronic species (e.g., AlF3 or MgF3-). The optimized geometry for each structure has been found starting from experimental coordinates of one of them (PDBID: 1WQ1). Direct comparison of the experimental and computed geometry configurations in the immediate vicinity of the active site suggests that MgF3- is the most likely candidate for the phosphate analog in the experimental structure.

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF, MP2, DFT, and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory.

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.

Theoretical studies on the hydrolysis of mono-phosphate and tri-phosphate in gas phase and aqueous solution.

Phosphate hydrolysis by GTPases plays an important role as a molecular switch in signal transduction and as an initiator of many other biological processes. Despite the centrality of this ubiquitous reaction, the mechanism is still poorly understood. As a first step to understand the mechanisms of this process, the nonenzymatic hydrolysis of mono-phosphate and tri-phosphate esters were systematically studied in gas phase and aqueous solution using hybrid density functional methods. The dielectric effect of the environment on the energetics of these processes was also explored. Theoretical results show that for mono-phosphate ester, the dissociative pathway is much more favorable than the associative pathway. However, the reaction barriers for the dissociative and associative pathways of tri-phosphate hydrolysis are very close in aqueous solution, though the dissociative pathway is more favorable in the gas phase. High dielectric solvents, such as water, significantly lower the activation barrier of the associative pathway due to the greater solvation energy of the associative transition states than that of the reactant complex. By contrast, the barrier of the dissociative pathway, with respect to the gas phase, is less sensitive to the surrounding dielectric. In the associative hydrolysis pathway of the tri-phosphate ester, negative charge is transferred from the gamma-phosphate to beta-phosphate through the bridging ester oxygen and results in Pgamma-O bond dissociation. No analogous charge transfer was observed in the dissociative pathway, where Pgamma-O bond dissociation resulted from proton transfer from the gamma-phosphate to the bridge oxygen. Finally, the active participation of local water molecules can significantly lower the activation energy of the dissociative pathway for both mono-phosphate and tri-phosphate.