Nanoscale structures and mechanics of peptide nucleic acids

DNA compaction and chromatin dynamics: The role of cationic polyamines and proteins

DNA compaction by polyaminic cations and proteins involves reversible condensation mechanisms. Polyamines, metal cations, and histone proteins are utilized to compact lengthy DNA chains. Chromatin organization begins with nucleosomal arrays, further compacted by linker histones. Various factors such as DNA methylation, histone modifications, and non-histone proteins influence chromatin structure. Posttranslational modifications like acetylation and methylation alter nucleosome shape. Polyamines induce significant phase transitions, while cationic surfactants drive conformational changes in DNA. In sperm cells, protamines replace histones, leading to dense DNA packing. Despite advances, unresolved aspects persist in understanding the dynamic regulation of chromatin structure, highlighting avenues for future research. An overview of current knowledge and cutting-edge discoveries in the field of reversible DNA compaction induced by charged polyamines and histone proteins is presented in this work, highlighting emerging mechanisms of chromatin compaction and their relevance to cellular function, disease, and potential therapeutic strategies.

Stretching and twisting of double-stranded RNA under forces: Unwinding mechanism and base-pair dependent elasticity

We used all-atom molecular dynamics simulations to investigate the mechanical response of double-stranded RNA (dsRNA) by applying various forces. We used the helical rise and helical twist, as well as a newly defined helical diameter, to characterize the stretching and twisting of dsRNA. The results indicate that dsRNA unwinds when stretched, accompanied by a linear increase in helical rise and helical diameter. Then, we utilized the normal modes, which are linear combinations of helical modes, to elucidate the underlying mechanism of dsRNA unwinding from an energetic perspective. On the other hand, we employed a stiffness matrix based on a rigid base pair model to examine the base-pair dependence of twist elasticity for dsRNA, as well as stretch elasticities with respect to the helical rise and helical diameter. The results show that the force induces variations in the local elasticities and their couplings of dsRNA, which are closely related to the distributions of base pairs. The mean stretch and twist elasticities can be considered as constants within the measurement uncertainties; however, their couplings demonstrate a slight linear dependency on applied force.

A review on salt-induced DNA compaction and charge inversion

This review delves into the reversible process of DNA compaction, vital for cellular functions like replication and transcription. The study highlights how various cations assist in the condensation of DNA chains, highlighting their specificity. The impact of the ionic environment on chromatin characteristics is discussed, emphasizing the roles of mono- and divalent cations in neutralizing DNA charge and promoting compaction. Trivalent ions induce significant compaction, while divalent ions also contribute, albeit less strongly. Charge inversion, facilitated by high concentrations of multivalent counterions, affects DNA condensation dynamics. Manipulating solution pH and dielectric constant can alter charge inversion bidirectionally. The hydrophobic effect driven by organic cations plays a crucial role in DNA compaction. The review underscores the implications of charge inversion, including macroscopic phase separation and DNA precipitation, driven by the binding of cationic micelles to DNA.

How Rigid Are Anthranilamide Molecular Electrets?

As important as molecular electrets are for electronic materials and devices, conformational fluctuations strongly impact their macrodipoles and intrinsic properties. Herein, we employ molecular dynamics (MD) simulations with the polarizable charge equilibrium (PQEq) method to investigate the persistence length (LP) of molecular electrets composed of anthranilamide (Aa) residues. The PQEq-MD dissipates the accepted static notions about Aa macromolecules, and LP represents the shortest Aa rigid segments. The classical model with a single LP value does not describe these oligomers. Introducing multiple LP values for the same macromolecule follows the observed trends and discerns the enhanced rigidity in their middle sections from the reduced stiffness at their terminal regions. Furthermore, LP distinctly depends on solvent polarity. The Aa oligomers maintain extended conformations in nonpolar solvents with LP exceeding 4 nm, while in polar media, increased conformational fluctuations reduce LP to about 2 nm. These characteristics set key guidelines about the utility of Aa conjugates for charge-transfer systems within organic electronics and energy engineering.

Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.

Dynamics of terminal fraying-peeling and hydrogen bonds dictates the sequential vs. cooperative melting pathways of nanoscale DNA and PNA triplexes

Peptide nucleic acids (PNAs) are charge-neutral synthetic DNA/RNA analogues. In many aspects of biology and biotechnology, the details of DNA and PNA melting reaction coordinates are crucial, and their associative/dissociative details remain inadequately understood. In the current study, we have attempted to gain insights into comparative melting pathways and binding affinity of iso-sequences of an 18-mer PNA-DNA-PNA triplex and the analogous DNA-DNA-DNA triplex, and DNA-DNA and PNA-DNA duplexes. It is intriguing that while the DNA-DNA-DNA triplex melts in two sequential steps, the PNA-DNA-PNA triplex melts in a single step and the mechanistic aspects for this difference are still not clear. We report an all-atom molecular dynamics simulation of both complexes in the temperature range of 300 to 500 K with 20 K intervals. Based on the trajectory analysis, we provide evidence that the association and dissociation are dictated by the differences in fraying-peeling effects from either terminus to the center in a zipper pattern among the PNA-DNA-PNA triplex and DNA-DNA-DNA triplexes. These are shown to be governed by the different characteristics of H-bonding, RMSD, and Free Energy Landscape (FEL) as analyzed by PCA, leading to the DNA-DNA-DNA triplex exhibiting sequential melting, while the PNA-DNA-PNA triplex shows cooperative melting of the whole fragment in a single-step. The PNA-DNA-PNA triplex base pairs are thermodynamically more stable than the DNA-DNA-DNA triplex, with the binding affinity of PNA-TFO to the PNA : DNA duplex being higher than that of DNA-TFO to the DNA : DNA duplex. The investigation of the association/dissociation of PNA-TFO to the PNA-DNA duplex has relevance and importance in the emerging effective applications of oligonucleotide therapy.

Effect of temperature on anisotropic bending elasticity of dsRNA: an all-atom molecular dynamics simulation

Employing all-atom molecular dynamics simulations, we examined the temperature-dependent behavior of bending elasticity in double-stranded RNA (dsRNA). Specifically, we focused on the bending persistence length and its constituent components, namely, the tilt and roll stiffness. Our results revealed a near-linear decrease in these stiffness components as a function of temperature, thereby highlighting the increased flexibility of dsRNA at elevated temperatures. Furthermore, our data revealed a significant anisotropy in dsRNA bending elasticity, which diminished with increasing temperature, attributable to marked disparities in tilt and roll stiffness components. We delineated the underlying biophysical mechanisms and corroborated our findings with extant literature. These observations offer salient implications for advancing our understanding of nucleic acid elasticity, and are pertinent to potential medical applications.

Systematic assessment of the flexibility of uracil damaged DNA

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3' base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5' T adjacent to the uracil strongly promotes flexibility, but other 5' bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping.Communicated by Ramaswamy H. Sarma.

DNA groove preference shift upon phosphorylation of a protamine-like cationic peptide

Protamines, arginine-rich DNA-binding proteins, are responsible for chromatin compaction in sperm cells, but their DNA groove preference, major or minor, is not clearly identified. We herein study the DNA groove preference of a short protamine-like cationic peptide before and after phosphorylation, using all-atom molecular dynamics and umbrella sampling simulations. According to various thermodynamic and structural analyses, a peptide in its non-phosphorylated native state prefers the minor groove over the major groove, but phosphorylation of the peptide bound to the minor groove not only reduces its binding affinity but also brings a serious deformation of the minor groove, eliminating the minor-groove preference. As protamines are heavily phosphorylated before binding to DNA, we expect that the structurally disordered phosphorylated protamines would prefer major grooves to enter into DNA during spermatogenesis.

Effect of phosphorylation of protamine-like cationic peptide on the binding affinity to DNA

Protamines are more arginine-rich and more basic than histones and are responsible for providing a highly compacted shape to the sperm heads in the testis. Phosphorylation and dephosphorylation are two events that occur in the late phase of spermatogenesis before the maturation of sperms. In this work, we have studied the effect of phosphorylation of protamine-like cationic peptides using all-atom molecular dynamics simulations. Through thermodynamic analyses, we found that phosphorylation reduces the binding efficiency of such cationic peptides on DNA duplexes. Peptide phosphorylation leads to a less efficient DNA condensation, due to a competition between DNA-peptide and peptide-peptide interactions. We hypothesize that the decrease of peptide bonds between DNA together with peptide self-assembly might allow an optimal re-organization of chromatin and an efficient condensation through subsequent peptide dephosphorylation. Based on the globular and compact conformations of phosphorylated peptides mediated by arginine-phosphoserine H-bonding, we furthermore postulate that phosphorylated protamines could more easily intrude into chromatin and participate to histone release through disruption of histone-histone and histone-DNA binding during spermatogenesis.

Diameter Dependent Melting and Softening of dsDNA Under Cylindrical Confinement

Carbon nanotubes (CNTs) are considered promising candidates for biomolecular confinement, including DNA encapsulation for gene delivery. Threshold values of diameters have been reported for double-stranded DNA (dsDNA) encapsulation inside CNTs. We have performed all-atom molecular dynamics (MD) simulations of dsDNAs confined inside single-walled CNTs (SWCNTs) at the physiologically relevant temperature of 300 K. We found that the dsDNA can be confined without being denatured only when the diameter of the SWCNT exceeds a threshold value. Below this threshold diameter, the dsDNA gets denatured and melts even at the temperature of 300 K. Our simulations using SWCNTs with chirality indices (20,20) to (30,30) at 300 K found the critical diameter to be 3.25 nm (corresponding to (24,24) chirality). Analyses of the hydrogen bonds (H-bonds), Van der Walls (VdW) energy, and other inter-base interactions show drastic reduction in the number of H-bonds, VdW energy, and electrostatic energies between the bases of dsDNA when it is confined in narrower SWCNTs (up to diameter of 3.12 nm). On the other hand, the higher interaction energy between the dsDNA and the SWCNT surface in narrower SWCNTs assists in the melting of the dsDNA. Electrostatic mapping and hydration status analyses show that the dsDNA is not adequately hydrated and the counter ion distribution is not uniform below the critical diameter of the SWCNT. As properly hydrated counter ions provide stability to the dsDNA, we infer that the inappropriate hydration of counter ions and their non-uniform distribution around the dsDNA cause the melting of the dsDNA inside SWCNTs of diameter below the critical value of 3.25 nm. For confined dsDNAs that do not get denatured, we computed their elastic properties. The persistence length of dsDNA was found to increase by a factor of about two and the torsional stiffness by a factor of 1.5 for confinement inside SWCNTs of diameters up to 3.79 nm, the stretch modulus also following nearly the same trend. Interestingly, for higher diameters of SWCNT, 3.79 nm and above, the dsDNA becomes more flexible, demonstrating that the mechanical properties of the dsDNA under cylindrical confinement depend non-monotonically on the confinement diameter.