We have been developing force fields designed for the eventual simulation of peptides and proteins using the Kirkwood-Buff (KB) theory of solutions as a guide. KB theory provides exact information on the relative distributions for each species present in solution. This information can also be obtained from computer simulations. Hence, one can use KB theory to help test and modify the parameters commonly used in biomolecular studies. A series of small molecule force fields representative of the fragments found in peptides and proteins have been developed. Since this approach is guided by the KB theory, our results provide a reasonable balance in the interactions between self-association of solutes and solute solvation. Here, we present our progress to date. In addition, our investigations have provided a wealth of data concerning the properties of solution mixtures, which is also summarized. Specific examples of the properties of aromatic (benzene, phenol, p-cresol) and sulfur compounds (methanethiol, dimethylsulfide, dimethyldisulfide) and their mixtures with methanol or toluene are provided as an illustration of this kind of approach.

No abstract

Fluctuation solution theory has provided an alternative view of many liquid mixture properties in terms of particle number fluctuations. The particle number fluctuations can also be related to integrals of the corresponding two body distribution functions between molecular pairs in order to provide a more physical picture of solution behavior and molecule affinities. Here, we extend this type of approach to provide expressions for higher order triplet and quadruplet fluctuations, and thereby integrals over the corresponding distribution functions, all of which can be obtained from available experimental thermodynamic data. The fluctuations and integrals are then determined using the International Association for the Properties of Water and Steam Formulation 1995 (IAPWS-95) equation of state for the liquid phase of pure water. The results indicate small, but significant, deviations from a Gaussian distribution for the molecules in this system. The pressure and temperature dependence of the fluctuations and integrals, as well as the limiting behavior as one approaches both the triple point and the critical point, are also examined. C 2015 AIP Publishing LLC. [http://dx

An extension of the traditional Kirkwood-Buff (KB) theory of solutions is outlined which provides additional fluctuating quantities that can be used to characterize and probe the behavior of solution mixtures. Particle-energy and energy-energy fluctuations for local regions of any multicomponent solution are expressed in terms of experimentally obtainable quantities, thereby supplementing the usual particle-particle fluctuations provided by the established KB inversion approach. The expressions are then used to analyze experimental data for pure water over a range of temperatures and pressures, a variety of pure liquids, and three binary solution mixtures - methanol and water, benzene and methanol, and aqueous sodium chloride. In addition to providing information on local properties of solutions it is argued that the particle-energy and energy-energy fluctuations can also be used to test and refine solute and solvent force fields for use in computer simulation studies.

In a continuation of our efforts to develop a united atom non-polarizable protein force field based upon the solution theory of Kirkwood and Buff i.e., the Kirkwood-Buff Force Field (KBFF) approach, we present KBFF models for the side chains of phenylalanine, tyrosine, tryptophan, and histidine, including both tautomers of neutral histidine and doubly-protonated histidine. The force fields were specifically designed to reproduce the thermodynamic properties of mixtures over the full composition range in an attempt to provide an improved description of intermolecular interactions. The models were developed by careful parameterization of the solution phase partial charges to reproduce the experimental Kirkwood-Buff integrals for mixtures of solutes representative of the amino acid sidechains in solution. The KBFF parameters and simulated thermodynamic and structural properties are presented for the following eleven binary mixtures: benzene + methanol, benzene + toluene, toluene + methanol, toluene + phenol, toluene + p-cresol, pyrrole + methanol, indole + methanol, pyridine + methanol, pyridine + water, histidine + water, and histidine hydrochloride + water. It is argued that the present approach and models provide a reasonable description of intermolecular interactions which ensures that the required balance between solute-solute, solute-solvent, and solvent-solvent distributions is obtained.

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.

Recent advances in fluctuation solution theory (FST) have provided access to information concerning triplet fluctuations and integrals, in addition to the established pair fluctuations and integrals, for liquids and liquid mixtures using both experimental and simulation data. Here, FST is used to investigate pair and triplet correlations for (i) pure water as provided by experiment and simulation using both polarizable and nonpolarizable water models, (ii) liquid mixtures of methanol and water as provided by experiment and simulation, and (iii) native and denatured states of proteins as provided by simulation. The last application is particularly powerful, as it provides exact equations for the volume, enthalpy, compressibility, thermal expansion, and heat capacity of a single protein form provided by a single simulation. In addition, a discussion of the quality of the integrals obtained from experiment and simulation is provided. The results clearly illustrate that FST can be a powerful tool for the analysis and interpretation of both experimental and simulation data in complex liquid mixtures, including biomolecular systems, and that current simulation protocols can provide reliable values for the pair and triplet correlations and integrals.

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. C 2015 AIP Publishing LLC. [http://dx