Molecular model for messenger RNA splicing.

A molecular model is presented for a messenger RNA (mRNA) "splice region." The model requires cation coordination to reduce backbone-backbone electrostatic repulsion and it allows for every base residue on the pre-mRNA to be stacked in A-form helical geometry with a recognition element on the intron or exon (or both) sides of the splice junction. The two nucleotides involved in the initial steps of the cleavage-ligation mechanism must adopt a non-A-form geometry, which ideally positions reactive groups on the pre-mRNA for the necessary catalytic chemistry. The model is also consistent with available biochemical data on splicing reactions.

Controlling the regiospecificity and coupling of cytochrome P450cam: T185F mutant increases coupling and abolishes 3-hydroxynorcamphor product.

Cytochrome P450cam (P450CIA1) catalyzes the hydroxylation of camphor and several substrate analogues such as norcamphor and 1-methyl-norcamphor. Hydroxylation was found experimentally at the 3, 5, and 6 positions of norcamphor, but only at the 5 and 6 positions of 1-methyl-norcamphor. In the catalytic cycle, the hydroxylation of substrate is coupled to the consumption of NADH. For camphor, the degree of coupling is 100%, but for both norcamphor and 1-methyl-norcamphor, the efficiency is dramatically lowered to 12% and 50%, respectively. Based on an examination of the active site of P450cam, it appeared that mutating position 185 might dramatically alter the product specificity and coupling of hydroxylation of norcamphor by P450cam. Analysis of molecular dynamics trajectories of norcamphor bound to the T185F mutant of cytochrome P450cam predicted that hydroxylation at the 3 position should be abolished and that the coupling should be dramatically increased. This mutant was constructed and the product profile and coupling experimentally determined. The coupling was doubled, and hydroxylation at the 3 position was essentially abolished. Both of these results are in agreement with the prediction.

Secondary structure in solution of two anti-HIV-1 hammerhead ribozymes as investigated by two-dimensional 1H 500 MHz NMR spectroscopy in water.

Two hammerhead chimeric RNA/DNA ribozymes (HRz) were synthesized in pure form. Both were 30 nucleotides long, and the sequences were such that they could be targeted to cleave the HIV-1 gag RNA. Named HRz-W and HRz-M, the former had its invariable core region conserved, the latter had a uridine in the invariable region replaced by a guanine. Their secodary structures were determined by 2D NOESY 1H 500 MHz NMR spectroscopy in 90% water and 10% D2(0), following the imino protons. The data show that both HRz-M and HRz-W form identical secondary structures with stem regions consisting of continuous stacks of AT and GT pairs. An energy minimized computer model of this stem region is provided. The results suggest that the loss of catalytic activity that is known to result when an invariant core residue is replaced is not related to the secondary structure of the ribozymes in the absence of substrate.

Rational approach to improving reductive catalysis by cytochrome P450cam.

Although halogenated hydrocarbons are noted for low chemical reactivity, small amounts are toxic to humans. Cytochromes P450 have been implicated in transforming these compounds to more reactive species, under anaerobic conditions, through reduction at the heme. A significant amount of effort has been directed toward turning this catalytic ability to our advantage by engineering P450 variants than can efficiently remediate these compounds in situ, before they come in contact with the human population. We have taken a 'rational' approach to this problem, in which a combination of theory and molecular modeling is applied to identify which properties of the enzyme have the greatest influence over reductive dehalogenation. Recent progress in this area is briefly reviewed. Two novel mutants, incorporating tryptophan (positions 87 and 396) and histidine (position 96, neutral and protonated) amino acid substitutions in the active site, are proposed and evaluated using molecular dynamics simulations. The upper bound on rate enhancement relative to wild-type is estimated in each mutant using electron transfer theory. The most significant rate enhancement is predicted for the His 96 mutant in the protonated state; while some His residues of certain proteins exhibit a pKa high enough to support a large protonated population, such information is not presently available for this proposed mutant.

Using molecular dynamics simulations on crambin to evaluate the suitability of different continuum dielectric and hydrogen atom models for protein simulations.

Molecular dynamics simulations of enzymes with enough explicit waters of solvation to realistically account for solute-solvent interactions can burden the computational resources required to perform the simulation by more than two orders of magnitude. Since enzyme simulations even with an implicit solvation model can be imposing for a supercomputer, it is important to assess the suitability of different continuum dielectric models for protein simulations. A series of 100-picosecond molecular dynamics simulations were performed on the X-ray crystal structure of the protein crambin to examine how well computed structures, obtained using seven continuum dielectric and two hydrogen atom models, agreed with the X-ray structure. The best level of agreement between computed and experimental structures was obtained using a constant dielectric of 2 and the all-hydrogen model. Continuum dielectric models of 1, 1r, and 2r also led to computed structures in reasonably good agreement with the X-ray structure. In all cases, the all-hydrogen model gave better agreement than the united atom model, although, in one case, the difference was not significant. Dielectric models of 4, 80, and 4r with either hydrogen model yielded significantly poorer fits. It is especially noteworthy that the observed trends did not semiquantitatively converge until about 50 picoseconds into the simulations, suggesting that validation studies for protein calculations based on energy minimizations or short simulations should be viewed with caution.

Ab initio quantum mechanics analysis of imidazole C-H...O water hydrogen bonding and a molecular mechanics forcefield correction.

While it is well established that classical hydrogen bonds play an important role in enzyme structure, function and dynamics, the role of weaker, but 'activated' C-H donor hydrogen bonds is poorly understood. The most important such case involves histidine which often plays a direct role in enzyme catalysis and possesses the most acidic C-H donor group of the standard amino acids. In the present study, we obtained optimized geometries and hydrogen bond interaction energies for C-H...O hydrogen bonded complexes between methane, ethylene, benzene, acetylene, and imidazole with water at the MP2-FC/6-31++G(2d,2p) and MP2-FC/aug-cc-pVDZ/MP2-FC/6-31++G(2d,2p) levels of theory. A strong linear relationship is obtained between the stability of the various hydrogen bonded complexes and both separation distances for H...O and C----O. In general, these calculations indicate that C-H...O interactions can be classified as hydrogen bonding interactions, albeit significantly weaker than the classical hydrogen bonds, but significantly stronger than just van der Waals interactions. For instance, while the electronic energy of stabilization at the MP2-FC/aug-cc-pVDZ/MP2-FC/6-31++G(2d,2p) level of theory of a water O-H...O water hydrogen bond is 4.36 kcal/mol more stable than the methane C-H...O water interaction, the water-water hydrogen bond is only 2.06 kcal/mol more stable than the imidazole Ce-H...O water hydrogen bond. Neglecting this latter hydrogen bonding interaction is obviously unacceptable. We next compare the potential energy surfaces for the imidazole Ce-H...O water and imidazole Na-H...O hydrogen bonded complexes computed at the MP2/6-31++G(2d,2p) level of theory with the potential energy surface computed using the AMBER molecular mechanics program and forcefields. While the Weiner et al and Cornell et al AMBER forcefields reasonably account for the imidazole N-H...O water interaction, these forcefields do not adequately account for the imidazole Ce-H...O water hydrogen bond. A forcefield modification is offered that results in excellent agreement between the ab initio and molecular mechanics geometry and energy for this C-H...O hydrogen bonded complex.

Catalytic role of the alpha-carboxylate of the Glu residue of glutathione in glutathione S-transferases.

The recently proposed catalytic role of the alpha-carboxylate of the Glu residue of glutathione in glutathione S-transferases (Widersten et al, Biochemistry 35, 7731-7742 (1996)) was examined. Based on structural considerations, it is clear that conformational changes in both glutathione and glutathione S-transferase are required. Recent kinetic studies by Ricci and coworkers (Ricci et al, J. Biol. Chem. 271, 16187-16192 (1996) and Caccuri et al, J. Biol. Chem. 271, 16193-16198 (1996)) may provide the missing evidence for these conformational changes. Possible ways to test this hypothesis are discussed.

DNA binding by TATA-box binding protein (TBP): a molecular dynamics computational study.

TATA-box binding protein (TBP) in a monomeric form and the complexes it forms with DNA have been elucidated with molecular dynamics simulations. Large TBP domain motions (bend and twist) are detected in the monomer as well as in the DNA complexes; these motions can be important for TBP binding of DNA. TBP interacts with guanine bases flanking the TATA element in the simulations of the complex; these interactions may explain the preference for guanine observed at these DNA positions. Side chains of some TBP residues at the binding interface display significant dynamic flexibility that results in 'flip-flop' contacts involving multiple base pairs of the DNA. We discuss the possible functional significance of these observations.

Molecular dynamics simulations indicate that F87W,T185F-cytochrome P450cam may reductively dehalogenate 1,1,1-trichloroethane.

Cytochrome P450cam is capable of reductively dehalogenating several chlorinated alkanes at low, but measurable, rates. In previous investigations of structure-function relationships in this enzyme using molecular dynamics simulations, we noticed that 1,1,1-trichloroethane (TCA) exhibits a very high degree of mobility in the active site due to its smaller molecular volume relative to the native substrate, camphor(1,2). Several amino acid sidechains lining the active site also exhibit significant dynamic fluctuations, possibly as a result of poor steric complementary to TCA. Guided by these results, we modeled double (F87W, T185F) and triple (F87W, T185F, V295I) mutants of P450cam, which provide additional bulk in the active site and increase the frequency of heme-substrate collision. Molecular dynamics simulations (300 ps on each protein) indicate that these mutants do not significantly perturb the three-dimensional fold of the enzyme, or local structure in the region of the active site. Both mutants bind the substrate more stably near the heme than the wild-type. Interestingly, however, the bulkier triple mutant seems to actually inhibit heme-substrate interactions relative to the double mutant. Over the final 200 ps of simulation, TCA is within 1 A of nonbonded contact with the heme 25% more often in the double mutant versus the wild-type. The triple mutant, on the other hand, binds TCA within 1 A of the heme only 15% as often as the wild-type. These results indicate that the double mutant may reductively dehalogenate TCA, a property not observed for the native protein. Implications for other experimentally measurable parameters are discussed.

Effect of warmup protocol and sampling time on convergence of molecular dynamics simulations of a DNA dodecamer using AMBER 4.1 and particle-mesh Ewald method.

This report describes one 3000 ps and two 1500 ps molecular dynamic simulations on a TATA box containing dodecamer DNA duplex in a periodic box of TIP3P water molecules, using the AMBER 4.1 implementation of the particle-mesh Ewald method. We compare the effect of warmup protocol and simulation time length on the root-mean square deviation (RMSD) parameter. For the longer simulation, the RMSD computed for the 500-1000 ps time interval is representative of longer time intervals, including 500-3000 ps. The various warmup protocols do not appear to have a significant effect on the simulation results. Based on the present results, DNA sequence-dependent differences in RMSD, or related properties, should exceed two standard deviations before being attributed to non-simulation factors, such as warmup protocol and sampling time effects; we recommend a minimum criterion of at least a three standard deviation difference with a sampling period of at least 500-1000 ps. In addition, while end effects appear negligible there is a consistent dependence of RMSD on DNA helix length.