Characterization of the NISTmAb Reference Material using small-angle scattering and molecular simulation : Part I: Dilute protein solutions.

Both conformational and colloidal stability of therapeutic proteins must be closely monitored and thoroughly characterized to assess the long-term viability of drug products. We characterized the IgG1 NISTmAb reference material in its histidine formulation buffer and report our findings on the higher order structure and interactions of NISTmAb under a range of conditions. In this paper we present the analysis of experimental small-angle scattering data with atomistic molecular simulations to characterize the monodisperse dilute solution of NISTmAb. In part II we describe the characterization of the NISTmAb at high protein concentration (Castellanos et al. 2018). The NISTmAb was found to be a flexible protein with a radius of gyration of 49.0 ± 1.2 Å in histidine formulation buffer using a variety of neutron and X-ray scattering measurements. Scattering data were then modeled using molecular simulation. After building and validating a starting NISTmAb structure from the Fc and Fab crystallographic coordinates, molecular dynamics and torsion-angle Monte Carlo simulations were performed to explore the configuration space sampled in the NISTmAb and obtain ensembles of structures with atomistic detail that are consistent with the experimental data. Our results indicate that the small-angle scattering profiles of the NISTmAb can be modeled using ensembles of flexible structures that explore a wide configuration space. The NISTmAb is flexible in solution with no single preferred orientation of Fc and Fab domains, but with some regions of configuration space that are more consistent with measured scattering profiles. Analysis of inter-domain atomistic contacts indicated that all ensembles contained configurations where residues between domains are ≤ 4 Å, although few contacts were observed for variable and C H 3 regions. Graphical Abstract Heavy atom self contact maps of the NISTmAb indicate a highly-flexible structure.

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. © 2016 Wiley Periodicals, Inc.

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.

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.

Electrostatically driven interactions between hybrid DNA-carbon nanotubes.

Single-stranded DNA is able to wrap around single-wall carbon nanotubes (CNT) and form stable DNA-CNT hybrids that are highly soluble in solution. Here we report quantitative measurements and analysis of the interactions between DNA-CNT hybrids at low salts. Condensation of DNA-CNT hybrids by neutral osmolytes leads to liquid crystalline phases, and varying the osmotic pressure modulates the interhybrid distance that is determined by x-ray diffraction. Thus obtained force-distance dependencies of DNA-CNT hybrids show a remarkable resemblance to that of double-stranded DNA with differences that can be largely accounted for by their different diameters. This establishes their common physical nature of electrostatically driven interactions. Quantitative modeling further reveals the roles of hydration in mediating the interhybrid forces within the last nanometer of surface separation. This study also suggests the utility of osmotic pressure to control DNA-CNT assemblies at subnanometer precision.

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.

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.

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.