Search NASA⌕ Search

SEARCH · Search NASA

Results for “RNA structure modeling”

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.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Structure-function analysis of microRNA 3'-end trimming by Nibbler

Nibbler (Nbr) is a 3'-to-5' exoribonuclease whose catalytic 3'-end trimming activity impacts microRNA (miRNA) and PIWI-interacting RNA (piRNA) biogenesis. Here, we report on structural and functional studies to decipher the contributions of Nbr’s N-terminal domain (NTD) and exonucleolytic domain (EXO) in miRNA 3'-end trimming. We have solved the crystal structures of the NTD core and EXO domains of Nbr, both in the apo-state. The NTD-core domain of Aedes aegypti Nbr adopts a HEAT-like repeat scaffold with basic patches constituting an RNA-binding surface exhibiting a preference for binding double-strand RNA (dsRNA) over single-strand RNA (ssRNA). Structure-guided functional assays in Drosophila S2 cells confirmed a principal role of the NTD in exonucleolytic miRNA trimming, which depends on basic surface patches. Gain-of-function experiments revealed a potential role of the NTD in recruiting Nbr to Argonaute-bound small RNA substrates. The EXO domain of A. aegypti and Drosophila melanogaster Nbr adopt a mixed α/β-scaffold with a deep pocket lined by a DEDDy catalytic cleavage motif. We demonstrate that Nbr’s EXO domain exhibits Mn2+-dependent ssRNA-specific 3'-to-5' exoribonuclease activity. Modeling of a 3' terminal Uridine into the catalytic pocket of Nbr EXO indicates that 2'-O-methylation of the 3'-U would result in a steric clash with a tryptophan side chain, suggesting that 2'-O-methylation protects small RNAs from Nbr-mediated trimming. Overall, our data establish that Nbr requires its NTD as a substrate recruitment platform to execute exonucleolytic miRNA maturation, catalyzed by the ribonuclease EXO domain.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Noonan syndrome patient-specific induced cardiomyocyte model carrying SOS1 gene variant c.1654A>G

Noonan syndrome (NS) is a dominant autosomal genetic disorder, associated with mutations in several genes that exhibit multisystem abnormal development including cardiac defects. NS associated with the Son of Sevenless homolog 1 (SOS1) gene mutation attributes to the development of cardiomyopathy and congenital heart defects. Since the treatment option for NS is very limited, an in vitro disease model with SOS1 gene mutation would be beneficial for exploring therapeutic possibilities for NS. We reprogrammed cardiac fibroblasts obtained from a NS patient and normal control skin fibroblasts (C-SF) into induced pluripotent stem cells (iPSCs). We identified NS-iPSCs carry a heterozygous single nucleotide variation in the SOS1 gene at the c.1654A > G. Furthermore, the control and NS-iPSCs were differentiated into induced cardiomyocytes (iCMCs), and the electron microscopic analysis showed that the sarcomeres of the NS-iCMCs were highly disorganized. FACS analysis showed that 47.5% of the NS-iCMCs co-expressed GATA4 and cardiac troponin T proteins, and the mRNA expression levels of many cardiac related genes, studied by qRT-PCR array, were significantly reduced when compared to the control C-iCMCs. We report for the first time that NS-iPSCs carry a single nucleotide variation in the SOS1 gene at the c.1654A>G were showing significantly reduced cardiac genes and proteins expression as well as structurally and functionally compromised when compared to C-iCMCs. These iPSCs and iCMCs can be used as a modeling platform to unravel the pathologic mechanisms and also the development of novel drug for the cardiomyopathy in patients with NS.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2006

Microgravity and Immunity: Changes in Lymphocyte Gene Expression

Earlier studies had shown that modeled and true microgravity (MG) cause multiple direct effects on human lymphocytes. MG inhibits lymphocyte locomotion, suppresses polyclonal and antigen-specific activation, affects signal transduction mechanisms, as well as activation-induced apoptosis. In this study we assessed changes in gene expression associated with lymphocyte exposure to microgravity in an attempt to identify microgravity-sensitive genes (MGSG) in general and specifically those genes that might be responsible for the functional and structural changes observed earlier. Two sets of experiments targeting different goals were conducted. In the first set, T-lymphocytes from normal donors were activated with antiCD3 and IL2 and then cultured in 1g (static) and modeled MG (MMG) conditions (Rotating Wall Vessel bioreactor) for 24 hours. This setting allowed searching for MGSG by comparison of gene expression patterns in zero and 1 g gravity. In the second set - activated T-cells after culturing for 24 hours in 1g and MMG were exposed three hours before harvesting to a secondary activation stimulus (PHA) thus triggering the apoptotic pathway. Total RNA was extracted using the RNeasy isolation kit (Qiagen, Valencia, CA). Affymetrix Gene Chips (U133A), allowing testing for 18,400 human genes, were used for microarray analysis. In the first set of experiments MMG exposure resulted in altered expression of 89 genes, 10 of them were up-regulated and 79 down-regulated. In the second set, changes in expression were revealed in 85 genes, 20 were up-regulated and 65 were down-regulated. The analysis revealed that significant numbers of MGS genes are associated with signal transduction and apoptotic pathways. Interestingly, the majority of genes that responded by up- or down-regulation in the alternative sets of experiments were not the same, possibly reflecting different functional states of the examined T-lymphocyte populations. The responder genes (MGSG) might play an essential role in adaptation to MG and/or be responsible for pathologic changes encountered in Space and thus represent potential targets for molecular-based countermeasures

Risin, D.↗
DOE OSTI2021

CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. Additionally, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. Finally, the capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

RNA nanotechnology to build a dodecahedral genome of single-stranded RNA virus

The quest for artificial RNA viral complexes with authentic structure while being non-replicative is on its way for the development of viral vaccines. RNA viruses contain capsid proteins that interact with the genome during morphogenesis. The sequence and properties of the protein and genome determine the structure of the virus. For example, the Pariacoto virus ssRNA genome assembles into a dodecahedron. Virus-inspired nanotechnology has progressed remarkably due to the unique structural and functional properties of viruses, which can inspire the design of novel nanomaterials. RNA is a programmable biopolymer able to self-assemble sophisticated 3D structures with rich functionalities. RNA dodecahedrons mimicking the Pariacoto virus quasi-icosahedral genome structures were constructed from both native and 2'-F modified RNA oligos. The RNA dodecahedron easily self-assembled using the stable pRNA three-way junction of bacteriophage phi29 as building blocks. The RNA dodecahedron cage was further characterized by cryo-electron microscopy and atomic force microscopy, confirming the spontaneous and homogenous formation of the RNA cage. The reported RNA dodecahedron cage will likely provide further studies on the mechanisms of interaction of the capsid protein with the viral genome while providing a template for further construction of the viral RNA scaffold to add capsid proteins for the assembly of the viral nucleocapsid as a model. Understanding the self-assembly and RNA folding of this RNA cage may offer new insights into the 3D organization of viral RNA genomes. Finally, the reported RNA cage also has the potential to be explored as a novel virus-inspired nanocarrier.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Deep structural clustering for single-cell RNA-seq data jointly through autoencoder and graph neural network

Abstract Single-cell RNA sequencing (scRNA-seq) permits researchers to study the complex mechanisms of cell heterogeneity and diversity. Unsupervised clustering is of central importance for the analysis of the scRNA-seq data, as it can be used to identify putative cell types. However, due to noise impacts, high dimensionality and pervasive dropout events, clustering analysis of scRNA-seq data remains a computational challenge. Here, we propose a new deep structural clustering method for scRNA-seq data, named scDSC, which integrate the structural information into deep clustering of single cells. The proposed scDSC consists of a Zero-Inflated Negative Binomial (ZINB) model-based autoencoder, a graph neural network (GNN) module and a mutual-supervised module. To learn the data representation from the sparse and zero-inflated scRNA-seq data, we add a ZINB model to the basic autoencoder. The GNN module is introduced to capture the structural information among cells. By joining the ZINB-based autoencoder with the GNN module, the model transfers the data representation learned by autoencoder to the corresponding GNN layer. Furthermore, we adopt a mutual supervised strategy to unify these two different deep neural architectures and to guide the clustering task. Extensive experimental results on six real scRNA-seq datasets demonstrate that scDSC outperforms state-of-the-art methods in terms of clustering accuracy and scalability. Our method scDSC is implemented in Python using the Pytorch machine-learning library, and it is freely available at https://github.com/DHUDBlab/scDSC.

Gan, Yanglan↗
NASA NTRS1993

Montmorillonite catalysis of RNA oligomer formation in aqueous solution. A model for the prebiotic formation of RNA

Oligomers of adenylic acid of up to the 11-mer in length are formed by the reaction of the phosphorimidazolide of adenosine (ImpA) in pH 8 aqueous solution at room temperature in the presence of Na(+)-montmorillonite. These oligomers are joined by phosphodiester bonds in which the 3',5'-linkage predominates over the 2',5'-linkage by a 2:1 ratio. Reaction of a 9:1 mixture of ImpA, A5'ppA results in the formation of oligomers with a 3:1 ratio of 3',5'- to 2',5'-linked phosphodiester bonds. A high proportion of these oligomers contain the A5'ppA grouping. A5'ppA reacts much more rapidly with ImpA than does 5'-ADP (ppA) or 5'-ATP (pppA). The exchangeable cation associated with the montmorillonite effects the observed catalysis with Li+, Na+, NH4+, and Ca2+ being the more effective while Mg2+ and Al3+ are almost ineffective catalysts. 2',5'-Linked oligomers, up to the tetramer in length, are formed using UO2(2+)-montmorillonite. The structure analysis of individual oligomer fractions was performed by selective enzymatic and KOH hydrolytic studies followed by HPLC analysis of the reaction products. It is concluded from the composition of the oligomers that the rate of addition ImpA to a 3'-terminus containing a 2',5'-linkage is slower than the addition to a nucleoside joined by a 3',5'-linked phosphodiester bond. The potential importance of mineral catalysis of the formation of RNA and other oligomers on primitive Earth is discussed.

Non-NASA Center↗
NASA NTRS2000

A DNA enzyme with N-glycosylase activity

In vitro evolution was used to develop a DNA enzyme that catalyzes the site-specific depurination of DNA with a catalytic rate enhancement of about 10(6)-fold. The reaction involves hydrolysis of the N-glycosidic bond of a particular deoxyguanosine residue, leading to DNA strand scission at the apurinic site. The DNA enzyme contains 93 nucleotides and is structurally complex. It has an absolute requirement for a divalent metal cation and exhibits optimal activity at about pH 5. The mechanism of the reaction was confirmed by analysis of the cleavage products by using HPLC and mass spectrometry. The isolation and characterization of an N-glycosylase DNA enzyme demonstrates that single-stranded DNA, like RNA and proteins, can form a complex tertiary structure and catalyze a difficult biochemical transformation. This DNA enzyme provides a new approach for the site-specific cleavage of DNA molecules.

Non-NASA Center↗
DOE OSTI2023

Structural basis of TFIIIC-dependent RNA polymerase III transcription initiation

RNA polymerase III (Pol III) is responsible for transcribing 5S ribosomal RNA (5S rRNA), tRNAs, and other short non-coding RNAs. Its recruitment to the 5S rRNA promoter requires transcription factors TFIIIA, TFIIIC, and TFIIIB. Here, we use cryoelectron microscopy (cryo-EM) to visualize the S. cerevisiae complex of TFIIIA and TFIIIC bound to the promoter. Gene-specific factor TFIIIA interacts with DNA and acts as an adaptor for TFIIIC-promoter interactions. We also visualize DNA binding of TFIIIB subunits, Brf1 and TBP (TATA-box binding protein), which results in the full-length 5S rRNA gene wrapping around the complex. Our smFRET study reveals that the DNA within the complex undergoes both sharp bending and partial dissociation on a slow timescale, consistent with the model predicted from our cryo-EM results. Our findings provide new insights into the transcription initiation complex assembly on the 5S rRNA promoter and allow us to directly compare Pol III and Pol II transcription adaptations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Dynamic basis of supercoiling-dependent DNA interrogation by Cas12a via R-loop intermediates

The sequence specificity and programmability of DNA binding and cleavage have enabled widespread applications of CRISPR-Cas12a in genetic engineering. As an RNA-guided CRISPR endonuclease, Cas12a engages a 20-base pair (bp) DNA segment by forming a three-stranded R-loop structure in which the guide RNA hybridizes to the DNA target. Here we use single-molecule torque spectroscopy to investigate the dynamics and mechanics of R-loop formation of two widely used Cas12a orthologs at base-pair resolution. We directly observe kinetic intermediates corresponding to a ~5 bp initial RNA-DNA hybridization and a ~17 bp intermediate preceding R-loop completion, followed by transient DNA unwinding that extends beyond the 20 bp R-loop. The complex multistate landscape of R-loop formation is ortholog-dependent and shaped by target sequence, mismatches, and DNA supercoiling. A four-state kinetic model captures essential features of Cas12a R-loop dynamics and provides a biophysical framework for understanding Cas12a activity and specificity.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The role of SAXS and molecular simulations in 3D structure elucidation of a DNA aptamer against lung cancer

Aptamers are short, single-stranded DNA or RNA oligonucleotide molecules that function as synthetic analogs of antibodies and bind to a target molecule with high specificity. Aptamer affinity entirely depends on its tertiary structure and charge distribution. Therefore, length and structure optimization are essential for increasing aptamer specificity and affinity. Here, we present a general optimization procedure for finding the most populated atomistic structures of DNA aptamers. Based on the existed aptamer LC-18 for lung adenocarcinoma, a new truncated LC-18 (LC-18t) aptamer LC-18t was developed. A three-dimensional (3D) shape of LC-18t was reported based on small-angle X-ray scattering (SAXS) experiments and molecular modeling by fragment molecular orbital or molecular dynamic methods. Molecular simulations revealed an ensemble of possible aptamer conformations in solution that were in close agreement with measured SAXS data. The aptamer LC-18t had stronger binding to cancerous cells in lung tumor tissues and shared the binding site with the original larger aptamer. The suggested approach reveals 3D shapes of aptamers and helps in designing better affinity probes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Structural basis for a highly conserved RNA-mediated enteroviral genome replication

Abstract Enteroviruses contain conserved RNA structures at the extreme 5′ end of their genomes that recruit essential proteins 3CD and PCBP2 to promote genome replication. However, the high-resolution structures and mechanisms of these replication-linked RNAs (REPLRs) are limited. Here, we determined the crystal structures of the coxsackievirus B3 and rhinoviruses B14 and C15 REPLRs at 1.54, 2.2 and 2.54 Å resolution, revealing a highly conserved H-type four-way junction fold with co-axially stacked sA-sD and sB-sC helices that are stabilized by a long-range A•C•U base-triple. Such conserved features observed in the crystal structures also allowed us to predict the models of several other enteroviral REPLRs using homology modeling, which generated models almost identical to the experimentally determined structures. Moreover, our structure-guided binding studies with recombinantly purified full-length human PCBP2 showed that two previously proposed binding sites, the sB-loop and 3′ spacer, reside proximally and bind a single PCBP2. Additionally, the DNA oligos complementary to the 3′ spacer, the high-affinity PCBP2 binding site, abrogated its interactions with enteroviral REPLRs, suggesting the critical roles of this single-stranded region in recruiting PCBP2 for enteroviral genome replication and illuminating the promising prospects of developing therapeutics against enteroviral infections targeting this replication platform.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Viroid-like “obelisk” agents are widespread in the ocean and exceed the abundance of RNA viruses in the prokaryotic fraction

Abstract “Obelisks” are recently discovered ribonucleic acid (RNA) viroid-like elements present in diverse environments with no phylogenetic similarity to any known biological agent. obelisks were first identified in the human gut and in a commensal bacterium acting as a replicative host. They have a circular ∼1 kb RNA genome, rod-like secondary structures, and the encoding of a protein superfamily called “Oblins”. We performed a large-scale search of obelisks in the ocean using the Pebblescout program and the transcriptomic Sequence Archive Read databases, revealing the biogeography and abundance of these viroid-like RNA elements. We detected 55 obelisk genomes resulting in 35 marine clusters at the species level. These obelisks were detected in the prokaryotic fraction and to a lesser extent in the eukaryotic fraction, and distributed across all the oceans from surface to mesopelagic including the Arctic, and even in the coldest seawater of Earth beneath the Antarctic Ross Ice Shelf. The obelisk hallmark protein Oblin-1 confirmed by 3D models was found in various marine samples. Some of the detected marine obelisks harbor hammerhead self-cleaving ribozymes in both polarities. In the prokaryotic, but not the eukaryotic, fraction of the Tara Ocean dataset, relative abundance of obelisks calculated by transcriptomic fragment recruitment indicated that they are abundant in marine samples, reaching or even exceeding the relative abundance of the previously discovered uncultured RNA viruses. In conclusion, obelisks are abundant and widespread viroid-like elements that should be included in ocean biogeochemical models.

Environmental Sciences & Ecology↗
DOE OSTI2024

Did the exposure of coacervate droplets to rain make them the first stable protocells?

Membraneless coacervate microdroplets have long been proposed as model protocells as they can grow, divide, and concentrate RNA by natural partitioning. However, the rapid exchange of RNA between these compartments, along with their rapid fusion, both within minutes, means that individual droplets would be unable to maintain their separate genetic identities. Hence, Darwinian evolution would not be possible, and the population would be vulnerable to collapse due to the rapid spread of parasitic RNAs. In this study, we show that distilled water, mimicking rain/freshwater, leads to the formation of electrostatic crosslinks on the interface of coacervate droplets that not only suppress droplet fusion indefinitely but also allow the spatiotemporal compartmentalization of RNA on a timescale of days depending on the length and structure of RNA. We suggest that these nonfusing membraneless droplets could potentially act as protocells with the capacity to evolve compartmentalized ribozymes in prebiotic environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1997

Equally parsimonious pathways through an RNA sequence space are not equally likely

An experimental system for determining the potential ability of sequences resembling 5S ribosomal RNA (rRNA) to perform as functional 5S rRNAs in vivo in the Escherichia coli cellular environment was devised previously. Presumably, the only 5S rRNA sequences that would have been fixed by ancestral populations are ones that were functionally valid, and hence the actual historical paths taken through RNA sequence space during 5S rRNA evolution would have most likely utilized valid sequences. Herein, we examine the potential validity of all sequence intermediates along alternative equally parsimonious trajectories through RNA sequence space which connect two pairs of sequences that had previously been shown to behave as valid 5S rRNAs in E. coli. The first trajectory requires a total of four changes. The 14 sequence intermediates provide 24 apparently equally parsimonious paths by which the transition could occur. The second trajectory involves three changes, six intermediate sequences, and six potentially equally parsimonious paths. In total, only eight of the 20 sequence intermediates were found to be clearly invalid. As a consequence of the position of these invalid intermediates in the sequence space, seven of the 30 possible paths consisted of exclusively valid sequences. In several cases, the apparent validity/invalidity of the intermediate sequences could not be anticipated on the basis of current knowledge of the 5S rRNA structure. This suggests that the interdependencies in RNA sequence space may be more complex than currently appreciated. If ancestral sequences predicted by parsimony are to be regarded as actual historical sequences, then the present results would suggest that they should also satisfy a validity requirement and that, in at least limited cases, this conjecture can be tested experimentally.

NASA Discipline Exobiology↗
DOE OSTI2022

Chromatin structure predicts survival in glioma patients

The pathological changes in epigenetics and gene regulation that accompany the progression of low-grade to high-grade gliomas are under-studied. The authors use a large set of paired atac-seq and RNA-seq data from surgically resected glioma specimens to infer gene regulatory relationships in glioma. Thirty-eight glioma patient samples underwent atac-seq sequencing and 16 samples underwent additional RNA-seq analysis. Using an atac-seq/RNA-seq correlation matrix, atac-seq peaks were paired with genes based on high correlation values (|r 2 | > 0.6). Samples clustered by IDH1 status but not by grade. Surprisingly there was a trend for IDH1 mutant samples to have more peaks. The majority of peaks are positively correlated with survival and positively correlated with gene expression. Constructing a model of the top six atac-seq peaks created a highly accurate survival prediction model (r 2 = 0.68). Four of these peaks were still significant after controlling for age, grade, pathology, IDH1 status and gender. Grade II, III, and IV (primary) samples have similar transcription factors and gene modules. However, grade IV (recurrent) samples have strikingly few peaks. Patient-derived glioma cultures showed decreased peak counts following radiation indicating that this may be radiation-induced. This study supports the notion that IDH1 mutant and IDH1 wildtype gliomas have different epigenetic landscapes and that accessible chromatin sites mapped by atac-seq peaks tend to be positively correlated with expression. The data in this study leads to a new model of treatment response wherein glioma cells respond to radiation therapy by closing open regions of DNA.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1998

Changing patterns of localization of the tobacco mosaic virus movement protein and replicase to the endoplasmic reticulum and microtubules during infection

Tobacco mosaic virus (TMV) derivatives that encode movement protein (MP) as a fusion to the green fluorescent protein (MP:GFP) were used in combination with antibody staining to identify host cell components to which MP and replicase accumulate in cells of infected Nicotiana benthamiana leaves and in infected BY-2 protoplasts. MP:GFP and replicase colocalized to the endoplasmic reticulum (ER; especially the cortical ER) and were present in large, irregularly shaped, ER-derived structures that may represent "viral factories." The ER-derived structures required an intact cytoskeleton, and microtubules appeared to redistribute MP:GFP from these sites during late stages of infection. In leaves, MP:GFP accumulated in plasmodesmata, whereas in protoplasts, the MP:GFP was targeted to distinct, punctate sites near the plasma membrane. Treating protoplasts with cytochalasin D and brefeldin A at the time of inoculation prevented the accumulation of MP:GFP at these sites. It is proposed that the punctate sites anchor the cortical ER to plasma membrane and are related to sites at which plasmodesmata form in walled cells. Hairlike structures containing MP:GFP appeared on the surface of some of the infected protoplasts and are reminiscent of similar structures induced by other plant viruses. We present a model that postulates the role of the ER and cytoskeleton in targeting the MP and viral ribonucleoprotein from sites of virus synthesis to the plasmodesmata through which infection is spread.

NASA Discipline Plant Biology↗
DOE OSTI2021

Peptide nucleic acid Hoogsteen strand linker design for major groove recognition of DNA thymine bases

Sequence-specific targeting of double-stranded DNA and non-coding RNA via triple-helix-forming peptide nucleic acids (PNAs) has attracted considerable attention in therapeutic, diagnostic and nanotechnological fields. An E-base (3-oxo-2,3-dihydropyridazine), attached to the polyamide backbone of a PNA Hoogsteen strand by a side-chain linker molecule, is typically used in the hydrogen bond recognition of the 4-oxo group of thymine and uracil nucleic acid bases in the major groove. We report on the application of quantum chemical computational methods, in conjunction with spatial constraints derived from the experimental structure of a homopyrimidine PNA·DNA-PNA hetero-triplex, to investigate the influence of linker flexibility on binding interactions of the E-base with thymine and uracil bases in geometry-optimised model systems. Hydrogen bond formation between the N2 E-base atom and target pyrimidine base 4-oxo groups in model systems containing a β-alanine linker (J Am Chem Soc 119:11116, 1997) was found to incur significant internal strain energy and the potential disruption of intra-stand aromatic base stacking interactions in an oligomeric context. Additionally, in geometry-optimised model systems containing a 3-trans olefin linker (Bioorg Med Chem Lett 14:1551, 2004) the E-base swung out away from the target pyrimidine bases into the solvent. These findings are in qualitative agreement with calorimetric measurements in hybridisation experiments at T–A and U–A inversion sites. In contrast, calculations on a novel 2-cis olefin linker design indicate that it could permit simultaneous E-base hydrogen bonding with the thymine 4-oxo group, circumvention and solvent screening of the thymine 5-methyl group, and maintenance of triplex intra-stand base stacking interactions.

59 BASIC BIOLOGICAL SCIENCES↗