Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

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.

The 8-hydroxyquinoline derivative, clioquinol, is an alpha-1 adrenoceptor antagonist

Clioquinol (5-chloro-7-iodo-8-hydroxyquinoline) is an antimicrobial agent whose actions as a zinc or copper ionophore and an iron chelator revived the interest in similar compounds for the treatment of fungal and bacterial infections, neurodegeneration and cancer. Recently, we reported zinc ionophores, including clioquinol, cause vasorelaxation in isolated arteries through mechanisms that involve sensory nerves, endothelium and vascular smooth muscle. Here, we report that clioquinol also uniquely acts as a competitive alpha-1 (α1) adrenoceptor antagonist. We employed ex vivo functional vascular contraction and pharmacological techniques in rat isolated mesenteric arteries, receptor binding assays using stabilized solubilized α1 receptor variants, or wild-type human α1-adrenoceptors transfected in COS-7 cells (African green monkey kidney fibroblast-like cells), and molecular dynamics homology modelling based on the recently published α1A adrenoceptor cryo-EM and α1B crystal structures. At higher concentrations, all ionophores including clioquinol cause a non-competitive antagonism of agonist-mediated contraction due to intracellular zinc delivery, as reported previously. However, at lower concentration ranges, clioquinol has an additional mechanism of competitively inhibiting α1-adrenoceptors that contributes to decreasing vascular contractility. Molecular dynamic simulation showed that clioquinol binds stably to the orthosteric binding site (Asp106) of the receptor, confirming the structural basis for competitive α1-adrenoceptor antagonism by clioquinol.

Curvature Matters: Modeling Calcium Binding to Neutral and Anionic Phospholipid Bilayers

In this work, the influence of membrane curvature on the Ca²⁺ binding to phospholipid bilayers is investigated by means of molecular dynamics simulations. In particular, we compared Ca²⁺ binding to flat, elastically buckled, or uniformly bent zwitterionic and anionic phospholipid bilayers. We demonstrate that Ca²⁺ ions bind preferably to the concave membrane surfaces in both types of bilayers. We also show that the membrane curvature leads to pronounced changes in Ca²⁺ binding including differences in free ion concentrations, lipid coordination distributions, and the patterns of ion binding to different chemical groups of lipids. Moreover, these effects differ substantially for the concave and convex membrane monolayers. Comparison between force fields with either full or scaled charges indicates that charge scaling results in reduction of the Ca²⁺ binding to curved phosphatidylserine bilayers, while for phosphatidylcholine membranes, calcium binds only weakly for both force fields.

Redox regulation of KV7 channels through EF3 hand of calmodulin

Neuronal KV7 channels, important regulators of cell excitability, are among the most sensitive proteins to reactive oxygen species. The S2S3 linker of the voltage sensor was reported as a site-mediating redox modulation of the channels. Recent structural insights reveal potential interactions between this linker and the Ca2+-binding loop of the third EF-hand of calmodulin (CaM), which embraces an antiparallel fork formed by the C-terminal helices A and B, constituting the calcium responsive domain (CRD). We found that precluding Ca2+ binding to the EF3 hand, but not to EF1, EF2, or EF4 hands, abolishes oxidation-induced enhancement of KV7.4 currents. Monitoring FRET (Fluorescence Resonance Energy Transfer) between helices A and B using purified CRDs tagged with fluorescent proteins, we observed that S2S3 peptides cause a reversal of the signal in the presence of Ca2+ but have no effect in the absence of this cation or if the peptide is oxidized. The capacity of loading EF3 with Ca2+ is essential for this reversal of the FRET signal, whereas the consequences of obliterating Ca2+ binding to EF1, EF2, or EF4 are negligible. Furthermore, we show that EF3 is critical for translating Ca2+ signals to reorient the AB fork. Our data are consistent with the proposal that oxidation of cysteine residues in the S2S3 loop relieves KV7 channels from a constitutive inhibition imposed by interactions between the EF3 hand of CaM which is crucial for this signaling.

Cumulative Millisecond-Long Sampling for a Comprehensive Energetic Evaluation of Aqueous Ionic Liquid Effects on Amino Acid Interactions

The interactions of amino acid side-chains confer diverse energetic contributions and physical properties to a protein's stability and function. Various computational tools estimate the effect of changing a given amino acid on the protein's stability based on parametrized (free) energy functions. When parametrized for the prediction of protein stability in water, such energy functions can lead to suboptimal results for other solvents, such as ionic liquids (IL), aqueous ionic liquids (aIL), or salt solutions. However, to our knowledge, no comprehensive data are available describing the energetic effects of aIL on intramolecular protein interactions. Here, we present the most comprehensive set of potential of mean force (PMF) profiles of pairwise protein-residue interactions to date, covering 50 relevant interactions in water, the two biotechnologically relevant aIL [BMIM/Cl] and [BMIM/TfO], and [Na/Cl]. These results are based on a cumulated simulation time of >1 ms. aIL and salt ions can weaken, but also strengthen, specific residue interactions by more than 3 kcal mol^-1, depending on the residue pair, residue-residue configuration, participating ions, and concentration, necessitating considering such interactions specifically. These changes originate from a complex interplay of competitive or cooperative noncovalent ion-residue interactions, changes in solvent structural dynamics, or unspecific charge screening effects and occur at the contact distance but also at larger, solvent-separated distances. This data provide explanations at the atomistic and energetic levels for complex IL effects on protein stability and should help improve the prediction accuracies of computational tools that estimate protein stability based on (free) energy functions.

Recalibrating the calcium trap in amino acid carboxyl groups via classical molecular dynamics simulations

In order to use classical molecular dynamics to complement experiments accurately, it is important to use robust descriptions of the system. The interactions between biomolecules, like aspartic and glutamic acid, and dissolved ions are often studied using standard biomolecular force-fields, where the interactions between biomolecules and cations are often not parameterized explicitly. In this study, we have employed metadynamics simulations to investigate different interactions of Ca with aspartic and glutamic acid and constructed the free energy profiles of Ca²⁺-carboxylate association. Starting from a generally accepted, AMBER-based force field, the association was substantially over and under-estimated, depending on the choice of water model (TIP3P and SPC/fw, respectively). To rectify this discrepancy, we have replaced the default calcium parameters. Additionally, we modified the σ_ij value in the hetero-atomic Lennard-Jones interaction by 0.5% to further improve the interaction between Ca and carboxylate, based on comparison with the experimentally determined association constant for Ca with the carboxylate group of L-aspartic acid. The corrected description retrieved the structural properties of the ion pair in agreement with the original biomolecule - Ca²⁺ interaction in AMBER, whilst also producing an association constant comparable to experimental observations. This refined force field was then used to investigate the interactions between amino acids, calcium and carbonate ions during biogenic and biomimetic calcium carbonate mineralisation.

Probing the Suitability of Different Ca2+ Parameters for Long Simulations of Diisopropyl Fluorophosphatase

Organophosphate hydrolases are promising as potential biotherapeutic agents to treat poisoning with pesticides or nerve gases. However, these enzymes often need to be further engineered in order to become useful in practice. One example of such enhancement is the alteration of enantioselectivity of diisopropyl fluorophosphatase (DFPase). Molecular modeling techniques offer a unique opportunity to address this task rationally by providing a physical description of the substrate-binding process. However, DFPase is a metalloenzyme, and correct modeling of metal cations is a challenging task generally coming with a tradeoff between simulation speed and accuracy. Here, we probe several molecular mechanical parameter combinations for their ability to empower long simulations needed to achieve a quantitative description of substrate binding. We demonstrate that a combination of the Amber19sb force field with the recently developed 12-6 Ca2+ models allows us to both correctly model DFPase and obtain new insights into the DFP binding process.

Computational Analyses of the AtTPC1 (Arabidopsis Two-Pore Channel 1) Permeation Pathway

Two Pore Channels (TPCs) are cation-selective voltage- and ligand-gated ion channels in membranes of intracellular organelles of eukaryotic cells. In plants, the TPC1 subtype forms the slowly activating vacuolar (SV) channel, the most dominant ion channel in the vacuolar membrane. Controversial reports about the permeability properties of plant SV channels fueled speculations about the physiological roles of this channel type. TPC1 is thought to have high Ca²⁺ permeability, a conclusion derived from relative permeability analyses using the Goldman–Hodgkin–Katz (GHK) equation. Here, we investigated in computational analyses the properties of the permeation pathway of TPC1 from Arabidopsis thaliana. Using the crystal structure of AtTPC1, protein modeling, molecular dynamics (MD) simulations, and free energy calculations, we identified a free energy minimum for Ca²⁺, but not for K⁺, at the luminal side next to the selectivity filter. Residues D269 and E637 coordinate in particular Ca²⁺ as demonstrated in in silico mutagenesis experiments. Such a Ca²⁺-specific coordination site in the pore explains contradicting data for the relative Ca²⁺/K⁺ permeability and strongly suggests that the Ca²⁺ permeability of SV channels is largely overestimated from relative permeability analyses. This conclusion was further supported by in silico electrophysiological studies showing a remarkable permeation of K⁺ but not Ca²⁺ through the open channel.

NEDD8 Deamidation Inhibits Cullin RING Ligase Dynamics

Cullin-RING ligases (CRLs) are a significant subset of Ubiquitin E3 ligases that regulate multiple cellular substrates involved in innate immunity, cytoskeleton modeling, and cell cycle. The glutamine deamidase Cycle inhibitory factor (Cif) from enteric bacteria inactivates CRLs to modulate these processes in the host cell. The covalent attachment of a Ubiquitin-like protein NEDD8 catalytically activates CRLs by driving conformational changes in the Cullin C-terminal domain (CTD). NEDDylation results in a shift from a compact to an open CTD conformation through non-covalent interactions between NEDD8 and the WHB subdomain of CTD, eliminating the latter’s inhibitory interactions with the RING E3 ligase-Rbx1/2. It is unknown whether the non-covalent interactions are sufficient to stabilize Cullin CTD’s catalytic conformation. We studied the dynamics of Cullin-CTD in the presence and absence of NEDD8 using atomistic molecular dynamics (MD) simulations. We uncovered that NEDD8 engages in non-covalent interactions with 4HB/αβ subdomains in Cullin-CTD to promote open conformations. Cif deamidates glutamine 40 in NEDD8 to inhibit the conformational change in CRLs by an unknown mechanism. We investigated the effect of glutamine deamidation on NEDD8 and its interaction with the WHB subdomain post-NEDDylation using MD simulations and NMR spectroscopy. Our results suggest that deamidation creates a new intramolecular salt bridge in NEDD8 to destabilize the NEDD8/WHB complex and reduce CRL activity.

Ground-State Structures of Hydrated Calcium Ion Clusters From Comprehensive Genetic Algorithm Search

We searched the lowest-energy structures of hydrated calcium ion clusters Ca2+(H2O)n (n = 10-18) in the whole potential energy surface by the comprehensive genetic algorithm (CGA). The lowest-energy structures of Ca2+(H2O)10-12 clusters show that Ca2+ is always surrounded by six H2O molecules in the first shell. The number of first-shell water molecules changes from six to eight at n = 12. In the range of n = 12-18, the number of first-shell water molecules fluctuates between seven and eight, meaning that the cluster could pack the water molecules in the outer shell even though the inner shell is not full. Meanwhile, the number of water molecules in the second shell and the total hydrogen bonds increase with an increase in the cluster size. The distance between Ca2+ and the adjacent water molecules increases, while the average adjacent O-O distance decreases as the cluster size increases, indicating that the interaction between Ca2+ and the adjacent water molecules becomes weaker and the interaction between water molecules becomes stronger. The interaction energy and natural bond orbital results show that the interaction between Ca2+ and the water molecules is mainly derived from the interaction between Ca2+ and the adjacent water molecules. The charge transfer from the lone pair electron orbital of adjacent oxygen atoms to the empty orbital of Ca2+ plays a leading role in the interaction between Ca2+ and water molecules.

Sulfation and Calcium Favor Compact Conformations of Chondroitin in Aqueous Solutions

The effects of sulfation and calcium cations (Ca²⁺) on the atomic-resolution conformational properties of chondroitin carbohydrate polymers in aqueous solutions are not well studied owing to experimental challenges. Here, we compare all-atom explicit-solvent molecular dynamics simulations results for pairs of O-type (nonsulfated) and A-type (GlcNAc 4-O-sulfated) chondroitin 20-mers in 140 mM NaCl with and without Ca²⁺ and find that both sulfation and Ca²⁺ favor more compact polymer conformations. We also show that subtle differences in force-field parametrization can have dramatic effects on Ca²⁺ binding to chondroitin carboxylate and sulfate functional groups and thereby determine Ca²⁺-mediated intra- and interstrand association. In addition to providing an atomic-resolution picture of the interaction of Ca²⁺ with sulfated and nonsulfated chondroitin polymers, the molecular dynamics data emphasize the importance of careful force-field parametrization to balance ion-water and ion-chondroitin interactions and suggest additional parametrization efforts to tune interactions involving sulfate.

The molecular basis for the pH-dependent calcium affinity of the pattern recognition receptor langerin

The C-type lectin receptor langerin plays a vital role in the mammalian defense against invading pathogens. Langerin requires a Ca²⁺ cofactor, the binding affinity of which is regulated by pH. Thus, Ca²⁺ is bound when langerin is on the membrane but released when langerin and its pathogen substrate traffic to the acidic endosome, allowing the substrate to be degraded. The change in pH is sensed by protonation of the allosteric pH sensor histidine H294. However, the mechanism by which Ca²⁺ is released from the buried binding site is not clear. We studied the structural consequences of protonating H294 by molecular dynamics simulations (total simulation time: about 120 μs) and Markov models. We discovered a relay mechanism in which a proton is moved into the vicinity of the Ca²⁺-binding site without transferring the initial proton from H294. Protonation of H294 unlocks a conformation in which a protonated lysine side chain forms a hydrogen bond with a Ca²⁺-coordinating aspartic acid. This destabilizes Ca²⁺ in the binding pocket, which we probed by steered molecular dynamics. After Ca²⁺ release, the proton is likely transferred to the aspartic acid and stabilized by a dyad with a nearby glutamic acid, triggering a conformational transition and thus preventing Ca²⁺ rebinding. These results show how pH regulation of a buried orthosteric binding site from a solvent-exposed allosteric pH sensor can be realized by information transfer through a specific chain of conformational arrangements.

A practical guide to biologically relevant molecular simulations with charge scaling for electronic polarization

Molecular simulations can elucidate atomistic-level mechanisms of key biological processes, which are often hardly accessible to experiment. However, the results of the simulations can only be as trustworthy as the underlying simulation model. In many of these processes, interactions between charged moieties play a critical role. Current empirical force fields tend to overestimate such interactions, often in a dramatic way, when polyvalent ions are involved. The source of this shortcoming is the missing electronic polarization in these models. Given the importance of such biomolecular systems, there is great interest in fixing this deficiency in a computationally inexpensive way without employing explicitly polarizable force fields. Here, we review the electronic continuum correction approach, which accounts for electronic polarization in a mean-field way, focusing on its charge scaling variant. We show that by pragmatically scaling only the charged molecular groups, we qualitatively improve the charge-charge interactions without extra computational costs and benefit from decades of force field development on biomolecular force fields.

Investigating dual Ca2+ modulation of the ryanodine receptor 1 by molecular dynamics simulation

The ryanodine receptors (RyR) are essential to calcium signaling in striated muscles. A deep understanding of the complex Ca²⁺ -activation/inhibition mechanism of RyRs requires detailed structural and dynamic information for RyRs in different functional states (eg, with Ca²⁺ bound to activating or inhibitory sites). Recently, high-resolution structures of the RyR isoform 1 (RyR1) were solved by cryo-electron microscopy, revealing the location of a Ca²⁺ binding site for activation. Toward elucidating the Ca²⁺ -modulation mechanism of RyR1, we performed extensive molecular dynamics simulation of the core RyR1 structure in the presence and absence of activating and solvent Ca²⁺ (total simulation time is >5 μs). In the presence of solvent Ca²⁺ , Ca²⁺ binding to the activating site enhanced dynamics of RyR1 with higher inter-subunit flexibility, asymmetric inter-subunit motions, outward domain motions and partial pore dilation, which may prime RyR1 for subsequent channel opening. In contrast, the solvent Ca²⁺ alone reduced dynamics of RyR1 and led to inward domain motions and pore contraction, which may cause inhibition. Combining our simulation with the map of disease mutation sites in RyR1, we constructed a wiring diagram of key domains coupled via specific hydrogen bonds involving the mutation sites, some of which were modulated by Ca²⁺ binding. The structural and dynamic information gained from this study will inform future mutational and functional studies of RyR1 activation and inhibition by Ca²⁺ .

Influence of Calcium Binding on Conformations and Motions of Anionic Polyamino Acids. Effect of Side Chain Length

Investigation of the effect of CaCl2 salt on conformations of two anionic poly(amino acids) with different side chain lengths, poly-(α-l glutamic acid) (PGA) and poly-(α-l aspartic acid) (PASA), was performed by atomistic molecular dynamics (MD) simulations. The simulations were performed using both unbiased MD and the Hamiltonian replica exchange (HRE) method. The results show that at low CaCl2 concentration adsorption of Ca2+ ions lead to a significant chain size reduction for both PGA and PASA. With the increase in concentration, the chains sizes partially recover due to electrostatic repulsion between the adsorbed Ca2+ ions. Here, the side chain length becomes important. Due to the longer side chain and its ability to distance the charged groups with adsorbed ions from both each other and the backbone, PGA remains longer in the collapsed state as the CaCl2 concentration is increased. The analysis of the distribution of the mineral ions suggests that both poly(amino acids) should induce the formation of mineral with the same structure of the crystal cell.

Divalent cations bind to phosphoinositides to induce ion and isomer specific propensities for nano-cluster initiation in bilayer membranes

We report all-atom molecular dynamics simulations of asymmetric bilayers containing phosphoinositides in the presence of monovalent and divalent cations. We have characterized the molecular mechanism by which these divalent cations interact with phosphoinositides. Ca²⁺ desolvates more readily, consistent with single-molecule calculations, and forms a network of ionic-like bonds that serve as a 'molecular glue' that allows a single ion to coordinate with up to three phosphatidylinositol-(4,5)-bisphosphate (PI(4, 5)P₂) lipids. The phosphatidylinositol-(3,5)-bisphosphate isomer shows no such effect and neither does PI(4, 5)P₂ in the presence of Mg²⁺. The resulting network of Ca²⁺-mediated lipid-lipid bonds grows to span the entire simulation space and therefore has implications for the lateral distribution of phosophoinositides in the bilayer. We observe context-specific differences in lipid diffusion rates, lipid surface densities and bilayer structure. The molecular-scale delineation of ion-lipid arrangements reported here provides insight into similar nanocluster formation induced by peripheral proteins to regulate the formation of functional signalling complexes on the membrane.

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

Classical molecular dynamics simulations reveal size-dependent trends of alkaline earth metal ions binding to DNA are due to ion size and hydration behavior.

Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction

The parametrization of classical nonbonded models of metal ions has been widely addressed in the recent years. Despite the continuous development of novel and more physically inspired functional forms, the 12-6 Lennard-Jones plus Coulomb potential is still the most adopted force field in molecular dynamics (MD) codes, owing to its simple form and easy implementation. However, due to the integer formal charge, unpolarizable force fields of ions may suffer from overestimated interatomic electrostatic interactions, leading to nonphysical clustering or repulsion between such full charges. The electronic continuum correction (ECC) can fix this problem through a simple inclusion of solvent polarization effects via ionic charge rescaling. In this work, the development of novel nonbonded models for mono, divalent, and highly charged metal ions is presented. For each metal species, the ionic charge has been scaled, according to the ECC. Lennard-Jones parameters have been optimized using experimental structural and thermodynamic properties as target quantities. Performances of the proposed models are discussed and compared with the literature data, while transferability attitudes among different and well-known water models are evaluated. © 2019 Wiley Periodicals, Inc.

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

A Kirkwood-Buff derived force field for alkaline earth halide salts

The activity and function of many macromolecules in cellular environments are coupled with the binding of divalent ions such as calcium or magnesium. In principle, computer simulations can be used to understand the molecular level aspects of how many important macromolecules interact with ions. However, most of the force fields currently available often fail to accurately reproduce the properties of divalent ions in aqueous environments. Here we develop classical non-polarizable force fields for the aqueous alkaline earth metal halides (MX₂), where M = Mg²⁺, Ca²⁺, Sr²⁺, Ba²⁺ and X = Cl^-, Br^-, I^-, which can be used in bimolecular simulations and which are compatible with the Simple Point Charge/Extended (SPC/E) water model. The force field parameters are specifically developed to reproduce the experimental Kirkwood-Buff integrals for aqueous solutions and thereby the experimental activity derivatives, partial molar volumes, and excess coordination numbers. This ensures that a reasonable balance between ion-ion, ion-water, and water-water distributions is obtained. However, this requires a scaling of the cation to water oxygen interaction strength in order to accurately reproduce the integrals. The scaling factors developed for chloride salts are successfully transferable to the bromide and iodide salts. Use of these new models leads to reasonable diffusion constants and dielectric decrements. However, the performance of the models decreases with increasing salt concentration (>4m), and simulations of the pure crystals exhibited unstable behavior.

Ion Permeability of a Microtubule in Neuron Environment

Microtubules, constituted by end-to-end negatively charged α- and β-tubulin dimers, are long, hollow, pseudohelical cylinders with internal and external diameters of about 16 and 26 nm, respectively, and widely exist in cell cytoplasm, neuron axons, and dendrites. Although their structural functions in physiological processes, such as cell mitosis, cell motility, and motor protein transport, have been widely accepted, their role in neuron activity remains attractively elusive. Here we show a new function of microtubules: they can generate instant response to a calcium pulse because of their specific permeability for ions. Our comprehensive simulations from all-atom molecular dynamics to potential of mean force and continuum modeling reveal that K⁺ and Na⁺ ions can permeate through the nanopores in the microtubule wall easily, while Ca²⁺ ions are blocked by the wall with a much higher free energy barrier. These cations are adsorbed to the surfaces of the wall with affinity decreasing in the sequence Ca²⁺, Na⁺, and K⁺. As a result, when the concentration of Ca²⁺ ions increases outside the microtubule during neuronal excitation, K⁺ and Na⁺ ions will be driven into the microtubule, triggering subsequent axial ion redistribution within the microtubule. The results shed light on the possibility of the ion-permeable microtubules being involved in neural signal processing.

New tricks for old dogs: improving the accuracy of biomolecular force fields by pair-specific corrections to non-bonded interactions

In contrast to ordinary polymers, the vast majority of biological macromolecules adopt highly ordered three-dimensional structures that define their functions. The key to folding of a biopolymer into a unique 3D structure or to an assembly of several biopolymers into a functional unit is a delicate balance between the attractive and repulsive forces that also makes such self-assembly reversible under physiological conditions. The all-atom molecular dynamics (MD) method has emerged as a powerful tool for studies of individual biomolecules and their functional assemblies, encompassing systems of ever increasing complexity. However, advances in parallel computing technology have outpaced the development of the underlying theoretical models-the molecular force fields, pushing the MD method into an untested territory. Recent tests of the MD method have found the most commonly used molecular force fields to be out of balance, overestimating attractive interactions between charged and hydrophobic groups, which can promote artificial aggregation in MD simulations of multi-component protein, nucleic acid, and lipid systems. One route towards improving the force fields is through the NBFIX corrections method, in which the intermolecular forces are calibrated against experimentally measured quantities such as osmotic pressure by making atom pair-specific adjustments to the non-bonded interactions. In this article, we review development of the NBFIX (Non-Bonded FIX) corrections to the AMBER and CHARMM force fields and discuss their implications for MD simulations of electrolyte solutions, dense DNA systems, Holliday junctions, protein folding, and lipid bilayer membranes.

Efficient Osmotic Pressure Calculations Using Coarse-Grained Molecular Simulations

Osmotic pressure data is increasingly used to parametrize all-atom simulation Force Fields (FFs), leading to large computational cost for larger molecules. Here, we show that the osmotic pressure can be calculated precisely using transferable coarse-grained FFs obtained from short atomistic simulations using an inhomogeneously regularized coarse-graining procedure. This is demonstrated for carbohydrates, where compared to the equivalent atomistic system, an increase of the computational efficiency by a factor of ≈500 is achieved.

Counting the ions surrounding nucleic acids

Nucleic acids are strongly negatively charged, and thus electrostatic interactions—screened by ions in solution—play an important role in governing their ability to fold and participate in biomolecular interactions. The negative charge creates a region, known as the ion atmosphere, in which cation and anion concentrations are perturbed from their bulk values. Ion counting experiments quantify the ion atmosphere by measuring the preferential ion interaction coefficient: the net total number of excess ions above, or below, the number expected due to the bulk concentration. The results of such studies provide important constraints on theories, which typically predict the full three-dimensional distribution of the screening cloud. This article reviews the state of nucleic acid ion counting measurements and critically analyzes their ability to test both analytical and simulation-based models.

Reparameterization of Solute-Solute Interactions for Amino Acid-Sugar Systems Using Isopiestic Osmotic Pressure Molecular Dynamics Simulations

AMBER/GLYCAM and CHARMM are popular force fields for simulations of amino acids and sugars. Here we report excessively attractive amino acid-sugar interactions in both force fields, and corrections to nonbonded interactions that match experimental osmotic pressures of mixed aqueous solutions of diglycine and sucrose. The modified parameters also improve the ΔGtrans of diglycine from water to aqueous sucrose and, with AMBERff99SB/GLYCAM06, eliminate a caging effect seen in previous simulations of the protein ubiquitin with glucose.

Probing the structure and in silico stability of cargo loaded DNA icosahedra using MD simulations

Platonic solids such as polyhedra based on DNA have been deployed for multifarious applications such as RNAi delivery, biological targeting and bioimaging. All of these applications hinge on the capability of DNA polyhedra for molecular display with high spatial precision. Therefore high resolution structural models of such polyhedra are critical to widen their applications in both materials and biology. Here, we present an atomistic model of a well-characterized DNA icosahedron, with demonstrated versatile functionalities in biological systems. We study the structure and dynamics of this DNA icosahedron using fully atomistic molecular dynamics (MD) simulation in explicit water and ions. The major modes of internal motion have been identified using principal component analysis. We provide a quantitative estimate of the radius of gyration (R_g), solvent accessible surface area (SASA) and volume of the icosahedron which is essential to estimate its maximal cargo carrying capacity. Importantly, our simulation of gold nanoparticles (AuNPs) encapsulated within DNA icosahedra revealed enhanced stability of the AuNP loaded DNA icosahedra compared to empty icosahedra. This is consistent with the experimental results that show high yields of cargo-encapsulated DNA icosahedra that have led to its diverse applications for precision targeting. These studies reveal that the stabilizing interactions between the cargo and the DNA scaffold powerfully position DNA polyhedra as targetable nanocapsules for payload delivery. These insights can be exploited for precise molecular display for diverse biological applications.

Refined Parameterization of Nonbonded Interactions Improves Conformational Sampling and Kinetics of Protein Folding Simulations

Recent advances in computational technology have enabled brute-force molecular dynamics (MD) simulations of protein folding using physics-based molecular force fields. The extensive sampling of protein conformations afforded by such simulations revealed, however, considerable compaction of the protein conformations in the unfolded state, which is inconsistent with experiment. Here, we show that a set of surgical corrections to nonbonded interactions between amine nitrogen-carboxylate oxygen and aliphatic carbon-carbon atom pairs can considerably improve the realism of protein folding simulations. Specifically, we show that employing our corrections in ∼500 μs all-atom replica-exchange MD simulations of the WW domain and villin head piece proteins increases the size of the denatured proteins' conformations and does not destabilize the native conformations of the proteins. In addition to making the folded conformations a global minimum of the respective free energy landscapes at room temperature, our corrections also make the free energy landscape smoother, considerably accelerating the folding kinetics and, hence, reducing the computational expense of a protein folding simulation.

Understanding Sodium Channel Function and Modulation Using Atomistic Simulations of Bacterial Channel Structures

Sodium channels are chief proteins involved in electrical signaling in the nervous system, enabling critical functions like heartbeat and brain activity. New high-resolution X-ray structures for bacterial sodium channels have created an opportunity to see how these proteins operate at the molecular level. An important challenge to overcome is establishing relationships between the structures and functions of mammalian and bacterial channels. Bacterial sodium channels are known to exhibit the main structural features of their mammalian counterparts, as well as several key functional characteristics, including selective ion conduction, voltage-dependent gating, pore-based inactivation and modulation by local anesthetic, antiarrhythmic and antiepileptic drugs. Simulations have begun to shed light on each of these features in the past few years. Despite deviations in selectivity signatures for bacterial and mammalian channels, simulations have uncovered the nature of the multiion conduction mechanism associated with Na(+) binding to a high-field strength site established by charged glutamate side chains. Simulations demonstrated a surprising level of flexibility of the protein, showing that these side chains are active participants in the permeation process. They have also uncovered changes in protein structure, leading to asymmetrical collapses of the activation gate that have been proposed to correspond to inactivated structures. These observations offer the potential to examine the mechanisms of state-dependent drug activity, focusing on pore-blocking and pore-based slow inactivation in bacterial channels, without the complexities of inactivation on multiple timescales seen in eukaryotic channels. Simulations have provided molecular views of the interactions of drugs, consistent with sites predicted in mammalian channels, as well as a wealth of other sites as potential new drug targets. In this chapter, we survey the new insights into sodium channel function that have emerged from studies of simpler bacterial channels, which provide an excellent learning platform, and promising avenues for mechanistic discovery and pharmacological development.

Molecular electrometer and binding of cations to phospholipid bilayers

The molecular electrometer – reorientation of lipid head due to bound charge – allows direct quantitative vetting of simulations against noninvasive NMR experiments; most simulation models overestimated lipid–cation affinities.