A Kirkwood-Buff force field for the aromatic amino acids

Kirkwood-Buff Analysis of Binary and Ternary Systems Consisting of Alcohols (Methanol, Ethanol, 1-Propanol, and 2-Propanol), Water, and n-Hexane to Understand the Formation of Surfactant-Free Microemulsions

Surfactant-free microemulsion (SFME) represents a class of fluid mixtures that can form microheterogeneous structures without detergents, offering an environmentally benign alternative to traditional microemulsions. However, the formation mechanism is still elusive. This work applies the Kirkwood-Buff theory to mixtures of alcohols, water, and n-hexane to elucidate the SFME formation mechanism. To ensure robust calculation of the Kirkwood-Buff integrals (KBIs), we construct a data set of densities and excess free energies of binary and ternary systems. Multiple excess Gibbs free energy models are assessed against this data set to select the most suitable model reproducing the experimental results. In addition, we introduce statistical methods to determine the optimal polynomial order of the Redlich-Kister correlation for the excess volume data. We first validate our methodology in binary systems. Then, we extend the calculation method to ternary mixtures. The KBI calculation results reveal that the alcohol-hexane and water-hexane interactions do not significantly affect SFME formation. In contrast, the interplay among water-water, water-alcohol, and alcohol-alcohol interactions critically influences the ability of a liquid mixture to form SFME structures. SFME systems exhibit the facile formation of water aggregates enveloped by alcohols, whereas non-SFME systems demonstrate homogeneous alcohol/water droplets dispersed in an oil continuous medium.

Impact of Cd(II) on the stability of human uracil DNA glycosylase enzyme; an implication of molecular dynamics trajectories on stability analysis

Uracil DNA glycosylase is a key enzyme that identifies and removes damaged bases from DNA in the base excision repair pathway. Experimentalists have identified the possibility of Cd(II) reducing the activity of human uracil DNA glycosylase (hUNG) by binding with the enzyme replacing the catalytic water molecule. The present study focus on the stability variation of the enzyme in the presence and absence of Cd(II) and confirms the reported results with the stability analysis done using molecular dynamic (MD) simulation trajectories. The CavityPlus web server identified seven cavities for the free enzyme as possible binding sites and a cavity containing the active site of the enzyme as the best binding cavity for a ligand. Based on the CavityPlus results and the previously reported work, a free hUNG system and two systems of the enzyme with Cd(II); one with Cd(II) replacing the catalytic water molecule in the active site of the enzyme and the other replacing a non-catalytic water molecule in the active site were generated for the simulation. The simulation trajectories were used for the structural stability analysis of the enzyme in all three systems. The binding free energy of the Cd(II) with the enzyme was calculated using molecular mechanics Poisson Boltzmann surface area method. The results showed that the enzyme achieves comparatively high stability with the removal of catalytic water of the enzyme by Cd(II). Therefore, this supports the previously reported idea that Cd(II) replaces catalytic water molecules and affects enzyme activity.

Simulating Polyproline II-Helix-Rich Peptides with the Latest Kirkwood-Buff Force Field: A Direct Comparison with AMBER and CHARMM

We simulated the dynamics of a set of peptides characterized by ensembles rich in PPII-helical content, to assess the ability of the most recent Kirkwood-Buff force field (KBFF20) to sample this conformational peculiarity. KBFF has been previously shown to capably reproduce experimental dimensions of disordered proteins, while being limited in confidently sampling structured proteins. Further development of the force field bridged this gap. It is however still unknown what are the main differences between KBFF and AMBER/CHARMM force fields. A direct comparison is now possible as both AMBER/CHARMM force fields have been used to sample peptides rich in PPII-helical content. We found that KBFF20 samples' PPII-helical content qualitatively matches both AMBER and CHARMM force fields, with the main difference being the KBFF ability to populate the α_R region of the Ramachandran plot in the set of simulated peptides. Overall, KBFF20 is a well-balanced force field, able to sample the dynamics of both structured and unstructured proteins.

Improving Force Field Accuracy by Training against Condensed-Phase Mixture Properties

Developing a sufficiently accurate classical force field representation of molecules is key to realizing the full potential of molecular simulations as a route to gaining a fundamental insight into a broad spectrum of chemical and biological phenomena. This is only possible, however, if the many complex interactions between molecules of different species in the system are accurately captured by the model. Historically, the intermolecular van der Waals (vdW) interactions have primarily been trained against densities and enthalpies of vaporization of pure (single-component) systems, with occasional usage of hydration free energies. In this study, we demonstrate how including physical property data of binary mixtures can better inform these parameters, encoding more information about the underlying physics of the system in complex chemical mixtures. To demonstrate this, we retrain a select number of Lennard-Jones parameters describing the vdW interactions of the OpenFF 1.0.0 (Parsley) fixed charge force field against training sets composed of densities and enthalpies of mixing for binary liquid mixtures as well as densities and enthalpies of vaporization of pure liquid systems and assess the performance of each of these combinations. We show that retraining against the mixture data improves the force field's ability to reproduce mixture properties, including solvation free energies, correcting some systematic errors that exist when training vdW interactions against properties of pure systems only.

Biomolecular Simulations with the Three-Dimensional Reference Interaction Site Model with the Kovalenko-Hirata Closure Molecular Solvation Theory

The statistical mechanics-based 3-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH) molecular solvation theory has proven to be an essential part of a multiscale modeling framework, covering a vast region of molecular simulation techniques. The successful application ranges from the small molecule solvation energy to the bulk phase behavior of polymers, macromolecules, etc. The 3D-RISM-KH successfully predicts and explains the molecular mechanisms of self-assembly and aggregation of proteins and peptides related to neurodegeneration, protein-ligand binding, and structure-function related solvation properties. Upon coupling the 3D-RISM-KH theory with a novel multiple time-step molecular dynamic (MD) of the solute biomolecule stabilized by the optimized isokinetic Nosé-Hoover chain thermostat driven by effective solvation forces obtained from 3D-RISM-KH and extrapolated forward by generalized solvation force extrapolation (GSFE), gigantic outer time-steps up to picoseconds to accurately calculate equilibrium properties were obtained in this new quasidynamics protocol. The multiscale OIN/GSFE/3D-RISM-KH algorithm was implemented in the Amber package and well documented for fully flexible model of alanine dipeptide, miniprotein 1L2Y, and protein G in aqueous solution, with a solvent sampling rate ~150 times faster than a standard MD simulation in explicit water. Further acceleration in computation can be achieved by modifying the extent of solvation layers considered in the calculation, as well as by modifying existing closure relations. This enhanced simulation technique has proven applications in protein-ligand binding energy calculations, ligand/solvent binding site prediction, molecular solvation energy calculations, etc. Applications of the RISM-KH theory in molecular simulation are discussed in this work.

Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20

A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.

Thermodynamic, structural and dynamic properties of selected non-associative neat liquids

Non-associative neat liquids benzene, acetone and carbon tetrachloride have been examined via molecular dynamics simulations. Several models of each neat liquid have been simulated, and selected thermodynamic and structural results are presented. However, the models have been compared mainly in terms of their dynamic quantities. Since models are rarely parametrized with the dynamic properties in mind, the principal goal of this work is to present quantities such as the power spectra, rotational correlation functions and relaxation times, diffusion coefficients and self and distinct parts of the van Hove functions in relation to available experimental data. The general trends of the calculated data provide a benchmark for the behavior of neat simple liquids which will be built upon in the cases of mixtures with associative liquids.

Coarse-Grained Molecular Model for the Glycosylphosphatidylinositol Anchor with and without Protein

Glycosylphosphatidylinositol (GPI) anchors are a unique class of complex glycolipids that anchor a great variety of proteins to the extracellular leaflet of plasma membranes of eukaryotic cells. These anchors can exist either with or without an attached protein called GPI-anchored protein (GPI-AP) both in vitro and in vivo. Although GPIs are known to participate in a broad range of cellular functions, it is to a large extent unknown how these are related to GPI structure and composition. Their conformational flexibility and microheterogeneity make it difficult to study them experimentally. Simplified atomistic models are amenable to all-atom computer simulations in small lipid bilayer patches but not suitable for studying their partitioning and trafficking in complex and heterogeneous membranes. Here, we present a coarse-grained model of the GPI anchor constructed with a modified version of the MARTINI force field that is suited for modeling carbohydrates, proteins, and lipids in an aqueous environment using MARTINI's polarizable water. The nonbonded interactions for sugars were reparametrized by calculating their partitioning free energies between polar and apolar phases. In addition, sugar-sugar interactions were optimized by adjusting the second virial coefficients of osmotic pressures for solutions of glucose, sucrose, and trehalose to match with experimental data. With respect to the conformational dynamics of GPI-anchored green fluorescent protein, the accessible time scales are now at least an order of magnitude larger than for the all-atom system. This is particularly important for fine-tuning the mutual interactions of lipids, carbohydrates, and amino acids when comparing to experimental results. We discuss the prospective use of the coarse-grained GPI model for studying protein-sorting and trafficking in membrane models.

Development of a Force Field for the Simulation of Single-Chain Proteins and Protein-Protein Complexes

The accuracy of atomistic physics-based force fields for the simulation of biological macromolecules has typically been benchmarked experimentally using biophysical data from simple, often single-chain systems. In the case of proteins, the careful refinement of force field parameters associated with torsion-angle potentials and the use of improved water models have enabled a great deal of progress toward the highly accurate simulation of such monomeric systems in both folded and, more recently, disordered states. In living organisms, however, proteins constantly interact with other macromolecules, such as proteins and nucleic acids, and these interactions are often essential for proper biological function. Here, we show that state-of-the-art force fields tuned to provide an accurate description of both ordered and disordered proteins can be limited in their ability to accurately describe protein-protein complexes. This observation prompted us to perform an extensive reparameterization of one variant of the Amber protein force field. Our objective involved refitting not only the parameters associated with torsion-angle potentials but also the parameters used to model nonbonded interactions, the specification of which is expected to be central to the accurate description of multicomponent systems. The resulting force field, which we call DES-Amber, allows for more accurate simulations of protein-protein complexes, while still providing a state-of-the-art description of both ordered and disordered single-chain proteins. Despite the improvements, calculated protein-protein association free energies still appear to deviate substantially from experiment, a result suggesting that more fundamental changes to the force field, such as the explicit treatment of polarization effects, may simultaneously further improve the modeling of single-chain proteins and protein-protein complexes.

Thermodynamic Analysis of n-Hexane-Ethanol Binary Mixtures Using the Kirkwood-Buff Theory

A complete thermodynamic analysis of mixtures consisting of molecules with complex chemical constitution can be rather demanding. The Kirkwood-Buff theory of solutions allows the estimation of thermodynamic properties, which cannot be directly extracted from atomistic simulations, such as the Gibbs energy of mixing (Δ_mix G). In this work, we perform molecular dynamics simulations of n-hexane-ethanol binary mixtures in the liquid state under two temperature-pressure conditions and at various mole fractions. On the basis of the recently published methodology of Galata [ Fluid Phase Equilib. 2018 , 470 , 25 - 37 ] , we first calculate the Kirkwood-Buff integrals in the isothermal-isobaric ( NpT) ensemble, identifying how system size affects their estimation. We then extract the activity coefficients, excess Gibbs energy, excess enthalpy, and excess entropy for the n-hexane-ethanol binary mixtures we simulate. We employ two approaches for quantifying composition fluctuations: one based on counting molecular centers of mass and a second one based on counting molecular segments. Results from the two approaches are practically indistinguishable. We compare our results against predictions of vapor-liquid equilibria obtained in a previous simulation work using the same force field, as well as with experimental data, and find very good agreement. In addition, we develop a simple methodology to identify the hydrogen bonds between ethanol molecules and analyze their effects on mixing properties.

Evaluating Force Field Performance in Thermodynamic Calculations of Cyclodextrin Host-Guest Binding: Water Models, Partial Charges, and Host Force Field Parameters

Computational prediction of noncovalent binding free energies with methods based on molecular mechanical force fields has become increasingly routine in drug discovery projects, where they promise to speed the discovery of small molecule ligands to bind targeted proteins with high affinity. Because the reliability of free energy methods still has significant room for improvement, new force fields, or modifications of existing ones, are regularly introduced with the aim of improving the accuracy of molecular simulations. However, comparatively little work has been done to systematically assess how well force fields perform, particularly in relation to the calculation of binding affinities. Hardware advances have made these calculations feasible, but comprehensive force field assessments for protein–ligand sized systems still remain costly. Here, we turn to cyclodextrin host–guest systems, which feature many hallmarks of protein–ligand binding interactions but are generally much more tractable due to their small size. We present absolute binding free energy and enthalpy calculations, using the attach-pull-release (APR) approach, on a set of 43 cyclodextrin-guest pairs for which experimental ITC data are available. The test set comprises both α- and β-cyclodextrin hosts binding a series of small organic guests, each with one of three functional groups: ammonium, alcohol, or carboxylate. Four water models are considered (TIP3P, TIP4Pew, SPC/E, and OPC), along with two partial charge assignment procedures (RESP and AM1-BCC) and two cyclodextrin host force fields. The results suggest a complex set of considerations when choosing a force field for biomolecular simulations. For example, some force field combinations clearly outperform others at the binding enthalpy calculations but not for the binding free energy. Additionally, a force field combination which we expected to be the worst performer gave the most accurate binding free energies – but the least accurate binding enthalpies. The results have implications for the development of improved force fields, and we propose this test set, and potential future elaborations of it, as a powerful validation suite to evaluate new force fields and help guide future force field development.

Lifshitz phase: the microscopic structure of aqueous and ethanol mixtures of 1,n-diols

We study binary mixtures of ethylene glycol and 1,3-propandiol with water or ethanol using computer simulations.

Osmotic Pressure Simulations of Amino Acids and Peptides Highlight Potential Routes to Protein Force Field Parameterization

Recent molecular dynamics (MD) simulations of proteins have suggested that common force fields overestimate the strength of amino acid interactions in aqueous solution. In an attempt to determine the causes of these effects, we have measured the osmotic coefficients of a number of amino acids using the AMBER ff99SB-ILDN force field with two popular water models, and compared the results with available experimental data. With TIP4P-Ew water, interactions between aliphatic residues agree well with experiment, but interactions of the polar residues serine and threonine are found to be excessively attractive. For all tested amino acids, the osmotic coefficients are lower when the TIP3P water model is used. Additional simulations performed on charged amino acids indicate that the osmotic coefficients are strongly dependent on the parameters assigned to the salt ions, with a reparameterization of the sodium/carboxylate interaction reported by the Aksimentiev group significantly improving description of the osmotic coefficient for glutamate. For five neutral amino acids, we also demonstrate a decrease in solute-solute attractions using the recently reported TIP4P-D water model and using the KBFF force field. Finally, we show that for four two-residue peptides improved agreement with experiment can be achieved by rederiving the partial charges for each peptide.

Optimizing Solute-Solute Interactions in the GLYCAM06 and CHARMM36 Carbohydrate Force Fields Using Osmotic Pressure Measurements

GLYCAM06 and CHARMM36 are successful force fields for modeling carbohydrates. To correct recently identified deficiencies with both force fields, we adjusted intersolute nonbonded parameters to reproduce the experimental osmotic coefficient of glucose at 1 M. The modified parameters improve behavior of glucose and sucrose up to 4 M and improve modeling of a dextran 55-mer. While the modified parameters may not be applicable to all carbohydrates, they highlight the use of osmotic simulations to optimize force fields.

From solutions to molecular emulsions

Concentration fluctuations play an important role in the statistical description of the stability of liquids, particularly in the neighborhood of phase transitions. Classical thermodynamics is blind to fluctuations, and statistical thermodynamics is required to fully understand quantities such as the isothermal compressibility or heat capacity, by linking them to fluctuations of appropriate statistical microscopic quantities and showing that they are response functions. This is illustrated by the seminal Kirkwood–Buff theory of solutions. However, the existence of micro-heterogeneous structures, particularly in aqueous mixtures, which leads to large Kirkwood–Buff integrals, suggest that micro-heterogeneity is a form of concentration fluctuation. This interpretation becomes difficult to accept when extrapolated to larger micro-heterogeneous structures such as micellar aggregates in micro-emulsions. By analyzing how different methods, experimental, computer experiments and theoretical approaches deal with the underlying duality behind these two physical manifestations, we put in evidence the need to reconsider the description of liquids by incorporating the description of emergent “objects”, such as the micro-heterogeneous structures from a molecular point of view. On this path, the concept of “molecular emulsion” allows to describe in a unified way all type of disordered liquids, from solutions to the organized liquids of soft matter.

Kirkwood-Buff Approach Rescues Overcollapse of a Disordered Protein in Canonical Protein Force Fields

Understanding the function of intrinsically disordered proteins is intimately related to our capacity to correctly sample their conformational dynamics. So far, a gap between experimentally and computationally derived ensembles exists, as simulations show overcompacted conformers. Increasing evidence suggests that the solvent plays a crucial role in shaping the ensembles of intrinsically disordered proteins and has led to several attempts to modify water parameters and thereby favor protein-water over protein-protein interactions. This study tackles the problem from a different perspective, which is the use of the Kirkwood-Buff theory of solutions to reproduce the correct conformational ensemble of intrinsically disordered proteins (IDPs). A protein force field recently developed on such a basis was found to be highly effective in reproducing ensembles for a fragment from the FG-rich nucleoporin 153, with dimensions matching experimental values obtained from small-angle X-ray scattering and single molecule FRET experiments. Kirkwood-Buff theory presents a complementary and fundamentally different approach to the recently developed four-site TIP4P-D water model, both of which can rescue the overcollapse observed in IDPs with canonical protein force fields. As such, our study provides a new route for tackling the deficiencies of current protein force fields in describing protein solvation.

Experimental triplet and quadruplet fluctuation densities and spatial distribution function integrals for liquid mixtures

Kirkwood-Buff or Fluctuation Solution Theory can be used to provide experimental pair fluctuations, and/or integrals over the pair distribution functions, from experimental thermodynamic data on liquid mixtures. Here, this type of approach is used to provide triplet and quadruplet fluctuations, and the corresponding integrals over the triplet and quadruplet distribution functions, in a purely thermodynamic manner that avoids the use of structure factors. The approach is then applied to binary mixtures of water + methanol and benzene + methanol over the full composition range under ambient conditions. The observed correlations between the different species vary significantly with composition. The magnitude of the fluctuations and integrals appears to increase as the number of the most polar molecule involved in the fluctuation or integral also increases. A simple physical picture of the fluctuations is provided to help rationalize some of these variations.

Simple and complex disorder in binary mixtures with benzene as a common solvent

Substituting benzene for water in computer simulations of binary mixtures, allows one to study the various forms of disorder, without the complications often encountered in aqueous mixtures.

A comparative Kirkwood-Buff study of aqueous methanol solutions modeled by the CHARMM additive and Drude polarizable force fields

A comparative study on aqueous methanol solutions modeled by the CHARMM additive and Drude polarizable force fields was carried out by employing Kirkwood-Buff analysis. It was shown that both models reproduced the experimental Kirkwood-Buff integrals and excess coordination numbers adequately well over the entire concentration range. The Drude model showed significant improvement over the additive model in solution densities, partial molar volumes, excess molar volumes, concentration-dependent diffusion constants, and dielectric constants. However, the additive model performed somewhat better than the Drude model in reproducing the activity derivative, excess molar Gibbs energy, and excess molar enthalpy of mixing. This is due to the additive achieving a better balance among solute-solute, solute-solvent, and solvent-solvent interactions, indicating the potential for improvements in the Drude polarizable alcohol model.

Kirkwood-Buff analysis of aqueous N-methylacetamide and acetamide solutions modeled by the CHARMM additive and Drude polarizable force fields

Kirkwood-Buff analysis was performed on aqueous solutions of N-methylacetamide and acetamide using the Chemistry at HARvard Molecular Mechanics additive and Drude polarizable all-atom force fields. Comparison of a range of properties with experimental results, including Kirkwood-Buff integrals, excess coordination numbers, solution densities, partial molar values, molar enthalpy of mixing, showed both models to be well behaved at higher solute concentrations with the Drude model showing systematic improvement at lower solution concentrations. However, both models showed difficulties reproducing experimental activity derivatives and the excess Gibbs energy, with the Drude model performing slightly better. At the molecular level, the improved agreement of the Drude model at low solute concentrations is due to increased structure in the solute-solute and solute-solvent interactions. The present results indicate that the explicit inclusion of electronic polarization leads to improved modeling of dilute solutions even when those properties are not included as target data during force field optimization.