McMillan–Mayer theory for solvent effects in inhomogeneous systems: Calculation of interaction pressure in aqueous electrical double layers

We demonstrate how to use the McMillan–Mayer theory to include solvent effects in effective solute–solute interactions for inhomogeneous systems, extending a recent derivation [S. Marčelja, Langmuir 16, 6081 (2000)] for symmetric planar double layers to the general case. In the exact treatment, the many-body potential of mean force between the solute molecules can be evaluated for an inhomogeneous reference system in equilibrium with pure bulk solvent. The reference system contains only solvent and a finite number, n, of fixed solute molecules and it has an external potential that in some cases is different from that of the original system. It is discussed how the n-body potential of mean force between the ions for the relevant cases of large n values can be approximated by a sum of effective singlet and pair interactions evaluated in the presence of, on average, all n ions, i.e., at finite concentration. In examples considered in this work we use effective interionic pair potentials evaluated from bulk electrolyte calculations at finite electrolyte concentrations. We calculate the contribution to the double layer interaction pressure arising from the interaction between ions dissolved in aqueous electrolyte. In cases of moderate or high surface charge, calculations show several new effects. At small surface separations one finds attractive and then strongly repulsive contributions. For surface charge density around one negative charge per 70 Å2 the full results for pressures resemble “secondary hydration force” measured in classical experiments in 1980s. When there is a tendency for ions to adsorb at the surfaces there is a marked change in behavior. The force is then oscillatory, reminiscent of results obtained with the surface force apparatus at low electrolyte concentration

Metal Ion-Induced Lateral Aggregation of Filamentous Viruses fd and M13

We report a detailed comparison between calculations of inter-filament interactions based on Monte-Carlo simulations and experimental features of lateral aggregation of bacteriophages fd and M13 induced by a number of divalent metal ions. The general findings are consistent with the polyelectrolyte nature of the virus filaments and confirm that the solution electrostatics account for most of the experimental features observed. One particularly interesting discovery is resolubilization for bundles of either fd or M13 viruses when the concentration of the bundle-inducing metal ion Mg2+ or Ca2+ is increased to large (\u3e100 mM) values. In the range of Mg2+ or Ca2+ concentrations where large bundles of the virus filaments are formed, the optimal attractive interaction energy between the virus filaments is estimated to be on the order of 0.01 kT per net charge on the virus surface when a recent analytical prediction to the experimentally defined conditions of resolubilization is applied. We also observed qualitatively distinct behavior between the alkali-earth metal ions and the divalent transition metal ions in their action on the charged viruses. The understanding of metal ions-induced reversible aggregation based on solution electrostatics may lead to potential applications in molecular biology and medicine

A multiscale model of protein adsorption on a nanoparticle surface

We present a methodology to quantify the essential interactions at the interface between inorganic solid nanoparticles (NPs) and biological molecules. Our model is based on pre-calculation of the repetitive contributions to the interaction from molecular segments, which allows us to efficiently scan a multitude of molecules and rank them by their adsorption affinity. The interaction between the biomolecular fragments and the nanomaterial are evaluated using a systematic coarse-graining scheme starting from all-atom molecular dynamics simulations. The NPs are modelled using a two-layer representation, where the outer layer is parameterized at the atomistic level and the core is treated at the continuum level using Lifshitz theory of dispersion forces. We demonstrate that the scheme reproduces the experimentally observed features of the NP protein coronas. To illustrate the use of the methodology, we compute the adsorption energies for human blood plasma proteins on gold NPs of different sizes as well as the preferred orientation of the molecules upon adsorption. The computed energies can be used for predicting the composition of the NP-protein corona for the corresponding material

Hydration interactions: aqueous solvent effects in electric double layers

A model for ionic solutions with an attractive short-range pair interaction between the ions is presented. The short-range interaction is accounted for by adding a quadratic non-local term to the Poisson-Boltzmann free energy. The model is used to study solvent effects in a planar electric double layer. The counter-ion density is found to increase near the charged surface, as compared with the Poisson-Boltzmann theory, and to decrease at larger distances. The ion density profile is studied analytically in the case where the ion distribution near the plate is dominated only by counter-ions. Further away from the plate the density distribution can be described using a Poisson-Boltzmann theory with an effective surface charge that is smaller than the actual one.Comment: 11 Figures in 13 files + LaTex file. 20 pages. Accepted to Phys. Rev. E. Corrected typos and reference

Effective interaction between helical bio-molecules

The effective interaction between two parallel strands of helical bio-molecules, such as deoxyribose nucleic acids (DNA), is calculated using computer simulations of the "primitive" model of electrolytes. In particular we study a simple model for B-DNA incorporating explicitly its charge pattern as a double-helix structure. The effective force and the effective torque exerted onto the molecules depend on the central distance and on the relative orientation. The contributions of nonlinear screening by monovalent counterions to these forces and torques are analyzed and calculated for different salt concentrations. As a result, we find that the sign of the force depends sensitively on the relative orientation. For intermolecular distances smaller than $6\AA$ it can be both attractive and repulsive. Furthermore we report a nonmonotonic behaviour of the effective force for increasing salt concentration. Both features cannot be described within linear screening theories. For large distances, on the other hand, the results agree with linear screening theories provided the charge of the bio-molecules is suitably renormalized.Comment: 18 pages, 18 figures included in text, 100 bibliog

Non-monotonic variation with salt concentration of the second virial coefficient in protein solutions

The osmotic virial coefficient $B_2$ of globular protein solutions is calculated as a function of added salt concentration at fixed pH by computer simulations of the ``primitive model''. The salt and counter-ions as well as a discrete charge pattern on the protein surface are explicitly incorporated. For parameters roughly corresponding to lysozyme, we find that $B_2$ first decreases with added salt concentration up to a threshold concentration, then increases to a maximum, and then decreases again upon further raising the ionic strength. Our studies demonstrate that the existence of a discrete charge pattern on the protein surface profoundly influences the effective interactions and that non-linear Poisson Boltzmann and Derjaguin-Landau-Verwey-Overbeek (DLVO) theory fail for large ionic strength. The observed non-monotonicity of $B_2$ is compared to experiments. Implications for protein crystallization are discussed.Comment: 43 pages, including 17 figure