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.

Elongation factor Tu (EF-Tu), the protein responsible for delivering aminoacyl-tRNAs (aa-tRNAs) to ribosomal A site during translation, belongs to the group of guanosine-nucleotide (GTP/GDP) binding proteins. Its active 'on'-state corresponds to the GTP-bound form, while the inactive 'off'-state corresponds to the GDP-bound form. In this work we focus on the chemical step, GTP + H 2 O → GDP + Pi, of the hydrolysis mechanism. We apply molecular modeling tools including molecular dynamics simulations and the combined quantum mechanical -molecular mechanical calculations for estimates of reaction energy profiles for two possible arrangements of switch II regions of EFTu. In the first case we presumably mimic binding of the ternary complex EF-Tu·GTP·aatRNA to the ribosome and allow the histidine (His85) side chain of the protein to approach the reaction active site. In the second case, corresponding to the GTP hydrolysis by EF-Tu alone, the side chain of His85 stays away from the active site, and the chemical reaction GTP + H 2 O → GDP + Pi proceeds without participation of the histidine but through water molecules. In agreement with the experimental observations which distinguish rate constants for the fast chemical reaction in EFTu·GTP·aatRNA·ribosome and the slow spontaneous GTP hydrolysis in EF-Tu, we show that the activation energy barrier for the first scenario is considerably lower compared to that of the second case.

Atomistic molecular dynamics simulations of diffusion of O2 from the hemes to the external solvent in the α- and β-subunits of the human hemoglobin (HbA) tetramer reveal transient gas tunnels that are not seen in crystal structures. We find here that the tunnel topology, which encompasses the reported experimental Xe binding cavities, is identical in HbA's T, R, and R2 quaternary states. However, the O2 population in the cavities and the preferred O2 escape portals vary significantly with quaternary structure. For example, most O2 molecules escape from the T β-subunit via the cavity at the center of the tetramer, but direct exit from the distal heme pocket dominates in the R2 β-subunit. To understand what triggers the quaternary-linked redistribution of O2 within its tunnels, we examined how the simulated tertiary structure and dynamics of each subunit differs among T, R, and R2 and report that minor adjustments in α-chain dynamics and β-heme position modulate O2 distribution and escape in HbA. Coupled to the β-heme position, residue βF71 undergoes quaternary-linked conformations that strongly regulate O2 migration between the β-subunit and HbA's central cavity. Remarkably, the distal histidine (HisE7) remains in a closed conformation near the α- and β-hemes in all states, but this does not prevent an average of 23, 31, and 46% of O2 escapes from the distal heme pockets of T, R, and R2, respectively, via several distal portals, with the balance of escapes occurring via the interior tunnels. Furthermore, preventing or restricting the access of O2 to selected cavities by mutating HisE7 and other heme pocket residues to tryptophan reveals how O2 migration adjusts to the bulky indole ring and sheds light on the experimental ligand binding kinetics of these variants. Overall, our simulations underscore the high gas porosity of HbA in its T, R, and R2 quaternary states and provide new mechanistic insights into why undergoing transitions among these states likely ensures effective O2 delivery by this tetrameric protein.

Hemoglobin transports O2 by binding the gas at its four hemes. Hydrogen bonding between the distal histidine (HisE7) and heme-bound O2 significantly increases the affinity of human hemoglobin (HbA) for this ligand. HisE7 is also proposed to regulate the release of O2 to the solvent via a transient E7 channel. To reveal the O2 escape routes controlled by HisE7 and to evaluate its role in gating heme access, we compare simulations of O2 diffusion from the distal heme pockets of the T and R states of HbA performed with HisE7 in its open (protonated) and closed (neutral) conformations. Irrespective of HisE7's conformation, we observe the same four or five escape routes leading directly from the α- or β-distal heme pockets to the solvent. Only 21-53% of O2 escapes occur via these routes, with the remainder escaping through routes that encompass multiple internal cavities in HbA. The conformation of the distal HisE7 controls the escape of O2 from the heme by altering the distal pocket architecture in a pH-dependent manner, not by gating the E7 channel. Removal of the HisE7 side chain in the GlyE7 variant exposes the distal pockets to the solvent, and the percentage of O2 escapes to the solvent directly from the α- or β-distal pockets of the mutant increases to 70-88%. In contrast to O2, the dominant water route from the bulk solvent is gated by HisE7 because protonation and opening of this residue dramatically increase the rate of influx of water into the empty distal heme pockets. The occupancy of the distal heme site by a water molecule, which functions as an additional nonprotein barrier to binding of the ligand to the heme, is also controlled by HisE7. Overall, analysis of gas and water diffusion routes in the subunits of HbA and its GlyE7 variant sheds light on the contribution of distal HisE7 in controlling polar and nonpolar ligand movement between the solvent and the hemes.

Tubulins are an ancient family of eukaryotic proteins characterized by an amino-terminal globular domain and disordered carboxyl terminus. These carboxyl termini play important roles in modulating the behavior of microtubules in living cells. However, the atomic-level basis of their function is not well understood. These regions contain multiple acidic residues and their overall charges are modulated in vivo by post-translational modifications, for example, phosphorylation. In this study, we describe an application of NMR and computer Monte Carlo simulations to investigate how the modification of local charge alters the conformational sampling of the γ-tubulin carboxyl terminus. We compared the dynamics of two 39-residue polypeptides corresponding to the carboxyl-terminus of yeast γ-tubulin. One polypeptide comprised the wild-type amino acid sequence while the second contained a Y > D mutation at Y11 in the polypeptide (Y445 in the full protein). This mutation introduces additional negative charge at a site that is phosphorylated in vivo and produces a phenotype with perturbed microtubule function. NMR relaxation measurements show that the Y11D mutation produces dramatic changes in the millisecond-timescale motions of the entire polypeptide. This observation is supported by Monte Carlo simulations that-similar to NMR-predict the WT γ-CT is largely unstructured and that the substitution of Tyr 11 with Asp causes the sampling of extended conformations that are unique to the Y11D polypeptide.