Steven C. Howell

DNA Methylation Regulated Nucleosome Dynamics

A strong correlation between nucleosome positioning and DNA methylation patterns has been reported in literature. However, the mechanistic model accounting for the correlation remains elusive. In this study, we evaluated the effects of specific DNA methylation patterns on modulating nucleosome conformation and stability using FRET and SAXS. CpG dinucleotide repeats at 10 bp intervals were found to play different roles in nucleosome stability dependent on their methylation states and their relative nucleosomal locations. An additional (CpG)5 stretch located in the nucleosomal central dyad does not alter the nucleosome conformation, but significant conformational differences were observed between the unmethylated and methylated nucleosomes. These findings suggest that the correlation between nucleosome positioning and DNA methylation patterns can arise from the variations in nucleosome stability dependent on their sequence and epigenetic content. This knowledge will help to reveal the detailed role of DNA methylation in regulating chromatin packaging and gene transcription.

Ion Competition in Condensed DNA Arrays in the Attractive Regime

Physical origin of DNA condensation by multivalent cations remains unsettled. Here, we report quantitative studies of how one DNA-condensing ion (Cobalt(3+) Hexammine, or Co(3+)Hex) and one nonDNA-condensing ion (Mg(2+)) compete within the interstitial space in spontaneously condensed DNA arrays. As the ion concentrations in the bath solution are systematically varied, the ion contents and DNA-DNA spacings of the DNA arrays are determined by atomic emission spectroscopy and x-ray diffraction, respectively. To gain quantitative insights, we first compare the experimentally determined ion contents with predictions from exact numerical calculations based on nonlinear Poisson-Boltzmann equations. Such calculations are shown to significantly underestimate the number of Co(3+)Hex ions, consistent with the deficiencies of nonlinear Poisson-Boltzmann approaches in describing multivalent cations. Upon increasing the concentration of Mg(2+), the Co(3+)Hex-condensed DNA array expands and eventually redissolves as a result of ion competition weakening DNA-DNA attraction. Although the DNA-DNA spacing depends on both Mg(2+) and Co(3+)Hex concentrations in the bath solution, it is observed that the spacing is largely determined by a single parameter of the DNA array, the fraction of DNA charges neutralized by Co(3+)Hex. It is also observed that only ∼20% DNA charge neutralization by Co(3+)Hex is necessary for spontaneous DNA condensation. We then show that the bath ion conditions can be reduced to one variable with a simplistic ion binding model, which is able to describe the variations of both ion contents and DNA-DNA spacings reasonably well. Finally, we discuss the implications on the nature of interstitial ions and cation-mediated DNA-DNA interactions.

Solution Scattering and FRET Studies on Nucleosomes Reveal DNA Unwrapping Effects of H3 and H4 Tail Removal

Using a combination of small-angle X-ray scattering (SAXS) and fluorescence resonance energy transfer (FRET) measurements we have determined the role of the H3 and H4 histone tails, independently, in stabilizing the nucleosome DNA terminal ends from unwrapping from the nucleosome core. We have performed solution scattering experiments on recombinant wild-type, H3 and H4 tail-removed mutants and fit all scattering data with predictions from PDB models and compared these experiments to complementary DNA-end FRET experiments. Based on these combined SAXS and FRET studies, we find that while all nucleosomes exhibited DNA unwrapping, the extent of this unwrapping is increased for nucleosomes with the H3 tails removed but, surprisingly, decreased in nucleosomes with the H4 tails removed. Studies of salt concentration effects show a minimum amount of DNA unwrapping for all complexes around 50-100mM of monovalent ions. These data exhibit opposite roles for the positively-charged nucleosome tails, with the ability to decrease access (in the case of the H3 histone) or increase access (in the case of the H4 histone) to the DNA surrounding the nucleosome. In the range of salt concentrations studied (0-200mM KCl), the data point to the H4 tail-removed mutant at physiological (50-100mM) monovalent salt concentration as the mononucleosome with the least amount of DNA unwrapping.

Elucidating Internucleosome Interactions and the Roles of Histone Tails

The nucleosome is the first level of genome organization and regulation in eukaryotes where negatively charged DNA is wrapped around largely positively charged histone proteins. Interaction between nucleosomes is dominated by electrostatics at long range and guided by specific contacts at short range, particularly involving their flexible histone tails. We have thus quantified how internucleosome interactions are modulated by salts (KCl, MgCl2) and histone tail deletions (H3, H4 N-terminal), using small-angle x-ray scattering and theoretical modeling. We found that measured effective charges at low salts are ∼1/5th of the theoretically predicted renormalized charges and that H4 tail deletion suppresses the attraction at high salts to a larger extent than H3 tail deletion.

Monte Carlo simulation algorithm for B-DNA

Understanding the structure–function relationship of biomolecules containing DNA has motivated experiments aimed at determining molecular structure using methods such as small‐angle X‐ray and neutron scattering (SAXS and SANS). SAXS and SANS are useful for determining macromolecular shape in solution, a process which benefits by using atomistic models that reproduce the scattering data. The variety of algorithms available for creating and modifying model DNA structures lack the ability to rapidly modify all‐atom models to generate structure ensembles. This article describes a Monte Carlo algorithm for simulating DNA, not with the goal of predicting an equilibrium structure, but rather to generate an ensemble of plausible structures which can be filtered using experimental results to identify a sub‐ensemble of conformations that reproduce the solution scattering of DNA macromolecules. The algorithm generates an ensemble of atomic structures through an iterative cycle in which B‐DNA is represented using a wormlike bead–rod model, new configurations are generated by sampling bend and twist moves, then atomic detail is recovered by back mapping from the final coarse‐grained configuration. Using this algorithm on commodity computing hardware, one can rapidly generate an ensemble of atomic level models, each model representing a physically realistic configuration that could be further studied using molecular dynamics.

Combined Monte Carlo/torsion-angle molecular dynamics for ensemble modeling of proteins, nucleic acids and carbohydrates

We describe a general method to use Monte Carlo simulation followed by torsion-angle molecular dynamics simulations to create ensembles of structures to model a wide variety of soft-matter biological systems. Our particular emphasis is focused on modeling low-resolution small-angle scattering and reflectivity structural data. We provide examples of this method applied to HIV-1 Gag protein and derived fragment proteins, TraI protein, linear B-DNA, a nucleosome core particle, and a glycosylated monoclonal antibody. This procedure will enable a large community of researchers to model low-resolution experimental data with greater accuracy by using robust physics based simulation and sampling methods which are a significant improvement over traditional methods used to interpret such data.

Contrast Variation SANS Study of DNA Packaging in Bacteriophage λ

The elegant physical structure of viruses has attracted wide-spread interests from structural biologists, biophysicists, and mathematicians. While the protein capsids of viruses posses the nearly universal quasi-equivalence of icosahedral symmetry, the configuration of viral DNA/RNA genome appears less structured and is still poorly understood. Here we report a contrast variation small angle neutron scattering (SANS) study of bacteriophage λ to probe the spatial organization of its dsDNA genome. Contrast variation SANS takes advantage of the different scattering lengths of protein and DNA, and can experimentally match the protein contrast so as to measure the DNA structure alone. Together with the known structure of the protein capsid, SANS data at different H/D ratios offer the opportunity for extensive structural modeling. Furthermore, knowledge of DNA configuration in bacteriophage λ impinges directly on the physics of pressurized DNA packaging known to drive the initial DNA ejection when infecting bacteria. Multivalent cations were used to modulate the internal DNA pressure, and the induced structural changes were measured by contrast variation SANS. We compare the various DNA packaging models proposed in the literature and discuss their consistencies and discrepancies with experimental data.