Base-pairing of uracil and 2,6-diaminopurine: from cocrystals to photoreactivity
Explore the source record for details and available documents.
Tomislav Stolar↗
SEARCH · Search NASA
Results for “Base Pairing”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Statistical and linguistic features of DNA sequences
We present evidence supporting the idea that the DNA sequence in genes containing noncoding regions is correlated, and that the correlation is remarkably long range--indeed, base pairs thousands of base pairs distant are correlated. We do not find such a long-range correlation in the coding regions of the gene. We resolve the problem of the "non-stationary" feature of the sequence of base pairs by applying a new algorithm called Detrended Fluctuation Analysis (DFA). We address the claim of Voss that there is no difference in the statistical properties of coding and noncoding regions of DNA by systematically applying the DFA algorithm, as well as standard FFT analysis, to all eukaryotic DNA sequences (33 301 coding and 29 453 noncoding) in the entire GenBank database. We describe a simple model to account for the presence of long-range power-law correlations which is based upon a generalization of the classic Levy walk. Finally, we describe briefly some recent work showing that the noncoding sequences have certain statistical features in common with natural languages. Specifically, we adapt to DNA the Zipf approach to analyzing linguistic texts, and the Shannon approach to quantifying the "redundancy" of a linguistic text in terms of a measurable entropy function. We suggest that noncoding regions in plants and invertebrates may display a smaller entropy and larger redundancy than coding regions, further supporting the possibility that noncoding regions of DNA may carry biological information.
Non-NASA Center↗
Structure of an anti-HIV-1 hammerhead ribozyme complex with a 17-mer DNA substrate analog of HIV-1 gag RNA and a mechanism for the cleavage reaction: 750 MHz NMR and computer experiments
The structure of an anti-HIV-1 ribozyme-DNA abortive substrate complex was investigated by 750 MHz NMR and computer modeling experiments. The ribozyme was a chimeric molecule with 30 residues-18 DNA nucleotides, and 12 RNA residues in the conserved core. The DNA substrate analog had 17 residues. The chimeric ribozyme and the DNA substrate formed a shortened ribozyme-abortive substrate complex of 47 nucleotides with two DNA stems (stems I and III) and a loop consisting of the conserved core residues. Circular dichroism spectra showed that the DNA stems assume A-family conformation at the NMR concentration and a temperature of 15 degrees C, contrary to the conventional wisdom that DNA duplexes in aqueous solution populate entirely in the B-form. It is proposed that the A-family RNA residues at the core expand the A-family initiated at the core into the DNA stems because of the large free energy requirement for the formation of A/B junctions. Assignments of the base H8/H6 protons and H1' of the 47 residues were made by a NOESY walk. In addition to the methyl groups of all T's, the imino resonances of stems I and III and AH2's were assigned from appropriate NOESY walks. The extracted NMR data along with available crystallographic data, were used to derive a structural model of the complex. Stems I and III of the final model displayed a remarkable similarity to the A form of DNA; in stem III, a GC base pair was found to be moving into the floor of the minor groove defined by flanking AT pairs; data suggest the formation of a buckled rhombic structure with the adjacent pair; in addition, the base pair at the interface of stem III and the loop region displayed deformed geometry. The loop with the catalytic core, and the immediate region of the stems displayed conformational multiplicity within the NMR time scale. A catalytic mechanism for ribozyme action based on the derived structure, and consistent with biochemical data in the literature, is proposed. The complex between the anti HIV-1 gag ribozyme and its abortive DNA substrate manifests in the detection of a continuous track of A.T base pairs; this suggests that the interaction between the ribozyme and its DNA substrate is stronger than the one observed in the case of the free ribozyme where the bases in stem I and stem III regions interact strongly with the ribozyme core region (Sarma, R. H., et al. FEBS Letters 375, 317-23, 1995). The complex formation provides certain guidelines in the design of suitable therapeutic ribozymes. If the residues in the ribozyme stem regions interact with the conserved core, it may either prevent or interfere with the formation of a catalytically active tertiary structure.
NASA Discipline Exobiology↗
Nonenzymatic RNA copying with a potentially primordial genetic alphabet
Nonenzymatic RNA copying is thought to have been responsible for the replication of genetic information during the origin of life. However, chemical copying with the canonical nucleotides (A, U, G, and C) strongly favors the incorporation of G and C and disfavors the incorporation of A and especially U because of the stronger G:C vs. A:U base pair and the weaker stacking interactions of U. Recent advances in prebiotic chemistry suggest that the 2-thiopyrimidines were precursors to the canonical pyrimidines, raising the possibility that they may have played an important early role in RNA copying chemistry. Furthermore, 2-thiouridine (s 2 U) and inosine (I) form by deamination of 2-thiocytidine (s 2 C) and A, respectively. We used thermodynamic and crystallographic analyses to compare the I:s 2 C and A:s 2 U base pairs. We find that the I:s 2 C base pair is isomorphic and isoenergetic with the A:s 2 U base pair. The I:s 2 C base pair is weaker than a canonical G:C base pair, while the A:s 2 U base pair is stronger than the canonical A:U base pair, so that a genetic alphabet consisting of s 2 U, s 2 C, I, and A generates RNA duplexes with uniform base pairing energies. Consistent with these results, kinetic analysis of nonenzymatic template-directed primer extension reactions reveals that s 2 C and s 2 U substrates bind similarly to I and A in the template, and vice versa. Our work supports the plausibility of a potentially primordial genetic alphabet consisting of s 2 U, s 2 C, I, and A and offers a potential solution to the long-standing problem of biased nucleotide incorporation during nonenzymatic template copying.
Science & Technology - Other Topics↗
Scaling features of noncoding DNA
We review evidence supporting the idea that the DNA sequence in genes containing noncoding regions is correlated, and that the correlation is remarkably long range--indeed, base pairs thousands of base pairs distant are correlated. We do not find such a long-range correlation in the coding regions of the gene, and utilize this fact to build a Coding Sequence Finder Algorithm, which uses statistical ideas to locate the coding regions of an unknown DNA sequence. Finally, we describe briefly some recent work adapting to DNA the Zipf approach to analyzing linguistic texts, and the Shannon approach to quantifying the "redundancy" of a linguistic text in terms of a measurable entropy function, and reporting that noncoding regions in eukaryotes display a larger redundancy than coding regions. Specifically, we consider the possibility that this result is solely a consequence of nucleotide concentration differences as first noted by Bonhoeffer and his collaborators. We find that cytosine-guanine (CG) concentration does have a strong "background" effect on redundancy. However, we find that for the purine-pyrimidine binary mapping rule, which is not affected by the difference in CG concentration, the Shannon redundancy for the set of analyzed sequences is larger for noncoding regions compared to coding regions.
Non-NASA Center↗
Targeting tRNA-synthetase interactions towards novel therapeutic discovery against eukaryotic pathogens
The development of chemotherapies against eukaryotic pathogens is especially challenging because of both the evolutionary conservation of drug targets between host and parasite, and the evolution of strain-dependent drug resistance. There is a strong need for new nontoxic drugs with broad-spectrum activity against trypanosome parasites such as Leishmania and Trypanosoma. A relatively untested approach is to target macromolecular interactions in parasites rather than small molecular interactions, under the hypothesis that the features specifying macromolecular interactions diverge more rapidly through coevolution. We computed tRNA Class-Informative Features in humans and independently in eight distinct clades of trypanosomes, identifying parasite-specific informative features, including base pairs and base mis-pairs, that are broadly conserved over approximately 250 million years of trypanosome evolution. Validating these observations, we demonstrated biochemically that tRNA:aminoacyl-tRNA synthetase (aaRS) interactions are a promising target for anti-trypanosomal drug discovery. From a marine natural products extract library, we identified several fractions with inhibitory activity toward Leishmania major alanyl-tRNA synthetase (AlaRS) but no activity against the human homolog. These marine natural products extracts showed cross-reactivity towards Trypanosoma cruzi AlaRS indicating the broad-spectrum potential of our network predictions. We also identified Leishmania major threonyl-tRNA synthetase (ThrRS) inhibitors from the same library. We discuss why chemotherapies targeting multiple aaRSs should be less prone to the evolution of resistance than monotherapeutic or synergistic combination chemotherapies targeting only one aaRS.
59 BASIC BIOLOGICAL SCIENCES↗
Identification of a TAAT-containing motif required for high level expression of the COL1A1 promoter in differentiated osteoblasts of transgenic mice
Our previous studies have shown that the 49-base pair region of promoter DNA between -1719 and -1670 base pairs is necessary for transcription of the rat COL1A1 gene in transgenic mouse calvariae. In this study, we further define this element to the 13-base pair region between -1683 and -1670. This element contains a TAAT motif that binds homeodomain-containing proteins. Site-directed mutagenesis of this element in the context of a COL1A1-chloramphenicol acetyltransferase construct extending to -3518 base pairs decreased the ratio of reporter gene activity in calvariae to tendon from 3:1 to 1:1, suggesting a preferential effect on activity in calvariae. Moreover, chloramphenicol acetyltransferase-specific immunofluorescence microscopy of transgenic calvariae showed that the mutation preferentially reduced levels of chloramphenicol acetyltransferase protein in differentiated osteoblasts. Gel mobility shift assays demonstrate that differentiated osteoblasts contain a nuclear factor that binds to this site. This binding activity is not present in undifferentiated osteoblasts. We show that Msx2, a homeodomain protein, binds to this motif; however, Northern blot analysis revealed that Msx2 mRNA is present in undifferentiated bone cells but not in fully differentiated osteoblasts. In addition, cotransfection studies in ROS 17/2.8 osteosarcoma cells using an Msx2 expression vector showed that Msx2 inhibits a COL1A1 promoter-chloramphenicol acetyltransferase construct. Our results suggest that high COL1A1 expression in bone is mediated by a protein that is induced during osteoblast differentiation. This protein may contain a homeodomain; however, it is distinct from homeodomain proteins reported previously to be present in bone.
NASA Discipline Cell Biology↗
Green’s function coupled cluster simulation of the near-valence ionizations of DNA-fragments
Accurate description of the ionization process in DNA is crucial to the understanding of the DNA damage under exposure to ionizing radiation and the exploration of the potential application of DNA strands in nanoelectronics. In this work, by employing our recently developed Green’s function coupled-cluster library on supercomputing facilities, we have studied the spectral functions of several guanine–cytosine (G–C) base pair structures ([G–C]n, n = 1–3) for the first time in a relatively broad near-valence regime ([−25.0, −5.0] eV) in the coupled-cluster with singles and doubles level. Our focus is to give a preliminary many-body coupled-cluster understanding and guideline of the vertical ionization energy (VIE), spectral profile, and ionization feature changes of these systems as the system size expands in this near-valence regime. The results show that, as the system size expands, even though the lowest VIEs keep decreasing, the changes of spectral function profile and the relative peak positions get unexpectedly smaller. Further analysis of the ionized states associated with the most intensive peak in the spectral functions reveals non-negligible |2h, 1p⟩’s in the ionized wave functions of the considered G–C base pair systems. The leading |2h, 1p⟩’s associated with the main ionizations from the cytosine part of the G–C base pairs feature a transition from the intra-base-pair cytosine π → π* excitation to the inter-base-pair electron excitation as the size of G–C base pairs expands, which also indicates the minimum quantum region in the many-body calculations of DNA systems.
Peng, Bo (ORCID:0000000242267294)↗
Structure of native four-repeat satellite III sequence with non-canonical base interactions
Abstract Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5′-ATGGA ATGGA ATGGA ATGGA-3′ using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05–0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G–G intercalation sandwiched by sheared A–G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A–T/T–A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions.
59 BASIC BIOLOGICAL SCIENCES↗
Revisiting the Hole Size in Double Helical DNA with Localized Orbital Scaling Corrections
The extent of electronic wave function delocalization for the charge carrier (electron or hole) in double helical DNA plays an important role in determining the DNA charge transfer mechanism and kinetics. The size of the charge carrier’s wave function delocalization is regulated by the solvation induced localization and the quantum delocalization among the π stacked base pairs at any instant of time. Using a newly developed localized orbital scaling correction (LOSC) density functional theory method, we accurately characterized the quantum delocalization of the hole wave function in double helical B-DNA. This approach can be used to diagnose the extent of delocalization in fluctuating DNA structures. Our studies indicate that the hole state tends to delocalize among 4 guanine–cytosine (GC) base pairs and among 3 adenine–thymine (AT) base pairs when these adjacent bases fluctuate into degeneracy. The relatively small delocalization in AT base pairs is caused by the weaker π–π interaction. This extent of delocalization has significant implications for assessing the role of coherent, incoherent, or flickering coherent carrier transport in DNA.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Crystal Structure of an i‐Motif from the HRAS Oncogene Promoter
Abstract An i‐motif is a non‐canonical DNA structure implicated in gene regulation and linked to cancers. The C‐rich strand of the HRAS oncogene, 5′‐CG CCC GTG CCC TGCG CCC GCAA CCC GA‐3′ (herein referred to as iHRAS), forms an i‐motif in vitro but its exact structure was unknown. HRAS is a member of the RAS proto‐oncogene family. About 19 % of US cancer patients carry mutations in RAS genes. We solved the structure of iHRAS at 1.77 Å resolution. The structure reveals that iHRAS folds into a double hairpin. The two double hairpins associate in an antiparallel fashion, forming an i‐motif dimer capped by two loops on each end and linked by a connecting region. Six C−C + base pairs form each i‐motif core, and the core regions are extended by a G−G base pair and a cytosine stacking. Extensive canonical and non‐canonical base pairing and stacking stabilizes the connecting region and loops. The iHRAS structure is the first atomic resolution structure of an i‐motif from a human oncogene. This structure sheds light on i‐motifs folding and function in the cell.
Li, Kevin S.↗
Crystal Structure of an i‐Motif from the HRAS Oncogene Promoter
Abstract An i‐motif is a non‐canonical DNA structure implicated in gene regulation and linked to cancers. The C‐rich strand of the HRAS oncogene, 5′‐CG CCC GTG CCC TGCG CCC GCAA CCC GA‐3′ (herein referred to as iHRAS), forms an i‐motif in vitro but its exact structure was unknown. HRAS is a member of the RAS proto‐oncogene family. About 19 % of US cancer patients carry mutations in RAS genes. We solved the structure of iHRAS at 1.77 Å resolution. The structure reveals that iHRAS folds into a double hairpin. The two double hairpins associate in an antiparallel fashion, forming an i‐motif dimer capped by two loops on each end and linked by a connecting region. Six C−C + base pairs form each i‐motif core, and the core regions are extended by a G−G base pair and a cytosine stacking. Extensive canonical and non‐canonical base pairing and stacking stabilizes the connecting region and loops. The iHRAS structure is the first atomic resolution structure of an i‐motif from a human oncogene. This structure sheds light on i‐motifs folding and function in the cell.
Li, Kevin S.↗
A Suppression-based STDP Rule Resilient to Jitter Noise in Spike Patterns for Neuromorphic Computing
Multi-spike models of synaptic plasticity, such as the triplet and suppression spike-timing-dependent plasticity (STDP) rules, exhibit better alignment with neurophysiological data in the brain compared to the pair-based STDP rule. Previous studies have empirically shown that the pair-based STDP rule can detect spatiotemporal spike patterns hidden in equally dense distractor spike trains in an unsupervised manner. However, it fails to detect spike patterns influenced by jitter noise. Given that spiking neural networks (SNNs) exhibit variability in generated spike trains in response to the same inputs, it becomes imperative to have learning rules capable of detecting spike patterns even in the presence of jitter noise. In this study, we introduce a simplified suppression-based STDP rule that demonstrates significantly enhanced tolerance to jitter in spike patterns compared to the pair-based STDP rule. Unlike the ideal suppression STDP rule, characterized by an exponential learning window and requiring high-resolution synapses, the simplified rule limits the synaptic efficacy update to a single bit at any given instant. Moreover, it employs 4-bit fixed-point synapses, facilitating straightforward implementation in neuromorphic hardware.
Gautam, Ashish [ORNL]↗
Accurate and efficient open-source implementation of domain-based local pair natural orbital (DLPNO) coupled-cluster theory using a t1-transformed Hamiltonian
We present an efficient, open-source formulation for coupled-cluster theory through perturbative triples with domain-based local pair natural orbitals [DLPNO-CCSD(T)]. Similar to the implementation of the DLPNO-CCSD(T) method found in the ORCA package, the most expensive integral generation and contraction steps associated with the CCSD(T) method are linear-scaling. In this work, we show that the t1-transformed Hamiltonian allows for a less complex algorithm when evaluating the local CCSD(T) energy without compromising efficiency or accuracy. Our algorithm yields sub-kJ mol−1 deviations for relative energies when compared with canonical CCSD(T), with typical errors being on the order of 0.1 kcal mol−1, using our TightPNO parameters. We extensively tested and optimized our algorithm and parameters for non-covalent interactions, which have been the most difficult interaction to model for orbital (PNO)-based methods historically. To highlight the capabilities of our code, we tested it on large water clusters, as well as insulin (787 atoms).
Chemistry↗
DNA dynamics in aqueous solution: opening the double helix
The opening of a DNA base pair is a simple reaction that is a prerequisite for replication, transcription, and other vital biological functions. Understanding the molecular mechanisms of biological reactions is crucial for predicting and, ultimately, controlling them. Realistic computer simulations of the reactions can provide the needed understanding. To model even the simplest reaction in aqueous solution requires hundreds of hours of supercomputing time. We have used molecular dynamics techniques to simulate fraying of the ends of a six base pair double strand of DNA, [TCGCGA]2, where the four bases of DNA are denoted by T (thymine), C (cytosine), G (guanine), and A (adenine), and to estimate the free energy barrier to this process. The calculations, in which the DNA was surrounded by 2,594 water molecules, required 50 hours of CRAY-2 CPU time for every simulated 100 picoseconds. A free energy barrier to fraying, which is mainly characterized by the movement of adenine away from thymine into aqueous environment, was estimated to be 4 kcal/mol. Another fraying pathway, which leads to stacking between terminal adenine and thymine, was also observed. These detailed pictures of the motions and energetics of DNA base pair opening in water are a first step toward understanding how DNA will interact with any molecule.
NASA Discipline Number 52-20↗
Monte Carlo approach in assessing damage in higher order structures of DNA
We have developed a computer monitor of nuclear DNA in the form of chromatin fibre. The fibres are modeled as a ideal solenoid consisting of twenty helical turns with six nucleosomes per turn. The chromatin model, in combination with are Monte Carlo theory of radiation damage induces by charged particles, based on general features of tack structure and stopping power theory, has been used to evaluate the influence of DNA structure on initial damage. An interesting has emerged from our calculations. Our calculated results predict the existence of strong spatial correlations in damage sites associated with the symmetries in the solenoidal model. We have calculated spectra of short fragments of double stranded DNA produced by multiple double strand breaks induced by both high and low LET radiation. The spectra exhibit peaks at multiples of approximately 85 base pairs (the nucleosome periodicity), and approximately 1000 base pairs (solenoid periodicity). Preliminary experiments to investigate the fragment distributions from irradiated DNA, made by B. Rydberg at Lawrence Berkeley Laboratory, confirm the existence of short DNA fragments and are in substantial agreement with the predictions of our theory.
Non-NASA Center↗
Self-replication of chemical systems based on recognition within a double or a triple helix - A realistic hypothesis
A scenario is proposed for the non-enzymatic self-replication of short RNA molecules. The self-replication of an oligopyrimidine strand is considered and the process of template-directed synthesis based on recognition within a double helix is discussed. Replication mechanisms are suggested for selected oligonucleotides. The mechanisms are based on Watson-Crick base pairing between complementary nucleotides as well as Hoogsteen base pairing between a duplex and the complementary third strand. It is suggested that self-replication based on these mechanisms may be accomplished but may result in a substantial amount of misinformation transfer when mixed oligonucleotides are used.
Kanavarioti, Anastassia↗
Crystal structure of a four-stranded intercalated DNA: d(C4)
The crystal structure of d(C4) solved at 2.3-A resolution reveals a four-stranded molecule composed of two interdigitated or intercalated duplexes. The duplexes are held together by hemiprotonated cytosine-cytosine base pairs and are parallel stranded, but the two duplexes point in opposite directions. The molecule has a slow right-handed twist of 12.4 degrees between covalently linked cytosine base pairs, and the base stacking distance is 3.1 A. This is in general agreement with the NMR studies. A biological role for DNA in this conformation is suggested.
Non-NASA Center↗