Robert Rein

Conformational Features of DNA Containing a Cis-Syn Photodimer

In an effort to understand the conformational and structural changes in DNA brought about by thymine photodimer, computer modeling and molecular mechanics energy calculations were performed on DNA hexamer and dodecamer duplexes containing a cis-syn photodimer. The conformation of the crystal structure of the cyanoethyl phosphate ester of the thymine dimer (Hruska et al., Biopolymers 25, 1399-1417 (1986)) was used in modeling the photodimer portion. Various starting conformations were used in the modeling procedure and the structures were minimized both retaining and later relaxing the crystallographic geometry of the cyclobutane ring. The results indicate that most of the deformation is restricted to the thymine dimer region, and that the conformational changes decrease rapidly on either side of the region containing the photodimer. The structural changes brought about by the introduction of the photodimer can be accommodated within six base paired duplex without significant bend in the DNA. More conformational changes are observed on the 5-side of the photodimer than on the 3-side. The conformational features, such as backbone torsion angles and sugar puckers, of the energy minimized structures are discussed in the context of the solution structures determined by NMR on a series of oligomers containing photodimers (Rycyna et al., Biochemistry 27, 3152-3163 (1988)).

Modeling DNA hydration: comparison of calculated and experimental hydration properties of nuclic acid bases.

Hydration properties of individual nucleic acid bases were calculated and compared with the available experimental data. Three sets of classical potential functions (PF) used in simulations of nucleic acid hydration were juxtaposed: (i) the PF developed by Poltev and Malenkov (PM), (ii) the PF of Weiner and Kollman (WK), which together with Jorgensens TIP3P water model are widely used in the AMBER program, and (iii) OPLS (optimized potentials for liquid simulations) developed by Jorgensen (J). The global minima of interaction energy of single water molecules with all the natural nucleic acid bases correspond to the formation of two water-base hydrogen bonds (water bridging of two hydrophilic atoms of the base). The energy values of these minima calculated via PM potentials are in somewhat better conformity with mass-spectrometric data than the values calculated via WK PF. OPLS gave much weaker water-base interactions for all compounds considered, thus these PF were not used in further computations. Monte Carlo simulations of the hydration of 9-methyladenine, 1-methyluracil and 1-methylthymine were performed in systems with 400 water molecules and periodic boundary conditions. Results of simulations with PM potentials give better agreement with experimental data on hydration energies than WK PF. Computations with PM PF of the hydration energy of keto and enol tautomers of 9-methylguanine can account for the shift in the tautomeric equilibrium of guanine in aqueous media to a dominance of the keto form in spite of nearly equal intrinsic stability of keto and enol tautomers. The results of guanine hydration computations are discussed in relation to mechanisms of base mispairing errors in nucleic acid biosynthesis. The data presented in this paper along with previous results on simulation of hydration shell structures in DNA duplex grooves provide ample evidence for the advantages of PM PF in studies of nucleic-acid hydration.

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.

Conformational requirements of suramin to target angiogenic growth factors

The conformational requirements of suramin for its binding to basic fibroblast growth factor (bFGF) and platelet-derived growth factor (PDGF) were examined by molecular modeling and docking simulations using the conformational features of suramin determined by the present proton nuclear magnetic resonance (1H-NMR) studies and the crystal structures of growth factors reported previously. The assignment of resonances of suramin to individual protons was accomplished by the combined analysis of the coupling constants, two-dimensional correlated spectroscopy (COSY) and nuclear Overhauser effect spectroscopy (NOESY). The NOESY data obtained for suramin were utilized in a conformational search algorithm with constraints to generate a family of conformers which were further refined by restrained energy minimization. A family of nine conformers generated by restrained modeling falls primarily into one of two categories, either the conformers two naphthyl rings are far apart, ∼28–30 Å, from one another or the conformers two naphthyl rings are relatively close, ∼16–18 Å. The NMR data provide evidence for the presence of more than one conformer in solution. The modeling and docking simulation studies suggest that suramin binds efficiently to bFGF and PDGF by an induced-fit mechanism, wherein suramin complements bFGF or PDGF by adjusting its conformational freedom around the two pairs of single bonds that link the middle phenyl rings to the secondary amide backbone. The interaction of suramin with bFGF or PDGF primarily involves ion-pair, hydrophobic and hydrogen bonding interactions, in addition to van der Waals contacts. The results indicate that suramin not only sterically blocks the receptor binding loop of the growth factors, but also competes for the binding sites of agonists such as heparin. The results suggest that suramins propensity to bind to several polypeptides of varying size and structure is due to its conformational flexibility. Collectively, the data emphasize that conformationally constrained suramin analogs that selectively and competitively target angiogenic growth factors could be designed to reduce non-specific binding and, accordingly, toxicity.

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

SummaryDistributed Point Charge Models (PCM) for CO, (H2O)2, and HS-SH molecules have been computed from analytical expressions using multicenter 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.

Molecular Modelling of Protein-Nucleic Acid Interactions

Computer modeling techniques to study the interaction of proteins with nucleic acids are presented. The methods utilize information from genetic and chemical modification experiments and macromolecular structural constraints. These techniques, in addition to computer model building procedures and theoretical energy calculations, are illustrated for the study of the lac and cro repressor-operator systems. Our predicted interactions between lac and its operator agree with those recently reported for lac based upon sequence alignment with the cro repressor. Several molecular models of the putative helical segment of cro interacting with its OR3 operator are presented. These models are reflective of intermediate conformations experienced by the repressor in recognition of the operator sequence. The results of our studies are further discussed in terms of the design of short peptides interacting with nucleic acid sequences and the evolutionary requirements in establishing these repressor interactions.

A Stereochemical Model of the Transpeptidation Complex

Molecular models are proposed to describe the relative arrangement of aminoacyl and peptidyl tRNAs when bound to their respective A and P sites on the ribosome. The crystallographically determined structures of tRNAasp and tRNAphe have served as the models for these bound structures, while the imposed steric constraints for the model complexes were based on the results of published experimental data. The constructed models satisfy the stereochemical requirements needed for codon-anticodon interaction and for peptide bond formation. In this paper, the results of the complex containing tRNAphe as the A and P site bound transfer RNAs, is compared to a similarly constructed model which uses tRNAasp as the ribosome-bound transfer RNAs. The models have the following three major features: 1) the aminoacyl and peptidyl transfer RNAs assume an angle of approximately 45 degrees relative to each other; 2) in providing the proper stereochemistry for peptide bond condensation, a significant kink must be present in the messenger RNA between the A site and P site codons; and 3) a comparison of the two model complexes indicates that structural variations between the tRNAs or any allosteric transitions of the transfer RNAs associated with codon-anticodon recognition may be accommodated in the model by way of freedom of rotation about the phosphate backbone bonds in the mRNA between consecutive codons.

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.

Structure of an Anti-HIV-1 Hammerhead Ribozyme Complex with a 17-mer DNA Substrate Analog of HIV-1 gag RNA and a Mechanism for the Cleavage Reaction: 750 MHz NMR and Computer Experiments

The structure of an anti-HIV-1 ribozyme-DNA abortive substrate complex was investigated by 750 MHz NMR and computer modeling experiments. The ribozyme was a chimeric molecule with 30 residues-18 DNA nucleotides, and 12 RNA residues in the conserved core. The DNA substrate analog had 17 residues. The chimeric ribozyme and the DNA substrate formed a shortened ribozyme-abortive substrate complex of 47 nucleotides with two DNA stems (stems I and III) and a loop consisting of the conserved core residues. Circular dichroism spectra showed that the DNA stems assume A-family conformation at the NMR concentration and a temperature of 15 degrees C, contrary to the conventional wisdom that DNA duplexes in aqueous solution populate entirely in the B-form. It is proposed that the A-family RNA residues at the core expand the A-family initiated at the core into the DNA stems because of the large free energy requirement for the formation of A/B junctions. Assignments of the base H8/H6 protons and H1 of the 47 residues were made by a NOESY walk. In addition to the methyl groups of all Ts, the imino resonances of stems I and III and AH2s were assigned from appropriate NOESY walks. The extracted NMR data along with available crystallographic data, were used to derive a structural model of the complex. Stems I and III of the final model displayed a remarkable similarity to the A form of DNA; in stem III, a GC base pair was found to be moving into the floor of the minor groove defined by flanking AT pairs; data suggest the formation of a buckled rhombic structure with the adjacent pair; in addition, the base pair at the interface of stem III and the loop region displayed deformed geometry. The loop with the catalytic core, and the immediate region of the stems displayed conformational multiplicity within the NMR time scale. A catalytic mechanism for ribozyme action based on the derived structure, and consistent with biochemical data in the literature, is proposed. The complex between the anti HIV-1 gag ribozyme and its abortive DNA substrate manifests in the detection of a continuous track of A.T base pairs; this suggests that the interaction between the ribozyme and its DNA substrate is stronger than the one observed in the case of the free ribozyme where the bases in stem I and stem III regions interact strongly with the ribozyme core region (Sarma, R. H., et al. FEBS Letters 375, 317-23, 1995). The complex formation provides certain guidelines in the design of suitable therapeutic ribozymes. If the residues in the ribozyme stem regions interact with the conserved core, it may either prevent or interfere with the formation of a catalytically active tertiary structure.

Comparing theoretical and experimental backbone-dependent sidechain conformational preferences for linear, branched, aromatic and polar residues

Abstract An Ecepp-3 conformational study based on a φ-ψ grid search with sidechain minimization was carried out on the N-acetyl N′-methyl amides of four representative amino acids: Met, Phe, Ile, and Ser, and the distribution of ξ 1 backbone-dependent rotamer preferences was compared with the similar distribution obtained from the backbone-dependent rotamer library for proteins developed by Dunbrack and Karplus (J. Mol. Biol. 230 (1993) 543). The experimental distribution is best reproduced theoretically in the case of the linear sidechain of Met, reasonably well for the bulky sidechains of the aromatic Phe and asymmetrically β-branched Ile, and only partially for the short polar sidechain of Ser. In the case of the Ser dipeptide the difference is accounted for by the missing H bonds.