Order and interactions in DNA arrays: Multiscale molecular dynamics simulation

Impact of Divalent Cations on In-Layer Positional Order of DNA-Based Liquid Crystals: Implications for DNA Condensation

The layered liquid crystalline phases formed by DNA molecules, which include rigid and flexible segments ("gapped DNA"), enable the study of both end-to-end stacking and side-to-side (helix-to-helix) lateral interactions, forming a model system to study such interactions at physiologically relevant DNA and ion concentrations. The observed layer structure exhibits long-range interlayer and in-layer positional correlations. In particular, the in-layer order has implications for DNA condensation, as it reflects whether these normally repulsive interactions become attractive under certain ionic conditions. Using synchrotron small-angle X-ray scattering measurements, we investigate the impact of divalent Mg2+ cations (in addition to a constant 150 mM Na+) on the stability of the inter- and in-layer DNA ordering as a function of temperature between 5 and 65 °C. DNA constructs with different terminal base pairings were created to mediate the strength of the attractive end-to-end stacking interactions between the blunt ends of the gapped DNA constructs. We demonstrate that the stabilities at a fixed DNA concentration of both interlayer and in-layer order are significantly enhanced even at a few mM Mg2+ concentration. The stabilities are even higher at 30 mM Mg2+; however, a marked decrease is observed at 100 mM Mg2+, suggesting a change in the nature of side-by-side interactions within this Mg2+ concentration range. We discuss the implications of these results in terms of counterion-mediated DNA-DNA attraction and DNA condensation.

Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media

Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.

Tuning the Dielectric Response of Water in Nanoconfinement through Surface Wettability

The tunable polarity of water can be exploited in emerging technologies including catalysis, gas storage, and green chemistry. Recent experimental and theoretical studies have shown that water can be rendered into an effectively apolar solvent under nanoconfinement. We furthermore demonstrate, through molecular simulations, that the static dielectric constant of water can be modified by changing the wettability of the confining material. We find the out-of-plane dielectric response to be highly sensitive to the level of confinement and can be reduced up to 40× , in accordance with experimental data. By altering the surface wettability from superhydrophilic to superhydrophobic, we observe a 36% increase for the out-of-plane and a 31% decrease for the in-plane dielectric constants. Our findings demonstrate the feasibility of tunable water polarity, a phenomenon with great potential for scientific and technological impact.

From adaptive resolution to molecular dynamics of open systems

ABSTRACT

We provide an overview of the Adaptive Resolution Simulation method (AdResS) based on discussing its basic principles and presenting its current numerical and theoretical developments. Examples of applications to systems of interest to soft matter, chemical physics, and condensed matter illustrate the method's advantages and limitations in its practical use and thus settle the challenge for further future numerical and theoretical developments.

Additive Modulation of DNA-DNA Interactions by Interstitial Ions

Quantitative understanding of biomolecular electrostatics, particularly involving multivalent ions and highly charged surfaces, remains lacking. Ion-modulated interactions between nucleic acids provide a model system in which electrostatics plays a dominant role. Using ordered DNA arrays neutralized by spherical cobalt³⁺ hexammine and Mg²⁺ ions, we investigate how the interstitial ions modulate DNA-DNA interactions. Using methods of ion counting, osmotic stress, and x-ray diffraction, we systematically determine thermodynamic quantities, including ion chemical potentials, ion partition, DNA osmotic pressure and force, and DNA-DNA spacing. Analyses of the multidimensional data provide quantitative insights into their interdependencies. The key finding of this study is that DNA-DNA forces are observed to linearly depend on the partition of interstitial ions, suggesting the dominant role of ion-DNA coupling. Further implications are discussed in light of physical theories of electrostatic interactions and like-charge attraction.

Liquid crystal ordering of nucleic acids

Several analytical calculations and computer simulations propose that cylindrical monodispersive rods having an aspect ratio (ratio of length to diameter) greater than 4 can exhibit liquid crystal (LC) ordering. But, recent experiments demonstrated the signature of LC ordering in systems of 4- to 20-base pair (bp) long nucleic acids (NAs) that do not satisfy the shape anisotropy criterion. Mechanisms of end-to-end adhesion and stacking have been proposed to explain this phenomenon. In this study, using all-atom molecular dynamics (MD) simulation, we explicitly verify the end-to-end stacking of double-stranded RNA (dsRNA) and demonstrate the LC ordering at the microscopic level. Using umbrella sampling (US) calculation, we quantify the potential of mean force (PMF) between two dsRNAs for various reaction coordinates (RCs) and compare our results with previously reported PMFs for double-stranded DNA (dsDNA). The PMF profiles demonstrate the anisotropic nature of inter-NA interaction. We find that, like dsDNA, dsRNA also prefers to stack on top of each other while repelling sideways, leading to the formation of supra-molecular-columns that undergo LC ordering at high NA volume fraction (φ). We also demonstrate and quantify the nematic ordering of the RNAs using several hundred nanosecond-long MD simulations that remain almost invariant for different initial configurations and under different external physiological conditions.

Multivalent counterions induced attraction between DNA polyelectrolytes

In this paper we study the electrostatic attraction between two parallel rodlike DNA polyelectrolytes induced by neutralizing multivalent counterions at the zero temperature limit. The counterions crystallize on the charged surfaces of DNA so that we can handle the system by using the Wigner crystal lattice model. We derived the 3D ground state configuration of counterions with minimized energy by use of the gradient descent method, and calculated the interaction between two DNA cylinders with divalent or trivalent counterions when they approach. The results show that the complex ground state configuration of counterions plays a key role in determining the caused attraction. The counterions form three-dimensional Wigner crystals on each cylinder at large separation. When the cylinders are brought together, some counterion lines will move towards the inner region and lead to strong attraction. The calculated interaction from our model is in good agreement with the simulation result, however, the single particle approximation considerably overestimates the attraction.

Counterions in the ground state partially condense in the intervening region of two approaching DNA cylinders and result in attraction.

Dynamical aspects of supercooled TIP3P-water in the grooves of DNA

We investigate by molecular dynamics simulations the mobility of the water located at the DNA minor and major grooves. We employ the TIP3P water model, and our system is analyzed for a range of temperatures 190-300 K. For high temperatures, the water at the grooves shows an Arrhenius behavior similar to that observed in the bulk water. At lower temperatures, a departure from the bulk behavior is observed. This slowing down in the dynamics is compared with the dynamics of the hydrogen of the DNA at the grooves and with the autocorrelation functions of the water hydrogen bonds. Our results indicate that the hydrogen bonds of the water at the minor grooves are highly correlated, which suggests that this is the mechanism for the slow dynamics at this high confinement.

SWINGER: a clustering algorithm for concurrent coupling of atomistic and supramolecular liquids

In this contribution, we review recent developments and applications of a dynamic clustering algorithm SWINGER tailored for the multiscale molecular simulations of biomolecular systems. The algorithm on-the-fly redistributes solvent molecules among supramolecular clusters. In particular, we focus on its applications in combination with the adaptive resolution scheme, which concurrently couples atomistic and coarse-grained molecular representations. We showcase the versatility of our multiscale approach on a few applications to biomolecular systems coupling atomistic and supramolecular water models such as the well-established MARTINI and dissipative particle dynamics models and provide an outlook for future work.

The physics of open systems for the simulation of complex molecular environments in soft matter

Molecular dynamics (MD) has become one of the most powerful tools of investigation in soft matter. Despite such success, simulations of large molecular environments are mostly run using the approximation of closed systems without the possibility of exchange of matter. Due to the molecular complexity of soft matter systems, an optimal simulation strategy would require the application of concurrent multiscale resolution approaches such that each part of a large system can be considered as an open subsystem at a high resolution embedded in a large coarser reservoir of energy and particles. This paper discusses the current capability and the future perspectives of multiscale adaptive resolution MD methods to satisfy the conceptual principles of open systems and to perform simulations of complex molecular environments in soft matter.

Splay-density coupling in semiflexible main-chain nematic polymers with hairpins

A main-chain nematic polymer melt/solution exhibits macroscopic orientational order of main polymer chains, i.e., a preferred (nematic) direction. It has long been known that in such polymeric liquid crystals spatial density/concentration variations and distortions of the nematic direction are coupled, obeying a vectorial continuity constraint whose rigidity increases with chain length. Its vectorial nature precludes the application to flexible chains, where backfolds (hairpins) are present and apolar nematic symmetry is manifest, which has been its puzzling feature from the beginning. We now establish a description of the splay-density coupling in the case of arbitrary backfolding, devising a continuity constraint for the "recovered" polar order of the chain tangents and introducing hairpins as its new type of sources. Performing detailed Monte Carlo simulations of nematic monodomain melts of "soft" worm-like chains with variable length and flexibility, we show via their structure factors that the weakening of the coupling due to the backfolding can be consistently quantified on the macroscopic level.

Adaptive resolution molecular dynamics technique: Down to the essential

We investigate the role of the thermodynamic (TD) force as an essential and sufficient technical ingredient for an efficient and accurate adaptive resolution algorithm. Such a force applied in the coupling region of an adaptive resolution molecular dynamics setup assures thermodynamic equilibrium between atomistically resolved and coarse-grained regions, allowing the proper exchange of molecules. We numerically prove that indeed for systems as relevant as liquid water and 1,3-dimethylimidazolium chloride ionic liquid, the combined action of the TD force and thermostat allows for computationally efficient and numerically accurate simulations, beyond the current capabilities of adaptive resolution setups, which employ switching functions in the coupling region.

Open-Boundary Molecular Dynamics of a DNA Molecule in a Hybrid Explicit/Implicit Salt Solution

The composition and electrolyte concentration of the aqueous bathing environment have important consequences for many biological processes and can profoundly affect the behavior of biomolecules. Nevertheless, because of computational limitations, many molecular simulations of biophysical systems can be performed only at specific ionic conditions: either at nominally zero salt concentration, i.e., including only counterions enforcing the system's electroneutrality, or at excessive salt concentrations. Here, we introduce an efficient molecular dynamics simulation approach for an atomistic DNA molecule at realistic physiological ionic conditions. The simulations are performed by employing the open-boundary molecular dynamics method that allows for simulation of open systems that can exchange mass and linear momentum with the environment. In our open-boundary molecular dynamics approach, the computational burden is drastically alleviated by embedding the DNA molecule in a mixed explicit/implicit salt-bathing solution. In the explicit domain, the water molecules and ions are both overtly present in the system, whereas in the implicit water domain, only the ions are explicitly present and the water is described as a continuous dielectric medium. Water molecules are inserted and deleted into/from the system in the intermediate buffer domain that acts as a water reservoir to the explicit domain, with both water molecules and ions free to enter or leave the explicit domain. Our approach is general and allows for efficient molecular simulations of biomolecules solvated in bathing salt solutions at any ionic strength condition.

Multiscale Simulation of Protein Hydration Using the SWINGER Dynamical Clustering Algorithm

To perform computationally efficient concurrent multiscale simulations of biological macromolecules in solution, where the all-atom (AT) models are coupled to supramolecular coarse-grained (SCG) solvent models, previous studies resorted to modified AT water models, such as the bundled-simple point charge (SPC) models, that use semiharmonic springs to restrict the relative movement of water molecules within a cluster. Those models can have a significant impact on the simulated biomolecules and can lead, for example, to a partial unfolding of a protein. In this work, we employ the recently developed alternative approach with a dynamical clustering algorithm, SWINGER, which enables a direct coupling of original unmodified AT and SCG water models. We perform an adaptive resolution molecular dynamics simulation of a Trp-Cage miniprotein in multiscale water, where the standard SPC water model is interfaced with the widely used MARTINI SCG model, and demonstrate that, compared to the corresponding full-blown AT simulations, the structural and dynamic properties of the solvated protein and surrounding solvent are well reproduced by our approach.