Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

The effect of machine learning predicted anharmonic frequencies on thermodynamic properties of fluid hydrogen fluoride

Anharmonic effects play a crucial role in determining thermochemical properties of liquids and gases. For such extended phases, the inclusion of anharmonicity in reliable electronic structure methods is computationally extremely demanding, and hence, anharmonic effects are often lacking in thermochemical calculations. In this study, we apply the quantum cluster equilibrium method to transfer density functional theory calculations at the cluster level to the macroscopic, liquid, and gaseous phase of hydrogen fluoride. This allows us to include anharmonicity, either via vibrational self-consistent field calculations for smaller clusters or using a regression model for larger clusters. We obtain the structural composition of the fluid phases in terms of the population of different clusters as well as isobaric heat capacities as an example for thermodynamic properties. We study the role of anharmonicities for these analyses and observe that, in particular, the dominating structural motifs are rather sensitive to the anharmonicity in vibrational frequencies. The regression model proves to be a promising way to get access to anharmonic features, and the extension to more sophisticated machine-learning models is promising.

Coarse-Grained Water Model Development for Accurate Dynamics and Structure Prediction

Several coarse-graining (CG) methods have been combined to develop a CG model of water capable of the accurate prediction of structure and dynamics properties. The multiscale coarse-graining (MS-CG) method based on force matching and the PDF-based coarse-graining method were used for accurate dynamics prediction. The iterative Boltzmann inversion (IBI) method was added for accurate structure representation. The approach is applied to bulk water, and the results show close reproduction of the CG structure when compared with the reference atomistic data. The combination of MS-CG and IBI methods facilitates the development of CG force fields at different temperatures based on a single MS-CG coarse-graining procedure. The dynamic properties of the CG water model closely match those obtained from the reference atomistic system. The general application of this approach to any existing coarse-graining methods is discussed.

Coarse-Grained Dynamically Accurate Simulations of Ionic Liquids: [pyr14][TFSI] and [EMIM][BF4]

In this work, coarse-grained (CG) models for two different sets of ionic liquids were developed from atomistic molecular dynamics (MD) reference systems, expanding their system size and time duration capabilities. The bonded force field of the CG systems was built using harmonic oscillator potential (HOP) fitting, while the nonbonded force field was generated with the multiscale coarse-graining (MS-CG) approach based on force matching. The dynamics of each system were corrected using the probability distribution function-based coarse-grained molecular dynamics (PDF-based CGMD) method. The structure and dynamics of each system were proven to match reference system data at two temperature scales. CG models and force fields for these liquids were developed to exemplify a general purpose methodology for producing MD results of ionic liquids and other fluids with accurate structural as well as dynamic properties. As an application, developed ionic liquids CG models were then applied to study vacuum-interface interaction. Density profile results of vacuum-interface exposure show significant deviation from bulk behavior. At the interface, multilayer ordering of ionic liquids is predicted to be similar to those observed from an experimental work. This ordering is intensified by decreasing temperature and use of the PDF-based CGMD method as opposed to conventional CG methods.

Hidden Halogen-Bonding Ability of Fluorine Manifesting in the Hydrogen-Bond Configurations of Hydrogen Fluoride

Elucidating how the intermolecular interactions of a covalently bonded fluorine atom are similar to and different from those of the other halogen atoms will be helpful for a better unified understanding of them. In the present study, the case of hydrogen fluoride is theoretically studied from this viewpoint by using the techniques of electron density analysis, molecular dynamics of liquid, and others. It is shown that the extra-point model, which locates an additional charge site on the line extended from (not within) the covalent bond and has been adopted for halogen-bonding systems as a key to the generation of proper stability and directionality, works well also in this case. A significantly bent hydrogen-bond configuration, which is characteristic of the intermolecular interactions of hydrogen fluoride, is reasonably well reproduced, meaning that it is a manifestation of the latent halogen-bonding ability, which is hidden by the strongly electronegative nature.

Computational and Experimental Approaches to Investigate Lipid Nanoparticles as Drug and Gene Delivery Systems

Lipid nanoparticles (LNPs) have been widely applied in drug and gene delivery. More than twenty years ago, DoxilTM was the first LNPs-based drug approved by the US Food and Drug Administration (FDA). Since then, with decades of research and development, more and more LNP-based therapeutics have been used to treat diverse diseases, which often offer the benefits of reduced toxicity and/or enhanced efficacy compared to the active ingredients alone. Here, we provide a review of recent advances in the development of efficient and robust LNPs for drug/gene delivery. We emphasize the importance of rationally combining experimental and computational approaches, especially those providing multiscale structural and functional information of LNPs, to the design of novel and powerful LNP-based delivery systems.

New Molecular-Mechanics Model for Simulations of Hydrogen Fluoride in Chemistry and Biology

Hydrogen fluoride (HF) is the most polar diatomic molecule and one of the simplest molecules capable of hydrogen-bonding. HF deviates from ideality both in the gas phase and in solution and is thus of great interest from a fundamental standpoint. Pure and aqueous HF solutions are broadly used in chemical and industrial processes, despite their high toxicity. HF is a stable species also in some biological conditions, because it does not readily dissociate in water unlike other hydrogen halides; yet, little is known about how HF interacts with biomolecules. Here, we set out to develop a molecular-mechanics model to enable computer simulations of HF in chemical and biological applications. This model is based on a comprehensive high-level ab initio quantum chemical investigation of the structure and energetics of the HF monomer and dimer; (HF)_n clusters, for n = 3-7; various clusters of HF and H₂O; and complexes of HF with analogs of all 20 amino acids and of several commonly occurring lipids, both neutral and ionized. This systematic analysis explains the unique properties of this molecule: for example, that interacting HF molecules favor nonlinear geometries despite being diatomic and that HF is a strong H-bond donor but a poor acceptor. The ab initio data also enables us to calibrate a three-site molecular-mechanics model, with which we investigate the structure and thermodynamic properties of gaseous, liquid, and supercritical HF in a wide range of temperatures and pressures; the solvation structure of HF in water and of H₂O in liquid HF; and the free diffusion of HF across a lipid bilayer, a key process underlying the high cytotoxicity of HF. Despite its inherent simplifications, the model presented significantly improves upon previous efforts to capture the properties of pure and aqueous HF fluids by molecular-mechanics methods and to our knowledge constitutes the first parameter set calibrated for biomolecular simulations.

Computational Characterization of a Cholesterol-Based Molecular Rotor in Lipid Membranes

Biophysical properties of cellular membranes critically depend on their content of cholesterol and its interaction with various other lipid species. Cholesterol-dependent friction at the nanoscale can be studied with molecular rotors, whose quantum yield depends on rotational dynamics of functional groups during their excited state lifetime. Here, we present a detailed computational analysis of a phenyl-BODIPY-linked cholesterol based molecular rotor in direct comparison with the well-known TopFluor-cholesterol. We describe a new parametrization strategy of force field parameters for the BODIPY moiety and carry out extensive molecular dynamics simulations of the probe in membranes in the absence or presence of cholesterol. Our study quantifies the extent of membrane perturbation by these probes, analyzes their tilting resistance in the bilayer and derives dynamic properties directly related to the rotor propensity. We show that phenyl-BODIPY-cholesterol bears potential as a cholesterol-dependent molecular rotor to report about microviscosity of sterol-containing model and cell membranes.

Development of the ReaxFF Methodology for Electrolyte-Water Systems

A new ReaxFF reactive force field has been developed for water-electrolyte systems including cations Li⁺, Na⁺, K⁺, and Cs⁺ and anions F^-, Cl^-, and I^-. The reactive force field parameters have been trained against quantum mechanical (QM) calculations related to water binding energies, hydration energies and energies of proton transfer. The new force field has been validated by applying it to molecular dynamics (MD) simulations of the ionization of different electrolytes in water and comparison of the results with experimental observations and thermodynamics. Radial distribution functions (RDF) determined for most of the atom pairs (cation or anion with oxygen and hydrogen of water) show a good agreement with the RDF values obtained from DFT calculations. On the basis of the applied force field, the ReaxFF simulations have described the diffusion constants for water and electrolyte ions in alkali metal hydroxide and chloride salt solutions as a function of composition and electrolyte concentration. The obtained results open opportunities to advance ReaxFF methodology to a wide range of applications involving electrolyte ions and solutions.

From classical to quantum and back: Hamiltonian adaptive resolution path integral, ring polymer, and centroid molecular dynamics

Path integral-based methodologies play a crucial role for the investigation of nuclear quantum effects by means of computer simulations. However, these techniques are significantly more demanding than corresponding classical simulations. To reduce this numerical effort, we recently proposed a method, based on a rigorous Hamiltonian formulation, which restricts the quantum modeling to a small but relevant spatial region within a larger reservoir where particles are treated classically. In this work, we extend this idea and show how it can be implemented along with state-of-the-art path integral simulation techniques, including path-integral molecular dynamics, which allows for the calculation of quantum statistical properties, and ring-polymer and centroid molecular dynamics, which allow the calculation of approximate quantum dynamical properties. To this end, we derive a new integration algorithm that also makes use of multiple time-stepping. The scheme is validated via adaptive classical-path-integral simulations of liquid water. Potential applications of the proposed multiresolution method are diverse and include efficient quantum simulations of interfaces as well as complex biomolecular systems such as membranes and proteins.

A Coarse-Grained Molecular Dynamics Study of Carbon Nanoparticle Aggregation

A multiscale coarse-graining procedure is used to study carbonaceous nanoparticle assembly. The computational methodology is applied to an ensemble of 10 000 nanoparticles (or effectively 2 million total carbon atoms) to simulate the agglomeration of carbonaceous nanoparticles using coarse-grained atomistic-scale information. In particular, with the coarse-graining approach, we are able to assess the influence of nanoparticle morphology and temperature on the agglomeration process. The coarse-graining of the interparticle force field is accomplished applying a force-matching procedure to data obtained from trajectories and forces from all-atom molecular dynamics simulation. The coarse-grained molecular dynamics results show rich and varied clustering behaviors for different particle morphologies. They are shown to reproduce accurately the structural properties of the nanoparticles systems studied, while allowing for molecular dynamics simulations of much larger self-assembled nanoparticles systems.

Multiscale Coarse-Graining of Mixed Phospholipid/Cholesterol Bilayers

Coarse-grained (CG) models for mixed dimyristoylphosphatidylcholine (DMPC)/cholesterol lipid bilayers are constructed using the recently developed multiscale coarse-graining (MS-CG) method. The MS-CG method permits a systematic fit of the bonded and nonbonded interactions and system pressure to trajectory and force data derived from an underlying reference all-atom molecular dynamics (MD) simulation. The CG sites for lipid and cholesterol molecules are associated with the centers-of-mass of atomic groups because of the simplicity in the evaluation of the forces acting on them from the atomistic MD data. Corresponding models with four-site and seven-site representations of the cholesterol molecule were also developed. The latter CG models differed by the bonding scheme of CG sites to represent intramolecular interactions. A one-site MS-CG model based on the TIP3P potential was used for water, with the interaction site placed at the molecular geometrical center, and the analytical fit of the model is presented. The MS-CG models were then used to conduct simulations in the constant NPT ensemble which reproduce accurately the structural properties as seen in the full all-atom MD simulation.

Efficient Multistate Reactive Molecular Dynamics Approach Based on Short-Range Effective Potentials

Nonbonded interactions between molecules usually include the van der Waals force and computationally expensive long-range electrostatic interactions. This article develops a more efficient approach: the effective-interaction multistate empirical-valence-bond (EI-MS-EVB) model. The EI-MS-EVB method relies on a mapping of all interactions onto a short-range and thus, computationally efficient effective potential. The effective potential is tabulated by matching its force to known trajectories obtained from the full-potential empirical multistate empirical-valence-bond (MS-EVB) model. The effective pairwise interaction depends on and is uniquely determined by the atomic configuration of the system, varying only with respect to the hydrogen-bonding topology. By comparing the EI-MS-EVB and full MS-EVB calculations of several equilibrium and dynamic properties important to hydrated excess proton solvation and transport, we show that the EI-MS-EVB model produces very accurate results for the specific system in which the tabulated potentials were generated. The EI-MS-EVB potential also transfers reasonably well to similar systems with different temperatures and box sizes. The EI-MS-EVB method also reduces the computational cost of the nonbonded interactions by about 1 order of magnitude in comparison with the full algorithm.

Quantum mechanical force field for hydrogen fluoride with explicit electronic polarization

The explicit polarization (X-Pol) theory is a fragment-based quantum chemical method that explicitly models the internal electronic polarization and intermolecular interactions of a chemical system. X-Pol theory provides a framework to construct a quantum mechanical force field, which we have extended to liquid hydrogen fluoride (HF) in this work. The parameterization, called XPHF, is built upon the same formalism introduced for the XP3P model of liquid water, which is based on the polarized molecular orbital (PMO) semiempirical quantum chemistry method and the dipole-preserving polarization consistent point charge model. We introduce a fluorine parameter set for PMO, and find good agreement for various gas-phase results of small HF clusters compared to experiments and ab initio calculations at the M06-2X/MG3S level of theory. In addition, the XPHF model shows reasonable agreement with experiments for a variety of structural and thermodynamic properties in the liquid state, including radial distribution functions, interaction energies, diffusion coefficients, and densities at various state points.

Developing accurate molecular mechanics force fields for conjugated molecular systems

A rapid method to parameterize the intramolecular component of classical force fields is proposed and applied to a molecular semiconductor, oligomers of conjugated polymers and a biological chromophore.

Generalized QM/MM Force Matching Approach Applied to the 11-cis Protonated Schiff Base Chromophore of Rhodopsin

We extended a previously developed force matching approach to systems with covalent QM/MM boundaries and describe its user-friendly implementation in the publicly available software package CPMD. We applied this approach to the challenging case of the retinal protonated Schiff base in dark state bovine rhodopsin. We were able to develop a highly accurate force field that is able to capture subtle structural changes within the chromophore that have a pronounced influence on the optical properties. The optical absorption spectrum calculated from configurations extracted from a MD trajectory using the new force field is in excellent agreement with QM/MM and experimental references.

Multiscale reactive molecular dynamics

Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data directly from condensed phase ab initio simulations to develop reactive molecular dynamics models that do not require predefined empirical functions. Instead, the interactions used in the reactive model are expressed as linear combinations of interpolating functions that are optimized by using a linear least-squares algorithm. One notable benefit of the procedure outlined here is the capability to minimize the number of parameters requiring nonlinear optimization. The method presented can be generally applied to multiscale problems and is demonstrated by generating reactive models for the hydrated excess proton and hydroxide ion based directly on condensed phase ab initio molecular dynamics simulations. The resulting models faithfully reproduce the water-ion structural properties and diffusion constants from the ab initio simulations. Additionally, the free energy profiles for proton transfer, which is sensitive to the structural diffusion of both ions in water, are reproduced. The high fidelity of these models to ab initio simulations will permit accurate modeling of general chemical reactions in condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus facilitating a proper statistical sampling of the coupling to slow, large-scale motions of the system.

Coarse-Grained Molecular Models of Water: A Review

Coarse-grained (CG) models have proven to be very effective tools in the study of phenomena or systems that involve large time- and length-scales. By decreasing the degrees of freedom in the system and using softer interactions than seen in atomistic models, larger timesteps can be used and much longer simulation times can be studied. CG simulations are widely used to study systems of biological importance that are beyond the reach of atomistic simulation, necessitating a computationally efficient and accurate CG model for water. In this review, we discuss the methods used for developing CG water models and the relative advantages and disadvantages of the resulting models. In general, CG water models differ with regards to how many waters each CG group or bead represents, whether analytical or tabular potentials have been used to describe the interactions, and how the model incorporates electrostatic interactions. Finally, how the models are parameterized depends on their application, so, while some are fitted to experimental properties such as surface tension and density, others are fitted to radial distribution functions extracted from atomistic simulations.

The curious case of the hydrated proton

Understanding the hydrated proton is a critically important problem that continues to engage the research efforts of chemists, physicists, and biologists because of its involvement in a wide array of phenomena. Only recently have several unique properties of the hydrated proton been unraveled through computer simulations. One such process is the detailed molecular mechanism by which protons hop between neighboring water molecules, thus giving rise to the anomalously high diffusion of protons relative to other simple cations. Termed Grotthuss shuttling, this process occurs over multiple time and length scales, presenting unique challenges for computer modeling and simulation. Because the hydrated proton is in reality a dynamical electronic charge defect that spans multiple water molecules, the simulation methodology must be able to dynamically readjust the chemical bonding topology. This reactive nature of the chemical process is automatically captured with ab initio molecular dynamics (AIMD) simulation methods, where the electronic degrees of freedom are treated explicitly. Unfortunately, these calculations can be prohibitively expensive for more complex proton solvation and transport phenomena in the condensed phase. These AIMD simulations remain extremely valuable, however, in validating empirical models, verifying results, and providing insight into molecular mechanisms. In this Account, we discuss recent progress in understanding the solvation and transport properties of the hydrated excess proton. The advances are based on results obtained from reactive molecular dynamics simulations using the multistate empirical valence bond (MS-EVB) methodology. This approach relies on a dynamic linear combination of chemical bond topologies to model charge delocalization and dynamic bonding environments. When parametrized via a variational force-matching algorithm from AIMD trajectories, the MS-EVB method can be viewed as a multiscale bridging of ab initio simulation results to a simpler and more efficient representation. The process allows sampling of longer time and length scales, which would normally be too computationally expensive with AIMD alone. With the MS-EVB methodology, the statistically important components of the excess proton solvation and hopping mechanisms in liquid water have been identified. The most likely solvation structure for the hydrated proton is a distorted Eigen-type complex (H(9)O(4)(+)). In this state, the excess proton charge defect rapidly resonates between three possible distorted Eigen-type structures until a successful proton hop occurs. This process, termed the "special-pair dance", serves as a kind of preparatory phase for the proton hopping while the neighboring water hydrogen-bonding network fluctuates and ultimately rearranges to facilitate a proton hop. The modifications of the solvation structure and transport properties of the excess proton in concentrated acid solutions were further investigated. The Eigen-type solvation structure also possesses both "hydrophilic" and "hydrophobic" sides, which accounts for the affinity of the hydrated proton for the air-water interface. This unusual "amphiphilic" character of the hydrated proton further leads to the metastable formation of contact ion pairs between two hydrated protons. It also engenders a surprisingly constant degree of solubility of hydrophobic species as a function of acid concentration, which contrasts with a markedly variable solubility as a function of salt (such as NaCl or KCl) concentration.

Communication: Hybrid ensembles for improved force matching

Force matching is a method for parameterizing empirical potentials in which the empirical parameters are fitted to a reference potential energy surface (PES). Typically, training data are sampled from a canonical ensemble generated with either the empirical potential or the reference PES. In this Communication, we show that sampling from either ensemble risks excluding critical regions of configuration space, leading to fitted potentials that deviate significantly from the reference PES. We present a hybrid ensemble which combines the Boltzmann probabilities of both potential surfaces into the fitting procedure, and we demonstrate that this technique improves the quality and stability of empirical potentials.

Molecular dynamics simulations of polyglutamine aggregation using solvent-free multiscale coarse-grained models

The multiscale coarse-graining (MS-CG) method is used to construct solvent-free CG models for polyglutamine peptides having various repeat lengths. Because the resulting CG models have fewer degrees of freedom than a corresponding all-atom simulations, they make it possible to study the self-assembly of polyglutamines at high concentrations for the first time by allowing for better equilibration and statistical sampling that is well beyond the range achievable by all-atom models. Molecular dynamics (MD) simulations performed with these models show that polyglutamine monomers with repeat lengths < or = 28 fluctuate between their folded and unfolded states. Monomers with 32 or more residues are stable and form alpha-helix solid structures. The degree of monomer compactness increases with chain length in both cases. CG MD simulations of equilibrium polyglutamine aggregates show that even at high concentrations, the system statistically fluctuates between heterogeneous and homogeneous configurations, rather than simply aggregates. The degree of aggregation and fluctuation increases with concentration and chain length. All of these phenomena are consistent with the experimental observations and may be explained by a mechanism that the collective nonbonded interactions between polyglutamine molecules in water solution are only weakly attractive. Finally, this work demonstrates that computer simulation of polypeptides self-assembly and aggregation, which is presently beyond the reach of all-atom MD simulations, is attainable using solvent-free MS-CG models.

The multiscale coarse-graining method. VI. Implementation of three-body coarse-grained potentials

Many methodologies have been proposed to build reliable and computationally fast coarse-grained potentials. Typically, these force fields rely on the assumption that the relevant properties of the system under examination can be reproduced using a pairwise decomposition of the effective coarse-grained forces. In this work it is shown that an extension of the multiscale coarse-graining technique can be employed to parameterize a certain class of two-body and three-body force fields from atomistic configurations. The use of explicit three-body potentials greatly improves the results over the more commonly used two-body approximation. The method proposed here is applied to develop accurate one-site coarse-grained water models.

Statistical thermodynamics of biomembranes

An overview of the major issues involved in the statistical thermodynamic treatment of phospholipid membranes at the atomistic level is summarized: thermodynamic ensembles, initial configuration (or the physical system being modeled), force field representation as well as the representation of long-range interactions. This is followed by a description of the various ways that the simulated ensembles can be analyzed: area of the lipid, mass density profiles, radial distribution functions (RDFs), water orientation profile, deuterium order parameter, free energy profiles and void (pore) formation; with particular focus on the results obtained from our recent molecular dynamic (MD) simulations of phospholipids interacting with dimethylsulfoxide (Me(2)SO), a commonly used cryoprotective agent (CPA).

Are many-body electronic polarization effects important in liquid water?

Many-body electronic polarization effects may be important for an accurate description of aqueous environments. As a result, numerous polarizable water models have been developed to include explicit polarization effects in intermolecular potential functions. In this paper, it is shown for liquid water at ambient conditions that such many-body polarization interactions can be decomposed into effective pairwise contributions using the force-matching (FM) method [Izvekov et al., J. Chem. Phys. 120, 10896 (2004)]. It is found that an effective pairwise water model obtained by the FM method can accurately reproduce various bulk structural and thermodynamic properties obtained from an accurate fully polarizable water model. In addition, the effective pairwise water model also provides a reasonable description of the water liquid-vapor interface, thus exhibiting a degree of transferability to heterogeneous environments. These results suggest that the role and importance of many-body electronic polarization effects in aqueous systems might be fruitfully explored relative to the best possible pairwise decomposable bulk phase model as the reference state.

Effective potentials from complex simulations: a potential-matching algorithm and remarks on coarse-grained potentials

The projection of complex interactions onto simple distance-dependent or angle-dependent classical mechanical functions is a long-standing theoretical challenge in the field of computational sciences concerning biomolecules, colloids, aggregates and simple systems as well. The construction of an effective potential may be based on theoretical assumptions, on the application of fitting procedures on experimental data and on the simplification of complex molecular simulations. Recently, a force-matching method was elaborated to project the data of Car-Parrinello ab initio molecular dynamics simulations onto two-particle classical interactions (Izvekov et al 2004 J. Chem. Phys. 120 10896). We have developed a potential-matching algorithm as a practical analogue of this force-matching method. The algorithm requires a large number of configurations (particle positions) and a single value of the potential energy for each configuration. We show the details of the algorithm and the test calculations on simple systems. The test calculation on water showed an example in which a similar structure was obtained for qualitatively different pair interactions. The application of the algorithm on reverse Monte Carlo configurations was tried as well. We detected inconsistencies in a part of our calculations. We found that the coarse graining of potentials cannot be performed perfectly both for the structural and the thermodynamic data. For example, if one applies an inverse method with an input of the pair-correlation function, it provides energetics data for the configurations uniquely. These energetics data can be different from the desired ones obtained by all atom simulations, as occurred in the testing of our potential-matching method.

Correlated atomic motions in liquid deuterium fluoride studied by coherent quasielastic neutron scattering

The collective dynamics of liquid deuterium fluoride are studied by means of high-resolution quasielastic and inelastic neutron scattering over a range of four decades in energy transfer. The spectra show a low-energy coherent quasielastic component which arises from correlated stochastic motions as well as a broad inelastic feature originating from overdamped density oscillations. While these results are at variance with previous works which report on the presence of propagating collective modes, they are fully consistent with neutron diffraction, nuclear magnetic resonance, and infrared/Raman experiments on this prototypical hydrogen-bonded fluid.

Multiscale coarse-graining of ionic liquids

A recently developed multiscale coarse-graining (MS-CG) approach for obtaining coarse-grained force fields from fully atomistic molecular dynamics simulation is applied to the challenging case of the EMIM+NO3- ionic liquid. The force-matching in the MS-CG methodology is accomplished with an explicit separation of bonded and nonbonded forces. While the nonbonded forces are adopted from this force-matching approach, the bonded forces are obtained from fitting the statistical configurational data from the atomistic simulations. The many-body electronic polarizability is also successfully broken into effective pair interactions. With a virial constraint fixing the system pressure, the MS-CG models rebuild satisfactory structural and thermodynamic properties for different temperatures. The MS-CG model developed from a modest atomistic simulation is therefore suitable for simulating much larger systems, because the coarse-grained models show significant time integration efficiency. This approach is expected to be general for coarse-graining other ionic liquids, as well as many other liquid-state systems. The limitations of the present coarse-graining procedure are also discussed.

Coarse-grained peptide modeling using a systematic multiscale approach

A systematic new approach to derive multiscale coarse-grained (MS-CG) models has been recently developed. The approach employs information from atomistically detailed simulations to derive CG forces and associated effective potentials. In this work, the MS-CG methodology is extended to study two peptides representing distinct structural motifs, alpha-helical polyalanine and the beta-hairpin V(5)PGV(5). These studies represent the first known application of this approach to peptide systems. Good agreement between the MS-CG and atomistic models is achieved for several structural properties including radial distribution functions, root mean-square deviation, and radius of gyration. The new MS-CG models are able to preserve the native states of these peptides within approximately 1 A backbone root mean-square deviation during CG simulations. The MS-CG approach, as with most coarse-grained models, has the potential to increase the length and timescales accessible to molecular simulations. However, it is also able to maintain a clear connection to the underlying atomistic-scale interactions.