Methylene blue binding to DNA with alternating GC base sequence: continuum treatment of salt effects.

Methylene blue (MB), an efficient singlet oxygen generating photoactive dye, binds to DNA and allows photosensitized reactions to be used for sequence-specific cleavage of the DNA backbone. Intercalation and groove binding are possible binding modes of the dye, depending on base sequences and environmental conditions. In a recent modeling study of methylene blue binding to a double stranded DNA decamer with an alternating GC sequence, six structural models for intercalation structures and for minor and major groove binding have been obtained. By estimating the binding energies (including electrostatic reaction field contributions of a salt-free aqueous solvent), symmetric intercalation at the 5'-CpG-3' and 5'-GpC-3' steps was found as the predominant binding mode, followed by a slightly weaker binding of the dye in the minor groove. In this study, the stability of the modeled structures has been analysed as a function of salt concentration. The results of finite difference numerical solutions of the non-linear Poisson-Boltzmann equation show that the stabilizing effect of salt is larger for free DNA than for the modeled MB-DNA complexes. Accordingly, the estimated binding energies decrease with increasing ionic strength. A slightly higher stabilization of the groove binding complexes results in comparable binding energies for symmetric intercalation and minor groove binding at high salt concentration. Both results are in qualitative agreement with experimental data.

Conformational deformability of RNA: a harmonic mode analysis.

The harmonic mode analysis method was used to characterize the conformational deformability of regular Watson-Crick paired, mismatch- and bulge-containing RNA. Good agreement between atomic Debye-Waller factors derived from x-ray crystallography of a regular RNA oligonucleotide and calculated atomic fluctuations was obtained. Calculated helical coordinate fluctuations showed a small sequence dependence of up to approximately 30-50%. A negative correlation between motions at a given base pair step and neighboring steps was found for most helical coordinates. Only very few calculated modes contribute significantly to global motions such as bending, twisting, and stretching of the RNA molecules. With respect to a local helical description of the RNA helix our calculations suggest that RNA bending is mostly due to a periodic change in the base pair step descriptors slide and roll. The presence of single guanine:uridine or guanine:adenine mismatches had little influence on the calculated RNA flexibility. In contrast, for tandem guanine:adenine base pairs the harmonic mode approach predicts a significantly reduced conformational flexibility in the case of a sheared arrangement and slightly enhanced flexibility for a face-to-face (imino proton) pairing relative to regular RNA. The presence of a single extra adenine bulge nucleotide stacked between flanking sequences resulted in an increased local atomic mobility around the bulge site (approximately 40%) and a slightly enhanced global bending flexibility. For an adenine bulge nucleotide in a looped-out conformation a strongly enhanced bulge nucleotide mobility but no increased bending flexibility compared to regular RNA was found.

Force field based conformational analysis of RNA structural motifs: GNRA tetraloops and their pyrimidine relatives.

The protocol of conformational analysis applied here to ribonucleotide oligomers combines conformational search in the space of torsion angles and energy minimization using the AMBER4.1 force field with a continuum treatment of electrostatic solute-solvent interactions. RNA fragments with 5'-GGGCGNNAGCCU-3' sequences commonly fold into hairpins with four-membered loops. The combinatorial search for acceptable conformations using the MC-SYM program was restricted to loop nucleotides and yielded roughly 1500 structures being compatible with a double-stranded stem. After energy minimization by the JUMNA program (without applying any experimental constraints), these structures converged into an ensemble of 74 different conformers including 26 structures which contained the sheared G-A base pair observed in experimental studies of GNRA tetraloops. Energetic analysis shows that inclusion of solvent electrostatic effects is critically important for the selection of conformers that agree with experimentally determined structures. The continuum model accounts for solvent polarization by means of the electrostatic reaction field. In the case of GNRA loop sequences, the contributions of the reaction field shift relative stabilities towards conformations showing most of the structural features derived from NMR studies. The agreement of computed conformations with the experimental structures of GAAA, GCAA, and GAGA tetraloops suggests that the continuum treatment of the solvent represents a definitive improvement over methods using simple damping models in electrostatic energy calculations. Application of the procedure described here to the evaluation of the relative stabilities of conformers resulting from searching the conformational space of RNA structural motifs provides some progress in (non-homology based) RNA 3D-structure prediction.

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model.

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A -form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions.

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.

Combining structural analysis of DNA with search routines for the detection of transcription regulatory elements.

MOTIVATION: Analysis of unannotated genomic sequences for regulatory regions depends on a reliable recognition of individual cis-acting elements. Since some of them have a very low conserved sequence pattern, additional criteria are required. RESULTS: Using molecular modelling techniques, we have created a complete database for the conversion of base sequences into profiles of structural parameters of DNA. On this basis, search routines can be developed that scan for profile matches. They may be used instead of or, probably most appropriate in most cases, in combination with conventional sequence pattern searches.

NMR studies and restrained-molecular-dynamics calculations of a long A+T-rich stretch in DNA. Effects of phosphate charge and solvent approximations.

The nonamer duplex d(GCAAAAACG).d(CGTTTTTGC) was studied by 1H-NMR at 500 MHz. With the exception of the H5' and H5" sugar protons, all protons were assigned by two-dimensional NMR experiments [NOE spectroscopy (NOESY), double-quantum-filtered J-correlated spectroscopy (DQF-COSY) and total correlation spectroscopy (TOCSY)]. The exchange kinetics of the imino protons of the Watson-Crick base pairing were studied at 15 degrees C by measuring inversion-recovery rates under conditions of extensive ammonia base catalysis. Extrapolation to infinite base concentration gave anomalous long lifetimes for the A-tract in accordance with previous results [Leroy, J.-L., Charettier, E., Kochoyan, M. & Guéron, M. (1988) Biochemistry 27, 8894-8898]. On average, 11 NOESY distance constraints/nucleotide were evaluated using the complete relaxation matrix approach. Deoxyribose coupling constants were obtained from simulations of the DQF-COSY cross-peaks, assuming a rapid two-state equilibrium between a C2'-endo and C3'-endo conformer. The sugars were found to be predominantly in the C2'-endo conformation. The NMR-derived distance and torsion constraints were implemented into three different restrained-molecular-dynamics (rMD) protocols, two in vacuo, with different charges on the phosphate group and the third with the solvent explicitly included. All protocols displayed good convergence from different starting structures. The structures derived from the three protocols satisfied experimental restraints equally well and had similar final energies. Although the overall pattern of sequence dependence of helical parameters shows some resemblance in all structures, we find that the absolute amplitudes of the parameters are largely dependent on the rMD protocols, particularly the twist parameters. The minor groove distance P(n + 2)-P(m + 2) varies from 0.7 nm to 1.2 nm in the three protocols. Still the NOESY-derived anomalously short distances AH2(n)-H1'(m + 1) and AH2(n)-H1'(n + 1), n and m denote complementary residues, which are assumed to be indicative of a compressed minor groove, are kept in all calculated structures.

Conformation of d(GGGATCCC)2 in crystals and in solution studied by X-ray diffraction, Raman spectroscopy and molecular modelling.

In the crystal, d(GGGATCCC)2 forms an A-DNA double helix as known from a single crystal X-ray diffraction study. Accordingly, in the Raman spectra of crystals the A-family marker bands at 664, 705, 807 and 1101 cm-1 and the spectral characteristics in the region 1200 to 1500 cm-1 clearly demonstrate the A-form as the dominant conformation. Bands at 691, 850, and 1080 cm-1, however, indicate that a minor fraction of the octamer molecules in the crystal is in an unusual, still not unequivocally identified conformation possibly belonging to the B-family. In solution, the octamer is in B-like conformation as shown by the presence of B-DNA Raman marker bands at 685, 837, 1094 and 1421 cm-1. Molecular modelling techniques lead to three structures with slightly different B-form geometries as the lowest energies models when a sigmoidal dielectric function with the bulk dielectric constant epsilon = 78 and the value q = -0.5e for the effective phosphate charges was used in the calculations. An A-form structure bearing a strong resemblance to the experimentally determined crystal structure becomes the lowest energy model structure when the electrostatic parameters are changed to epsilon = 30 and q = -0.25e, respectively.

Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI.

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.

Defining the structure of irregular nucleic acids: conventions and principles.

The algorithm "Curves", that we have recently presented in this journal (J. Biolmol. Str. Dynam. 6, 63-91 (1988], is updated to take into account the conventions developed at the Cambridge meeting on DNA curvature (September 1988) and extended to the calculation of local parameters. In addition, the principles which govern the choices made in establishing the Curves algorithm are compared with the approaches adopted by other authors.

Conformational and helicoidal analysis of 30 PS of molecular dynamics on the d(CGCGAATTCGCG) double helix: "curves", dials and windows.

A new procedure for the analysis of the structure and molecular dynamics of duplex DNA is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the values of the conformational and helicoidal parameters for each structure entering the analysis using the method "Curves" developed by Lavery and Sklenar, J. Biomol. Str. Dyn. 6, 63 (1988), followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels, or "dials", for each IUPAC torsional parameter and "windows" on the range of values assumed by the linear and angular helicoidal parameters, and is presented in a form isomorphous with the structure per se. The complete time evolution of the conformational and helicoidal parameters of a DNA double helix can then be depicted in a set of six composite figures. Dynamical aspects of helix bending are also subsumed in this analysis. The procedure is illustrated with an analysis of the structures of canonical A and B forms of DNA and the 300 degrees K native dodecamer duplex d(CGCGAATTCGCG). The "dials and windows" are then used for a comprehensive analysis of 30 psec of molecular dynamics on the dodecamer in the vicinity of a canonical B-DNA energy minimum. This involves presentation of the time evolution of 206 conformational and 230 helicoidal parameters for the dodecamer. A number of interesting structural features can be recognized in the analysis, including crankshaft motions, BI - BII transitions, sugar repuckerings, and a description of spontaneous helix bending at what corresponds to the 1 degrees and 2 degrees "hinge points" indicated in the crystal structure. Our approach is expected to be directly useful for critical analysis of the effects of various assumptions about force field parameters, hydration and electrostatic effects and thus contribute to the development of reliable simulation protocols for nucleic acid systems. Extension of the method to present differential changes in conformational and helicoidal parameters is expected to be valuable for the analysis of structural and molecular dynamics studies of the reorganization and adaptation of DNA on complexation with various drugs and regulatory proteins.

Describing protein structure: a general algorithm yielding complete helicoidal parameters and a unique overall axis.

We present a general and mathematically rigorous algorithm which allows the helicoidal structure of a protein to be calculated starting from the atomic coordinates of its peptide backbone. This algorithm yields a unique curved axis which quantifies the folding of the backbone and a full set of helicoidal parameters describing the location of each peptide unit. The parameters obtained form a complete and independent set and can therefore be used for analyzing, comparing, or reconstructing protein backbone geometry. This algorithm has been implemented in a computer program named P-Curve. Several examples of its possible applications are discussed.

The definition of generalized helicoidal parameters and of axis curvature for irregular nucleic acids.

An algorithm is presented which solves the problem of obtaining a rigorous helicoidal description of an irregular nucleic acid segment. Central to this approach is the definition of a function describing simultaneously the curvature of the nucleic acid segment in question and the corresponding stepwise variation of helicoidal parameters along the segment. Minimisation of this function leads to an optimal distribution of the conformational irregularity of the segment between these two components. Further, it is shown that this approach can be applied equally easily to single or double stranded nucleic acids. The results of this analysis yield both the absolute helicoidal parameters of individual bases/base pairs and the relative helicoidal parameters between successive bases/base pairs as well as the overall locus of the helical axis. The possibilities of this mathematical approach are demonstrated with the help of a computer program termed "Curves" which is applied to the study of a number of different nucleic acid structures.

The flexibility of the nucleic acids: (III). The interaction of an aliphatic diamine, putrescine, with flexible B-DNA.

A theoretical modelling of the interaction of putrescine (H3+N-(CH2)4-(+NH3) with DNA is carried out, introducing two new features which make the simulation of this interaction considerably more realistic. Firstly, the DNA to which putrescine is bound is fully flexible and thus able to respond to the distorting influence of the ligand. Secondly, the effect of changing the ratio of DNA base pairs per bound ligand is explicitly modelled. In this way, we have been able to confirm the experimentally known preference of putrescine binding with AT base pairs in B-DNA, but we also show, through the new features introduced, that the nature of the binding site of the ligand and the resulting impact on DNA conformation is strongly modified by the ligand binding density.

The flexibility of the nucleic acids: (I). "SIR", a novel approach to the variation of polymer geometry in constrained systems.

A novel and powerful methodology is developed which allows the alteration of molecular structures subjected to constraints and its application to polynucleotides with mononucleotide repeat symmetry, including the treatment of the flexible sugar rings is described. In contrast to procedures proposed by other authors, the constraints are formulated as differential equations which are linear with respect to the differentials of the geometrical variables. These equations can be solved easily by stepwise numerical integration involving sucessive infinitesimal rotations (SIR). Moreover, these equations define a set of independent curvilinear coordinates which can be used directly as the parameters of the energy functional in an energy minimisation procedure. This methodology allows the scanning of the full configurational space of a complex macromolecule, with direct access to the helicoidal variables in the case of periodic systems. Through this approach many problems involving biomacromolecular conformation, which would otherwise be intractable, may be studied with considerable ease.

The flexibility of the nucleic acids: (II). The calculation of internal energy and applications to mononucleotide repeat DNA.

Results concerning the flexibility of mononucleotide repeat DNA are presented using a novel methodology, denoted "SIR", to describe continuous changes in the structure of the nucleic acid. This methodology, combined with internal energy calculations and analytical energy gradients allows us to determine optimal conformations of poly(dG).poly(dC) and poly (dA).poly(dT) in both the A and B forms, taking into account the influence of the solvent medium and explicit counterions. Subsequently, several different types of distorsion of these optimal structures are investigated. It is shown that excellent correlation with experimental results concerning coupled changes in structural variables is obtained and several new correlations are also detected.

A structural model for the interaction of haem with unsaturated fatty acids explaining its quasi-lipoxygenase activity. Quantum chemical calculations.

The quasi-lipoxygenase activity of haemoglobin differs in many respects from the well-known haemin-catalyzed lipid peroxidation (1-4), especially in its high substrate specificity for unsaturated fatty acids containing one 1,4-pentadiene system (dienoic fatty acids). In this report a structural model for the fatty acid haem interaction based on quantum-chemical calculations is presented which show that only dienoic fatty acids are bound to the haem in such a way that the initial hydrogen abstraction that is necessary for the over-all reaction, is favoured sterically and energetically.