A combined experimental and theoretical study of ion solvation in liquid N-methylacetamide

Balancing Group I Monatomic Ion-Polar Compound Interactions for Condensed Phase Simulation in the Polarizable Drude Force Field

Molecular dynamics (MD) simulations are a commonly used method for investigating molecular behavior at the atomic level. Achieving reliable MD simulation results necessitates the use of an accurate force field. In the present work, we present a protocol to enhance the quality of group 1 monatomic ions (specifically Li+, Na+, K+, Rb+, and Cs+) with respect to their interactions with common polar model compounds in biomolecules in condensed phases in the context of the Drude polarizable force field. Instead of adjusting preexisting individual parameters for ions, model compounds, and water, we employ atom-pair specific Lennard-Jones (LJ) (known as NBFIX in CHARMM) and through-space Thole dipole screening (NBTHOLE) terms to fine-tune the balance of ion-model compound, ion-water, and model compound-water interactions. This involved establishing a protocol for the optimization of NBFIX and NBTHOLE parameters targeting the difference between molecular mechanical (MM) and quantum mechanical (QM) potential energy scans (PES). It is shown that targeting PES involving complexes that include multiple model compounds and/or ions as trimers and tetramers yields parameters that produce condensed phase properties in agreement with experimental data. Validation of this protocol involved the reproduction of experimental thermodynamic benchmarks, including solvation free energies of ions in methanol and N-methylacetamide, osmotic pressures, ionic conductivities, and diffusion coefficients within the condensed phase. These results show the importance of including more complex ion-model compound complexes beyond dimers in the QM target data to account for many-body effects during parameter fitting. The presented parameters represent a significant refinement of the Drude polarizable force field, which will lead to improved accuracy for modeling ion-biomolecular interactions.

Insight into the K channel's selectivity from binding of K+, Na+ and water to N-methylacetamide

In potassium channels that conduct K+ selectively over Na+, which sites are occupied by K+ or water and the mechanism of selectivity are unresolved questions. The combination of the energetics and the constraints imposed by the protein structure yield the selective permeation and occupancy. To gain insight into the combination of structure and energetics, we performed density functional theory (DFT) calculations of multiple N-methyl acetamide (NMA) ligands binding to K+ and Na+, relative to hydrated K+ and Na+. NMA is an analogue of the amino acid backbone and provides the carbonyl binding to the ions that occurs in most binding sites of the K+ channel. Unconstrained optimal structures are obtained through geometry optimization calculations of the NMA ligand binding. The complexes formed by 8 NMA binding to the cations have the O atoms positioned in nearly identical locations as the O atoms in the selectivity filter. The transfer free energies between bulk water and K+ or Na+ bound to 8 NMA are almost identical, implying there is no selectivity by a single site. For water optimized with 8 NMA, binding is weak and O atoms are not positioned as in the K+ channel selectivity filter, suggesting that the ions are much more favored than water. Optimal structures of 8 NMA binding with two cations (K+ or Na+) are stable and have lower binding free energy than the optimal structures with just one cation. However, in the Na+ case, the optimal structure deforms and does not match the K+ channel; that is, two bound Na+ are destabilizing. In contrast, the two K+ structure is stabilized and the selectivity free energy favors K+. Overall, this study shows that binding site occupancy and the mechanism for K+ selectivity involves multiple K+ binding in multiple neighboring layers or sites of the K+ channel selectivity filter.

Interactions between selectivity filter and pore helix control filter gating in the MthK channel

K+ channel activity can be limited by C-type inactivation, which is likely initiated in part by dissociation of K+ ions from the selectivity filter and modulated by the side chains that surround it. While crystallographic and computational studies have linked inactivation to a "collapsed" selectivity filter conformation in the KcsA channel, the structural basis for selectivity filter gating in other K+ channels is less clear. Here, we combined electrophysiological recordings with molecular dynamics simulations, to study selectivity filter gating in the model potassium channel MthK and its V55E mutant (analogous to KcsA E71) in the pore-helix. We found that MthK V55E has a lower open probability than the WT channel, due to decreased stability of the open state, as well as a lower unitary conductance. Simulations account for both of these variables on the atomistic scale, showing that ion permeation in V55E is altered by two distinct orientations of the E55 side chain. In the "vertical" orientation, in which E55 forms a hydrogen bond with D64 (as in KcsA WT channels), the filter displays reduced conductance compared to MthK WT. In contrast, in the "horizontal" orientation, K+ conductance is closer to that of MthK WT; although selectivity filter stability is lowered, resulting in more frequent inactivation. Surprisingly, inactivation in MthK WT and V55E is associated with a widening of the selectivity filter, unlike what is observed for KcsA and reminisces recent structures of inactivated channels, suggesting a conserved inactivation pathway across the potassium channel family.

The Persistent Question of Potassium Channel Permeation Mechanisms

Potassium channels play critical roles in many physiological processes, providing a selective permeation route for K⁺ ions in and out of a cell, by employing a carefully designed selectivity filter, evolutionarily conserved from viruses to mammals. The structure of the selectivity filter was determined at atomic resolution by x-ray crystallography, showing a tight coordination of desolvated K⁺ ions by the channel. However, the molecular mechanism of K⁺ ions permeation through potassium channels remains unclear, with structural, functional and computational studies often providing conflicting data and interpretations. In this review, we will present the proposed mechanisms, discuss their origins, and will critically assess them against all available data. General properties shared by all potassium channels are introduced first, followed by the introduction of two main mechanisms of ion permeation: soft and direct knock-on. Then, we will discuss critical computational and experimental studies that shaped the field. We will especially focus on molecular dynamics (MD) simulations, that provided mechanistic and energetic aspects of K⁺ permeation, but at the same time created long-standing controversies. Further challenges and possible solutions are presented as well.

Polarisable force fields: what do they add in biomolecular simulations?

The quality of biomolecular simulations critically depends on the accuracy of the force field used to calculate the potential energy of the molecular configurations. Currently, most simulations employ non-polarisable force fields, which describe electrostatic interactions as the sum of Coulombic interactions between fixed atomic charges. Polarisation of these charge distributions is incorporated only in a mean-field manner. In the past decade, extensive efforts have been devoted to developing simple, efficient, and yet generally applicable polarisable force fields for biomolecular simulations. In this review, we summarise the latest developments in accounting for key biomolecular interactions with polarisable force fields and applications to address challenging biological questions. In the end, we provide an outlook for future development in polarisable force fields.

Challenges and advances in atomistic simulations of potassium and sodium ion channel gating and permeation

Ion channels are implicated in many essential physiological events such as electrical signal propagation and cellular communication. The advent of K⁺ and Na⁺ ion channel structure determination has facilitated numerous investigations of molecular determinants of their behaviour. At the same time, rapid development of computer hardware and molecular simulation methodologies has made computational studies of large biological molecules in all-atom representation tractable. The concurrent evolution of experimental structural biology with biomolecular computer modelling has yielded mechanistic details of fundamental processes unavailable through experiments alone, such as ion conduction and ion channel gating. This review is a short survey of the atomistic computational investigations of K⁺ and Na⁺ ion channels, focusing on KcsA and several voltage-gated channels from the K_V and Na_V families, which have garnered many successes and engendered several long-standing controversies regarding the nature of their structure-function relationship. We review the latest advancements and challenges facing the field of molecular modelling and simulation regarding the structural and energetic determinants of ion channel function and their agreement with experimental observations.

LiCl solvation in N-methyl-acetamide (NMA) as a model for understanding Li(+) binding to an amide plane

The thermodynamics of ion solvation in non-aqueous solvents remains of great significance for understanding cellular transport and ion homeostasis for the design of novel ion-selective materials and applications in molecular pharmacology. Molecular simulations play pivotal roles in connecting experimental measurements to the microscopic structures of liquids. One of the most useful and versatile mimetic systems for understanding biological ion transport is N-methyl-acetamide (NMA). A plethora of theoretical studies for ion solvation in NMA have appeared recently, but further progress is limited by two factors. One is an apparent lack of experimental data on solubility and thermodynamics of solvation for a broad panel of 1 : 1 salts over an appropriate temperature and concentration range. The second concern is more substantial and has to do with the limitations hardwired in the additive (fixed charge) approximations used for most of the existing force-fields. In this submission, we report on the experimental evaluation of LiCl solvation in NMA over a broad range of concentrations and temperatures and compare the results with those of MD simulations with several additive and one polarizable force-field (Drude). By comparing our simulations and experimental results to density functional theory computations, we discuss the limiting factors in existing potential functions. To evaluate the possible implications of explicit and implicit polarizability treatments on ion permeation across biological channels, we performed potential of mean force (PMF) computations for Li(+) transport through a model narrow ion channel with additive and polarizable force-fields.

Representation of Ion-Protein Interactions Using the Drude Polarizable Force-Field

Small metal ions play critical roles in numerous biological processes. Of particular interest is how metalloenzymes are allosterically regulated by the binding of specific ions. Understanding how ion binding affects these biological processes requires atomic models that accurately treat the microscopic interactions with the protein ligands. Theoretical approaches at different levels of sophistication can contribute to a deeper understanding of these systems, although computational models must strike a balance between accuracy and efficiency in order to enable long molecular dynamics simulations. In this study, we present a systematic effort to optimize the parameters of a polarizable force field based on classical Drude oscillators to accurately represent the interactions between ions (K⁺, Na⁺, Ca²⁺, and Cl^–) and coordinating amino-acid residues for a set of 30 biologically important proteins. By combining ab initio calculations and experimental thermodynamic data, we derive a polarizable force field that is consistent with a wide range of properties, including the geometries and interaction energies of gas-phase ion/protein-like model compound clusters, and the experimental solvation free-energies of the cations in liquids. The resulting models display significant improvements relative to the fixed-atomic-charge additive CHARMM C36 force field, particularly in their ability to reproduce the many-body electrostatic nonadditivity effects estimated from ab initio calculations. The analysis clarifies the fundamental limitations of the pairwise additivity assumption inherent in classical fixed-charge force fields, and shows its dramatic failures in the case of Ca²⁺ binding sites. These optimized polarizable models, amenable to computationally efficient large-scale MD simulations, set a firm foundation and offer a powerful avenue to study the roles of the ions in soluble and membrane transport proteins.

Anharmonic vibrations of N-H in Cl(-)(N-methylacetamide)1(H2O)(0-2)Ar2 cluster ions. Combined IRPD experiments and BOMD simulations

Infrared Predissociation (IRPD) spectra of Cl(-)(NMA)1(H2O)0-2Ar2 combined with Born-Oppenheimer Molecular Dynamics (BOMD) IR spectra have been acquired, providing the structure and dynamics of these systems. We show that the chloride ion is bound to the hydrogen of the amide N-H group, forming a strong ionic hydrogen bond, weakening the N-H stretch, and shifting it to lower frequency. The presence of water molecules enhances the ionic hydrogen bond by binding to the amide carbonyl oxygen of NMA and shifts the N-H stretch further to lower frequency. The BOMD IR spectra can recapture all, but about 100 cm(-1), of the 600 to 700 cm(-1) shifts due to the strong N-H stretch anharmonicities observed in experiments. This residual error was found to be due to the lack of zero point energy in the classical treatment of motion in the BOMD method.

Relative Free Energies for Hydration of Monovalent Ions from QM and QM/MM Simulations

Methods directly evaluating the hydration structure and thermodynamics of physiologically relevant cations (Na(+), K(+), Cl(-), etc.) have wide ranging applications in the fields of inorganic, physical, and biological chemistry. All-atom simulations based on accurate potential energy surfaces appear to offer a viable option for assessing the chemistry of ion solvation. Although MD and free energy simulations of ion solvation with classical force fields have proven their usefulness, a number of challenges still remain. One of them is the difficulty of force field benchmarking and validation against structural and thermodynamic data obtained for a condensed phase. Hybrid quantum mechanical/molecular mechanical (QM/MM) models combined with sampling algorithms have the potential to provide an accurate solvation model and to incorporate the effects from the surrounding, which is often missing in gas-phase ab initio computations. Herein, we report the results from QM/MM free energy simulations of Na(+)/K(+) and Cl(-)/Br(-) hydration where we simultaneously characterized the relative thermodynamics of ion solvation and changes in the solvation structure. The Flexible Inner Region Ensemble Separator (FIRES) method was used to impose a spatial separation between QM region and the outer sphere of solvent molecules treated with the CHARMM27 force field. FEP calculations based on QM/MM simulations utilizing the CHARMM/deMon2k interface were performed with different basis set combinations for K(+)/Na(+) and Cl(-)/Br(-) perturbations to establish the dependence of the computed free energies on the basis set level. The dependence of the computed relative free energies on the size of the QM and MM regions is discussed. The current methodology offers an accurate description of structural and thermodynamic aspects of the hydration of alkali and halide ions in neat solvents and can be used to obtain thermodynamic data on ion solvation in condensed phase along with underlying structural properties of the ion-solvent system.

AUTOMATED FORCE FIELD PARAMETERIZATION FOR NON-POLARIZABLE AND POLARIZABLE ATOMIC MODELS BASED ON AB INITIO TARGET DATA

Classical molecular dynamics (MD) simulations based on atomistic models are increasingly used to study a wide range of biological systems. A prerequisite for meaningful results from such simulations is an accurate molecular mechanical force field. Most biomolecular simulations are currently based on the widely used AMBER and CHARMM force fields, which were parameterized and optimized to cover a small set of basic compounds corresponding to the natural amino acids and nucleic acid bases. Atomic models of additional compounds are commonly generated by analogy to the parameter set of a given force field. While this procedure yields models that are internally consistent, the accuracy of the resulting models can be limited. In this work, we propose a method, General Automated Atomic Model Parameterization (GAAMP), for generating automatically the parameters of atomic models of small molecules using the results from ab initio quantum mechanical (QM) calculations as target data. Force fields that were previously developed for a wide range of model compounds serve as initial guess, although any of the final parameter can be optimized. The electrostatic parameters (partial charges, polarizabilities and shielding) are optimized on the basis of QM electrostatic potential (ESP) and, if applicable, the interaction energies between the compound and water molecules. The soft dihedrals are automatically identified and parameterized by targeting QM dihedral scans as well as the energies of stable conformers. To validate the approach, the solvation free energy is calculated for more than 200 small molecules and MD simulations of 3 different proteins are carried out.

Charging free energy calculations using the Generalized Solvent Boundary Potential (GSBP) and periodic boundary condition: a comparative analysis using ion solvation and oxidation free energy in proteins

Free energy simulations using a finite sphere boundary condition rather than a periodic boundary condition (PBC) are attractive in the study of very large biomolecular systems. To understand the quantitative impact of various approximations in such simulations, we compare charging free energies in both solution and protein systems calculated in a linear response framework with the Generalized Solvent Boundary Potential (GSBP) and PBC simulations. For simple ions in solution, we find good agreements between GSBP and PBC charging free energies, once the relevant correction terms are taken into consideration. For PBC simulations with the particle-mesh-Ewald for long-range electrostatics, the contribution (ΔG(P-M)) due to the use of a particle rather than molecule based summation scheme in real space is found to be significant, as pointed out by Hünenberger and co-workers. For GSBP, when the inner region is close to be charge neutral, the key correction is the overpolarization of water molecules at the inner/outer dielectric boundary; the magnitude of the correction (ΔG(s-pol)), however, is relatively small. For charging (oxidation) free energy in proteins, the situation is more complex, although good agreement between GSBP and PBC can still be obtained when care is exercised. The smooth dielectric boundary approximation inherent to GSBP tends to make significant errors when the inner region is featured with a high net charge. However, the error can be corrected with Poisson-Boltzmann calculations using snapshots from GSBP simulations in a straightforward and robust manner. Because of the more complex charge and solvent distributions, the magnitudes of ΔG(P-M) and ΔG(s-pol) in protein simulations appear to be different from those derived for solution simulations, leading to uncertainty in directly comparing absolute charging free energies from PBC and GSBP simulations for protein systems. The relative charging/oxidation free energies, however, are robust. With the linear response approximation, for the specific protein system (CueR) studied, the effect of freezing the protein structure in the outer region is found to be small, unless a very small (8 Å) inner region is used; even in the latter case, the result is substantially improved when the nearby metal binding loop is allowed to respond to metal oxidation. The implications of these results to the applicability of GSBP to complex biomolecules and in ab initio QM/MM simulations are discussed.

Six-site polarizable model of water based on the classical Drude oscillator

A polarizable water model, SWM6, was developed and optimized for liquid phase simulations under ambient conditions. Building upon the previously developed SWM4-NDP model, additional sites representing oxygen lone-pairs were introduced. The geometry of the sites is assumed to be rigid. Considering the large number of adjustable parameters, simulated annealing together with polynomial fitting was used to facilitate model optimization. The new water model was shown to yield the correct self-diffusion coefficient after taking the system size effect into account, and the dimer geometry is better reproduced than in the SWM4 models. Moreover, the experimental oxygen-oxygen radial distribution is better reproduced, indicating that the new model more accurately describes the local hydrogen bonding structure of bulk phase water. This was further validated by its ability to reproduce the experimental nuclear magnetic shielding and related chemical shift of the water hydrogen in the bulk phase, a property sensitive to the local hydrogen bonding structure. In addition, comparison of the liquid properties of the SWM6 model is made with those of a number of widely used additive and polarizable models. Overall, improved balance between the description of monomer, dimer, clustered, and bulk phase water is obtained with the new model compared to its SWM4-NDP polarizable predecessor, though application of the model requires an approximately twofold increase on computational resources.

How does overcoordination create ion selectivity?

Some biological molecules can distinguish between ions of similar nature, which may be achieved by enforcing specific coordination numbers on ions in the binding site. It is suggested that when this number is favourable for one ion type, but too large for another, this creates ion selectivity through the proposed mechanism of 'overcoordination'. Much debate has occurred about the role overcoordination plays, and suggestions made as to how molecules can enforce particular coordination numbers, but there has not been an examination of the microscopic underpinning of ion selectivity by overcoordination. Here we use molecular-dynamics to systematically investigate how the number of ligands affects the ion-ligand and ligand-ligand interaction energies, and thus the thermodynamic ion selectivity, of a combination of model systems: three ions (Li(+)/Na(+)/K(+)) with three different ligands (water/formaldehyde/formamide). We find that the ligand-ligand repulsion controls the changes in geometry of each system with changing ligand number. Ion selectivity by overcoordination is achieved as smaller ions exhibit anomalous geometrical changes with the addition of extra ligands, whilst larger ions do not.

SIZE DEPENDENCE OF SOLVATION STRUCTURE AND DYNAMICS OF IONS IN LIQUID N-METHYLACETAMIDE: A MOLECULAR DYNAMICS SIMULATION STUDY

The solvation structure and dynamics of alkali metal (Li+, Na+, K+, Rb+, Cs+) and halide (F-, Cl-, Br-, I-) ions in liquid N -methylacetamide (NMA) are calculated at two different temperatures T = 313 K and 453 K, by using classical molecular dynamics simulations. We have also considered [Formula: see text] and some larger cations such as I +, Me 4 N +, and Et4N+ in this study to investigate the size dependence solvation structure and dynamics of ions in liquid NMA. With the increase of ion size, the self-diffusion coefficients of cations are found to increase and the maximum is observed at Me4N+ , whereas for halide ions the increase of diffusion coefficient with ion size continues up to I- and no maximum is observed. However, the relative increase of the diffusion coefficients of larger ion compared to those of Li+ and F+ are found to be significantly higher at low temperature. Results are very good in agreement with experimental observation.

Molecular Mechanisms of K+ Selectivity in Na/K Pump

The sodium–potassium (Na/K) pump plays an essential role in maintaining cell volume and secondary active transport of other solutes by establishing the Na+ and K+ concentration gradients across the plasma membrane of animal cells. The recently determined crystal structures of the Na/K pump to atomic resolution provide a new impetus to investigate molecular determinants governing the binding of Na+ and K+ ions and conformational transitions during the functional cycle. The pump cycle is generally described by the alternating access mechanism, in which the pump toggles between different conformational states, where ions can bind from either the intracellular or the extracellular side. However, important issues concerning the selectivity of the Na/K pump remain to be addressed. In particular, two out of the three binding sites are shared between Na+ and K+ and it is not clear how the protein is able to select K+ over Na+ when it is in the outwardly facing phosphorylated conformation (E2P), and Na+ over K+ when it is in the inwardly facing conformation (E1). In this review article, we will first briefly review the recent advancement in understanding the microscopic mechanism of K+ selectivity in the Na/K pump at the E2·Pi state and then outline the remaining challenges to be addressed about ion selectivity.

Specific interactions of sodium salts with alanine dipeptide and tetrapeptide in water: insights from molecular dynamics

We examine computationally the dipeptide and tetrapeptide of alanine in pure water and solutions of sodium chloride (NaCl) and iodide (NaI), with salt concentrations up to 3 M. Enhanced sampling of the configuration space is achieved by the replica exchange method. In agreement with other works, we observe preferential sodium interactions with the peptide carbonyl groups, which are enhanced in the NaI solutions due to the increased affinity of the less hydrophilic iodide anion for the peptide methyl side-chains and terminal blocking groups. These interactions have been associated with a decrease in the helicities of more complex peptides. In our simulations, both salts have a small effect on the dipeptide, but consistently stabilize the intramolecular hydrogen-bonding interactions and "α-helical" conformations of the tetrapeptide. This behavior, and an analysis of the intermolecular interaction energies show that ion-peptide interactions, or changes in the peptide hydration due to salts, are not sufficient determining factors of the peptide conformational preferences. Additional simulations suggest that the observed stabilizing effect is not due to the employed force-field, and that it is maintained in short peptides but is reversed in longer peptides. Thus, the peptide conformational preferences are determined by an interplay of energetic and entropic factors, arising from the peptide sequence and length and the composition of the solution.