-
A closed-form, analytical approximation for apparent surface charge and electric field of molecules
Authors:
Dan Folescu,
Alexey V. Onufriev
Abstract:
Closed-form, analytical approximations for electrostatic properties of molecules are of unique value, as these can provide computational speed, versatility, and physical insight. Here, we derive a simple, closed-form formula for the apparent surface charge (ASC), as well as for the electric field, generated by a molecular charge distribution in aqueous solution. The approximation, with no fitted p…
▽ More
Closed-form, analytical approximations for electrostatic properties of molecules are of unique value, as these can provide computational speed, versatility, and physical insight. Here, we derive a simple, closed-form formula for the apparent surface charge (ASC), as well as for the electric field, generated by a molecular charge distribution in aqueous solution. The approximation, with no fitted parameters, is tested against numerical solutions of the Poisson equation, where it yields a significant speed-up. For neutral small molecules, the hydration free energies estimated from the closed-form ASC formula are within 0.8 kcal/mol RMSD from the numerical Poisson reference; the electric field at the surface is in quantitative agreement with the reference. Performance of the approximation is also tested on larger structures, including a protein, a DNA fragment, and a viral receptor-target complex. For all structures tested, a near quantitative agreement with the numerical Poisson reference is achieved, except in regions of high negative curvature, where the new approximation is still qualitatively correct. A unique efficiency feature of the proposed "source-based" closed-form approximation is that the ASC and electric field can be estimated individually at any point or surface patch, without the need to obtain the full global solution. An open source software implementation of the method is available: http://people.cs.vt.edu/~onufriev/CODES/aasc.zip
△ Less
Submitted 10 March, 2022; v1 submitted 22 June, 2021;
originally announced June 2021.
-
Strongly bent double-stranded DNA: reconciling theory and experiment
Authors:
Aleksander V. Drozdetski,
Abhishek Mukhopadhyay,
Alexey V. Onufriev
Abstract:
The strong bending of polymers is poorly understood. We propose a general quantitative framework of polymer bending that includes both the weak and strong bending regimes on the same footing, based on a single general physical principle. As the bending deformation increases beyond a certain (polymer-specific) point, the change in the convexity properties of the effective bending energy of the poly…
▽ More
The strong bending of polymers is poorly understood. We propose a general quantitative framework of polymer bending that includes both the weak and strong bending regimes on the same footing, based on a single general physical principle. As the bending deformation increases beyond a certain (polymer-specific) point, the change in the convexity properties of the effective bending energy of the polymer makes the harmonic deformation energetically unfavorable: in this strong bending regime the energy of the polymer varies linearly with the average bending angle as the system follows the convex hull of the deformation energy function. For double-stranded DNA, the effective bending deformation energy becomes non-convex for bends greater than about 2 degrees per base-pair, equivalent to the curvature of a closed circular loop of about 160 base pairs. A simple equation is derived for the polymer loop energy that covers both the weak and strong bending regimes. The theory shows quantitative agreement with recent DNA cyclization experiments on short DNA fragments, while maintaining the expected agreement with experiment in the weak bending regime. Counter-intuitively, cyclization probability (j-factor) of very short DNA loops is predicted to increase with decreasing loop length; the j-factor reaches its minimum for loops of 45 base pairs. Atomistic simulations reveal that the attractive component of the short-range Lennard-Jones interaction between the backbone atoms can explain the underlying non-convexity of the DNA effective bending energy, leading to the linear bending regime. Applicability of the theory to protein-DNA complexes, including the nucleosome, is discussed.
△ Less
Submitted 2 July, 2019;
originally announced July 2019.
-
Multi-shell model of ion-induced nucleic acid condensation
Authors:
Igor S. Tolokh,
Aleksander Drozdetski,
Lois Pollack,
Nathan A. Baker,
Alexey V. Onufriev
Abstract:
We present a semi-quantitative model of condensation of short nucleic acid (NA) duplexes induced by tri-valent cobalt(III) hexammine (CoHex) ions. The model is based on partitioning of bound counterion distribution around singleNA duplex into "external" and "internal" ion binding shells distinguished by the proximity to duplex helical axis. In the aggregated phase the shells overlap, which leads t…
▽ More
We present a semi-quantitative model of condensation of short nucleic acid (NA) duplexes induced by tri-valent cobalt(III) hexammine (CoHex) ions. The model is based on partitioning of bound counterion distribution around singleNA duplex into "external" and "internal" ion binding shells distinguished by the proximity to duplex helical axis. In the aggregated phase the shells overlap, which leads to significantly increased attraction of CoHex ions in these overlaps with the neighboring duplexes. The duplex aggregation free energy is decomposed into attractive and repulsive components in such a way that they can be represented by simple analytical expressions with parameters derived from molecular dynamic (MD) simulations and numerical solutions of Poisson equation. The short-range interactions described by the attractive term depend on the fractions of bound ions in the overlapping shells and affinity of CoHex to the "external" shell of nearly neutralized duplex. The repulsive components of the free energy are duplex configurational entropy loss upon the aggregation and the electrostatic repulsion of the duplexes that remains after neutralization by bound CoHex ions. The estimates of the aggregation free energy are consistent with the experimental range of NA duplex condensation propensities, including the unusually poor condensation of RNA structures and subtle sequence effects upon DNA condensation. The model predicts that, in contrast to DNA, RNA duplexes may condense into tighter packed aggregates with a higher degree of duplex neutralization. The model also predicts that longer NA fragments will condense more readily than shorter ones. The ability of this model to explain experimentally observed trends in NA condensation, lends support to proposed NA condensation picture based on the multivalent "ion binding shells".
△ Less
Submitted 6 April, 2016;
originally announced April 2016.
-
Understanding Nucleic Acid Structural Changes by Comparing Wide-Angle X-ray Scattering (WAXS) Experiments to Molecular Dynamics Simulations
Authors:
Suzette A. Pabit,
Andrea M. Katz,
Igor S. Tolokh,
Aleksander Drozdetski,
Nathan Baker,
Alexey V. Onufriev,
Lois Pollack
Abstract:
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 W…
▽ More
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 in 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.
△ Less
Submitted 15 May, 2016; v1 submitted 26 December, 2015;
originally announced December 2015.
-
The role of correlation and solvation in ion interactions with B-DNA
Authors:
Maria L. Sushko,
Dennis G. Thomas,
Suzette A. Pabit,
Lois Pollack,
Alexey V. Onufriev,
Nathan A. Baker
Abstract:
The ionic atmospheres around nucleic acids play important roles in biological function. Large-scale explicit solvent simulations coupled to experimental assays such as anomalous small-angle X-ray scattering (ASAXS) can provide important insights into the structure and energetics of such atmospheres but are time- and resource-intensive. In this paper, we use classical density functional theory (cDF…
▽ More
The ionic atmospheres around nucleic acids play important roles in biological function. Large-scale explicit solvent simulations coupled to experimental assays such as anomalous small-angle X-ray scattering (ASAXS) can provide important insights into the structure and energetics of such atmospheres but are time- and resource-intensive. In this paper, we use classical density functional theory (cDFT) to explore the balance between ion-DNA, ion-water, and ion-ion interactions in ionic atmospheres of RbCl, SrCl$_2$, and CoHexCl$_3$ (cobalt hexammine chloride) around a B-form DNA molecule. The accuracy of the cDFT calculations was assessed by comparison between simulated and experimental ASAXS curves, demonstrating that an accurate model should take into account ion-ion correlation and ion hydration forces, DNA topology, and the discrete distribution of charges on DNA strands. As expected, these calculations revealed significant differences between monovalent, divalent, and trivalent cation distributions around DNA. About half of the DNA-bound Rb$^+$ ions penetrate into the minor groove of the DNA and half adsorb on the DNA strands. The fraction of cations in the minor groove decreases for the larger Sr$^{2+}$ ions and becomes zero for CoHex$^{3+}$ ions, which all adsorb on the DNA strands. The distribution of CoHex$^{3+}$ ions is mainly determined by Coulomb and steric interactions, while ion-correlation forces play a central role in the monovalent Rb$^+$ distribution and a combination of ion-correlation and hydration forces affect the Sr$^{2+}$ distribution around DNA.
△ Less
Submitted 6 December, 2015; v1 submitted 25 June, 2015;
originally announced June 2015.
-
Accurate Evaluation of Charge Asymmetry in Aqueous Solvation
Authors:
Abhishek Mukhopadhyay,
Igor S. Tolokh,
Alexey V. Onufriev
Abstract:
Charge hydration asymmetry (CHA)--a characteristic dependence of hydration free energy on the sign of the solute charge--quantifies the asymmetric response of water to electric field at microscopic level. Accurate estimates of CHA are critical for understanding hydration effects ubiquitous in chemistry and biology. However, measuring hydration energies of charged species is fraught with significan…
▽ More
Charge hydration asymmetry (CHA)--a characteristic dependence of hydration free energy on the sign of the solute charge--quantifies the asymmetric response of water to electric field at microscopic level. Accurate estimates of CHA are critical for understanding hydration effects ubiquitous in chemistry and biology. However, measuring hydration energies of charged species is fraught with significant difficulties, which lead to unacceptably large (up to 300%) variation in the available estimates of the CHA effect. We circumvent these difficulties by developing a framework which allows us to extract and accurately estimate the intrinsic propensity of water to exhibit CHA from accurate experimental hydration free energies of neutral polar molecules. Specifically, from a set of 504 small molecules we identify two pairs that are analogous, with respect to CHA, to the K+/F- pair--a classical probe for the effect. We use these "CHA-conjugate" molecule pairs to quantify the intrinsic charge-asymmetric response of water to the microscopic charge perturbations: the asymmetry of the response is strong, ~50% of the average hydration free energy of these molecules. The ability of widely used classical water models to predict hydration energies of small molecules correlates with their ability to predict CHA.
△ Less
Submitted 11 May, 2015; v1 submitted 1 April, 2015;
originally announced April 2015.