Unveiling the solution structure of a DNA duplex with continuous silver-modified Watson-Crick base pairs
Not provided.
SEARCH · Search NASA
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.
Not provided.
RNA self-copying using an activated adenosine substrate is challenging. Our structure studies reveal the binding motifs of adenosine substrates. The obtained structures explain the disfavorable conformations for nucleophilic reactions.
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.
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.
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.
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.
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.
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.
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.
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.
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).
Nutrient limitation has been shown to reduce bacterial genome size and influence nucleotide composition; however, much of this work has been conducted in marine systems and the factors which shape soil bacterial genomic traits remain largely unknown. Here, for this work, we determined average genome size, GC content, codon usage, and amino acid content from 398 soil metagenomes across a broad geographic range and used machine-learning to determine the environmental parameters that most strongly explain the distribution of these traits. We found that genomic trait averages were most related to pH, which we suggest is primarily due to the correlation of pH with several environmental parameters, particularly soil carbon content. Low pH soils had higher carbon to nitrogen ratios (C:N) and tended to have communities with lower GC content and larger genomes, potentially a response to increased physiological stress and a requirement for metabolic diversity. Conversely, communities in high pH and low soil C:N had smaller genomes and higher GC content—indicating potential resource driven selection against AT base pairs, which have a higher C:N than GC base pairs. Similarly, we found that nutrient conservation also applied to amino acid stoichiometry, where bacteria in soils with low C:N ratios tended to code for amino acids with lower C:N. Together, these relationships point towards fundamental mechanisms that underpin genome size, and nucleotide and amino acid selection in soil bacteria.
The interaction energy of base–acid plays a key role in acid retention of phosphoric acid (PA)-doped polymer electrolytes under fuel cell operating conditions. Here, we investigate the energetics of proton-accepting and hydroxide-donating organic bases using density functional theory calculations. Because of their weak basicity, proton-accepting organic bases such as benzimidazole have relatively low interaction energy with the acid in the absence of water (15.3–28.0 kcal mol –1 ). Energetics of the proton-accepting base–PA complex increases by adding water, indicating that the interactions in the base–acid complex strengthen in the presence of water. On the other hand, hydroxide-donating organic bases, such as tetramethylammonium hydroxide, have high interaction energy with PA (~110 kcal mol –1 ), which remains high in the presence of water. The chemical shifts of 31P NMR support the energetics of the base–acid complexes. Overall, this study further discusses the benefit of incorporating hydroxide-donating organic bases into the polymeric structure over proton-accepting bases as a way to increase acid retention.
PrimPol is a human DNA polymerase-primase that localizes to mitochondria and nucleus and bypasses the major oxidative lesion 7,8-dihydro-8-oxoguanine (oxoG) via translesion synthesis, in mostly error-free manner. We present structures of PrimPol insertion complexes with a DNA template-primer and correct dCTP or erroneous dATP opposite the lesion, as well as extension complexes with C or A as a 3'–terminal primer base. We show that during the insertion of C and extension from it, the active site is unperturbed, reflecting the readiness of PrimPol to accommodate oxoG(anti). The misinsertion of A opposite oxoG(syn) also does not alter the active site, and is likely less favorable due to lower thermodynamic stability of the oxoG(syn)•A base-pair. During the extension step, oxoG(syn) induces an opening of its base-pair with A or misalignment of the 3'-A primer terminus. Together, the structures show how PrimPol accurately synthesizes DNA opposite oxidatively damaged DNA in human cells.
The CCCTC-binding factor (CTCF) binds tens of thousands of enhancers and promoters on mammalian chromosomes by means of its 11 tandem zinc finger (ZF) DNA-binding domain. In addition to the 12–15-bp CORE sequence, some of the CTCF binding sites contain 5' upstream and/or 3' downstream motifs. Here, we describe two structures for overlapping portions of human CTCF, respectively, including ZF1–ZF7 and ZF3–ZF11 in complex with DNA that incorporates the CORE sequence together with either 3' downstream or 5' upstream motifs. Like conventional tandem ZF array proteins, ZF1–ZF7 follow the right-handed twist of the DNA, with each finger occupying and recognizing one triplet of three base pairs in the DNA major groove. ZF8 plays a unique role, acting as a spacer across the DNA minor groove and positioning ZF9–ZF11 to make cross-strand contacts with DNA. We ascribe the difference between the two subgroups of ZF1–ZF7 and ZF8–ZF11 to residues at the two positions -6 and -5 within each finger, with small residues for ZF1–ZF7 and bulkier and polar/charged residues for ZF8–ZF11. ZF8 is also uniquely rich in basic amino acids, which allows salt bridges to DNA phosphates in the minor groove. Highly specific arginine–guanine and glutamine–adenine interactions, used to recognize G:C or A:T base pairs at conventional base-interacting positions of ZFs, also apply to the cross-strand interactions adopted by ZF9–ZF11. The differences between ZF1–ZF7 and ZF8–ZF11 can be rationalized structurally and may contribute to recognition of high-affinity CTCF binding sites.
Missing linkers induce coordinatively unsaturated Lewis acid-base pair sites in metal–organic frameworks such as UiO-66 (Zr). The Lewis acid-base pair is induced by replacing the benzene dicarboxylic acid linker with formic acid, which subsequently is thermally decomposed. The concentration of coordinatively unsaturated Zr 4+ cations was titrated with acetone. In parallel to their concentration, the catalytic activity increased while the selectivity for aldol condensation between furfural and acetone remained constant. The missing linkers generate a complex ensemble consisting of an (at least partly) accessible Zr 4+ cation surrounded by oxygen, catalyzing the cross-aldol condensation. As a result, density functional theory calculations, matching the experimental findings provide a more detailed understanding of the elementary steps during the condensation reaction.
Dynamic programming algorithms within the NUPACK software suite enable analysis of equilibrium base-pairing properties for complex and test tube ensembles containing arbitrary numbers of interacting nucleic acid strands. Currently, calculations are limited to single-material systems that are either all-RNA or all-DNA. Here, to enable analysis of mixed-material systems that are critical for modern applications in vitro, in situ, and in vivo, we develop physical models and dynamic programming algorithms that allow the material of the system to be specified at nucleotide resolution. Free energy parameter sets are constructed for both RNA/DNA and RNA/2'OMe-RNA mixed-material systems by combining available empirical mixed-material parameters with single-material parameter sets to enable treatment of the full complex and test tube ensembles. New dynamic programming recursions account for the material of each nucleotide throughout the recursive process. For a complex with N nucleotides, the mixed-material dynamic programming algorithms maintain the O(N 3 ) time complexity of the single-material algorithms, enabling efficient calculation of diverse physical quantities over complex and test tube ensembles (e.g., complex partition function, equilibrium complex concentrations, equilibrium base-pairing probabilities, minimum free energy secondary structure(s), and Boltzmann-sampled secondary structures) at a cost increase of roughly 2.0-3.5×. The results of existing single-material algorithms are exactly reproduced when applying the new mixed-material algorithms to single-material systems. Accuracy is significantly enhanced using mixed-material models and algorithms to predict RNA/DNA and RNA/2'OMe-RNA duplex melting temperatures from the experimental literature as well as RNA/DNA melt profiles from new experiments. In conclusion, mixed-material analyses can be performed online using the NUPACK web app (www.nupack.org) or locally using the NUPACK Python module.
Understanding DNA-ligand binding interactions requires ligand screening, crystallization, and structure determination. In order to obtain insights into the amyloid peptide precursor (APP) gene-Thioflavin T (ThT) interaction, single crystals of two DNA sequences 5'-GCCCACCACGGC-3' (PDB 8ASK) and d(CCGGGGTACCCCGG) 2 (PDB 8ASH) were grown in the presence of ThT or its analogue 2-((4-(dimethylamino)benzylidene)amino)-3,6-dimethylbenzo[d]thiazol-3-ium iodide (XRB). Both structures were solved by molecular replacement. In the case of 8ASK, the space group was H3 with unit cell dimensions of a = b = 64.49 Å, c = 46.19 Å. Phases were obtained using a model generated by X3DNA. The novel 12-base-pair B-DNA structure did not have extra density for the ThT ligand. The 14-base-pair A-DNA structure with bound ThT analog XRB was isomorphous with previously the obtained apo-DNA structure 5WV7 (space group was P 4 1 2 1 2 with unit cell dimensions a = b = 41.76 Å, c = 88.96 Å). Binding of XRB to DNA slightly changes the DNA's buckle parameters at the CpG regions. Comparison of the two conformations of the XRB molecule: alone and bound to DNA indicates that the binding results from the freedom of rotation of the two aromatic rings.