Computational DNA binding studies of (–)-epigallocatechin-3-gallate

<p>The catechin family of molecules that are present in the leaves of green tea has been under investigation since the antioxidant and anti-inflammatory properties of tea were discovered. Among multiple proposed therapeutic targets of these molecules, the direct interaction with nucleic acids has been proposed and experimentally observed but without clear knowledge about the potential binding modes between these ligands and DNA. One of these catechin structures, (–)-epigallocatechin gallate (EGCG), has three aromatic rings that could interact with double-stranded DNA via terminal base-pair stacking, intercalation, or through groove binding. Using enhanced sampling techniques and molecular dynamics simulations, we have found a stable complex between the EGCG ligand and DNA through intercalation of the trihydroxybenzoate aromatic ring and an ApC step. Moreover, we have calculated the absorption spectra of four possible binding modes and compared these to absorption profiles reported in the literature, and explored the possible DNA sequence preference for the EGCG ligand to bind. Our results suggest that an intercalative mode of interaction through the major groove is possible between the EGCG ligands and DNA with apparently very little DNA sequence selectivity.</p

Investigating the ion dependence of the first unfolding step of GTPase-Associating Center ribosomal RNA

<p>The interactions in the tertiary structure of a ribosomal RNA fragment in the GTPase Associating Center (<i>GAC</i>) have been experimentally studied, but the roles of the bound and diffuse cations in its folding pathway have not yet been fully elucidated. Melting experiments have shown that the temperature of the first of the two distinguishable transitions in the unfolding pathway of the <i>GAC</i> RNA can be regulated by altering the magnesium concentration, yet the physical interpretation of such ion-dependent effects on folding have not been clearly understood in spite of the availability of crystal structures that depict many <i>GAC</i> RNA–ion interactions. Here, we use umbrella sampling and molecular dynamics (MD) simulations to provide a physical description for the first transition in this unfolding pathway, with a focus on the role of a chelated magnesium ion. Our results indicate that the presence of cations mediating the local interaction of two loops stabilizes the folded state relative to the unfolded or partially folded states. Also, our findings suggest that a bridging magnesium ion between the two loops improves the stabilizing effect. This is consistent with the multistep unfolding pathway proposed for the <i>GAC</i> RNA and highlights the importance of ions in the first unfolding step. The results suggest how MD simulations can provide insight into RNA unfolding pathways as a complementary approach to experiments.</p

Insight into G-DNA Structural Polymorphism and Folding from Sequence and Loop Connectivity through Free Energy Analysis

The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence–structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT₂, dT₃, dTTA, and dT₄] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K⁺ solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure

Exploring potentially alternative non-canonical DNA duplex structures through simulation

<p>Hopkins proposed an alternative and chirally distinct family of double-stranded DNA (dsDNA) models that have antiparallel chains with 5′→3′ senses opposite to those of the right-handed Watson–Crick (WC) family. Termed configuration II, this family of dsDNA models contains both right-handed (II-R) and left-handed (II-L) forms, with Z-DNA as an example of the latter. Relative interstrand binding energies for six DNA duplex models, two each of configuration I-R (standard WC canonical B-DNA), II-R, and II-L for the duplex d(CGCGAATTCGCG), have been estimated under identical conditions using MM-PBSA analysis from molecular dynamics trajectories using three different AMBER force fields. These simulations support the stereo chemical soundness of configuration II dsDNA forms. Recent force fields (Barcelona Supercomputing Center [BSC] [bsc1] and Olomouc 2015 [OL15]) successfully render stable II-L structures, whereas the previous force field, bsc0, generated stable II-R structures, although with an energy difference between II-R and II-L of ∼30 kcal/mol.</p> <p>Communicated by Ramaswamy H. Sarma</p