Search PubMed⌕ Search

Biomedical subjects

R L Ornstein

Publications and source records attributed to R L Ornstein.

At least 55 records · Page 3Linked to original sources
PubMed1992

Cumulative atomic multipole moments complement any atomic charge model to obtain more accurate electrostatic properties.

The quality of several atomic charge models based on different definitions has been analyzed using cumulative atomic multipole moments (CAMM). This formalism can generate higher atomic moments starting from any atomic charges, while preserving the corresponding molecular moments. The atomic charge contribution to the higher molecular moments, as well as to the electrostatic potentials, has been examined for CO and HCN molecules at several different levels of theory. The results clearly show that the electrostatic potential obtained from CAMM expansion is convergent up to R-5 term for all atomic charge models used. This illustrates that higher atomic moments can be used to supplement any atomic charge model to obtain more accurate description of electrostatic properties.

Biopolymers↗
PubMed1992

A method for determining the positions of polar hydrogens added to a protein structure that maximizes protein hydrogen bonding.

An automated method for the optimal placement of polar hydrogens in a protein structure is described. This method treats the polar, side chain hydrogens of lysine, serine, threonine, and tyrosine and the amino terminus of a protein. The program, called NETWORK, divides the potential hydrogen-bonding pairs of a protein into groups of interacting donors and acceptors. A search is conducted on each of the local groups to find an arrangement which forms the most hydrogen bonds. If two or more arrangements have the same number of hydrogen bonds, the arrangement with the shortest set of hydrogen bonds is selected. The polar hydrogens of the histidyl side chain are specifically treated, and the ionization state of this residue is allowed to change, if this change results in additional hydrogen bonds for the local group. The program will accept Protein Data Bank as well as Biosym-format coordinate files. Input and output routines can be easily modified to accept other coordinate file formats. The predictions from this method are compared to known hydrogen positions for bovine pancreatic trypsin inhibitor, insulin, RNase-A, and trypsin for which the neutron diffraction structures have been determined. The usefulness of this program is further demonstrated by a comparison of molecular dynamics simulations for the enzyme cytochrome P-450cam with and without using NETWORK.

Animals↗
PubMed1992

Substrate mobility in a deeply buried active site: analysis of norcamphor bound to cytochrome P-450cam as determined by a 201-psec molecular dynamics simulation.

While cytochrome P-450cam catalyzes the hydroxylation of camphor to 5-exo-hydroxycamphor with 100% stereospecificity, norcamphor is hydroxylated by this enzyme yielding 45% 5-exo-, 47% 6-exo-, and 8% 3-exo-hydroxynorcamphor (Atkins, W.M., Sligar, S.G., J. Am. Chem. Soc. 109:3754-3760, 1987). The present study describes a 201-psec molecular dynamics (MD) stimulation of norcamphorbound cytochrome P-450cam to elucidate the relationship between substrate conformational mobility and formation of alternative products. First, these data suggest that the product specificity is, at least partially, due to the mobility of the substrate within the active site. Second, the high mobility of norcamphor in the active site leads to an average increase in separation between the heme iron and the substrate of about 1.0 A; this increase in separation may be the cause of the uncoupling of electron transfer when norcamphor is the substrate. Third, the active site water located in the norcamphorbound crystal structure possesses mobility that correlates well with the spin-state equilibrium of this enzyme-substrate complex.

Binding Sites↗
PubMed1992

Predicting the product specificity and coupling of cytochrome P450cam.

We present an analysis of several molecular dynamics trajectories of substrate-bound cytochrome P450cam. Trajectories were calculated for the native substrate, camphor, as well as for the alternative substrates, norcamphor and thiocamphor. The system modeled consisted of the crystallographically resolved amino acids, the heme group with a single oxygen atom as the distal ligand, the bound substrate, and the crystallographic waters. These trajectories of the presumptive ferryl oxygen intermediate were used to predict regiospecificity of hydroxylation and coupling between NADH consumption and product formation. Simple geometric criteria in combination with electronic considerations were used to calculate the probability of hydroxylation at specific sites on the substrate. We found that for all the cases examined, the predicted product ratios were in good agreement with the experimentally observed values. We also determined that these simple geometric criteria can be used to predict the degree of coupling between NADH consumption and product formation for a given substrate, which was in good agreement with the experimental values.

Camphor↗
PubMed1992

3 Nsec molecular dynamics simulation of the protein ubiquitin and comparison with X-ray crystal and solution NMR structures.

Mainly due to computational limitations, past protein molecular dynamics simulations have rarely been extended to 300 psec; we are not aware of any published results beyond 350 psec. The present work compares a 3000 psec simulation of the protein ubiquitin with the available x-ray crystallographic and solution NMR structures. Aside from experimental structure availability, ubiquitin was studied because of its relatively small size (76 amino acids) and lack of disulfide bridges. An implicit solvent model was used except for explicit treatment of waters of crystallization. We found that the simulated average structure retains most of the character of the starting x-ray crystal structure. In two highly surface accessible regions, the simulation was not in agreement with the x-ray structure. In addition, there are six backbone-backbone hydrogen bonds that are in conflict between the solution NMR and x-ray crystallographic structures; two are bonds that the NMR does not locate, and four are ones that the two methods disagree upon the donor. Concerning these six backbone-backbone hydrogen bonds, the present simulation agrees with the solution NMR structure in five out-of-the six cases, in that if a hydrogen bond is present in the x-ray structure and not in the NMR structure, the bond breaks within 700 psec. Of the two hydrogen bonds that are found in the NMR structure and not in the x-ray structure, one forms at 1400 psec and the other forms rarely. The present results suggest that relatively long molecular dynamics simulations, that use protein x-ray crystal coordinates for the starting structure and a computationally efficient solvent representation, may be used to gain an understanding of conformational and dynamic differences between the solid-crystal and dilute-solution states.

Computer Simulation↗
PubMed1992

A molecular dynamics simulation of bacteriophage T4 lysozyme.

An analysis of a 400 ps molecular dynamics simulation of the 164 amino acid enzyme T4 lysozyme is presented. The simulation was carried out with all hydrogen atoms modeled explicitly, the inclusion of all 152 crystallographic waters and at a temperature of 300 K. Temporal analysis of the trajectory versus energy, hydrogen bond stability, r.m.s. deviation from the starting crystal structure and radius of gyration, demonstrates that the simulation was both stable and representative of the average experimental structure. Average structural properties were calculated from the enzyme trajectory and compared with the crystal structure. The mean value of the C alpha displacements of the average simulated structure from the X-ray structure was 1.1 +/- 0.1 A; differences of the backbone phi and psi angles between the average simulated structure and the crystal structure were also examined. Thermal-B factors were calculated from the simulation for heavy and backbone atoms and both were in good agreement with experimental values. Relationships between protein secondary structure elements and internal motions were studied by examining the positional fluctuations of individual helix, sheet and turn structures. The structural integrity in the secondary structure units was preserved throughout the simulation; however, the A helix did show some unusually high atomic fluctuations. The largest backbone atom r.m.s. fluctuations were found in non-secondary structure regions; similar results were observed for r.m.s. fluctuations of non-secondary structure phi and psi angles. In general, the calculated values of r.m.s. fluctuations were quite small for the secondary structure elements. In contrast, surface loops and turns exhibited much larger values, being able to sample larger regions of conformational space. The C alpha difference distance matrix and super-positioning analyses comparing the X-ray structure with the average dynamics structure suggest that a 'hinge-bending' motion occurs between the N- and C-terminal domains.

Bacteriophage T4↗
PubMed1991

A 175-psec molecular dynamics simulation of camphor-bound cytochrome P-450cam.

The structure and internal motions of cytochrome P-450cam, a monooxygenase heme enzyme with 414 amino acid residues, with camphor bound at the active site have been evaluated on the basis of a 175-psec molecular dynamics simulation carried out at 300 K. All hydrogen atoms were explicitly modeled, and 204 crystallographic waters were included in the simulation. Based on an analysis of the time course of the trajectory versus potential energy, root mean square deviation, radius of gyration, and hydrogen bonding, the simulation was judged to be stable and representative of the average experimental structure. The averaged structural properties of the enzyme were evaluated from the final 135 psec of the simulation. The average atomic displacement from the X-ray structure was 1.39 A for all heavy atoms and 1.17 A for just C-alpha atoms. The average root-mean-square (rms) fluctuations of all heavy atoms and backbone atoms were 0.42 and 0.37 A, respectively. The computed rms fluctuations were in reasonable agreement with the experimentally determined temperature factors. All 13 segments of alpha-helix and 5 segments of beta-sheet were well preserved with the exception of the N-terminal half of helix F which alternated between an alpha-helix and a 310-helix. In addition there were in general only small variations in the relative orientation of adjacent alpha-helices. The rms fluctuations of the backbone dihedral angles in the secondary structure elements were almost uniformly smaller, with the fluctuation in alpha-helices and beta-sheets, 31 and 10% less, respectively, than those in nonsecondary structure regions. The reported crystal structure contains kinks in both helices C and I. In the simulation, both of these regions showed high mobility and large deviations from their starting positions. Since the kink in the I helix is at the oxygen binding site, these motions may have mechanistic implications.

Binding Sites↗
PubMed1991

Analysis of active site motions from a 175 picosecond molecular dynamics simulation of camphor-bound cytochrome P450cam.

The structure and internal motions of the active site residues of camphor-bound cytochrome P450cam have been evaluated on the basis of a 175 psec molecular dynamics simulation. The active site residues generally show very small deviations away from their starting crystal positions. These residues also generally show much smaller fluctuations than for the enzyme as a whole. Phe 87 is dynamically very unusual and is suggested to play a role in substrate movement into and/or out of the active site. The average distance between the heme iron and atoms C5, C6, and C3 of camphor is 5.3, 6.0, and 7.0 A, respectively. This trend is consistent with the experimentally observed stereospecificity of the hydroxylation reaction. On the basis of distance and angle criteria, both 5-exo and 5-endo hydrogen abstraction are predicted to occur during the hydroxylation reaction; although the 5-exo pathway is expected to be 3-fold more likely.

Binding Sites↗
PubMed1990

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.

Computer Simulation↗
PubMed1984

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.

Adenosine Triphosphate↗
PubMed1983

Correlation of crystallographically determined and computationally predicted hydrogen-bonded pairing configurations of nucleic acid bases.

Crystals of pairs of H-bonded nucleic acid bases are generally grown from nonaqueous solutions. We have been able to predict the H-bonded configuration of most of the base pairs in such crystals by using an empirical-potential function we recently developed for calculating the energetics of such interactions in chloroform solution. The following configurations were computationally predicted to predominate and are those observed in crystal structures: the Watson-Crick G.C configuration instead of two competing configurations; the Hoogsteen-type configurations for A.T, A.U, and A.br5U instead of Watson-Crick-type configurations; the Watson-Crick-type configurations for 2-aminopurine.br5U instead of the purine N3-type configuration; the Watson-Crick-type configurations for 8-bromo-2,6-diaminopurine.T instead of the Hoogsteen or purine-N3-type configurations; the syn-anti configuration for br8A.br8I instead of the anti-anti configuration; the Watson-Crick-type configurations for br8A.br5U instead of the Hoogsteen-type configurations; and the Hoogsteen-type configurations for me8A.T instead of the Watson-Crick configurations. In addition, the H-bonded base triplet br5U.2,6-diaminopurine.br5U was calculated to have Hoogsteen and Watson-Crick-type configurations but not the purine N3-type configuration. Apparently, lattice forces and chance nucleation of a minor base pairing configuration are not significant when the stability difference between the preferred and alternative configurations exceeds a relatively small value. In one case, in order to correctly predict the base pairing configuration in the crystal, it was necessary to include a contribution due to a C--H...O bond, suggesting that this type of H bond can make a significant contribution to base pair stability.

Base Composition↗
PubMed1980

Molecular models of induced DNA premutational damage and mutational pathways for the carcinogen 4-nitroquinoline 1-oxide and its metabolites.

The carcinogen 4-nitroquinoline 1-oxide (4NQO) and its metabolites undergo intercalative or covalent binding with DNA. Recent evidence indicates that the latter binding pattern is probably facilitated by an initial weaker intercalative interaction that can align potentially reactive sites on a 4NQO-metabolite and adjacent stacked bases. In the present study, we have proposed numerous possible covalent reaction products between 4NQO and its metabolites with DNA mini-helices based on chemical properties and key 'short-contacts' after energy-minimization in 21 different intercalative-like complexes. It is known from numerous experimental studies that 90% of the quinoline-bound DNAs in vivo involve guanine with the remaining 10% apparently involving adenine residues. The results of the present study suggest that this trend is not due to the greater affinity of the quinolines for guanine, but instead results from secondary processes involving the preferential formation of apurinic sites at aralkyl-adenine residues over that of aralkyl-guanine residues. In addition, observed mutational patterns can be rationalized in terms of the proposed reaction-products. The role of DNA repair mechanisms in the removal and correction of the different proposed reaction products are discussed. The binding pattern of several other aromatic carcinogens are similar to those depicted in the present work for the 4NQO-metabolites; hence the present study may be of some general significance.

4-Nitroquinoline-1-oxide↗