Scaling molecular dynamics beyond 100,000 processor cores for large-scale biophysical simulations.

The growing interest in the complexity of biological interactions is continuously driving the need to increase system size in biophysical simulations, requiring not only powerful and advanced hardware but adaptable software that can accommodate a large number of atoms interacting through complex forcefields. To address this, we developed and implemented strategies in the GENESIS molecular dynamics package designed for large numbers of processors. Long-range electrostatic interactions were parallelized by minimizing the number of processes involved in communication. A novel algorithm was implemented for nonbonded interactions to increase single instruction multiple data (SIMD) performance, reducing memory usage for ultra large systems. Memory usage for neighbor searches in real-space nonbonded interactions was reduced by approximately 80%, leading to significant speedup. Using experimental data describing physical 3D chromatin interactions, we constructed the first atomistic model of an entire gene locus (GATA4). Taken together, these developments enabled the first billion-atom simulation of an intact biomolecular complex, achieving scaling to 65,000 processes (130,000 processor cores) with 1 ns/day performance. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

Mesoscale modeling reveals formation of an epigenetically driven HOXC gene hub.

Gene expression is orchestrated at the structural level by nucleosome positioning, histone tail acetylation, and linker histone (LH) binding. Here, we integrate available data on nucleosome positioning, nucleosome-free regions (NFRs), acetylation islands, and LH binding sites to "fold" in silico the 55-kb HOXC gene cluster and investigate the role of each feature on the gene's folding. The gene cluster spontaneously forms a dynamic connection hub, characterized by hierarchical loops which accommodate multiple contacts simultaneously and decrease the average distance between promoters by ∼100 nm. Contact probability matrices exhibit "stripes" near promoter regions, a feature associated with transcriptional regulation. Interestingly, while LH proteins alone decrease long-range contacts and acetylation alone increases transient contacts, combined LH and acetylation produce long-range contacts. Thus, our work emphasizes how chromatin architecture is coordinated strongly by epigenetic factors and opens the way for nucleosome resolution models incorporating epigenetic modifications to understand and predict gene activity.

Chromatin Fiber Folding Directed by Cooperative Histone Tail Acetylation and Linker Histone Binding.

In eukaryotic chromatin, islands of histone tail acetylation are found near transcription start sites and enhancers, whereas linker histones (LHs) are localized in intergenic regions with wild-type (WT) histone tails. However, the structural mechanisms by which acetylation, in combination with LH binding, modulates chromatin compaction and hence transcription regulation are unknown. To explore the folding propensity by which these features may govern gene expression, we analyze 20 kb fibers that contain regularly spaced acetylation islands of two sizes (2 or 5 kb) with various LH levels by mesoscale modeling. Specifically, we investigate the effect of acetylating each histone tail (H3, H4, H2A, and H2B) individually, in combination (H3 and H4, or all tails), and adding LH to WT regions. We find that fibers with acetylated H4 tails lose local contacts (<1 kb) and fibers with all tails acetylated have decreased long-range contacts in those regions. Tail interaction plots show that this opening of the fiber is driven by the loss of tail-tail interactions in favor of tail-parent core interactions and/or increase in free tails. When adding LH to WT regions, the fibers undergo hierarchical looping, enriching long-range contacts between WT and acetylated domains. For reference, adding LH to the entire fiber results in local condensation and loss of overall long-range contacts. These findings highlight the cooperation between histone tail acetylation and regulatory proteins like LH in directing folding and structural heterogeneity of chromatin fibers. The results advance our understanding of chromatin contact domains, which represent a pivotal part of the cell cycle, diseased states, and differentiation states in eukaryotic cells.

Cohesin Loss Eliminates All Loop Domains.

The human genome folds to create thousands of intervals, called "contact domains," that exhibit enhanced contact frequency within themselves. "Loop domains" form because of tethering between two loci-almost always bound by CTCF and cohesin-lying on the same chromosome. "Compartment domains" form when genomic intervals with similar histone marks co-segregate. Here, we explore the effects of degrading cohesin. All loop domains are eliminated, but neither compartment domains nor histone marks are affected. Loss of loop domains does not lead to widespread ectopic gene activation but does affect a significant minority of active genes. In particular, cohesin loss causes superenhancers to co-localize, forming hundreds of links within and across chromosomes and affecting the regulation of nearby genes. We then restore cohesin and monitor the re-formation of each loop. Although re-formation rates vary greatly, many megabase-sized loops recovered in under an hour, consistent with a model where loop extrusion is rapid.

Kilobase Pair Chromatin Fiber Contacts Promoted by Living-System-Like DNA Linker Length Distributions and Nucleosome Depletion.

Nucleosome placement, or DNA linker length patterns, are believed to yield specific spatial features in chromatin fibers, but details are unknown. Here we examine by mesoscale modeling how kilobase (kb) range contacts and fiber looping depend on linker lengths ranging from 18 to 45 bp, with values modeled after living systems, including nucleosome free regions (NFRs) and gene encoding segments. We also compare artificial constructs with alternating versus randomly distributed linker lengths in the range of 18-72 bp. We show that nonuniform distributions with NFRs enhance flexibility and encourage kb-range contacts. NFRs between neighboring gene segments diminish short-range contacts between flanking nucleosomes, while enhancing kb-range contacts via hierarchical looping. We also demonstrate that variances in linker lengths enhance such contacts. In particular, moderate sized variations in fiber linker lengths (∼27 bp) encourage long-range contacts in randomly distributed linker length fibers. Our work underscores the importance of linker length patterns, alongside bound proteins, in biological regulation. Contacts formed by kb-range chromatin folding are crucial to gene activity. Because we find that special linker length distributions in living systems promote kb contacts, our work suggests ways to manipulate these patterns for regulation of gene activity.

Linking Chromatin Fibers to Gene Folding by Hierarchical Looping.

While much is known about DNA structure on the basepair level, this scale represents only a fraction of the structural levels involved in folding the genomic material. With recent advances in experimental and theoretical techniques, a variety of structures have been observed on the fiber, gene, and chromosome levels of genome organization. Here we view chromatin architecture from nucleosomes and fibers to genes and chromosomes, highlighting the rich structural diversity and fiber fluidity emerging from both experimental and theoretical techniques. In this context, we discuss our recently proposed folding mechanism, which we call "hierarchical looping", similar to rope flaking used in mountain climbing, where 10-nm zigzag chromatin fibers are compacted laterally into self-associating loops which then stack and fold in space. We propose that hierarchical looping may act as a bridge between fibers and genes as well as provide a mechanism to relate key features of interphase and metaphase chromosome architecture to genome structural changes. This motif emerged by analysis of ultrastructural internucleosome contact data by electron microscopy-assisted nucleosome interaction capture cross-linking experiments, in combination with mesoscale modeling. We suggest that while the local folding of chromatin can be regulated at the fiber level by adjustment of internal factors such as linker-histone binding affinities, linker DNA lengths, and divalent ion levels, hierarchical looping on the gene level can additionally be controlled by posttranslational modifications and external factors such as polycomb group proteins. From a combination of 3C data and mesoscale modeling, we suggest that hierarchical looping could also play a role in epigenetic gene silencing, as stacked loops may occlude access to transcription start sites. With advances in crystallography, single-molecule in vitro biochemistry, in vivo imaging techniques, and genome-wide contact data experiments, various modeling approaches are allowing for previously unavailable structural interpretation of these data at multiple spatial and temporal scales. An unprecedented level of productivity and opportunity is on the horizon for the chromatin structure field.

Mesoscale Modeling Reveals Hierarchical Looping of Chromatin Fibers Near Gene Regulatory Elements.

While it is well-recognized that chromatin loops play an important role in gene regulation, structural details regarding higher order chromatin loops are only emerging. Here we present a systematic study of restrained chromatin loops ranging from 25 to 427 nucleosomes (fibers of 5-80 Kb DNA in length), mimicking gene elements studied by 3C contact data. We find that hierarchical looping represents a stable configuration that can effectively bring distant regions of the GATA-4 gene together, satisfying connections reported by 3C experiments. Additionally, we find that restrained chromatin fibers larger than 100 nucleosomes (∼20Kb) form closed plectonemes, whereas fibers shorter than 100 nucleosomes form simple hairpin loops. By studying the dependence of loop structures on internal parameters, we show that loop features are sensitive to linker histone concentration, loop length, divalent ions, and DNA linker length. Specifically, increasing loop length, linker histone concentration, and divalent ion concentration are associated with increased persistence length (or decreased bending), while varying DNA linker length in a manner similar to experimentally observed "nucleosome free regions" (found near transcription start sites) disrupts intertwining and leads to loop opening and increased persistence length in linker histone depleted (-LH) fibers. Chromatin fiber structure sensitivity to these parameters, all of which vary throughout the cell cycle, tissue type, and species, suggests that caution is warranted when using uniform polymer models to fit chromatin conformation capture genome-wide data. Furthermore, the folding geometry we observe near the transcription initiation site of the GATA-4 gene suggests that hierarchical looping provides a structural mechanism for gene inhibition, and offers tunable parameters for design of gene regulation elements.

Hierarchical looping of zigzag nucleosome chains in metaphase chromosomes.

The architecture of higher-order chromatin in eukaryotic cell nuclei is largely unknown. Here, we use electron microscopy-assisted nucleosome interaction capture (EMANIC) cross-linking experiments in combination with mesoscale chromatin modeling of 96-nucleosome arrays to investigate the internal organization of condensed chromatin in interphase cell nuclei and metaphase chromosomes at nucleosomal resolution. The combined data suggest a novel hierarchical looping model for chromatin higher-order folding, similar to rope flaking used in mountain climbing and rappelling. Not only does such packing help to avoid tangling and self-crossing, it also facilitates rope unraveling. Hierarchical looping is characterized by an increased frequency of higher-order internucleosome contacts for metaphase chromosomes compared with chromatin fibers in vitro and interphase chromatin, with preservation of a dominant two-start zigzag organization associated with the 30-nm fiber. Moreover, the strong dependence of looping on linker histone concentration suggests a hierarchical self-association mechanism of relaxed nucleosome zigzag chains rather than longitudinal compaction as seen in 30-nm fibers. Specifically, concentrations lower than one linker histone per nucleosome promote self-associations and formation of these looped networks of zigzag fibers. The combined experimental and modeling evidence for condensed metaphase chromatin as hierarchical loops and bundles of relaxed zigzag nucleosomal chains rather than randomly coiled threads or straight and stiff helical fibers reconciles aspects of other models for higher-order chromatin structure; it constitutes not only an efficient storage form for the genomic material, consistent with other genome-wide chromosome conformation studies that emphasize looping, but also a convenient organization for local DNA unraveling and genome access.

Single-Walled Carbon Nanotubes Modulate the B- to A-DNA Transition.

We study the conformational equilibrium between B-to-A forms of ds-DNA adsorbed onto a single-walled carbon nanotube (SWNT) using free energy profile calculations based on all-atom molecular dynamics simulations. The potential of mean force (PMF) of the B-to-A transition of ds-DNA in the presence of an uncharged (10,0) carbon nanotube for two dodecamers with poly-AT or poly-GC sequences is calculated as a function of a root-mean-square-distance (ΔRMSD) difference metric for the B-to-A transition. The calculations reveal that in the presence of a SWNT DNA favors B-form DNA significantly in both poly-GC and poly-AT sequences. Furthermore, the poly-AT DNA:SWNT complex shows a higher energy penalty for adopting an A-like conformation than poly-GC DNA:SWNT by several kcal/mol. The presence of a SWNT on either poly-AT or poly-GC DNA affects the PMF of the transition such that the B form is favored by as much as 10 kcal/mol. In agreement with published data, we find a potential energy minimum between A and B-form DNA at ΔRMSD ≈ -1.5 Å and that the presence of the SWNT moves this minimum by as much as ΔRMSD = 3 Å.

Design, synthesis, and activity of a series of arylpyrid-3-ylmethanones as type I positive allosteric modulators of α7 nicotinic acetylcholine receptors.

A series of novel arylpyrid-3-ylmethanones (7a-aa) were designed as modulators of α7 nicotinic acetylcholine receptors (nAChRs). The methanones were found to be type I positive allosteric modulators (PAMs) of human α7 nAChRs expressed in Xenopus ooctyes. Structure-activity relationship (SAR) studies resulted in the identification of compound 7v as a potent and efficacious type I PAM with maximum modulation of a nicotine EC5 response of 1200% and EC50 = 0.18 µM. Compound 7z was active in reversing the effect of scopolamine in the novel object recognition (NOR) paradigm with a minimum effective ip dose of 1.0 mg/kg (2.7 µmol/kg). This effect was blocked by the selective α7 nAChR antagonist methyllycaconitine (MLA). These compounds are potent type I positive allosteric modulators of α7 nAChRs that may have therapeutic value in restoring impaired sensory gating and cognitive deficits in schizophrenia and Alzheimer's disease.

A general method for constructing atomic-resolution RNA ensembles using NMR residual dipolar couplings: the basis for interhelical motions revealed.

The ability to modulate alignment and measure multiple independent sets of NMR residual dipolar couplings (RDCs) has made it possible to characterize internal motions in proteins at atomic resolution and with time scale sensitivity ranging from picoseconds up to milliseconds. The application of such methods to the study of RNA dynamics, however, remains fundamentally limited by the inability to modulate alignment and by strong couplings between internal and overall motions that complicate the quantitative interpretation of RDCs. Here, we address this problem by showing that RNA alignment can be generally modulated, in a controlled manner, by variable elongation of A-form helices and that the information contained within the measured RDCs can be extracted even in the presence of strong couplings between motions and overall alignment via structure-based prediction of alignment. Using this approach, four RDC data sets, and a broad conformational pool obtained from a 8.2 µs molecular dynamics simulation, we successfully construct and validate an atomic resolution ensemble of human immunodeficiency virus type I transactivation response element RNA. This ensemble reveals local motions in and around the bulge involving changes in stacking and hydrogen-bonding interactions, which are undetectable by traditional spin relaxation and drive global changes in interhelical orientation. This new approach broadens the scope of using RDCs in characterizing the dynamics of nucleic acids.

Probing sequence-specific DNA flexibility in a-tracts and pyrimidine-purine steps by nuclear magnetic resonance (13)C relaxation and molecular dynamics simulations.

Sequence-specific DNA flexibility plays a key role in a variety of cellular interactions that are critical for gene packaging, expression, and regulation, yet few studies have experimentally explored the sequence dependence of DNA dynamics that occur on biologically relevant time scales. Here, we use nuclear magnetic resonance (NMR) carbon spin relaxation combined with molecular dynamics (MD) simulations to examine the picosecond to nanosecond dynamics in a variety of dinucleotide steps as well as in varying length homopolymeric A(n)·T(n) repeats (A(n)-tracts, where n = 2, 4, or 6) that exhibit unusual structural and mechanical properties. We extend the NMR spin relaxation time scale sensitivity deeper into the nanosecond regime by using glycerol and a longer DNA duplex to slow overall tumbling. Our studies reveal a structurally unique A-tract core (for n > 3) that is uniformly rigid, flanked by junction steps that show increasing sugar flexibility with A-tract length. High sugar mobility is observed at pyrimidine residues at the A-tract junctions, which is encoded at the dinucleotide level (CA, TG, and CG steps) and increases with A-tract length. The MD simulations reproduce many of these trends, particularly the overall rigidity of A-tract base and sugar sites, and suggest that the sugar-backbone dynamics could involve transitions in sugar pucker and phosphate backbone BI ↔ BII equilibria. Our results reinforce an emerging view that sequence-specific DNA flexibility can be imprinted in dynamics occurring deep within the nanosecond time regime that is difficult to characterize experimentally at the atomic level. Such large-amplitude sequence-dependent backbone fluctuations might flag the genome for specific DNA recognition.