Effect of Protein-Polarized Ligand Charges on Relative Protein Ligand Binding Affinities

A major challenge in computer-aided drug design is predicting relative binding energies of different molecules to a target protein using fast and accurate free-energy calculation methods. Free-energy calculations are primarily computed by utilizing classical molecular dynamics simulations based on all-atom force fields (FF) to model the interactions in the system. The present standard classical all-atom FFs contain fixed partial charges on the atoms, and hence electrostatic interactions are modeled between them. The parametrization process to determine these partial charges usually relies on quantum mechanics or semiempirical calculations of the molecule in the gas phase or homogeneous water surrounding. These present standard parametrization schemes of the partial charges neglect, therefore, polarization effects from the protein surrounding. The absence of protein polarization effects can lead to significant errors in free-energy calculations in proteins. We present a parametrization scheme for the partial charges of ligands, named protein-induced polarization (PIP) charges, which account for the electrostatic polarization due to the protein surrounding. The scheme involves single-point quantum mechanics/molecular mechanics calculations of the ligand charges in the protein/water surrounding. Using PIP ligand partial charges, we have calculated the relative binding free energies (RBFEs) of well-studied protein-ligand systems. We show here that RBFEs computed with PIP charges are either significantly improved or at least comparable to those computed with nonpolarized standard GAFF charges. Overall, we present a simple-to-use parametrization scheme to include protein polarization in any type of binding free-energy calculations. The parametrization scheme increases the accuracy in RBFE calculations, while it does not add significant computation time to standard parametrization procedures.

Effective Inclusion of Electronic Polarization Improves the Description of Electrostatic Interactions: The prosECCo75 Biomolecular Force Field

prosECCo75 is an optimized force field effectively incorporating electronic polarization via charge scaling. It aims to enhance the accuracy of nominally nonpolarizable molecular dynamics simulations for interactions in biologically relevant systems involving water, ions, proteins, lipids, and saccharides. Recognizing the inherent limitations of nonpolarizable force fields in precisely modeling electrostatic interactions essential for various biological processes, we mitigate these shortcomings by accounting for electronic polarizability in a physically rigorous mean-field way that does not add to computational costs. With this scaling of (both integer and partial) charges within the CHARMM36 framework, prosECCo75 addresses overbinding artifacts. This improves agreement with experimental ion binding data across a broad spectrum of systems─lipid membranes, proteins (including peptides and amino acids), and saccharides─without compromising their biomolecular structures. prosECCo75 thus emerges as a computationally efficient tool providing enhanced accuracy and broader applicability in simulating the complex interplay of interactions between ions and biomolecules, pivotal for improving our understanding of many biological processes.

Developing and Benchmarking Sulfate and Sulfamate Force Field Parameters via Ab Initio Molecular Dynamics Simulations To Accurately Model Glycosaminoglycan Electrostatic Interactions

Glycosaminoglycans (GAGs) are negatively charged polysaccharides found on cell surfaces, where they regulate transport pathways of foreign molecules toward the cell. The structural and functional diversity of GAGs is largely attributed to varied sulfation patterns along the polymer chains, which makes understanding their molecular recognition mechanisms crucial. Molecular dynamics (MD) simulations, thanks to their unmatched microscopic resolution, have the potential to be a reference tool for exploring the patterns responsible for biologically relevant interactions. However, the capability of molecular dynamics force fields used in biosimulations to accurately capture sulfation-specific interactions is not well established, partly due to the intrinsic properties of GAGs that pose challenges for most experimental techniques. In this work, we evaluate the performance of molecular dynamics force fields for sulfated GAGs by studying ion pairing of Ca2+ to sulfated moieties─N-methylsulfamate and methylsulfate─that resemble N- and O-sulfation found in GAGs, respectively. We tested available nonpolarizable (CHARMM36 and GLYCAM06) and explicitly polarizable (Drude and AMOEBA) force fields, and derived new implicitly polarizable models through charge scaling (prosECCo75 and GLYCAM-ECC75) that are consistent with our developed "charge-scaling" framework. The calcium-sulfamate/sulfate interaction free energy profiles obtained with the tested force fields were compared against reference ab initio molecular dynamics (AIMD) simulations, which serve as a robust alternative to experiments. AIMD simulations indicate that the preferential Ca2+ binding mode to sulfated GAG groups is solvent-shared pairing. Only our scaled-charge models agree satisfactorily with the AIMD data, while all other force fields exhibit poorer agreement, sometimes even qualitatively. Surprisingly, even explicitly polarizable force fields display a notable disagreement with the AIMD data, likely attributed to difficulties in their optimization and possible inherent limitations in depicting high-charge-density ion interactions accurately. Finally, the underperforming force fields lead to unrealistic aggregation of sulfated saccharides, which qualitatively disagrees with our understanding of the soft glycocalyx environment. Our results highlight the importance of accurately treating electronic polarization in MD simulations of sulfated GAGs and caution against over-reliance on currently available models without thorough validation and optimization.

Dressing a Nonpolarizable Force Field for OH- in TIP4P/2005 Aqueous Solutions with Corrected Hirshfeld Charges

We present a rigid model for the OH- ion parametrized for binary mixtures with TIP4P/2005-type water molecules. Li+, Na+ and K+ were selected as counterions, hence mimicking the important and widely used solutions of soluble alkaline hydroxides. The optimized atomic charge distributions were obtained by scaling in a factor of 0.85 those derived from the atomic dipole corrected Hirshfeld approach. The agreement between experimental and Molecular Dynamics simulation results is remarkable for a set of properties, namely, the dependence of the density of the solutions on the hydroxide concentration and on temperature, the structure (i.e., positions of the atom-to-atom radial distribution functions and coordination numbers), the viscosity coefficients, the surface tension, or the freezing point depression. The proposed optimized potential parameters for OH- thus enlarge the set of models comprised within the Madrid-2019 force field and widen the potential applicability of the TIP4P/2005 water model in basic environments.

Structure and Physicochemical Properties of Solutions of Lithium Polysulfides in Tetrasolvate of Lithium Perchlorate with Sulfolane Molecular Dynamics Modeling

Lithium polysulfides are intermediate products formed during the discharge and charge of lithium-sulfur batteries and have good solubility in electrolyte solutions. Therefore, the properties and structure of solutions of lithium polysulfides in electrolyte solutions affect the energy characteristics of the lithium-sulfur battery. In this work, the structure and physicochemical properties (density, ionic conductivity, and self-diffusion coefficients) of solutions of lithium disulfide, tetrasulfide, and octasulfide in lithium perchlorate tetrasolvate with sulfolane in a wide range of concentrations (0.2-12 M) were studied using the molecular dynamics method. With increasing concentrations of lithium polysulfides, the proportion of sulfur atoms in the coordination sphere of the lithium cation increases and the proportion of sulfolane molecules decreases. It has been established that the ability of the polysulfide anion to act as a bridging ligand leads to the formation of clusters, including lithium perchlorate, lithium polysulfides, and sulfolane. It has been shown that the tendency to form clusters increases with an increasing number of sulfur atoms in the polysulfide anion. At high concentrations of lithium polysulfides, regardless of their size, the electrolyte system becomes a single cluster. With an increase in the concentration of lithium polysulfides in the electrolyte system, its density increases and the ionic conductivity and diffusion coefficients of the system components decrease.

Diffusion and structure of propylene carbonate-metal salt electrolyte solutions for post-lithium-ion batteries: From experiment to simulation

The diffusion of cations in organic solvent solutions is important for the performance of metal-ion batteries. In this article, pulsed field gradient nuclear magnetic resonance experiments and fully atomistic molecular dynamic simulations were employed to study the temperature-dependent diffusive behavior of various liquid electrolytes representing 1M propylene carbonate solutions of metal salts with bis(trifluoromethylsulfonyl)imide (TFSI-) or hexafluorophosphate (PF6-) anions commonly used in lithium-ion batteries and beyond. The experimental studies revealed the temperature dependence of the diffusion coefficients for the propylene carbonate (PC) solvent and for the anions following an Arrhenius type of behavior. It was observed that the PC molecules are the faster species. For the monovalent cations (Li+, Na+, K+), the PC solvent diffusion was enhanced as the cation size increased, while for the divalent cations (Mg2+, Ca2+, Sr2+, Ba2+), the opposite trend was observed, i.e., the diffusion coefficients decreased as the cation size increased. The anion diffusion in LiTFSI and NaTFSI solutions was found to be similar, while in electrolytes with divalent cations, a decrease in anion diffusion with increasing cation size was observed. It was shown that non-polarizable charge-scaled force fields could correspond perfectly to the experimental values of the anion and PC solvent diffusion coefficients in salt solutions of both monovalent (Li+, Na+, K+) and divalent (Mg2+, Ca2+, Sr2+, Ba2+) cations at a range of operational temperatures. Finally, after calculating the radial distribution functions between cations, anions, and solvent molecules, the increase in the PC diffusion coefficient established with the increase in cation size for monovalent cations was clearly explained by the large hydration shell of small Li+ cations, due to their strong interaction with the PC solvent. In solutions with larger monovalent cations, such as Na+, and with a smaller solvation shell of PC, the PC diffusion is faster due to more liberated solvent molecules. In the salt solutions with divalent cations, both the anion and the PC diffusion coefficients decreased as the cation size increased due to an enhanced cation-anion coordination, which was accompanied by an increase in the amount of PC in the cation solvation shell due to the presence of anions.

Searching for proton transfer channels in respiratory complex I

We have explored a strategy to identify potential proton transfer channels using computational analysis of a protein structure based on Voronoi partitioning and applied it for the analysis of proton transfer pathways in redox-driven proton-pumping respiratory complex I. The analysis results in a network of connected voids/channels, which represent the dual structure of the protein; we then hydrated the identified channels using our water placement program Dowser++. Many theoretical water molecules found in the channels perfectly match the observed experimental water molecules in the structure; some other predicted water molecules have not been resolved in the experiments. The channels are of varying cross sections. Some channels are big enough to accommodate water molecules that are suitable to conduct protons; others are too narrow to hold water but require only minor conformational changes to accommodate proton transfer. We provide a preliminary analysis of the proton conductivity of the network channels, classifying the proton transfer channels as open, closed, and partially open, and discuss possible conformational changes that can modulate, i.e., open and close, the channels.

The temperature of maximum in density of aqueous solutions of nitrate and ammonium salts: Testing the Madrid-2019 force field

The shift in the temperature of maximum in density (TMD) at room pressure of aqueous solutions of a set of five salts containing NO3- and/or NH4+ groups is studied both through experiments and through molecular dynamics simulations using the Madrid-2019 force field for ions and the TIP4P/2005 model for water. The experiments demonstrate the potential transferability and limitations of the Madrid-2019 force field for nitrate and ammonium ions recently developed by our group at different temperatures and add updated information to the reported datasets of TMDs for strong electrolytes. By using the Despretz law, individual ion contributions are extracted for predictive purposes from the experimental values of the shift in the TMD. Interesting findings for the behavior of the shift in the TMD in nitrate salts expose that this property might be particularly challenging for modelization approaches when dealing with polyatomic species.

Accuracy limit of non-polarizable four-point water models: TIP4P/2005 vs OPC. Should water models reproduce the experimental dielectric constant?

The last generation of four center non-polarizable models of water can be divided into two groups: those reproducing the dielectric constant of water, as OPC, and those significantly underestimating its value, as TIP4P/2005. To evaluate the global performance of OPC and TIP4P/2005, we shall follow the test proposed by Vega and Abascal in 2011 evaluating about 40 properties to fairly address this comparison. The liquid-vapor and liquid-solid equilibria are computed, as well as the heat capacities, isothermal compressibilities, surface tensions, densities of different ice polymorphs, the density maximum, equations of state at high pressures, and transport properties. General aspects of the phase diagram are considered by comparing the ratios of different temperatures (namely, the temperature of maximum density, the melting temperature of hexagonal ice, and the critical temperature). The final scores are 7.2 for TIP4P/2005 and 6.3 for OPC. The results of this work strongly suggest that we have reached the limit of what can be achieved with non-polarizable models of water and that the attempt to reproduce the experimental dielectric constant deteriorates the global performance of the water force field. The reason is that the dielectric constant depends on two surfaces (potential energy and dipole moment surfaces), whereas in the absence of an electric field, all properties can be determined simply from just one surface (the potential energy surface). The consequences of the choice of the water model in the modeling of electrolytes in water are also discussed.

In Silico View of Crowding: Biomolecular Processes to Nanomaterial Design

It is widely accepted that deciphering biomolecular structure and function requires going beyond the single-molecule or single-complex paradigm. The densely packed macromolecules, cosolutes, and metabolites in the living cell impose crowding effects on the biomolecular structure and dynamics that need to be accounted for. Molecular simulations have proven to be a powerful tool to advance the current molecular-level understanding of such a highly concentrated, complex milieu. This Mini-Review focuses on summarizing the understanding achieved so far for the effects of crowding on biomolecular processes using computational methods, along with highlighting a new direction in employing crowding as a tool for tunable nanomaterial design. The two schools of thought that form the pillars of the current understanding of crowding effects are discussed. The investigation of crowded solutions using physics-based models that encompass different time and length scales to mimic the intracellular environment are described. The limitations and challenges faced by the current models and simulation methods are addressed, highlighting the gaps to be filled for better agreement with experiments. Crowding can also act as an effective tool to modulate the structure-property-function relationships of nanomaterials, leading to the development of novel functional materials. A few recent studies, mostly experimental, have been summarized in this direction. The Mini-Review concludes with an outlook for future developments in this field in order to enable accurate mimicking of the intracellular environment using simulations and to bridge the gap between biological processes and nanomaterial design.

Effect of ethylene oxide groups on calcite wettability reversal by nonionic surfactants: An experimental and molecular dynamics simulation investigation

HYPOTHESIS

Ethoxylated nonionic surfactants are promising candidates for enhanced oil recovery (EOR) from oil-wet carbonate reservoirs due to their ability to reverse the mineral wettability. The wettability-reversal efficiency increases with the number of the ethoxy (EO) groups in the surfactant molecule.

METHODOLOGY

Contact angle measurements, scanning electron microscopy (SEM) and molecular dynamics (MD) simulations were combined to investigate the wettability reversal of an oil-wet calcite by three ethoxylated nonionic surfactants with 1, 4 and 8 EO groups, respectively, to directly probe the role of the EO groups and to uncover the molecular mechanism responsible for the wettability reversal.

FINDINGS

Both experiments and simulations consistently show a clear correlation between the number of EO groups and the wettability reversal efficiency of the surfactants, whereby the higher number of EO groups results in greater degree of wettability reversal. This is due to 1) the more hydrophilic surfactant headgroup weakening the carboxylate interactions with the surface by expanding the surface-adjacent water layer, and 2) the physically larger surfactant molecule attracting the carboxylates more strongly, thus aiding in their removal from the surface.

Building polarization into protein-inhibitor binding dynamics in rational drug design for rheumatoid arthritis

Standard force field-based simulations to accomplish structure-based evaluations of lead molecules is a powerful tool. Combining protein fragmentation into tractable sub-systems with continuum solvation method is envisaged to enable quantum mechanics-based electronic structure calculations of macromolecules in their realistic environment. This along with incorporation of many-body polarization effect in molecular dynamics simulations may augment an accurate description of electrostatics of protein-inhibitor systems for effective drug design. Rheumatoid arthritis (RA) is a complex autoimmune disorder plagued by the ceiling effect of current targeted therapies, encouraging identification of new druggable targets and corresponding drug design to tackle the refractory form of disease. In this study, polarization-inclusive force field approach has been used to model protein solvation and ligand binding for 'Mitogen-activated protein kinase' (MAP3K8), a regulatory node of notable pharmacological relevance in RA synovial biology. For MAP3K8 inhibitors belonging to different scaffold series, the calculations illustrated differential electrostatic contribution to their relative binding affinities and successfully explained examples from available structure-activity relationship studies. Results from this study exemplified i) the advantage of this approach in reliably ranking inhibitors having close nanomolar range activities for the same target; and ii) its prospective application in lead molecule identification aiding drug discovery efforts in RA.Communicated by Ramaswamy H. Sarma.

Sodium Triflate Water-in-Salt Electrolytes in Advanced Battery Applications: A First-Principles-Based Molecular Dynamics Study

Offering a compelling combination of safety and cost-effectiveness, water-in-salt (WiS) electrolytes have emerged as promising frontiers in energy storage technology. Still, there is a strong demand for research and development efforts to make these electrolytes ripe for commercialization. Here, we present a first-principles-based molecular dynamics (MD) study addressing in detail the properties of a sodium triflate WiS electrolyte for Na-ion batteries. We have developed a workflow based on a machine learning (ML) potential derived from ab initio MD simulations. As ML potentials are typically restricted to the interpolation of the data points of the training set and have hardly any predictive properties, we subsequently optimize a classical force field based on physics principles to ensure broad applicability and high performance. Performing and analyzing detailed MD simulations, we identify several very promising properties of the sodium triflate as a WiS electrolyte but also indicate some potential stability challenges associated with its use as a battery electrolyte.

Inclusion of Water Multipoles into the Implicit Solvation Framework Leads to Accuracy Gains

The current practical "workhorses" of the atomistic implicit solvation─the Poisson-Boltzmann (PB) and generalized Born (GB) models─face fundamental accuracy limitations. Here, we propose a computationally efficient implicit solvation framework, the Implicit Water Multipole GB (IWM-GB) model, that systematically incorporates the effects of multipole moments of water molecules in the first hydration shell of a solute, beyond the dipole water polarization already present at the PB/GB level. The framework explicitly accounts for coupling between polar and nonpolar contributions to the total solvation energy, which is missing from many implicit solvation models. An implementation of the framework, utilizing the GAFF force field and AM1-BCC atomic partial charges model, is parametrized and tested against the experimental hydration free energies of small molecules from the FreeSolv database. The resulting accuracy on the test set (RMSE ∼ 0.9 kcal/mol) is 12% better than that of the explicit solvation (TIP3P) treatment, which is orders of magnitude slower. We also find that the coupling between polar and nonpolar parts of the solvation free energy is essential to ensuring that several features of the IWM-GB model are physically meaningful, including the sign of the nonpolar contributions.

Effect of counterions on the structure and dynamics of water near a negatively charged surfactant: a theoretical vibrational sum frequency generation study

Charged aqueous interfaces are of paramount importance in electrochemical, biological and environmental sciences. The properties of aqueous interfaces with ionic surfactants can be influenced by the presence of counterions. Earlier experiments involving vibrational sum frequency generation (VSFG) spectroscopy of aqueous interfaces with negatively charged sodium dodecyl sulfate (Na+DS- or SDS) surfactants revealed that the hydrogen bonding strength of the interfacial water molecules follows a certain order when salts of monovalent and divalent cations are added. It is known that cations do not directly participate in hydrogen bonding with water molecules, rather they only influence the hydrogen bonded network through their electrostatic fields. In the current work, we have simulated the aqueous interfacial systems of sodium dodecyl sulfate in the presence of chloride salts of mono and divalent countercations. The electronic polarization effects on the ions are considered at a mean-field level within the electronic continuum correction model. Our calculations of the VSFG spectra show a blue shift in the presence of added countercations whose origin is traced to different relative contributions of water molecules from the solvation shells of the surfactant headgroups and the remaining water molecules in the presence of countercations. Furthermore, the cations shield the electric fields of the surfactant headgroups, which in turn influences the contributions of water molecules to the total VSFG spectrum. This shielding effect becomes more significant when divalent countercations are present. The dynamics of water molecules is found to be slower at the interface in comparison to the bulk. The interfacial depth dependence of various dynamical quantities shows that the interface is structurally and dynamically more heterogeneous at the microscopic level.

Additive-Driven Nanoscale Architecture of Solid Electrolyte Interphase Revealed by Cryogenic Transmission Electron Microscopy

In Li metal batteries (LMBs), which boast the highest theoretical capacity, the chemical structure of the solid electrolyte interphase (SEI) serves as the key component that governs the growth of reactive Li. Various types of additives have been developed for electrolyte optimization, representing one of the most effective strategies to enhance the SEI properties for stable Li plating. However, as advanced electrolyte systems become more chemically complicated, the use of additives is empirically optimized. Indeed, their role in SEI formation and the resulting cycle life of LMBs are not well-understood. In this study, we employed cryogenic transmission electron microscopy combined with Raman spectroscopy, theoretical studies including molecular dynamics (MD) simulations and density functional theory (DFT) calculations, and electrochemical measurements to explore the nanoscale architecture of SEI modified by the most representative additives, lithium nitrate (LiNO3) and vinylene carbonate (VC), applied in a localized high-concentration electrolyte. We found that LiNO3 and VC play distinct roles in forming the SEI, governing the solvation structure, and influencing the kinetics of electrochemical reduction. Their collaboration leads to the desired SEI, ensuring prolonged cycle performance for LMBs. Moreover, we propose mechanisms for different Li growth and cycling behaviors that are determined by the physicochemical properties of SEI, such as uniformity, elasticity, and ionic conductivity. Our findings provide critical insights into the appropriate use of additives, particularly regarding their chemical compatibility.

A locally solvent-tethered polymer electrolyte for long-life lithium metal batteries

Solid polymer electrolytes exhibit enhanced Li+ conductivity when plasticized with highly dielectric solvents such as N,N-dimethylformamide (DMF). However, the application of DMF-containing electrolytes in solid-state batteries is hindered by poor cycle life caused by continuous DMF degradation at the anode surface and the resulting unstable solid-electrolyte interphase. Here we report a composite polymer electrolyte with a rationally designed Hofmann-DMF coordination complex to address this issue. DMF is engineered on Hofmann frameworks as tethered ligands to construct a locally DMF-rich interface which promotes Li+ conduction through a ligand-assisted transport mechanism. A high ionic conductivity of 6.5 × 10-4 S cm-1 is achieved at room temperature. We demonstrate that the composite electrolyte effectively reduces the free shuttling and subsequent decomposition of DMF. The locally solvent-tethered electrolyte cycles stably for over 6000 h at 0.1 mA cm-2 in Li | |Li symmetric cell. When paired with sulfurized polyacrylonitrile cathodes, the full cell exhibits a prolonged cycle life of 1000 cycles at 1 C. This work will facilitate the development of practical polymer-based electrolytes with high ionic conductivity and long cycle life.

Coordination Behavior of a Confined Ionic Liquid in Carbon Nanotubes from Molecular Dynamics Simulations

To understand the behavior of ionic liquids (ILs) at carbon material, i.e., carbon nanotube (CNT)-containing pores, we simulated different systems and analyzed their structural─in particular their coordination─behavior as well as their velocity distribution. The extension of our analysis tool CONAN presented here allowed us to study the coordination behavior as a function of the distance to the carbon material. Our systems were composed of three different CNTs combined with either the neat IL 1-ethyl-3-methylimidazolium tetrafluoroborate or with their NaBF4 salt mixtures. We investigated the impact of the force field charge scaling. As previously detected, the neat IL assumed radial layers within the confinement, with the radial density distribution depending strongly on the pore size. For the salt mixtures, the sodium cation remained in the bulk and was observed only once inside a tube. In all systems, the ions showed an overall decreased coordination behavior for regions in the bulk phase close to the carbon pore and within the confinement. The coordination number was always reduced with scaled charges. For charge scaling, higher dynamics was observed also in confinement. Interestingly, the average velocity of the atoms near the surface inside the confined space was higher than that in the center of the pore.

Tensorial Properties via the Neuroevolution Potential Framework: Fast Simulation of Infrared and Raman Spectra

Infrared and Raman spectroscopy are widely used for the characterization of gases, liquids, and solids, as the spectra contain a wealth of information concerning, in particular, the dynamics of these systems. Atomic scale simulations can be used to predict such spectra but are often severely limited due to high computational cost or the need for strong approximations that limit the application range and reliability. Here, we introduce a machine learning (ML) accelerated approach that addresses these shortcomings and provides a significant performance boost in terms of data and computational efficiency compared with earlier ML schemes. To this end, we generalize the neuroevolution potential approach to enable the prediction of rank one and two tensors to obtain the tensorial neuroevolution potential (TNEP) scheme. We apply the resulting framework to construct models for the dipole moment, polarizability, and susceptibility of molecules, liquids, and solids and show that our approach compares favorably with several ML models from the literature with respect to accuracy and computational efficiency. Finally, we demonstrate the application of the TNEP approach to the prediction of infrared and Raman spectra of liquid water, a molecule (PTAF-), and a prototypical perovskite with strong anharmonicity (BaZrO3). The TNEP approach is implemented in the free and open source software package gpumd, which makes this methodology readily available to the scientific community.

Water Potential from Adaptive Force Matching for Ice and Liquid with Revised Dispersion Predicts Supercooled Liquid Anomalies in Good Agreement with Two Independent Experimental Fits

A revised version of the Water potential from Adaptive force matching for Ice and Liquid (WAIL) was developed by using the previous data set for fitting the WAIL model but with a dispersion term calculated using symmetry adapted perturbation theory (SAPT). The model has no adjustable parameters and relies solely on fitting first-principles information. The new model, named revised WAIL (rWAIL), shows improved predictions of most properties of water when compared to the previously published WAIL model. The rWAIL model also compares favorably to other first-principles-derived water models, such as MB-Pol, at only a fraction of the computational cost. The rWAIL model is used to study the properties of supercooled water. The model shows evidence of a liquid-liquid phase transition (LLPT) in the supercooled regimes with the liquid-liquid critical point (LLCP) at 203 K and 90 MPa. This estimate is in good agreement with a recent polynomial fit to the experimental density of water. Also, the fit to the surface tension of supercooled water based on the rWAIL model shows excellent agreement with the corresponding fit to the experimental data. Consistent with previously published molecular dynamics and experimental data, the surface tension of water exhibits exponential growth in the supercooled regime, which is likely a result of the emergence of a low-density liquid form of water. The simulation thus unites two separate experimental fits with one first-principles-based model, lending strong evidence of an LLPT in real water.

Force Fields for High Concentration Aqueous KOH Solutions and Zincate Ions

Motivated by increasing interest in electrochemical devices that include highly alkaline electrolytes, we investigated two force fields for potassium hydroxide (KOH) at high concentrations in water. The "FNB" model uses the SPC/E water model, while the "FHM" model uses the TIP4P/2005 water model. We also developed parameters to describe zincate ions in these solutions. The density and viscosity of KOH using the FHM model are in better agreement with experiment than the values from the FNB model. Comparing the properties of the zincate solutions to the available experimental data, we find that both force fields agree reasonably well, although the FHM parameters give a better prediction of the viscosity. The developed force field parameters can be used in future simulations of zincate/KOH solutions in combination with other species of interest.

Adsorption Isotherm and Mechanism of Ca2+ Binding to Polyelectrolyte

Polyelectrolytes, such as poly(acrylic acid) (PAA), can effectively mitigate CaCO3 scale formation. Despite their success as antiscalants, the underlying mechanism of binding of Ca2+ to polyelectrolyte chains remains unresolved. Through all-atom molecular dynamics simulations, we constructed an adsorption isotherm of Ca2+ binding to sodium polyacrylate (NaPAA) and investigated the associated binding mechanism. We find that the number of calcium ions adsorbed [Ca2+]ads to the polymer saturates at moderately high concentrations of free calcium ions [Ca2+]aq in the solution. This saturation value is intricately connected with the binding modes accessible to Ca2+ ions when they bind to the polyelectrolyte chain. We identify two dominant binding modes: the first involves binding to at most two carboxylate oxygens on a polyacrylate chain, and the second, termed the high binding mode, involves binding to four or more carboxylate oxygens. As the concentration of free calcium ions [Ca2+]aq increases from low to moderate levels, the polyelectrolyte chain undergoes a conformational transition from an extended coil to a hairpin-like structure, enhancing the accessibility to the high binding mode. At moderate concentrations of [Ca2+]aq, the high binding mode accounts for at least one-third of all binding events. The chain's conformational change and its consequent access to the high binding mode are found to increase the overall Ca2+ ion binding capacity of the polyelectrolyte chain.

Charge of a transmembrane peptide alters its interaction with lipid membranes

Transmembrane (TM) proteins interact closely with the surrounding membrane lipids. Lipids in the vicinity of TM proteins were reported to have hindered mobility, which has been associated with lipids being caught up in the rough surface of the TM domains. These reports, however, neglect one important factor that largely influences the membrane behavior - electrostatics of the TM peptides that are usually positively charged at their cytosolic end. Here, we study on the example of a neutral and a positively charged WALP peptide, how the charge of a TM peptide influences the membrane. We investigate both its dynamics and mechanics by: (i) time dependent fluorescent shift in combination with classical and FRET generalized polarization to evaluate the mobility of lipids at short and long-range distance from the peptide, (ii) atomic force microscopy to observe the mechanical stability of the peptide-containing membranes, and (iii) molecular dynamics simulations to analyze the peptide-lipid interactions. We show that both TM peptides lower lipid mobility in their closest surroundings. The peptides cause lateral heterogeneity in lipid mobility, which in turn prevents free lipid rearrangement and lowers the membrane ability to seal ruptures after mechanical indentations. Introduction of a positive charge to the peptide largely enhances these effects, affecting the whole membrane. We thus highlight that unspecific peptide-lipid interactions, especially the electrostatics, should not be overlooked as they have a great impact on the mechanics and dynamics of the whole membrane.

Spontaneous Appearance of Triiodide Covering the Topmost Layer of the Iodide Solution Interface Without Photo-Oxidation

Ions containing iodine atoms at the vapor-aqueous solution interfaces critically affect aerosol growth and atmospheric chemistry due to their complex chemical nature and multivalency. While the surface propensity of iodide ions has been intensely discussed in the context of the Hofmeister series, the stability of various ions containing iodine atoms at the vapor-water interface has been debated. Here, we combine surface-specific sum-frequency generation (SFG) vibrational spectroscopy with ab initio molecular dynamics simulations to examine the extent to which iodide ions cover the aqueous surface. The SFG probe of the free O-D stretch mode of heavy water indicates that the free O-D group density decreases drastically at the interface when the bulk NaI concentration exceeds ∼2 M. The decrease in the free O-D group density is attributed to the spontaneous appearance of triiodide that covers the topmost interface rather than to the surface adsorption of iodide. This finding demonstrates that iodide is not surface-active, yet the highly surface-active triiodide is generated spontaneously at the water-air interface, even under dark and oxygen-free conditions. Our study provides an important first step toward clarifying iodine chemistry and pathways for aerosol formation.

Effect of the Ion, Solvent, and Thermal Interaction Coefficients on Battery Voltage

In order to increase the adoption of batteries for sustainable transport and energy storage, improved charging and discharging capabilities of lithium-ion batteries are necessary. To achieve this, accurate data that describe the internal state of the cells are essential. Several models have been derived, and transport coefficients have been reported for use in these models. We report for the first time a complete set of transport coefficients to model the concentration and temperature polarization in a lithium-ion battery ternary electrolyte, allowing us to test common assumptions. We include effects due to gradients in chemical potentials and in temperature. We find that the voltage contributions due to salt and solvent polarization are of the same order of magnitude as the ohmic loss and must be taken into account for more accurate modeling and understanding of battery performance. We report new Soret and Seebeck coefficients and find thermal polarization to be significant in cases relevant to battery research. The analysis is suitable for electrochemical systems, in general.

Ion transport mechanisms in pectin-containing EC-LiTFSI electrolytes

Using all-atom molecular dynamics simulations, we report the structure and ion transport characteristics of a new class of solid polymer electrolytes that contain the biodegradable and mechanically stable biopolymer pectin. We used highly conducting ethylene carbonate (EC) as a solvent for simulating lithium-trifluoromethanesulfonimide (LiTFSI) salt containing different weight percentages of pectin. Our simulations reveal that the pectin chains reduce the coordination number of lithium ions around their counterions (and vice versa) because of stronger lithium-pectin interactions compared to lithium-TFSI interactions. Furthermore, the pectin is found to promote smaller ionic aggregates over larger ones, in contrast to the results typically reported for liquid and polymer electrolytes. We observed that the loading of pectin in EC-LiTFSI electrolytes increases their viscosity (η) and relaxation timescales (τc), indicating higher mechanical stability, and, consequently, a decrease of the mean squared displacement, diffusion coefficient (D), and Nernst-Einstein conductivity (σNE). Interestingly, while the lithium diffusivities are related to the ion-pair relaxation timescales as D+ ∼ τc-3.1, the TFSI- diffusivities exhibit excellent correlations with ion-pair relaxation timescales as D- ∼ τc-0.95. On the other hand, the NE conductivities are dictated by distinct transport mechanisms and scales with ion-pair relaxation timescales as σNE ∼ τc-1.85.

Hyaluronan-arginine enhanced and dynamic interaction emerges from distinctive molecular signature due to electrostatics and side-chain specificity

Hyaluronan is a natural carbohydrate polymer with a negative charge that fosters gel-like conditions crucial for its cellular functions and industrial applications. As a recognized ligand for proteins, understanding their mutual interactions provides solid ground to tune hyaluronan's gel properties using biocompatible peptides. This work employs NMR and molecular dynamics simulations to identify molecular motifs relevant to hyaluronan-peptide interactions using arginine, lysine, and glycine oligopeptides. Arginine-rich peptides exhibit the strongest binding to hyaluronan according to chemical shift perturbation measurements, followed distantly by the similarly charged lysine. This difference highlights the significance of electrostatics and the peculiarities of the guanidinium side chain in arginine, capable of non-polar interactions that further stabilize the binding. Additional nuclear Overhauser effect measurements do not show stable interaction partners, precluding strong and well-defined complexes. Finally, molecular simulations support our findings and show an extended but significant interaction region, especially for arginine, responsible for the observed enhanced binding, which can also promote cross-linking of hyaluronan polymers. Our findings pave the way for optimizing biocompatible peptides to alter hyaluronan gels' properties efficiently and also explain why hyaluronan-protein interaction typically involves positively charged arginine-rich regions also capable of forming hydrogen bonds and non-polar interactions.

Modelling structure and ionic diffusion in a class of ionic liquid crystal-based solid electrolytes

Next-generation high-efficiency Li-ion batteries require an electrolyte that is both safe and thermally stable. A possible choice for high performance all-solid-state Li-ion batteries is a liquid crystal, which possesses properties in-between crystalline solids and isotropic liquids. By employing molecular dynamics simulations together with various experimental techniques, we have designed and analyzed a novel liquid crystal electrolyte composed of rigid naphthalene-based moieties as mesogenic units, grafted to flexible alkyl chains of different lengths. We have synthesized novel highly ordered lamellar phase liquid crystal electrolytes at 99% purity and have evaluated the effect of alkyl chain length variation on ionic conduction. We find that the conductivity of the liquid crystal electrolytes is directly dependent on the extent of the nanochannels formed by molecule self-organization, which itself depends non-monotonously on the size of the alkyl chains. In addition, we show that the ion pair interaction between the anionic center of the liquid crystal molecules and the Li+ ions plays a crucial role in the overall conductivity. Based on our results, we suggest that further improvement of the ionic conductivity performance is possible, making this novel family of liquid crystal electrolytes a promising option for the design of entirely solid-state Li+ ion batteries.

Hydration of biologically relevant tetramethylammonium cation by neutron scattering and molecular dynamics

Neutron scattering and molecular dynamics studies were performed on a concentrated aqueous tetramethylammonium (TMA) chloride solution to gain insight into the hydration shell structure of TMA, which is relevant for understanding its behavior in biological contexts of, e.g., properties of phospholipid membrane headgroups or interactions between DNA and histones. Specifically, neutron diffraction with isotopic substitution experiments were performed on TMA and water hydrogens to extract the specific correlation between hydrogens in TMA (HTMA) and hydrogens in water (HW). Classical molecular dynamics simulations were performed to help interpret the experimental neutron scattering data. Comparison of the hydration structure and simulated neutron signals obtained with various force field flavors (e.g. overall charge, charge distribution, polarity of the CH bonds and geometry) allowed us to gain insight into how sensitive the TMA hydration structure is to such changes and how much the neutron signal can capture them. We show that certain aspects of the hydration, such as the correlation of the hydrogen on TMA to hydrogen on water, showed little dependence on the force field. In contrast, other correlations, such as the ion-ion interactions, showed more marked changes. Strikingly, the neutron scattering signal cannot discriminate between different hydration patterns. Finally, ab initio molecular dynamics was used to examine the three-dimensional hydration structure and thus to benchmark force field simulations. Overall, while neutron scattering has been previously successfully used to improve force fields, in the particular case of TMA we show that it has only limited value to fully determine the hydration structure, with other techniques such as ab initio MD being of a significant help.

Engineering the Dielectric Constants of Polymers: From Molecular to Mesoscopic Scales

Polymers are essential components of modern-day materials and are widely used in various fields. The dielectric constant, a key physical parameter, plays a fundamental role in the light-, electricity-, and magnetism-related applications of polymers, such as dielectric and electrical insulation, battery and photovoltaic fabrication, sensing and electrical contact, and signal transmission and communication. Over the past few decades, numerous efforts have been devoted to engineering the intrinsic dielectric constant of polymers, particularly by tailoring the induced and orientational polarization modes and ferroelectric domain engineering. Investigations into these methods have guided the rational design and on-demand preparation of polymers with desired dielectric constants. This review article exhaustively summarizes the dielectric constant engineering of polymers from molecular to mesoscopic scales, with emphasis on application-driven design and on-demand polymer synthesis rooted in polymer chemistry principles. Additionally, it explores the key polymer applications that can benefit from dielectric constant regulation and outlines the future prospects of this field.

Osmotic Force Balance Evaluation of Aqueous Electrolyte Osmotic Pressures and Chemical Potentials

Concentrated aqueous salt solutions are ubiquitous in problems of biological and environmental relevance. The development of accurate force fields that capture the interactions between dissolved species in solution is crucial to simulating these systems to gain molecular insights into the underlying processes under saline conditions. The osmotic pressure is a relatively simple thermodynamic property connecting the experimental and simulation measurements of the associative properties of the ions in solution. Milner [C. Gillespie and S. T. Milner, Soft Matter, 16, 9816 (2020)] proposed a simulation approach to evaluate the osmotic pressures of salts in solution by applying a restraint potential to the ions alone in solution and determining the resulting pressure required to balance that potential, referred to here as the osmotic force balance. Here, we expand Milner's approach, demonstrating that the chemical potentials of the salts in solution as a function of concentration can be fitted to the concentration profiles determined from simulation, additionally providing an analytical expression for the osmotic pressure. This approach is used to determine the osmotic pressures of 15 alkali halide salts in water from simulations. The cross interactions between cations and anions in solution are subsequently optimized to capture their experimental osmotic pressures.

Molecular Dynamics Simulation on the Charge Transport Properties in a Salt-in-Ionic Liquid Electrolyte

A clear picture of charge transport properties in salt-in-ionic liquid electrolyte (SILE) is indispensable for the applications in lithium-ion batteries. In this study, we applied molecular dynamics (MD) simulations on a typical SILE system, composed of lithium bis(fluorosulfonyl)imide (LiFSI) with a molar fraction of 0.3 doped in 1-ethyl-3-methylimidazolium bis(fluorosulfonyl)imide (EMIMFSI). Based on the MD simulations, we calculated conductivity spectra from 108 Hz to 1014 Hz, charge current correlation functions, and charge mean square displacements, based on the center-of-mass (COM) velocities of the ions. The conductivity spectra show a bimodal feature between 1012 Hz and 1013 Hz, attributed to the interionic vibrations of the EMIM+-FSI- and Li+-FSI- contact ion pairs, respectively. Structural relaxation is observed between 109 Hz and 1012 Hz, and a flat plateau below 109 Hz, attributed to the direct current (DC) conductivity. For this SILE composed of three constituent ions, i.e., Li+, EMIM+, and FSI-, the above transport properties are further partitioned to the contributions of the individual constituent ions, including self, distinct contribution of the same constituent ions, and also the cross correlation between them. Detailed analyses on the individual contributions reveal strongly correlated motions in this complex ionic system.

In silico discovery of food-derived phytochemicals against asialoglycoprotein receptor 1 for treatment of hypercholesterolemia: Pharmacophore modeling, molecular docking and molecular dynamics simulation approach

Hypercholesterolemia is a significant risk factor for atherosclerotic cardiovascular disease (ASCVD). Successful management of cholesterol metabolism disorders can prevent these ASCVD effectively. Asialoglycoprotein receptor 1 (ASGR1) is the main subtype of sialoglycoprotein receptor, which is specifically expressed in the liver and mediates the endocytosis of blood asialoglycoprotein to lysosome degradation. Recently, ASGR1 has been reported as a new therapeutic target for the treatment of hypercholesterolemia. In this study, the main aim was to identify natural ASGR1 inhibitors from plant food chemicals library through pharmacophore and docking based virtual screening. Total 14 phytochemicals of potential ASGR1 inhibitors were identified, which presented docking affinity higher than control compound through docking based virtual screening. The docking pose showed the top three hits interacted residues were located at active pocket of ASGR1 with hydrogen bonds, hydrophobic interactions and electrostatic interactions. The top three hits (ZINC85664954, ZINC169372863, and ZINC195764535) were then subjected to 200 ns molecular dynamics simulation to evaluate the stability of docked complexes. These results showed that selected phytochemicals bound to ASGR1 with higher stability than control compound. Binding free energy of each docked complex was calculated by the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) method. The binding free energy of ZINC85664954, ZINC169372863, ZINC195764535, and control-ASGR1 docked complexes were -18.359, -13.303, -14.389, and -6.229 kcal/mol, respectively. This indicated that selected hits bound to ASGR1 with higher affinity than control compound. Network pharmacology analysis shows that these phytochemicals have obvious multiple-effects and can regulate various biochemical pathways related to hypercholesterolemia. Besides, selected phytochemicals have suitable pharmacokinetics properties, suggesting that these compounds may be potential drug candidates. This study may be contributed to rational design of novel ASGR1 inhibitors for treatment of hypercholesterolemia.

Novel Computational Chemistry Infrastructure for Simulating Astatide in Water: From Basis Sets to Force Fields Using Particle Swarm Optimization

Using the example of astatine, the heaviest naturally occurring halogen whose isotope At-211 has promising medical applications, we propose a new infrastructure for large-scale computational models of heavy elements with strong relativistic effects. In particular, we focus on developing an accurate force field for At- in water based on reliable relativistic density functional theory (DFT) calculations. To ensure the reliability of such calculations, we design novel basis sets for relativistic DFT, via the particle swarm optimization algorithm to optimize the coefficients of the new basis sets and the polarization-consistent basis set idea's extension to heavy elements to eliminate the basis set error from DFT calculations. The resulting basis sets enable the well-grounded evaluation of relativistic DFT against "gold-standard" CCSD(T) results. Accounting for strong relativistic effects, including spin-orbit interaction, via our redesigned infrastructure, we elucidate a noticeable dissimilarity between At- and I- in halide-water force field parameters, radial distribution functions, diffusion coefficients, and hydration energies. This work establishes the framework for the systematic development of polarization-consistent basis sets for relativistic DFT and accurate force fields for molecular dynamics simulations to be used in large-scale models of complex molecular systems with elements from the bottom of the periodic table, including actinides and even superheavy elements.

In-Cell Dynamics: The Next Focus of All-Atom Simulations

The cell is a crowded space where large biomolecules and metabolites are in continuous motion. Great strides have been made in in vitro studies of protein dynamics, folding, and protein-protein interactions, and much new data are emerging of how they differ in the cell. In this Perspective, we highlight the current progress in atomistic modeling of in-cell environments, both bacteria and mammals, with emphasis on classical all-atom molecular dynamics simulations. These simulations have been recently used to capture and characterize functional and non-functional protein-protein interactions, protein folding dynamics of small proteins with varied topologies, and dynamics of metabolites. We further discuss the challenges and efforts for updating modern force fields critical to the progress of cellular environment simulations. We also briefly summarize developments in relevant state-of-the-art experimental techniques. As computational and experimental methodologies continue to progress and produce more directly comparable data, we are poised to capture the complex atomistic picture of the cell.

Effects of ion type and concentration on the structure and aggregation of the amyloid peptide A β 16 - 22 $$ {\boldsymbol{\beta}}_{16-22} $$

Among the various factors controlling the amyloid aggregation process, the influences of ions on the aggregation rate and the resulting structures are important aspects to consider, which can be studied by molecular simulations. There is a wide variety of protein force fields and ion models, raising the question of which model to use in such studies. To address this question, we perform molecular dynamics simulations of Aβ16-22 , a fragment of the Alzheimer's amyloid β peptide, using different protein force fields, AMBER99SB-disp (A99-d) and CHARMM36m (C36m), and different ion parameters. The influences of NaCl and CaCl2 at various concentrations are studied and compared with the systems without the addition of ions. Our results indicate a sensitivity of the peptide-ion interactions to the different ion models. In particular, we observe a strong binding of Ca2+ to residue E22 with C36m and also with the Åqvist ion model used together with A99-d, which slightly affects the monomeric Aβ16-22 structures and the aggregation rate, but significantly affects the oligomer structures formed in the aggregation simulations. For example, at high Ca2+ concentrations, there was a switch from an antiparallel to a parallel β-sheet. Such ionic influences are of biological relevance because local ion concentrations can change in vivo and could help explain the polymorphism of amyloid fibrils.

Binding of carboxylate and water to monovalent cations

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. To gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li+, Na+ and K+. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

Why do polyarginines adsorb at neutral phospholipid bilayers and polylysines do not? An insight from density functional theory calculations and molecular dynamics simulations

Adsorption of cell-penetrating peptides (CPPs) at cellular membranes is the first and necessary step for their subsequent translocation across cellular membranes into the cytosol. It has been experimentally shown that CPPs rich in arginine (Arg) amino acid penetrate across phospholipid bilayers more effectively than their lysine (Lys) rich counterparts. In this work, we aim to understand the differences in the first translocation step, adsorption of Arg9 and Lys9 peptides at fully hydrated neutral phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid bilayers and evaluate in detail the energetics of the process using molecular dynamics (MD) simulations and free energy calculations of adsorption of the single peptide. We show that the adsorption of Arg9 is energetically feasible, with the free energy of adsorption being ∼-5.0 kcal mol-1 at PC and ∼-5.5 kcal mol-1 at PE bilayers. In contrast, adsorption of Lys9 is not observed at PC bilayers, and their adsorption at PE bilayers is very weak, being ∼-0.5 kcal mol-1. We show by energy decomposition and analysis of peptide hydration along the membrane that significantly stronger electrostatic interactions of Arg9 with lipid phosphate groups, together with the greater loss of peptide hydration (and in turn stronger hydrophobic interactions) along the membrane translocation path, are the main driving factors governing the adsorption of Arg-rich peptides at neutral lipid bilayers in contrast to Lys-rich peptides. Finally, we also compare the energetics in lipid/bilayer systems with the density functional theory (DFT) calculations of the corresponding model systems in the continuum water model and reveal the energetic differences in different environments.

Evolution of the electrical double layer with electrolyte concentration probed by second harmonic scattering

Investigating the electrical double layer (EDL) structure has been a long-standing challenge and has seen the emergence of several sophisticated techniques able to probe selectively the few molecular layers of a solid/water interface. While a qualitative estimation of the thickness of the EDL can be obtained using simple theoretical models, following experimentally its evolution is not straightforward and can be even more complicated in nano- or microscale systems, particularly when changing the ionic concentration by several orders of magnitude. Here, we bring insight into the structure of the EDL of SiO2 nanoparticle suspensions and its evolution with increasing ionic concentration using angle-resolved second harmonic scattering (AR-SHS). Below millimolar salt concentrations, we can successively characterize inner-sphere adsorption, diffuse layer formation, and outer-sphere adsorption. Moreover, we show for the first time that, by appropriately selecting the nanoparticle size, it is possible to retrieve information also in the millimolar range. There, we observe a decrease in the magnitude of the surface potential corresponding to a compression in the EDL thickness, which agrees with the results of several other electroanalytical and optical techniques. Molecular dynamics simulations suggest that the EDL compression mainly results from the diffuse layer compression rather than outer-sphere ions (Stern plane) moving closer to the surface.

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

We develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

Nonaqueous Ion Pairing Exemplifies the Case for Including Electronic Polarization in Molecular Dynamics Simulations

The inclusion of electronic polarization is of crucial importance in molecular simulations of systems containing charged moieties. When neglected, as often done in force field simulations, charge-charge interactions in solution may become severely overestimated, leading to unrealistically strong bindings of ions to biomolecules. The electronic continuum correction introduces electronic polarization in a mean-field way via scaling of charges by the reciprocal of the square root of the high-frequency dielectric constant of the solvent environment. Here, we use ab initio molecular dynamics simulations to quantify the effect of electronic polarization on pairs of like-charged ions in a model nonaqueous environment where electronic polarization is the only dielectric response. Our findings confirm the conceptual validity of this approach, underlining its applicability to complex aqueous biomolecular systems. Simultaneously, the results presented here justify the potential employment of weaker charge scaling factors in force field development.

In silico discovery of potential sodium-glucose cotransporter-2 inhibitors from natural products for treatment of heart failure via molecular docking and molecular dynamics simulation approach

Heart failure (HF) is the end stage of cardiovascular disease. Because of its complex condition and poor prognosis, HF has become an important public health problem in the world. Sodium-glucose cotransporter-2 (SGLT2) is a member of the glucose transporter family. Recently, SGLT2 inhibitors have been applied to treat HF. In this study, the main aim was to discover natural SGLT2 inhibitor from Chinese herbs through docking-based virtual screening. Totally 113 natural compounds of potential SGLT2 inhibitor were identified, which displayed docking affinity higher than six approved inhibitors (dapagliflozin (IC50 = 4.9 nM), canagliflozin (IC50 = 4.4 nM 6.7), ipragliflozin (IC50 = 7.4 nM), empagliflozin (IC50 = 3.1 nM), tofogliflozin (IC50 = 4 nM) and luseogliflozin (IC50 = 2.3 nM)) through docking-based virtual screening. Then, the top three hits (ZINC70455591, ZINC85594065 and ZINC14588133) and six known inhibitors were selected for molecular dynamics simulation and the binding free energy calculation using molecular mechanics Poisson-Boltzmann surface area to demonstrate the stability and affinity of docked complexes. These results pointed out that the three docked complexes were stabilized and the chosen compounds were tightly adhering to the binding site of SGLT2. Besides, pharmacokinetic properties of the selected compounds showed those natural compounds may be potential drug candidates. This study may be contributed to further in vitro and in vivo validation and the development of novel SGLT2 inhibitor for treating HF.Communicated by Ramaswamy H. Sarma.

Accurate Calculation of Solvation Properties of Lithium Ions in Nonaqueous Solutions

We perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

Molecular Modeling of Water-in-Salt Electrolytes: A Comprehensive Analysis of Polarization Effects and Force Field Parameters in Molecular Dynamics Simulations

Accurate modeling of highly concentrated aqueous solutions, such as water-in-salt (WiS) electrolytes in battery applications, requires proper consideration of polarization contributions to atomic interactions. Within the force field molecular dynamics (MD) simulations, the atomic polarization can be accounted for at various levels. Nonpolarizable force fields implicitly account for polarization effects by incorporating them into their van der Waals interaction parameters. They can additionally mimic electron polarization within a mean-field approximation through ionic charge scaling. Alternatively, explicit polarization description methods, such as the Drude oscillator model, can be selectively applied to either a subset of polarizable atoms or all polarizable atoms to enhance simulation accuracy. The trade-off between simulation accuracy and computational efficiency highlights the importance of determining an optimal level of accounting for atomic polarization. In this study, we analyze different approaches to include polarization effects in MD simulations of WiS electrolytes, with an example of a Na-OTF solution. These approaches range from a nonpolarizable to a fully polarizable force field. After careful examination of computational costs, simulation stability, and feasibility of controlling the electrolyte properties, we identify an efficient combination of force fields: the Drude polarizable force field for salt ions and non-polarizable models for water. This cost-effective combination is sufficiently flexible to reproduce a broad range of electrolyte properties, while ensuring simulation stability over a relatively wide range of force field parameters. Furthermore, we conduct a thorough evaluation of the influence of various force field parameters on both the simulation results and technical requirements, with the aim of establishing a general framework for force field optimization and facilitating parametrization of similar systems.

Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property

In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.

Polarizability Plays a Decisive Role in Modulating Association between Molecular Cations and Anions

Electrostatic interactions involving proteins depend on not only the ionic charges involved but also their chemical identities. Here we examine the origins of incompletely understood differences in the strength of association of different pairs of monovalent molecular ions that are relevant to protein-protein and protein-ligand interactions. Cationic analogues of the basic amino acid side chains are simulated, along with oxyanionic analogues of cation-exchange ligands and acidic amino acids. Experimentally observed association trends with respect to the cations, but not anions, are captured by a nonpolarizable model. An effective continuum correction to account for electronic polarizability can capture both trends better but at the expense of fidelity to the underlying free energy landscape for ion-pair association. A polarizable model proves decisive in capturing experimentally suggested trends with respect to both cations and anions; critically, the free energy landscape for ion-pair association is itself altered, thus altering configurational sampling.

Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions

Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.

RNA-protein complexes and force field polarizability

Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.

Surface Affinity of Tetramethylammonium Iodide in Aqueous Solutions: A Combined Experimental and Computer Simulation Study

The surface affinity of tetramethylammonium iodide (TMAI) in aqueous solutions is investigated by surface tension measurements and molecular dynamics computer simulations. Experiments, performed in the entire composition range of solubility using the pendant drop method with two different setups, clearly reveal that TMAI is a weakly capillary active salt. Computer simulations performed with the AMBER force field reproduce the experimental data very well, while two other major force fields (i.e., CHARMM and OPLS) can still reproduce the experimental trend qualitatively; however, even qualitative reproduction of the experimental trend requires scaling down the ion charges according to the Leontyev-Stuchebrukhov correction. On the other hand, the GROMOS force field fails in reproducing the experimentally confirmed capillary activity of TMAI. Molecular dynamics simulation results show that, among the two ions, iodide has a clearly larger surface affinity than tetramethylammonium (TMA⁺). Further, the adsorption of the I^- anions is strictly limited to the first molecular layer beneath the liquid-vapor interface, which is followed by several layers of their depletion. On the other hand, the net negative charge of the surface layer, caused by the excess amount of I^- with respect to TMA⁺, is compensated by a diffuse layer of adsorbed TMA⁺ cations, extending to or beyond the fourth molecular layer beneath the liquid surface.