Effective pair potentials in charged molecular systems: Water mediated interactions

On the evolution of the exclusion zone produced by hydrophilic surfaces: A contracted description

After exciting scientific debates about its nature, the development of the exclusion zone, a region near hydrophilic surfaces from which charged colloidal particles are strongly expelled, has been finally traced back to the diffusiophoresis produced by unbalanced ion gradients. This was done by numerically solving the coupled Poisson equation for electrostatics, the two stationary Stokes equations for low Reynolds numbers in incompressible fluids, and the Nernst-Planck equation for mass transport. Recently, it has also been claimed that the leading mechanism behind the diffusiophoretic phenomenon is electrophoresis [Esplandiu et al., Soft Matter 16, 3717 (2020)]. In this paper, we analyze the evolution of the exclusion zone based on a one-component interaction model at the Langevin equation level, which leads to simple analytical expressions instead of the complex numerical scheme of previous works, yet being consistent with it. We manage to reproduce the evolution of the exclusion zone width and the mean-square displacements of colloidal particles we measure near Nafion, a perfluorinated polymer membrane material, along with all characteristic time regimes, in a unified way. Our findings are also strongly supported by complementary experiments using two parallel planar conductors kept at a fixed voltage, mimicking the hydrophilic surfaces, and some computer simulations.

Excluded volume and ion-ion correlation effects on the ionic atmosphere around B-DNA: theory, simulations, and experiments

The ionic atmosphere around a nucleic acid regulates its stability in aqueous salt solutions. One major source of complexity in biological activities involving nucleic acids arises from the strong influence of the surrounding ions and water molecules on their structural and thermodynamic properties. Here, we implement a classical density functional theory for cylindrical polyelectrolytes embedded in aqueous electrolytes containing explicit (neutral hard sphere) water molecules at experimental solvent concentrations. Our approach allows us to include ion correlations as well as solvent and ion excluded volume effects for studying the structural and thermodynamic properties of highly charged cylindrical polyelectrolytes. Several models of size and charge asymmetric mixtures of aqueous electrolytes at physiological concentrations are studied. Our results are in good agreement with Monte Carlo simulations. Our numerical calculations display significant differences in the ion density profiles for the different aqueous electrolyte models studied. However, similar results regarding the excess number of ions adsorbed to the B-DNA molecule are predicted by our theoretical approach for different aqueous electrolyte models. These findings suggest that ion counting experimental data should not be used alone to validate the performance of aqueous DNA-electrolyte models.

Three-dimensional structures of RNA obtained by means of knowledge-based interaction potentials

We derive a set of effective potentials describing the interaction between pairs of nucleotides that belong to an RNA molecule. Such interaction potentials are then used as the main constituents of a simplified simulation model, which is tested in the description of small secondary structure motifs. Our simulated RNA hairpins are consistent with the experimental structures obtained by NMR.

Asymmetric charge renormalization for nanoparticles in aqueous media

The effective renormalized charge of nanoparticles in an aqueous electrolyte is essential to determine their solubility. By using a molecular model for the supporting aqueous electrolyte, we find that the effective renormalized charge of the nanoparticles is strongly dependent on the sign of the bare charge. Negatively charged nanoparticles have a lower effective renormalized charge than positively charged nanoparticles. The degree of asymmetry is a nonmonotonic function of the bare charge of the nanoparticle. We show that the effect is due to the asymmetric charge distribution of the water molecules, which we model using a simple three-site molecular structure of point charges.

Protein tyrosine phosphatase oligomerization studied by a combination of 15N NMR relaxation and 129Xe NMR. Effect of buffer containing arginine and glutamic acid

15N NMR relaxation and 129Xe NMR chemical shift measurements offer complementary information to study weak protein-protein interactions. They have been applied to study the oligomerization equilibrium of a low-molecular-weight protein tyrosine phosphatase in the presence of 50 mM arginine and 50 mM glutamic acid. These experimental conditions are shown to enhance specific protein-protein interactions while decreasing nonspecific aggregation. In addition, 129Xe NMR chemical shifts become selective reporters of one particular oligomer in the presence of arginine and glutamic acid, indicating that a specific Xe binding site is created in the oligomerization process. It is suggested that the multiple effects of arginine and glutamic acid are related to their effective excluded volume that favors specific protein association and the destabilization of partially unfolded forms that preferentially interact with xenon and are responsible for nonspecific protein aggregation.

A simple phenomenological fix for the dielectric constant within the reference interaction site model approach

The effective interactions among ions immersed in water are studied by means of the effective pair potentials (EPPs) [J. Chem. Phys. 2002, 117, 6133] obtained after contracting (integrating out) the degrees of freedom of the solvent molecules. The dressed interaction site theory (DIST) leads to a simple way of adjusting the effective dielectric constant of the model solvent to its experimental value at standard conditions. The molecular structure of the solvent is mirrored in the structure of the short-ranged component of the induced EPPs, with noticeable differences between the cases with trivial (ideal gas) and nontrivial (experimental) values of the dielectric constant. The shape of these EPPs remains almost invariant over the whole range of salt concentrations considered here. The asymptotic behavior of the EPP between two macroions obtained after contracting the supporting electrolyte (water molecules plus small ions) is also briefly discussed.

Two-Component Polypeptides Modeled with Effective Pair Potentials

We present Monte Carlo simulations performed within a model based on a set of distance-dependent effective potentials which are used to describe the interactions between a pair of distinct amino acids. These effective potentials are extracted from experimental correlation functions through the Ornstein-Zernike equations and adequate closure approximations. We focus our attention on the sequences of two specific residues, namely, alanine and glycine. The studied sequences are (a) (Ala)(12)-(Gly)(4)-(Ala)(12) and (b) three interacting chains of alternating alanines and glycines (with five residues per chain). The resulting structures are combinations of known secondary structures. More importantly, we verify that our simulated structures are in thermodynamic equilibrium by means of an estimation of the density of states.

Effective multipoles and Yukawa electrostatics in dressed molecule theory

In this paper we derive the multipolar expansion of the screened Coulomb potential in electrolyte solutions with molecular solvent. The solute and solvent molecules can have arbitrary sizes, shapes, and internal charge distributions. We use the exact statistical mechanical definition of renormalized charge distributions coming from "dressed molecule theory" to determine the effective multipoles of a molecule immersed in an electrolyte. The effects of many-body correlations are fully included in our formally exact theory. We restrict ourselves to sufficiently dilute solutions so the screened Coulomb potential decays for large distances like a Yukawa function, exp(-kappa r)/r, where r is the distance and 1/kappa is the decay length (it is normally different from the Debye length). The resulting "Yukawa electrostatics" differ in many respects from ordinary, unscreened electrostatics. The "Yukawa charge" of a molecule (the lowest order moment in the multipolar expansion) is in general not equal to its Coulombic charge and it is not the integral of the renormalized charge distribution of the molecule. Moreover, as shown in this paper, the multipolar expansion of the Yukawa potential does not correspond, contrary to the case of the Coulomb potential, to its asymptotic expansion for large r. As a consequence, the charge term in the multipolar expansion is not the leading term in the asymptotic expansion. Instead, for large r values, multipoles of all orders contribute to the leading asymptotic term. Thus, the electrostatic potential from, for example, an electroneutral solvent molecule in an electrolyte solution has generally the same range as that from an ion. The proper asymptotic expansion for electrostatic interactions in electrolytes is derived. It is briefly shown how the multipole expansion formalism can also be applied in the Poisson-Boltzmann approximation for primitive model electrolytes.

A reference interaction site model approach to depletion forces induced by hard rodlike particles

We study the effective interactions among large hard spherical colloidal particles induced by small hard rodlike particles and compare them with those induced by small hard spherical particles to highlight the specific effects due to the anisotropic shape of the former. This is done by determining the effective pair potentials within the framework of the reference interaction site model approach. The rodlike particles are modeled as N nonoverlapping spherical units arranged in a straight line, so that their total length is N times their transversal diameter. These results are compared against those obtained in the Asakura-Oosawa limit.

Polypeptide foldings obtained with effective pair potentials

We present a model of protein folding which is based on a potential function that describes the effective interaction between two amino acids (alanines, in this case). Our model is consistent with the formation of two important secondary structures, namely, an alpha-helix and a beta-ladder. In each case, we estimate the density of states using a random walk in energy space. This function allows the direct calculation of certain thermodynamic properties. By means of the configurational temperature, we also verify that the obtained polypeptides are in their native state.

Molecular multivalent electrolytes: microstructure and screening lengths

We study small rod-like molecular electrolytes solutions with their corresponding atomic counterions. The asymptotic length scales (decay length and wavelength) of the structural correlations are analyzed using the formalism of the dressed interaction site theory (DIST). The correlation functions are determined using the reference interaction site model equation complemented with a mixed approach in which the hypernetted-chain closure is used for the repulsive interactions, and the mean spherical approximation is used for the attractive interactions. The results from this scheme are in good agreement with the Monte Carlo computer simulations reported here. The asymptotic properties of the correlation functions of this molecular system are compared against those corresponding to two related simple (atomic) electrolyte models. The main conclusion is that the molecular structure of the ions lowers by two orders of magnitude the concentration at which the transition from monotonic to oscillatory decay occurs.

Effective pair potentials between protein amino acids

We present an effective potential describing the interaction between pairs of residues (alanines, in this case) that belong to a protein. The effective potential is extracted from an experimental correlation function, by means of the Ornstein-Zernike equation together with a closure approximation. It is found that the most relevant features of the effective potential are consistent with the formation of two different secondary structures of proteins.

Polyelectrolytes in the presence of multivalent ions: gelation versus segregation

We analyze solutions of strongly charged chains bridged by linkers such as multivalent ions. The gelation induced by the strong short range electrostatic attractions is dramatically suppressed by the long range electrostatic correlations due to the charge along the non-cross-linked monomers and ions. A modified Debye-Hückel approach of cross-linked clusters of charged chains is used to determine the mean field gelation transition self-consistently. Highly dilute polyelectrolyte solutions tend to segregate macroscopically. Semidilute solutions can form gels if the Bjerrum length l(B) and the distance between neighboring charged monomers along the chain b are both greater than the ion size a.