Role of glutamine-61 in the hydrolysis of GTP by p21H-ras: an experimental and theoretical study.

The active GTP-bound form of p21ras is converted to the biologically inactive GDP-bound form by enzymatic hydrolysis and this function serves to regulate the wild-type ras protein. The side chain of the amino acid at position 61 may play a key role in this hydrolysis of GTP by p21. Experimental studies that define properties of the Q61E mutant of p21H-ras are presented along with supporting molecular dynamics simulations. We find that under saturating concentrations of GTP the Q61E mutant of p21H-ras has a 20-fold greater rate of intrinsic hydrolysis (kcat = 0.57 min-1) than the wild type. The affinity of the Q61E variant for GTP (Kd = 115 microM) is much lower than that of the wild type. GTPase activating protein does not activate the variant. From molecular dynamics simulations, we find that both the wild type and Q61E mutant have the residue 61 side chain in transient contact with a water molecule that is well-positioned for hydrolytic attack on the gamma phosphate. Thr-35 also is found to form a transient hydrogen bond with this critical water. These elements may define the catalytic complex for hydrolysis of the GTP [Pai et al. (1990) EMBO J. 9, 2351]. Similarly, the G12P mutant, which also has an intrinsic hydrolysis rate similar to the wild type, is found to form the same complex in simulation. In contrast, molecular dynamics analysis of the mutants G12R, G12V, and Q61L, which have much lower intrinsic rates than the wild-type p21, do not show this complex.(ABSTRACT TRUNCATED AT 250 WORDS)

Simulation of the solution structure of the H-ras p21-GTP complex.

An unconstrained simulation of the GTP-bound form of the H-ras protein p21 is performed in an aqueous environment with charge-neutralizing counterions. The simulation is compared to the 1.35-A structure of Pai et al. [(1990) EMBO J. 9, 2351] and a proposed alternate structure, in which the loop at residues 60-65 is modeled into a form which may activate a water molecule for the GTP hydrolysis. The simulation suggests that some protein intermolecular H-bond contacts which are present in the crystal structure are lost in the solvation process and this loss may lead to localized refolding of the molecule. For instance, we find that the gamma-phosphate of the GTP has somewhat weaker contact with the protein in the simulation structure. The antiparallel beta-sheet (residues 38-57) partially melts. The 60-65 loop, which is hypervariable in the X-ray study, is initially relatively distant from the gamma-phosphate region. However, this loop moves so as to sample the space around the gamma-phosphate. For a significant fraction of the simulation time, forms similar to the alternate structure are observed, and a water molecule is localized near the hydrolytic site. The molecular dynamics simulations of p21-GTP in solution support a postulated hydrolysis mechanism for the biological inactivation of the nucleotide complex based on crystallographic data.

Theoretical and experimental measures of DNA helix stability and their relation to sequence specific repair of O6-ethylguanine lesions.

Recent work (Breslauer et al. (1986) Proc. Natl. Acad. Sci. (U.S.A.), 83, 3746) has provided a method for calculating empirical thermodynamic quantities for helix to coil transitions from the base sequence of any oligomer. It is shown in this work that the DNA helix binding energy, calculated with the AMBER force field, for 9-mers of the type 5'-GGGXGeYGGG-3', where X and Y are any base and the central Ge is O6-ethylguanine, correlates well with the empirical delta G for helix to strand transitions. The mutation spectrum of ethane methylsulfonate (EMS) in the lacI gene of Escherichia coli can be modeled using the calculated local binding energy but the empirical free energies, enthalpies and melting temperatures predict these levels of repair less well. The relation of the binding energy to the mutation spectrum can be somewhat improved by including entropic effects in a theoretical free energy of binding as given by delta G theoretical identical to delta E binding - T delta S.