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

Ion-Ion Association in Bulk Mixed Electrolytes Using Global and Local Electroneutrality Constraints

Ion atmospheres play a critical role in modulating the interactions between charged components in solutions. However, a detailed description of the nature of ion atmospheres remains elusive. Here, we use Kirkwood-Buff theory, an exact theory of solution mixtures, together with a series of local and bulk electroneutrality constraints to provide relationships between all the net ion-ion distributions in bulk electrolyte mixtures. The validity of the underlying relationships is then confirmed using classical explicit solvent molecular simulations of a range of electrolyte mixtures. Further analysis indicates the ion distributions can be separated into two contributions, one resulting in charge neutralization, for which each ion contributes in proportion to its ionic strength, and the other accounting for all the solution thermodynamics. The relationships hold for atomic and molecular ions of any size and valency regardless of ionic strength, temperature, or pressure, in any solvent system.

Structural Study of Aqueous Electrolyte Solution by MeV Liquid Electron Scattering

The impact of ions on water has long been a subject of great interest, as it is closely tied to the hydration structure, dynamics, and properties of electrolyte solutions. Over centuries of investigation, the influence of ions on water's structure remains highly debated. Prevailing techniques, such as neutron and X-ray scattering, primarily focus on the microscopic structure of salt solutions at very high concentrations, mostly above 1 mol/L. In this study, we measured the structure of aqueous potassium iodide (KI) and potassium chloride (KCl) solutions using MeV liquid electron scattering (MeV-LES) across a concentration range of 0.10 to 0.75 mol/L. The obtained results provide detailed insights into the variations in ion-oxygen and oxygen-oxygen correlations as a function of concentration. The observed structural differences between KI and KCl solutions are in line with the structure maker/breaker theory, which suggests that iodide ions exert a more pronounced effect than chloride ions on disrupting the water shell. This work demonstrates the potency of MeV-LES for investigating the atomic structure in liquids, augmenting the modern analytical toolbox.

Preferential Solvation Phenomena as Solute-Induced Structure-Making/Breaking Processes: Linking Thermodynamic Preferential Interaction Parameters to Fundamental Structure Making/Breaking Functions

In this work, we identify the explicit macroscopic-to-microscopic rigorous links between existing thermodynamic preferential interaction parameters Γ Q α Q β ( χ i ) and microstructural descriptors based on total correlation function integrals, leading to their unambiguous characterization in terms of fundamental structure making/breaking functions S α β . First, we provide the statistics mechanical framework to identify a universal molecular-based signature for the preferential solvation P S phenomenon involving solutes at infinite dilution in mixed-solvent environments and discuss its fundamental properties. Then, we characterize the S α β functions relevant to the P S process, identify the microscopic markers for the existing preferential interaction parameters Γ Q α Q β ( χ i ) in terms of the S α β functions, and test their compliance with a pair of essential microstructural constraints linked to the properties of the universal P S signature. Moreover, we illustrate the analysis by probing the behavior of a representative ternary system comprising the solubility of methane in aqueous 1,4-dioxane mixed-solvent environments under ambient conditions. Finally, we discuss some relevant issues surrounding the statistical mechanical (microstructural) interpretation of the thermodynamic (macroscopic) preferential interaction parameters, review some pitfalls in their evaluation from molecular simulation, and provide an outlook.

Toward a general neural network force field for protein simulations: Refining the intramolecular interaction in protein

Molecular dynamics (MD) is an extremely powerful, highly effective, and widely used approach to understanding the nature of chemical processes in atomic details for proteins. The accuracy of results from MD simulations is highly dependent on force fields. Currently, molecular mechanical (MM) force fields are mainly utilized in MD simulations because of their low computational cost. Quantum mechanical (QM) calculation has high accuracy, but it is exceedingly time consuming for protein simulations. Machine learning (ML) provides the capability for generating accurate potential at the QM level without increasing much computational effort for specific systems that can be studied at the QM level. However, the construction of general machine learned force fields, needed for broad applications and large and complex systems, is still challenging. Here, general and transferable neural network (NN) force fields based on CHARMM force fields, named CHARMM-NN, are constructed for proteins by training NN models on 27 fragments partitioned from the residue-based systematic molecular fragmentation (rSMF) method. The NN for each fragment is based on atom types and uses new input features that are similar to MM inputs, including bonds, angles, dihedrals, and non-bonded terms, which enhance the compatibility of CHARMM-NN to MM MD and enable the implementation of CHARMM-NN force fields in different MD programs. While the main part of the energy of the protein is based on rSMF and NN, the nonbonded interactions between the fragments and with water are taken from the CHARMM force field through mechanical embedding. The validations of the method for dipeptides on geometric data, relative potential energies, and structural reorganization energies demonstrate that the CHARMM-NN local minima on the potential energy surface are very accurate approximations to QM, showing the success of CHARMM-NN for bonded interactions. However, the MD simulations on peptides and proteins indicate that more accurate methods to represent protein-water interactions in fragments and non-bonded interactions between fragments should be considered in the future improvement of CHARMM-NN, which can increase the accuracy of approximation beyond the current mechanical embedding QM/MM level.

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.

Protein assembly and crowding simulations

Proteins encounter frequent molecular interactions in biological environments. Computer simulations have become an increasingly important tool in providing mechanistic insights into how such interactions in vivo relate to their biological function. The review here focuses on simulations describing protein assembly and molecular crowding effects as two important aspects that are distinguished mainly by how specific and long-lived protein contacts are. On the topic of crowding, recent simulations have provided new insights into how crowding affects protein folding and stability, modulates enzyme activity, and affects diffusive properties. Recent studies of assembly processes focus on assembly pathways, especially for virus capsids, amyloid aggregation pathways, and the role of multivalent interactions leading to phase separation. Also, discussed are technical challenges in achieving increasingly realistic simulations of interactions in cellular environments.