Empty liquid state and re-entrant phase behavior of the patchy colloids confined in porous media.

Patchy colloids with three and four equivalent patches, confined in an attractive random porous medium, undergo re-entrant gas-liquid phase separation with the liquid phase density approaching zero at low temperatures. The (bonding) colloid-colloid interaction causes the liquid-gas phase separation, which is modulated by the presence of the randomly distributed hard-sphere obstacles, attracting the colloids via Yukawa potential. Due to this interaction, a layer of mutually bonded colloids around the obstacles is formed. The network becomes nonuniform, with colloid particles locally centered on the obstacles. Features described in this article may open possibilities to produce equilibrium gels with predefined nonuniform distribution of particles and indicate how complicated the phase behavior of biological macromolecules in a crowded environment may be.

Modeling the depletion effect caused by an addition of polymer to monoclonal antibody solutions.

We present a theoretical study of colloidal stability of the model mixtures of monoclonal antibody molecules and non-adsorbing (no polymer-protein attraction) polymers. The antibodies are pictured as an assembly of seven hard spheres assuming a Y-like shape. Polymers present in the mixture are modeled as chain-like molecules having from 32 up to 128 monomers represented as hard spheres. We use Wertheim's thermodynamic perturbation theory to construct the two molecular species and to calculate measurable properties. The calculations are performed in the osmotic ensemble. In view that no direct attractive interaction is present in the model Hamiltonian, we only account for the entropic contribution to the phase equilibrium. We calculate chemical potentials and the equation of state for the model mixture to determine the liquid-liquid part of the phase diagram. We investigate how the critical antibody number density depends on the degree of polymerization and the bead size ratio of the polymer and protein components. The model mixture qualitatively correctly predicts some basic features of real systems. The effects of the model 'protein' geometry, that is the difference in results for the flexible Y-shaped protein versus the rigid spherical one, are also examined.

Model for a mixture of macroions, counterions, and co-ions in a waterlike fluid.

We propose an integral equation theory for a mixture of macroions, counterions, and co-ions in a waterlike fluid in which all the components are accounted for explicitly. The macroions can carry positive and negative surface charges simultaneously, mimicking in this way the situation occurring in protein solutions. To solve this complex model numerically, we utilize the associative mean spherical approximation, developed earlier for low-molecular-mass charge-symmetric electrolyte solutions. To illustrate the potential of this approach, we present numerical results for various experimental conditions. Among the measurable properties we choose to calculate the osmotic coefficient, a quantity that reflects the stability of the solution. We show that the osmotic coefficient depends not only on the magnitude of the net charge on the macroion but also on its sign, as well as on the nature of the low-molecular-mass electrolyte present. These specific ion effects are the consequence of differences in hydration between the ions in solution and charged groups on the macroion.

Computer simulations of ionenes, hydrophobic ions with unusual solution thermodynamic properties. The ion-specific effects.

Ionenes are alkyl polymer chains in which different numbers of methylene groups separate quaternary ammonium groups. They are ideal molecules for studying the balance between hydrophobic and charge effects in water. Implicit-solvent models predict osmotic coefficients that are too high (too low water vapor pressures), compared to experiments. We present a molecular dynamics simulation, in explicit SPC/E water, of a solution of aliphatic 6,6 ionene oligocations with sodium co-ions and fluorine, chlorine, bromine, or iodine counterions. In the 6,6 ionene solution, the latter polyion has more hydrophobic groups than its 3,3 counterpart, the waters are displaced more from the oligoion surface. Also, we find that the large ions, such as iodine, act like hydrophobic groups insofar as they bind to ionene's methylene groups. The water-mediated attraction between fluorine ions is enhanced in presence of weakly charged 6,6 ionene molecules. This effect may additionally reduce the osmotic pressure in such systems. Our results can explain some experimental trends in ionene solutions and weakly charged polyelectrolytes in general.

Explicit-water molecular dynamics study of a short-chain 3,3 ionene in solutions with sodium halides.

Ionenes are alkyl polymer chains in which hydrophobic groups are separated by ionic charges. They are useful for studying the properties of water as a solvent because they demonstrate a sufficiently complex combination of hydrophobicity, charge interactions, and specific-ion effects that some properties cannot be predicted by implicit-solvation theories. On the other hand, they are simple enough that their molecular structures can be varied and controlled in systematic experiments. In particular, implicit-solvent models predict that all such solutes will have negative enthalpies of dilution, whereas experiments show that enthalpies of dilution are positive for the chaotropic counterions. Here, we study ionenes that are short chains (six monomer units) in solutions of different counterions, with sodium as the coion by molecular dynamics simulations in explicit water. We explore the pair distributions of various atoms within the system at three different temperatures: T=278, 298, and 318 K. We find (i) that the molecular dynamics simulations are consistent with the experimental trends for the osmotic coefficients and enthalpies of dilution, (ii) that the fluorine-nitrogen and fluorine-carbon correlations decrease with decreasing temperature, (iii) while the opposite behavior is found for iodine ions, and (iv) that in the counterion-Na(+) pair distributions, too, fluorine ions behave oppositely to iodine ions upon temperature increase.

Ionic correlations in the inhomogeneous atmosphere surrounding cylindrical polyions: catalytic effects of polyions.

The structural properties of linear polyelectrolyte solutions in the presence of a salt as evidenced through ionic correlations in the inhomogeneous atmosphere around a polyion and their consequence such as the catalytic potential are studied by using Monte Carlo simulation techniques. The simulations are performed on the cylindrical cell model where a uniformly charged hard cylinder mimics the linear polyion, which is caged in its own cylindrical cell containing counterions and salt. The cell (volume) average of the interionic correlations is presented as a function of the polyion and salt concentrations and ion radius. These results are utilized to study the catalytic effects of polyions as manifested through the changes in the collision frequency between ions in the double layer surrounding the polyion relative to that in the pure electrolyte solution. The reported results suggest a strong influence of the added salt/polyelectrolyte concentration ratio on the structural properties of the solution and hence on ion-ion collision frequency. The machine simulations are supplemented by nonlinear Poisson-Boltzmann results. Fair agreement between two different theoretical methods of calculating the collision frequency is obtained.

Theory for the solvation of nonpolar solutes in water.

We recently developed an angle-dependent Wertheim integral equation theory (IET) of the Mercedes-Benz (MB) model of pure water [Silverstein et al., J. Am. Chem. Soc. 120, 3166 (1998)]. Our approach treats explicitly the coupled orientational constraints within water molecules. The analytical theory offers the advantage of being less computationally expensive than Monte Carlo simulations by two orders of magnitude. Here we apply the angle-dependent IET to studying the hydrophobic effect, the transfer of a nonpolar solute into MB water. We find that the theory reproduces the Monte Carlo results qualitatively for cold water and quantitatively for hot water.

An improved thermodynamic perturbation theory for Mercedes-Benz water.

We previously applied Wertheim's thermodynamic perturbation theory for associative fluids to the simple Mercedes-Benz model of water. We found that the theory reproduced well the physical properties of hot water, but was less successful in capturing the more structured hydrogen bonding that occurs in cold water. Here, we propose an improved version of the thermodynamic perturbation theory in which the effective density of the reference system is calculated self-consistently. The new theory is a significant improvement, giving good agreement with Monte Carlo simulations of the model, and predicting key anomalies of cold water, such as minima in the molar volume and large heat capacity, in addition to giving good agreement with the isothermal compressibility and thermal expansion coefficient.

Counterion-counterion correlation in the double layer around cylindrical polyions: counterion size and valency effects.

Monte Carlo simulation and Poisson-Boltzmann results on some aspects of structure and thermodynamics of aqueous polyelectrolyte solutions are presented. The polyelectrolyte solution is described by an infinitely long cylindrical polyion surrounded by counterions modeled as rigid ions moving in a continuum dielectric. Ion-ion correlations in the form of volume average of the counterion-counterion distribution function in the double layer surrounding the polyion are reported for mono- and divalent counterions and for a range of polyion concentrations and charge density parameters in each case. These results confirm again strong influence of the charge density parameter of polyions on properties of polyelectrolyte solutions. The structural information is supplemented by the calculated thermodynamic properties such as osmotic coefficients and heats of dilutions; the latter quantity has not been examined yet in detail by computer simulations. The results are discussed in view of the existing experimental data from the literature for these properties.

Computer simulations and theoretical aspects of the depletion interaction in protein-oligomer mixtures.

The depletion interaction between proteins caused by addition of either uncharged or partially charged oligomers was studied using the canonical Monte Carlo simulation technique and the integral equation theory. A protein molecule was modeled in two different ways: either as (i) a hard sphere of diameter 30.0 A with net charge 0, or +5, or (ii) as a hard sphere with discrete charges (depending on the pH of solution) of diameter 45.4 A. The oligomers were pictured as tangentially jointed, uncharged, or partially charged, hard spheres. The ions of a simple electrolyte present in solution were represented by charged hard spheres distributed in the dielectric continuum. In this study we were particularly interested in changes of the protein-protein pair-distribution function, caused by addition of the oligomer component. In agreement with previous studies we found that addition of a nonadsorbing oligomer reduces the phase stability of solution, which is reflected in the shape of the protein-protein pair-distribution function. The value of this function in protein-protein contact increases with increasing oligomer concentration, and is larger for charged oligomers. The range of the depletion interaction and its strength also depend on the length (number of monomer units) of the oligomer chain. The integral equation theory, based on the Wertheim Ornstein-Zernike approach applied in this study, was found to be in fair agreement with Monte Carlo results only for very short oligomers. The computer simulations for a model mimicking the lysozyme molecule (ii) are in qualitative agreement with small-angle neutron experiments for lysozyme-dextran mixtures.

Theoretical aspects and computer simulations of flexible charged oligomers in salt-free solutions.

The structural and thermodynamic properties of a model solution containing flexible charged oligomers and an equivalent number of counterions were studied by means of the canonical Monte Carlo simulation and integral equation theory. The oligomers were represented as freely jointed chains of charged hard spheres. In accordance with the primitive model of electrolyte solutions, the counterions were modeled as charged hard spheres and the solvent as a dielectric continuum. Simulations were performed for a set of model parameters, independently varying the chain length and concentration of the oligomers. Structural properties in the form of pair distribution functions were calculated as functions of model parameters. In addition, thermodynamic properties such as the excess energy of solution and the excess chemical potential of counterions were obtained. These properties were correlated with the conformational averages of oligomers as reflected in the end-to-end distances and radii of gyration obtained from the simulations. The relation with the experimental data for heats of dilution and for the activity coefficient is discussed. Finally, theories based on Wertheim's integral equation approach (product reactant Ornstein-Zernike approach) [J. Stat. Phys. 42, 477 (1986)] in the so-called polymer mean spherical and polymer hypernetted chain approximations were tested against the new and existing computer simulations. For the values of parameters examined in this study, the integral equation theory yields semiquantitative agreement with computer simulations.

Analysis of osmotic pressure data for aqueous protein solutions via a multicomponent model.

Integral equation theories and Monte Carlo simulations were used to study the Donnan equilibrium, which is established by an equilibrium distribution of a simple electrolyte between an aqueous protein-electrolyte mixture and an aqueous solution of the same simple electrolyte, when these two phases are separated by a semipermeable membrane. In order to describe the unusually low osmotic pressure found in many experiments we assumed that protein molecules can form dimers. The model solution contains proteins in a monomeric form, represented as charged hard spheres, or in a dimerized form, modeled as fused charged hard spheres. The counterions and coions were also modeled as charged hard spheres but of a much smaller size. The associative mean spherical and hypernetted-chain approximations were applied to this model. In addition, Monte Carlo computer simulations were performed for the same model system mimicking a lysozyme solution in the presence of 0.1 M sodium chloride. Theory and simulations were found to be in reasonably good agreement for the thermodynamic properties such as chemical potential and osmotic pressure under these conditions. Using the theoretical approaches mentioned above, we analyzed the experimental data for the osmotic pressure of bovine serum albumin in 0.15 M sodium chloride, human serum albumin solution (HSA) in 0.1 M phosphate buffer, and lysozyme in sulphate and phosphate buffers. The theoretically determined osmotic coefficients were fitted to the existing experimental data in order to obtain the fraction of dimers in solution. Our analysis indicated that there was relatively small self-association of protein molecules for bovine serum albumin solutions at pH=5.4 and 7.3, with the fraction of dimers smaller than 10%, while at pH=4.5 the dimer fraction was equal to 50%. In the case of HSA solutions, strong negative deviations from the ideal value were found and at pH=8.0 a reasonably good agreement between the theory and experiment is obtained by assuming full dimerization. For HSA solution at pH=5.4, the best fit to the experimental results was obtained for a fraction of dimers equal to 80%.

Confined water: a Mercedes-Benz model study.

We study water that is confined within small geometric spaces. We use the Mercedes-Benz (MB) model of water, in NVT and muVT Monte Carlo computer simulations. For MB water molecules between two planes separated by a distance d, we explore the structures, hydrogen bond networks, and thermodynamics as a function of d, temperature T, and water chemical potential mu. We find that squeezing the planes close enough together leads to a vaporization of waters out of the cavity. This vaporization transition has a corresponding peak in the heat capacity of the water. We also find that, in small pores, hydrogen bonding is not isotropic but, rather, it preferentially forms chains along the axis of the cavity. This may be relevant for fast proton transport in pores. Our simulations show oscillations in the forces between the inert plates, due to water structure, even for plate separations of 5-10 water diameters, consistent with experiments by Israelachvili et al. [Nature 1983, 306, 249]. Finally, we find that confinement affects water's heat capacity, consistent with recent experiments of Tombari et al. on Vycor nanopores [J. Chem. Phys. 2005, 122, 104712].

Theoretical study of catalytic effects in micellar solutions.

The catalytic effect of charged micelles as manifested through the increased collision frequency between the counterions of an electrolyte in the presence of such micelles is explored by the Monte Carlo simulation technique and various theoretical approaches. The micelles and ions are pictured as charged hard spheres embedded in a dielectric continuum with the properties of water at 298 K with the charge on micelles varying from zero to z(m) = 50 negative elementary charges. Analytical theories such as (i) the symmetric Poisson-Boltzmann theory, (ii) the modified Poisson-Boltzmann theory, and (iii) the hypernetted-chain integral equation are applied and tested against the Monte Carlo data for micellar ions (m) with up to 50 negative charges in aqueous solution with monovalent counterions (c; z(c) = +1) and co-ions (co; z(co) = -1). The results for the counterion-counterion pair correlation function at contact, g(cc)(sigma(cc)), are calculated in a micellar concentration range from c(m) = 5 x 10(-)(6) to 0.1 mol/dm(3) with an added +1:-1 electrolyte concentration of 0.005 mol/dm(3) (for most cases), and for various model parameters. Our computations indicate that even a small concentration of a highly charged polyelectrolyte added to a +1:-1 electrolyte solution strongly increases the probability of finding two counterions in contact. This result is in agreement with experimental data. For low charge on the micelles (z(m) below -8), all the theories are in qualitative agreement with the new computer simulations. For highly charged micelles, the theories either fail to converge (the hypernetted-chain theory) or, alternatively, yield poor agreement with computer data (the symmetric Poisson-Boltzmann and modified Poisson-Boltzmann theories). The nonlinear Poisson-Boltzmann cell model results yield reasonably good agreement with computer simulations for this system.

Clustering of macroions in solutions of highly asymmetric electrolytes.

In this paper, we present results of computer simulations for a primitive model of asymmetric electrolyte solutions containing macroions, counterions and in a few cases, also co-ions. The results show that the valency of counterions plays an important role in shaping the net interaction between the macroions. For solutions with monovalent counterions, the macroions are distributed at larger distances, and in solutions with divalent counterions, the macroions come closer to each other and share a layer of counterions, whereas, in solutions with trivalent counterions, the macroions form clusters. These clusters dissolve upon dilution or addition of a simple electrolyte. These findings suggest a mechanism whereby the nonuniform distribution of macroions observed experimentally in charged systems may occur.

Ionic effects beyond Poisson-Boltzmann theory.

Polyelectrolytes are electrolytes asymmetric both in charge and size. Their properties in solution are dominated by Coulombic forces, and without a detailed understanding of these interactions, no interpretation of experimental data is possible. This paper is a review of recent developments in the theory of highly asymmetric electrolytes of spherical shape resembling surfactant micelles. Three different models are discussed: (a) the cell model, which is focused on the small ion-macroion interaction; (b) the model that treats the solution as an effective one-component fluid of macroions; and (c) the isotropic model, where the solution is represented as a mixture of charged spheres. Traditionally, the electrostatic interactions are accounted for via the solution of the Poisson-Boltzmann equation. This theory, however, ignores the fluctuations around the most probable distribution and may yield poor results for systems with multivalent ions. This paper focuses on developments beyond the Poisson-Boltzmann theory; the results of computer simulations and integral equation theories represent the major part of the review.

Liquid-Vapor Coexistence in the Screened Coulomb (Yukawa) Hard Sphere Binary Mixture in Disordered Porous Media: The Mean Spherical Approximation.

The thermodynamics of a two-component fluid with a hard core interaction and screened Coulomb (Yukawa) interaction between particles, similar to the primitive model of an electrolyte solution, adsorbed in a disordered matrix of hard spheres, is studied by using replica Ornstein-Zernike integral equations and the mean spherical approximation (MSA). The gas-liquid transition is localized. The coexistence curve is investigated dependent on the range of interaction between fluid species, on matrix density, and on fluid-matrix attraction. We have observed shrinking of the coexistence envelope with increasing matrix density. The critical temperature of adsorbed mixture decreases with increasing matrix density. The critical density is less affected; however, it also decreases slightly. The critical temperature is sensitive to the fluid species-matrix attraction and depends nonmonotonously on their strength. For a given matrix microporosity, it increases slightly and then decreases with augmenting strength of fluid-matrix attraction. The critical density is less affected by this attraction. However, it decreases for the model with a sufficiently long-range tail of fluid-matrix attraction. Copyright 1998 Academic Press.

Distribution of ions around DNA, probed by energy transfer.

Measurements of the effect of DNA on rates of bimolecular energy transfer between ions provide a direct indication of how cations cluster in regions near DNA and how anions are repelled from the same regions. Energy transfer from luminescent lanthanide ions (in the "rapid-diffusion" limit) probes collision frequencies that are dependent on the equilibrium spatial distributions of ions. The addition of 1 mM DNA (phosphate) to a 2 mM salt solution increases the overall collision frequency between monovalent cations by a factor of 6 +/- 1.5; it increases the divalent-monovalent cation collision frequency by a factor of 29 +/- 3; and it decreases the divalent cation-monovalent anion collision frequency by a factor of 0.24 +/- 0.03. Comparisons are made with the changes in collision frequencies predicted by several different theoretical descriptions of ion distributions. The closest agreement with experimental results for monovalent ions at 1 mM DNA is obtained with a static accessibility-modified discrete charge calculation, based on a detailed molecular model of B-DNA. At high DNA concentration (10 mM), the best results are obtained by numerical solutions of the Poisson-Boltzmann equation for a "soft-rod" model of DNA. Poisson-Boltzmann calculations for a "hard-rod" model greatly overestimate the effects of DNA on collision frequencies, as does a calculation based on counterion-condensation theory.

Thermodynamics of the isothermal interaction of human immunoglobulin G with guanidinium chloride.

A thermodynamic study of the isothermal interaction of human immunoglobulin G with guanidinium chloride, a strong denaturant, has been performed. Free energies of interaction were calculated using preferential binding data obtained by measuring densities at constant chemical potential and constant composition, respectively. Enthalpies of interaction were determined calorimetrically. The values of both thermodynamic parameters as well as those of entropies of interaction have been found to depend crucially on the extent of denaturant binding.

Thermodynamics of the isothermal interaction of beta-lactoglobulin with guanidinium chloride and urea.

A thermodynamic study of the isothermal interaction of beta-lactoglobulin with guanidinium chloride and urea has been performed. Enthalpies of interaction of the two denaturants have been obtained by calorimetric measurements, and the free energy of interaction calculated from previously determined preferential binding of denaturants. In separate dilatometric experiments the volume changes accompanying the interaction of guanidinium chloride with beta-lactoglobulin have also been determined.