Density-Nematic Coupling in Isotropic Solution of DNA: Multiscale Model.

Monte Carlo simulations of isotropic solutions of double-stranded DNA (deoxyribonucleic acid) are performed using the well-established oxDNA model. By comparing the fluctuation amplitudes with theoretical predictions, the parameters of a generic macroscopic model of an isotropic linear polymer solution/melt are determined. A multiscale continuum field model is thus obtained, corresponding to the full specificity of the isotropic phase of double-stranded DNA in the usual B-form as perceived at the macroscopic level. Present research is particularly focused on the coupling between spatial concentration/density variations of the polymer and the emerging nematic orientation order of the chains. This rather unfamiliar, only recently described phenomenon, inherent to linear polymers, is outlined and interpreted. Quantitative predictions are provided for the degree of nematic order induced by concentration gradients in isotropic solutions of double-stranded DNA.

Lessons learned from urgent computing in Europe: Tackling the COVID-19 pandemic.

PRACE (Partnership for Advanced Computing in Europe), an international not-for-profit association that brings together the five largest European supercomputing centers and involves 26 European countries, has allocated more than half a billion core hours to computer simulations to fight the COVID-19 pandemic. Alongside experiments, these simulations are a pillar of research to assess the risks of different scenarios and investigate mitigation strategies. While the world deals with the subsequent waves of the pandemic, we present a reflection on the use of urgent supercomputing for global societal challenges and crisis management.

Sub-THz acoustic excitation of protein motion.

The application of terahertz radiation has been shown to affect both protein structure and cellular function. As the key to such structural changes lies in the dynamic response of a protein, we focus on the susceptibility of the protein's internal dynamics to mechanical stress induced by acoustic pressure waves. We use the open-boundary molecular dynamics method, which allows us to simulate the protein exposed to acoustic waves. By analyzing the dynamic fluctuations of the protein subunits, we demonstrate that the protein is highly susceptible to acoustic waves with frequencies corresponding to those of the internal protein vibrations. This is confirmed by changes in the compactness of the protein structure. As the amplitude of the pressure wave increases, even larger deviations from average positions and larger changes in protein compactness are observed. Furthermore, performing the mode-projection analysis, we show that the breathing-like character of collective modes is enhanced at frequencies corresponding to those used to excite the protein.

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.

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.

Adaptive resolution simulations of biomolecular systems.

In this review article, we discuss and analyze some recently developed hybrid atomistic-mesoscopic solvent models for multiscale biomolecular simulations. We focus on the biomolecular applications of the adaptive resolution scheme (AdResS), which allows solvent molecules to change their resolution back and forth between atomistic and coarse-grained representations according to their positions in the system. First, we discuss coupling of atomistic and coarse-grained models of salt solution using a 1-to-1 molecular mapping-i.e., one coarse-grained bead represents one water molecule-for development of a multiscale salt solution model. In order to make use of coarse-grained molecular models that are compatible with the MARTINI force field, one has to resort to a supramolecular mapping, in particular to a 4-to-1 mapping, where four water molecules are represented with one coarse-grained bead. To this end, bundled atomistic water models are employed, i.e., the relative movement of water molecules that are mapped to the same coarse-grained bead is restricted by employing harmonic springs. Supramolecular coupling has recently also been extended to polarizable coarse-grained water models with explicit charges. Since these coarse-grained models consist of several interaction sites, orientational degrees of freedom of the atomistic and coarse-grained representations are coupled via a harmonic energy penalty term. The latter aligns the dipole moments of both representations. The reviewed multiscale solvent models are ready to be used in biomolecular simulations, as illustrated in a few examples.

Deciphering the dynamics of star molecules in shear flow.

This work analyses the rotation of star polymers under shear flow, in melts, and in good solvent dilute solution. The latter is modeled by single molecule Brownian hydrodynamics, while melts are modeled using non-equilibrium molecular dynamics in closed (periodic) boxes and in open boundaries. A Dissipative Particle Dynamics (DPD) thermostat introduces pairwise monomer friction in melts at will, in directions normal and tangent to the monomer-monomer vectors. Although tangential friction is seldom modeled, we show that it is essential to control hydrodynamic effects in melts. We analyze the different sources of molecular angular momentum in solution and melts and distinguish three dynamic regimes as the shear rate [small gamma, Greek, dot above] is increased. These dynamic regimes are related with the disruption of the different relaxation mechanisms of the star in equilibrium. Although strong differences are found between harmonic springs and finitely extensible bonds, above a critical shear rate the star molecule has a "breathing" mode with successive elongations and contractions in the flow direction with frequency Ω. The force balance in the flow direction unveils a relation between Ω and the orientation angle. Using literature results for the tumbling of rings and linear chains, either in melt or in solution, we show that the relation is general. A different "tank-treading" dynamics determines the rotation of monomers around the center of mass of the molecule. We show that the tank-treading frequency does not saturate but keeps increasing with [small gamma, Greek, dot above]. This is at odds with previous studies which erroneously calculated the molecular angular frequency, used as a proxy for tank-treading.

Application of the Eckart frame to soft matter: rotation of star polymers under shear flow.

The Eckart co-rotating frame is used to analyze the dynamics of star polymers under shear flow, either in melt or solution and with different types of bonds. This formalism is compared with the standard approach used in many previous studies on polymer dynamics, where an apparent angular velocity ω is obtained from the relation between the tensor of inertia and angular momentum. A common mistake is to interpret ω as the molecular rotation frequency, which is only valid for rigid-body rotation. The Eckart frame, originally formulated to analyze the infrared spectra of small molecules, dissects different kinds of displacements: vibrations without angular momentum, pure rotation, and vibrational angular momentum (leading to a Coriolis cross-term). The Eckart frame co-rotates with the molecule with an angular frequency Ω obtained from the Eckart condition for minimal coupling between rotation and vibration. The standard and Eckart approaches are compared with a straight description of the star's dynamics taken from the time autocorrelation of the monomer positions moving around the molecule's center of mass. This is an underdamped oscillatory signal, which can be described by a rotation frequency ωR and a decorrelation rate Γ. We consistently find that Ω coincides with ωR, which determines the characteristic tank-treading rotation of the star. By contrast, the apparent angular velocity ω < Ω does not discern between pure rotation and molecular vibrations. We believe that the Eckart frame will be useful to unveil the dynamics of semiflexible molecules where rotation and deformations are entangled, including tumbling, tank-treading motions and breathing modes.

Adaptive resolution simulations coupling atomistic water to dissipative particle dynamics.

Multiscale methods are the most efficient way to address the interlinked spatiotemporal scales encountered in soft matter and molecular liquids. In the literature reported hybrid approaches span from quantum to atomistic, coarse-grained, and continuum length scales. In this article, we present the hybrid coupling of the molecular dynamics (MD) and dissipative particle dynamics (DPD) methods, bridging the micro- and mesoscopic descriptions. The interfacing is performed within the adaptive resolution scheme (AdResS), which is a linear momentum conserving coupling technique. Our methodology is hence suitable to simulate fluids on the micro/mesoscopic scale, where hydrodynamics plays an important role. The presented approach is showcased for water at ambient conditions. The supramolecular coupling is enabled by a recently developed clustering algorithm SWINGER that assembles, disassembles, and reassembles clusters as needed during the course of the simulation. This allows for a seamless coupling between standard atomistic MD and DPD models. The developed framework can be readily applied to various applications in the fields of materials and life sciences, e.g., simulations of phospholipids and polymer melts, or to study the red blood cells behavior in normal and disease states.

Open boundary molecular dynamics of sheared star-polymer melts.

Open boundary molecular dynamics (OBMD) simulations of a sheared star polymer melt under isothermal conditions are performed to study the rheology and molecular structure of the melt under a fixed normal load. Comparison is made with the standard molecular dynamics (MD) in periodic (closed) boxes at a fixed shear rate (using the SLLOD dynamics). The OBMD system exchanges mass and momentum with adjacent reservoirs (buffers) where the external pressure tensor is imposed. Insertion of molecules in the buffers is made feasible by implementing there a low resolution model (blob-molecules with soft effective interactions) and then using the adaptive resolution scheme (AdResS) to connect with the bulk MD. Straining with increasing shear stress induces melt expansion and a significantly different redistribution of pressure compared with the closed case. In the open sample, the shear viscosity is also a bit lowered but more stable against the viscous heating. At a given Weissenberg number, molecular deformations and material properties (recoverable shear strain and normal stress ratio) are found to be similar in both setups. We also study the modelling effect of normal and tangential friction between monomers implemented in a dissipative particle dynamics (DPD) thermostat. Interestingly, the tangential friction substantially enhances the elastic response of the melt due to a reduction of the kinetic stress viscous contribution.

STOCK: structure mapper and online coarse-graining kit for molecular simulations.

We present a web toolkit STructure mapper and Online Coarse-graining Kit for setting up coarse-grained molecular simulations. The kit consists of two tools: structure mapping and Boltzmann inversion tools. The aim of the first tool is to define a molecular mapping from high, for example, all-atom, to low, that is, coarse-grained, resolution. Using a graphical user interface it generates input files, which are compatible with standard coarse-graining packages, for example, Versatile Object-oriented Toolkit for Coarse-graining Applications and DL_CGMAP. Our second tool generates effective potentials for coarse-grained simulations preserving the structural properties, for example, radial distribution functions, of the underlying higher resolution model. The required distribution functions can be provided by any simulation package. Simulations are performed on a local machine and only the distributions are uploaded to the server. The applicability of the toolkit is validated by mapping atomistic pentane and polyalanine molecules to a coarse-grained representation. Effective potentials are derived for systems of TIP3P (transferable intermolecular potential 3 point) water molecules and salt solution. The presented coarse-graining web toolkit is available at http://stock.cmm.ki.si.

Adaptive resolution simulation of polarizable supramolecular coarse-grained water models.

Multiscale simulations methods, such as adaptive resolution scheme, are becoming increasingly popular due to their significant computational advantages with respect to conventional atomistic simulations. For these kind of simulations, it is essential to develop accurate multiscale water models that can be used to solvate biophysical systems of interest. Recently, a 4-to-1 mapping was used to couple the bundled-simple point charge water with the MARTINI model. Here, we extend the supramolecular mapping to coarse-grained models with explicit charges. In particular, the two tested models are the polarizable water and big multiple water models associated with the MARTINI force field. As corresponding coarse-grained representations consist of several interaction sites, we couple orientational degrees of freedom of the atomistic and coarse-grained representations via a harmonic energy penalty term. This additional energy term aligns the dipole moments of both representations. We test this coupling by studying the system under applied static external electric field. We show that our approach leads to the correct reproduction of the relevant structural and dynamical properties.

Adaptive resolution simulation of an atomistic protein in MARTINI water.

We present an adaptive resolution simulation of protein G in multiscale water. We couple atomistic water around the protein with mesoscopic water, where four water molecules are represented with one coarse-grained bead, farther away. We circumvent the difficulties that arise from coupling to the coarse-grained model via a 4-to-1 molecule coarse-grain mapping by using bundled water models, i.e., we restrict the relative movement of water molecules that are mapped to the same coarse-grained bead employing harmonic springs. The water molecules change their resolution from four molecules to one coarse-grained particle and vice versa adaptively on-the-fly. Having performed 15 ns long molecular dynamics simulations, we observe within our error bars no differences between structural (e.g., root-mean-squared deviation and fluctuations of backbone atoms, radius of gyration, the stability of native contacts and secondary structure, and the solvent accessible surface area) and dynamical properties of the protein in the adaptive resolution approach compared to the fully atomistically solvated model. Our multiscale model is compatible with the widely used MARTINI force field and will therefore significantly enhance the scope of biomolecular simulations.

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.

Rotational Dynamics of a Protein under Shear Flow Studied by the Eckart Frame Formalism.

Proteins are natural polymers that play an essential role in both living organisms and biotechnological applications. During certain bioprocessing steps, they can be exposed to significant mechanical stress induced by, for example, shear flow or sonication, resulting in reduced therapeutic efficacy, aggregation, or even a loss of activity. For this reason, there is a need to understand and determine the susceptibility of the protein activity to the experienced mechanical stress. To acquire this knowledge, it is necessary to study the rotational dynamics of the protein. Commonly, the rotational dynamics of soft molecules is interpreted based on a theoretical analysis performed in an inertial laboratory frame. However, the obtained angular velocity mixes pure rotations and vibrations with angular momentum, consequently lacking a clear dynamical interpretation. On the other hand, the use of the noninertial internal Eckart frame allows the determination of pure angular velocity as it minimizes the coupling between the rotational and vibrational degrees of freedom. In the present work, by conducting open-boundary molecular dynamics simulations and exploiting the Eckart frame formalism, we study the rotational dynamics of a small protein under the shear flow of various strengths. Our results show that the angular velocity increases nonlinearly with increasing shear rate. Furthermore, the protein gains vibrational angular momentum at higher shear rates, which is reflected in the higher angular velocity computed by employing the Eckart frame formalism and confirmed by analysis of the contributions to the total kinetic energy of the biomolecule.

Dissipative Particle Dynamics Simulation of Ultrasound Propagation through Liquid Water.

Ultrasound is widely used as a noninvasive method in therapeutic and diagnostic applications. These can be further optimized by computational approaches, as they allow for controlled testing and rational optimization of the ultrasound parameters, such as frequency and amplitude. Usually, continuum numerical methods are used to simulate ultrasound propagating through different tissue types. In contrast, ultrasound simulations using particle description are less common, as the implementation is challenging. In this work, a dissipative particle dynamics model is used to perform ultrasound simulations in liquid water. The effects of frequency and thermostat parameters are studied and discussed. We show that frequency and thermostat parameters affect not only the attenuation but also the computed speed of sound. The present study paves the way for development and optimization of a virtual ultrasound machine for large-scale biomolecular simulations.

Accelerated Simulations of Molecular Systems through Learning of Effective Dynamics.

Simulations are vital for understanding and predicting the evolution of complex molecular systems. However, despite advances in algorithms and special purpose hardware, accessing the time scales necessary to capture the structural evolution of biomolecules remains a daunting task. In this work, we present a novel framework to advance simulation time scales by up to 3 orders of magnitude by learning the effective dynamics (LED) of molecular systems. LED augments the equation-free methodology by employing a probabilistic mapping between coarse and fine scales using mixture density network (MDN) autoencoders and evolves the non-Markovian latent dynamics using long short-term memory MDNs. We demonstrate the effectiveness of LED in the Müller-Brown potential, the Trp cage protein, and the alanine dipeptide. LED identifies explainable reduced-order representations, i.e., collective variables, and can generate, at any instant, all-atom molecular trajectories consistent with the collective variables. We believe that the proposed framework provides a dramatic increase to simulation capabilities and opens new horizons for the effective modeling of complex molecular systems.

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.