The effect of structural water molecules on the normal mode spectrum of dTn . dAn x dTn DNA

In this paper we present a theoretical treatment of triplex B type DNA hydration using normal mode calculation techniques. Discrete solvent is added as spines of hydration in the Watson-Crick and Crick-Hoogsteen grooves as well as water bridges between the Phosphate groups. The effect of binding the discrete structural waters on the normal mode of vibration of the system was studied by introducing a parameter, Xw, that is proportional to the degree of water binding and inversely proportional to the relative humidity (RH) of the system. We examined the variation of the dipole moments of characteristic modes with Xw. The results show that there is a direct relationship between the degree of binding of the water molecules to the atoms in the triple helix, the relative humidity of the system and the conformation and stability of the triple helix. At high RH and Xw = 0:0 the triple helix has mostly B type conformation characteristics, with C'2 -endo sugars. The emergence of normal modes of vibration characteristic to the A type conformation (C'3 - endo sugars) at Xw = 0:4 and 60% RH indicates a conformational shift towards A-type for some of the sugars between Xw = 0.2 (80% RH) and Xw = 0.4 (60% RH). These results are in agreement with the "economy of hydration hypothesis" of Saenger (Saenger et al., 1986) which maintains that the main difference in the hydration of A- and B- forms of DNA is the presence of water bridges between adjacent Phosphate groups in the low-hydration A-form but not in the B- form. Free energy calculations for the triplex DNA with structural waters show that there is a minimum of the free energy at Xw = 0.2 and the free energy increases with Xw and becomes larger than the free energy of the B conformation without structural waters for Xw equal to and larger than 0.4. This result indicates that the B conformation is more stable with bound structural water molecules (for degrees of water binding that are not over 20% higher than the degree of binding between bulk water molecules). The structural water molecules are bound much tighter in the A conformation than in the B conformation. The model predicts that the B to A transition occurs at higher relative humidities in D2O than in H2O. Part of these results (Dadarlat, 1997) have been subsequently confirmed by the experimental work and MD simulations of Ouali (Ouali et al., 1997). The experimental results showed that the N-type sugars corresponding to the A conformation are clearly detected below 75% RH.     

One left-handed strand in DNA-oligonucleotide complexes?

A single strand of oligonucleotide can bind to double helical DNA under certain conditions. This must involve some unwinding of the original double helix in a process leading to the formation of a three-stranded region. The free energy for such an entropically unlikely reaction may come from a change in the degree of supercoiling of the original DNA. The conformation of the triple strand is investigated here using computer graphics and molecular mechanics calculations. It is suggested that on binding the oligonucleotide (strand 3) to two paired strands (1 and 2) in a supercoiled DNA molecule, strand 2 might adopt a left-handed conformation whilst strand 1 and strand 3 pair in the normal Watson-Crick B-configuration.     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

Effective energy function for proteins in solution

A Gaussian solvent-exclusion model for the solvation free energy is developed. It is based on theoretical considerations and parametrized with experimental data. When combined with the CHARMM 19 polar hydrogen energy function, it provides an effective energy function (EEF1) for proteins in solution. The solvation model assumes that the solvation free energy of a protein molecule is a sum of group contributions, which are determined from values for small model compounds. For charged groups, the self-energy contribution is accounted for primarily by the exclusion model. Ionic side-chains are neutralized, and a distance-dependent dielectric constant is used to approximate the charge-charge interactions in solution. The resulting EEF1 is subjected to a number of tests. Molecular dynamics simulations at room temperature of several proteins in their native conformation are performed, and stable trajectories are obtained. The deviations from the experimental structures are similar to those observed in explicit water simulations. The calculated enthalpy of unfolding of a polyalanine helix is found to be in good agreement with experimental data. Results reported elsewhere show that EEF1 clearly distinguishes correctly from incorrectly folded proteins, both in static energy evaluations and in molecular dynamics simulations and that unfolding pathways obtained by high-temperature molecular dynamics simulations agree with those obtained by explicit water simulations. Thus, this energy function appears to provide a realistic first approximation to the effective energy hypersurface of proteins.     

Acoustic modes and nonbonded interactions of the double helix

Observations of acoustic velocities in DNA fibers have been used to refine nonbonded force constants for the DNA double helix. Long-range forces are found to be needed for A conformation and are likely to dominate in B conformation as well. The acoustic dispersion curves are described and calculated. A correction due to the effects of water is calculated. The effect of nonbonded interaction on other vibrational modes is calculated.     

Folding of a DNA hairpin loop structure in explicit solvent using replica-exchange molecular dynamics simulations

Hairpin loop structures are common motifs in folded nucleic acids. The 5'-GCGCAGC sequence in DNA forms a characteristic and stable trinucleotide hairpin loop flanked by a two basepair stem helix. To better understand the structure formation of this hairpin loop motif in atomic detail, we employed replica-exchange molecular dynamics (RexMD) simulations starting from a single-stranded DNA conformation. In two independent 36 ns RexMD simulations, conformations in very close agreement with the experimental hairpin structure were sampled as dominant conformations (lowest free energy state) during the final phase of the RexMDs ( approximately 35% at the lowest temperature replica). Simultaneous compaction and accumulation of folded structures were observed. Comparison of the GCA trinucleotides from early stages of the simulations with the folded topology indicated a variety of central loop conformations, but arrangements close to experiment that are sampled before the fully folded structure also appeared. Most of these intermediates included a stacking of the C(2) and G(3) bases, which was further stabilized by hydrogen bonding to the A(5) base and a strongly bound water molecule bridging the C(2) and A(5) in the DNA minor groove. The simulations suggest a folding mechanism where these intermediates can rapidly proceed toward the fully folded hairpin and emphasize the importance of loop and stem nucleotide interactions for hairpin folding. In one simulation, a loop motif with G(3) in syn conformation (dihedral flip at N-glycosidic bond) accumulated, resulting in a misfolded hairpin. Such conformations may correspond to long-lived trapped states that have been postulated to account for the folding kinetics of nucleic acid hairpins that are slower than expected for a semiflexible polymer of the same size.     

Polyproline II helix is the preferred conformation for unfolded polyalanine in water

Does aqueous solvent discriminate among peptide conformers? To address this question, we computed the solvation free energy of a blocked, 12‐residue polyalanyl‐peptide in explicit water and analyzed its solvent structure. The peptide was modeled in each of 4 conformers: α‐helix, antiparallel β‐strand, parallel β‐strand, and polyproline II helix (P_II). Monte Carlo simulations in the canonical ensemble were performed at 300 K using the CHARMM 22 forcefield with TIP3P water. The simulations indicate that the solvation free energy of P_II is favored over that of other conformers for reasons that defy conventional explanation. Specifically, in these 4 conformers, an almost perfect correlation is found between a residue's solvent‐accessible surface area and the volume of its first solvent shell, but neither quantity is correlated with the observed differences in solvation free energy. Instead, solvation free energy tracks with the interaction energy between the peptide and its first‐shell water. An additional, previously unrecognized contribution involves the conformation‐dependent perturbation of first‐shell solvent organization. Unlike P_II, β‐strands induce formation of entropically disfavored peptide:water bridges that order vicinal water in a manner reminiscent of the hydrophobic effect. The use of explicit water allows us to capture and characterize these dynamic water bridges that form and dissolve during our simulations. Proteins 2004. © 2004 Wiley‐Liss, Inc.     

Neutron fibre diffraction study of DNA hydration

Interactions with water are crucial to the conformation assumed by the DNA double helix. The location of water around the D conformation has been investigated in a neutron fibre diffraction study which shows that water is ordered in the minor groove of the DNA. The D conformation is important since its occurrence is limited to specific DNA base pair sequences which have been identified as functionally significant. This study is of particular interest because the D conformation has not been reported in single crystal studies of oligonucleotides.     

Molecular dynamics simulations of a nucleosome and free DNA

All atom molecular dynamics simulations (10ns) of a nucleosome and of its 146 basepairs of DNA free in solution have been conducted. DNA helical parameters (Roll, Tilt, Twist, Shift, Slide, Rise) were extracted from each trajectory to compare the conformation, effective force constants, persistence length measures, and fluctuations of nucleosomal DNA to free DNA. The conformation of DNA in the nucleosome, as determined by helical parameters, is found to be largely within the range of thermally accessible values obtained for free DNA. DNA is found to be less flexible on the nucleosome than when free in solution, however such measures are length scale dependent. A method for disassembling and reconstructing the conformation and dynamics of the nucleosome using Fourier analysis is presented. Long length variations in the conformation of nucleosomal DNA are identified other than those associated with helix repeat. These variations are required to create a proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns of the nucleosome's superhelix. Reconstruction of free DNA using selected long wavelength variations in conformation can produce either a left-handed or a right-handed superhelix. The long wavelength variations suggest 146 basepairs is a natural length of DNA to wrap around the histone core.     

Prediction of conformation of rat galanin in the presence and absence of water with the use of monte Carlo methods and the ECEPP/3 force field

The conformation of the 29-residue rat galanin neuropeptide was studied using the Monte Carlo with energy minimization (MCM) and electrostatically driven Monte Carlo (EDMC) methods. According to a previously elaborated procedure, the polypeptide chain was first treated in a united-residue approximation, in order to enable extensive exploration of the conformational space to be carried out (with the use of MCM), Then the low-energy united-residue conformations were converted to the all-atom representations, and EDMC simulations were carried out for the all-atom polypeptide chains, using the ECEPP/3 force field with hydration included. In order to estimate the effect of environment on galanin conformation, the low-energy conformations obtained as a result of these simulations were taken as starting structures for further EDMC runs that did not include hydration. The lowest-energy conformation obtained in aqueous solution calculations had a nonhelical N-terminal part packed against the nonpolar face of a residual helix that extended from Pro¹³ toward the C-terminus. One next lowest-energy structure was a nearly-all-helical conformation, but with a markedly higher energy. In contrast, all of the low-energy conformations in the absence of water were all-helical differing only by the extent to which the helix was kinked around Pro¹³. These results are in qualitative agreement with the available NMR and CD data of galanin in aqueous and nonaqueous solvents.     

Computer simulations of a tumor surface octapeptide epitope

Molecular dynamics using CHARMM and GEMM programs with the Star Technologies ST 100 array processor functioning at the speed of super computers was used as a searching algorithm for conformational exploration of the octapeptide Gly-Asn-Thr-Ile-Val-Ala-Glu. This poorly soluble octapeptide is the N-terminal epitope of an 11 KD glycoprotein antigen residing on human ductal carcinoma (breast) cells. Very long (nanoseconds) simulations were required. Both an alpha-helix and the N-acetyl-N1-methylamide derived minimized starting structures gave the same lowest potential energy conformation with simulations at 600 K. The same conformation was found only when using the latter starting conformation with simulations at 300 K. The lowest potential energy conformation was stabilized by 4 hydrophobic contacts and 13 H bonds completing one turn of a left-handed helix.     

Phosphorus-31 nuclear magnetic resonance of double- and triple-helical nucleic acids. Phosphorus-31 chemical shifts as a probe of phosphorus-oxygen ester bond torsional angles

The temperature dependence to the 31P NMR spectra of poly[d(GC)] . poly [d(GC)],d(GC)4, phenylalanine tRNA (yeast) and mixtures of poly(A) + oligo(U) is presented. The 31P NMR spectra of mixtures of complementary RNA and of the poly d(GC) self-complementary DNA provide torsional information on the phosphate ester conformation in the double, triple, and "Z" helix. The increasing downfield shift with temperature of the single-strand nucleic acids provides a measure of the change in the phosphate ester conformation in the single helix to coil conversion. A separate upfield peak (20-60% of the total phosphates) is observed at lower temperatures in the oligo(U) . poly(A) mixtures which is assigned to the double helix/triple helix. Proton NMR and UV spectra confirm the presence of the multistrand forms. The 31P chemical shift for the double helix/triple helix is 0.2-0.5 ppm upfield from the chemical shift for the single helix which in turn is 1.0 ppm upfield from the chemical shift for the random coil conformation.     

Mechanical stability of single DNA molecules

Using a modified atomic force microscope (AFM), individual double-stranded (ds) DNA molecules attached to an AFM tip and a gold surface were overstretched, and the mechanical stability of the DNA double helix was investigated. In lambda-phage DNA the previously reported B-S transition at 65 piconewtons (pN) is followed by a second conformational transition, during which the DNA double helix melts into two single strands. Unlike the B-S transition, the melting transition exhibits a pronounced force-loading-rate dependence and a marked hysteresis, characteristic of a nonequilibrium conformational transition. The kinetics of force-induced melting of the double helix, its reannealing kinetics, as well as the influence of ionic strength, temperature, and DNA sequence on the mechanical stability of the double helix were investigated. As expected, the DNA double helix is considerably destabilized under low salt buffer conditions (相似文献   

An effective mesoscopic model of double-stranded DNA

Watson and Crick’s epochal presentation of the double helix structure in 1953 has paved the way to intense exploration of DNA’s vital functions in cells. Also, recent advances of single molecule techniques have made it possible to probe structures and mechanics of constrained DNA at length scales ranging from nanometers to microns. There have been a number of atomistic scale quantum chemical calculations or molecular level simulations, but they are too computationally demanding or analytically unfeasible to describe the DNA conformation and mechanics at mesoscopic levels. At micron scales, on the other hand, the wormlike chain model has been very instrumental in describing analytically the DNA mechanics but lacks certain molecular details that are essential in describing the hybridization, nano-scale confinement, and local denaturation. To fill this fundamental gap, we present a workable and predictive mesoscopic model of double-stranded DNA where the nucleotides beads constitute the basic degrees of freedom. With the inter-strand stacking given by an interaction between diagonally opposed monomers, the model explains with analytical simplicity the helix formation and produces a generalized wormlike chain model with the concomitant large bending modulus given in terms of the helical structure and stiffness. It also explains how the helical conformation undergoes overstretch transition to the ladder-like conformation at a force plateau, in agreement with the experiment.     

Water-DNA interactions as studied by X-ray and neutron fibre diffraction

X-ray fibre-diffraction studies indicate a high degree of stereochemical specificity in interactions between water and the DNA double helix. Evidence for this comes from data that show that the molecular conformations assumed by DNA in fibres are highly reproducible and that the hydration-driven transitions between these conformations are fully reversible. These conformational transitions are induced by varying the relative humidity of the fibre environment and hence its water content. Further evidence for stereochemical specificity comes from the observed dependence of the conformation assumed on the ionic content of the fibre and the nucleotide sequence of the DNA. For some transitions, information on stereochemical pathways has come from real-time X-ray fibre diffraction using synchrotron radiation; information on the location of water with respect to the double helix for a number of DNA conformations has come from neutron fibre diffraction. This structural information from fibre-diffraction studies of DNA is complemented by information from X-ray single-crystal studies of oligonucleotides. If the biochemical processes involving DNA have evolved to exploit the structural features observed in DNA fibres and oligonucleotide single crystals, the challenges in developing alternatives to a water environment can be expected to be very severe.     

Molecular Dynamics Simulations of a Nucleosome and Free DNA

Abstract

All atom molecular dynamics simulations (10ns) of a nucleosome and of its 146 basepairs of DNA free in solution have been conducted. DNA helical parameters (Roll, Tilt, Twist, Shift, Slide, Rise) were extracted from each trajectory to compare the conformation, effective force constants, persistence length measures, and fluctuations of nucleosomal DNA to free DNA. The conformation of DNA in the nucleosome, as determined by helical parameters, is found to be largely within the range of thermally accessible values obtained for free DNA. DNA is found to be less flexible on the nucleosome than when free in solution, however such measures are length scale dependent. A method for disassembling and reconstructing the conformation and dynamics of the nucleosome using Fourier analysis is presented. Long length variations in the conformation of nucleosomal DNA are identified other than those associated with helix repeat. These variations are required to create a proposed tetrasome conformation or to qualitatively reconstruct the 1.75 turns of the nucleosome's superhelix. Reconstruction of free DNA using selected long wavelength variations in conformation can produce either a left-handed or a right-handed superhelix. The long wavelength variations suggest 146 basepairs is a natural length of DNA to wrap around the histone core.     

Conformation of endothelin in aqueous ethylene glycol determined by 1H-NMR and molecular dynamics simulations

The solution conformation of a 21-residue vasoconstrictor peptide endothelin-1 (ET-1) in water-ethylene glycol has been determined by two-dimensional 1H-NMR spectroscopy and constrained molecular dynamics simulations. The N-terminus (residues 1-4) appears to undergo conformational averaging and no single structure consistent with the NMR constraints could be found for this region. Residues 5-8 form a turn, and residues 9-16 exist in a helical conformation. A flexible 'hinge' between residues 8-9 allows various orientations of the turn relative to the helix. Another 'hinge' at residue 17 connects the extended C-terminus to the bicyclic core region (residues 1-15). Residues important for binding and biological activity form a contiguous surface on one side of the helix, with the two disulfides extending from the other side of the helix.     

Conformational flexibility of B-DNA at 0.74 A resolution: d(CCAGTACTGG)(2)

The affinity and specificity of a ligand for its DNA site is a function of the conformational changes between the isolated and complexed states. Although the structures of a hydroxypyrrole-imidazole-pyrrole polyamide dimer with 5'-CCAGTACTGG-3' and the trp repressor recognizing the sequence 5'-GTACT-3' are known, the baseline conformation of the DNA site would contribute to our understanding of DNA recognition by these ligands. The 0.74 A resolution structure of a B-DNA double helix, 5'-CCAGTACTGG-3', has been determined by X-ray crystallography. Six of the nine phosphates, two of four bound calcium ions and networks of water molecules hydrating the oligonucleotide have alternate conformations. By contrast, nine of the ten bases have a single, unique conformation with hydrogen atoms visible in most cases. The polyamide molecules alter the geometry of the phosphodiester backbone, and the water molecules mediating contacts in the trp repressor/operator complex are conserved in the unliganded DNA. Furthermore, the multiple conformational states, ions and hydration revealed by this ultrahigh resolution structure of a B-form oligonucleotide are potentially general considerations for understanding DNA-binding affinity and specificity by ligands.     

Conformation of the Flavin Adenine Dinucleotide Cofactor FAD in DNA-Photolyase: A Molecular Dynamics Study

In order to gain insight into the light-driven repair of DNA by the enzyme DNA photolyase, the conformation of the photoactive cofactor FAD, a flavin adenine dinucleotide, has been studied by molecular dynamic simulations. In contrast to FAD in the gas phase and in water where the MD procedure yields various "open" I-shaped as well as "closed" U-shaped conformations, the calculations of FAD binding to the enzyme show essentially a single U-shaped conformation of this cofactor which, so far, is unique among FAD-carrying proteins. It is characteristic for this U-shaped conformation that the FAD components occupy opposite sides of the pocket in the surface of the protein which provides the binding site for the defect pyrimidine dimer structure on DNA. In fact, the calculated U-shaped conformation is very close to the one revealed by the X-ray structure analysis of DNA photolyase. Moreover, the simulations yield details on the binding of the photoactive isoalloxazine moiety and the dynamics of the amino acids forming the binding cavity of the enzyme.