Electric field of DNA in solution: who is in charge?

In solution, DNA is a highly charged macromolecule which bears a unit of negative charge on each phosphate of its sugar-phosphate backbone. Although partially compensated by counterions adsorbed at or condensed near it, DNA still produces a substantial electric field in its vicinity, which is screened by buffer electrolyte at longer distances from the DNA. Such field has been explored so far predominantly within the scope of a primitive model of the electrolytic solution, not considering more complicated structural effects of the water solvent. In this paper we investigate the distribution of electric field around DNA using linear response nonlocal electrostatic theory, applied here for helix-specific charge distributions, and compare the predictions of such theory with specially performed fully atomistic large scale molecular dynamics simulations. The main finding of this study is that oscillations in the electrostatic potential distribution are present around DNA, caused by the overscreening effect of structured water. Surprisingly, electrolyte ions at physiological concentrations do not strongly disrupt these oscillations, and rather distribute according to these oscillating patterns, indicating that water structural effects dominate the short-range electrostatics. We also show that (i) structured water adsorbed in the grooves of DNA lead to a positive electrostatic potential core, (ii) the Debye length some 10 {\AA} away from the DNA is reduced, effectively renormalised by the helical pitch of the DNA, and (iii) Lorentzian contributions to the nonlocal dielectric function of water, effectively reducing the dielectric constant close to the DNA, enhances the overall electric field. The impressive agreement between the atomistic simulations and the developed theory substantiates the use of nonlocal electrostatics when considering solvent effects in molecular processes in biology.

The structure and physical properties of a packaged bacteriophage particle.

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

Length-dependent Intramolecular Coil-to-Globule Transition in Poly(ADP-ribose) Induced by Cations.

Poly(ADP-ribose) (PAR), as part of a post-translational modification, serves as a flexible scaffold for noncovalent protein binding. Such binding is influenced by PAR chain length through a mechanism yet to be elucidated. Structural insights have been elusive, partly due to the difficulties associated with synthesizing PAR chains of defined lengths. Here, we employ an integrated approach combining molecular dynamics (MD) simulations with small-angle X-ray scattering (SAXS) experiments, enabling us to identify highly heterogeneous ensembles of PAR conformers at two different, physiologically relevant lengths: PAR 15 and PAR 22 . Our findings reveal that numerous factors including backbone conformation, base stacking, and chain length contribute to determining the structural ensembles. We also observe length-dependent compaction of PAR upon the addition of small amounts of Mg 2+ ions, with the 22-mer exhibiting ADP-ribose bundles formed through local intramolecular coil-to-globule transitions. This study illuminates how such bundling could be instrumental in deciphering the length-dependent action of PAR.