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.

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 (MX2), where M = Mg2+, Ca2+, Sr2+, Ba2+ 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.

A Kirkwood-Buff Derived Force Field for Aqueous Alkali Halides.

A classical nonpolarizable force field is presented for the simulation of aqueous alkali halide solutions (MX), where M = Li(+), Na(+), K(+), Rb(+), Cs(+) and X = F(-), Cl(-), Br(-), I(-), and their interactions with biomolecules. The models are specifically designed to reproduce the experimental Kirkwood-Buff integrals, and thereby the solution salt activities, as a function of salt concentration. Additionally, we demonstrate that these models reasonably reproduce other experimental properties including ion diffusion constants, dielectric decrements, and the excess heats of mixing. The parameters are developed by considering the properties of aqueous NaX and MCl solutions using a previously established model for NaCl. Transferability of the parameters to other salts is then established by the successful simulation of additional aqueous salt solutions, KI and CsBr, not originally included in the parameterization procedure.

Kirkwood-Buff integrals for ideal solutions.

The Kirkwood-Buff (KB) theory of solutions is a rigorous theory of solution mixtures which relates the molecular distributions between the solution components to the thermodynamic properties of the mixture. Ideal solutions represent a useful reference for understanding the properties of real solutions. Here, we derive expressions for the KB integrals, the central components of KB theory, in ideal solutions of any number of components corresponding to the three main concentration scales. The results are illustrated by use of molecular dynamics simulations for two binary solutions mixtures, benzene with toluene, and methanethiol with dimethylsulfide, which closely approach ideal behavior, and a binary mixture of benzene and methanol which is nonideal. Simulations of a quaternary mixture containing benzene, toluene, methanethiol, and dimethylsulfide suggest this system displays ideal behavior and that ideal behavior is not limited to mixtures containing a small number of components.

Developing Force Fields from the Microscopic Structure of Solutions.

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.

A Kirkwood-Buff derived force field for thiols, sulfides, and disulfides.

A force field has been developed for molecular simulations of methanethiol, dimethyl sulfide, and dimethyl disulfide mixtures. The force field specifically attempts to balance the solvation and self-association of these solutes in solution mixtures with methanol. The force field is based on the Kirkwood-Buff (KB) theory of solutions and is parametrized using the KB integrals obtained from the experimental activity coefficients for the solution mixtures. The transferability of the force field was tested and confirmed by the accurate prediction of the activity coefficients for methanethiol/dimethyl sulfide solutions, which were not used in the initial parametrization of the force fields. The ideality of this latter solution is excellently reproduced. The applicability of the force field to simulations in water was corroborated with a reasonably accurate prediction for the low solubility of dimethyl sulfide in water. The aggregation of methanol molecules at low methanol mole fractions displayed by all the mixtures is reproduced and further analyzed.

Droplet deformation in dc electric fields: the extended leaky dielectric model.

The leaky dielectric model (LDM) was extended to large droplet distortions in dc electric fields. The resulting extended LDM (ELDM) reduces to the LDM for small droplet aspect ratios and to the pure dielectric model when the ratio of droplet and matrix conductivities equals the inverse ratio of their permittivities. The ELDM distinguishes between two types of phenomena possible at high electric fields: continuous deformation and hysteresis. For droplets deforming parallel to the electric field, the relationship that distinguishes between the two phenomena is a function of the droplet and matrix conductivities and viscosities but not of their permittivities. For droplets deforming perpendicular to the electric field, the relationship is a function of the droplet permittivities and conductivities but depends only slightly on the ratio of their viscosities. Some of the predictions of the LDM and the ELDM were compared with our own data and with data from the literature. For the systems that deformed parallel to the field direction, the ELDM not only predicted the data qualitatively but also predicted the data quantitatively when the experimental errors in its input parameters were taken into account, whereas the older LDM did not even predict the qualitative trend of the data. For the systems that deformed perpendicular to the field direction, however, the ELDM predicted the observed the aspect ratios in only one out of the four systems examined. In the other three systems, the LDM appeared to give reasonable predictions when either the ratio of the matrix/droplet viscosities was relatively small or the value of total charge relaxation time was relatively large. Thus, the applicability of the ELDM, as presently formulated, appears to be limited in the case of deformations perpendicular to the electric field.