Statistical Mechanical Approach for Predicting the Transition to Non-B DNA Structures in Supercoiled DNA

Abstract

Supercoiling causes global twist of DNA structure and the supercoiled state has wide influence on conformational transition. A statistical mechanical approach was made for prediction of the transition probability to non-B DNA structures under torsional stress. A conditional partition function was defined as the sum over all possible states of the DNA sequence with basepair 1 and basepair n being in B-form helix and a recurrence formula was developed which expressed the partition function for basepair n with those for less number of pairs. This new definition permits a quick enumeration of every configuration of secondary structures. Energetic parameters of all conformations concerned, involving B-form, interior loop, cruciform and Z-form, were included in the equation. The probability of transition to each non-B conformation could be derived from these conditional partition functions. For treatment of effects of superhelicity, supercoiling energy was considered, and a twist of each conformation was determined to minimize the supercoiling energy. As the twist itself affects the transition probability, the whole scheme of equations was solved by renormalization technique. The present method permits a simultaneous treatment of serveral types of conformations under a common torsional stress.

A set of energetic parameters of DNA secondary structures has been chosen for calculation. Some DNA sequences were submitted to the calculation, and all the sequences that we submitted gave stable convergence. Some of them have been investigated the critical supercoil density for the transition to non-B DNA structures. Even though the reliability of the set of parameters was not enough, the prediction of secondary structure transition showed good agreement with reported observation. Hence, the present algorithm can estimate the probability of local conformational change of DNA under a given supercoil density, and also be employed to predict some specific sequences in which conformational change is sensitive to superhelicity.     

Mitochondrial DNA deletions are associated with non-B DNA conformations

Mitochondrial DNA (mtDNA) deletions are a primary cause of mitochondrial disease and are believed to contribute to the aging process and to various neurodegenerative diseases. Despite strong observational and experimental evidence, the molecular basis of the deletion process remains obscure. In this study, we test the hypothesis that the primary cause of mtDNA vulnerability to breakage resides in the formation of non-B DNA conformations, namely hairpin, cruciform and cloverleaf-like elements. Using the largest database of human mtDNA deletions built thus far (753 different cases), we show that site-specific breakage hotspots exist in the mtDNA. Furthermore, we discover that the most frequent deletion breakpoints occur within or near predicted structures, a result that is supported by data from transgenic mice with mitochondrial disease. There is also a significant association between the folding energy of an mtDNA region and the number of breakpoints that it harbours. In particular, two clusters of hairpins (near the D-loop 3'-terminus and the L-strand origin of replication) are hotspots for mtDNA breakage. Consistent with our hypothesis, the highest number of 5'- and 3'-breakpoints per base is found in the highly structured tRNA genes. Overall, the data presented in this study suggest that non-B DNA conformations are a key element of the mtDNA deletion process.     

Base Dependence of B-DNA Sugar Conformation in Solution and in the Solid State

Abstract

Analysis of 1H-NOESY solution data for eight short DNA duplexes has revealed pronounced differences between the sugar conformations of purine and pyrimidine nucleotides. It was found that the H1′-H4′ interproton distance is less than ca. 3.0 A. in pyrimidine sugars, while in purine sugars it is more than ca. 3.OA. This difference has been analyzed by comparison with the sugar conformations of highly resolved B-DNA crystal structures and model sugar conformations. The conclusion can be drawn that the deoxyribose conformation is of the general C2′-endo type but pyrimidine sugars are characterized by smaller phase angles of pseudorotation P (90°<P<150°), while purine sugars have larger P values that are greater than ca. 140° (140°<P<180°). There is no such clear base dependence of sugar conformation in highly resolved B-DNA crystal structures; however the similar trend can be seen as in the solution studies. Based on B-type DNA crystal structures, I-coupling constants have been calculated, and the applicability of experimental coupling measurements to the determination of sugar conformation is discussed.     

Sequence-specific transitions of the torsion angle gamma change the polar-hydrophobic profile of the DNA grooves: implication for indirect protein–DNA recognition

Variations of the shape and polarity of the DNA grooves caused by changes of the DNA conformation play an important role in the DNA readout. Despite the fact that non-canonical trans and gauche- conformations of the DNA backbone angle γ (O5′–C5′–C4′–C3′) are frequently found in the DNA crystal structures, their possible role in the DNA recognition has not been studied systematically. In order to fill in this gap, we analyze the available high-resolution crystal structures of the naked and complexed DNA. The analysis shows that the non-canonical γ angle conformations are present both in the naked and bound DNA, more often in the bound vs. naked DNA, and in the nucleotides with the A-like vs. the B-like sugar pucker. The alternative angle γ torsions are more frequently observed in the purines with the A-like sugar pucker and in the pyrimidines with the B-like sugar conformation. The minor groove of the nucleotides with non-canonical γ angle conformation is more polar, while the major groove is more hydrophobic than in the nucleotides with the classical γ torsions due to variations in exposure of the polar and hydrophobic groups of the DNA backbone. The propensity of the nucleotides with different γ angle conformations to participate in the protein–nucleic acid contacts in the minor and major grooves is connected with their sugar pucker and sequence-specific. Our findings imply that the angle γ transitions contribute to the process of the protein–DNA recognition due to modification of the polar/hydrophobic profile of the DNA grooves.     

Molecular Dynamics Simulations of Chlorambucil/DNA Adducts. A Structural Basis for the 5′ -GNC Interstrand DNA Crosslink Formed by Nitrogen Mustards

Abstract

The alkylation of DNA by chlorambucil has been studied using a computational approach. Molecular dynamics simulations were performed on the fully solvated non-covalent complex, two monoadducts and a crosslinked diadduct of chlorambucil with the d(CGG₃G₂CGC).- d(GCG₁CCCG) duplex, in which the N7 atoms of G₁, G₂ and G₃ are potential alkylation sites. The results provide a structural basis for the preference of nitrogen mustards to crosslink DNA duplexes at a 5′-GNC site (a 1,3 crosslink, G₁ -G₃) rather than at a 5′-GC sites (a 1,2 crosslink, G₁ -G₂).

In the non-covalent complex simulation the drug reoriented from a non-interstrand crosslinking location to a position favorable for G₁ -G₃ diadduct formation. It proved possible to construct a G₁ -G₃ diadduct from a structure from the non-covalent simulation, and continue the molecular dynamics calculation without further disruption of the DNA structure. A crosslinked diadduct developed with four B_II conformations on the 3′ side of each alkylated guanine and of their respective complementary cytosine. In the first monoadduct simulation the starting point was the same DNA conformation used in the crosslinked diadduct simulation with alkylation at G₁. In this simulation the DNA deformation was reduced, with the helix returning to a more canonical form. A second monoadduct simulation was started from a canonical DNA conformation alkylated at G₃. Here, no significant motion towards a potential crosslinking conformation occurred. Collectively, the results suggest that crosslink formation is dependent upon the drug orientation prior to alkylation and the required deformation of the DNA to permit 1,3 crosslinking can largely be achieved in the non-covalent complex.     

Non-B DNA conformations formed by long repeating tracts of myotonic dystrophy type 1, myotonic dystrophy type 2, and Friedreich's ataxia genes, not the sequences per se, promote mutagenesis in flanking regions

The expansions of long repeating tracts of CTG.CAG, CCTG.CAGG, and GAA.TTC are integral to the etiology of myotonic dystrophy type 1 (DM1), myotonic dystrophy type 2 (DM2), and Friedreich's ataxia (FRDA). Essentially all studies on the molecular mechanisms of this expansion process invoke an important role for non-B DNA conformations which may be adopted by these repeat sequences. We have directly evaluated the role(s) of the repeating sequences per se, or of the non-B DNA conformations formed by these sequences, in the mutagenic process. Studies in Escherichia coli and three types of mammalian (COS-7, CV-1, and HEK-293) fibroblast-like cells revealed that conditions which promoted the formation of the non-B DNA structures enhanced the genetic instabilities, both within the repeat sequences and in the flanking sequences of up to approximately 4 kbp. The three strategies utilized included: the in vivo modulation of global negative supercoil density using topA and gyrB mutant E. coli strains; the in vivo cleavage of hairpin loops, which are an obligate consequence of slipped-strand structures, cruciforms, and intramolecular triplexes, by inactivation of the SbcC protein; and by genetic instability studies with plasmids containing long repeating sequence inserts that do, and do not, adopt non-B DNA structures in vitro. Hence, non-B DNA conformations are critical for these mutagenesis mechanisms.     

Non-B DNA conformations, mutagenesis and disease

Recent discoveries have revealed that simple repeating DNA sequences, which are known to adopt non-B DNA conformations (such as triplexes, cruciforms, slipped structures, left-handed Z-DNA and tetraplexes), are mutagenic. The mutagenesis is due to the non-B DNA conformation rather than to the DNA sequence per se in the orthodox right-handed Watson-Crick B-form. The human genetic consequences of these non-B structures are approximately 20 neurological diseases, approximately 50 genomic disorders (caused by gross deletions, inversions, duplications and translocations), and several psychiatric diseases involving polymorphisms in simple repeating sequences. Thus, the convergence of biochemical, genetic and genomic studies has demonstrated a new paradigm implicating the non-B DNA conformations as the mutagenesis specificity determinants, not the sequences as such.     

On the Bioactive Conformation of a Small Peptide and its Set of Thermodynamically Accessible Conformations

Abstract

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics framework using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.     

Modulation of DNA structure formation using small molecules

Genome integrity is essential for proper cell function such that genetic instability can result in cellular dysfunction and disease. Mutations in the human genome are not random, and occur more frequently at “hotspot” regions that often co-localize with sequences that have the capacity to adopt alternative (i.e. non-B) DNA structures. Non-B DNA-forming sequences are mutagenic, can stimulate the formation of DNA double-strand breaks, and are highly enriched at mutation hotspots in human cancer genomes. Thus, small molecules that can modulate the conformations of these structure-forming sequences may prove beneficial in the prevention and/or treatment of genetic diseases. Further, the development of molecular probes to interrogate the roles of non-B DNA structures in modulating DNA function, such as genetic instability in cancer etiology are warranted. Here, we discuss reported non-B DNA stabilizers, destabilizers, and probes, recent assays to identify ligands, and the potential biological applications of these DNA structure-modulating molecules.     

Dependence of Deformability of Geometries and Characteristics of Intramolecular Hydrogen Bonds in Canonical 2′-Deoxyribonucleotides on DNA Conformations

Abstract

The molecular structure and deformability (with respect to average geometry) of methyl ethers of canonical 2′-deoxyribonucleotides thymidine-5′-phosphate (mTMP), 2-deoxycytidine-5′- phosphate (mCMP), 2-deoxyadenosine-5′-phosphate (mAMP) and 2′-deoxyguanosine-5′- phosphate (mGMP) in different types of DNA have been calculated using B3LYP/cc-pvdz method. Comparison of energy at equilibrium conformations of nucleotides and conformations with torsion angles of backbone fixed to average values for different types of DNA reveals that incorporation of nucleotides to A-DNA macromolecules requires the minimum amount of deformation energy. Therefore, this type of DNA should be the least strained from viewpoint of intramolecular deformations of monomers. Modeling of environmental effects within the PCM approach reveals that the immersion of nucleotides in polar medium results in significant decrease of energy differences between anti conformers of all DNTs and syn conformers of mGMP This also leads to reduction by almost a half nucleotides' deformation energy facilitating formation of DNA macromolecule. Change of DNTs conformation causes switch between different types of intramolecular H bonds. Every type of DNA possesses unique set of intramolar hydrogen bonds in nucleotides.     

Explicit Solvent Molecular Dynamics Simulation of Duplex Formed by the Modified Oligonucleotide with Alternating Phosphate/Phosphonate Internucleoside Linkages and its Natural Counterpart

Abstract

Impact of the internucleoside linkage modification by inserting a methylene group on the ability of the modified oligonucleotide to hybridize with a natural DNA strand was studied by fully solvated molecular dynamics (MD) simulations. Three undecamer complexes were analyzed: natural dT₁₁.dA₁₁ duplex as a reference and two its analogs with alternating modified and natural linkages in the deoxyadenosine chain. The isopolar, non-isosteric modified linkages were of 5′-O-PO₂-CH₂-O-3′ (5′PC3′) or 5′-O-CH₂-PO₂-O-3′ (5′CP3′) type. Simulations were performed by using the AMBER 5.0 software package with the force field completed by a set of parameters needed to model the modified segments. Both modifications were found to lead to double helical complexes, in which the thymidine strand as well as deoxyriboses and unmodified linkages in the adenosine strand adopted conformations typical for the B-type structure. For each of the two conformational richer modified linkages two stable conformations were found at 300 K: the -ggt and ggt for the 5′PC3′ and ggg, tgg for the 5′CP3′, respectively. Both modified chains adopted helical conformations with heightened values of the inclination parameter but without affecting the Watson-Crick hydrogen bonds.     

A Molecular Mechanics and Dynamics Study of Alternate Triple-Helices Involving the Integrase-Binding Site of the HIV-1 Virus and Oligonucleotides Having a 3′-3′ Internucleotide Junction

Abstract

Triple helix formation by oligonucleotides can be extended beyond polypurine tracts with the help of specially designed linkers. In this paper we focus our attention on the integrase- binding site of the HIV-1 virus located on the U5 LTR end which contains two adjacent purine tracts on opposite strands. Two alternate triple helices with a 3′-3′ junction in the third strand are considered:

5′-GGTTTTp3′-3′pTGTGT-5′ and 5′-GGAAAAp3′-3′pAGAGA-5′

The structural plausibility of these triplexes is investigated using molecular mechanics and dynamics simulations, both in vacuo and in aqua.

The non-isomorphism of the triplets in the GpT steps in the first sequence, gives rise to non canonical conformations in the torsion angles, hydration appears to be crucial for this triplex. Sugar puckers are predominantly South during in vacuo simulations while they tum East in aqua. In the simulation in aqua the triplexes are shrouded by an hydration shell, however, we have not been able to detect any permanent hydrogen bond bridge between DNA and water.

The solvation of ions as well as their radial distribution, appear to be relatively well behaved despite the artifacts known to be generated by the simulation procedure. The experimental feasibility of these structures is discussed.     

Non-B DNA structure-induced genetic instability

Repetitive DNA sequences are abundant in eukaryotic genomes, and many of these sequences have the potential to adopt non-B DNA conformations. Genes harboring non-B DNA structure-forming sequences increase the risk of genetic instability and thus are associated with human diseases. In this review, we discuss putative mechanisms responsible for genetic instability events occurring at these non-B DNA structures, with a focus on hairpins, left-handed Z-DNA, and intramolecular triplexes or H-DNA. Slippage and misalignment are the most common events leading to DNA structure-induced mutagenesis. However, a number of other mechanisms of genetic instability have been proposed based on the finding that these structures not only induce expansions and deletions, but can also induce DNA strand breaks and rearrangements. The available data implicate a variety of proteins, such as mismatch repair proteins, nucleotide excision repair proteins, topoisomerases, and structure specific-nucleases in the processing of these mutagenic DNA structures. The potential mechanisms of genetic instability induced by these structures and their contribution to human diseases are discussed.     

Molecular Dynamics Study of the Opening Mechanism for DNA Polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme∶DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.     

Machine Learning Estimates of Natural Product Conformational Energies

Machine learning has been used for estimation of potential energy surfaces to speed up molecular dynamics simulations of small systems. We demonstrate that this approach is feasible for significantly larger, structurally complex molecules, taking the natural product Archazolid A, a potent inhibitor of vacuolar-type ATPase, from the myxobacterium Archangium gephyra as an example. Our model estimates energies of new conformations by exploiting information from previous calculations via Gaussian process regression. Predictive variance is used to assess whether a conformation is in the interpolation region, allowing a controlled trade-off between prediction accuracy and computational speed-up. For energies of relaxed conformations at the density functional level of theory (implicit solvent, DFT/BLYP-disp3/def2-TZVP), mean absolute errors of less than 1 kcal/mol were achieved. The study demonstrates that predictive machine learning models can be developed for structurally complex, pharmaceutically relevant compounds, potentially enabling considerable speed-ups in simulations of larger molecular structures.     

Molecular Conformation of 2′-Deoxy-2′-methylidene-cytidine: A Potent Antineoplastic Nucleoside

Abstract

2′-Deoxy-2′-methylidenecytidine (DMDC), a potent inhibitor of the growth of tumor cells, was crystallized with two different forms. One is dihydrated (DMDC·2H₂O) and the other is its hydrochloride salt (DMDC·HCLl). Both crystal and molecular structures have been determined by the X-ray diffraction method. In both forms the glycosidic and sugar conformations are anti and C(4′)-exo, respectively, whereas the conformation about the exocyclic bond is trans for DMDC·2H₂O and gauche ⁺ for DMDC·HCl. Proton nuclear magnetic resonance data of DMDC indicate a preference for the anti C(4′)-exo conformation found in the solid state. These molecular conformations were compared with the related pyrimidine nucleosides. When the cytosine bases are brought into coincidence, DMDC displays the exocyclic C(4′)-C(5′) bond located on the very close position to those of pyrimidine nucleosides with typical overall conformations. On the other hand, the hydroxyl O(3′)-H groups are separated by ca. 3 Å in the cases of DMDC and other pyrimidine nucleosides which have the C(2′)-endo sugar conformation. This result may be useful for the implication about the mechanism of the biological activity of DMDC.     

The Molecular Determinants of the Increased Reduction Potential of the Rubredoxin Domain of Rubrerythrin Relative to Rubredoxin

Based on the crystal structures, three possible sequence determinants have been suggested as the cause of a 285 mV increase in reduction potential of the rubredoxin domain of rubrerythrin over rubredoxin by modulating the polar environment around the redox site. Here, electrostatic calculations of crystal structures of rubredoxin and rubrerythrin and molecular dynamics simulations of rubredoxin wild-type and mutants are used to elucidate the contributions to the increased reduction potential. Asn¹⁶⁰ and His¹⁷⁹ in rubrerythrin versus valines in rubredoxins are predicted to be the major contributors, as the polar side chains contribute significantly to the electrostatic potential in the redox site region. The mutant simulations show both side chains rotating on a nanosecond timescale between two conformations with different electrostatic contributions. Reduction also causes a change in the reduction energy that is consistent with a linear response due to the interesting mechanism of shifting the relative populations of the two conformations. In addition to this, a simulation of a triple mutant indicates the side-chain rotations are approximately anticorrelated so whereas one is in the high potential conformation, the other is in the low potential conformation. However, Ala¹⁷⁶ in rubrerythrin versus a leucine in rubredoxin is not predicted to be a large contributor, because the solvent accessibility increases only slightly in mutant simulations and because it is buried in the interface of the rubrerythrin homodimer.