On base flipping

Adenosine deaminases that act on RNA, then and now

In this article, I recount my memories of key experiments that led to my entry into the RNA editing/modification field. I highlight initial observations made by the pioneers in the ADAR field, and how they fit into our current understanding of this family of enzymes. I discuss early mysteries that have now been solved, as well as those that still linger. Finally, I discuss important, outstanding questions and acknowledge my hope for the future of the RNA editing/modification field.

Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond

Damaged or mismatched DNA bases are normally thought to be able to flip out of the helical stack, providing enzymes with access to the faulty genetic information otherwise hidden inside the helix. Thymine glycol (Tg) is one of the most common products of nucleic acid damage. However, the static and dynamic structures of DNA duplexes affected by 5R-Tg epimers are still not clearly understood, including the ability of these to undergo spontaneous base flipping. Structural effects of the 5R-Tg epimers on the duplex DNA are herein studied using molecular dynamics together with reliable DFT based calculations. In comparison with the corresponding intact DNA, the cis-5R,6S-Tg epimer base causes little perturbation to the duplex DNA, and a barrier of 4.9 kcal mol^–1 is obtained by meta-eABF for cis-5R,6S-Tg base flipping out of the duplex DNA, comparable to the 5.4 kcal mol^–1 obtained for the corresponding thymine flipping in intact DNA. For the trans-5R,6R-Tg epimer, three stable local structures were identified, of which the most stable disrupts the Watson–Crick hydrogen-bonded G5/C20 base pair, leading to conformational distortion of the duplex. Interestingly, the relative barrier height of the 5R-Tg flipping is only 1.0 kcal mol^–1 for one of these trans-5R,6R-Tg epimers. Water bridge interactions were identified to be essential for 5R-Tg flipping. The study clearly demonstrates the occurrence of partial trans-5R,6R-Tg epimer flipping in solution.

Structural Basis for the Interactions of the Colibactin Resistance Gene Product ClbS with DNA

The natural product colibactin, along with its associated biosynthetic gene cluster, is an example system for the role microbially derived small molecules play in the human microbiome. This is particularly relevant in the human gut, where host microbiota is involved in various disorders, including colorectal cancer pathogenesis. Bacteria harboring the colibactin gene cluster induce alkylation of nucleobases in host DNA, forming interstrand cross-links both in vivo and in vitro. These lesions can lead to deleterious double-strand breaks and have been identified as the primary mechanism of colibactin-induced cytotoxicity. The gene product ClbS is one of several mechanisms utilized by the producing bacteria to maintain genome integrity. ClbS catalyzes hydrolytic inactivation of colibactin and has been shown to bind DNA, incurring self-resistance. Presented is the molecular basis for ClbS bound to a DNA oligonucleotide. The structure shows the interaction of the protein with the ends of a DNA duplex with terminal nucleotides flipped to the enzyme active site. The structure suggests an additional function for ClbS, the binding to damaged DNA followed by repair. Additionally, our study provides general insight into the function of the widely distributed and largely uncharacterized DUF1706 protein family.

Duplexes Formed by G4C2 Repeats Contain Alternate Slow- and Fast-Flipping G·G Base Pairs

Aberrant expansion of the hexanucleotide GGGGCC (or G4C2) repeat in the human C9ORF72 gene is the most common genetic factor found behind amyotrophic lateral sclerosis and frontotemporal dementia. The hypothesized pathways, through which the repeat expansions contribute to the pathology, involve one or more secondary structural forms of the DNA and/or RNA sequences, such as G-quadruplexes, duplexes, and hairpins. Here, we study the structures of DNA and RNA duplexes formed by G4C2 repeats, which contain G(syn)·G(anti) base pairs flanked by either G·C or C·G base pairs. We show that duplexes formed by G4C2 repeats contain alternately two types of G·G pair contexts exhibiting different syn-anti base flipping dynamics (∼100 ms vs ∼2 ms for DNA and ∼50 ms vs ∼20 ms for RNA at 10 °C, respectively) depending on the flanking bases, with the slow-flipping G·G pairs being flanked by a guanine at the 5'-end and the fast-flipping G·G pairs being flanked by a cytosine at the 5'-end. Our findings on the structures and dynamics of G·G base pairs in DNA and RNA duplexes formed by G4C2 repeats provide a foundation for further studies of the functions and targeting of such biologically relevant motifs.

Binding Dynamics of siRNA with Selected Lipopeptides: A Computer-Aided Study of the Effect of Lipopeptides' Functional Groups and Stereoisomerism

The engineering issues pertaining to nanoparticle systems toward targeted gene therapies have not been fully probed. Recent experiments have identified specific structural characteristics of a novel class of lipopeptides (LP) that may lead to potent nanocarriers intended as RNAi therapeutics, albeit the molecular mechanism that underlies their performance remains unexplored. We conducted molecular dynamics simulations in atomistic detail coupled with free energy computations to study the dynamics and thermodynamics of an acrylate- and an epoxide-derived LP, members of the aforesaid class, upon their binding to siRNA in aqueous solution aiming at examining structure-potency relations. We found that the entropic part of the free energy of binding predominates; moreover, the first LP class tends to disrupt the Watson-Crick base pairing of siRNA, whereas the latter leaves the double helix intact. Moreover, the identified tug-of-war effect between LP-water and LP-siRNA hydrogen bonding in the supramolecular complex can underpin synthesis routes toward tuning the association dynamics. Our simulations on two diastereomers of the epoxide-derived LP showed significant structural and energetics differences upon binding, as a result of steric effects imposed by the different absolute configurations at their chiral centers. These findings may serve as crucial design parameters toward modulating the interplay between complex stability and ease of releasing the nucleic acid drug into the cell.

Sulfur-substitution-induced base flipping in the DNA duplex

Base flipping is widely observed in a number of important biological processes. The genetic codes deposited inside the DNA duplex become accessible to external agents upon base flipping. The sulfur substitution of guanine leads to thioguanine, which alters the thermodynamic stability of the GC base pairs and the GT mismatches. Experimental studies conclude that the sulfur substitution decreases the lifetime of the GC base pair. In this work, under three AMBER force fields for nucleotide systems, we firstly performed equilibrium and nonequilibrium free energy simulations to investigate the variation of the thermodynamic profiles in base flipping upon sulfur substitution. It is found that the bsc0 modification, the bsc1 modification and the OL15 modification of AMBER force fields are able to qualitatively describe the sulfur-substitution dependent behavior of the thermodynamics. However, only the two last-generation AMBER force fields are able to provide quantitatively correct predictions. The second computational study on the sulfur substitutions focused on the relative stability of the S6G-C base pair and the S6G-T mismatch. Two conflicting experimental observations were reported by the same authors. One suggested that the S6G-C base pair was more stable, while the other concludes that the S6G-T mismatch was more stable. We answered this question by constructing the free energy profiles along the base flipping pathway computationally.

A guanine-flipping and sequestration mechanism for G-quadruplex unwinding by RecQ helicases

Homeostatic regulation of G-quadruplexes (G4s), four-stranded structures that can form in guanine-rich nucleic acids, requires G4 unwinding helicases. The mechanisms that mediate G4 unwinding remain unknown. We report the structure of a bacterial RecQ DNA helicase bound to resolved G4 DNA. Unexpectedly, a guanine base from the unwound G4 is sequestered within a guanine-specific binding pocket. Disruption of the pocket in RecQ blocks G4 unwinding, but not G4 binding or duplex DNA unwinding, indicating its essential role in structure-specific G4 resolution. A novel guanine-flipping and sequestration model that may be applicable to other G4-resolving helicases emerges from these studies.

IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Studies of Sodium Cationized Thymidine and 5-Methyluridine: Kinetic Trapping During the ESI Desolvation Process Preserves the Solution Structure of [Thd+Na]<sup/>

Thymidine (dThd) is a fundamental building block of DNA nucleic acids, whereas 5-methyluridine (Thd) is a common modified nucleoside found in tRNA. In order to determine the conformations of the sodium cationized thymine nucleosides [dThd+Na]⁺ and [Thd+Na]⁺ produced by electrospray ionization, their infrared multiple photon dissociation (IRMPD) action spectra are measured. Complementary electronic structure calculations are performed to determine the stable low-energy conformations of these complexes. Geometry optimizations and frequency analyses are performed at the B3LYP/6-311+G(d,p) level of theory, whereas energies are calculated at the B3LYP/6-311+G(2d,2p) level of theory. As protonation preferentially stabilizes minor tautomers of dThd and Thd, tautomerization facilitated by Na⁺ binding is also considered. Comparisons of the measured IRMPD and computed IR spectra find that [dThd+Na]⁺ prefers tridentate (O2,O4',O5') coordination to the canonical 2,4-diketo form of dThd with thymine in a syn orientation. In contrast, [Thd+Na]⁺ prefers bidentate (O2,O2') coordination to the canonical 2,4-diketo tautomer of Thd with thymine in an anti orientation. Although 2,4-dihydroxy tautomers and O2 protonated thymine nucleosides coexist in the gas phase, no evidence for minor tautomers is observed for the sodium cationized species. Consistent with experimental observations, the computational results confirm that the sodium cationized thymine nucleosides exhibit a strong preference for the canonical form of the thymine nucleobase. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized dThd and Thd follow the order [dThd+H]⁺ < [Thd+H]⁺ < [dThd+Na]⁺ < [Thd+Na]⁺. Graphical Abstract ᅟ.

Evaluating the role of coherent delocalized phonon-like modes in DNA cyclization

The innate flexibility of a DNA sequence is quantified by the Jacobson-Stockmayer's J-factor, which measures the propensity for DNA loop formation. Recent studies of ultra-short DNA sequences revealed a discrepancy of up to six orders of magnitude between experimentally measured and theoretically predicted J-factors. These large differences suggest that, in addition to the elastic moduli of the double helix, other factors contribute to loop formation. Here, we develop a new theoretical model that explores how coherent delocalized phonon-like modes in DNA provide single-stranded "flexible hinges" to assist in loop formation. We combine the Czapla-Swigon-Olson structural model of DNA with our extended Peyrard-Bishop-Dauxois model and, without changing any of the parameters of the two models, apply this new computational framework to 86 experimentally characterized DNA sequences. Our results demonstrate that the new computational framework can predict J-factors within an order of magnitude of experimental measurements for most ultra-short DNA sequences, while continuing to accurately describe the J-factors of longer sequences. Further, we demonstrate that our computational framework can be used to describe the cyclization of DNA sequences that contain a base pair mismatch. Overall, our results support the conclusion that coherent delocalized phonon-like modes play an important role in DNA cyclization.

Transient Kinetic Methods for Mechanistic Characterization of DNA Binding and Nucleotide Flipping

Enzymes that modify nucleobases in double-stranded genomic DNA, either as part of a DNA repair pathway or as an epigenetic modifying pathway, adopt a multistep pathway to locate target sites and reconfigure the DNA to gain access. Work on several different enzymes has shown that in almost all cases base flipping, also known as nucleotide flipping, is a key feature of specific site recognition. In this chapter, we discuss some of the strategies that can be used to perform a kinetic characterization for DNA binding and nucleotide flipping. The resulting kinetic and thermodynamic framework provides a platform for understanding substrate specificity, mechanisms of inhibition, and the roles of important amino acids. We use a human DNA repair glycosylase called alkyladenine DNA glycosylase as a case study, because this is one of the best-characterized nucleotide-flipping enzymes. However, the approaches that are described can be readily adapted to study other enzymes, and future studies are needed to understand the mechanism of substrate recognition in each individual case. As more enzymes are characterized, we can hope to uncover which features of DNA searching and nucleotide flipping are fundamental features shared by many different families of DNA modifying enzymes, and which features are specific to a particular enzyme. Such an understanding provides reasonable models for less characterized enzymes that are important for epigenetic DNA modification and DNA repair pathways.

OGG1-DNA interactions facilitate NF-κB binding to DNA targets

DNA repair protein counteracting oxidative promoter lesions may modulate gene expression. Oxidative DNA bases modified by reactive oxygen species (ROS), primarily as 7, 8-dihydro-8-oxo-2′-deoxyguanosine (8-oxoG), which is repaired by 8-oxoguanine DNA glycosylase1 (OGG1) during base excision repair (BER) pathway. Because cellular response to oxidative challenge is accompanied by DNA damage repair, we tested whether the repair by OGG1 is compatible with transcription factor binding and gene expression. We performed electrophoretic mobility shift assay (EMSA) using wild-type sequence deriving from Cxcl2 gene promoter and the same sequence bearing a single synthetic 8-oxoG at defined 5′ or 3′ guanine in runs of guanines to mimic oxidative effects. We showed that DNA occupancy of NF-κB present in nuclear extracts from tumour necrosis factor alpha (TNFα) exposed cells is OGG1 and 8-oxoG position dependent, importantly, OGG1 counteracting 8-oxoG outside consensus motif had a profound influence on purified NF-κB binding to DNA. Furthermore, OGG1 is essential for NF-κB dependent gene expression, prior to 8-oxoG excised from DNA. These observations imply that pre-excision step(s) during OGG1 initiated BER evoked by ROS facilitates NF-κB DNA occupancy and gene expression.

Biosynthesis and Function of Modified Bases in Bacteria and Their Viruses

Naturally occurring modification of the canonical A, G, C, and T bases can be found in the DNA of cellular organisms and viruses from all domains of life. Bacterial viruses (bacteriophages) are a particularly rich but still underexploited source of such modified variant nucleotides. The modifications conserve the coding and base-pairing functions of DNA, but add regulatory and protective functions. In prokaryotes, modified bases appear primarily to be part of an arms race between bacteriophages (and other genomic parasites) and their hosts, although, as in eukaryotes, some modifications have been adapted to convey epigenetic information. The first half of this review catalogs the identification and diversity of DNA modifications found in bacteria and bacteriophages. What is known about the biogenesis, context, and function of these modifications are also described. The second part of the review places these DNA modifications in the context of the arms race between bacteria and bacteriophages. It focuses particularly on the defense and counter-defense strategies that turn on direct recognition of the presence of a modified base. Where modification has been shown to affect other DNA transactions, such as expression and chromosome segregation, that is summarized, with reference to recent reviews.

Twenty Years of Delila and Molecular Information Theory: The Altenberg-Austin Workshop in Theoretical Biology Biological Information, Beyond Metaphor: Causality, Explanation, and Unification Altenberg, Austria, 11-14 July 2002

A brief personal history is given about how information theory can be applied to binding sites of genetic control molecules on nucleic acids. The primary example used is ribosome binding sites in Escherichia coli. Once the sites are aligned, the information needed to describe the sites can be computed using Claude Shannon's method. This is displayed by a computer graphic called a sequence logo. The logo represents an average binding site, and the mathematics easily allows one to determine the components of this average. That is, given a set of binding sites, the information for individual binding sites can also be computed. One can go further and predict the information of sites that are not in the original data set. Information theory also allows one to model the flexibility of ribosome binding sites, and this led us to a simple model for ribosome translational initiation in which the molecular components fit together only when the ribosome is at a good ribosome binding site. Since information theory is general, the same mathematics applies to human splice junctions, where we can predict the effect of sequence changes that cause human genetic diseases and cancer. The second example given is the Pribnow 'box' which, when viewed by the information theory method, reveals a mechanism for initiation of both transcription and DNA replication. Replication, transcription, splicing, and translation into protein represent the central dogma, so these examples show how molecular information theory is contributing to our knowledge of basic biology.

Reversibly locked thionucleobase pairs in DNA to study base flipping enzymes

Covalently interstrand cross-linked DNA is an interesting tool to study DNA binding proteins that locally open up the DNA duplex by flipping single bases out of the DNA helix or melting whole stretches of base pairs to perform their function. The ideal DNA cross-link to study protein–DNA interactions should be specific and easy to synthesize, be stable during protein binding experiments, have a short covalent linker to avoid steric hindrance of protein binding, and should be available as a mimic for both A/T and G/C base pairs to cover all possible binding specificities. Several covalent interstrand cross-links have been described in the literature, but most of them fall short of at least one of the above criteria. We developed an efficient method to site-specifically and reversibly cross-link thionucleoside base pairs in synthetic duplex oligodeoxynucleotides by bisalkylation with 1,2-diiodoethane resulting in an ethylene-bridged base pair. Both linked A/T and G/C base pair analogs can conveniently be prepared which allows studying any base pair-opening enzyme regardless of its sequence specificity. The cross-link is stable in the absence of reducing agents but the linker can be quickly and tracelessly removed by the addition of thiol reagents like dithiothreitol. This property makes the cross-linking reaction fully reversible and allows for a switching of the linked base pair from locked to unlocked during biochemical experiments. Using the DNA methyltransferase from Thermus aquaticus (M.TaqI) as example, we demonstrate that the presented cross-linked DNA with an ethylene-linked A/T base pair analog at the target position is a useful tool to determine the base-flipping equilibrium constant of a base-flipping enzyme which lies mostly on the extrahelical side for M.TaqI.

Hidden Conformation Events in DNA Base Extrusions: A Generalized Ensemble Path Optimization and Equilibrium Simulation Study

DNA base extrusion is a crucial component of many biomolecular processes. Elucidating how bases are selectively extruded from the interiors of double-strand DNAs is pivotal to accurately understanding and efficiently sampling this general type of conformational transitions. In this work, the on-the-path random walk (OTPRW) method, which is the first generalized ensemble sampling scheme designed for finite-temperature-string path optimizations, was improved and applied to obtain the minimum free energy path (MFEP) and the free energy profile of a classical B-DNA major-groove base extrusion pathway. Along the MFEP, an intermediate state and the corresponding transition state were located and characterized. The MFEP result suggests that a base-plane-elongation event rather than the commonly focused base-flipping event is dominant in the transition state formation portion of the pathway; and the energetic penalty at the transition state is mainly introduced by the stretching of the Watson-Crick base pair. Moreover to facilitate the essential base-plane-elongation dynamics, the surrounding environment of the flipped base needs to be intimately involved. Further taking the advantage of the extended-dynamics nature of the OTPRW Hamiltonian, an equilibrium generalized ensemble simulation was performed along the optimized path; and based on the collected samples, several base-flipping (opening) angle collective variables were evaluated. In consistence with the MFEP result, the collective variable analysis result reveals that none of these commonly employed flipping (opening) angles alone can adequately represent the base extrusion pathway, especially in the pre-transition-state portion. As further revealed by the collective variable analysis, the base-pairing partner of the extrusion target undergoes a series of in-plane rotations to facilitate the base-plane-elongation dynamics. A base-plane rotation angle is identified to be a possible reaction coordinate to represent these in-plane rotations. Notably, these in-plane rotation motions may play a pivotal role in determining the base extrusion selectivity.

Bullied no more: when and how DNA shoves proteins around

The predominant protein-centric perspective in protein-DNA-binding studies assumes that the protein drives the interaction. Research focuses on protein structural motifs, electrostatic surfaces and contact potentials, while DNA is often ignored as a passive polymer to be manipulated. Recent studies of DNA topology, the supercoiling, knotting, and linking of the helices, have shown that DNA has the capability to be an active participant in its transactions. DNA topology-induced structural and geometric changes can drive, or at least strongly influence, the interactions between protein and DNA. Deformations of the B-form structure arise from both the considerable elastic energy arising from supercoiling and from the electrostatic energy. Here, we discuss how these energies are harnessed for topology-driven, sequence-specific deformations that can allow DNA to direct its own metabolism.

Early stage intercalation of doxorubicin to DNA fragments observed in molecular dynamics binding simulations

The intercalation mode between doxorubicin (an anticancer drug) and two 6-base-pair DNA model fragments (d(CGATCG)₂ and d(CGTACG)₂) has been well studied by X-ray crystallography and NMR experimental methods. Yet, the detailed intercalation pathway at molecular level remains elusive. In this study, we conducted molecular dynamics binding simulations of these two systems using AMBER DNA (parmbsc0) and drug (GAFF) force fields starting from the unbound state. We observed outside binding (minor groove binding or end-binding) in all six independent binding simulations (three for each DNA fragment), followed by the complete intercalation of a drug molecule in two simulations (one for each DNA fragment). First, our data directly supported that the minor groove binding is the dominant pre-intercalation step. Second, we observed that the opening and flipping of a local base pair (A3-T10 for d(CGATCG)₂ and C1-G12 for d(CGTACG)₂) in the two intercalation trajectories. This locally cooperative flipping-intercalation mechanism was different from the previously proposed rise-insertion mechanism by which the distance between two neighboring intact base pairs increases to create a space for the drug insertion. Third, our simulations provided the first set of data to support the applicability of the AMBER DNA and drug force fields in drug-DNA atomistic binding simulations. Implications on the kinetics pathway and drug action are also discussed.

Conformational diversity of short DNA duplex

Two 25 base-pair cDNA strands are encapsulated within an optically trapped nanodroplet, and we observe the kinetics of their hybridization in dynamic equilibrium via single-molecule fluorescence resonance energy transfer (FRET) as a function of temperature and of the solution's NaCl concentration. We have observed the duplex unfolding and refolding, and we have observed quasistable partially unfolded states under low salinity conditions. Furthermore, our measurements reveal that, even in conditions under which the duplex is stable, it undergoes conformational fluctuations in solution.

Opening of DNA double helix at room temperature: Application of alpha-cyclodextrin self-aggregates

Self-aggregation of alpha-cyclodextrin (alpha-CD) can induce DNA opening at room temperature (25 degrees C) owing to the hydroxyl groups on the surface of the spherical aggregates of alpha-CD, which promote hydrogen bonding with the flipped-out bases in DNA duplex prohibiting them from reverting back.

Fluorescence properties of pyrimidopyrimidoindole nucleoside dC(PPI) incorporated into oligodeoxynucleotides

A series of oligodeoxynucleotides labeled by a pyrimidopyrimidoindole deoxynucleoside (1a: dC(PPI)) and its derivatives 2a and 3a substituted with electron-donating and -withdrawing groups, respectively, were synthesized according to the phosphoramidite approach. The photophysical properties and quenching efficiencies of oligonucleotides incorporating dC(PPI) derivatives were studied in detail. The thermal denaturation experiments and molecular dynamics simulation of DNA duplexes incorporating dC(PPI) suggested that a modified base of dC(PPI) could form base pairs with guanine and adenine in canonical Watson-Crick and reverse-wobble geometries, respectively. The fluorescence of oligonucleotides incorporating dC(PPI) derivatives increased upon binding to the counter strands, except when dC(PPI) and guanine formed a base pair. It was revealed that dGMP quenched the fluorescence of the cyano derivative 3a most effectively, whereas it affected that of the methoxy derivative 2a least effectively. The involvement of the electron transfer from guanine to the dC(PPI) derivatives in the fluorescence quenching was supported by energy considerations.

In the absence of writhe, DNA relieves torsional stress with localized, sequence-dependent structural failure to preserve B-form

To understand how underwinding and overwinding the DNA helix affects its structure, we simulated 19 independent DNA systems with fixed degrees of twist using molecular dynamics in a system that does not allow writhe. Underwinding DNA induced spontaneous, sequence-dependent base flipping and local denaturation, while overwinding DNA induced the formation of Pauling-like DNA (P-DNA). The winding resulted in a bimodal state simultaneously including local structural failure and B-form DNA for both underwinding and extreme overwinding. Our simulations suggest that base flipping and local denaturation may provide a landscape influencing protein recognition of DNA sequence to affect, for examples, replication, transcription and recombination. Additionally, our findings help explain results from single-molecule experiments and demonstrate that elastic rod models are strictly valid on average only for unstressed or overwound DNA up to P-DNA formation. Finally, our data support a model in which base flipping can result from torsional stress.

Damage detection and base flipping in direct DNA alkylation repair

THE FOREIGN LESION: The mechanistic questions for DNA base damage detection by repair proteins are discussed in this Minireview. Repair proteins could either probe and locate a weakened base pair that results from base damage, or passively capture an extrahelical base lesion in the first step of damage searching on double-stranded DNA. How some repair proteins, such as AGT (see figure), locate base lesions in DNA is still not fully understood.To remove a few damaged bases efficiently from the context of the entire genome, the DNA base repair proteins rely on remarkably specific detection mechanisms to locate base lesions. This efficient molecular recognition event inside cells has been extensively studied with various structural and biochemical tools. These studies suggest that DNA base damage can be located by repair proteins by using two mechanisms: a repair protein can probe and detect a weakened base pair that results from mutagenic or cytotoxic base damage; alternatively, a protein can passively capture and stabilize an extrahelical base lesion. Our chemical and structural studies on the direct DNA repair proteins hAGT, C-Ada and ABH2 suggest that these proteins search for weakened base pairs in their first step of damage searching. We have also discovered a very unique base-flipping mechanism used by the DNA repair protein AlkB. This protein distorts DNA and favors single stranded DNA (ssDNA) substrates over double-stranded (dsDNA) ones. Potentially, it locates base lesions in dsDNA by imposing a constraint that targets less rigid regions of the duplex DNA. The exact mechanism of how AlkB and related proteins search for damage in ssDNA and dsDNA still awaits further studies.

Structural basis for recognition of hemi-methylated DNA by the SRA domain of human UHRF1

Epigenetic inheritance in mammals is characterized by high-fidelity replication of CpG methylation patterns during development. UHRF1 (also known as ICBP90 in humans and Np95 in mouse) is an E3 ligase important for the maintenance of global and local DNA methylation in vivo. The preferential affinity of UHRF1 for hemi-methylated DNA over symmetrically methylated DNA by means of its SET and RING-associated (SRA) domain and its association with the maintenance DNA methyltransferase 1 (DNMT1) suggests a role in replication of the epigenetic code. Here we report the 1.7 A crystal structure of the apo SRA domain of human UHRF1 and a 2.2 A structure of its complex with hemi-methylated DNA, revealing a previously unknown reading mechanism for methylated CpG sites (mCpG). The SRA-DNA complex has several notable structural features including a binding pocket that accommodates the 5-methylcytosine that is flipped out of the duplex DNA. Two specialized loops reach through the resulting gap in the DNA from both the major and the minor grooves to read the other three bases of the CpG duplex. The major groove loop confers both specificity for the CpG dinucleotide and discrimination against methylation of deoxycytidine of the complementary strand. The structure, along with mutagenesis data, suggests how UHRF1 acts as a key factor for DNMT1 maintenance methylation through recognition of a fundamental unit of epigenetic inheritance, mCpG.

Molecular dynamics of a kappaB DNA element: base flipping via cross-strand intercalative stacking in a microsecond-scale simulation

The sequence-dependent structural variability and conformational dynamics of DNA play pivotal roles in many biological milieus, such as in the site-specific binding of transcription factors to target regulatory elements. To better understand DNA structure, function, and dynamics in general, and protein···DNA recognition in the ‘κB’ family of genetic regulatory elements in particular, we performed molecular dynamics simulations of a 20-bp DNA encompassing a cognate κB site recognized by the proto-oncogenic ‘c-Rel’ subfamily of NF-κB transcription factors. Simulations of the κB DNA in explicit water were extended to microsecond duration, providing a broad, atomically detailed glimpse into the structural and dynamical behavior of double helical DNA over many timescales. Of particular note, novel (and structurally plausible) conformations of DNA developed only at the long times sampled in this simulation—including a peculiar state arising at ≈0.7 μs and characterized by cross-strand intercalative stacking of nucleotides within a longitudinally sheared base pair, followed (at ≈1 μs) by spontaneous base flipping of a neighboring thymine within the A-rich duplex. Results and predictions from the microsecond-scale simulation include implications for a dynamical NF-κB recognition motif, and are amenable to testing and further exploration via specific experimental approaches that are suggested herein.

Germ-line mutations, DNA damage, and global hypermethylation in mice exposed to particulate air pollution in an urban/industrial location

Particulate air pollution is widespread, yet we have little understanding of the long-term health implications associated with exposure. We investigated DNA damage, mutation, and methylation in gametes of male mice exposed to particulate air pollution in an industrial/urban environment. C57BL/CBA mice were exposed in situ to ambient air near two integrated steel mills and a major highway, alongside control mice breathing high-efficiency air particulate (HEPA) filtered ambient air. PCR analysis of an expanded simple tandem repeat (ESTR) locus revealed a 1.6-fold increase in sperm mutation frequency in mice exposed to ambient air for 10 wks, followed by a 6-wk break, compared with HEPA-filtered air, indicating that mutations were induced in spermatogonial stem cells. DNA collected after 3 or 10 wks of exposure did not exhibit increased mutation frequency. Bulky DNA adducts were below the detection threshold in testes samples, suggesting that DNA reactive chemicals do not reach the germ line and cause ESTR mutation. In contrast, DNA strand breaks were elevated at 3 and 10 wks, possibly resulting from oxidative stress arising from exposure to particles and associated airborne pollutants. Sperm DNA was hypermethylated in mice breathing ambient relative to HEPA-filtered air and this change persisted following removal from the environmental exposure. Increased germ-line DNA mutation frequencies may cause population-level changes in genetic composition and disease. Changes in methylation can have widespread repercussions for chromatin structure, gene expression and genome stability. Potential health effects warrant extensive further investigation.

6MAP, a fluorescent adenine analogue, is a probe of base flipping by DNA photolyase

Cyclobutylpyrimidine dimers (CPDs) are formed between adjacent pyrimidines in DNA when it absorbs ultraviolet light. CPDs can be directly repaired by DNA photolyase (PL) in the presence of visible light. How PL recognizes and binds its substrate is still not well understood. Fluorescent nucleic acid base analogues are powerful probes of DNA structure. We have used the fluorescent adenine analogue 6MAP, a pteridone, to probe the local double helical structure of the CPD substrate when bound by photolyase. Duplex melting temperatures were obtained by both UV-vis absorption and fluorescence spectroscopies to ascertain the effect of the probe and the CPD on DNA stability. Steady-state fluorescence measurements of 6MAP-containing single-stranded and doubled-stranded oligos with and without protein show that the local region around the CPD is significantly disrupted. 6MAP shows a different quenching pattern compared to 2-aminopurine, another important adenine analogue, although both probes show that the structure of the complementary strand opposing the 5'-side of the CPD lesion is more destacked than that opposing the 3'-side in substrate/protein complexes. We also show that 6MAP/CPD duplexes are substrates for PL. Vertical excitation energies and transition dipole moment directions for 6MAP were calculated using time-dependent density functional theory. Using these results, the Förster resonance energy transfer efficiency between the individual adenine analogues and the oxidized flavin cofactor was calculated to account for the observed intensity pattern. These calculations suggest that energy transfer is highly efficient for the 6MAP probe and less so for the 2Ap probe. However, no experimental evidence for this process was observed in the steady-state emission spectra.

Base-flipping dynamics in a DNA hairpin processing reaction

Many enzymes that repair or modify bases in double-stranded DNA gain access to their substrates by base flipping. Although crystal structures provide stunning snap shots, biochemical approaches addressing the dynamics have proven difficult, particularly in complicated multi-step reactions. Here, we use protein-DNA crosslinking and potassium permanganate reactivity to explore the base-flipping step in Tn5 transposition. We present a model to suggest that base flipping is driven by a combination of factors including DNA bending and the intrusion of a probe residue. The forces are postulated to act early in the reaction to create a state of tension, relieved by base flipping after cleavage of the first strand of DNA at the transposon end. Elimination of the probe residue retards the kinetics of nicking and reduces base flipping by 50%. Unexpectedly, the probe residue is even more important during the hairpin resolution step. Overall, base flipping is pivotal to the hairpin processing reaction because it performs two opposite but closely related functions. On one hand it disrupts the double helix, providing the necessary strand separation and steric freedom. While on the other, transposase appears to position the second DNA strand in the active site for cleavage using the flipped base as a handle.

Differential distortion of substrate occurs when it binds to DNA photolyase: a 2-aminopurine study

Cyclobutylpyrimidine dimers (CPDs) are formed between adjacent pyrimidines in DNA when it is exposed to ultraviolet light. CPDs can be directly repaired by DNA photolyase (PL) upon absorption of blue-green light. We have used the fluorescent adenine analogue 2-aminopurine (2Ap) to probe the local double-helical structure of the DNA substrate when it binds to the protein. Duplex melting temperatures and van't Hoff enthalpies were obtained by both UV-vis absorption and fluorescence spectroscopies to ascertain the effect of the probe and CPD on DNA stability. Steady-state fluorescence measurements of the single- and double-stranded oligos showed that the local region around the 5'-side of the CPD lesion was more disrupted and destacked than the 3'-side in substrate-protein complexes. These results were compared with those of a protein-substrate crystal structure, demonstrating that the crystal structure and solution-state studies are in agreement with regard to the differential distortions of the target DNA at the active site of the protein.

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.

Sequence-dependent base pair opening in DNA double helix

Preservation of genetic information in DNA relies on shielding the nucleobases from damage within the double helix. Thermal fluctuations lead to infrequent events of the Watson-Crick basepair opening, or DNA "breathing", thus making normally buried groups available for modification and interaction with proteins. Fluctuational basepair opening implies the disruption of hydrogen bonds between the complementary bases and flipping of the base out of the helical stack. Prediction of sequence-dependent basepair opening probabilities in DNA is based on separation of the two major contributions to the stability of the double helix: lateral pairing between the complementary bases and stacking of the pairs along the helical axis. The partition function calculates the basepair opening probability at every position based on the loss of two stacking interactions and one base-pairing. Our model also includes a term accounting for the unfavorable positioning of the exposed base, which proceeds through a formation of a highly constrained small loop, or a ring. Quantitatively, the ring factor is found as an adjustable parameter from the comparison of the theoretical basepair opening probabilities and the experimental data on short DNA duplexes measured by NMR spectroscopy. We find that these thermodynamic parameters suggest nonobvious sequence dependent basepair opening probabilities.

NMR imino proton exchange experiments on duplex DNA primarily monitor the opening of purine bases

Molecular dynamics simulations were performed to investigate GC and AT base opening events in DNA. Calculated equilibrium constants between the base open (or flipped) and closed states were shown to be in good agreement with experimental data from NMR imino proton exchange experiments. Analysis of the computed results indicates that the equilbrium constants are dominated by the opening of the A and G bases in the AT and GC base pairs, respectively. Thus, the present results predict that NMR imino proton exchange experiments of base opening are primarily monitoring the opening of purine bases.

Crystal structure of thymine DNA glycosylase conjugated to SUMO-1

Members of the small ubiquitin-like modifier (SUMO) family can be covalently attached to the lysine residue of a target protein through an enzymatic pathway similar to that used in ubiquitin conjugation, and are involved in various cellular events that do not rely on degradative signalling via the proteasome or lysosome. However, little is known about the molecular mechanisms of SUMO-modification-induced protein functional transfer. During DNA mismatch repair, SUMO conjugation of the uracil/thymine DNA glycosylase TDG promotes the release of TDG from the abasic (AP) site created after base excision, and coordinates its transfer to AP endonuclease 1, which catalyses the next step in the repair pathway. Here we report the crystal structure of the central region of human TDG conjugated to SUMO-1 at 2.1 A resolution. The structure reveals a helix protruding from the protein surface, which presumably interferes with the product DNA and thus promotes the dissociation of TDG from the DNA molecule. This helix is formed by covalent and non-covalent contacts between TDG and SUMO-1. The non-covalent contacts are also essential for release from the product DNA, as verified by mutagenesis.

Excitation, Ionization, and Fragmentation of Chiral Molecules in Asymmetric Microenvironments: A Mass-Resolved R2PI Spectroscopic Study

One- and two-color, mass selected R2PI spectra of the S(1) <-- S(0) transitions in the bare (R)-(+)-1-phenyl-1-propanol and its complexes with bidentate solvent molecules, like the (R)-(-)- and (S)-(+)-3-hydroxytetrahydrofuran enantiomers, have been recorded after a supersonic molecular beam expansion. The one-color R2PI excitation spectra of the diastereomeric complexes are characterized by three main peaks, one red-shifted and the other two blue-shifted relative to the band origin of the most stable anti conformer of the bare chromophore. The opposite direction of these spectral shifts is ascribed to the occurrence of three different hydrogen bonded isomeric structures for each individual complex, while their different magnitude depends on the configuration of the bidentate solvent molecule as well as its specific hydrogen bond interaction center, whether the ethereal oxygen atom or the hydroxyl group. The same factors play a major role in determining the magnitude of the phenomenological activation barriers for the loss of an ethyl radical from the ionized diastereomeric complexes.

The Two-Photon Excitation Cross Section of 6MAP, a Fluorescent Adenine Analogue

6MAP is a fluorescent analogue of adenine that undergoes Watson-Crick base pairing and base stacking in double-stranded DNA. The one-photon absorption spectrum of 6MAP is characterized by a maximum around 330 nm with moderate quantum yield fluorescence centered at about 420 nm. To take advantage of this probe for confocal and single-molecule microscopy, it would be advantageous to be able to excite the analogue via two photons. We report the first determination of the two-photon excitation cross section and spectrum for 6MAP from 614 to 700 nm. The power dependence of the fluorescence indicates that emission results from the absorption of two photons. The one-photon and two-photon emission line shapes are identical within experimental error. A study of the concentration dependence of the fluorescence yield for one-photon excitation shows no measurable quenching up to about 5 microM. The maximum in the two-photon excitation spectrum gives a two-photon cross section, delta(TPE), of 3.4 +/- 0.1 Goeppert-Mayer (G.M.) at 659 nm, which correlates well with the one-photon absorption maximum. This compares quite favorably with cross sections of various naturally fluorescent biological molecules such as flavins and nicotiamide. In addition, we have also obtained the two-photon-induced fluorescence emission spectrum of quinine sulfate. It is approximately the same as that for one-photon excitation, suggesting that two-photon excitation of quinine sulfate may be used for calibration purposes.

Dynamic opening of DNA during the enzymatic search for a damaged base

Uracil DNA glycosylase (UDG) removes uracil from U.A or U.G base pairs in genomic DNA by extruding the aberrant uracil from the DNA base stack. A question in enzymatic DNA repair is whether UDG and related glycosylases also use an extrahelical recognition mechanism to inspect the integrity of undamaged base pairs. Using NMR imino proton exchange measurements we find that UDG substantially increases the equilibrium constant for opening of T-A base pairs by almost two orders of magnitude relative to free B-DNA. This increase is brought about by enzymatic stabilization of an open state of the base pair without increasing the rate constant for spontaneous base pair opening. These findings indicate a passive search mechanism in which UDG uses the spontaneous opening dynamics of DNA to inspect normal base pairs in a rapid genome-wide search for uracil in DNA.

Atomistic view of base flipping in DNA

Base flipping is essential for the enzyme-catalysed methylation of DNA. In our previous studies, the flipping of bases out of duplex DNA in DNA alone and when bound to the (cytosine-C5)-methyltransferase from HhaI (M.HhaI) were investigated via potential of mean force calculations. Insights into various experimental observations were obtained. In the present paper we present an overview of previous computational studies of base flipping along with new detailed structural and energetic analysis on atomic events that contribute to the free energy surfaces. The contributions from different intrinsic and environmental effects to the base-flipping process are explored, and experimental data derived from a variety of methods are reconciled. A detailed protein-facilitated base-flipping mechanism is proposed. Ground-state destabilization is achieved via disruption of the target base Watson-Crick interactions by substitution with favourable DNA-protein interactions. In addition, specific DNA-protein interactions and favourable solvation effects further promote target base flipping along the major groove through the protein matrix, and maximal interactions occur between the DNA and the protein upon reaching the fully flipped state. Other DNA binding proteins that involve distortion of DNA's conformation may use a similar mechanism to that by which M.HhaI facilitates base flipping.

Asynchronous basepair openings in transcription initiation: CRP enhances the rate-limiting step

The mechanism of isomerization (basepair openings) during transcription initiation by RNA polymerase at the galP1 promoter of Escherichia coli was investigated by 2-aminopurine (2,AP) fluorescence. The fluorescence of 2,AP is quenched in DNA duplex and enhanced when the basepair is distorted or deformed. The increase of 2,AP fluorescence was used to monitor basepair distortion at several individual positions in the promoter. We observed that basepair distortions during isomerization are a multi-step process. Three distinct hitherto unresolved steps in kinetic terms were observed, where significant fluorescence change occurs: a fast step with a half-life of around 1 s, which is followed by two slower steps occurring with a half-life in the range of minutes at 25 degrees C. Contrary to commonly held expectations, basepairs at different positions opened by 2,AP assays without any obvious pattern, suggesting that basepair opening is an asynchronous multi-step process. cAMP.CRP, which activates transcription at galP1, enhanced the rate-limiting step.

Mechanisms of DNA damage, DNA hypomethylation, and tumor progression in the folate/methyl-deficient rat model of hepatocarcinogenesis

Using the folate/methyl-deficient rat model of hepatocarcinogenesis, we obtained evidence that may provide new insights into a major unresolved paradox in DNA methylation and cancer research: the mechanistic basis for genome-wide hypomethylation despite an increase in DNA methyltransferase activity and gene-specific regional hypermethylation. Previous studies revealed that the methyltransferase binds with higher affinity to DNA strand breaks, gaps, abasic sites, and uracil than it does to its cognate hemimethylated CpG sites, consistent with its ancestral function as a DNA repair enzyme. These same DNA lesions are an early occurrence in models of folate and methyl deficiency and are often present in human preneoplastic cells. We hypothesized that the high-affinity binding of the maintenance DNA methyltransferase to unrepaired lesions in DNA could sequester available enzyme away from the replication fork and promote passive replication-dependent demethylation. In support of this possibility, we found that lesion-containing DNA is less efficiently methylated than lesion-free DNA from folate/methyl-deficient rats and that an increase in DNA strand breaks precedes DNA hypomethylation. Despite an adaptive increase in DNA methyltransferase activity, hemimethylated DNA from folate/methyl-deficient rats is progressively replaced by double-stranded unmethylated DNA that is resistant to remethylation with dietary methyl repletion. In promoter regions, the inappropriate binding of the DNA methyltransferase to unrepaired lesions or mispairs may promote local histone deacetylation, methylation, and regional hypermethylation associated with tumor suppressor gene silencing. These insights in an experimental model are consistent with the possibility that DNA lesions may be a necessary prerequisite for the disruption of normal DNA methylation patterns in preneoplastic and neoplastic cells.

Speculations about alternative DNA structures

An alternative model to the Watson & Crick (W&C) double DNA-spiral and the Pauling & Corey (P&C) triple spiral is presented. In this model: (1). the rotation axis of the polynucleotide chain is in the ribose ring; (2). there is a -H- bond or direct covalent bond between the O2 (PO(4) and C2(') (in ribose) which makes the nucleic acid strands 'stiff'; (3). when there is a covalent bond between O2 and C2('), the unit of the DNA is the ribonucleoside 2('), 3(')-cyclic monophosphate, an intermediate form between DNA and RNA; (4). the bases point outwards from the rotation axis and may interact with each other to connect 2-4 strands together through complementary base pairs; (5). two strands may, but do not necessarily, form a helical structure and if they do, the interacting strands do not turn around each other. The architecture of this model, termed the Homulus DNA model is open (in contrast to the inverted W&C model) and using it might help us to understand the nature of some specific DNA-protein interactions, ordered chromatin formation (coiling and de-coiling), specific gene-to-gene interaction (gene targeting). It is possible that a small portion of the total DNA, the transcribed, 'working DNA', might be built by this way.

Theoretical ab initio study of the effects of methylation on structure and stability of G:C Watson-Crick base pair

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6-31G(d) basis set was used. The single-point calculations have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.

Detection and analysis of enzymatic DNA methylation of oligonucleotide substrates by matrix-assisted laser desorption ionization time-of-flight mass spectrometry

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF) mass spectrometry was employed to analyze DNA methylation carried out by the Escherichia coli dam DNA methyltransferase using oligonucleotide substrates with molecular masses of 5000-10,000 Da per strand. The mass spectrometry assay offers several advantages: (i) it directly shows the methylation as the increase in the mass of the substrate DNA, (ii) it is nonradioactive, (iii) it is quantitative, and (iv) it can be automated for high-throughput applications. Since unmethylated and methylated DNA are detected, the ratio of methylation can be determined directly and accurately. Furthermore, the assay allows detection individually of the methylation of several substrates in competition, offering an ideal setup to analyze the specificity of DNA interacting with enzymes. We could not identify methylation at any noncanonical site, indicating that the dam MTase is a very specific enzyme. Finally, MALDI-TOF mass spectrometry permitted assessment of the number of methyl groups incorporated into each DNA strand, thereby, allowing study of mechanistic details such as the processivity of the methylation reaction. We provide evidence that the dam MTase modifies DNA in a processive reaction, confirming earlier findings.