Igor S. Tolokh

Charge Hydration Asymmetry: The Basic Principle and How to Use It to Test and Improve Water Models

Charge hydration asymmetry (CHA) manifests itself in the experimentally observed strong dependence of free energy of ion hydration on the sign of the ion charge. This asymmetry is not consistently accounted for by popular models of solvation; its magnitude varies greatly between the models. While it is clear that CHA is somehow related to charge distribution within a water molecule, the exact nature of this relationship is unknown. We propose a simple, yet general and rigorous criterion that relates rotational and charge inversion properties of a water molecules charge distribution with its ability to cause CHA. We show which electric multipole components of a water molecule are key to explain its ability for asymmetric charge hydration. We then test several popular water models and explain why specific models show none, little, or strong CHA in simulations. We use the gained insight to derive an analogue of the Born equation that includes the missing physics necessary to account for CHA and does not rely on redefining the continuum dielectric boundary. The proposed formula is as simple as the original, does not contain any fitting parameters, and predicts hydration free energies and entropies of spherical cations and anions within experimental uncertainty. Our findings suggest that the gap between the practical continuum electrostatics framework and the more fundamental explicit solvent treatment may be reduced considerably by explicitly introducing CHA into the existing continuum framework.

Why double-stranded RNA resists condensation

The addition of small amounts of multivalent cations to solutions containing double-stranded DNA leads to inter-DNA attraction and eventual condensation. Surprisingly, the condensation is suppressed in double-stranded RNA, which carries the same negative charge as DNA, but assumes a different double helical form. Here, we combine experiment and atomistic simulations to propose a mechanism that explains the variations in condensation of short (25 base-pairs) nucleic acid (NA) duplexes, from B-like form of homopolymeric DNA, to mixed sequence DNA, to DNA:RNA hybrid, to A-like RNA. Circular dichroism measurements suggest that duplex helical geometry is not the fundamental property that ultimately determines the observed differences in condensation. Instead, these differences are governed by the spatial variation of cobalt hexammine (CoHex) binding to NA. There are two major NA-CoHex binding modes—internal and external—distinguished by the proximity of bound CoHex to the helical axis. We find a significant difference, up to 5-fold, in the fraction of ions bound to the external surfaces of the different NA constructs studied. NA condensation propensity is determined by the fraction of CoHex ions in the external binding mode.

Conduction through the Inward Rectifier Potassium Channel, Kir2.1, Is Increased by Negatively Charged Extracellular Residues

Ion channel conductance can be influenced by electrostatic effects originating from fixed “surface” charges that are remote from the selectivity filter. To explore whether surface charges contribute to the conductance properties of Kir2.1 channels, unitary conductance was measured in cell-attached recordings of Chinese hamster ovary (CHO) cells transfected with Kir2.1 channels over a range of K+ activities (4.6–293.5 mM) using single-channel measurements as well as nonstationary fluctuation analysis for low K+ activities. K+ ion concentrations were shown to equilibrate across the cell membrane in our studies using the voltage-sensitive dye DiBAC4(5). The dependence of γ on the K+ activity (aK) was fit well by a modified Langmuir binding isotherm, with a nonzero intercept as aK approaches 0 mM, suggesting electrostatic surface charge effects. Following the addition of 100 mM N-methyl-d-glucamine (NMG+), a nonpermeant, nonblocking cation or following pretreatment with 50 mM trimethyloxonium (TMO), a carboxylic acid esterifying agent, the γ–aK relationship did not show nonzero intercepts, suggesting the presence of surface charges formed by glutamate or aspartate residues. Consistent with surface charges in Kir2.1 channels, the rates of current decay induced by Ba2+ block were slowed with the addition of NMG or TMO. Using a molecular model of Kir2.1 channels, three candidate negatively charged residues were identified near the extracellular mouth of the pore and mutated to cysteine (E125C, D152C, and E153C). E153C channels, but not E125C or D152C channels, showed hyperbolic γ–aK relationships going through the origin. Moreover, the addition of MTSES to restore the negative charges in E53C channels reestablished wild-type conductance properties. Our results demonstrate that E153 contributes to the conductance properties of Kir2.1 channels by acting as a surface charge.

Introducing Charge Hydration Asymmetry into the Generalized Born Model.

The effect of charge hydration asymmetry (CHA)—non-invariance of solvation free energy upon solute charge inversion—is missing from the standard linear response continuum electrostatics. The proposed charge hydration asymmetric–generalized Born (CHA–GB) approximation introduces this effect into the popular generalized Born (GB) model. The CHA is added to the GB equation via an analytical correction that quantifies the specific propensity of CHA of a given water model; the latter is determined by the charge distribution within the water model. Significant variations in CHA seen in explicit water (TIP3P, TIP4P-Ew, and TIP5P-E) free energy calculations on charge-inverted “molecular bracelets” are closely reproduced by CHA–GB, with the accuracy similar to models such as SEA and 3D-RISM that go beyond the linear response. Compared against reference explicit (TIP3P) electrostatic solvation free energies, CHA–GB shows about a 40% improvement in accuracy over the canonical GB, tested on a diverse set of 248 rigid small neutral molecules (root mean square error, rmse = 0.88 kcal/mol for CHA–GB vs 1.24 kcal/mol for GB) and 48 conformations of amino acid analogs (rmse = 0.81 kcal/mol vs 1.26 kcal/mol). CHA–GB employs a novel definition of the dielectric boundary that does not subsume the CHA effects into the intrinsic atomic radii. The strategy leads to finding a new set of intrinsic atomic radii optimized for CHA–GB; these radii show physically meaningful variation with the atom type, in contrast to the radii set optimized for GB. Compared to several popular radii sets used with the original GB model, the new radii set shows better transferability between different classes of molecules.

Opposing Effects of Multivalent Ions on the Flexibility of DNA and RNA.

Increasing the concentration of counterions (salt) is known to reduce the bending persistence length of DNA. Here we use atomistic molecular dynamics simulations to predict that multivalent counterions have the opposite effect on double-stranded RNA, increasing its bending rigidity by at least 30%. This counterintuitive effect is observed for various tri- and tetravalent ions alike, and is robust to methodological details and the RNA sequence. In contrast to DNA, multivalent counterions bind inside the RNA major groove, causing significant contraction of the molecule along its helical axis-as a result, its further deformation due to bending becomes energetically more expensive compared to bending without bound multivalent ions. Thus, the relationship between mechanical properties of a charged polymer and its ionic atmosphere may be richer than previously thought.

Hydration of uracil and thymine methylderivatives: a Monte Carlo simulation.

The simulation performed shows that under methylation of uracil and thymine NH-groups the interaction energy between a base and water (Uwb) is increased. It is also detected that the increase in this energy was observed in the 1st and the 3rd sectors. These conclusions do not confirm the assumption made in the literature on the character of an interaction between methylated bases and water. According to this assumption, when the NH-groups are methylated, the energy of Uwb in these sectors decreases as a result of the van der Waals interactions between a methyl group and water, whose energy compensates the increase in the Uwb energy due to the breaking of an H-bond. Regularity of water molecules near a hydrophobic group under the hydration of polar molecules is detected for the first time.

Nature of the stacking of nucleic acid bases in water : a Monte Carlo study

Abstract The results of a Monte Carlo simulation of the hydration of uracil and thymine molecules, their stacked dimers and hydrogen-bonded base pairs are presented. Simulations have been performed in a cluster approximation. The semiempirical atom-atom potential functions have been used (cluster consisting of 200 water molecules). It has been shown that stacking interactions of uracil and thymine molecules in water arise mainly due to the increase in the water—water interaction during the transition from monomers to dimer. It has been found that stacked base associates are more favourable than base pairs in water. This is mainly due to the energetically more preferable structure of water around stacks.

Understanding nucleic acid structural changes by comparing wide-angle x-ray scattering (WAXS) experiments to molecular dynamics simulations

Wide-angle x-ray scattering (WAXS) is emerging as a powerful tool for increasing the resolution of solution structure measurements of biomolecules. Compared to its better known complement, small angle x-ray scattering (SAXS), WAXS targets higher scattering angles and can enhance structural studies of molecules by accessing finer details of solution structures. Although the extension from SAXS to WAXS is easy to implement experimentally, the computational tools required to fully harness the power of WAXS are still under development. Currently, WAXS is employed to study structural changes and ligand binding in proteins; however, the methods are not as fully developed for nucleic acids. Here, we show how WAXS can qualitatively characterize nucleic acid structures as well as the small but significant structural changes driven by the addition of multivalent ions. We show the potential of WAXS to test all-atom molecular dynamics (MD) simulations and to provide insight into understanding how the trivalent ion cobalt(III) hexammine (CoHex) affects the structure of RNA and DNA helices. We find that MD simulations capture the RNA structural change that occurs due to addition of CoHex.