Polarizable molecular dynamics simulations of ionic liquids: Influence of temperature control.
Ontology highlight
ABSTRACT:
SUBMITTER: Heid E
PROVIDER: S-EPMC7610910 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Developing advanced ionic electroactive devices such as ionic actuators and supercapacitors requires the understanding of ionic diffusion and drifting processes, which depend on the distances over which the ions travel, in these systems. The charge dynamics of [C(4)mim][PF(6)] ionic liquid films and Aquivion membranes with 40 wt% [C(2)mim][TfO] were investigated over a broad film thickness (d) range. It was found that the double layer charging time ?(DL) follows the classic model ?(DL) = ?(D)d/(2D) very well, where D is the diffusion coefficient and ?(D) the Debye length. In the longer time regimes (t ? ?(DL)) where diffusion dominates, the charge dynamics become voltage dependent. For low applied voltage, the later stage charge process seems to follow the d(2) dependence. However, at high voltages (> 0.5 V) in which significant device responses occur, the charging process does not show d(2) dependence so that ?(diff) = d(2)/(4D), corresponding to the ion diffusion from the bulk region, was not observed.
| S-EPMC3298447 | biostudies-literature
Project description:We present, for the first time, scalable polarizable molecular dynamics (MD) simulations within a polarizable continuum solvent with molecular shape cavities and exact solution of the mutual polarization. The key ingredients are a very efficient algorithm for solving the equations associated with the polarizable continuum, in particular, the domain decomposition Conductor-like Screening Model (ddCOSMO), which involves a rigorous coupling of the continuum with the polarizable force field achieved through a robust variational formulation and an effective strategy to solve the coupled equations. The coupling of ddCOSMO with nonvariational force fields, including AMOEBA, is also addressed. The MD simulations are feasible, for real-life systems, on standard cluster nodes; a scalable parallel implementation allows for further acceleration in the context of a newly developed module in Tinker, named Tinker-HP. NVE simulations are stable, and long-term energy conservation can be achieved. This paper is focused on the methodological developments, the analysis of the algorithm, and the stability of the simulations; a proof-of-concept application is also presented to attest to the possibilities of this newly developed technique.
| S-EPMC4621208 | biostudies-other
Project description:Additive force fields are designed to account for induced electronic polarization in a mean-field average way, using effective empirical fixed charges. The limitation of this approximation is cause for serious concerns, particularly in the case of lipid membranes, where the molecular environment undergoes dramatic variations over microscopic length scales. A polarizable force field based on the classical Drude oscillator offers a practical and computationally efficient framework for an improved representation of electrostatic interactions in molecular simulations. Building on the first-generation Drude polarizable force field for the dipalmitoylphosphatidylcholine 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) molecule, the present effort was undertaken to improve this initial model and expand the force field to a wider range of phospholipid molecules. New lipids parametrized include dimyristoylphosphatidylcholine (DMPC), dilauroylphosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), dipalmitoylphosphatidylethanolamine (DPPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE). The iterative optimization protocol employed in this effort led to lipid models that achieve a good balance between reproducing quantum mechanical data on model compound representative of phospholipids and reproducing a range of experimental condensed phase properties of bilayers. A parametrization strategy based on a restrained ensemble-maximum entropy methodology was used to help accurately match the experimental NMR order parameters in the polar headgroup region. All the parameters were developed to be compatible with the remainder of the Drude polarizable force field, which includes water, ions, proteins, DNA, and selected carbohydrates.
| S-EPMC5595662 | biostudies-other
Project description:A computational method is developed to allow molecular dynamics simulations of biomembrane systems under realistic ionic gradients and asymmetric salt concentrations while maintaining the conventional periodic boundary conditions required to minimize finite-size effects in an all-atom explicit solvent representation. The method, which consists of introducing a nonperiodic energy step acting on the ionic species at the edge of the simulation cell, is first tested with illustrative applications to a simple membrane slab model and a phospholipid membrane bilayer. The nonperiodic energy-step method is then used to calculate the reversal potential of the bacterial porin OmpF, a large cation-specific ?-barrel channel, by simulating the I-V curve under an asymmetric 10:1 KCl concentration gradient. The calculated reversal potential of 28.6 mV is found to be in excellent agreement with the values of 26-27 mV measured from lipid bilayer experiments, thereby demonstrating that the method allows realistic simulations of nonequilibrium membrane transport with quantitative accuracy. As a final example, the pore domain of Kv1.2, a highly selective voltage-activated K(+) channel, is simulated in a lipid bilayer under conditions that recreate, for the first time, the physiological K(+) and Na(+) concentration gradients and the electrostatic potential difference of living cells.
| S-EPMC3787774 | biostudies-literature
Project description:Today's standard molecular dynamics simulations of moderately sized biomolecular systems at full atomic resolution are typically limited to the nanosecond timescale and therefore suffer from limited conformational sampling. Efficient ensemble-preserving algorithms like replica exchange (REX) may alleviate this problem somewhat but are still computationally prohibitive due to the large number of degrees of freedom involved. Aiming at increased sampling efficiency, we present a novel simulation method combining the ideas of essential dynamics and REX. Unlike standard REX, in each replica only a selection of essential collective modes of a subsystem of interest (essential subspace) is coupled to a higher temperature, with the remainder of the system staying at a reference temperature, T(0). This selective excitation along with the replica framework permits efficient approximate ensemble-preserving conformational sampling and allows much larger temperature differences between replicas, thereby considerably enhancing sampling efficiency. Ensemble properties and sampling performance of the method are discussed using dialanine and guanylin test systems, with multi-microsecond molecular dynamics simulations of these test systems serving as references.
| S-EPMC1877756 | biostudies-literature
Project description:In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.
| S-EPMC6017617 | biostudies-literature
Project description:DNA and RNA sequences rich in guanine can fold into noncanonical structures called G-quadruplexes (GQs), which exhibit a common stem structure of Hoogsteen hydrogen-bonded guanine tetrads and diverse loop structures. GQ sequence motifs are overrepresented in promoters, origins of replication, telomeres, and untranslated regions in mRNA, suggesting roles in modulating gene expression and preserving genomic integrity. Given these roles and unique aspects of different structures, GQs are attractive targets for drug design, but greater insight into GQ folding pathways and the interactions stabilizing them is required. Here, we performed molecular dynamics simulations to study two bimolecular GQs, a telomeric DNA GQ and the analogous telomeric repeat-containing RNA (TERRA) GQ. We applied the Drude polarizable force field, which we show outperforms the additive CHARMM36 force field in both ion retention and maintenance of the GQ folds. The polarizable simulations reveal that the GQs bind bulk K+ ions differently, and that the TERRA GQ accumulates more K+ ions, suggesting different ion interactions stabilize these structures. Nucleobase dipole moments vary as a function of position and also contribute to ion binding. Finally, we show that the TERRA GQ is more sensitive than the telomeric DNA GQ to water-mediated modulation of ion-induced dipole-dipole interactions.
| S-EPMC6954416 | biostudies-literature
Project description:The accelerated molecular dynamics (aMD) method has recently been shown to enhance the sampling of biomolecules in molecular dynamics (MD) simulations, often by several orders of magnitude. Here, we describe an implementation of the aMD method for the OpenMM application layer that takes full advantage of graphics processing units (GPUs) computing. The aMD method is shown to work in combination with the AMOEBA polarizable force field (AMOEBA-aMD), allowing the simulation of long time-scale events with a polarizable force field. Benchmarks are provided to show that the AMOEBA-aMD method is efficiently implemented and produces accurate results in its standard parametrization. For the BPTI protein, we demonstrate that the protein structure described with AMOEBA remains stable even on the extended time scales accessed at high levels of accelerations. For the DNA repair metalloenzyme endonuclease IV, we show that the use of the AMOEBA force field is a significant improvement over fixed charged models for describing the enzyme active-site. The new AMOEBA-aMD method is publicly available (http://wiki.simtk.org/openmm/VirtualRepository) and promises to be interesting for studying complex systems that can benefit from both the use of a polarizable force field and enhanced sampling.
| S-EPMC3948463 | biostudies-literature
Project description:In this work, we have explored the thermophysical properties of tetraalkylammonium hydroxide ionic liquids (ILs) such as tetrapropylammonium hydroxide (TPAH) and tetrabutylammonium hydroxide (TBAH) with isomers of butanol (1-butanol, 2-butanol and 2-methyl-2-propanol) within the temperature range 293.15-313.15 K, with interval of 5 K and over the varied concentration range of ILs. The molecular interactions between ILs and butanol isomers are essential for understanding the function of ILs in related measures and excess functions are sensitive probe for the molecular interactions. Therefore, we calculated the excess molar volume (V(E) ) and the deviation in isentropic compressibility (??s ) using the experimental values such as densities (?) and ultrasonic sound velocities (u) that are measured over the whole compositions range at five different temperatures (293.15, 298.15, 303.15, 308.15 and 313.15 K) and atmospheric pressure. These excess functions were adequately correlated by using the Redlich-Kister polynomial equation. It was observed that for all studied systems, the V(E) and ??s values are negative for the whole composition range at 293.15 K. And, the excess function follows the sequence: 2-butanol>1-butanol>2-methyl-2-propanol, which reveals that (primary or secondary or tertiary) position of hydroxyl group influence the magnitude of interactions with ILs. The negative values of excess functions are contributions from the ion-dipole interaction, hydrogen bonding and packing efficiency between the ILs and butanol isomers. Hence, the position of hydroxyl group plays an important role in the interactions with ILs. The hydrogen bonding features between ILs and alcohols were analysed using molecular modelling program by using HyperChem 7.
| S-EPMC3906063 | biostudies-literature
Project description:Ionic Liquids are a broad group of salts with low melting points that can be specifically tuned for a broad range of applications. Despite being initially considered “green” solvents, their better environmental friendliness compared to traditional solvents has been increasingly challenged. In this study, we aimed to investigate the molecular effects of ILs exposure by using RNA-sequencing to study differential gene expression patterns. Thus, we exposed Daphnia magna to 1-ethyl-3-methylimidazolium chloride ([C2mim]Cl), 1-dodecyl chloride-3-methylimidazolium ([C12mim]Cl) and cholinium chloride ([Chol]Cl). Results suggest that the three ILs share several mechanisms of toxicity, including cellular membrane and cytoskeleton damage, oxidative stress, inhibition of antioxidant enzymes, mitochondrial affectation, changes in protein biosynthesis and energy production, DNA damage, and ultimately, programmed cell death and disease initiation. Overall, the dataset revealed that [C2mim]Cl and [C12mim]Cl were, respectively, the least and the most toxic ILs at the transcriptional level. Also, it is reinforced that [Chol]Cl is not devoid of environmental hazardous potential. Unique gene expression signatures could also be identified for each IL.
2021-08-13 | GSE156769 | GEO