Communication: Asymmetrical cation movements through G-quadruplex DNA

Stability of Two-Quartet G-Quadruplexes and Their Dimers in Atomistic Simulations

G-quadruplexes (GQs) are four-stranded noncanonical DNA and RNA architectures that can be formed by guanine-rich sequences. The stability of GQs increases with the number of G-quartets, and three G-quartets generally form stable GQs. However, the stability of two-quartet GQs is an open issue. To understand the intrinsic stability of two-quartet GQ stems, we have carried out a series of unbiased molecular dynamics (MD) simulations (505 μs in total) of two- and four-quartet DNA and RNA GQs, with attention paid mainly to parallel-stranded arrangements. We used AMBER DNA parmOL15 and RNA parmOL3 force fields and tested different ion and water models. Two-quartet parallel-stranded DNA GQs unfolded in all the simulations, while the equivalent RNA GQ was stable in most of the simulations. GQs composed of two stacked units of two-quartet GQs were stable for both DNA and RNA. The simulations suggest that a minimum of three quartets are needed to form an intrinsically stable all-anti parallel-stranded DNA GQ. Parallel two-quartet DNA GQ may exist if substantially stabilized by another molecule or structural element, including multimerization. On the other hand, we predict that isolated RNA two-quartet parallel GQs may form, albeit being weakly stable. We also show that ionic parameters and water models should be chosen with caution because some parameter combinations can cause spurious instability of GQ stems. Some in-so-far unnoticed limitations of force-field description of multiple ions inside the GQs are discussed, which compromise the capability of simulations to fully capture the effect of increase in the number of quartets on the GQ stability.

Cationic porphyrins as destabilizer of a G-quadruplex located at the promoter of human MYH7 β gene

G-quadruplex (GQ) architecture is adopted by guanine rich sequences, present throughout the eukaryotic genome including promoter locations and telomeric ends. The in vivo presence indicates their involvement and role in various biological processes. Various small ligands have been developed to interact and stabilize/destabilize G-quadruplex structures. Cationic porphyrins are among the most studied ligands, reported to bind and stabilize G-quadruplexes. Herein, we report the recognition and destabilization of a parallel G-quadruplex by porphyrins (TMPyP3 and TMPyP4). This G-quadruplex forming 23-nt G-rich sequence is in the promoter region of Human Myosin Heavy Chain β gene (MYH7β). Presence of various putative regulatory sequence elements (TATA Box, CCAAT, SP-1) located in the vicinity of this quadruplex motif, highlight its regulatory implications. Biophysical methods as Circular Dichroism Spectroscopy, UV-Absorption Spectroscopy, UV-Thermal Denaturation and Fluorescence Spectroscopy (steady as well as Time Resolved) have been used for studying the interaction and binding parameters. It is proposed that porphyrins have a destabilizing effect on the G-quadruplexes with parallel topology and a stronger binding specifically via intercalation mode is needed to cause destabilization. The study deals with better understanding and insights of DNA-Drug interactions in biological systems.Communicated by Ramaswamy H. Sarma.

Radial distribution function of semiflexible oligomers with stretching flexibility

The radial distribution of the end-to-end distance R_ee is crucial for quantifying the global size and flexibility of a linear polymer. For semiflexible polymers, several analytical formulas have been derived for the radial distribution of R_ee ignoring the stretching flexibility. However, for semiflexible oligomers, such as DNA or RNA, the stretching flexibility can be rather pronounced and can significantly affect the radial distribution of R_ee. In this study, we obtained an extended formula that includes the stretch modulus to describe the distribution of R_ee for semiflexible oligomers on the basis of previous formulas for semiflexible polymers without stretching flexibility. The extended formula was validated by extensive Monte Carlo simulations over wide ranges of the stretch modulus and persistence length, as well as all-atom molecular dynamics simulations of short DNAs and RNAs. Additionally, our analyses showed that the effect of stretching flexibility on the distribution of R_ee becomes negligible for DNAs longer than ∼130 base pairs and RNAs longer than ∼240 base pairs.

Quantum mechanical investigation of G-quartet systems of DNA

Minima of the electric field and positions of K⁺ and Na⁺ (zero of the x-coordinate is the center of the cavity).