Thomas C. Bishop

A systematic molecular dynamics study of nearest-neighbor effects on base pair and base pair step conformations and fluctuations in B-DNA

It is well recognized that base sequence exerts a significant influence on the properties of DNA and plays a significant role in protein–DNA interactions vital for cellular processes. Understanding and predicting base sequence effects requires an extensive structural and dynamic dataset which is currently unavailable from experiment. A consortium of laboratories was consequently formed to obtain this information using molecular simulations. This article describes results providing information not only on all 10 unique base pair steps, but also on all possible nearest-neighbor effects on these steps. These results are derived from simulations of 50–100 ns on 39 different DNA oligomers in explicit solvent and using a physiological salt concentration. We demonstrate that the simulations are converged in terms of helical and backbone parameters. The results show that nearest-neighbor effects on base pair steps are very significant, implying that dinucleotide models are insufficient for predicting sequence-dependent behavior. Flanking base sequences can notably lead to base pair step parameters in dynamic equilibrium between two conformational sub-states. Although this study only provides limited data on next-nearest-neighbor effects, we suggest that such effects should be analyzed before attempting to predict the sequence-dependent behavior of DNA.

Difficulties with multiple time stepping and fast multipole algorithm in molecular dynamics

Numerical experiments are performed on a 36,000‐atom protein–DNA–water simulation to ascertain the effectiveness of two devices for reducing the time spent computing long‐range electrostatics interactions. It is shown for Verlet‐I/r‐RESPA multiple time stepping, which is based on approximating long‐range forces as widely separated impulses, that a long time step of 5 fs results in a dramatic energy drift and that this is reduced by using an even larger long time step. It is also shown that the use of as many as six terms in a fast multipole algorithm approximation to long‐range electrostatics still fails to prevent significant energy drift even though four digits of accuracy is obtained. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1785–1791, 1997

μABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA.

We present the results of microsecond molecular dynamics simulations carried out by the ABC group of laboratories on a set of B-DNA oligomers containing the 136 distinct tetranucleotide base sequences. We demonstrate that the resulting trajectories have extensively sampled the conformational space accessible to B-DNA at room temperature. We confirm that base sequence effects depend strongly not only on the specific base pair step, but also on the specific base pairs that flank each step. Beyond sequence effects on average helical parameters and conformational fluctuations, we also identify tetranucleotide sequences that oscillate between several distinct conformational substates. By analyzing the conformation of the phosphodiester backbones, it is possible to understand for which sequences these substates will arise, and what impact they will have on specific helical parameters.

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.

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 nucleosomes 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.

Atomistic simulations of nucleosomes

Nearly a dozen all‐atom molecular dynamics (MD) simulations of the nucleosome have been performed. Collectively, these simulations provide insights into the structure and dynamics of the biomolecular complex that serves as the fundamental folding unit of chromatin. Nucleosomes contain 146 base pairs of DNA wrapped in a left‐handed superhelix around a core of eight histones. This review provides a survey of what has been learned about DNA, histones, and solvent interactions based on all‐atom MD studies of the nucleosome. The longest simulations to date are on the order of 100 nanoseconds. On this time scale, nucleosomes are quite stable. DNA kinks, the histone tails, solvent, and ions are highly dynamic and can be readily investigated using equilibrium dynamics methods. Steered MD is required to observe large‐scale structural changes. The need for explicit solvent techniques is underscored by the inability of continuum solvent methods to properly describe the ion‐nucleosome radial distribution functions. The atomistic techniques reviewed here are deemed necessary for exploration of the near infinite variations in atomic composition that exists even in the canonical nucleosome octamer. Continued development of these nascent simulation efforts will enable experimentalists to utilize rational design strategies in their efforts to investigate nucleosomes and chromatin.

Geometry of the Nucleosomal DNA Superhelix

Nucleosome stability is largely an indirect measure of DNA sequence based on the material properties of DNA and the ability of a sequence to assume the required left-handed superhelical conformation. Here we focus attention only on the geometry of the superhelix and present two distinct mathematical expressions that rely on the DNA helical parameters (Shift, Slide, Rise, Tilt, Roll, Twist). One representation requires torsion for superhelix formation; the other requires shear. To compare these mathematical expressions to experimental data we develop a strategy for Fourier-filtering the helical parameters that identifies necessary and sufficient conditions to achieve a high-resolution model of the nucleosome superhelix. We apply this filtering strategy to 24 high-resolution structures of the nucleosome and demonstrate that all structures have a highly conserved distribution of Roll, Slide and Twist that involves two length scales. One length scale spans the entire length of nucleosomal DNA. The other is associated with the helix repeat. Our strategy also enables us to identify ground state or simple nucleosomes and altered nucleosome structures. These results form a basis for characterizing structural variations in the emerging family of nucleosome structures and a method for further developing structure-based models of nucleosome stability.

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° 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.

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 nucleosomes 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.

Evaluation of Elastic Rod Models with Long Range Interactions for Predicting Nucleosome Stability

Abstract The ability of a dinucleotide-step based elastic-rod model of DNA to predict nucleosome binding free energies is investigated using four available sets of elastic parameters. We compare the predicted free energies to experimental values derived from nucleosome reconstitution experiments for 84 DNA sequences. Elastic parameters (conformation and stiffnessess) obtained from MD simulations are shown to be the most reliable predictors, as compared to those obtained from analysis of base-pair step melting temperatures, or from analysis of x-ray structures. We have also studied the effect of varying the folded conformation of nucleosomal DNA by means of our Fourier filtering knock-out and knock-in procedure. This study confirmed the above ranking of elastic parameters, and helped to reveal problems inherent in models using only a local elastic energy function. Long-range interactions were added to the elastic-rod model in an effort to improve its predictive ability. For this purpose a Debye-Huckel energy term with a single, homogenous point charge per base- pair was introduced. This term contains only three parameters,—its weight relative to the elastic energy, the Debye screening length, and a minimum sequence distance for including pairwise interactions between charges. After optimization of these parameters, our Debye-Huckel term is attractive, and yields the same level of correlation with experiment (R = 0.75) as was achieved merely by varying the nucleosomal shape in the elastic-rod model. We suggest this result indicates a linker DNA—histone attraction or, possibly, entropic effects, that lead to a stabilization of a nucleosome away from the ends of DNA segments longer than 147 bp. Such effects are not accounted for by a localized elastic energy model.