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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Identifying candidate structured RNAs in CRISPR operons

Noncoding RNAs with secondary structures play important roles in CRISPR-Cas systems. Many of these structures likely remain undiscovered. We used a large-scale comparative genomics approach to predict 156 novel candidate structured RNAs from 36,111 CRISPR-Cas systems. A number of these were found to overlap with coding genes, including palindromic candidates that overlapped with a variety of Cas genes in type I and III systems. Among these 156 candidates, we identified 46 new models of CRISPR direct repeats and 1 tracrRNA. This tracrRNA model occasionally overlapped with predicted cas9 coding regions, emphasizing the importance of expanding our search windows for novel structure RNAs in coding regions. We also demonstrated that the antirepeat sequence in this tracrRNA model can be used to accurately assign thousands of predicted CRISPR arrays to type II-C systems. This study highlights the importance of unbiased identification of candidate structured RNAs across CRISPR-Cas systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Improved Biofuel Production through Discovery and Engineering of Terpene Metabolism in Switchgrass

Project Objectives - Of the myriad specialized metabolites that plants deploy to adapt to environmental challenges, terpenes form the largest group. In many major crops, unique terpene blends serve as key stress defenses that directly impact plant fitness and yield. In addition, terpenes, such as bisabolene and pinene, are used for producing renewable biofuels. Essential to advancing a broader use of terpenes for biofuel feedstock engineering is a system-wide knowledge of the diverse biosynthetic machinery and defensive potential of often species-specific terpene blends. The proposed project would merge genome-wide enzyme discovery with comparative –omics, protein structural and plant microbiome studies to define the biosynthesis and stress-defensive functions of the switchgrass (Panicum virgatum) terpene network. These insights would be combined with developing and applying non-transgenic genome editing tools to design plants with desirable terpene blends for higher productivity and biofuel production on marginal lands. As a dedicated lignocellulosic feedstock for U.S. biofuel production with high net energy yield, stress tolerance, and available genome resources, switchgrass is well-suited for devising new avenues for biofuel production. Project Description – The diversity of plant terpene defenses is governed by species-specific families of terpene synthase (TPS) and cytochrome P450 monooxygenase (P450) enzymes. Mining of the switchgrass genome (genotype Alamo) identified ~100 TPS and P450 candidate genes, and combinatorial biochemical analysis of synthesized TPSs and P450s revealed more than a dozen enzymes with common and novel activities. In addition, several identified terpene metabolites and the corresponding transcripts were up-regulated in response to abiotic stressors. These findings demonstrate a unique switchgrass terpene network with probable importance to abiotic stress tolerance, thus providing a large chemical portfolio for optimizing crop resistance, yield, and biofuel composition. Leveraging these preliminary data, we propose to generate a genome-wide map of the switchgrass terpene metabolic network through multi-gene co-expression analyses that allow the efficient cross-validation of TPS and P450 functions. Key enzymes would further be applied to structure-function studies via X-ray protein crystallography, homology modeling and site-directed mutagenesis to gain mechanistic insight into the catalytic specificity of switchgrass terpene metabolism and provide gene and amino acid targets for genome editing. In tandem with terpene pathway discovery, system-wide metabolomics, transcriptomics and proteomics studies in switchgrass accessions of contrasting drought tolerance would define the role of switchgrass terpene metabolism in conferring abiotic stress resilience. Metabolic changes would be assessed in a combined approach of targeted (terpenes) and untargeted metabolite profiling using a high-resolution LC-MS/MS approach, differential gene expression analyses through multiplexed Illumina RNA sequencing, and quantitative analysis of high-priority pathway enzymes using multiple reaction monitoring (MRM). Drawing on these insights, knock-down/out mutants of stress-associated pathway nodes would be generated by optimizing transient virus-induced gene silencing (VIGS) and CRISPR/Cas9 systems under control of the Tobacco Rattle Virus (TRV). The resulting mutant lines would then be analyzed for stress susceptibility and the impact on the root microbiome to define gene functions in planta. Knowledge of terpene pathways, enzyme mechanisms and bioactivities would be applied to enhance switchgrass stress resilience and to tailor-make terpene blends for biofuel production. Here, TRV-enabled CRISPR/Cas9 genome editing, including allele-specific knock-out of redundant genes, engineering of enzyme specificity via structure-guided point mutations, and overexpression of terpene genes relevant to stress-protection or biofuel production, would be used to increase metabolic flux toward desired pathways. Broader Impacts - Integrating the system-wide discovery, mechanistic analysis and non-transgenic genome engineering of the switchgrass terpene network aligns the required steps to unlock the chemical potential of this important metabolite class to generate crops that are more resistant to stress and provide advanced biofuel production in light of rising climate pressures as foreseeable challenges for bioenergy crop cultivation. The proposed project would further offer interdisciplinary student training through active involvement in the project and integration of research concepts and outcomes into newly-developed graduate and undergraduate courses on Plant Biotechnology.

09 BIOMASS FUELS↗
DOE OSTI2024

Biallelic variants in SNUPN cause a limb girdle muscular dystrophy with myofibrillar-like features

Abstract Alterations in RNA-splicing are a molecular hallmark of several neurological diseases, including muscular dystrophies, where mutations in genes involved in RNA metabolism or characterized by alterations in RNA splicing have been described. Here, we present five patients from two unrelated families with a limb-girdle muscular dystrophy (LGMD) phenotype carrying a biallelic variant in SNUPN gene. Snurportin-1, the protein encoded by SNUPN, plays an important role in the nuclear transport of small nuclear ribonucleoproteins (snRNPs), essential components of the spliceosome. We combine deep phenotyping, including clinical features, histopathology and muscle MRI, with functional studies in patient-derived cells and muscle biopsies to demonstrate that variants in SNUPN are the cause of a new type of LGMD according to current definition. Moreover, an in vivo model in Drosophila melanogaster further supports the relevance of Snurportin-1 in muscle. SNUPN patients show a similar phenotype characterized by proximal weakness starting in childhood, restrictive respiratory dysfunction and prominent contractures, although inter-individual variability in terms of severity even in individuals from the same family was found. Muscle biopsy showed myofibrillar-like features consisting of myotilin deposits and Z-disc disorganization. MRI showed predominant impairment of paravertebral, vasti, sartorius, gracilis, peroneal and medial gastrocnemius muscles. Conservation and structural analyses of Snurportin-1 p.Ile309Ser variant suggest an effect in nuclear-cytosol snRNP trafficking. In patient-derived fibroblasts and muscle, cytoplasmic accumulation of snRNP components is observed, while total expression of Snurportin-1 and snRNPs remains unchanged, which demonstrates a functional impact of SNUPN variant in snRNP metabolism. Furthermore, RNA-splicing analysis in patients’ muscle showed widespread splicing deregulation, in particular in genes relevant for muscle development and splicing factors that participate in the early steps of spliceosome assembly. In conclusion, we report that SNUPN variants are a new cause of limb girdle muscular dystrophy with specific clinical, histopathological and imaging features, supporting SNUPN as a new gene to be included in genetic testing of myopathies. These results further support the relevance of splicing-related proteins in muscle disorders.

Iruzubieta, Pablo (ORCID:0000000303316222)↗
DOE OSTI2020

Experimental and theoretical rationalization for the base pairing abilities of inosine, guanosine, adenosine, and their corresponding 8-oxo -7,8-dihydropurine, and 8-bromopurine analogues within A-form duplexes of RNA

Inosine is an important RNA modification, furthermore RNA oxidation has gained interest due, in part, to its potential role in the development/progression of disease as well as on its impact on RNA structure and function. In this report we established the base pairing abilities of purine nucleobases G, I, A, as well as their corresponding, 8-oxo-7,8-dihydropurine (common products of oxidation at the C8-position of purines), and 8-bromopurine (as probes to explore conformational changes), derivatives, namely 8-oxoG, 8-oxoI, 8-oxoA, 8-BrG, and 8-BrI. Dodecamers of RNA were obtained using standard phosphoramidite chemistry via solid-phase synthesis, and used as models to establish the impact that each of these nucleobases have on the thermal stability of duplexes, when base pairing to canonical and noncanonical nucleobases. Thermal stabilities were obtained from thermal denaturation transition (Tm) measurements, via circular dichroism (CD). The results were then rationalized using models of base pairs between two monomers, via density functional theory (DFT), that allowed us to better understand potential contributions from H-bonding patterns arising from distinct conformations. Overall, some of the important results indicate that: (a) an anti-I:syn-A base pair provides thermal stability, due to the absence of the exocyclic amine; (b) 8-oxoG base pairs like U, and does not induce destabilization within the duplex when compared to the pyrimidine ring; (c) a U:G wobble-pair is only stabilized by G; and (d) 8-oxoA displays an inherited base pairing promiscuity in this sequence context. Gaining a better understanding of how this oxidatively generated lesions potentially base pair with other nucleobases will be useful to predict various biological outcomes, as well as in the design of biomaterials and/or nucleotide derivatives with biological potential.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Structural basis of RNA cap modification by SARS-CoV-2

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), the causative agent of COVID-19 illness, has caused millions of infections worldwide. In SARS coronaviruses, the non-structural protein 16 (nsp16), in conjunction with nsp10, methylates the 5'-end of virally encoded mRNAs to mimic cellular mRNAs, thus protecting the virus from host innate immune restriction. We report here the high-resolution structure of a ternary complex of SARS-CoV-2 nsp16 and nsp10 in the presence of cognate RNA substrate analogue and methyl donor, S-adenosyl methionine (SAM). The nsp16/nsp10 heterodimer is captured in the act of 2'-O methylation of the ribose sugar of the first nucleotide of SARS-CoV-2 mRNA. We observe large conformational changes associated with substrate binding as the enzyme transitions from a binary to a ternary state. This induced fit model provides mechanistic insights into the 2'-O methylation of the viral mRNA cap. We also discover a distant (25 Å) ligand-binding site unique to SARS-CoV-2, which can alternatively be targeted, in addition to RNA cap and SAM pockets, for antiviral development.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Functional anatomy of zinc finger antiviral protein complexes

Abstract ZAP is an antiviral protein that binds to and depletes viral RNA, which is often distinguished from vertebrate host RNA by its elevated CpG content. Two ZAP cofactors, TRIM25 and KHNYN, have activities that are poorly understood. Here, we show that functional interactions between ZAP, TRIM25 and KHNYN involve multiple domains of each protein, and that the ability of TRIM25 to multimerize via its RING domain augments ZAP activity and specificity. We show that KHNYN is an active nuclease that acts in a partly redundant manner with its homolog N4BP1. The ZAP N-terminal RNA binding domain constitutes a minimal core that is essential for antiviral complex activity, and we present a crystal structure of this domain that reveals contacts with the functionally required KHNYN C-terminal domain. These contacts are remote from the ZAP CpG binding site and would not interfere with RNA binding. Based on our dissection of ZAP, TRIM25 and KHNYN functional anatomy, we could design artificial chimeric antiviral proteins that reconstitute the antiviral function of the intact authentic proteins, but in the absence of protein domains that are otherwise required for activity. Together, these results suggest a model for the RNA recognition and action of ZAP-containing antiviral protein complexes.

Science & Technology - Other Topics↗
DOE OSTI2023

Downregulation of barley ferulate 5-hydroxylase dramatically alters straw lignin structure without impact on mechanical properties

Barley is a major cereal crop for temperate climates, and a diploid genetic model for polyploid wheat. Cereal straw biomass is an attractive source of feedstock for green technologies but lignin, a key determinant of feedstock recalcitrance, complicates bio-conversion processes. However, manipulating lignin content to improve the conversion process could negatively affect agronomic traits. An alternative approach is to manipulate lignin composition which influences the physical and chemical properties of straw. This study validates the function of a barley ferulate 5-hydroxylase gene and demonstrates that its downregulation using the RNA-interference approach substantially impacts lignin composition. We identified five barley genes having putative ferulate 5-hydroxylase activity. Downregulation of HvF5H1 substantially reduced the lignin syringyl/guaiacyl (S/G) ratio in straw while the lignin content, straw mechanical properties, plant growth habit, and grain characteristics all remained unaffected. Metabolic profiling revealed significant changes in the abundance of 173 features in the HvF5H1-RNAi lines. The drastic changes in the lignin polymer of transgenic lines highlight the plasticity of barley lignification processes and the associated potential for manipulating and improving lignocellulosic biomass as a feedstock for green technologies. On the other hand, our results highlight some differences between the lignin biosynthetic pathway in barley, a temperate climate grass, and the warm climate grass, rice, and underscore potential diversity in the lignin biosynthetic pathways in grasses.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Cryo-EM structure of the diapause chaperone artemin

The protein artemin acts as both an RNA and protein chaperone and constitutes over 10% of all protein in Artemia cysts during diapause. However, its mechanistic details remain elusive since no high-resolution structure of artemin exists. Here we report the full-length structure of artemin at 2.04 Å resolution. The cryo-EM map contains density for an intramolecular disulfide bond between Cys22-Cys61 and resolves the entire C-terminus extending into the core of the assembled protein cage but in a different configuration than previously hypothesized with molecular modeling. We also provide data supporting the role of C-terminal helix F towards stabilizing the dimer form that is believed to be important for its chaperoning activity. We were able to destabilize this effect by placing a tag at the C-terminus to fully pack the internal cavity and cause limited steric hindrance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Two related families of metal transferases, ZNG1 and ZNG2, are involved in acclimation to poor Zn nutrition in Arabidopsis

Metal homeostasis has evolved to tightly modulate the availability of metals within the cell, avoiding cytotoxic interactions due to excess and protein inactivity due to deficiency. Even in the presence of homeostatic processes, however, low bioavailability of these essential metal nutrients in soils can negatively impact crop health and yield. While research has largely focused on how plants assimilate metals, acclimation to metal-limited environments requires a suite of strategies that are not necessarily involved in metal transport across membranes. The identification of these mechanisms provides a new opportunity to improve metal-use efficiency and develop plant foodstuffs with increased concentrations of bioavailable metal nutrients. Here, we investigate the function of two distinct subfamilies of the nucleotide-dependent metallochaperones (NMCs), named ZNG1 and ZNG2, that are found in plants, using Arabidopsis thaliana as a reference organism. AtZNG1 (AT1G26520) is an ortholog of human and fungal ZNG1, and like its previously characterized eukaryotic relatives, localizes to the cytosol and physically interacts with methionine aminopeptidase type I (AtMAP1A). Analysis of At ZNG1 , At MAP1A , At MAP2A , and At MAP2B transgenic mutants are consistent with the role of Arabidopsis ZNG1 as a Zn transferase for AtMAP1A, as previously described in yeast and zebrafish. Structural modeling reveals a flexible cysteine-rich loop that we hypothesize enables direct transfer of Zn from AtZNG1 to AtMAP1A during GTP hydrolysis. Based on proteomics and transcriptomics, loss of this ancient and conserved mechanism has pleiotropic consequences impacting the expression of hundreds of genes, including those involved in photosynthesis and vesicle transport. Members of the plant-specific family of NMCs, ZNG2A1 (AT1G80480) and ZNG2A2 (AT1G15730), are also required during Zn deficiency, but their target protein(s) remain to be discovered. In conclusion, RNA-seq analyses reveal wide-ranging impacts across the cell when the genes encoding these plastid-localized NMCs are disrupted.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Spontaneous Reorganization of DNA-Based Polymers in Higher Ordered Structures Fueled by RNA

We demonstrate a strategy that allows for the spontaneous reconfiguration of self-assembled DNA polymers exploiting RNA as chemical fuel. To do this, we have rationally designed orthogonally addressable DNA building blocks that can be transiently deactivated by RNA fuels and subtracted temporarily from participation in the self-assembly process. Through a fine modulation of the rate at which the building blocks are reactivated we can carefully control the final composition of the polymer and convert a disordered polymer in a higher order polymer, which is disfavored from a thermodynamic point of view. We measure the dynamic reconfiguration via fluorescent signals and confocal microscopy, and we derive a kinetic model that captures the experimental results. Our approach suggests a novel route toward the development of biomolecular materials in which engineered chemical reactions support the autonomous spatial reorganization of multiple components.

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

Targeted inhibition of MASTL kinase activity induces apoptosis in breast cancer

Microtubule-associated serine/threonine kinase-like (MASTL) (or Greatwall kinase (GWL)) is an important cell cycle regulating kinase that regulates the G2-M transition. Uncontrolled MASTL activity is implicated in breast cancer progression. To date, very few inhibitors have been reported against this protein. Here, structure-based computational modeling indicates that the natural product flavopiridol (FLV) binds strongly to MASTL and these results are validated using molecular dynamics simulation studies. Further, an in vitro kinase assay reveals an EC 50 (effective concentration) value of FLV to be 82.1 nM and a better IC 50 compared to the positive reference compound, staurosporine. FLV is found to inhibit MASTL kinase activity, arresting the cell growth in the G1 phase and inducing apoptosis in breast cancer cells. Consistent with these results differential gene expression obtained using RNA sequencing studies, and validated by RT PCR and immunoblot analysis, indicate that MASTL inhibition induces cell cycle arrest and apoptotic-related genes. Furthermore, metastasis- and inflammation- related genes are downregulated. Thus, the deregulation of MASTL signaling pathways on targeted inhibition of its kinase activity is revealed. This study lays a strong foundation for investigating FLV as a lead compound in breast cancer therapeutics.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Seed coat transcriptomic profiling of 5-593, a genotype important for genetic studies of seed coat color and patterning in common bean ( Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) market classes have distinct seed coat colors, which are directly related to the diverse flavonoids found in the mature seed coat. To understand and elucidate the molecular mechanisms underlying the regulation of seed coat color, RNA-Seq data was collected from the black bean 5-593 and used for a differential gene expression and enrichment analysis from four different seed coat color development stages. 5-593 carries dominant alleles for 10 of the 11 major genes that control seed coat color and expression and has historically been used to develop introgression lines used for seed coat genetic analysis. Pairwise comparison among the four stages identified 6,294 differentially expressed genes (DEGs) varying from 508 to 5,780 DEGs depending on the compared stages. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that phenylpropanoid biosynthesis, flavonoid biosynthesis, and plant hormone signal transduction comprised the principal pathways expressed during bean seed coat pigment development. Transcriptome analysis suggested that most structural genes for flavonoid biosynthesis and some potential regulatory genes were significantly differentially expressed. Further studies detected 29 DEGs as important candidate genes governing the key enzymatic flavonoid biosynthetic pathways for common bean seed coat color development. Additionally, four gene models, Pv5-593.02G016100, 593.02G078700, Pv5-593.02G090900, and Pv5-593.06G121300, encode MYB-like transcription factor family protein were identified as strong candidate regulatory genes in anthocyanin biosynthesis which could regulate the expression levels of some important structural genes in flavonoid biosynthesis pathway. These findings provide a framework to draw new insights into the molecular networks underlying common bean seed coat pigment development.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Structured sequences emerge from random pool when replicated by templated ligation

The central question in the origin of life is to understand how structure can emerge from randomness. The Eigen theory of replication states, for sequences that are copied one base at a time, that the replication fidelity has to surpass an error threshold to avoid that replicated specific sequences become random because of the incorporated replication errors [M. Eigen, Naturwissenschaften 58 (10), 465–523 (1971)]. Here, we showed that linking short oligomers from a random sequence pool in a templated ligation reaction reduced the sequence space of product strands. We started from 12-mer oligonucleotides with two bases in all possible combinations and triggered enzymatic ligation under temperature cycles. Surprisingly, we found the robust creation of long, highly structured sequences with low entropy. At the ligation site, complementary and alternating sequence patterns developed. However, between the ligation sites, we found either an A-rich or a T-rich sequence within a single oligonucleotide. Our modeling suggests that avoidance of hairpins was the likely cause for these two complementary sequence pools. What emerged was a network of complementary sequences that acted both as templates and substrates of the reaction. This self-selecting ligation reaction could be restarted by only a few majority sequences. The findings showed that replication by random templated ligation from a random sequence input will lead to a highly structured, long, and nonrandom sequence pool. This is a favorable starting point for a subsequent Darwinian evolution searching for higher catalytic functions in an RNA world scenario.

60 APPLIED LIFE SCIENCES↗
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↗
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↗