Triplex hydration: nanosecond molecular dynamics simulation of the solvated triplex formed by mixed sequences

Influence of hydrogen bonds on the reaction of guanine and hydroxyl radical: DFT calculations in C(H+)GC motif

A comprehensive theoretical investigation was performed to illuminate the influence of hydrogen bonds (H-bonds) on the obscure reaction of a hydroxyl radical (HO˙) and guanine (G) by selecting the building block of parallel triplex DNA, C(H+)GC, as the model. By mapping the energy profiles for addition and hydrogen abstraction reactions, the favorable pathway is predicted. The results reveal that in the C(H+)GC context, barrierless hydrogen abstraction from N2 of G leading to a neutral radical G(N2-H)˙ appears to become significant, but electrophilic attack by HO˙ on C8 of G resulting in 8-oxoG is the most thermodynamically favorable course. This shows a strong structural dependence due to the context constrained by the H-bond, which is dramatically different from the situation in unencumbered G. More interestingly, it proves that the stability order of resulting adduct radicals is not altered by H-bonding, but the activity for possible sites of the hydroxylation reaction changes. The significant influence of the H-bond on elementary reactions involved in the reaction is emphasized in the C(H+)GC context but is not restricted to the H-abstraction reaction. It is greatly anticipated that the present study could provide thoughtful insights into the vague hydroxyl radical-induced oxidation chemistry.

Selective Preference of Parallel DNA Triplexes Is Due to the Disruption of Hoogsteen Hydrogen Bonds Caused by the Severe Nonisostericity between the G*GC and T*AT Triplets

Implications of DNA, RNA and RNA.DNA hybrid triplexes in diverse biological functions, diseases and therapeutic applications call for a thorough understanding of their structure-function relationships. Despite exhaustive studies mechanistic rationale for the discriminatory preference of parallel DNA triplexes with G*GC & T*AT triplets still remains elusive. Here, we show that the highest nonisostericity between the G*GC & T*AT triplets imposes extensive stereochemical rearrangements contributing to context dependent triplex destabilisation through selective disruption of Hoogsteen scheme of hydrogen bonds. MD simulations of nineteen DNA triplexes with an assortment of sequence milieu reveal for the first time fresh insights into the nature and extent of destabilization from a single (non-overlapping), double (overlapping) and multiple pairs of nonisosteric base triplets (NIBTs). It is found that a solitary pair of NIBTs, feasible either at a G*GC/T*AT or T*AT/G*GC triplex junction, does not impinge significantly on triplex stability. But two overlapping pairs of NIBTs resulting from either a T*AT or a G*GC interruption disrupt Hoogsteen pair to a noncanonical mismatch destabilizing the triplex by ~10 to 14 kcal/mol, implying that their frequent incidence in multiples, especially, in short sequences could even hinder triplex formation. The results provide (i) an unambiguous and generalised mechanistic rationale for the discriminatory trait of parallel triplexes, including those studied experimentally (ii) clarity for the prevalence of antiparallel triplexes and (iii) comprehensive perspectives on the sequence dependent influence of nonisosteric base triplets useful in the rational design of TFO's against potential triplex target sites.

Multiscale analytic continuation approach to nanosystem simulation: applications to virus electrostatics

Electrostatic effects in nanosystems are understood via a physical picture built on their multiscale character and the distinct behavior of mobile ions versus charge groups fixed to the nanostructure. The Poisson-Boltzmann equation is nondimensionalized to introduce a factor lambda that measures the density of mobile ion charge versus that due to fixed charges; the diffusive smearing and volume exclusion effects of the former tend to diminish its value relative to that from the fixed charges. We introduce the ratio sigma of the average nearest-neighbor atom distance to the characteristic size of the features of the nanostructure of interest (e.g., a viral capsomer). We show that a unified treatment (i.e., lambda proportional to sigma) and a perturbation expansion around sigma=0 yields, through analytic continuation, an approximation to the electrostatic potential of high accuracy and computational efficiency. The approach was analyzed via Padé approximants and demonstrated on viral system electrostatics; it can be generalized to accommodate extended Poisson-Boltzmann models, and has wider applicability to nonequilibrium electrodiffusion and many-particle quantum systems.

Modulation of psoralen DNA crosslinking kinetics associated with a triplex-forming oligonucleotide

A triplex-forming oligonucleotide (TFO), HPRT3, conjugated to a psoralen derivative, was designed to target a psoralen reaction site within the HPRT gene. HPRT3 bound with high affinity to a synthetic duplex target sequence. At a uniform UVA radiation dose, the ratio of psoralen monoadducts (MA) to interstrand crosslinks decreased and inverted with increasing TFO concentration. As the TFO concentration increased from 10 nm to 10 microm, the efficiency of psoralen MA formation remained relatively constant but the efficiency of interstrand crosslink formation increased several-fold. Neither shortening the TFO to reduce its dissociation constant nor altering the DNA sequences flanking the TFO binding site altered the concentration dependence of MA and crosslink yields. The psoralen photokinetics associated with 10 nm HPRT3 converted to those associated with 10 microm HPRT3 with the addition of other unrelated TFOs at 10 microm that do not specifically interact with the HPRT3 target sequence. Glycerol at concentrations of 0.5% (vol/vol) or higher also mimicked high TFO concentrations in enhancing crosslink formation. These results demonstrate that while psoralen may be targeted to react at a particular sequence by TFOs, photoreactivity associated with triplex formation is also modulated by sequence-independent factors that may affect the local macromolecular environment.

Hydration and conformational transitions in DNA, RNA, and mixed DNA-RNA triplexes studied by gravimetry and FTIR spectroscopy

We have studied by gravimetric measurements and FTIR spectroscopy the hydration of duplexes and triplexes formed by combinations of dA(n), dT(n), rA(n), and rU(n) strands. Results obtained on hydrated films show important differences in their hydration and in the structural transitions which can be induced by varying the water content of the samples. The number of water molecules per nucleotide (w/n) measured at high relative humidity (98% R.H.) is found to be 21 for dA(n).dT(n) and 15 for rA(n).rU(n). Addition of a third rU(n) strand does not change the number of water molecules per nucleotide: w/n=21 for rU(n)*dA(n).dT(n) and w/n=15 for rU(n)*rA(n).rU(n). On the contrary, the addition of a third dT(n) strand changes the water content but in a different way, depending whether the duplex is DNA or RNA. Thus, a loss of four water molecules per nucleotide is measured for dT(n)*dA(n).dT(n) while an increase of two water molecules per nucleotide is observed for dT(n)*rA(n).rU(n). The final hydration is the same for both triplexes (w/n=17). The desorption profiles obtained by gravimetry and FTIR spectroscopy are similar for the rA(n).rU(n) duplex and the rU(n)*rA(n).rU(n) triplex. On the contrary, the desorption profiles of the dA(n).dT(n) duplex and the triplexes formed with it (rU(n)*dA(n).dT(n) and dT(n)*dA(n).dT(n)) are different from each other. This is correlated with conformational transitions induced by varying the hydration content of the different structures, as shown by FTIR spectroscopy. Modifications of the phosphate group hydration and of the sugar conformation (S to N type repuckering) induced by decrease of the water content are observed in the case of triplexes formed on the dA(n).dT(n) duplex.

Towards an understanding of DNA recognition by the methyl-CpG binding domain 1

CpG methylation determines a variety of biological functions of DNA. The methylation signal is interpreted by proteins containing a methyl-CpG binding domain (MBDs). Based on the NMR structure of MBD1 complexed with methylated DNA we analysed the recognition mode by means of molecular dynamics simulations. As the protein is monomeric and recognizes a symmetrically methylated CpG step, the recognition mode is an asymmetric one. We find that the two methyl groups do not contribute equally to the binding energy. One methyl group is associated with the major part of the binding energy and the other one nearly does not contribute at all. The contribution of the two cytosine methyl groups to binding energy is calculated to be -3.6 kcal/mol. This implies a contribution of greater than two orders of magnitude to the binding constant. The conserved amino acid Asp32 is known to be essential for DNA binding by MBD1, but so far no direct contact with DNA has been observed. We detected a direct DNA base contact to Asp32. This could be the main reason for the importance of this amino acid. MBD contacts DNA exclusively in the major groove, the minor groove is reserved for histone contacts. We found a deformation of the minor groove shape due to complexation by MBD1, which indicates an information transfer between the major and the minor groove.