Search PubMed⌕ Search

Biomedical subjects

R L Ornstein

Publications and source records attributed to R L Ornstein.

At least 37 records · Page 2Linked to original sources
PubMed1995

Dramatic differences in the motions of the mouth of open and closed cytochrome P450BM-3 by molecular dynamics simulations.

Molecular dynamics trajectories were calculated separately for each of the two molecules in the asymmetric unit of the crystal structure of the hemoprotein domain of cytochrome P450BM-3. Each simulation was 200 ps in length and included a 10 A layer of explicit solvent. The simulated time-average structure of each P450BM-3 molecule is closer to its crystal structure than the two molecular dynamics time-averaged structures are to each other. In the crystal structure, molecule 2 has a more accessible substrate binding pocket than molecule 1, and this difference is maintained throughout the simulations presented here. In particular, the substrate docking regions of molecule 1 and molecule 2 diverge in the solution state simulations. The mouth of the substrate binding pocket is significantly more mobile in the simulation of molecule 2 than in the simulation of molecule 1. For molecule 1, the width of the mouth is only slightly larger than its X-ray value of 8.7 A and undergoes fluctuations of about 1 A. However, in molecule 2, the mouth of the substrate binding pocket is dramatically more open in the time-average molecular dynamics structure (14.7 A) than in the X-ray structure (10.9 A). Furthermore, this region of the protein undergoes large amplitude motions during the trajectory that are not seen in the trajectory of molecule 1, repeatedly opening and closing up to 7 A. Presumably, the binding of different substrates will induce the mouth region to adopt different conformations from within the wide range of structures that are accessible.

Binding Sites↗
PubMed1995

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.

Binding Sites↗
PubMed1995

Modeling of a possible conformational change associated with the catalytic mechanism in the hammerhead ribozyme.

Here we describe a possible model of the cleavage mechanism in the hammerhead ribozyme. In this model, the 2' hydroxyl of C17 is moved into an appropriate orientation for an in-line attack on the G1.1 phosphate through a change in its sugar pucker from C3' endo to C2' endo. This conformational change in the active site is caused by a change in the uridine turn placing the N2 and N3 atoms of G5 of the conserved core in hydrogen bonding geometry with the N3 and N2 atoms on the conserved G16.2 residue. The observed conformational change in the uridine turn suggests an explanation for the conservation of G5. In the crystal structure of H.M. Pley et al., Nature 372, 68-74 (1994), G5 is situated 5.3A away from G16.2. However, the uridine turn is sufficiently flexible to allow this conformational change with relatively modest changes in the backbone torsion angles (average change of 14.2 degrees). Two magnesium ions were modeled into the active site with positions analogous to those described in the functionally similar Klenow fragment 3'-5' exonuclease (L.S. Beese and T.A. Steitz, EMBO J. 10, 25-33 (1991)), the Group I intron (T.A. Steitz and J.A. Steitz, P.N.A.S. U.S.A. 90, 6498-6502 (1993); R.F. Setlik et al., J. Biomol. Str. Dyn. 10, 945-972 (1993)) and other phosphotransferases. Comparison of this model with one in which the uridine turn conformation was not changed showed that although the changes in the C17 sugar pucker could be modeled, insufficient space existed for the magnesium ions in the active site.

Base Sequence↗
PubMed1995

Enzyme-catalyzed dehalogenation of pentachloroethane: why F87W-cytochrome P450cam is faster than wild type.

Under anaerobic conditions, cytochromes P450 can reductively dehalogenate heavily halogenated hydrocarbons, such as one- and two-carbon organic solvents. This catalytic capacity has drawn attention to the potential use of engineered forms of P450s in the remediation of contaminated deep subsurface ecosystems. Loida (1994, PhD Thesis, University of Illinois at Urbana-Champaign, IL) and S.G. Sligar (personal communication) have observed recently that an active-site variant of cytochrome P450cam (F87W) dechlorinates pentachloroethane approximately three times faster than the wild-type enzyme. Molecular dynamics simulations have revealed that the mutant enzyme binding pocket remains smaller, and that pentachloroethane assumes configurations closer to the heme-Fe in the F87W mutant twice as often as in the wild-type enzyme. This result is consistent with a collisional model of dehalogenation, which agrees with experimental observations [Li and Wackett (1993) Biochemistry, 32, 9355-9361] that solutions containing wild-type P450cam dehalogenate pentachloroethane 100 times faster than those containing free heme. The simulations suggest that it is unlikely that Trp87 significantly stabilizes the developing negative charge on the substrate during carbon-halogen bond reduction. The design of improved microbial enzymes that incorporate both steric and electronic effects continues for use in remediating halogenated contaminants in situ.

Binding Sites↗
PubMed1995

Finding the global minimum: a fuzzy end elimination implementation.

The 'fuzzy end elimination theorem' (FEE) is a mathematically proven theorem that identifies rotameric states in proteins which are incompatible with the global minimum energy conformation. While implementing the FEE we noticed two different aspects that directly affected the final results at convergence. First, the identification of a single dead-ending rotameric state can trigger a 'domino effect' that initiates the identification of additional rotameric states which become dead-ending. A recursive check for dead-ending rotameric states is therefore necessary every time a dead-ending rotameric state is identified. It is shown that, if the recursive check is omitted, it is possible to miss the identification of some dead-ending rotameric states causing a premature termination of the elimination process. Second, we examined the effects of removing dead-ending rotameric states from further considerations at different moments of time. Two different methods of rotameric state removal were examined for an order dependence. In one case, each rotamer found to be incompatible with the global minimum energy conformation was removed immediately following its identification. In the other, dead-ending rotamers were marked for deletion but retained during the search, so that they influenced the evaluation of other rotameric states. When the search was completed, all marked rotamers were removed simultaneously. In addition, to expand further the usefulness of the FEE, a novel method is presented that allows for further reduction in the remaining set of conformations at the FEE convergence. In this method, called a tree-based search, each dead-ending pair of rotamers which does not lead to the direct removal of either rotameric state is used to reduce significantly the number of remaining conformations. In the future this method can also be expanded to triplet and quadruplet sets of rotameric states. We tested our implementation of the FEE by exhaustively searching ten protein segments and found that the FEE identified the global minimum every time. For each segment, the global minimum was exhaustively searched in two different environments: (i) the segments were extracted from the protein and exhaustively searched in the absence of the surrounding residues; (ii) the segments were exhaustively searched in the presence of the remaining residues fixed at crystal structure conformations. We also evaluated the performance of the method for accurately predicting side chain conformations. We examined the influence of factors such as type and accuracy of backbone template used, and the restrictions imposed by the choice of potential function, parameterization and rotamer database. Conclusions are drawn on these results and future prospects are given.

Algorithms↗
PubMed1994

An evaluation of implicit and explicit solvent model systems for the molecular dynamics simulation of bacteriophage T4 lysozyme.

In this report we examine several solvent models for use in molecular dynamics simulations of protein molecules with the Discover program from Biosym Technologies. Our goal was to find a solvent system which strikes a reasonable balance among theoretical rigor, computational efficiency, and experimental reality. We chose phage T4 lysozyme as our model protein and analyzed 14 simulations using different solvent models. We tested both implicit and explicit solvent models using either a linear distance-dependent dielectric or a constant dielectric. Use of a linear distance-dependent dielectric with implicit solvent significantly diminished atomic fluctuations in the protein and kept the protein close to the starting crystal structure. In systems using a constant dielectric and explicit solvent, atomic fluctuations were much greater and the protein was able to sample a larger portion of conformational space. A series of nonbonded cutoff distances (9.0, 11.5, 15.0, 20.0 A) using both abrupt and smooth truncation of the nonbonded cutoff distances were tested. The method of dual cutoffs was also tested. We found that a minimum nonbonded cutoff distance of 15.0 A was needed in order to properly couple solvent and solute. Distances shorter than 15.0 A resulted in a significant temperature gradient between the solvent and solute. In all trajectories using the proprietary Discover switching function, we found significant denaturation in the protein backbone; we were able to run successful trajectories only in those simulations that used no switching function. We were able to significantly reduce the computational burden by using dual cutoffs and still calculate a quality trajectory. In this method, we found that an outer cutoff distance of 15.0 A and an inner cutoff distance of 11.5 worked well. While a 10 A shell of explicit water yielded the best results, a 6 A shell of water yielded satisfactory results with nearly a 40% reduction in computational cost.

Bacteriophage T4↗
PubMed1994

An evaluation of molecular models of the cytochrome P450 Streptomyces griseolus enzymes P450SU1 and P450SU2.

P450SU1 and P450SU2 are herbicide-inducible bacterial cytochrome P450 enzymes from Streptomyces griseolus. They have two of the highest sequence indentities to camphor hydroxylase (P450cam from Pseudomonas putida), the cytochrome P450 with the first known crystal structure. We have built several models of these two proteins to investigate the variability in the structures that can occur from using different modeling protocols. We looked at variability due to alignment methods, backbone loop conformations and refinement methods. We have constructed two models for each protein using two alignment algorithms, and then an additional model using an identical alignment but different loop conformations for both buried and surface loops. The alignments used to build the models were created using the Needleman-Wunsch method, adapted for multiple sequences, and a manual method that utilized both a dot-matrix search matrix and the Needleman-Wunsch method. After constructing the initial models, several energy minimization methods were used to explore the variability in the final models caused by the choice of minimization techniques. Features of cytochrome P450cam and the cytochrome P450 superfamily, such as the ferredoxin binding site, the heme binding site and the substrate binding site were used to evaluate the validity of the models. Although the final structures were very similar between the models with different alignments, active-site residues were found to be dependent on the conformations of buried loops and early stages of energy minimization. We show which regions of the active site are the most dependent on the particular methods used, and which parts of the structures seem to be independent of the methods.

Amino Acid Sequence↗
PubMed1994

Active-site mobility inhibits reductive dehalogenation of 1,1,1-trichloroethane by cytochrome P450cam.

Recent studies by Wackett and co-workers have shown that cytochrome P450cam is capable of reductively dehalogenating hexachloroethane at a significant rate, but that no appreciable dehalogenation of 1,1,1-trichloroethane is observed. A growing body of evidence indicates that differences in intrinsic reactivity can not completely explain this observation. We therefore explored the possible role of differences in preferred binding orientation and in active-site mobility. A detailed analysis of molecular dynamics trajectories with each of these substrates bound at the active site of P450cam is presented. While the dynamics and overall time-average structure calculated for the protein are similar in the two trajectories, the two substrates behave quite differently. The smaller substrate, 1,1,1-trichloroethane, is significantly more mobile than hexachloroethane and has a preferred orientation in which the substituted carbon is generally far from the heme iron. In contrast, for hexachloroethane, one of the chlorine atoms is nearly always in van der Waals contact with the heme iron, which should favor the initial electron transfer step.

Binding Sites↗
PubMed1994

Investigation of domain motions in bacteriophage T4 lysozyme.

Hinge-bending in T4 lysozyme has been inferred from single amino acid mutant crystalline allomorphs by Matthews and coworkers. This raises an important question: are the different conformers in the unit cell artifacts of crystal packing forces, or do they represent different solution state structures? The objective of this theoretical study is to determine whether domain motions and hinge-bending could be simulated in T4 lysozyme using molecular dynamics. An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. Molecular dynamics calculations were computed using the Discover software package (Biosym Technologies). All hydrogen atoms were modeled explicitly with the inclusion of all 152 crystallographic waters at a temperature of 300 K. The native T4 lysozyme molecular dynamics simulation demonstrated hinge-bending in the protein. Relative domain motions between the N-terminal and C-terminal domains were evident. The enzyme hinge bending sites resulted from small changes in backbone atom conformations over several residues rather than rotation about a single bound. Two hinge foci were found in the simulation. One locus comprises residues 8-14 near the C-terminal of the A helix; the other site, residues 77-83 near the C-terminal of the C helix. Comparison of several snapshot structures from the dynamics trajectory clearly illustrates domain motions between the two lysozyme lobes. Time correlated atomic motions in the protein were analyzed using a dynamical cross-correlation map. We found a high degree of correlated atomic motions in each of the domains and, to a lesser extent, anticorrelated motions between the two domains. We also found that the hairpin loop in the N-terminal lobe (residues 19-24) acted as a mobile 'flap' and exhibited highly correlated dynamic motions across the cleft of the active site, especially with residue 142.

Bacteriophage T4↗
PubMed1994

A full-coordinate model of the polymerase domain of HIV-1 reverse transcriptase and its interaction with a nucleic acid substrate.

We present a full-coordinate model of residues 1-319 of the polymerase domain of HIV-I reverse transcriptase. This model was constructed from the x-ray crystallographic structure of Jacobo-Molina et al. (Jacobo-Molina et al., P.N.A.S. USA 90, 6320-6324 (1993)) which is currently available to the degree of C-coordinates. The backbone and side-chain atoms were constructed using the MAXSPROUT suite of programs (L. Holm and C. Sander, J. Mol. Biol. 218, 183-194 (1991)) and refined through molecular modeling. A seven base pair A-form dsDNA was positioned in the nucleic acid binding cleft to represent the template-primer complex. The orientation of the template-primer complex in the nucleic acid binding cleft was guided by the positions of phosphorus atoms in the crystal structure.

Amino Acid Sequence↗
PubMed1994

Application of 3-dimensional homology modeling of cytochrome P450 2B1 for interpretation of site-directed mutagenesis results.

Three-dimensional structures of cytochrome P450 2B1 were modeled based on the crystallographic structure of P450cam. The effect of the alignment, loop choice, and minimization with or without water was assessed. Although final models were similar in overall structure, the identity of active site residues depended upon the alignment. An example is Phe-206, which may or may not form part of the active site. The choice of the loop conformation had a lesser effect, while including water in the final minimization step was essential for preserving the shape and size of the active site. The best model (model 2) was in good agreement with the data from site-directed mutagenesis studies, and correctly predicted the effect of substitutions at 9 out of 10 amino acid positions. Thus, residues important for P450 2B1 activity, such as Ile-114, Phe-206, Ile-290, Thr-302, Val-363, and Gly-478, constitute part of the active site and are able to interact with the substrate androstenedione through hydrophobic interactions. On the other hand, Ser-303, Ser-360 and Lys-473 are far from the active site and/or cannot interact with the substrate, in agreement with experimental data. The model indicates other residues likely to be important for enzyme function, such as Tyr-111, Leu-209, Ile-477, and Ile-480, which can be tested experimentally. The substrate may assume numerous binding orientations consistent with observed patterns of hydroxylation at C15 and C16. The replacement in the model of certain amino acid residues to mimic residue substitutions from site-directed mutagenesis studies and docking of the substrate into the modified active site allowed a plausible explanation for alterations in regio- and stereospecificities of some mutants of P450 2B1, such as Gly-478-->Ala or Val-363-->Ala.

Amino Acid Sequence↗
PubMed1993

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides.

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.

Chemistry, Physical↗
PubMed1993

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.

Binding Sites↗
PubMed1993

Point charge representation of multicenter multipole moments in calculation of electrostatic properties.

Distributed Point Charge Models (PCM) for CO, (H2O)2, and HS-SH molecules have been computed from analytical expressions using multi-center multipole moments. The point charges (set of charges including both atomic and non-atomic positions) exactly reproduce both molecular and segmental multipole moments, thus constituting an accurate representation of the local anisotropy of electrostatic properties. In contrast to other known point charge models, PCM can be used to calculate not only intermolecular, but also intramolecular interactions. Comparison of these results with more accurate calculations demonstrated that PCM can correctly represent both weak and strong (intramolecular) interactions, thus indicating the merit of extending PCM to obtain improved potentials for molecular mechanics and molecular dynamics computational methods.

Carbon Monoxide↗
PubMed1993

Applicability of PM3 to transphosphorylation reaction path: toward designing a minimal ribozyme.

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on the system under investigation. In the present study we assess the validity of the MNDO/PM3 method on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.

Models, Biological↗
PubMed1993

Modeling study on the cleavage step of the self-splicing reaction in group I introns.

A three-dimensional model of the Tetrahymena thermophila group I intron is used to further explore the catalytic mechanism of the transphosphorylation reaction of the cleavage step. Based on the coordinates of the catalytic core model proposed by Michel and Westhof (Michel, F., Westhof, E. J. Mol. Biol. 216, 585-610 (1990)), we first converted their ligation step model into a model of the cleavage step by the substitution of several bases and the removal of helix P9. Next, an attempt to place a trigonal bipyramidal transition state model in the active site revealed that this modified model for the cleavage step could not accommodate the transition state due to insufficient space. A lowering of P1 helix relative to surrounding helices provided the additional space required. Simultaneously, it provided a better starting geometry to model the molecular contacts proposed by Pyle et al. (Pyle, A. M., Murphy, F. L., Cech, T. R. Nature 358, 123-128. (1992)), based on mutational studies involving the J8/7 segment. Two hydrated Mg2+ complexes were placed in the active site of the ribozyme model, using the crystal structure of the functionally similar Klenow fragment (Beese, L.S., Steitz, T.A. EMBO J. 10, 25-33 (1991)) as a guide. The presence of two metal ions in the active site of the intron differs from previous models, which incorporate one metal ion in the catalytic site to fulfill the postulated roles of Mg2+ in catalysis. The reaction profile is simulated based on a trigonal bipyramidal transition state, and the role of the hydrated Mg2+ complexes in catalysis is further explored using molecular orbital calculations.

Animals↗
PubMed1993

Substrate mobility in thiocamphor-bound cytochrome P450cam: an explanation of the conflict between the observed product profile and the X-ray structure.

Thiocamphor is an unusual substrate for P450cam in that in the X-ray structure it binds in the active site pocket in two distinct orientations and neither of these orientations are consistent with the 5-alcohol being the primary product. Other camphor analogs such as norcamphor or camphane bind in a single orientation consistent with the 5-alcohol being a major product. We present an analysis of four 175 ps molecular dynamics trajectories of thiocamphor-bound cytochrome P450cam. The first two trajectories were calculated for cytochrome P450cam with thiocamphor bound in both its major and minor crystallographic orientations. In the second set of simulations, a single oxygen atom was added as a distal ligand to the heme group in order to model the putative ferryl oxygen reaction intermediate. Trajectories were again calculated starting with thiocamphor in its major and minor orientations. While the protein dynamics were quite similar in all four trajectories, the substrate showed distinctly different motions in each of the trajectories. In particular, the preferred substrate orientations were very different in the presence of the ferryl oxygen than in the absence of that oxygen. The preferred orientations in the absence of the distal oxygen were consistent with the 3-alcohol being the major product, while the preferred orientations in the presence of the distal oxygen were consistent with the 5-alcohol being a major product. These simulations offer an explanation for the inconsistency between the X-ray data and the product profile.

Binding Sites↗
PubMed1992

Ethylbenzene hydroxylation by cytochrome P450cam.

The metabolism of ethylbenzene by cytochrome P450cam was analyzed by experimental and theoretical methods. The present experiments indicate that ethylbenzene is hydroxylated almost exclusively at the secondary ethyl carbon with about a 2:1 ratio of R:S product. Several molecular dynamics trajectories were performed with different starting conformations of ethylbenzene in the active site of P450cam. The stereochemistry of hydroxylation predicted from the molecular dynamics simulations was found to be in good agreement with the observed products.

Benzene Derivatives↗