Search

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 T's, the imino resonances of stems I and III and AH2's 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.

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.

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.

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.

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.

The review is presented of experimental and computational data on the influence of genotoxic modification of bases (deamination, alkylation, oxidation) on the structure and biological functioning of nucleic acids. Pathways are discussed for the influence of modification on coding properties of bases, on possible errors of nucleic acid biosynthesis, and on configurations of nucleotide mispairs. The atomic structure of nucleic acid fragments with modified bases and the role of base damages in mutagenesis and carcinogenesis are considered.

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.

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.

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.

The origin of torsional potentials in H3CSSCH3, H3CSSH, and HOOH and the anisotropy of the local charge distribution has been analyzed in terms of atomic multipoles calculated from the ab initio LCAO-MO-SCF wave function in the 6-31G* basis set. The results indicate that for longer -S-S-bonds the major contribution to these torsional barriers are electrostatic interactions of the atomic multipoles located on two atoms forming the rotated bond. This finding demonstrates the important role of electrostatic 1-2 interatomic interactions, usually neglected in conformational studies. It also opens the possibility to derive directly from accurate ab initio wave functions a simple nonempirical torsional potential involving atomic multipoles of two bonded atoms defining the torsional angle. For shorter -O-O- bonds, use of more precise models and inclusion of 1-3 interactions seems to be necessary.

A systematic analysis is presented of the algorithm for converting a virtual-bond chain, defined by the coordinates of the alpha-carbons of a given protein, into a complete polypeptide backbone. An alternative algorithm, based upon the same set of geometric parameters used in the Purisima-Scheraga algorithm but with a different "linkage map" of the algorithmic procedures, is proposed. The global virtual-bond chain geometric constraints are more easily separable from the loal peptide geometric and energetic constraints derived from, for example, the Ramachandran criterion, within the framework of this approach.

The phospholipid and Ca2+ dependent protein kinase (PKC) plays an essential role in a variety of cellular events. Inhibition of PKC was shown to arrest growth in tumor cell cultures making it a target for possible antitumor therapy. Calphostins are potent inhibitors of PKC with high affinity for the enzyme regulatory site. Structural characteristics of calphostins, which confer the inhibitory activity, are investigated by comparing their optimized structures with the existing models for PKC activation. The resulting model of inhibitory activity assumes interaction with two out of the three electrostatic interaction sites postulated for activators. The model shows two sites of hydrophobic interaction and enables the inhibitory activity of gossypol to be accounted for.

Two methods of side-chain placement on a modeled protein have been examined. Two molecular models of calmodulin were constructed that differ in the treatment of side chains prior to optimization of the molecule. A virtual bond analysis program developed by Purisima and Scheraga was used to determine the backbone conformation based on 2.2 angstroms resolution C alpha coordinates for the molecules. In the first model, side chains were initially constructed in an extended conformation. In the second model, a conformational grid search technique was employed. Calcium ions were treated explicitly during energy optimization using CHARMM. The models are compared to a recently published refined crystal structure of calmodulin. The results indicate that the initial choices for side-chains, but also significant effects on the main-chain conformation and supersecondary structure. The conformational differences are discussed. Analysis of these and other methods makes possible the formulation of a methodology for more appropriate side-chain placement in modeled proteins.

Recent discoveries have established the fact that RNA is capable of acting as an enzyme. In this study two different types of molecular orbital calculations, INDO and ab initio, were used in an attempt to assess the structural/functional role of the Mg2+ hydrated complex in ribozyme reactions. Preliminary studies indicate that the reaction is multistep and that the Mg2+ complex exerts a stabilizing effect on the intermediate or midpoint of the reaction.

Analysis of C-alpha atom positions from cysteines involved in disulphide bridges in protein crystals shows that their geometric characteristics are unique with respect to other Cys-Cys, non-bridging pairs. They may be used for predicting disulphide connections in incompletely determined protein structures, such as low resolution crystallography or theoretical folding experiments. The basic unit for analysis and prediction is the 3 x 3 distance matrix for Cx positions of residues (i - 1), Cys(i), (i +1) with (j - 1), Cys(j), (j + 1). In each of its columns, row and diagonal vector--outer distances are larger than the central distance. This analysis is compared with some analytical models.