Accurate Sequence-Dependent Coarse-Grained Model for Conformational and Elastic Properties of Double-Stranded DNA

CGeNArateWeb: a web server for the atomistic study of the structure and dynamics of chromatin fibers

We present CGeNArateWeb, a new web tool for the three-dimensional simulation of naked DNA and protein-bound chromatin fibers. The server allows the user to obtain a dynamic representation of long segments of linear, circular, or protein-DNA segments thanks to a Langevin dynamics coarse-grained (CG) model working with a machine-learning (ML) fitted C1'-resolution Hamiltonian. The CG trajectories can be back-mapped to atomistic resolution using another ML algorithm trained on a large database of molecular dynamics (MD) simulations. The method allows the user to get structural and dynamic information on large (kilobase range) portions of both protein-bound and free DNA, to transform conceptual cartoons into structural and dynamical models. Trajectories are analyzed using an extensive set of nucleic acid-specific analysis tools, and the results are displayed using a powerful and flexible graphic interface. The web tool uses state-of-the-art technologies such as (i) Docker components orchestrated by Docker Swarm, with containers deployed on demand for computations, (ii) WebGL-programmed NGL molecular viewer and the JavaScript plotly library for interactive plots, and (iii) noSQL-MongoDB for storage. The server is accessible at https://mmb.irbbarcelona.org/CGNAW/. The web tool is free and open to all users, and there are no login requirements.

Mechanisms of DNA Damage Recognition by UDG and PARP1 in the Nucleosome

The DNA base-excision repair (BER) pathway shares the second part of its enzymatic chain with the single-strand break (SSB) repair pathway. BER is initiated by a glycosylase, such as UDG, while SSBR is initiated by the multifunctional enzyme PARP1. The very early steps in the identification of the DNA damage are crucial to the correct initiation of the repair chains, and become even more complex when considering the realistic environment of damage to the DNA in the nucleosome. We performed molecular dynamics computer simulations of the interaction between the glycosylase UDG and a mutated uracil (as could result from oxidative deamination of cytosine), and between the Zn1-Zn2 fragment of PARP1 and a simulated SSB. The model system is a whole nucleosome in which DNA damage is inserted at various typical positions along the 145-bp sequence. It is shown that damage recognition by the enzymes requires very strict conditions, unlikely to be matched by pure random search along the DNA. We propose that mechanical deformation of the DNA around the defective sites may help signaling the presence of the defect, accelerating the search process.

Lipidic drug delivery systems are responsive to the human microbiome

In vitro and in vivo tests for therapeutic agents are typically conducted in sterile environments, but many target areas for drug delivery are home to thousands of microbial species. Here, we examine the behaviour of lipidic nanomaterials after exposure to representative strains of four bacterial species found in the gastrointestinal tract and skin. Small angle X-ray scattering measurements show that the nanostructure of monoolein cubic and inverse hexagonal phases are transformed, respectively, into inverse hexagonal and inverse micellar cubic phases upon exposure to a strain of live Staphylococcus aureus often present on skin and mucosa. Further investigation demonstrates that enzymatic hydrolysis and cell membrane lipid transfer are both likely responsible for this effect. The structural responses to S. aureus are rapid and significantly reduce the rate of drug release from monoolein-based nanomaterials. These findings are the first to demonstrate how a key species in the live human microbiome can trigger changes in the structure and drug release properties of lipidic nanomaterials. The effect appears to be strain specific, varies from patient to patient and body region to body region, and is anticipated to affect the bioapplication of monoglyceride-based formulations.

Mechanical response of double-stranded DNA: Bend, twist, and overwind

We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.

Shape and size tunability of sheets of interlocked ring copolymers

Mechanically bonded membranes of interlocked ring polymers are a significant generalization of conventional elastic sheets, where connectivity is provided by covalent bonding, and represent a promising class of topological meta-materials. In this context, two open questions regard the large-scale reverberations of the heterogeneous composition of the rings and the inequivalent modes of interlocking neighboring rings. We address these questions with Langevin dynamics simulations of chainmails with honeycomb-lattice connectivity, where the rings are block copolymers with two segments of different rigidity. We considered various combinations of the relative lengths of the two segments and the patterns of the over- and under-passes linking neighboring rings. We find that varying ring composition and linking patterns have independent and complementary effects. While the former sets the overall size of the chainmail, the latter defines the shape, enabling the selection of starkly different conformation types. Notably, one of the considered linking patterns favors saddle-shaped membranes, providing a first example of spontaneous negative Gaussian curvature in mechanically bonded sheets. The results help establish the extent to which mechanically bonded membranes can differ from conventional elastic ones, particularly for the achievable shape and size tunability.

CGeNArate: a sequence-dependent coarse-grained model of DNA for accurate atomistic MD simulations of kb-long duplexes

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles.

Skipping events impose repeated binding attempts: profound kinetic implications of protein-DNA conformational changes

The kinetics of protein-DNA recognition, along with its thermodynamic properties, including affinity and specificity, play a central role in shaping biological function. Protein-DNA recognition kinetics are characterized by two key elements: the time taken to locate the target site amid various nonspecific alternatives; and the kinetics involved in the recognition process, which may necessitate overcoming an energetic barrier. In this study, we developed a coarse-grained (CG) model to investigate interactions between a transcription factor called the sex-determining region Y (SRY) protein and DNA, in order to probe how DNA conformational changes affect SRY-DNA recognition and binding kinetics. We find that, not only does a requirement for such a conformational DNA transition correspond to a higher energetic barrier for binding and therefore slower kinetics, it may further impede the recognition kinetics by increasing unsuccessful binding events (skipping events) where the protein partially binds its DNA target site but fails to form the specific protein-DNA complex. Such skipping events impose the need for additional cycles protein search of nonspecific DNA sites, thus significantly extending the overall recognition time. Our results highlight a trade-off between the speed with which the protein scans nonspecific DNA and the rate at which the protein recognizes its specific target site. Finally, we examine molecular approaches potentially adopted by natural systems to enhance protein-DNA recognition despite its intrinsically slow kinetics.

Insights into elastic properties of coarse-grained DNA models: q-stiffness of cgDNA vs cgDNA

Coarse-grained models have emerged as valuable tools to simulate long DNA molecules while maintaining computational efficiency. These models aim at preserving interactions among coarse-grained variables in a manner that mirrors the underlying atomistic description. We explore here a method for testing coarse-grained vs all-atom models using stiffness matrices in Fourier space (q-stiffnesses), which are particularly suited to probe DNA elasticity at different length scales. We focus on a class of coarse-grained rigid base DNA models known as cgDNA and its most recent version, cgDNA+. Our analysis shows that while cgDNA+ closely follows the q-stiffnesses of the all-atom model, the original cgDNA shows some deviations for twist and bending variables, which are rather strong in the q → 0 (long length scale) limit. The consequence is that while both cgDNA and cgDNA+ give a suitable description of local elastic behavior, the former misses some effects that manifest themselves at longer length scales. In particular, cgDNA performs poorly on twist stiffness, with a value much lower than expected for long DNA molecules. Conversely, the all-atom and cgDNA+ twist are strongly length scale dependent: DNA is torsionally soft at a few base pair distances but becomes more rigid at distances of a few dozen base pairs. Our analysis shows that the bending persistence length in all-atom and cgDNA+ is somewhat overestimated.

Systematic Comparison of Atomistic Force Fields for the Mechanical Properties of Double-Stranded DNA

The response of double-stranded DNA to external mechanical stress plays a central role in its interactions with the protein machinery in the cell. Modern atomistic force fields have been shown to provide highly accurate predictions for the fine structural features of the duplex. In contrast, and despite their pivotal function, less attention has been devoted to the accuracy of the prediction of the elastic parameters. Several reports have addressed the flexibility of double-stranded DNA via all-atom molecular dynamics, yet the collected information is insufficient to have a clear understanding of the relative performance of the various force fields. In this work, we fill this gap by performing a systematic study in which several systems, characterized by different sequence contexts, are simulated with the most popular force fields within the AMBER family, bcs1 and OL15, as well as with CHARMM36. Analysis of our results, together with their comparison with previous work focused on bsc0, allows us to unveil the differences in the predicted rigidity between the newest force fields and suggests a roadmap to test their performance against experiments. In the case of the stretch modulus, we reconcile these differences, showing that a single mapping between sequence-dependent conformation and elasticity via the crookedness parameter captures simultaneously the results of all force fields, supporting the key role of crookedness in the mechanical response of double-stranded DNA.

Influence of temperature on bend, twist and twist-bend coupling of dsDNA

The temperature-dependent bend and twist elasticities of dsDNA, as well as their couplings, were explored through all-atom molecular dynamics simulations. Three rotational parameters, tilt, roll, and twist, were employed to assess the bend and twist elasticities through their stiffness matrix. Our analysis indicates that the bend and twist stiffnesses decrease as the temperature rises, primarily owing to entropic influences stemming from thermodynamic fluctuations. Furthermore, the couplings between these rotational parameters also exhibit a decline with increasing temperature, although the roll-twist coupling displays greater strength than the tilt-roll and tilt-twist couplings, attributed to its more robust correction component. We elucidated the influence of temperature on bend and twist elasticities based on the comparisons between various models and existing data.

Bi-directional nucleosome sliding by the Chd1 chromatin remodeler integrates intrinsic sequence-dependent and ATP-dependent nucleosome positioning

Chromatin remodelers use a helicase-type ATPase motor to shift DNA around the histone core. Although not directly reading out the DNA sequence, some chromatin remodelers exhibit a sequence-dependent bias in nucleosome positioning, which presumably reflects properties of the DNA duplex. Here, we show how nucleosome positioning by the Chd1 remodeler is influenced by local DNA perturbations throughout the nucleosome footprint. Using site-specific DNA cleavage coupled with next-generation sequencing, we show that nucleosomes shifted by Chd1 can preferentially localize DNA perturbations - poly(dA:dT) tracts, DNA mismatches, and single-nucleotide insertions - about a helical turn outside the Chd1 motor domain binding site, super helix location 2 (SHL2). This phenomenon occurs with both the Widom 601 positioning sequence and the natural +1 nucleosome sequence from the Saccharomyces cerevisiae SWH1 gene. Our modeling indicates that localization of DNA perturbations about a helical turn outward from SHL2 results from back-and-forth sliding due to remodeler action on both sides of the nucleosome. Our results also show that barrier effects from DNA perturbations can be extended by the strong phasing of nucleosome positioning sequences.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Force-dependent elasticity of nucleic acids

The functioning of double-stranded (ds) nucleic acids (NAs) in cellular processes is strongly mediated by their elastic response. These processes involve proteins that interact with dsDNA or dsRNA and distort their structures. The perturbation of the elasticity of NAs arising from these deformations is not properly considered by most theoretical frameworks. In this work, we introduce a novel method to assess the impact of mechanical stress on the elastic response of dsDNA and dsRNA through the analysis of the fluctuations of the double helix. Application of this approach to atomistic simulations reveals qualitative differences in the force dependence of the mechanical properties of dsDNA with respect to those of dsRNA, which we relate to structural features of these molecules by means of physically-sound minimalistic models.

Temperature dependence of DNA elasticity: An all-atom molecular dynamics simulation study

We used all-atom molecular dynamics simulation to investigate the elastic properties of double-stranded DNA (dsDNA). We focused on the influences of temperature on the stretch, bend, and twist elasticities, as well as the twist-stretch coupling, of the dsDNA over a wide range of temperature. The results showed that the bending and twist persistence lengths, together with the stretch and twist moduli, decrease linearly with temperature. However, the twist-stretch coupling behaves in a positive correction and enhances as the temperature increases. The potential mechanisms of how temperature affects dsDNA elasticity and coupling were investigated by using the trajectories from atomistic simulation, in which thermal fluctuations in structural parameters were analyzed in detail. We analyzed the simulation results by comparing them with previous simulation and experimental data, which are in good agreement. The prediction about the temperature dependence of dsDNA elastic properties provides a deeper understanding of DNA elasticities in biological environments and potentially helps in the further development of DNA nanotechnology.

Molecular dynamics simulations for the study of chromatin biology

The organization of Eukaryotic DNA into chromatin has profound implications for the processing of genetic information. In the past years, molecular dynamics (MD) simulations proved to be a powerful tool to investigate the mechanistic basis of chromatin biology. We review recent all-atom and coarse-grained MD studies revealing how the structure and dynamics of chromatin underlie its biological functions. We describe the latest method developments; the structural fluctuations of nucleosomes and the various factors affecting them; the organization of chromatin fibers, with particular emphasis on its liquid-like character; the interactions and dynamics of transcription factors on chromatin; and how chromatin organization is modulated by molecular motors acting on DNA.

Ab initio predictions for 3D structure and stability of single- and double-stranded DNAs in ion solutions

The three-dimensional (3D) structure and stability of DNA are essential to understand/control their biological functions and aid the development of novel materials. In this work, we present a coarse-grained (CG) model for DNA based on the RNA CG model proposed by us, to predict 3D structures and stability for both dsDNA and ssDNA from the sequence. Combined with a Monte Carlo simulated annealing algorithm and CG force fields involving the sequence-dependent base-pairing/stacking interactions and an implicit electrostatic potential, the present model successfully folds 20 dsDNAs (≤52nt) and 20 ssDNAs (≤74nt) into the corresponding native-like structures just from their sequences, with an overall mean RMSD of 3.4Å from the experimental structures. For DNAs with various lengths and sequences, the present model can make reliable predictions on stability, e.g., for 27 dsDNAs with/without bulge/internal loops and 24 ssDNAs including pseudoknot, the mean deviation of predicted melting temperatures from the corresponding experimental data is only ~2.0°C. Furthermore, the model also quantificationally predicts the effects of monovalent or divalent ions on the structure stability of ssDNAs/dsDNAs.

To determine 3D structures and quantify stability of single- (ss) and double-stranded (ds) DNAs is essential to unveil the mechanisms of their functions and to further guide the production and development of novel materials. Although many DNA models have been proposed to reproduce the basic structural, mechanical, or thermodynamic properties of dsDNAs based on the secondary structure information or preset constraints, there are very few models can be used to investigate the ssDNA folding or dsDNA assembly from the sequence. Furthermore, due to the polyanionic nature of DNAs, metal ions (e.g., Na⁺ and Mg²⁺) in solutions can play an essential role in DNA folding and dynamics. Nevertheless, ab initio predictions for DNA folding in ion solutions are still an unresolved problem. In this work, we developed a novel coarse-grained model to predict 3D structures and thermodynamic stabilities for both ssDNAs and dsDNAs in monovalent/divalent ion solutions from their sequences. As compared with the extensive experimental data and available existing models, we showed that the present model can successfully fold simple DNAs into their native-like structures, and can also accurately reproduce the effects of sequence and monovalent/divalent ions on structure stability for ssDNAs including pseudoknot and dsDNAs with/without bulge/internal loops.