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.

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.

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.

Tying up the Loose Ends: A Mathematically Knotted Protein.

Knots have attracted scientists in mathematics, physics, biology, and engineering. Long flexible thin strings easily knot and tangle as experienced in our daily life. Similarly, long polymer chains inevitably tend to get trapped into knots. Little is known about their formation or function in proteins despite >1,000 knotted proteins identified in nature. However, these protein knots are not mathematical knots with their backbone polypeptide chains because of their open termini, and the presence of a "knot" depends on the algorithm used to create path closure. Furthermore, it is generally not possible to control the topology of the unfolded states of proteins, therefore making it challenging to characterize functional and physicochemical properties of knotting in any polymer. Covalently linking the amino and carboxyl termini of the deeply trefoil-knotted YibK from Pseudomonas aeruginosa allowed us to create the truly backbone knotted protein by enzymatic peptide ligation. Moreover, we produced and investigated backbone cyclized YibK without any knotted structure. Thus, we could directly probe the effect of the backbone knot and the decrease in conformational entropy on protein folding. The backbone cyclization did not perturb the native structure and its cofactor binding affinity, but it substantially increased the thermal stability and reduced the aggregation propensity. The enhanced stability of a backbone knotted YibK could be mainly originated from an increased ruggedness of its free energy landscape and the destabilization of the denatured state by backbone cyclization with little contribution from a knot structure. Despite the heterogeneity in the side-chain compositions, the chemically unfolded cyclized YibK exhibited several macroscopic physico-chemical attributes that agree with theoretical predictions derived from polymer physics.

Synchronous RNA conformational changes trigger ordered phase transitions in crystals.

Time-resolved studies of biomacromolecular crystals have been limited to systems involving only minute conformational changes within the same lattice. Ligand-induced changes greater than several angstroms, however, are likely to result in solid-solid phase transitions, which require a detailed understanding of the mechanistic interplay between conformational and lattice transitions. Here we report the synchronous behavior of the adenine riboswitch aptamer RNA in crystal during ligand-triggered isothermal phase transitions. Direct visualization using polarized video microscopy and atomic force microscopy shows that the RNA molecules undergo cooperative rearrangements that maintain lattice order, whose cell parameters change distinctly as a function of time. The bulk lattice order throughout the transition is further supported by time-resolved diffraction data from crystals using an X-ray free electron laser. The synchronous molecular rearrangements in crystal provide the physical basis for studying large conformational changes using time-resolved crystallography and micro/nanocrystals.

Topologically knotted deubiquitinases exhibit unprecedented mechanostability to withstand the proteolysis by an AAA+ protease.

More than one thousand knotted protein structures have been identified so far, but the functional roles of these knots remain elusive. It has been postulated that backbone entanglement may provide additional mechanostability. Here, we employed a bacterial proteasome, ClpXP, to mechanically unfold 52-knotted human ubiquitin C-terminal hydrolase (UCH) paralogs from their C-termini, followed by processive translocation into the proteolytic chamber for degradation. Our results revealed unprecedentedly slow kinetics of ClpXP-mediated proteolysis for the proteasome-associated UCHL5: ten thousand times slower than that of a green fluorescence protein (GFP), which has a comparable size to the UCH domain but much higher chemical and thermal stabilities. The ClpXP-dependent mechanostability positively correlates with the intrinsic unfolding rates of the substrates, spanning over several orders of magnitude for the UCHs. The broad range of mechanostability within the same protein family may be associated with the functional requirements for their differential malleabilities.

A Natively Monomeric Deubiquitinase UCH-L1 Forms Highly Dynamic but Defined Metastable Oligomeric Folding Intermediates.

Oligomerization of misfolded protein species is implicated in many human disorders. Here we showed by size-exclusion chromatography-coupled multiangle light scattering (SEC-MALS) and small-angle X-ray scattering (SEC-SAXS) that urea-induced folding intermediate of human ubiquitin C-terminal hydrolase, UCH-L1, can form well-defined dimers and tetramers under denaturing conditions despite being highly disordered. Introduction of a Parkinson disease-associated mutation, I93M, resulted in increased aggregation propensity and formation of irreversible precipitants in the presence of a moderate amount of urea. Since UCH-L1 exhibits highly populated partially unfolded forms under native conditions that resemble urea-induced folding intermediates, it is likely that these metastable dimers and tetramers can form under physiological conditions. Our findings highlighted the unique strength of integrated SEC-MALS/SAXS in quantitative analyses of the structure and dynamics of oligomeric folding intermediates that enabled us to extract information that is inaccessible to conventional biophysical techniques.

Entropic stabilization of a deubiquitinase provides conformational plasticity and slow unfolding kinetics beneficial for functioning on the proteasome.

Human ubiquitin C-terminal hydrolyase UCH-L5 is a topologically knotted deubiquitinase that is activated upon binding to the proteasome subunit Rpn13. The length of its intrinsically disordered cross-over loop is essential for substrate recognition. Here, we showed that the catalytic domain of UCH-L5 exhibits higher equilibrium folding stability with an unfolding rate on the scale of 10-8 s-1, over four orders of magnitudes slower than its paralogs, namely UCH-L1 and -L3, which have shorter cross-over loops. NMR relaxation dynamics analysis confirmed the intrinsic disorder of the cross-over loop. Hydrogen deuterium exchange analysis further revealed a positive correlation between the length of the cross-over loop and the degree of local fluctuations, despite UCH-L5 being thermodynamically and kinetically more stable than the shorter UCHs. Considering the role of UCH-L5 in removing K48-linked ubiquitin to prevent proteasomal degradation of ubiquitinated substrates, our findings offered mechanistic insights into the evolution of UCH-L5. Compared to its paralogs, it is entropically stabilized to withstand mechanical unfolding by the proteasome while maintaining structural plasticity. It can therefore accommodate a broad range of substrate geometries at the cost of unfavourable entropic loss.

Familial Mutations and Post-translational Modifications of UCH-L1 in Parkinson's Disease and Neurodegenerative Disorders.

Parkinson's disease (PD) is one of the most common progressive neurodegenerative disorders in modern society. The disease involves many genetic risk factors as well as a sporadic pathogenesis that is age- and environment-dependent. Of particular interest is the formation of intra-neural fibrillar aggregates, namely Lewy bodies (LBs), the histological hallmark of PD, which results from aberrant protein homeostasis or misfolding that results in neurotoxicity. A better understanding of the molecular mechanism and composition of these cellular inclusions will help shed light on the progression of misfolding-associated neurodegenerative disorders. Ubiquitin carboxyl-terminal hydrolase L1 (UCH-L1) is found to co-aggregate with α-synuclein (αS), the major component of LBs. Several familial mutations of UCH-L1, namely p.Ile93Met (p.I93M), p.Glu7Ala (p.E7A), and p.Ser18Tyr (p.S18Y), are associated with PD and other neurodegenerative disorders. Here, we review recent progress and recapitulate the impact of PD-associated mutations of UCH-L1 in the context of their biological functions gleaned from biochemical and biophysical studies. Finally, we summarize the effect of these genetic mutations and post-translational modifications on the association of UCHL1 and PD in terms of loss of cellular functions or gain of cellular toxicity.

The Knotted Protein UCH-L1 Exhibits Partially Unfolded Forms under Native Conditions that Share Common Structural Features with Its Kinetic Folding Intermediates.

The human ubiquitin C-terminal hydrolase, UCH-L1, is an abundant neuronal deubiquitinase that is associated with Parkinson's disease. It contains a complex Gordian knot topology formed by the polypeptide chain alone. Using a combination of fluorescence-based kinetic measurements, we show that UCH-L1 has two distinct kinetic folding intermediates that are transiently populated on parallel pathways between the denatured and native states. NMR hydrogen-deuterium exchange (HDX) experiments indicate the presence of partially unfolded forms (PUFs) of UCH-L1 under native conditions. HDX measurements as a function of urea concentration were used to establish the structure of the PUFs and pulse-labelled HDX NMR was used to show that the PUFs and the folding intermediates are likely the same species. In both cases, a similar stable core encompassing most of the central ß-sheet is highly structured and α-helix 3, which is partially formed, packs against it. In contrast to the stable ß-sheet core, the peripheral α-helices display significant local fluctuations leading to rapid exchange. The results also suggest that the main difference between the two kinetic intermediates is structure and packing of α-helices 3 and 7 and the degree of structure in ß-strand 5. Together, the fluorescence and NMR results establish that UCH-L1 neither folds through a continuum of pathways nor by a single discrete pathway. Its folding is complex, the ß-sheet core forms early and is present in both intermediate states, and the rate-limiting step which is likely to involve the threading of the chain to form the 52-knot occurs late on the folding pathway.

Random-Coil Behavior of Chemically Denatured Topologically Knotted Proteins Revealed by Small-Angle X-ray Scattering.

Recent studies on the mechanisms by which topologically knotted proteins attain their natively knotted structures have intrigued theoretical and experimental biophysicists. Of particular interest is the finding that YibK and YbeA, two small trefoil knotted proteins, remain topologically knotted in their chemically denatured states. Using small-angle X-ray scattering (SAXS), we examine whether these chemically denatured knotted proteins are different from typical random coils. By revisiting the scaling law of radius of gyration (Rg) as a function of polypeptide chain length for chemically denatured proteins and natively folded proteins, we find that the chemically denatured knotted proteins in fact follow the same random coil-like behavior, suggesting that the formation of topological protein knots do not necessarily require global compaction while the loosely knotted polypeptide chains are capable of maintaining the correct chirality without defined secondary or tertiary structures.

A streamlined method for preparing split intein for NMR study.

A protein ligase, intein, mediates a protein-splicing reaction. It can be split into two complementary fragments and reconstituted as a whole intein scaffold to perform protein trans-splicing. To understand the association of intein fragments and the splicing mechanism, it is necessary to produce a large quantity of split intein for structural study. Conventionally, two fragments are prepared separately and assembled in solution, but severe aggregation of intein fragments occurs, and precise control of the relative concentration of each fragment is difficult. Here, we present a streamlined method to incorporate a circular permutation concept into the production of split intein. By circular permutation of the native split Nostoc punctiforme DnaE intein (NpuInt), a new backbone opening is relocated to the native split site at residue 102. As the protein splicing activity is preserved, the expressed NpuInt can immediately self-cleave into a two-piece split NpuInt. Because of a tight association between the two complementary fragments, split NpuInt can be purified in one step. The idea is simple and applicable to other split inteins. Employing the new preparation, we use NMR spectra to assign the backbone and side chain resonances for the native split NpuInt.

Circular permutation prediction reveals a viable backbone disconnection for split proteins: an approach in identifying a new functional split intein.

Split-protein systems have emerged as a powerful tool for detecting biomolecular interactions and reporting biological reactions. However, reliable methods for identifying viable split sites are still unavailable. In this study, we demonstrated the feasibility that valid circular permutation (CP) sites in proteins have the potential to act as split sites and that CP prediction can be used to search for internal permissive sites for creating new split proteins. Using a protein ligase, intein, as a model, CP predictor facilitated the creation of circular permutants in which backbone opening imposes the least detrimental effects on intein folding. We screened a series of predicted intein CPs and identified stable and native-fold CPs. When the valid CP sites were introduced as split sites, there was a reduction in folding enthalpy caused by the new backbone opening; however, the coincident loss in entropy was sufficient to be compensated, yielding a favorable free energy for self-association. Since split intein is exploited in protein semi-synthesis, we tested the related protein trans-splicing (PTS) activities of the corresponding split inteins. Notably, a novel functional split intein composed of the N-terminal 36 residues combined with the remaining C-terminal fragment was identified. Its PTS activity was shown to be better than current reported two-piece intein with a short N-terminal segment. Thus, the incorporation of in silico CP prediction facilitated the design of split intein as well as circular permutants.