Equilibrium dialysis data and the relationships between preferential interaction parameters for biological systems in terms of Kirkwood-Buff integrals

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.

Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences

Interplay between divalent cations (Mg2+ and Ca2+) and single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA), as well as stacking interactions, is important in nucleosome stability and phase separation in nucleic acids. Quantitative techniques accounting for ion-DNA interactions are needed to obtain insights into these and related problems. Toward this end, we created a sequence-dependent computational TIS-ION model that explicitly accounts for monovalent and divalent ions. Simulations of the rigid 24 base-pair (bp) dsDNA and flexible ssDNA sequences, dT30 and dA30, with varying amounts of the divalent cations show that the calculated excess number of ions around the dsDNA and ssDNA agree quantitatively with ion-counting experiments. Using an ensemble of all-atom structures generated from coarse-grained simulations, we calculated the small-angle X-ray scattering profiles, which are in excellent agreement with experiments. Although ion-counting experiments mask the differences between Mg2+ and Ca2+, we find that Mg2+ binds to the minor grooves and phosphate groups, whereas Ca2+ binds specifically to the minor groove. Both Mg2+ and Ca2+ exhibit a tendency to bind to the minor groove of DNA as opposed to the major groove. The dA30 conformations are dominated by stacking interactions, resulting in structures with considerable helical order. The near cancellation of the favorable stacking and unfavorable electrostatic interactions leads to dT30 populating an ensemble of heterogeneous conformations. The successful applications of the TIS-ION model are poised to confront many problems in DNA biophysics.

Influence of salt and temperature on the self-assembly of cyclic peptides in water: a molecular dynamics study

It is found in the literature that cyclic peptides (CPs) are able to self-assemble in water to form cyclic peptide nanotubes (CPNTs) and are used extensively in the field of nanotechnology. Several factors influence the formation and stability of these nanotubes in water. However, an extensive study of the contribution of several important factors is still lacking. The purpose of this study is to explore the effect of temperature and salt (NaCl) on the association tendency of CPs. Furthermore, the self-association behavior of CPs in aqueous solutions at various temperatures is also thoroughly discussed. Cyclo-[(Asp-D-Leu-Lys-D-Leu)₂] is considered for this study and a series of classical molecular dynamics (MD) simulations at three different temperatures, viz. 280 K, 300 K, and 320 K, both in pure water and in NaCl solutions of different concentrations are carried out. The calculations of radial distribution functions, preferential interaction parameters, cluster formation and hydrogen bonding properties suggest a strong influence of NaCl concentration on the association propensity of CPs. Low NaCl concentration hinders CP association while high NaCl concentration facilitates the association of CPs. Besides this, the association of CPs is found to be enhanced at low temperature. Furthermore, the thermodynamics of CP association is predominantly found to be enthalpy driven in both the presence and absence of salt. No crossover between enthalpy and entropy in CP association is observed. In addition, the MM-GBSA method is used to investigate the binding free energies of the CP rings that self-assembled to form nanotube like structures at all three temperatures.

The chaotropic effect of ions on the self-aggregating propensity of Whitlock's molecular tweezers

Molecular tweezers feature the first class of artificial receptors to pique the interest of researchers and emerge as an effective therapeutic candidate. The exceptional structure and exquisite binding specificity of tweezers establish this overall class of receptors as a promising tool, with abundant applications. However, their inclination to self-aggregate by mutual π-π stacking interactions of their aromatic arms diminishes their efficacy as a therapeutic candidate. Therefore, following up on sporadic studies, since the discovery of the Hofmeister series, on the ability of ions to either solvate (salting-in) or induce aggregation (salting-out) of hydrophobic solutes, the notions of ion-specificity effects are utilized on tweezer moieties. The impacts of three different aluminum salts bearing anions Cl^-, ClO₄^- and SCN^- on the self-association propensity of Whitlock's caffeine-pincered molecular tweezers are investigated, with a specific emphasis placed on elucidating the varied behavior of the ions on the hydration ability of tweezers. The comparative investigation is conducted employing a series of all-atom molecular dynamics simulations of five tweezer molecules in pure water and three salt solutions, at two different concentrations each, maintaining a temperature of 300 K and a pressure of 1 atm, respectively. Radial distribution functions, coordination numbers, and SASA calculations display a steady reduction in the aggregation proclivity of the receptor molecules with an increase in salt concentration, as progressed along the Hofmeister series. Orientational preferences between the tweezer arms reveal a disruptive effect in the regular π-π stacking interactions, in the presence of high concentrations of ClO₄^- and SCN^- ions, while preferential interactions and tetrahedral order parameters unveil the underlying mechanism, by which the anions alter the solubility of the hydrophobic molecules. Overall, it is observed that SCN^- exhibits the highest salting-in effect, followed by ClO₄^-, with both anions inhibiting tweezer aggregation through different mechanisms. ClO₄^- ions impart an effect by moderately interacting with the solute molecules as well as modifying the water structure of the bulk solution promoting solvation, whereas, SCN^- ions engage entirely in interaction with specific tweezer sites. Cl^- being the most charge-dense of the three anionic species experiences stronger hydration and therefore, imparts a very negligible salting-in effect.

Accelerating the 3D reference interaction site model theory of molecular solvation with treecode summation and cut-offs

The 3D reference interaction site model (3D-RISM) of molecular solvation is a powerful tool for computing the equilibrium thermodynamics and density distributions of solvents, such as water and co-ions, around solute molecules. However, 3D-RISM solutions can be expensive to calculate, especially for proteins and other large molecules where calculating the potential energy between solute and solvent requires more than half the computation time. To address this problem, we have developed and implemented treecode summation for long-range interactions and analytically corrected cut-offs for short-range interactions to accelerate the potential energy and long-range asymptotics calculations in non-periodic 3D-RISM in the AmberTools molecular modeling suite. For the largest single protein considered in this work, tubulin, the total computation time was reduced by a factor of 4. In addition, parallel calculations with these new methods scale almost linearly and the iterative solver remains the largest impediment to parallel scaling. To demonstrate the utility of our approach for large systems, we used 3D-RISM to calculate the solvation thermodynamics and density distribution of 7-ring microtubule, consisting of 910 tubulin dimers, over 1.2 million atoms.

Phase separation property of a hydrophobic deep eutectic solvent-water binary mixture: A molecular dynamics simulation study

Over the past decade, deep eutectic solvents (DESs) have earned applicability in numerous fields as non-flammable, non-volatile, and greener alternatives to conventional organic solvents. In a first of its kind, a hydrophobic DES composed of a 1:1 mixture of oleic acid and lidocaine was recently reported, possessing a lower critical solution temperature in water. The thermoreversible phase property of this DES-water system was utilized to sequester out dye molecules from their aqueous solutions. In this article, we explore the phase separation phenomena for this particular DES in its aqueous solution using an all-atom molecular dynamics simulation. A 50 wt. % solution of the DES in water was studied at three different temperatures (253, 293, and 313 K) to understand the various molecular interactions that dictate the phase segregation property of these systems. In this work, we have elaborated on the importance of hydrogen bonding interactions and the non-bonding interactions between the components and the competition between the two that leads to phase separation. Overall, we observe that the increase in unfavorable interaction between the DES components and water with increasing temperature determines the phase separation behavior. We have also studied the modification in the dynamical properties of water molecules close to the phase boundary. Such molecular insights would be beneficial for designing novel solvent systems that can be used as extraction-based media in industries.

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

Here, we perform structural, thermodynamic, and kinetics tests of the Kirkwood-Buff-derived force field, KBFF20, for peptides and proteins developed in the previous article. The physical/structural tests measure the ability of KBFF20 to capture the experimental J-couplings for small peptides, to keep globular monomeric and oligomeric proteins folded, and to produce the experimentally relevant expanded conformational ensembles of intrinsically disordered proteins. The thermodynamic-based tests probe KBFF20's ability to quantify the preferential interactions of sodium chloride around native β-lactoglobulin and urea around native lysozyme, to reproduce the melting curves for small helix- and sheet-based peptides, and to fold the small proteins Trp-cage and Villin. The kinetics-based tests quantify how well KBFF20 can match the experimental contact formation rates of small, repeat-sequence peptides of variable lengths and the rotational diffusion coefficients of globular proteins. The results suggest that KBFF20 is naturally able to reproduce properties of both folded and disordered proteins, which we attribute to the use of the Kirkwood-Buff theory as the foundation of the force field's development. However, we show that KBFF20 tends to lose some well-defined secondary structural elements and increases the percentage of coil regions, indicating that the perfect balance of all interactions remains elusive. Nevertheless, we argue that KBFF20 is an improvement over recently modified force fields that require ad hoc interventions to prevent the collapse of intrinsically disordered proteins.

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.

Protein Preferential Solvation in Water:Glycerol Mixtures

For proteins in solvent mixtures, the relative abundances of each solvent in their solvation shell have a critical impact on their properties. Preferential solvation of a series of proteins in water-glycerol mixtures is studied here over a broad range of solvent compositions via classical molecular dynamics simulations. Our simulation results reveal that the differences between shell and bulk compositions exhibit dramatic changes with solvent composition, temperature, and protein nature. In contrast with the simple and widely used picture where glycerol is completely excluded from the protein interface, we show that for aqueous solutions with less than 50% glycerol in volume, protein solvation shells have approximately the same composition as the bulk solvent and proteins are in direct contact with glycerol. We further demonstrate that at high glycerol concentration, glycerol depletion from the solvation shell is due to an entropic factor arising from the reduced accessibility of bulky glycerol molecules in protein cavities. The resulting molecular picture is important to understand protein activity and cryopreservation in mixed aqueous solvents.

Theory and simulations for RNA folding in mixtures of monovalent and divalent cations

RNA molecules cannot fold in the absence of counterions. Experiments are typically performed in the presence of monovalent and divalent cations. How to treat the impact of a solution containing a mixture of both ion types on RNA folding has remained a challenging problem for decades. By exploiting the large concentration difference between divalent and monovalent ions used in experiments, we develop a theory based on the reference interaction site model (RISM), which allows us to treat divalent cations explicitly while keeping the implicit screening effect due to monovalent ions. Our theory captures both the inner shell and outer shell coordination of divalent cations to phosphate groups, which we demonstrate is crucial for an accurate calculation of RNA folding thermodynamics. The RISM theory for ion-phosphate interactions when combined with simulations based on a transferable coarse-grained model allows us to predict accurately the folding of several RNA molecules in a mixture containing monovalent and divalent ions. The calculated folding free energies and ion-preferential coefficients for RNA molecules (pseudoknots, a fragment of the rRNA, and the aptamer domain of the adenine riboswitch) are in excellent agreement with experiments over a wide range of monovalent and divalent ion concentrations. Because the theory is general, it can be readily used to investigate ion and sequence effects on DNA properties.

Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force Fields

Many applications in chemistry, biology, and energy storage/conversion research rely on molecular simulations to provide fundamental insight into structural and transport properties of materials with high ionic concentrations. Whether the system is comprised entirely of ions, like ionic liquids, or is a mixture of a polar solvent with a salt, e.g., liquid electrolytes for battery applications, the presence of ions in these materials results in strong local electric fields polarizing solvent molecules and large ions. To predict properties of such systems from molecular simulations often requires either explicit or mean-field inclusion of the influence of polarization on electrostatic interactions. In this manuscript, we review the pros and cons of different treatments of polarization ranging from the mean-field approaches to the most popular explicit polarization models in molecular dynamics simulations of ionic materials. For each method, we discuss their advantages and disadvantages and emphasize key assumptions as well as their adjustable parameters. Strategies for the development of polarizable models are presented with a specific focus on extracting atomic polarizabilities. Finally, we compare simulations using polarizable and nonpolarizable models for several classes of ionic systems, discussing the underlying physics that each approach includes or ignores, implications for implementation and computational efficiency, and the accuracy of properties predicted by these methods compared to experiments.

Understanding of Structure and Thermodynamics of Melamine Association in Aqueous Solution from a Unified Theoretical and Experimental Approach

The aggregation propensity of melamine molecules in aqueous solutions in a range of melamine concentrations is investigated by means of a combination of theoretical and experimental approaches. It is observed that the hydrogen bonding interactions of sp³ nitrogen atoms of one melamine with sp² nitrogen atoms of another melamine play a major role in the melamine association. This finding is complemented by the observed favorable electrostatic energies between melamine molecules. The estimation of the orientational probability of melamine aromatic ring rules out any role of π-π interaction in melamine association. Further, the quantum chemical calculations suggest that a melamine molecule prefers to bind with another like molecule with a dihedral angle ranging from 36° to 46°. We have also determined the dimer existence autocorrelation functions to investigate the melamine-dimer stability with time in aqueous solution. Our results are well validated by the experimental findings (Chapman, R. P.; Averell, P. R.; Harris, R. R. Solubility of Melamine in Water. Ind. Eng. Chem. 1943, 35, 137-138. Ahromi, A. J.; Moosheimer, U. Oxygen Barrier Coatings Based on Supramolecular Assembly of Melamine. Macromolecules 2000, 33, 7582-7587. Yang, C.; Liu. Y. Studying on the Steady-State and Time-Resolved Fluorescence Characteristics of Melamine. Spectrochim. Acta A 2010, 75, 1329-1332. Mircescu, N. E.; Oltean, M.; Chis, V.; Leopold, N. FTIR, FT-Raman, SERS and DFT study on Melamine. Vib. Spectrosc. 2012, 62, 165-171. Makowski. S. J.; Lacher. M.; Schnick. W. Supramolecular Hydrogenbonded Structures between Melamine and N-Heterocycles. J. Mol. Struct. 2012, 1013, 19-25. Li, Z.; Chen, G.; Xu, Y.; Wang, X.; Wang, Z. Study of the Structural and Spectral Characteristics of C₃N₃(NH₂)₃(n = 1-4) Clusters. J. Phys. Chem. A 2013, 117, 12511-12518. Li, P.; Arman, D. H.; Wang, H.; Weng, L.; Alfooty, K.; Angawi, R. F.; Chen. B. Solvent Dependent Structures of Melamine: Porous or Non-porous. Cryst. Growth Des. 2015, 15, 1871-1875. Chen, J.; Lei, X.; Peng, B. Study on the Fluorescence Spectra of Melamine in Pure Milk. J. Opt. 2017, 46, 183-186.). Moreover, the thermodynamics of melamine association reveals that the association process is essentially driven by enthalpy, and this enthalpy-driven phenomenon is also confirmed by the experimental isothermal titration calorimetry measurements.

Ion hydration: linking self-diffusion and reorientational motion to water structure

A link between water dynamics and the “water structure” has been established through the combination of the extended jump model, transition state theory and the Kirkwood-Buff theory.

Kirkwood-Buff integrals for hard-core Yukawa fluids

The Kirkwood-Buff (KB) theory of solution is employed to investigate several macroscopic properties of the one-component hard-core Yukawa (HCY) fluid, where the key physical quantities are the KB integrals (KBIs). For both repulsive and attractive HCY fluids, the radial distribution functions are calculated by using the classical density functional theory, and then the corresponding KBIs are carried out. Since the local structure and global properties of a fluid can be related by KBI, we presented the isothermal compressibility and the derivative of the chemical potential with respect to bulk density for both repulsive and attractive HCY fluids. It is found that a transition of the affinity of particles in an attractive HCY fluid exists. The corresponding phase diagrams on the affinity are illustrated, which consist of repulsive and attractive regions with the boundary line of KBIs being zero. These results show that the aggregated structure of a HCY fluid can be effectively regulated by the screening parameter, bulk density and interaction energy, while KBIs can provide a quantitative reliable description on the properties of HCY fluids.

How Osmolytes Counteract Pressure Denaturation on a Molecular Scale

Life in the deep sea exposes enzymes to high hydrostatic pressure, which decreases their stability. For survival, deep sea organisms tend to accumulate various osmolytes, most notably trimethylamine N-oxide used by fish, to counteract pressure denaturation. However, exactly how these osmolytes work remains unclear. Here, a rigorous statistical thermodynamics approach is used to clarify the mechanism of osmoprotection. It is shown that the weak, nonspecific, and dynamic interactions of water and osmolytes with proteins can be characterized only statistically, and that the competition between protein-osmolyte and protein-water interactions is crucial in determining conformational stability. Osmoprotection is driven by a stronger exclusion of osmolytes from the denatured protein than from the native conformation, and water distribution has no significant effect on these changes at low osmolyte concentrations.

Ammonium based stabilizers effectively counteract urea-induced denaturation in a small protein: insights from molecular dynamics simulations

Room temperature ionic liquids (IL) and deep eutectic solvents (DES) are known to aid the conformational stability and activity of proteins and enzymes in aqueous solutions.

Competitive interaction of monovalent cations with DNA from 3D-RISM

The composition of the ion atmosphere surrounding nucleic acids affects their folding, condensation and binding to other molecules. It is thus of fundamental importance to gain predictive insight into the formation of the ion atmosphere and thermodynamic consequences when varying ionic conditions. An early step toward this goal is to benchmark computational models against quantitative experimental measurements. Herein, we test the ability of the three dimensional reference interaction site model (3D-RISM) to reproduce preferential interaction parameters determined from ion counting (IC) experiments for mixed alkali chlorides and dsDNA. Calculations agree well with experiment with slight deviations for salt concentrations >200 mM and capture the observed trend where the extent of cation accumulation around the DNA varies inversely with its ionic size. Ion distributions indicate that the smaller, more competitive cations accumulate to a greater extent near the phosphoryl groups, penetrating deeper into the grooves. In accord with experiment, calculated IC profiles do not vary with sequence, although the predicted ion distributions in the grooves are sequence and ion size dependent. Calculations on other nucleic acid conformations predict that the variation in linear charge density has a minor effect on the extent of cation competition.

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.

Synergy in protein-osmolyte mixtures

Virtually all taxa use osmolytes to protect cells against biochemical stress. Osmolytes often occur in mixtures, such as the classical combination of urea with TMAO (trimethylamine N-oxide) in cartilaginous fish or the cocktail of at least six different osmolytes in the kidney. The concentration patterns of osmolyte mixtures found in vivo make it likely that synergy between them plays an important role. Using statistical mechanical n-component Kirkwood-Buff theory, we show from first principles that synergy in protein-osmolyte systems can arise from two separable sources: (1) mutual alteration of protein surface solvation and (2) effects mediated through bulk osmolyte chemical activities. We illustrate both effects in a four-component system with the experimental example of the unfolding of a notch ankyrin domain in urea-TMAO mixtures, which make urea a less effective denaturant and TMAO a more effective stabilizer. Protein surface effects are primarily responsible for this synergy. The specific patterns of surface solvation point to denatured state expansion as the main factor, as opposed to direct competition.

Hydrotrope accumulation around the drug: the driving force for solubilization and minimum hydrotrope concentration for nicotinamide and urea

A rigorous statistical thermodynamic theory explains how urea and nicotinamide can solubilize hydrophobic drugs in water.

Hydrotropy: monomer-micelle equilibrium and minimum hydrotrope concentration

Drug molecules with low aqueous solubility can be solubilized by a class of cosolvents, known as hydrotropes. Their action has often been explained by an analogy with micelle formation, which exhibits critical micelle concentration (CMC). Indeed, hydrotropes also exhibit "minimum hydrotrope concentration" (MHC), a threshold concentration for solubilization. However, MHC is observed even for nonaggregating monomeric hydrotropes (such as urea); this raises questions over the validity of this analogy. Here we clarify the effect of micellization on hydrotropy, as well as the origin of MHC when micellization is not accompanied. On the basis of the rigorous Kirkwood-Buff (KB) theory of solutions, we show that (i) micellar hydrotropy is explained also from preferential drug-hydrotrope interaction; (ii) yet micelle formation reduces solubilization effeciency per hydrotrope molecule; (iii) MHC is caused by hydrotrope-hydrotrope self-association induced by the solute (drug) molecule; and (iv) MHC is prevented by hydrotrope self-aggregation in the bulk solution. We thus need a departure from the traditional view; the structure of hydrotrope-water mixture around the drug molecule, not the structure of the aqueous hydrotrope solutions in the bulk phase, is the true key toward understanding the origin of MHC.

Ion counting from explicit-solvent simulations and 3D-RISM

The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20-25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na(+) binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.

Mechanism of hydrophobic drug solubilization by small molecule hydrotropes

Drugs that are poorly soluble in water can be solubilized by the addition of hydrotropes. Albeit known for almost a century, how they work at a molecular basis is still controversial due to the lack of a rigorous theoretical basis. To clear up this situation, a combination of experimental data and Fluctuation Theory of Solutions (FTS) has been employed; information on the interactions between all the molecular species present in the solution has been evaluated directly. FTS has identified two major factors of hydrotrope-induced solubilization: preferential hydrotrope-solute interaction and water activity depression. The former is dominated by hydrotrope-solute association, and the latter is enhanced by ionic dissociation and hindered by the self-aggregation of the hydrotropes. Moreover, in stark contrast to previous hypotheses, neither the change of solute hydration nor the water structure accounts for hydrotropy. Indeed, the rigorous FTS poses serious doubts over the other common hypothesis: self-aggregation of the hydrotrope hinders, rather than promotes, solubilization.

Theory and Simulation of Multicomponent Osmotic Systems

Most cellular processes occur in systems containing a variety of components many of which are open to material exchange. However, computer simulations of biological systems are almost exclusively performed in systems closed to material exchange. In principle, the behavior of biomolecules in open and closed systems will be different. Here, we provide a rigorous framework for the analysis of experimental and simulation data concerning open and closed multicomponent systems using the Kirkwood-Buff (KB) theory of solutions. The results are illustrated using computer simulations for various concentrations of the solutes Gly, Gly(2) and Gly(3) in both open and closed systems, and in the absence or presence of NaCl as a cosolvent. In addition, KB theory is used to help rationalize the aggregation properties of the solutes. Here one observes that the picture of solute association described by the KB integrals, which are directly related to the solution thermodynamics, and that provided by more physical clustering approaches are different. It is argued that the combination of KB theory and simulation data provides a simple and powerful tool for the analysis of complex multicomponent open and closed systems.

Fluctuation theory of molecular association and conformational equilibria

General expressions relating the effects of pressure, temperature, and composition on solute association and conformational equilibria using the fluctuation theory of solutions are provided. The expressions are exact and can be used to interpret experimental or computer simulation data for any multicomponent mixture involving molecules of any size and character at any composition. The relationships involve particle-particle, particle-energy, and energy-energy correlations within local regions in the vicinity of each species involved in the equilibrium. In particular, it is demonstrated that the results can be used to study peptide and protein association or aggregation, protein denaturation, and protein-ligand binding. Exactly how the relevant fluctuating properties may be obtained from experimental or computer simulation data are also outlined. It is shown that the enthalpy, heat capacity, and compressibility differences associated with the equilibrium process can, in principle, be obtained from a single simulation. Fluctuation based expressions for partial molar heat capacities, thermal expansions, and isothermal compressibilities are also provided.

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.

The effect of urea on the morphology of NaCl crystals: A combined theoretical and simulation study

It has been known for over a century that the presence of cosolvents such as urea and formamide can alter the morphology of NaCl crystals grown from solution. To help understand this effect we have been developing a theoretical approach based on the Kirkwood-Buff (KB) theory of solutions, and have combined this with computer simulations of the interation of urea with different crystal faces of NaCl. In this way one can predict the effect of urea on the thermodynamic stability of different NaCl faces, with atomic level detail provided by the simulations. We observe that urea is preferentially excluded from 100 and 111 crystal faces, but is less excluded from 111 faces which present chloride ions at the surface. The results indicate that the 111 face is stabilized in urea solutions and promotes the formation of octahedral over cubic NaCl crystals. The approach is totally general and can be applied to understand a variety of interfacial properties. Furthermore, we apply KB theory to study several other issues regarding the simulation of crystal growth.

Kirkwood-Buff theory of four and higher component mixtures

Explicit expressions are developed for the chemical potential derivatives, partial molar volumes, and isothermal compressibility of solution mixtures involving four components at finite concentrations using the Kirkwood-Buff theory of solutions. In addition, a general recursion relationship is provided which can be used to generate the chemical potential derivatives for higher component solutions.

On the theory of solute solubility in mixed solvents

A series of equations are developed for the study of the effects of cosolvents on the solubility of a solute in mixed solutions where the solute displays a finite solubility. The equations differ depending on the scale used for the solute (and cosolvent) concentrations. The expressions use Kirkwood-Buff integrals to relate the changes in solubility to changes in the local solution composition around the solute and can be applied to study any type of ternary system including electrolyte cosolvents. The expressions provided here differ from previous approaches because of the use of a semi-open ensemble and the extension to finite solute solubilities.

Recent applications of Kirkwood-Buff theory to biological systems

The effect of cosolvents on biomolecular equilibria has traditionally been rationalized using simple binding models. More recently, a renewed interest in the use of Kirkwood-Buff (KB) theory to analyze solution mixtures has provided new information on the effects of osmolytes and denaturants and their interactions with biomolecules. Here we review the status of KB theory as applied to biological systems. In particular, the existing models of denaturation are analyzed in terms of KB theory, and the use of KB theory to interpret computer simulation data for these systems is discussed.

An analysis of the molecular origin of osmolyte-dependent protein stability

Protein solvation is the key determinant for isothermal, concentration-dependent effects on protein equilibria, such as folding. The required solvation information can be extracted from experimental thermodynamic data using Kirkwood-Buff theory. Here we derive and discuss general properties of proteins and osmolytes that are pertinent to their biochemical behavior. We find that hydration depends very little on osmolyte concentration and type. Strong dependencies on both osmolyte concentration and type are found for osmolyte self-solvation and protein-osmolyte solvation changes upon unfolding. However, solvation in osmolyte solutions does not involve complex concentration dependencies as found in organic molecules that are not used as osmolytes in nature. It is argued that the simple solvation behavior of naturally occurring osmolytes is a prerequisite for their usefulness in osmotic regulation in vivo.

The temperature dependence of salt-protein association is sequence specific

Molecular dynamics (MD) simulations are used to probe the origin of the unexpected temperature dependence of salt accumulation in the C-terminal region of the protein human lymphotactin. As in previous MD simulations, sodium ions accumulate in an enhanced manner near the C-terminal helix at the lower temperature, while the temperature dependence of chloride accumulation is much weaker and slightly positive. In a designed mutant in which all positively charged residues in the C-terminal helix are replaced with neutral polar groups (Ser), the unexpected temperature dependence of the sodium ions is no longer observed. Therefore, these simulations convincingly verified the previous hypothesis that the temperature dependence of ion-protein association is sensitive to the local sequence. This is explained qualitatively in terms of the entropy of association between charged species in solution. These findings have general implications for the interpretation of thermodynamic quantities associated with binding events where ion release is important, such as protein-DNA interactions.

Chemical potential derivatives and preferential interaction parameters in biological systems from Kirkwood-Buff theory

New expressions for chemical potential derivatives and preferential interaction parameters for ternary mixtures are derived for open, semiopen, and closed ensembles in terms of Kirkwood-Buff integrals, where all three components are present at finite concentrations. This is achieved using a simple approach that avoids the use of the general matrix formulation of Kirkwood-Buff theory. The resulting expressions provide a rigorous foundation for the analysis of experimental and simulation data. Using the results, a simple model is developed and used to investigate the possible effects of finite protein concentrations on the corresponding cosolvent dependent chemical potential and denaturation thermodynamics.