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

Crystal structure of LacI member, PurR, bound to DNA: minor groove binding by alpha helices

The three-dimensional structure of a ternary complex of the purine repressor, PurR, bound to both its corepressor, hypoxanthine, and the 16-base pair purF operator site has been solved at 2.7 A resolution by x-ray crystallography. The bipartite structure of PurR consists of an amino-terminal DNA-binding domain and a larger carboxyl-terminal corepressor binding and dimerization domain that is similar to that of the bacterial periplasmic binding proteins. The DNA-binding domain contains a helix-turn-helix motif that makes base-specific contacts in the major groove of the DNA. Base contacts are also made by residues of symmetry-related alpha helices, the "hinge" helices, which bind deeply in the minor groove. Critical to hinge helix-minor groove binding is the intercalation of the side chains of Leu54 and its symmetry-related mate, Leu54', into the central CpG-base pair step. These residues thereby act as "leucine levers" to pry open the minor groove and kink the purF operator by 45 degrees.

Stacking-unstacking of the dinucleoside monophosphate guanylyl-3',5'-uridine studied with molecular dynamics

Molecular dynamics simulations were carried out on two conformations of the dinucleoside monophosphate guanylyl-3',5'-uridine (GpU) in aqueous solution with one sodium counterion. One stacked conformation and one with the C3'-O3'-P-O5' backbone torsion angle twisted 180 degrees to create an unstacked conformation. We observed a relatively stable behavior of the stacked conformation, which remained stacked throughout the simulation, whereas the unstacked conformation showed major changes in the backbone torsion and glycosidic angles. During the simulation the unstacked conformation transformed into a more stacked form and then back again to an unstacked one. The calculated correlation times for rotational diffusion from the molecular dynamics simulations are in agreement with fluorescence anisotropy and nuclear magnetic resonance data. As expected, the correlation times for rotational diffusion of the unstacked conformation were observed to be longer than for the stacked conformation. The 2'OH group may contribute in stabilizing the stacked conformation, where the O2'-H...O4' hydrogen bond occurred in 82.7% of the simulation.

Structure-function studies of DNA damage using ab initio quantum mechanics and molecular dynamics simulation

Studies of ring-saturated pyrimidine base lesions are used to illustrate an integrated modeling approach that combines quantum-chemical calculations with molecular dynamics simulation. Electronic structure calculations on the lesions in isolation reveal strong conformational preferences due to interactions between equatorial substituents to the pyrimidine ring. Large distortions of DNA should result when these interactions force the methyl group of thymine to assume an axial orientation, as is the case for thymine glycol but not for dihydrothymine. Molecular dynamics simulations of the dodecamer d(CGCGAATTCGCG)2 with and without a ring-saturated thymine lesion at position T7 support this conclusion. Implications of these studies for recognition of thymine lesions by endonuclease III are also discussed.

Molecular dynamics simulations of oligonucleotides in solution: visualization of intrinsic curvature

We have undertaken molecular dynamics simulations on the d(CGCAAAAAAGCG).d(CGCTTTTTTGCG) dodecamer in solution. In this study, we focus on aspects of conformation and dynamics, including the possibility of cross-strand hydrogen bonds. We compare our results with those from crystallography as well as infrared, Raman and NMR spectroscopy and cyclization kinetics. Our method of analysis allows us to visualise the curvature of the helix as a function of time during the simulation. We find that the major distortions of the helix axis path occur at the junctions between the (essentially straight) A-tract and the CG- and GC-tracts, although at one junction this is due to hyperflexibility (i.e., regions of high flexibility with no preferred direction of curvature), while at the other junction a static curvature is found (i.e., a preferred, sustained direction of curvature).

Measuring the geometry of DNA grooves

We present a new method for measuring the widths and depths of the grooves formed within DNA helices. This method overcomes the limitations of simply measuring interstrand phosphate-phosphate distances and has the advantage of yielding continuous values for groove geometry along a DNA fragment. In the case of oligonucleotides, it also clearly indicates the zones in which grooves exist, bounded by two phosphodiester backbones. The methodology has been developed within the Curves algorithm for studying irregular DNA geometries and is based on the optimal, and generally curved, helical axis obtained by this analysis.

DNA recognition by beta-sheets in the Arc repressor-operator crystal structure

Transcription of the ant gene during lytic growth of bacteriophage P22 (ref. 1) is regulated by the cooperative binding of two Arc repressor dimers to a 21-base-pair operator site. Here we report the co-crystal structure of this Arc tetramer-operator complex at 2.6 A resolution. As expected from genetic and structural studies and from the co-crystal structure of the homologous Escherichia coli MetJ repressor, each Arc dimer uses an antiparallel beta-sheet to recognize bases in the major groove. However, the Arc and MetJ complexes differ in several important ways: the beta-sheet-DNA interactions of Arc are far less symmetrical; DNA binding by Arc is accompanied by important conformational changes in the beta-sheet; and Arc uses a different part of its protein surface for dimer-dimer interactions.

Solution structure of the octamer motif in immunoglobulin genes via restrained molecular dynamics calculations

The solution structure of the DNA decamer d(CATTTGCATC)-d(GATGCAAATG), comprising the octamer motif of immunoglobulin genes, is determined by restrained molecular dynamics (rMD) simulations. The restraint data set includes interproton distances and torsion angles for the deoxyribose sugar ring which were previously obtained by a complete relaxation matrix analysis of the two-dimensional nuclear Overhauser enhancement (2D NOE) intensities and by the quantitative simulation of cross-peaks in double-quantum-filtered correlated (2QF-COSY) spectra. The influence of torsion angles and the number of experimental distance restraints on the structural refinement has been systematically examined. Omitting part of the experimental NOE-derived distances results in reduced restraint violations and lower R factors but impairs structural convergence in the rMD refinement. Eight separate restrained molecular dynamics simulations were carried out for 20 ps each, starting from either energy-minimized A- or B-DNA. Mutual atomic root-mean-square (rms) differences among the refined structures are well below 1 A and comparable to the rms fluctuations of the atoms about their average position, indicating convergence to essentially identical structures. The average refined structure was subjected to an additional 100 ps of rMD simulations and analyzed in terms of average torsion angles and helical parameters. The B-type duplex exhibits clear sequence-dependent variations in its geometry with a narrow minor groove at the T3.A3 tract and a large positive roll at the subsequent TG.CA step. This is accompanied by a noticeable bend of the global helix axis into the major groove. There is also evidence of significant flexibility of the sugar-phosphate backbone with rapid interconversion among different conformers.

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates

beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.

Molecular dynamics simulations of poly(dA).poly(dT): comparisons between implicit and explicit solvent representations

The program AMBER 3.0 has been used to generate molecular dynamics trajectories of a poly(dA).poly(dT) decamer. The simulations were performed using different methods to treat solvent effects. Results of a simulation including 18 counterions NH4+ and 4109 water molecules under (N, P, T) conditions were compared to simulation runs with implicit solvent representation in which solvent screening effects were represented by the use of a sigmoidal distance-dependent dielectric function. In the latter case, the system was simulated under microcanonical (N, V, E) and canonical (N, V, T) conditions. For the fully hydrated system simulation, a preequilibration protocol was developed since it was observed that long and progressive periods of heating and equilibration on the overall system were necessary in order to avoid energetic collisions between the solute and the solvent molecules, leading to severe irreversible deformation of the solute. A detailed analysis of DNA conformations, sugar puckers, and stability of the hydrogen bonds, Watson-Crick and three-center H bonds, is reported. The results show that DNA remains essentially in the B conformer with a tendency in the hydrated model to adopt a slightly distorted, unwound, and stretched conformation in comparison to standard B-DNA. Concerning sugar puckers, the mean pseudorotation phases of the adenine residues are systematically higher than those of the thymine residues, except in the case of the hydrated model for which a articular behavior is observed for the adenine strand. In this case, the terminal bases oscillate between C2'-endo and O4'-endo and the central ones stay in the C3'-endo domain. The mean lifetimes of the internal Watson-Crick H-bond (A) HN6...O4(T) are also dependent on the base pairs included in the calculation, excepted for the implicit solvent simulation at constant temperature. The three-center H bonds have very small mean lifetimes in all three cases of MD simulation. In the minor groove of the hydrated model, a spine of hydration is found as observed by x-ray crystallography and other theoretical simulations. On the basis of the rms deviations, it appears that the fully hydrated simulation has not reached a plateau at the end of the run, while the implicit simulation at constant energy seems to have converged. At constant temperature, very large oscillations in rms deviations are observed.

Molecular dynamics investigations of DNA triple helical models: unique features of the Watson-Crick duplex

We have built computer models of triple helical structures with a third poly(dT) strand Hoogsteen base paired to the major groove of a poly(dA).poly(dT) Watson-Crick (WC) base-paired duplex in the canonical A-DNA as well as B-DNA. For the A-DNA form, the sugar-phosphate backbone of the third strand intertwines and clashes with the poly(dA) strand requiring a radical alteration of the duplex to access the hydrogen bonding sites in the major groove. In contrast, when the duplex was in the canonical B-DNA form, the third strand was readily accommodated in the major groove without perturbing the duplex. The triple helical model, with the duplex in the B-DNA form, was equilibrated for 400ps using molecular dynamics simulations including water molecules and counter-ions. During the entire simulations, the deoxyriboses of the adenine strand oscillate between the S-type and E-type conformations. However, 30% of the sugars of the thymine strands-II & III switch to the N-type conformation early in the simulations but return to the S-type conformation after 200ps. In the equilibrium structure, the WC duplex portion of the triplex is unique and its geometry differs from both the A- or B-DNA. the deoxyriboses of the three strands predominantly exhibit S-type conformation. Besides the sugar pucker, the major groove width and the base-tilt are analogous to B-DNA, while the X-displacement and helical twist resemble A-DNA, giving a unique structure to the triplex and the Watson & Crick and Hoogsteen duplexes.

SUBCUR: visualization of structural differences between DNA duplexes

A computer program, SUBCUR, is described which permits analysis and rapid identification of geometrical differences and patterns of variance between two DNA duplexes. The program is compatible with the CURVES 3.1 package and allows graphical visualization of the structural differences. Examples are provided which illustrate the applicability of the program in analyzing the different backbone conformations of two helices and the different curvatures of two helices.

Molecular dynamics study of the base pair opening process in the self-complementary octanucleotide d(CTGATCAG)

We report an analysis of a 200 ps Molecular Dynamics simulation of the double stranded oligonucleotide d(CTGATCAG) in the presence of 1534 water molecules and 14 Na+ ions. We focus on the opening process of Thymine 5, by analyzing in detail the glycosidic bond rotational motion about the helix axis. The present analysis is mainly based on autocorrelation functions and on mean square displacements. We show that the opening of the base has a Brownian character and we find a rotational diffusion coefficient of 4.7 rad2s-1. Furthermore we estimate the DNA torsional constant to be about 0.5 10(-18) J.rad-2 and the RMS of the angular displacement to be 8.3 degrees. All these values are in fair agreement with those determined experimentally by fluorescence polarization of DNA-Ethidium bromide complexes. This shows that the rotational motions of the bases detected in the range 10(-9)-10(-7) s. by fluorescence techniques are the same as those analyzed in the present study (10(-12)-2 10(-10) s).

Solution structure of the exocyclic 1,N2-propanodeoxyguanosine adduct opposite deoxyadenosine in a DNA nonamer duplex at pH 8.9. model of pH-dependent conformational transition

The solution structure of the complementary d(C1-A2-T3-G4-X5-G6-T7-A8-C9).d(G10-T11-A12-C13-A14-C15-A16-T17-G18) DNA duplex (designated X.A 9-mer), which contains a 1,N2-propanodeoxyguanosine exocyclic adduct X5 opposite deoxyadenosine A14 at the center, is pH dependent [Kouchakdjian, M., Eisenberg, M., Live, D., Marinelli, E., Grollman, A., & Patel, D.J. (1990) Biochemistry 29, 4456-4465]. In our previous paper [Huang, P., & Eisenberg, M. (1992) Biochemistry 31, 6518-6532] we established the three-dimensional structure of this X.A 9-mer duplex at pH 5.8 by use of restrained molecular dynamics followed by NOE-based back-calculation refinement. The present paper discusses the structure at pH 8.9 and the pH-dependent conformational transition between the structures at pH 5.8 and at pH 8.9. The structure at pH 8.9 is calculated starting from five different conformations. The final structures converge and agree well with the experimental NOE intensities. These structures are essentially B-type DNA (with X5 and A14 in the BII conformation while the other residues are in the most commonly described BI conformation) and display an approximate 27 degrees kink at the center of the helix. At the kink site, X5 is positioned in the major groove with the exocyclic ring directed toward the G6.C13 base pair, unstacked from the flanking base G6 and exposed to the solvent. A14, opposite the lesion, remains stacked with its neighbor C15, but not with C13. The kinked helix can accommodate the rotation of the bulky X5 about its glycosidic bond. We propose here a model for the pH-dependent transition. Our model explains the conformational change, which includes the anti and syn rotation of the bulky adduct around its glycosidic bond, with a minimal energy barrier and with an overall kink of the DNA helix. These new findings, fully consistent with the NMR experimental data, were revealed only after restrained dynamics refinement. Distance-restrained energy minimization by itself was insufficient, as shown by the previous NMR study.

Molecular dynamics simulations of B-DNA: an analysis of the role of initial molecular configuration, randomly assigned velocity distribution, long integration times, and nonconstrained termini

Molecular dynamics simulations of three DNA sequences using the AMBER 3.0 force field were performed with implicit inclusion of water through a distance-dependent dielectric constant and solvated counterions. Simulations of the self-complementary DNA dodecamer d(CGCGAATTCGCG) were started from a regular B-DNA structure and the x-ray single crystal B-DNA structure. Although mean convergence during the 89-ps calculation was confirmed, localized differences in backbone torsionals and base-pair helicoidals were observed. A nanosecond simulation of the nonself-complementary 14 base-pair DNA d(GGCGGAATTGGCGG) indicates that most structural parameters stabilize within the first 100-200 ps, while isolated features show low-frequency oscillations throughout the calculation. The lack of harmonic constraints on the ends of the molecules was shown not to perturb the structural dynamics of the internal oligonucleotide beyond the external 2 base pairs. Comparison of three simulations of the nonself-complementary 14 base-pair DNA d(GGCGAAATTCGCGG), identical in all respects other than the assignment of initial Maxwellian atomic velocity distributions, revealed the inherent systematic variability. The three calculations result in nearly superimposable global structures, with localized variations in torsionals and helicoidals. Our results provide a basis for performing a comparative analysis of the effect of DNA sequence on localized structure.

Persistence analysis of the static and dynamical helix deformations of DNA oligonucleotides: application to the crystal structure and molecular dynamics simulation of d(CGCGAATTCGCG)2

A theory and graphical presentation for the analysis of helix structure and deformations in oligonucleotides is presented. The parameters "persistence" and "flexibility" as defined in the configurational statistics of polymers of infinite length are reformulated at the oligonucleotide level in an extension of J. A. Schellman's method [(1974) Biopolymers, Vol. 17, pp. 217-226], and used as a basis for a systematic "Persistence Analysis" of the helix deformation properties for all possible subsequences in the structure. The basis for the analysis is a set of link vectors referenced to individual base pairs, and is limited to sequences exhibiting only perturbed rod-like behavior, i.e., below the threshold for supercoiling. The present application of the method is concerned with a physical model for the angular component of bending, so the link vectors are defined as the unit components of a global helix axis obtained by the procedure "Curves" of R. Lavery and H. Sklenar [(1988) J. Biomol. Struct. Dynam., Vol. 6, pp. 63-91; (1989) ibid., Vol. 6, pp. 655-667]. A discussion of the relationship between global bending and relative orientation of base pairs is provided. Our approach is illustrated by analysis of some model oligonucleotide structures with intrinsic kinks, the crystal structure of the dodecamer d(CGCGAATTCGCG)2, and the results of two molecular dynamics simulations on this dodecamer using two variations of the GROMOS force field. The results indicate that essentially all aspects of curvature in short oligonucleotides can be determined, such as the position and orientation of each bend, the sharpness or smoothness, and the location and linearity of subsequences. In the case of molecular dynamics simulations, where a Boltzmann ensemble of structures is analyzed, the spatial extent of the deformations (flexibility) is also considered.

A molecular dynamics study of conformational changes and hydration of left-handed d(CGCGCGCGCGCG)2 in a nonsalt solution

Twelve dinucleotides (one complete turn) of left-handed, flexible, double-helix poly(dG-dC) Z-DNA have been simulated in aqueous solution with K+ counterions for 70 ps. Most of the d(GpC) phosphates have rotated in accordance with a ZI----ZII transition. The ZII conformation was probably partly stabilized by counterions, which coordinate one of the anionic oxygens and the guanine-N7 of the next (5'----3' direction) base. The presence of base-coordinating ions close to the helical axis rotated and pulled about half of the d(CpG) phosphates further into the groove. These ions also gave rise to rather large deviations from the crystal structure (ZI) with their tendency of pulling the bases closer toward the helical axis. A flipping of the orientation about the glycosyl bond from the +sc to the -sc region was observed for one guanosine, also leading to deviations from the crystal structure. Many bridges containing one or two water molecules were found, with a dominance for the latter. They essentially formed a network of intra- and interstrand bridges between anionic and esterified phosphate oxygens. A "spine" of water molecules could be distinguished as a dark zig-zag pattern in the water density map. The lifetime of a bridge containing one water was about twice as long as that of a two-water bridge and it lasted 5-15 times longer than a hydrogen bond in water. The lifetimes were also calculated for a selection of bridge types, in order of decreasing stability: O1P/O2P ... W ... O'4 much greater than O1P/O2P ... W ... guanine-N2 greater than O1P/O2P ... W ... O1P/O2P. The reorientational motion of water molecules in the first hydration shell around selected groups was slowed down considerably compared to bulk water and the decreasing order of correlation times was guanine-N2 greater than O'4 greater than O'3/O'5 greater than O1P/O2P.

A molecular dynamics simulation of a polyamine-induced conformational change of DNA. A possible mechanism for the B to Z transition

A 75ps molecular dynamics simulation has been performed on a fully solvated complex of spermine with the B DNA decamer (dGdC)5.(dGdC)5. The simulation indicates a possible mechanism by which polyamines might induce the formation of a left-handed helix, the B to Z transition. Spermine was initially located in the major groove, hydrogen bonded to the helix. During the simulation the ligand migrates deeper into the DNA, maintaining strong hydrogen bonding to the central guanine bases and destroying the Watson-Crick base pairing with their respective cytosines. Significant rotation of these and other cytosine bases was observed, in part due to interactions of the helix with the aminopropyl chains of spermine. An intermediate BII conformation might be of importance in this process.

Conformation and dynamics of drug-DNA intercalation

Molecular dynamics simulations have been undertaken for a B-form dodecanucleotide duplex in solution with and without an intercalated proflavine molecule between the central C.G base pairs. The introduction of this simple intercalator affects both the conformational features and dynamic properties of the oligonucleotide double helix. Changes are seen in the rms atomic fluctuations and anisotropy of phosphate, sugar and base atoms. The backbone conformation is slightly changed on average and more sugars adopt the C3' endo conformation in the simulation of the complex compared with the simulation of the oligonucleotide alone. Both major and minor grooves becomes wider on average with the addition of the intercalating drug. Flanking A.T base pairs on both sides of the intercalation site have undergone an increase in flexibility, with the base pairs, especially at the 5' side, having the N1...N3 hydrogen bonds being broken.

Molecular dynamics of HIV-1 protease

Molecular dynamics simulations have been carried out based on the GROMOS force field on the aspartyl protease (PR) of the human immunodeficiency virus HIV-1. The principal simulation treats the HIV-1 PR dimer and 6990 water molecules in a hexagonal prism cell under periodic boundary conditions and was carried out for a trajectory of 100 psec. Corresponding in vacuo simulations, i.e., treating the isolated protein without solvent, were carried out to study the influence of solvent on the simulation. The results indicate that including waters explicitly in the simulation results in a model considerably closer to the crystal structure than when solvent is neglected. Detailed conformational and helicoidal analysis was performed on the solvated form to determine the exact nature of the dynamical model and the exact points of agreement and disagreement with the crystal structure. The calculated dynamical model was further elucidated by means of studies of the time evolution of the cross-correlation coefficients for atomic displacements of the atoms comprising the protein backbone. The cross-correlation analysis revealed significant aspects of structure originating uniquely in the dynamical motions of the molecule. In particular, an unanticipated through-space, domain-domain correlation was found between the mobile flap region covering the active site and a remote regions of the structure, which collectively act somewhat like a molecular cantilever. The significance of these results is discussed with respect to the inactivation of the protease by site-specific mutagenesis, and in the design of inhibitors.

Toward a dynamical structure of DNA: comparison of theoretical and experimental NOE intensities

Comparisons of experimental and calculated interproton nuclear Overhauser effect (NOE) buildup curves for duplex d(CGCGAATTCGCG)2 have been made. The calculated NOEs are based on molecular dynamics simulations including counterions and water and on the single-structure canonical A, B, and crystal forms. The calculated NOE effects include consideration of the motions of individual interproton vectors and the anisotropic tumbling of the DNA. The effects due to inclusion of anisotropic tumbling are much larger than those due to the local motion, and both improve the agreement between calculated and experimental results. The predictions based on the dynamical models agree significantly better with experiment than those based on either of the canonical forms or the crystal structure.

Molecular dynamics simulations of dinucleoside and dinucleoside-drug crystal hydrates

Molecular dynamics simulations have been performed on the dinucleoside monophosphates rGpC and dCpG, the latter in its intercalation complex with the acridine drug proflavine. The simulations were performed on the crystal structures, with crystallographically-located solvent molecules. It was found that satisfactory results were best obtained with restraints placed on the movements of the water molecules. Motions of individual atoms have been examined in terms of rms fluctuations and anisotropy and correlation functions. Relative motions of groups (phosphates, sugars, bases and proflavine molecules) have been analysed.

Crystal structure of a CAP-DNA complex: the DNA is bent by 90 degrees

The 3 angstrom resolution crystal structure of the Escherichia coli catabolite gene activator protein (CAP) complexed with a 30-base pair DNA sequence shows that the DNA is bent by 90 degrees. This bend results almost entirely from two 40 degrees kinks that occur between TG/CA base pairs at positions 5 and 6 on each side of the dyad axis of the complex. DNA sequence discrimination by CAP derives both from sequence-dependent distortion of the DNA helix and from direct hydrogen-bonding interactions between three protein side chains and the exposed edges of three base pairs in the major groove of the DNA. The structure of this transcription factor--DNA complex provides insights into possible mechanisms of transcription activation.

Are time-averaged restraints necessary for nuclear magnetic resonance refinement? A model study for DNA

A recently suggested method for refinement of structural data obtained from two-dimensional nuclear magnetic resonance experiments using molecular dynamics (MD) is explored. In this method, the time-averaged values of the appropriate internal co-ordinates of the molecule, calculated from the MD trajectory, are driven by restraints towards the experimental target values. This contrasts with most refinement procedures currently in use, where restraints are applied based on the instantaneous values of the appropriate co-ordinates. Both refinement methods are applied to the EcoRI restriction site DNA hexamer d(GAATTC)2, using target nuclear Overhauser enhancement distances derived from a one nanosecond unrestrained MD simulation of this structure. The resulting refined structures are compared to the results of the unrestrained MD trajectory, which serves as our "experimental" data. We show that although both methods can yield an average structure with the correct gross morphology, the new method allows both a much more realistic picture of inherent flexibility, and reproduces fine conformational detail better, such as sequence dependency. We also analyze the very long MD trajectory generated here (longer than any previously reported for a DNA oligomer), and find that significantly shorter simulations, typical of those frequently performed, may not yield acceptably reliable values for certain structural parameters.

100ps molecular dynamic simulation of d(TATCACC)2

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage lambda is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H-bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.

A molecular dynamics study of the effect of G.T mispairs on the conformation of DNA in solution

The effect of G.T mispair incorporation into a double-helical environment was examined by molecular dynamics simulation. The 60-ps simulations performed on the two hexanucleotide duplexes d (G3C3)2 and d(G3TC2)2 included 10 Na+ counterions and first hydration shell waters. The resulting backbone torsional angle trajectories were analyzed to select time spans representative of conformational domains. The average backbone angles and helical parameters of the last time span for both duplexes are reported. During the simulation the hexamers retained B-type DNA structures that differed from typical A- or B-DNA forms. The overall helical structures for the two duplexes are vary similar. The presence of G.T mispairs did not alter the overall helical structure of the oligonucleotide duplex. Large propeller twist and buckle angles were obtained for both duplexes. The purine/pyrimidine crossover step showed a large decrease in propeller twist in the normal duplex but not in the mismatch duplex. Upon the formation of wobble mispairs in the mismatched duplex, the guanines moved into the minor groove and the thymines moved into the major groove. This helped prevent purine/purine clash and created a deformation in the relative orientation of the glycosidic bonds. It also exposed the free O4 of the thymines in the major groove and N2 of the guanines in the minor groove to interactions with solvent and counterions. These factors seemed to contribute to the apparently higher rigidity of the mismatched duplex during the simulation.

Molecular-mechanics modelling of drug-DNA structures; the effects of differing dielectric treatment on helix parameters and comparison with a fully solvated structural model

This study analyses the influence that the nature of the dielectric constant has on the final structures obtained from in vacuo molecular mechanics calculations on a drug-DNA complex and compares these structures with the energy minimised complex including explicit solvent molecules. Minimisations have been performed on a proflavine-decanucleotide structure, where the drug was intercalated at the d(CpG) site of the d(GpApTpApCpGpApTpApC) decamer duplex, using two expressions for the dielectric constant: a distance-independent, epsilon ij = EPS, and a distance-dependent, epsilon ij = EPS*Rij, form and for values of EPS from 1 to 8. Significantly different structures are obtained for the distance-independent and the distance-dependent expressions of the dielectric constant. The use of a distance-independent dielectric constant leads to distorted structures, which are very sensitive to slight changes in the value of EPS. The use of a distance-dependent dielectric constant leads to less distorted and more stable structures. The effects on helical parameters are analysed in detail. The structures obtained for different values of EPS (within the distance-dependent formalism) seem to converge for values of EPS equal to 4 or greater. Based on these results a distance-dependent form of the dielectric with an EPS value of 4 is recommended in order to produce reliable refined nucleic acid structures by molecular mechanics. These conclusions have been supported by molecular-mechanics minimisation of the same structure with the inclusion of explicit water molecules and counter-ions.

A note on graphing helical parameters of dynamics structure of DNA

A graphical procedure for analysis of helical parameters in dynamic structure of DNA is described. The performance of the procedure is illustrated by analysis of a 20 ps dynamics simulation of the non-self-complementary ninemer, 5'CAAACAGGA:5'TCCTGTTTG, which is a part of DNA from lacI gene. The dynamics trajectory of the duplex shows sequence-dependent fluctuations of helical parameters.

Molecular dynamics of an in vacuo model of duplex d(CGCGAATTCGCG) in the B-form based on the amber 3.0 force field

The characteristics of 100 ps of molecular dynamics (MD) on the DNA dodecamer d(CGCGAATTCGCG) at 300 K are described and investigated. The simulation is based on an in vacuo model of the oligomer and the AMBER 3.0 force field configured in the manner of Singh, U. C., S. J. Weiner, and P. A. Kollman, (1985, Proc. Natl. Acad. Sci. USA. 82:755-759). The analysis of the results was carried out using the "curves, dials, and windows" procedure (Ravishanker, G., S. Swaminathan, D. L. Beveridge, R. Lavery, and H. Sklenar. 1989. J. Biomol. Struct. Dyn. 6:669-699). The results indicate this dynamical model to be a provisionally stable double helix which lies at approximately 3.2 A rms deviation from the canonical B-form. There is, however, a persistent nonplanarity in the base pair orientations which resemble that observed in canonical A-DNA. The major groove width is seen to narrow during the course of the simulation and the minor groove expands, contravariant to the alterations in groove width seen in the crystal structure of the native dodecamer (Drew, H. R., R. M. Wing, T. Takano, C. Broka, S. Tanaka, I. Itakura, and R. E. Dickerson, 1981. Proc. Natl. Acad. Sci. USA. 78:2179-2183). The propeller twist in the bases, the sequence dependence of the base pair roll and aspects of bending in the helix axis are in some degree of agreement with the crystal structure. The patterns in DNA bending are observed to follow Zhurkin theory (Zhurkin, V. B. 1985. J. Biomol. Struct. Dyn. 2:785-804.). The relationship between the dynamical model and structure in solution is discussed.

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.

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.