Structural basis of MALAT1 RNA maturation and mascRNA biogenesis.

The metastasis-associated lung adenocarcinoma transcript 1 (MALAT1) long noncoding RNA (lncRNA) has key roles in regulating transcription, splicing, tumorigenesis, etc. Its maturation and stabilization require precise processing by RNase P, which simultaneously initiates the biogenesis of a 3' cytoplasmic MALAT1-associated small cytoplasmic RNA (mascRNA). mascRNA was proposed to fold into a transfer RNA (tRNA)-like secondary structure but lacks eight conserved linking residues required by the canonical tRNA fold. Here we report crystal structures of human mascRNA before and after processing, which reveal an ultracompact, quasi-tRNA-like structure. Despite lacking all linker residues, mascRNA faithfully recreates the characteristic 'elbow' feature of tRNAs to recruit RNase P and ElaC homolog protein 2 (ELAC2) for processing, which exhibit distinct substrate specificities. Rotation and repositioning of the D-stem and anticodon regions preclude mascRNA from aminoacylation, avoiding interference with translation. Therefore, a class of metazoan lncRNA loci uses a previously unrecognized, unusually streamlined quasi-tRNA architecture to recruit select tRNA-processing enzymes while excluding others to drive bespoke RNA biogenesis, processing and maturation.

Overexpression of a BR inactivating enzyme gene GhPAG1 impacts eggplant fruit development and anthocyanin accumulation mainly by altering hormone homeostasis.

Brassinosteroids (BRs) function importantly in plant growth and development, but the roles in regulating fruit development and anthocyanin pigmentation remain unclear. Eggplant (Solanum melongena L.) is an important Solanaceae vegetable crop rich in anthocyanins. The fruit size and coloration are important agronomic traits for eggplant breeding. In this study, transgenic eggplant exhibiting endogenous BRs deficiency was created by overexpressing a heterologous BRs-inactivating enzyme gene GhPAG1 driven by CaMV 35 S promoter. 35 S::GhPAG1 eggplant exhibited severe dwarfism, reduced fruit size, and less anthocyanin accumulation. Microscopic observation showed that the cell size of 35 S::GhPAG1 eggplant was significantly reduced compared to WT. Furthermore, the levels of IAA, ME-IAA, and active JAs (JA, JA-ILE, and H2JA) all decreased in 35 S::GhPAG1 eggplant fruit. RNA-Seq analyses showed a decrease in the expression of genes involved in cell elongation, auxin signaling, and JA signaling. Besides, overexpression of GhPAG1 significantly downregulated anthocyanin biosynthetic genes and associated transcription regulators. Altogether, these results strongly suggest that endogenous brassinosteroid deficiency arising from GhPAG1 overexpression impacts eggplant fruit development and anthocyanin coloration mainly by altering hormone homeostasis.

Direct observation of tRNA-chaperoned folding of a dynamic mRNA ensemble.

T-box riboswitches are multi-domain noncoding RNAs that surveil individual amino acid availabilities in most Gram-positive bacteria. T-boxes directly bind specific tRNAs, query their aminoacylation status to detect starvation, and feedback control the transcription or translation of downstream amino-acid metabolic genes. Most T-boxes rapidly recruit their cognate tRNA ligands through an intricate three-way stem I-stem II-tRNA interaction, whose establishment is not understood. Using single-molecule FRET, SAXS, and time-resolved fluorescence, we find that the free T-box RNA assumes a broad distribution of open, semi-open, and closed conformations that only slowly interconvert. tRNA directly binds all three conformers with distinct kinetics, triggers nearly instantaneous collapses of the open conformations, and returns the T-box RNA to their pre-binding conformations upon dissociation. This scissors-like dynamic behavior is enabled by a hinge-like pseudoknot domain which poises the T-box for rapid tRNA-induced domain closure. This study reveals tRNA-chaperoned folding of flexible, multi-domain mRNAs through a Venus flytrap-like mechanism.

Capturing heterogeneous conformers of cobalamin riboswitch by cryo-EM.

RNA conformational heterogeneity often hampers its high-resolution structure determination, especially for large and flexible RNAs devoid of stabilizing proteins or ligands. The adenosylcobalamin riboswitch exhibits heterogeneous conformations under 1 mM Mg2+ concentration and ligand binding reduces conformational flexibility. Among all conformers, we determined one apo (5.3 Å) and four holo cryo-electron microscopy structures (overall 3.0-3.5 Å, binding pocket 2.9-3.2 Å). The holo dimers exhibit global motions of helical twisting and bending around the dimer interface. A backbone comparison of the apo and holo states reveals a large structural difference in the P6 extension position. The central strand of the binding pocket, junction 6/3, changes from an 'S'- to a 'U'-shaped conformation to accommodate ligand. Furthermore, the binding pocket can partially form under 1 mM Mg2+ and fully form under 10 mM Mg2+ within the bound-like structure in the absence of ligand. Our results not only demonstrate the stabilizing ligand-induced conformational changes in and around the binding pocket but may also provide further insight into the role of the P6 extension in ligand binding and selectivity.

Determining structures of individual RNA conformers using atomic force microscopy images and deep neural networks.

The vast percentage of the human genome is transcribed into RNA, many of which contain various structural elements and are important for functions. RNA molecules are conformationally heterogeneous and functionally dyanmics1, even when they are structured and well-folded2, which limit the applicability of methods such as NMR, crystallography, or cryo-EM. Moreover, because of the lack of a large structure RNA database, and no clear correlation between sequence and structure, approaches like AlphaFold3 for protein structure prediction, do not apply to RNA. Therefore determining the structures of heterogeneous RNA is an unmet challenge. Here we report a novel method of determining RNA three-dimensional topological structures using deep neural networks and atomic force microscopy (AFM) images of individual RNA molecules in solution. Owing to the high signal-to-noise ratio of AFM, our method is ideal for capturing structures of individual conformationally heterogeneous RNA. We show that our method can determine 3D topological structures of any large folded RNA conformers, from ~ 200 to ~ 420 residues, the size range that most functional RNA structures or structural elements fall into. Thus our method addresses one of the major challenges in frontier RNA structural biology and may impact our fundamental understanding of RNA structure.

Architectural basis for cylindrical self-assembly governing Plk4-mediated centriole duplication in human cells.

Proper organization of intracellular assemblies is fundamental for efficient promotion of biochemical processes and optimal assembly functionality. Although advances in imaging technologies have shed light on how the centrosome is organized, how its constituent proteins are coherently architected to elicit downstream events remains poorly understood. Using multidisciplinary approaches, we showed that two long coiled-coil proteins, Cep63 and Cep152, form a heterotetrameric building block that undergoes a stepwise formation into higher molecular weight complexes, ultimately generating a cylindrical architecture around a centriole. Mutants defective in Cep63•Cep152 heterotetramer formation displayed crippled pericentriolar Cep152 organization, polo-like kinase 4 (Plk4) relocalization to the procentriole assembly site, and Plk4-mediated centriole duplication. Given that the organization of pericentriolar materials (PCM) is evolutionarily conserved, this work could serve as a model for investigating the structure and function of PCM in other species, while offering a new direction in probing the organizational defects of PCM-related human diseases.

Prognostic Analysis and Biomarkers Identification of Immune Infiltration in Early and Late Stage Hepatocellular Carcinoma Based on TCGA Data.

Background: Hepatocellular carcinoma (HCC) is a major cause of cancer death in the world. The aim of this study was to establish a new model to predict the prognosis of HCC. Materials and Methods: The mRNA, miRNA and lncRNA expression profiles of early (stage I-II) and late (stage III-IV) stage HCC patients were acquired from The Cancer Genome Atlas (TCGA) database. The differentially expressed mRNAs (DEmRNAs), miRNAs (DEmiRNAs) and lncRNAs (DElncRNAs) were identified between early and late stage HCC. Key molecules associated with the prognosis, and important immune cell types in HCC were identified. The nomogram based on incorporating age, gender, stage, and all important factors was constructed to predict the survival of HCC. Results: A total of 1516 DEmRNAs, 97 DEmiRNAs and 87 DElncRNAs were identified. A DElncRNA-DEmiRNA-DEmRNA regulatory network including 78 mRNAs, 50 miRNAs and 1 lncRNA was established. Among the regulatory network, 11 molecules were significantly correlated with the prognosis of HCC based on Lasso regression analysis. Then, Preadipocytes and 3 survival-associated DEmRNAs were identified as crucial biomarkers. Subsequently, a nomogram with a differentiation degree of 0.758, including 1 immune cell, 11 mRNAs and 3 miRNAs, was generated. Conclusion: Our study constructed a model by incorporating clinical information, significant biomarkers and immune cells to predict the survival of HCC, which achieved a good performance.

Crystal structure of Escherichia coli thiamine pyrophosphate-sensing riboswitch in the apo state.

The thiamine pyrophosphate (TPP)-sensing riboswitch is one of the earliest discovered and most widespread riboswitches. Numerous structural studies have been reported for this riboswitch bound with various ligands. However, the ligand-free (apo) structure remains unknown. Here, we report a 3.1 Å resolution crystal structure of Escherichia coli TPP riboswitch in the apo state, which exhibits an extended, Y-shaped conformation further supported by small-angle X-ray scattering data and driven molecular dynamics simulations. The loss of ligand interactions results in helical uncoiling of P5 and disruption of the key tertiary interaction between the sensory domains. Opening of the aptamer propagates to the gene-regulatory P1 helix and generates the key conformational flexibility needed for the switching behavior. Much of the ligand-binding site at the three-way junction is unaltered, thereby maintaining a partially preformed pocket. Together, these results paint a dynamic picture of the ligand-induced conformational changes in TPP riboswitches that confer conditional gene regulation.

SmCIP7, a COP1 interactive protein, positively regulates anthocyanin accumulation and fruit size in eggplant.

In higher plants, COP1 (Constitutively Photomorphogenic 1) acts as a central regulator of light-signaling networks and globally conditions the target proteins via the ubiquitin-proteasome pathway. However, the function of COP1-interacting proteins in light-regulated fruit coloration and development remains unknown in Solanaceous plants. Here, a COP1-interacting protein-encoding gene, SmCIP7, expressed specifically in the eggplant (Solanum melongena L.) fruit, was isolated. Gene-specific silencing of SmCIP7 using RNA interference (RNAi) significantly altered fruit coloration, fruit size, flesh browning, and seed yield. SmCIP7-RNAi fruits showed evident repression of the accumulation of anthocyanins and chlorophyll, indicating functional similarities between SmCIP7 and AtCIP7. However, the reduced fruit size and seed yield indicated SmCIP7 had evolved a distinctly new function. With the comprehensive application of HPLC-MS, RNA-seq, qRT-PCR, Y2H, BiFC, LCI, and dual-luciferase reporter system (DLR™), it was found that SmCIP7, a COP1 interactive protein in light signaling promoted anthocyanin accumulation, probably by regulating the transcription of SmTT8. Additionally, the drastic up-regulation of SmYABBY1, a homologous gene of SlFAS, might account for the strongly retarded fruit growth in SmCIP7-RNAi eggplant. Altogether, this study proved that SmCIP7 is an essential regulatory gene to modulate fruit coloration and development, serving as a key gene locus in eggplant molecular breeding.

Visualizing RNA conformational and architectural heterogeneity in solution.

RNA flexibility is reflected in its heterogeneous conformation. Through direct visualization using atomic force microscopy (AFM) and the adenosylcobalamin riboswitch aptamer domain as an example, we show that a single RNA sequence folds into conformationally and architecturally heterogeneous structures under near-physiological solution conditions. Recapitulated 3D topological structures from AFM molecular surfaces reveal that all conformers share the same secondary structural elements. Only a population-weighted cohort, not any single conformer, including the crystal structure, can account for the ensemble behaviors observed by small-angle X-ray scattering (SAXS). All conformers except one are functionally active in terms of ligand binding. Our findings provide direct visual evidence that the sequence-structure relationship of RNA under physiologically relevant solution conditions is more complex than the one-to-one relationship for well-structured proteins. The direct visualization of conformational and architectural ensembles at the single-molecule level in solution may suggest new approaches to RNA structural analyses.

Combining Biophysical Methods for Structure-Function Analyses of RNA in Solution.

RNA-level regulation by riboswitches relies on the specific binding of small metabolites to the aptamer domain to trigger substantial conformational changes that affect transcription or translation. Although several biophysical methods have been employed to study such RNAs, the utility of any one single method is limited. Hybrid approaches, therefore, are essential to better characterize these intrinsically dynamic molecules and elucidate their regulatory mechanisms driven by ligand-induced conformational changes. This chapter outlines procedures for biochemical and biophysical characterization of RNA that employs a combination of solution-based methods: isothermal titration calorimetry (ITC), small-angle X-ray scattering (SAXS), and atomic force microscopy (AFM). Collectively, these tools provide a semi-quantitative assessment of the thermodynamics associated with ligand binding and subsequent conformational changes.

Mix-and-Inject Serial Femtosecond Crystallography to Capture RNA Riboswitch Intermediates.

Time-resolved structure determination of macromolecular conformations and ligand-bound intermediates is extremely challenging, particularly for RNA. With rapid technological advances in both microfluidic liquid injection and X-ray free electron lasers (XFEL), a new frontier has emerged in time-resolved crystallography whereby crystals can be mixed with ligand and then probed with X-rays (mix-and-inject) in real time and at room temperature. This chapter outlines the basic setup and procedures for mix-and-inject experiments for recording time-resolved crystallographic data of riboswitch RNA reaction states using serial femtosecond crystallography (SFX) and an XFEL.

Deriving RNA topological structure from SAXS.

Structures of well-folded RNA molecules can be determined with atomic resolution by either X-ray crystallography, cryo-EM, or NMR spectroscopy, but those of conformationally-flexible RNAs often are difficult to study with these methods. However, flexible RNAs have biological relevance and likely represent the majority of the RNA conformational space. Due to the high electron density of the phosphate-sugar backbone, RNA is very sensitive to small-angle X-ray scattering (SAXS), and SAXS data can be recorded with sub-µM concentrations and under near-physiological solution conditions without the need for labeling. For these reasons, SAXS has significant advantages over other techniques for obtaining global structural information of flexible RNAs in the form of molecular envelopes or low-resolution topological structural models. The SAXS-derived information is extremely valuable for bridging secondary structure data, often determined by other techniques, with a three-dimensional structure description. In this chapter, we present a detailed account of the principle, algorithms, and experimental and computational protocols for topological structure determination of RNA molecules in solution. To illustrate the applications of the methodology, we provide several case studies that cover a broad spectrum of the RNA conformational landscape.

Developing methods to study conformational changes in RNA crystals using a photocaged ligand.

Crystallographic observation of structural changes in real time requires that those changes be uniform both spatially and temporally. A primary challenge with time-resolved ligand-mixing diffraction experiments is asynchrony caused by variable factors, such as efficiency of mixing, rate of diffusion, crystal size, and subsequently, conformational heterogeneity. One method of minimizing such variability is use of a photolabile caged ligand, which can fully saturate the crystal environment (spatially), and whose photoactivation can rapidly (temporally) trigger the reaction in a controlled manner. Our recently published results on a ligand-mixing experiment using time-resolved X-ray crystallography (TRX) with an X-ray free electron laser (XFEL) demonstrated that large conformational changes upon ligand binding resulted in a solid-to-solid phase transition (SSPT), while maintaining Bragg diffraction. Here we investigate this SSPT by polarized video microscopy (PVM) after light-triggered release of a photo-caged adenine (pcADE). In general, the mean transition times and transition widths of the SSPT were less dependent on crystal size than what was observed in previous PVM studies with direct ADE mixing. Instead, the photo-induced transition appears to be heavily influenced by the equilibrium between caged and uncaged ADE due to relatively low sample exposure and uncaging efficiency. Nevertheless, we successfully demonstrate a method for the characterization of phase transitions in RNA crystals that are inducible with a photocaged ligand. The transition data for three crystals of different sizes were then applied to kinetic analysis by fitting to the known four-state model associated with ligand-induced conformational changes, revealing an apparent concentration of uncaged ADE in crystal of 0.43-0.46 mM. These results provide further insight into approaches to study time-resolved ligand-induced conformational changes in crystals, and in particular, highlight the feasibility of triggering phase transitions using a light-inducible system. Developing such approaches may be paramount for the rapidly emerging field of time-resolved crystallography.

Taf2 mediates DNA binding of Taf14.

The assembly and function of the yeast general transcription factor TFIID complex requires specific contacts between its Taf14 and Taf2 subunits, however, the mechanism underlying these contacts remains unclear. Here, we determined the molecular and structural basis by which the YEATS and ET domains of Taf14 bind to the C-terminal tail of Taf2 and identified a unique DNA-binding activity of the linker region connecting the two domains. We show that in the absence of ligands the linker region of Taf14 is occluded by the surrounding domains, and therefore the DNA binding function of Taf14 is autoinhibited. Binding of Taf2 promotes a conformational rearrangement in Taf14, resulting in a release of the linker for the engagement with DNA and the nucleosome. Genetic in vivo data indicate that the association of Taf14 with both Taf2 and DNA is essential for transcriptional regulation. Our findings provide a basis for deciphering the role of individual TFIID subunits in mediating gene transcription.

CYP26A1 Is a Novel Cancer Biomarker of Pancreatic Carcinoma: Evidence from Integration Analysis and In Vitro Experiments.

Background: CYP26A1 has been reported in multiple cancers. However, the role of CYP26A1 in pancreatic cancer (PC) has not been explored. Method: The public data used for this study was obtained from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and Cancer Cell Line Encyclopedia (CCLE) cell lines. CCK8, colony formation, and EdU assay were used to assess the proliferation ability of cancer cells. Transwell and wound healing assays were used to evaluate the invasion and migration ability of cancer cells. qRT-PCR and western blot assays were used to analyze the RNA and protein level of genes. Survival package was used for prognosis analysis. Gene Set Enrichment Analysis (GSEA) was used to identify biological pathway differences between two groups. ssGSEA analysis was used to quantify the immune microenvironment in PC tissue. GDSC and TIDE analyses were used for sensitivity analysis of chemotherapy and immunotherapy. Results: Our results showed that CYP26A1 was overexpressed in PC tissue and cell lines. Meanwhile, metastatic PC cell lines tend to have a higher CYP26A1 level compared with the primary PC cell lines based on CCLE data. Moreover, CYP26A1 was associated with worse clinical features. Also, we found that CYP26A1 had a satisfactory efficiency in predicting overall survival, disease-specific survival, and progression-free interval of PC patients, independent of other clinical features. In vitro experiments indicated that CYP26A1 could significantly facilitate the proliferation, invasion, and migration ability of PC cells. GSEA showed that the pathways of angiogenesis, E2F target, MYC target, mTORC signaling, G2M checkpoint, and epithelial-mesenchymal transition were activated in high CYP26A1 patients. Immune infiltration analysis showed that CYP26A1 was positively correlated with macrophages, Th1 cells, and Treg cells, but negatively correlated with Th17 cells. TIDE analysis showed that non_responder patients had a higher CYP26A1 level compared with predicted responder patients of immunotherapy. Drug sensitivity analysis and assay showed that CYP26A1 could increase the chemotherapy sensitivity of gemcitabine. Conclusions: In summary, CYP26A1 promotes PC progression and is a novel biomarker of PC, with potential for clinical application.

Relative domain orientation of the L289K HIV-1 reverse transcriptase monomer.

HIV-1 reverse transcriptase (RT) is a heterodimer comprised p66 and p51 subunits (p66/p51). Several single amino acid substitutions in RT, including L289K, decrease p66/p51 dimer affinity, and reduce enzymatic functioning. Here, small-angle X-ray scattering (SAXS) with proton paramagnetic relaxation enhancement (PRE), 19 F site-specific NMR, and size exclusion chromatography (SEC) were performed for the p66 monomer with the L289K mutation, p66L289K . NMR and SAXS experiments clearly elucidated that the thumb and RNH domains in the monomer do not rigidly interact with each other but are spatially close to the RNH domain. Based on this structural model of the monomer, p66L289K and p51 were predicted to form a heterodimer while p66 and p51L289K not. We tested this hypothesis by SEC analysis of p66 and p51 containing L289K in different combinations and clearly demonstrated that L289K substitution in the p51 subunit, but not in the p66 subunit, reduces p66/p51 formation. Based on the derived monomer model and the importance of the inter-subunit RNH-thumb domain interaction in p66/p51, validated by SEC, the mechanism of p66 homodimer formation was discussed.

The low-resolution structural models of hepatitis C virus RNA subdomain 5BSL3.2 and its distal complex with domain 3'X point to conserved regulatory mechanisms within the Flaviviridae family.

Subdomain 5BSL3.2 of hepatitis C virus RNA lies at the core of a network of distal RNA-RNA contacts that connect the 5' and 3' regions of the viral genome and regulate the translation and replication stages of the viral cycle. Using small-angle X-ray scattering and NMR spectroscopy experiments, we have determined at low resolution the structural models of this subdomain and its distal complex with domain 3'X, located at the 3'-terminus of the viral RNA chain. 5BSL3.2 adopts a characteristic 'L' shape in solution, whereas the 5BSL3.2-3'X distal complex forms a highly unusual 'Y'-shaped kissing junction that blocks the dimer linkage sequence of domain 3'X and promotes translation. The structure of this complex may impede an effective association of the viral polymerase with 5BSL3.2 and 3'X to start negative-strand RNA synthesis, contributing to explain the likely mechanism used by these sequences to regulate viral replication and translation. In addition, sequence and shape features of 5BSL3.2 are present in functional RNA motifs of flaviviruses, suggesting conserved regulatory processes within the Flaviviridae family.

Dependence of phase transition uniformity on crystal sizes characterized using birefringence.

Solid-solid phase transitions (SSPTs) have been widely observed in crystals of organic or inorganic small-molecules. Although SSPTs in macromolecular crystals have been reported, the majority involve local atomic changes, such as those induced by changes in hydration. SSPTs driven by large conformational changes, however, can be more difficult to characterize since they often significantly disrupt lattice packing interactions. Such drastic changes make the cooperativity of molecular motion at the atomic level less easily achieved and more dependent on intrinsic properties of the crystal that define lattice order. Here, we investigate the effect of crystal size on the uniformity of SSPT in thin plate-like crystals of the adenine riboswitch aptamer RNA (riboA) by monitoring changes in crystal birefringence upon the diffusion of adenine ligand. The birefringence intensity is directly related to molecular order and the concurrent changes to polarizability of molecules that results from structural changes throughout the phase transition. The riboA crystals were loosely grouped into three categories (small, medium, and large) based on the surface area of the crystal plates. The time width of transition increased as a function of crystal size, ranging from ∼13 s for small crystals to ∼40 s for the largest crystal. Whereas the transitions in small crystals (<10 µm2) were mostly uniform throughout, the medium and large crystals exhibited large variations in the time and width of the transition peak depending on the region of the crystal being analyzed. Our study provides insight into the spatiotemporal behavior of phase transitions in crystals of biological molecules and is of general interest to time-resolved crystallographic studies, where the kinetics of conformational changes may be governed by the kinetics of an associated SSPT.