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

Effect of cytarabine on the NMR structure of a model okazaki fragment from the SV40 genome

Okazaki fragments occur as intermediates during lagging strand DNA replication. Alterations in Okazaki fragment structure may contribute to the anticancer activities of nucleoside analogues such as cytarabine, a potent anti-leukemic agent that inhibits lagging strand replication. We have determined the solution structures for two model Okazaki fragments, [OKA] and [ARAC]. These sequences are derived from a frequent initiation site for primase during replication of the SV-40 viral genome. The sequence of [ARAC] differs from [OKA] only by substitution of cytarabine for one deoxycytidine. The structure of each model Okazaki fragment was elucidated using NMR spectroscopy and restrained molecular dynamics simulations. The solution structures of [OKA] and [ARAC] each consist of two distinct domains: a DNA duplex region (DDR) and an RNA-DNA hybrid duplex region (HDR). The DDR of [OKA] adopts geometry similar to B-form except for variations in helical parameters, especially twist and roll, which occur in the purine tract, increasing base overlap among the five consecutive purines. The helical axes for the DDR and HDR of [OKA] are bent 22 degrees relative to one another. Although the local structures for the DDR and HDR of [ARAC] are similar to those in [OKA] (root-mean-square deviation (rmsd) approximately 0.8, 1.7 A), the bending at the junction is different (41 degrees for [ARAC] vs 22 degrees for [OKA]). Increased helical bending of cytarabine-substituted Okazaki fragments may contribute to the propensity of cytarabine to inhibit elongation of the lagging strand during DNA replication, and in effecting anticancer activity.

Molecular dynamic simulations of environment and sequence dependent DNA conformations: the development of the BMS nucleic acid force field and comparison with experimental results

Molecular dynamic (MD) simulations using the BMS nucleic acid force field produce environment and sequence dependent DNA conformations that closely mimic experimentally derived structures. The parameters were initially developed to reproduce the potential energy surface, as defined by quantum mechanics, for a set of small molecules that can be used as the building blocks for nucleic acid macromolecules (dimethyl phosphate, cyclopentane, tetrahydrofuran, etc.). Then the dihedral parameters were fine tuned using a series of condensed phase MD simulations of DNA and RNA (in zero added salt, 4M NaCl, and 75% ethanol solutions). In the tuning process the free energy surface for each dihedral was derived from the MD ensemble and fitted to the conformational distributions and populations observed in 87 A- and B-DNA x-ray and 17 B-DNA NMR structures. Over 41 nanoseconds of MD simulations are presented which demonstrate that the force field is capable of producing stable trajectories, in the correct environments, of A-DNA, double stranded A-form RNA, B-DNA, Z-DNA, and a netropsin-DNA complex that closely reproduce the experimentally determined and/or canonical DNA conformations. Frequently the MD averaged structure is closer to the experimentally determined structure than to the canonical DNA conformation. MD simulations of A- to B- and B- to A-DNA transitions are also shown. A-DNA simulations in a low salt environment cleanly convert into the B-DNA conformation and converge into the RMS space sampled by a low salt simulation of the same sequence starting from B-DNA. In MD simulations using the BMS force field the B-form of d(GGGCCC)2 in a 75% ethanol solution converts into the A-form. Using the same methodology, parameters, and conditions the A-form of d(AAATTT)2 correctly converts into the B-DNA conformation. These studies demonstrate that the force field is capable of reproducing both environment and sequence dependent DNA structures. The 41 nanoseconds (nsec) of MD simulations presented in this paper paint a global picture which suggests that the DNA structures observed in low salt solutions are largely due to the favorable internal energy brought about by the nearly uniform screening of the DNA electrostatics. While the conformations sampled in high salt or mixed solvent environments occur from selective and asymmetric screening of the phosphate groups and DNA grooves, respectively, brought about by sequence induced ion and solvent packing.

Molecular dynamics simulations of DNA in solutions with different counter-ions

Molecular dynamics simulations of the [d(ATGCAGTCAG]2 fragment of DNA, in water and in the presence of three different counter-ions (Li+, Na+ and Cs+) are reported. Three-dimensional hydration structure and ion distribution have been calculated using spatial distribution functions for a detailed picture of local concentrations of ions and water molecules around DNA. According to the simulations, Cs+ ions bind directly to the bases in the minor groove, Na+ ions bind prevailing to the bases in the minor groove through one water molecule, whereas Li+ ions bind directly to the phosphate oxygens. The different behavior of the counter-ions is explained by specific hydration structures around the DNA and the ions. It is proposed how the observed differences in the ion binding to DNA may explain different conformational behavior of DNA. Calculated self-diffusion coefficients for the ions agree well with the available NMR data.

Selective binding of the TATA box-binding protein to the TATA box-containing promoter: analysis of structural and energetic factors

We report the results of an energy-based exploration of the components of selective recognition of the TATA box-binding protein (TBP) to a TATA box sequence that includes 1) the interaction between the hydrophobic Leu, Pro, and Phe residues of TBP with the TA, AT, AA, TT, and CG steps, by ab initio quantum mechanical calculations; and 2) the free energy penalty, calculated from molecular dynamics/potential of mean force simulations, for the conformational transition from A-DNA and B-DNA into the TA-DNA form of DNA observed in a complex with TBP. The GTAT, GATT, GAAT, and GTTT tetramers were explored. The results show that 1) the discrimination of TA, AT, AA, TT, or CG steps by TBP cannot rest on their interaction with the inserting Phe side chains; 2) the steric clash between the bulky and hydrophobic Pro and Leu residues and the protruding -NH2 group of guanine is responsible for the observed selectivity against any Gua-containing basepair; 3) the Pro and Leu residues cannot selectively discriminate among TA, AT, AA, or TT steps; and 4) the calculated energy required to achieve the TA-DNA conformation of DNA that is observed in the complex with TBP appears to be a key determinant for the observed selectivity against the AT, AA, and TT steps. The simulations also indicate that only the TA step can form a very efficient interbase hydrogen bond network in the TA-DNA conformation. Such an energetically stabilizing network is not achievable in the AA and TT steps. While it is viable in the AT step, structural constraints render the hydrogen bonding network energetically ineffective there.

Inherent DNA curvature and flexibility correlate with TATA box functionality

Four 1.5 ns molecular dynamics (MD) simulations were performed on the d(GCTATAAAAGGG).d(CCCTTTTATAGC) double helix dodecamer bearing the Adenovirus major late promoter TATA element and three iso-composition mutants for which physical and biochemical data are available from the same laboratory. Three of these DNA sequences experimentally induce tight binding with the TATA box binding protein (TBP) and induce high transcription rates; the other DNA sequence induces much lower TBP binding and transcription. The x-ray crystal structures have previously shown that the duplex DNA in DNA-TBP complexes are highly bent. We performed and analyzed MD simulations for these four DNAs, whose experimental structures are not available, in order to address the issue of whether inherent DNA structure and flexibility play a role in establishing these observed preferences. A comparison of the experimental and simulated results demonstrated that DNA duplex sequence-dependent curvature and flexibility play a significant role in TBP recognition, binding, and transcriptional activation.

Water molecules in DNA recognition II: a molecular dynamics view of the structure and hydration of the trp operator

The structure and hydration of the DNA duplex d-(AGCGTACTAGTACGCT)2 corresponding to the trp operator fragment used in the crystal structure of the half site complex (PDB entry 1TRR) was studied by a 1.4 ns molecular dynamics simulation in water. The simulation, starting from a B-DNA conformation, used a non-bonded cutoff of 1.4 nm with a reaction field correction and resulted in a stable trajectory. The average DNA conformation obtained was closer to the ones found in the crystal structures of the complexes (PDB entries 1TRO and 1TRR) than to the crystal structure of unbound trp operator (Nucleic Acid Database entry BDJ061). The DNA hydration was characterized in terms of hydrogen bond percentages and corresponding residence times. The residence times of water molecules within 0.35 nm of the DNA non-exchangeable protons were calculated for comparison with NMR measurements of intermolecular water-DNA NOEs and nuclear magnetic relaxation dispersion measurements. No significant difference was found between major and minor groove hydration. The DNA donors and acceptors were hydrogen bonded to water molecules for 77(+/-19)% of the time on average. The average residence time of the hydrogen bonded water molecules was 11(+/-11) ps with a maximum of 223 ps. When all water molecules within NOE distance (0.35 nm) of non-exchangeable protons were considered, the average residence times increased to an average of 100(+/-4) ps and a maximum of 608 ps. These results agree with the experimental NMR results of Sunnerhagen et al. which did not show any evidence for water molecules bound with more than 1 ns residence time on the DNA surface. The exchange of hydration water from the DNA occurred in the major groove primarily through direct exchange with the bulk solvent, while access to and from the minor groove frequently proceeded via pathways involving ribose O3' and O4' and phosphate O2P oxygen atoms. The most common water diffusion pathways in the minor groove were perpendicular to the groove direction. In general, water molecules visited only a limited number of sites in the DNA grooves before exiting. The hydrogen bonding sites, where hydrogen bonds could be formed with donor and acceptor groups of the DNA, were filled with water molecules with an average B-factor value of 0.58 mn2. No special values were observed at any of the sites, where water molecules were observed both in the trp repressor/operator co-crystals and in the crystal structure of unbound DNA.

Molecular dynamics and continuum solvent studies of the stability of polyG-polyC and polyA-polyT DNA duplexes in solution

Molecular dynamics simulation in explicit solvent and continuum solvent models are applied to investigate the relative stability of A- and B-form helices for two DNA sequences, dA10-dT10 and dG10-dC10 in three structural forms. One structural form is based on an unrestrained molecular dynamics (MD) trajectory starting from a canonical B-DNA structure, the second is based on a MD trajectory starting in a canonical B-DNA structure with the sugars constrained to be C2'-endo and the third simulation started from a canonical A-DNA structure with the sugars constrained to C3'-endo puckers. For the energetic analysis, structures were taken as snapshots from nanosecond length molecular dynamics simulations computed in a consistent fashion in explicit solvent, applying the particle mesh Ewald method and the Cornell et al. force field. The electrostatic contributions to solvation free energies are computed using both a finite-difference Poisson-Boltzmann model and a pairwise Generalized Born model. The non-electrostatic contributions to the solvation free energies are estimated with a solvent accessible surface area dependent term. To estimate the gas phase component of the relative free energy between the various structures, the mean solute internal energies (determined with the Cornell et al. molecular mechanics potential including all pairwise interactions within the solute) and estimates of the solute entropy (using a harmonic approximation) were used. Consistent with experiment, the polyG-polyC (GC) structures are found to be much more A-phillic than the polyA-polyT (AT) structures, the latter being quite A-phobic. The dominant energy components responsible for this difference comes from the internal and van der Waal energies. A perhaps less appreciated difference between the GC and AT rich sequences is suggested by the calculated salt dependence which demonstrates a significantly enhanced ability to drive GC rich sequences towards an A-form structure compared to AT rich sequences. In addition to being A-phobic, the AT structure also has a noticably larger helical repeat than GC and other mixed sequence duplexes, consistent with experiment. Analysis of the average solvent density from the trajectories shows hydration patterns in qualitative agreement with experiment and previous theoretical treatments.

Molecular dynamics simulations of an oligonucleotide duplex with adenine tracts phased by a full helix turn

A theoretical model of a DNA oligonucleotide duplex featuring A-tracts phased by a full helix turn is developed based on molecular dynamics computer simulation. The extent to which this model agrees with relevant experimental data on axis bending and the relationship of A-tracts to bending and other aspects of helix morphology is investigated. Specifically, a series of nanosecond-level molecular dynamics (MD) simulations have been carried out for the 25 bp duplex d(ATAGGCAAAAAATAGGCAAAAATGG) at various concentrations of saline solution. A 30 base-pair sequence composed of three 10 bp repeats of the BamHI recognition sequence ligated together, d(CGGGATCCCG. CGGGATCCCG.CGGGATCCCG), was simulated as a control. The MD was carried out using the AMBER 4.1 suite of programs, and utilized the Cornell et al. force-field with the electrostatic boundary conditions treated by the particle-mesh Ewald summation protocol. The MD results show that at a concentration of 60 mM KCl, 10 mM MgCl2 added salt plus minimal neutralizing cations, the MD model exhibits concerted axis bending to the extent of 15.5 degrees per A-tract. This compares favorably with the bending per turn of 17 to 21 degrees inferred from cyclization experiments. The MD model also exhibits a progressive 5' to 3' narrowing of the minor-groove region of A-tracts, a feature inferred from DNA footprinting experiments. Analysis of the dynamic structure of the MD models shows that the origin of the bending follows a junction-type bending model with an admixture of mixed sequence effects, with A-tracts relatively straight, as in oligonucleotide crystal structures of sequences containing A-tracts. The results are shown to be sensitive to environmental conditions: MD on d(ATAGGCAAAAAATAGGCAAAAATGG) in neutralizing Na+ buffer results in markedly reduced curvature, and the removal of Mg2+ measurably affects bending. Carrying out the simulations at experimental salt conditions appears to be essential to obtain an accurate account of the experimentally observed bending.

Nuclear magnetic resonance structure of d(GCATATGATAG). d(CTATCATATGC): a consensus sequence for promoters recognized by sigma K RNA polymerase

The three-dimensional structure of d(GCATATGATAG).d(CTATCATATGC), from the promoter region of a gene regulating sporulation in Bacillus subtilis mother cells, was determined utilizing two-dimensional nuclear Overhauser effect (2D NOE) and double-quantum-filtered COSY (2QF-COSY) spectra. To minimize the effect of methods used to obtain restraints and refine structure, several variables were studied. Interproton distance bounds were calculated very conservatively by running the complete relaxation matrix program MARDIGRAS hundreds of times using 2D NOE spectra for exchangeable and for nonexchangeable protons at different mixing times, assuming different overall correlation times and different starting structures. The 435 distance restraints were used with two different structural refinement methods: restrained molecular dynamics (rMD) and restrained Monte Carlo calculations (rMC). Refinement using different procedures and starting structures resulted in essentially the same structure (<0.8 A rmsd), indicating that the structure is defined by experimental restraints and not the refinement method or variables used. R factors indicate the structures fit the experimental NOE data very well. Some helical parameters, notably large negative X displacement, are characteristic of A-DNA, but others are characteristic of B-DNA. As with TG.CA steps in other duplex DNA sequences studied in our laboratory, the two TG.CA steps have a positive roll, with T6-G7 exhibiting the largest, and consequently a bent helix axis. The converged structure represents a time-averaged structure. However, multiple conformations, especially in deoxyriboses, were evident from vicinal coupling constants obtained from quantitative simulations of 2QF-COSY cross-peaks and from persistent inconsistencies in experimental distances due to nonlinear conformational averaging.

Structure of the metal-ion-activated diphtheria toxin repressor/tox operator complex

The virulent phenotype of the pathogenic bacterium Corynebacterium diphtheriae is conferred by diphtheria toxin, whose expression is an adaptive response to low concentrations of iron. The expression of the toxin gene (tox) is regulated by the repressor DtxR, which is activated by transition metal ions. X-ray crystal structures of DtxR with and without (apo-form) its coordinated transition metal ion have established the general architecture of the repressor, identified the location of the metal-binding sites, and revealed a metal-ion-triggered subunit-subunit 'caliper-like' conformational change. Here we report the three-dimensional crystal structure of the complex between a biologically active Ni(II)-bound DtxR(C102D) mutant, in which a cysteine is replaced by an aspartate at residue 102, and a 33-base-pair DNA segment containing the toxin operator toxO. This structure shows that DNA interacts with two dimeric repressor proteins bound to opposite sides of the tox operator. We propose that a metal-ion-induced helix-to-coil structural transition in the amino-terminal region of the protein is partly responsible for the unique mode of repressor activation by transition metal ions.

Crystal structure of lambda-Cro bound to a consensus operator at 3.0 A resolution

The structure of the Cro protein from bacteriophage lambda in complex with a 19 base-pair DNA duplex that includes the 17 base-pair consensus operator has been determined at 3.0 A resolution. The structure confirms the large changes in the protein and DNA seen previously in a crystallographically distinct low-resolution structure of the complex and, for the first time, reveals the detailed interactions between the side-chains of the protein and the base-pairs of the operator. Relative to the crystal structure of the free protein, the subunits of Cro rotate 53 degrees with respect to each other on binding DNA. At the same time the DNA is bent by 40 degrees through the 19 base-pairs. The intersubunit connection includes a region within the protein core that is structurally reminiscent of the "ball and socket" motif seen in the immunoglobulins and T-cell receptors. The crystal structure of the Cro complex is consistent with virtually all available biochemical and related data. Some of the interactions between Cro and DNA proposed on the basis of model-building are now seen to be correct, but many are different. Tests of the original model by mutagenesis and biochemical analysis corrected some but not all of the errors. Within the limitations of the crystallographic resolution it appears that operator recognition is achieved almost entirely by direct hydrogen-bonding and van der Waals contacts between the protein and the exposed bases within the major groove of the DNA. The discrimination of Cro between the operators OR3 and OR1, which differ in sequence at just three positions, is inferred to result from a combination of small differences, both favorable and unfavorable. A van der Waals contact at one of the positions is of primary importance, while the other two provide smaller, indirect effects. Direct hydrogen bonding is not utilized in this distinction.

Progressive DNA bending is made possible by gradual changes in the torsion angle of the glycosyl bond

Structural comparisons have led to the suggestion that the conformational rearrangement that would be required to change A-DNA into the TA-DNA form of DNA observed in the complex with the TATA box binding protein (TBP) could be completed by modifying only the value of the glycosyl bond chi by approximately 45 degrees. The lack of a high number of crystal structures of this type makes it difficult to conclude whether a smooth transition from A-DNA to TA-DNA can occur without disrupting at any point either the Watson-Crick base pairing or the A-DNA conformation of the backbone. To explore the possibility of such a smooth transition, constrained molecular dynamics simulations were carried out for the double-stranded dodecamer d(GGTATATAAAAC), in which a transition from A-DNA to TA-DNA was induced by modifying only the chi angle values. The results demonstrate the feasibility of a continuous path in the A-DNA to TA-DNA transition. Varying extents of DNA curvature are also attainable, by maintaining the A-DNA backbone structure and Watson-Crick hydrogen bonding while changing the chi angle value smoothly from that in A-DNA to one corresponding to B-DNA.

The intrinsic curvature of a 51 bp K-DNA fragment of Leishmania tarentolae: a molecular model

DNA intrinsic structure and curvature is a subject of debate because of the importance of these attributes in processes such as DNA packaging, transcription, and gene regulation. X-ray crystallography of DNA single crystals has provided a wealth of information about the local, short range conformational features of DNA. On the other hand, gel electrophoresis analysis of DNA has not only uncovered the macroscopic curvature of DNA but it also provides most of the available data on DNA intrinsic curvature. However, gel electrophoresis can not identify features of DNA structure at the nucleotide or atomic level. In order to address the problem of DNA intrinsic curvature in an attempt to bridge the gap between X-ray crystallography and gel electrophoresis, we use the computational method of molecular dynamics (MD). In this study, we report the results of 2.0 ns MD simulations on a 51 bp fragment of the K-DNA of Leishmania tarentolae containing several A-tracts. The K-DNA double helix is very stable and remains in an intermediate state between the canonical A and B forms of the duplex. The magnitude of global curvature (75 degrees) agrees well with the experimental estimate (72 degrees) available. Analysis of local (every base triplet) and sublocal (every helix turn) curvature shows that the 51 bp K-DNA fragment has curvature features also present in the Wedge, Junction and Calladine's models of DNA intrinsic curvature. We further characterize the flexibility of individual nucleotides in the molecule and find the sugar flexibility within the A-tracts to be strongly correlated with the pattern of A-tract cleavage by the hydroxyl radical. Differential curvature and flexibility at the 5' and 3'junctions between A-tracts and general-sequence DNA are found to modulate the global curvature of the K-DNA fragment.

Solution structure of a highly stable DNA duplex conjugated to a minor groove binder

The tripeptide 1,2-dihydro-(3 H )-pyrrolo[3,2- e ]indole-7-carboxylate (CDPI3) binds to the minor groove of DNA with high affinity. When this minor groove binder is conjugated to the 5'-end of short oligonucleotides the conjugates form unusually stable hybrids with complementary DNA and thus may have useful diagnostic and/or therapeutic applications. In order to gain an understanding of the structural interactions between the CDPI3minor groove binding moiety and the DNA, we have determined and compared the solution structure of a duplex consisting of oligodeoxyribonucleotide 5'-TGATTATCTG-3' conjugated at the 5'-end to CDPI3 and its complementary strand to an unmodified control duplex of the same sequence using nuclear magnetic resonance techniques. Thermal denaturation studies indicated that the hybrid of this conjugate with its complementary strand had a melting temperature that was 30 degrees C higher compared with the unmodified control duplex. Following restrained molecular dynamics and relaxation matrix refinement, the solution structure of the CDPI3-conjugated DNA duplex demonstrated that the overall shape of the duplex was that of a straight B-type helix and that the CDPI3moiety was bound snugly in the minor groove, where it was stabilized by extensive van der Waal's interactions.

Structure prediction of a complex between the chromosomal protein HMG-D and DNA

Non-histone chromosomal proteins are an important part of nuclear structure and function due to their ability to interact with DNA to form and modulate chromatin structure and regulate gene expression. However, the understanding of the function of chromosomal proteins at the molecular level has been hampered by the lack of structures of chromosomal protein-DNA complexes. We have carried out a molecular dynamics modeling study to provide insight into the mode of DNA binding to the chromosomal HMG-domain protein, HMG-D. Three models of a complex of HMG-D bound to DNA were derived through docking the protein to two different DNA fragments of known structure. Molecular dynamics simulations of the complexes provided data indicating the most favorable model. This model was further refined by molecular dynamics simulation and extensively analyzed. The structure of the corresponding HMG-D-DNA complex exhibits many features seen in the NMR structures of the sequence-specific HMG-domain-DNA complexes, lymphoid enhancer factor 1 (LEF-1) and testis determining factor (SRY). The model reveals differences from these known structures that suggest how chromosomal proteins bind to many different DNA sequences with comparable affinity.

Observation of an A-DNA to B-DNA transition in a nonhelical nucleic acid hairpin molecule using molecular dynamics

One of the truly challenging problems for molecular dynamics (MD) simulations is demonstrating that the trajectories can sample not only in the vicinity of an experimentally determined structure, but also that the trajectories can find the correct experimental structure starting from some other structure. Frequently these transitions to the correct structure require that the simulations overcome energetic barriers to conformational change. Here we present unrestrained molecular dynamics simulations of the DNA analogs of the RNA 5'-GGACUUCGGUCC-3' hairpin tetraloop. In one simulation we have used deoxyuracil residues, and in the other we have used the native DNA deoxythymine residues. We demonstrate that, on a nanosecond time scale, MD is able to simulate the transitions of both of the A-DNA stems to B-DNA stems within the constraints imposed by the four-base loop that caps the helix. These results suggest that we are now in a position to use MD to address the nature of sequence-dependent structural effects in nonduplex DNA structures.

A 5-nanosecond molecular dynamics trajectory for B-DNA: analysis of structure, motions, and solvation

We report the results of four new molecular dynamics (MD) simulations on the DNA duplex of sequence d(CGCGAATTCGCG)2, including explicit consideration of solvent water, and a sufficient number of Na+ counterions to provide electroneutrality to the system. Our simulations are configured particularly to characterize the latest MD models of DNA, and to provide a basis for examining the sensitivity of MD results to the treatment of boundary conditions, electrostatics, initial placement of solvent, and run lengths. The trajectories employ the AMBER 4.1 force field. The simulations use particle mesh Ewald summation for boundary conditions, and range in length from 500 ps to 5.0 ns. Analysis of the results is carried out by means of time series for conformationalm, helicoidal parameters, newly developed indices of DNA axis bending, and groove widths. The results support a dynamically stable model of B-DNA for d(CGCGAATTCGCG)2 over the entire length of the trajectory. The MD results are compared with corresponding crystallographic and NMR studies on the d(CGCGAATTCGCG)2 duplex, and placed in the context of observed behavior of B-DNA by comparisons with the complete crystallographic data base of B-form structures. The calculated distributions of mobile solvent molecules, both water and counterions, are displayed. The calculated solvent structure of the primary solvation shell is compared with the location of ordered solvent positions in the corresponding crystal structure. The results indicate that ordered solvent positions in crystals are roughly twice as structured as bulk water. Detailed analysis of the solvent dynamics reveals evidence of the incorporation of ions in the primary solvation of the minor groove B-form DNA. The idea of localized complexation of otherwise mobile counterions in electronegative pockets in the grooves of DNA helices introduces an additional source of sequence-dependent effects on local conformational, helicoidal, and morphological structure, and may have important implications for understanding the functional energetics and specificity of the interactions of DNA and RNA with regulatory proteins, pharmaceutical agents, and other ligands.

Insight into the stabilization of A-DNA by specific ion association: spontaneous B-DNA to A-DNA transitions observed in molecular dynamics simulations of d[ACCCGCGGGT]2 in the presence of hexaamminecobalt(III)

BACKGROUND

Duplex DNA is more than a simple information carrier. The sequence-dependent structure and its inherent deformability, in concert with the subtle modulating effects of the environment, play a crucial role in the regulation and packaging of DNA. Recent advances in force field and simulation methodologies allow molecular dynamics simulations to now represent the specific effects of the environment. An understanding of the environmental dependence of DNA structure gives insight into how histones are able to package DNA, how various proteins are able to bind and modulate nucleic acid structure and will ultimately aid the design of molecules to package DNA for more effective gene therapy.

RESULTS

Molecular dynamics simulations of d[ACCCGCGGGT]2 in solution in the presence of hexaamminecobalt(III) [Co(NH3)6(3+)] show stabilization of A-DNA and spontaneous B-DNA to A-DNA transitions, which is consistent with experimental results from NMR and Raman spectroscopic and X-ray crystallographic studies. In the absence of Co(NH3)6(3+), A-DNA to B-DNA transitions are observed instead. In addition to their interaction with the guanines in the major groove, Co(NH3)6(3+) ions bridge opposing strands in the bend across the major groove, probably stabilizing A-DNA.

CONCLUSIONS

The simulation methods and force fields have advanced to a sufficient level that some representation of the environment can be seen in nanosecond length molecular dynamics simulations. These simulations suggest that, in addition to the general explanation of A-DNA stabilization by dehydration, hydration and ion association in the major groove stabilize A-DNA.

The crystal structure of an RNA oligomer incorporating tandem adenosine-inosine mismatches

The X-ray crystallographic structure of the RNA duplex [r(CGCAIGCG)]2 has been refined to 2.5 A. It shows a symmetric internal loop of two non-Watson-Crick base pairs which form in the middle of the duplex. The tandem A-I/I-A pairs are related by a crystallographic two-fold axis. Both A(anti)-I(anti) mismatches are in a head-to-head conformation forming hydrogen bonds using the Watson-Crick positions. The octamer duplexes stack above one another in the cell forming a pseudo-infinite helix throughout the crystal. A hydrated calcium ion bridges between the 3'-terminal of one molecule and the backbone of another. The tandem A-I mismatches are incorporated with only minor distortion to the backbone. This is in contrast to the large helical perturbations often produced by sheared G-A pairs in RNA oligonucleotides.

Molecular dynamics simulation study of DNA dodecamer d(CGCGAATTCGCG) in solution: conformation and hydration

A molecular dynamics simulation of the dodecamer duplex d(CGCGAATTCGCG) using the particle mesh Ewald sum assumed a B-conformation remarkably close to the observed X-ray structure. The Ewald summation method effectively eliminates the usual "cut-off" of long-range interactions and allowed us to evaluate the full effect of the electrostatic forces. This simulation showed remarkable agreement with the Dickerson X-ray structure in both average structure and B-factors; within the EcoRI site itself, the rms deviation between the average theoretical and observed structures was 1.1 A. The width of the minor groove fluctuated between a wide and narrow configuration with the latter corresponding closely to the X-ray structure. The simulation also suggested a strong sequence-dependent signature on the minor groove width in both wide and narrow conformers. Hydration shells in both the major and minor grooves were observed. The "spine of hydration" in the minor groove was clear. In the major groove the first hydration shell appears to be a ribbon-like structure that reproduces the principal features of observed X-ray structures; subtle variations of this hydration pattern suggest sequence dependencies. Sequence-dependent features were also examined for helical and other geometric parameters. The successful reproduction of many experimentally observed fine structural features shows that the Ewald summation significantly improves the fidelity of the calculations.

A molecular level picture of the stabilization of A-DNA in mixed ethanol-water solutions

Advances in computer power, methodology, and empirical force fields now allow routine "stable" nanosecond-length molecular dynamics simulations of DNA in water. The accurate representation of environmental influences on structure remains a major, unresolved issue. In contrast to simulations of A-DNA in water (where an A-DNA to B-DNA transition is observed) and in pure ethanol (where disruption of the structure is observed), A-DNA in approximately 85% ethanol solution remains in a canonical A-DNA geometry as expected. The stabilization of A-DNA by ethanol is likely due to disruption of the spine of hydration in the minor groove and the presence of ion-mediated interhelical bonds and extensive hydration across the major groove.

Does TATA matter? A structural exploration of the selectivity determinants in its complexes with TATA box-binding protein

The binding of the TATA box-binding protein (TBP) to a TATA sequence in DNA is essential for eukaryotic basal transcription. TBP binds in the minor groove of DNA, causing a large distortion of the DNA helix. Given the apparent stereochemical equivalence of AT and TA basepairs in the minor groove, DNA deformability must play a significant role in binding site selection, because not all AT-rich sequences are bound effectively by TBP. To gain insight into the precise role that the properties of the TATA sequence have in determining the specificity of the DNA substrates of TBP, the solution structure and dynamics of seven DNA dodecamers have been studied by using molecular dynamics simulations. The analysis of the structural properties of basepair steps in these TATA sequences suggests a reason for the preference for alternating pyrimidine-purine (YR) sequences, but indicates that these properties cannot be the sole determinant of the sequence specificity of TBP. Rather, recognition depends on the interplay between the inherent deformability of the DNA and steric complementarity at the molecular interface.

Binding of the estrogen receptor to DNA. The role of waters

Molecular dynamics simulations are carried out to investigate the binding of the estrogen receptor, a member of the nuclear hormone receptor family, to specific and non-specific DNA. Two systems have been simulated, each based on the crystallographic structure of a complex of a dimer of the estrogen receptor DNA binding domain with DNA. One structure includes the dimer and a consensus segment of DNA, ds(CCAGGTCACAGTGACCTGG); the other structure includes the dimer and a nonconsensus segment of DNA, ds(CCAGAACACAGTGACCTGG). The simulations involve an atomic model of the protein-DNA complex, counterions, and a sphere of explicit water with a radius of 45 A. The molecular dynamics package NAMD was used to obtain 100 ps of dynamics for each system with complete long-range electrostatic interactions. Analysis of the simulations revealed differences in the protein-DNA interactions for consensus and nonconsensus sequences, a bending and unwinding of the DNA, a slight rearrangement of several amino acid side chains, and inclusion of water molecules at the protein-DNA interface region. Our results indicate that binding specificity and stability is conferred by a network of direct and water mediated protein-DNA hydrogen bonds. For the consensus sequence, the network involves three water molecules, residues Glu-25, Lys-28, Lys-32, Arg-33, and bases of the DNA. The binding differs for the nonconsensus DNA sequence in which case the fluctuating network of hydrogen bonds allows water molecules to enter the protein-DNA interface. We conclude that water plays a role in furnishing DNA binding specificity to nuclear hormone receptors.

Theoretical studies of an exceptionally stable RNA tetraloop: observation of convergence from an incorrect NMR structure to the correct one using unrestrained molecular dynamics

We report on the results of five independent and unrestrained molecular dynamics simulations of an RNA tetraloop, r(GGACUUCGGUCC), and its related structures with the loop UUCG sugars changed to deoxyribose. Two separate NMR structures have been reported for the loop portion of this molecule, with the second refinement resulting in a slightly different and more accurate conformation for the loop. The root-mean-square deviation (RMSd) between the two NMR structures, for the loop portions only, is 2.5 A. Our simulations, starting from the two NMR structures, demonstrate that this tetraloop is a very stable and rigid structure with both nanosecond length simulations staying very close to the initial structures. Additionally, both simulations preserved most, if not all, of the NMR-derived interactions and violated very few of the nuclear Overhauser effect (NOE)-derived distances used in the structure refinements. However, when the two NMR structures were simulated with deoxyriboses in the loops instead of the native riboses, the flexibility of the systems increased and we observed a conversion from the incorrect to the correct loop conformation in the simulation which started in the incorrect loop conformation. When the riboses were subsequently re-introduced back into the structure which underwent the conversion, the agreement between this simulation and the one starting from the correct NMR structure was a remarkably low 0.5 A, demonstrating an almost complete convergence from the incorrect to the correct structure using unrestrained molecular dynamics.

Solution structures of 5-fluorouracil-substituted RNA duplexes containing G-U wobble base pairs

The structures and stabilities of three RNA duplexes that differed only in the position of 5-fluorouridine (FUrd) substitution were elucidated using NMR spectroscopy and UV hyperchromicity studies to determine if FUrd substitution altered the structure or stability of RNA duplexes that contained G-U base pairs. The duplexes investigated corresponded to the region of the U4-U6 snRNA complex that contained the 5' terminus of U4 snRNA. The control duplex contained a G-U wobble base pair and also a G-A mismatched base pair. FUrd was substituted in one duplex at the G-U wobble base pair and in the second duplex at an A-U base pair adjacent to the wobble base pair. FUrd substitution slightly destabilized the duplex that contained a G-FU base pair but stabilized the duplex that contained an A-FU base pair. NOESY spectra were used to determine interproton distances, and these distance constraints were used in a restrained molecular dynamics protocol to determine the three-dimensional structures of these RNA duplexes. Analyses of helical parameters, backbone torsion angles, and rms deviations between the final structures revealed no systematic differences due to FUrd substitution in RNA duplexes that contained G-U base pairs. The G-FU base pair adopted wobble geometry, while the G-A mismatch formed a sheared base pair. NOESY spectra in H2O solution revealed the imino 1H from FUrd exchanged more rapidly with solvent than did the Urd imino 1H but did not show the G-FU base pair adopted an ionized structure. Reduced stacking occurred for the G-FU base pair relative to the G-U base pair in the time-averaged structure, and this, rather than ionization of the base pair, was responsible for the slight destabilization of the duplex that contained the G-FU base pair.

Solution structure of a DNA decamer containing the antiviral drug ganciclovir: combined use of NMR, restrained molecular dynamics, and full relaxation matrix refinement

The nucleoside analog 9-[(1,3-dihydroxy-2-propoxy)methyl]guanine (ganciclovir, DHPG) is an antiviral drug that is used in the treatment of a variety of herpes viruses in immunocompromised patients and in a gene therapy protocol that has shown promising activity for the treatment of cancer. To probe the structural effects of ganciclovir when incorporated into DNA, we determined and compared the solution structure of a modified ganciclovir-containing decamer duplex [d(CTG)(ganciclovir)d(ATCCAG)]2 and a control duplex d[(CTGGATCCAG)]2 using nuclear magnetic resonance techniques. 1H and 31P resonances in both duplexes were assigned using a combination of 2-D 1H and 31P NMR experiments. Proton-proton distances determined from NOESY data and dihedral angles determined from DQF-COSY data were used in restrained molecular dynamics simulations starting from canonical A- and B-form DNA models. Both the control and ganciclovir sets of simulations converged to B-type structures. These structures were subjected to full relaxation matrix refinement to produce final structures that were in excellent agreement with the observed NOE intensities. Examination of the final ganciclovir-containing structures reveals that the base of the ganciclovir residue is hydrogen bonded to its complementary dC and is stacked in the helix; in fact, the base of ganciclovir exhibits increased stacking with the 5' base relative to the control. Interestingly, some of the most significant distortions in the structures occur 3' to the lesion site, including a noticeable kink in the sugar-phosphate backbone at this position. Further examination reveals that the backbone conformation, sugar pucker, and glycosidic torsion angle of the residue 3' to the lesion site all indicate an A-type conformation at this position. A possible correlation of these structural findings with results obtained from earlier biochemical studies will be discussed.

How hormone receptor-DNA binding affects nucleosomal DNA: the role of symmetry

Molecular dynamics simulations have been employed to determine the optimal conformation of an estrogen receptor DNA binding domain dimer bound to a consensus response element, ds(AGGTCACAGTGACCT), and to a nonconsensus response element, ds(AGAACACAGTGACCT). The structures simulated were derived from a crystallographic structure and solvated by a sphere (45-A radius) of explicit water and counterions. Long-range electrostatic interactions were accounted for during 100-ps simulations by means of a fast multipole expansion algorithm combined with a multiple time-step scheme in the molecular dynamics package NAMD. The simulations demonstrate that the dimer induces a bent and underwound (10.7 bp/turn) conformation in the DNA. The bending reflects the dyad symmetry of the receptor dimer and can be described as an S-shaped curve in the helical axis of DNA when projected onto a plane. A similar bent and underwound conformation is observed for nucleosomal DNA near the nucleosome's dyad axis that reflects the symmetry of the histone octamer. We propose that when a receptor dimer binds to a nucleosome, the most favorable dimer-DNA and histone-DNA interactions are achieved if the respective symmetry axes are aligned. Such positioning of a receptor dimer over the dyad of nucleosome B in the mouse mammary tumor virus promoter is in agreement with experiment.

Calculated distortions of duplex DNA by a cis, syn cyclobutane thymine dimer are unaffected by a 3' TpA step

Molecular dynamics simulations were performed on the duplex DNA dodecamers d(CGCGAA TT CGCG): d(CGCGAATTCGCG) and d(GCACGAA TT AAG): d(CTTAATTCGTGC), where TT denotes a cis, syn cyclobutane thymine dimer. The constant temperature and pressure algorithm of the AMBER 4.1 molecular-modeling package was used with explicit water and counterions, periodic boundary conditions and electrostatic interactions evaluated by the particle-mesh Ewald method. Results were analyzed by the CURVES algorithm and its implementation in DIALS and WINDOWS. Calculated distortions of DNA structure by the thymine dimer were qualitatively and quantitatively similar for the two sequences. Despite the enhanced flexibility of the native TpA dinucleotide step, major deviations from the B-DNA values of helicoidal parameters were found only at the Ap and p dinucleotide steps in both sequences. Only the AT base pairs of the two sequences that contain the 5' thymine of the dimers exhibited weakened Watson-Crick hydrogen bonds and anomalous stretching. Hence, we conclude that the pattern of structural perturbations responsible for recognition of cis, syn thymine dimers by repair enzymes is not sensitive to their sequence context.

Solution structures of 5-fluorouracil-substituted DNA and RNA decamer duplexes

The structures in solution of eight oligonucleotide duplexes each containing either zero, one, or two 5-fluorodeoxyuridine (FdUrd) or 5-fluorouridine (FUrd) nucleosides were determined by the combined use of NMR spectroscopy, restrained molecular dynamics, and full relaxation matrix refinement to determine how FdUrd and FUrd substitution affects the structure of duplex DNA and RNA and to establish whether structural differences due to FdUrd and FUrd substitution in nucleic acids may be responsible, in part, for the biological effects of the anticancer drug 5-fluorouracil (FUra). The nucleic acid directed effects of FUra include induction of single-strand breaks in duplex DNA and altered processing of pre-mRNA and rRNA. Four self-complementary oligodeoxyribonucleotide sequences were prepared and studied as duplexes in aqueous solution: (5' dGCGAAUUCGC)2, (5' dGCGAAUFCGC)2, (5' dGCGAAFUCGC)2, and (5' dGCGAAFFCGC)2. The corresponding oligoribonucleotide sequences (5' rGCGAAUUCGC)2, (5' rGCGAAUFCGC)2, (5' rGCGAAFUCGC)2, and (5' rGCGAAFFCGC)2 were also prepared and studied. The helical parameters for the structures of these eight duplexes were analyzed to determine how substitution of FdUrd and FUrd affects the three-dimensional structures of duplex DNA and RNA. FdUrd substitution affects the base roll angle at the site of FdUrd substitution, causing the helical axis of FdUrd-substituted DNA duplexes to be bent compared to the nonsubstituted duplex. A-FUrd base pairs show substantial RMS deviations from A-Urd base pairs in all three of the RNA duplexes substituted with FUrd. Bending of the helical axis due to FdUrd substitution may contribute to the occurrence of single-strand breaks in duplex DNA while the altered structures of A-FUrd base pairs may affect RNA-RNA and RNA-protein recognition.

Effect of adenine methylation on the structure and dynamics of TpA steps in DNA: NMR structure determination of [d(CGAGGTTTAAACCTCG)]2 and its A9-methylated derivative at 750 MHz

At TpA steps in DNA, the adenine base experiences exceptionally large amplitude (20 degrees-50 degrees) and slow (10 ms-1 microsecond) motion [Kennedy et al. (1993) Biochemistry 32, 8022-8035] which has been correlated with transitions between multiple conformational states [Lefevre et al. (1985) FEBS Lett. 190, 37-40]. The base dynamics can be detected in one-dimensional 1H NMR spectra as excess line width of the aromatic proton resonances. The magnitude of the excess line width is temperature dependent and reaches a maximum at some temperature. In order to better understand the origin of the dynamics, we have studied the effect of N6-methylation of the TpA adenine on both the line widths and its local structure. Here, solution-state 500 and 750 MHz 1H NMR data collected on [d(CGAGGTTTAAACCTCG)]2 show that the excess line width of the TpA adenine-H2 is diminished when the TpA adenine is N6-methylated and that the line width no longer experiences a maximum as the temperature is varied. The resonances sharpen upon methylation because the amplitude of base motion is restricted due to steric effects and due to other structural changes at the TpA site. Additionally, both the TpA adenine-H8 and the exchangeable imino resonance of thymine at the TpA step were also found to have excess line width that is diminished upon N6-methylation. In order to elucidate the structural features responsible for TpA base dynamics, solution-state NMR structures of [d(CGAGGTTTAAACCTCG)]2 and its A9 N6-methylated derivative were determined at 750 MHz. Comparison of the structures shows that poor base stacking at the TpA step may contribute to, or be the origin of, its base dynamics.

Molecular dynamics study of glucocorticoid receptor-DNA binding

Molecular dynamics simulations have been conducted to investigate the binding of the glucocorticoid receptor (GR) dimer to DNA. For this purpose simulations of the complex formed by a DNA segment and a dimer of GR-DNA binding domains (GR-DBD) have been carried out, employing an available X-ray structure. A second set of simulations was based on this structure as well, except that the DNA segment was altered to the consensus glucocorticoid response element (GRE). Simulations of a single GR-DBD and of the uncomplexed GRE served as controls. For the simulations, each system was encapsulated in an ellipsoid of water. Protein-DNA interactions, dimer interactions, and DNA structural parameters were analyzed for each system and compared. The consensus GRE is found to yield more favorable and symmetric interactions between the GR-DBDs and the GRE, explaining the ability of the GR dimer to recognize this DNA segment. Further analysis focused on deformations of the DNA that are induced by the binding of GR. The deformations observed involve a 35 degree bend of the DNA, an unwinding, and a displacement of the helical axis. These deformations are consistent with a mechanism for transcriptional regulation that involves a change of nucleosome packing upon GR binding. Significant protein-protein and protein-DNA interactions, both direct and water mediated, develop due to the deformations of the GRE and are indicative of an increased recognition achieved through DNA deformation. The interactions include direct interactions between the GRE and glycine-458 and serine-459, side groups which differentiate GR from other members of the nuclear hormone receptor family.

Potential of mean force calculations of the stacking-unstacking process in single-stranded deoxyribodinucleoside monophosphates

The free energy of the stacking-unstacking process of deoxyribodinucleoside monophosphates in aqueous solution has been investigated by potential of mean force calculations along a reaction coordinate, defined by the distance between the glycosidic nitrogen atoms of the bases. The stacking-unstacking process of a ribodinucleoside monophosphate was observed to be well characterized by this coordinate, which has the advantage that it allows for a dynamical backbone and flexible bases. All 16 naturally occurring DNA dimers composed of the adenine, cytosine, guanine, or thymine bases in both the 5' and the 3' positions were studied. From the free-energy profiles we observed the deepest minima for the stacked states of the purine-purine dimers, but good stacking was also observed for the purine-pyrimidine and pyrimidine-purine dimers. Substantial stacking ability was found for the dimers composed of a thymine base and a purine base and also for the deoxythymidylyl-3',5'-deoxythymidine dimer. Very poor stacking was observed for the dCpdC dimer. Conformational properties and solvent accessibility are discussed for the stacked and unstacked dimers. The potential of mean force profiles of the stacking-unstacking process for the DNA dimers are compared with the RNA dimers.

Analysis of local helix bending in crystal structures of DNA oligonucleotides and DNA-protein complexes

Sequence-dependent bending of the helical axes in 112 oligonucleotide duplex crystal structures resident in the Nucleic Acid Database have been analyzed and compared with the use of bending dials, a computer graphics tool. Our analysis includes structures of both A and B forms of DNA and considers both uncomplexed forms of the double helix as well as those bound to drugs and proteins. The patterns in bending preferences in the crystal structures are analyzed by base pair steps, and emerging trends are noted. Analysis of the 66 B-form structures in the Nucleic Acid Database indicates that uniform trends within all pyrimidine-purine and purine-pyrimidine steps are not necessarily observed but are found particularly at CG and GC steps of dodecamers. The results support the idea that AA steps are relatively straight and that larger roll bends occur at or near the junctions of these A-tracts with their flanking sequences. The data on 16 available crystal structures of protein-DNA complexes indicate that the majority of the DNA bends induced via protein binding are sharp localized kinks. The analysis of the 30 available A-form DNA structures indicates that these structures are also bent and show a definitive preference for bending into the deep major groove over the shallow minor groove.

Structure and dynamics of a DNA.RNA hybrid duplex with a chiral phosphorothioate moiety: NMR and molecular dynamics with conventional and time-averaged restraints

The three-dimensional structure of two thiophosphate-modified DNA.RNA hybrid duplexes d(GCTATAApsTGG).r(CCAUUAUAGC), one with R-thiophosphate chirality and one with S-thiophosphate chirality, have been determined by restrained molecular dynamics simulations (rMD). As the two yielded almost identical results, a description of results can be presented in the singular. The conformational flexibility of this hybrid has been investigated by employing time-averaged constraints during the molecular dynamics simulations (MD-tar). A set of structural restraints, comprising 322 precise interproton distance constraints obtained by a complete relaxation matrix analysis of the 2D NOE intensities as well as J coupling constants obtained from quantitative simulations of DQF-COSY cross-peaks in deoxyriboses, was reported in our previous paper [González, C., Stec, W., Kobylanska, A., Hogrefe, R. I., Reynolds, M., & James, T. L. (1994) Biochemistry 33, 11062-11072]. Multiple conformations of the deoxyribose moieties were evident from the scalar coupling constant analysis. Accurate distance constraints, obtained from complete relaxation matrix analysis, yielded a time-averaged solution structure via conventional restrained molecular dynamics which is not compatible with the experimental J coupling constants (root-mean-square deviation in J value approximately 2 Hz). However, vicinal coupling constant information can be reproduced when time-averaged constraints are used during the molecular dynamics calculations instead of the conventional restraints (Jrms approximately 0.6 Hz). MD-tar simulations also improve the NMR R factors. This improvement is more evident in the DNA than in the RNA strand, where no indication of conformational flexibility had been obtained. Analysis of the MD-tar trajectories confirms that deoxyriboses undergo pucker transitions between the S and N domain, with the major conformer in the S domain. The ribose moieties in the RNA strand, however, remain in the N domain during the entire simulation. Conformations of deoxyriboses in the intermediate domain near O4'-endo are obtained when the average structure is calculated with conventional NMR restraints. Since these conformations cannot account for the experimental J coupling information, and they only appear in a very low population in the MD-tar ensemble, we conclude that intermediate E sugar puckers are artifacts produced by the attempt to fit all the structural constraints simultaneously when in reality more than one conformer is present. Most structural features of the duplex remain the same in the average structure and in the MD-tar ensemble, e.g., the minor groove width, exhibiting an intermediate value compared with those of canonical A- and B-like structures.(ABSTRACT TRUNCATED AT 250 WORDS)

Modifying the helical structure of DNA by design: recruitment of an architecture-specific protein to an enforced DNA bend

BACKGROUND

Proteins can force DNA to adopt distorted helical structures that are rarely if ever observed in naked DNA. The ability to synthesize DNA that contains defined helical aberrations would offer a new avenue for exploring the structural and energetic plasticity of DNA. Here we report a strategy for the enforcement of non-canonical helical structures through disulfide cross-linking; this approach is exemplified by the design and synthesis of an oligonucleotide containing a pronounced bend.

RESULTS

A localized bend was site-specifically introduced into DNA by the formation of a disulfide cross-link between the 5' adenines of a 5'-AATT-3' region in complementary strands of DNA. The DNA bend was characterized by high-resolution NMR structure determination of a cross-linked dodecamer and electrophoretic mobility assays on phased multimers, which together indicate that the cross-linked tetranucleotide induces a helical bend of approximately 30 degrees and a modest degree of unwinding. The enforced bend was found to stimulate dramatically the binding of an architecture-specific protein, HMG-D, to the DNA. DNase I foot-printing analysis revealed that the protein is recruited to the section of DNA that is bent.

CONCLUSIONS

The present study reports a novel approach for the investigation of non-canonical DNA structures and their recognition by architecture-specific proteins. The mode of DNA bending induced by disulfide cross-linking resembles that observed in structures of protein-DNA complexes. The results reveal common elements in the DNA-binding mode employed by sequence-specific and architecture-specific HMG proteins.

Crystal structure of a paired domain-DNA complex at 2.5 A resolution reveals structural basis for Pax developmental mutations

The 2.5 A resolution structure of a cocrystal containing the paired domain from the Drosophila paired (prd) protein and a 15 bp site shows structurally independent N-terminal and C-terminal subdomains. Each of these domains contains a helical region resembling the homeodomain and the Hin recombinase. The N-terminal domain makes extensive DNA contacts, using a novel beta turn motif that binds in the minor groove and a helix-turn-helix unit with a docking arrangement surprisingly similar to that of the lambda repressor. The C-terminal domain is not essential for prd binding and does not contact the optimized site. All known developmental missense mutations in the paired box of mammalian Pax genes map to the N-terminal subdomain, and most of them are found at the protein-DNA interface.

Molecular dynamics simulations of the glucocorticoid receptor DNA-binding domain in complex with DNA and free in solution

Molecular dynamics simulations have been performed on the glucocorticoid receptor DNA binding domain (GR DBD) in aqueous solution as a dimer in complex with DNA and as a free monomer. In the simulated complex, we find a slightly increased bending of the DNA helix axis compared with the crystal structure in the spacer region of DNA between the two half-sites that are recognized by GR DBD. The bend is mainly caused by an increased number of interactions between DNA and the N-terminal extended region of the sequence specifically bound monomer. The recognition helices of GR DBD are pulled further into the DNA major groove leading to a weakening of the intrahelical hydrogen bonds in the middle of the helices. Many ordered water molecules with long residence times are found at the intermolecular interfaces of the complex. The hydrogen-bonding networks (including water bridges) on either side of the DNA major groove involve residues that are highly conserved within the family of nuclear receptors. Very similar hydrogen-bonding networks are found in the estrogen receptor (ER) DBD in complex with DNA, which suggests that this is a common feature for proper positioning of the recognition helix in ER DBD and GR DBD.

Molecular dynamics simulations of chlorambucil/DNA adducts. A structural basis for the 5'-GNC interstrand DNA crosslink formed by nitrogen mustards

The alkylation of DNA by chlorambucil has been studied using a computational approach. Molecular dynamics simulations were performed on the fully solvated non-covalent complex, two monoadducts and a crosslinked diadduct of chlorambucil with the d(CGG3G2CGC).-d(GCG1CCCG) duplex, in which the N7 atoms of G1, G2 and G3 are potential alkylation sites. The results provide a structural basis for the preference of nitrogen mustards to crosslink DNA duplexes at a 5'-GNC site (a 1,3 crosslink, G1-G3) rather than at a 5'-GC sites (a 1,2 crosslink, G1-G2). In the non-covalent complex simulation the drug reoriented from a non-interstrand crosslinking location to a position favorable for G1-G3 diadduct formation. It proved possible to construct a G1-G3 diadduct from a structure from the non-covalent simulation, and continue the molecular dynamics calculation without further disruption of the DNA structure. A crosslinked diadduct developed with four BII conformations on the 3' side of each alkylated guanine and of their respective complementary cytosine. In the first monoadduct simulation the starting point was the same DNA conformation used in the crosslinked diadduct simulation with alkylation at G1. In this simulation the DNA deformation was reduced, with the helix returning to a more canonical form. A second monoadduct simulation was started from a canonical DNA conformation alkylated at G3. Here, no significant motion towards a potential crosslinking conformation occurred. Collectively, the results suggest that crosslink formation is dependent upon the drug orientation prior to alkylation and the required deformation of the DNA to permit 1,3 crosslinking can largely be achieved in the non-covalent complex.

Molecular dynamics simulations of the effects of ring-saturated thymine lesions on DNA structure

The effect of thymine lesions produced by radiation or oxidative damage on DNA structure was studied by molecular dynamics simulations of native and damaged DNA. Thymine in position 7 of native dodecamer d(CGCGAATTCGCG)2 was replaced by one of the four thymine lesions 5-hydroxy-5,6-dihydrothymine, 6-hydroxy-5,6-dihydrothymine (thymine photohydrate), 5,6-dihydroxy-5,6-dihydro-thymine (thymine glycol), and 5,6-dihydrothymine. Simulations were performed with Assisted Model Building with Energy Refinement force field. Solvent was represented by a rectangular box of water with periodic boundary conditions applied. A constant temperature and constant volume protocol was used. The observed level of distortions of DNA structure depends on the specific nature of the lesion. The 5,6-dihydrothymine does not cause distinguishable perturbations to DNA. Other lesions produce a dramatic increase in the rise parameter between the lesion and the 5' adjacent adenine. These changes are accompanied by weakening of Watson-Crick hydrogen bonds in the A6-T19 base pair on the 5' side of the lesion. The lesioned bases also show negative values of inclination relative to the helical axis. No changes in the pattern of backbone torsional angles are observed with any of the lesions incorporated into DNA. The structural distortions in DNA correlate well with known biological effects of 5,6-dihydrothymine and thymine glycol on such processes as polymerase action or recognition by repair enzymes.

Refinement of the solution structure of a branched DNA three-way junction

We have refined the structure of the DNA Three-Way Junction complex, TWJ-TC, described in the companion paper by quantitative analysis of two 2D NOESY spectra (mixing times 60 and 200 ms) obtained in D2O solution. NOESY crosspeak intensities extracted from these spectra were used in two kinds of refinement procedure: 1) distance-restrained energy minimization (EM) and molecular dynamics (MD) and 2) full relaxation matrix back calculation refinement. The global geometry of the refined model is very similar to that of a published, preliminary model (Leontis, 1993). Two of the helical arms of the junction are stacked. These are Helix 1, defined by basepairs S1-G1/S3-C12 through S1-C5/S3-G8 and Helix 2, which comprises basepairs S1-C6/S2-G5 through S1-G10/S2-G1. The third helical arm (Helix 3), comprised of basepairs S2-C6/S3-G5 through S2-C10/S3-G1 extends almost perpendicularly from the axis defined by Helices 1 and 2. The bases S1-C5 and S1-C6 of Strand 1 are continuously stacked across the junction region. The conformation of this strand is close to that of B-form DNA along its entire length, including the S1-C5 to S1-C6 dinucleotide step at the junction. The two unpaired bases S3-T6 and S3-C7 lie outside of the junction along the minor groove of Helix 1 and largely exposed to solvent. Analysis of the refined structure reveals that the glycosidic bond of S3-T6 exists in the syn conformation, allowing the methyl group of this residue to contact the hydrophobic surface of the minor groove of Helix 1, at S3-G11. The helical parameters of the three helical arms of the structure exhibit only weak deviations from typical values for right-handed B-form DNA. Unusual dihedral angles are only observed for the sugarphosphate backbone joining the "hinge" residues, S2-G5 and S2-C6, and S3-G5 through S3-G8. The glycosidic bond of S3-G8also lies within the syn range, allowing favorable Watson-Crick base-pairing interactions with Si -C5. The stability of this structure was checked in 39 ps molecular dynamic simulation at 330 K in water. The structure of TWJ-TC retained the geometrical features mentioned above at the end of the simulation period. The final R(1/6)-factor of the refined structure is 5%.

Structure determination and analysis of local bending in an A-tract DNA duplex: comparison of results from crystallography, nuclear magnetic resonance, and molecular dynamics simulation on d(CGCAAAAATGCG)

We have presented a detailed analysis for structure determinations for the DNA duplex d(CGCAAAAATGCG) obtained from X-ray crystallography, nuclear magnetic resonance, and molecular dynamics simulation. Each of the structures for the duplex deviates from the structure of the canonical form of B-DNA in a number of observable characteristics. Specifically, the three determinations all contain DNA axis deflections at the junctions of the A-tract with the flanking sequences. The analysis provided shows that the general characteristics of the structures obtained for d(CGCAAAAATGCG) from X-ray, NMR, and MD methods turn out to be quite similar. The extent to which this result can be generalized remains to be established by consideration of similar cross-comparisons on diverse oligonucleotide structures.

Molecular dynamics simulations suggest that the Eco RI kink is an example of molecular strain

The energy surface in the vicinity of the "Eco RI kink" was investigated by conducting both in vacuo molecular dynamics simulations as well as a simulation with explicit solvent. The in vacuo simulations used the "all atom" AMBER 3.0 force field with a distant dependent dielectric function and "hydrated" counter ions while the simulation with explicit solvent used the AMBER 4.0 force field, fully charged phosphates and counter ions and a dielectric constant of 1.0. The thrust of the simulations was to discriminate between two models of the energy surface of the deformed DNA as found in the recognition complex with Eco RI endonuclease. In the intrinsic model, the kinked DNA is a local minimum of the energy surface intrinsic to the DNA itself while in the strained model there is no significant energy barrier separating kinked and regular B-DNA. The two models have significant implications for theories of indirect recognition of DNA based on sequence-dependent deformability. The simulations suggest that the Eco RI-kinked structure is an example of molecular strain because it is not near a minimum of any of the potential energy functions examined. The simulations leave the question of an energy barrier somewhat open and raise the possibility that the Eco RI kink is at (or near) a point of dynamic instability of the energy surface (either a true maximum or a saddle point).