RNA sectors and allosteric function within the ribosome

Coevolving residues distant from the ligand binding site are involved in GAF domain function

Ligand binding to GAF domains regulates the activity of associated catalytic domains in various proteins, such as the cGMP-hydrolyzing catalytic domain of phosphodiesterase 5 (PDE5) activated by cGMP binding to GAFa domain. However, the specific residues involved and the mechanism of GAF domain function remain unclear. Here, we combine computational and experimental approaches to demonstrate that two highly coevolving residues, L267 and F295, distant from the ligand binding site, play a critical role in GAF domain allostery. Statistical Coupling Analysis (SCA) of GAF domain sequences identified these residues, and molecular dynamics (MD) simulations of both apo and holo forms of wild-type and mutant (L267A, F295A) PDE5 GAFa domains revealed significant changes in structural dynamics and cGMP interaction. Mutational incorporation into a Bioluminescence Resonance Energy Transfer (BRET)-based biosensors, which detects ligand-induced conformational changes, showed altered GAF domain conformation and increased EC50 for cGMP-induced conformational changes. Similar effects were observed in full-length PDE5 and the GAF domain fluorescent protein, miRFP670nano3. Structural analysis of conformers observed in MD simulations suggested a mechanism by which these coevolving residues influence GAF domain allostery. Our findings provide insight into the role of distant residues in GAF domain function and may enhance understanding of allostery in proteins.

Recognition Mechanism of RNA by TLR13: Structural Insights and Implications for Immune Activation

RNA serves as a distinctive pathogen-associated molecular pattern (PAMP) that plays a critical role in innate immunity. However, the specific mechanisms of RNA recognition remain largely unexplored, especially given RNA's vulnerability to degradation and the absence of sequence specificity in most RNA recognition receptors. Notably, Toll-like receptor 13 (TLR13) is capable of detecting a conserved RNA sequence, RNA15 (2054-2068, ACG GAA AGA CCC CGU), within bacterial 23S rRNA, thereby triggering an immune response. To unravel the exact mechanism by which TLR13 recognizes RNA15, we combined experimental approaches with molecular dynamics simulations. Our results suggest that RNA15 adopts a stable hairpin structure in solution, protected from nuclease degradation by intramolecular interactions. TLR13 specifically recognizes this hairpin structure, leading to the dimerization of TLR13. This interaction further induces RNA15 to transition into a stem-loop-like conformation, thereby activating TLR13 downstream signaling. Additionally, our study indicates that TLR13 can form stable dimers in the membrane independently of ligand binding. Although the hairpin structure is the predominant form of RNA15 in solution, the temporary stem-loop-like structure can spontaneously bind to dimeric TLR13, initiating the immune response. These insights deepen our understanding of the complex recognition process of RNA15 by TLR13 and explore the complicated mechanisms governing innate immune system function.

The role of ribosomal protein networks in ribosome dynamics

Accurate protein synthesis requires ribosomes to integrate signals from distant functional sites and execute complex dynamics. Despite advances in understanding ribosome structure and function, two key questions remain: how information is transmitted between these distant sites, and how ribosomal movements are synchronized? We recently highlighted the existence of ribosomal protein networks, likely evolved to participate in ribosome signaling. Here, we investigate the relationship between ribosomal protein networks and ribosome dynamics. Our findings show that major motion centers in the bacterial ribosome interact specifically with r-proteins, and that ribosomal RNA exhibits high mobility around each r-protein. This suggests that periodic electrostatic changes in the context of negatively charged residues (Glu and Asp) induce RNA-protein 'distance-approach' cycles, controlling key ribosomal movements during translocation. These charged residues play a critical role in modulating electrostatic repulsion between RNA and proteins, thus coordinating ribosomal dynamics. We propose that r-protein networks synchronize ribosomal dynamics through an 'electrostatic domino' effect, extending the concept of allostery to the regulation of movements within supramolecular assemblies.

GPCR Signaling: A Study of the Interplay Between Structure, Energy, and Function

G protein-coupled receptors (GPCRs) exemplify sophisticated allosteric communication, transducing extracellular signals through ligand-induced structural rearrangements that resonate through the molecular scaffold. Despite extensive study, the biophysical underpinnings of how conformational changes spread remain unclear. This work employs a novel physics-based framework to characterize the role of energy dissipation in directing intramolecular signaling pathways. By modeling each residue as a network of coupled oscillators, we generate a localization landscape depicting the vibrational energy distribution throughout the protein scaffold. Quantifying directional energy flux between residues reveals distinct pathways for energy and information transfer, illuminating sequences of allosteric communication. Our analysis of CB1 and CCR5 crystal structures unveils an anisotropic pattern of energy dissipation aligning with key functional dynamics, such as activation-related conformational changes. These anisotropic patterns of vibrational energy flow constitute pre-configured channels for allosteric signaling. Elucidating the relationship between structural topology and energy dissipation patterns provides key insights into the thermodynamic drivers of conformational signaling. This methodology significantly advances our mechanistic understanding of allostery in GPCRs and presents a broadly applicable approach for rationally dissecting allosteric communication pathways, with potential implications for structure-based drug design targeting these critical receptors.

Distant Ribose 2'-O-Methylation of 23S rRNA Helix 69 Pre-Orders the Capreomycin Drug Binding Pocket at the Ribosome Subunit Interface

Loss of ribosomal RNA (rRNA) modifications incorporated by the intrinsic methyltransferase TlyA results in reduced sensitivity to tuberactinomycin antibiotics such as capreomycin. However, the mechanism by which rRNA methylation alters drug binding, particularly at the distant but functionally more important site in 23S rRNA Helix 69 (H69), is currently unknown. We determined high-resolution cryo-electron microscopy structures of the Mycolicibacterium smegmatis 70S ribosome with or without the two ribose 2'-O-methyl modifications incorporated by TlyA. In the unmodified ribosome, the tip of H69 adopts a more compact conformation, positioning two key nucleotides (A2137 and C2138) such that interactions with capreomycin would be lost and the binding pocket partially occluded. In contrast, methylation of 23S rRNA nucleotide C2144 results in conformational changes that propagate from the site of modification to the H69 tip, resulting in its movement away from h44, a more favorable positioning of C2138 and adoption of a more open conformation to enable capreomycin binding. Methylation of h44 also results in structural rearrangements at the H69-h44 interface that further support antibiotic binding. These structures thus reveal the effect and regulation of distant rRNA methylation on ribosome-targeting antibiotic binding.

Ribosome Structural Changes Dynamically Affect Ribosome Function

Ribosomes were known to be multicomponent complexes as early as the 1960s. Nonetheless, the prevailing view for decades considered active ribosomes to be a monolithic population, in which all ribosomes are identical in composition and function. This implied that ribosomes themselves did not actively contribute to the regulation of protein synthesis. In this perspective, I review evidence for a different model, based on results showing that ribosomes can harbor different types of ribosomal RNA (rRNA) and ribosomal proteins (r-proteins) and, furthermore, need not contain a complete set of r-proteins. I also summarize recent results favoring the notion that such distinct types of ribosomes have different affinities for specific messenger RNAs and may execute the translation process differently. Thus, ribosomes should be considered active contributors to the regulation of protein synthesis.

A parameterized two-domain thermodynamic model explains diverse mutational effects on protein allostery

New experimental findings continue to challenge our understanding of protein allostery. Recent deep mutational scanning study showed that allosteric hotspots in the tetracycline repressor (TetR) and its homologous transcriptional factors are broadly distributed rather than spanning well-defined structural pathways as often assumed. Moreover, hotspot mutation-induced allostery loss was rescued by distributed additional mutations in a degenerate fashion. Here, we develop a two-domain thermodynamic model for TetR, which readily rationalizes these intriguing observations. The model accurately captures the in vivo activities of various mutants with changes in physically transparent parameters, allowing the data-based quantification of mutational effects using statistical inference. Our analysis reveals the intrinsic connection of intra- and inter-domain properties for allosteric regulation and illustrate epistatic interactions that are consistent with structural features of the protein. The insights gained from this study into the nature of two-domain allostery are expected to have broader implications for other multi-domain allosteric proteins.

A parametrized two-domain thermodynamic model explains diverse mutational effects on protein allostery

New experimental findings continue to challenge our understanding of protein allostery. Recent deep mutational scanning study showed that allosteric hotspots in the tetracycline repressor (TetR) and its homologous transcriptional factors are broadly distributed rather than spanning well-defined structural pathways as often assumed. Moreover, hotspot mutation-induced allostery loss was rescued by distributed additional mutations in a degenerate fashion. Here, we develop a two-domain thermodynamic model for TetR, which readily rationalizes these intriguing observations. The model accurately captures the in vivo activities of various mutants with changes in physically transparent parameters, allowing the data-based quantification of mutational effects using statistical inference. Our analysis reveals the intrinsic connection of intra- and inter-domain properties for allosteric regulation and illustrate epistatic interactions that are consistent with structural features of the protein. The insights gained from this study into the nature of two-domain allostery are expected to have broader implications for other multidomain allosteric proteins.

Modulation of Allostery with Multiple Mechanisms by Hotspot Mutations in TetR

Modulating allosteric coupling offers unique opportunities for biomedical applications. Such efforts can benefit from efficient prediction and evaluation of allostery hotspot residues that dictate the degree of cooperativity between distant sites. We demonstrate that effects of allostery hotspot mutations can be evaluated qualitatively and semiquantitatively by molecular dynamics simulations in a bacterial tetracycline repressor (TetR). The simulations recapitulate the effects of these mutations on abolishing the induction function of TetR and provide a rationale for the different rescuabilities observed to restore allosteric coupling of the hotspot mutations. We demonstrate that the same noninducible phenotype could be the result of perturbations in distinct structural and energetic properties of TetR. Our work underscores the value of explicitly computing the functional free energy landscapes to effectively evaluate and rank hotspot mutations despite the prevalence of compensatory interactions and therefore provides quantitative guidance to allostery modulation for therapeutic and engineering applications.

Modulation of Allostery with Multiple Mechanisms by Hotspot Mutations in TetR

Modulating allosteric coupling offers unique opportunities for biomedical applications. Such efforts can benefit from efficient prediction and evaluation of allostery hotspot residues that dictate the degree of co-operativity between distant sites. We demonstrate that effects of allostery hotspot mutations can be evaluated qualitatively and semi-quantitatively by molecular dynamics simulations in a bacterial tetracycline repressor (TetR). The simulations recapitulate the effects of these mutations on abolishing the induction function of TetR and provide a rationale for the different degrees of rescuability observed to restore allosteric coupling of the hotspot mutations. We demonstrate that the same non-inducible phenotype could be the result of perturbations in distinct structural and energetic properties of TetR. Our work underscore the value of explicitly computing the functional free energy landscapes to effectively evaluate and rank hotspot mutations despite the prevalence of compensatory interactions, and therefore provide quantitative guidance to allostery modulation for therapeutic and engineering applications.

Aminobenzoic Acid Derivatives Obstruct Induced Fit in the Catalytic Center of the Ribosome

The Escherichia coli (E. coli) ribosome can incorporate a variety of non-l-α-amino acid monomers into polypeptide chains in vitro but with poor efficiency. Although these monomers span a diverse set of compounds, there exists no high-resolution structural information regarding their positioning within the catalytic center of the ribosome, the peptidyl transferase center (PTC). Thus, details regarding the mechanism of amide bond formation and the structural basis for differences and defects in incorporation efficiency remain unknown. Within a set of three aminobenzoic acid derivatives-3-aminopyridine-4-carboxylic acid (Apy), ortho-aminobenzoic acid (oABZ), and meta-aminobenzoic acid (mABZ)-the ribosome incorporates Apy into polypeptide chains with the highest efficiency, followed by oABZ and then mABZ, a trend that does not track with the nucleophilicity of the reactive amines. Here, we report high-resolution cryo-EM structures of the ribosome with each of these three aminobenzoic acid derivatives charged on tRNA bound in the aminoacyl-tRNA site (A-site). The structures reveal how the aromatic ring of each monomer sterically blocks the positioning of nucleotide U2506, thereby preventing rearrangement of nucleotide U2585 and the resulting induced fit in the PTC required for efficient amide bond formation. They also reveal disruptions to the bound water network that is believed to facilitate formation and breakdown of the tetrahedral intermediate. Together, the cryo-EM structures reported here provide a mechanistic rationale for differences in reactivity of aminobenzoic acid derivatives relative to l-α-amino acids and each other and identify stereochemical constraints on the size and geometry of non-monomers that can be accepted efficiently by wild-type ribosomes.

Multiple Non-Canonical Base-Stacking Interactions as One of the Major Determinants of RNA Tertiary Structure Organization

Stacking interactions of heterocyclic bases of ribonucleotides are one of the most important factors in the organization of RNA secondary and tertiary structure. Most of these (canonical) interactions are formed between adjacent residues in RNA polynucleotide chains. However, with the accumulation of data on the atomic tertiary structures of various RNAs and their complexes with proteins, it has become clear that nucleotide residues that are not adjacent in the polynucleotide chains and are sometimes separated in the RNA primary structure by tens or hundreds of nucleotides can interact via (non-canonical) base stacking. This paper presents an exhaustive database of such nonadjacent base-stacking elements (NA-BSEs) and their environment in the macromolecules of natural and synthetic RNAs. Analysis of these data showed that NA-BSE-forming nucleotides, on average, account for about a quarter of all nucleotides in a particular RNA and, therefore, should be considered as bona fide motifs of the RNA tertiary structure. We also classified NA-BSEs by their location in RNA macromolecules. It was shown that the structure-forming role of NA-BSEs involves compact folding of single-stranded RNA loops, transformation of double-stranded bulges into imperfect helices, and binding of RNA regions distant in the primary and secondary RNA structure.

Undersampling and the inference of coevolution in proteins

Protein structure, function, and evolution depend on local and collective epistatic interactions between amino acids. A powerful approach to defining these interactions is to construct models of couplings between amino acids that reproduce the empirical statistics (frequencies and correlations) observed in sequences comprising a protein family. The top couplings are then interpreted. Here, we show that as currently implemented, this inference unequally represents epistatic interactions, a problem that fundamentally arises from limited sampling of sequences in the context of distinct scales at which epistasis occurs in proteins. We show that these issues explain the ability of current approaches to predict tertiary contacts between amino acids and the inability to obviously expose larger networks of functionally relevant, collectively evolving residues called sectors. This work provides a necessary foundation for more deeply understanding and improving evolution-based models of proteins.

Binding of the peptide deformylase on the ribosome surface modulates the exit tunnel interior

Proteosynthesis on ribosomes is regulated at many levels. Conformational changes of the ribosome, possibly induced by external factors, may transfer over large distances and contribute to the regulation. The molecular principles of this long-distance allostery within the ribosome remain poorly understood. Here, we use structural analysis and atomistic molecular dynamics simulations to investigate peptide deformylase (PDF), an enzyme that binds to the ribosome surface near the ribosomal protein uL22 during translation and chemically modifies the emerging nascent peptide. Our simulations of the entire ribosome-PDF complex reveal that the PDF undergoes a swaying motion on the ribosome surface at the submicrosecond timescale. We show that the PDF affects the conformational dynamics of parts of the ribosome over distances of more than 5 nm. Using a supervised-learning algorithm, we demonstrate that the exit tunnel is influenced by the presence or absence of PDF. Our findings suggest a possible effect of the PDF on the nascent peptide translocation through the ribosome exit tunnel.

Biological innovation in the functional landscape of a model regulator, or the lactose operon repressor

The operon model was proposed six decades ago. And yet, despite all this time, the lactose operon repressor, LacI, remains a subject of major interest. While it is well established that LacI can exist in two functional forms, one that renders the operon inactive via binding of LacI to DNA and another, bound to an inducer that does not allow repression, how it switches from one to the other is still not well understood. The construction of a library of several tens of thousands of LacI mutants has revealed some unexpected features. In particular, the transition implemented in some of them reveals a new type of transcription regulation: band-pass (OFF/ON/OFF) and band-stop (ON/OFF/ON) filters. This makes it natural to think that it is the network of hydrogen bonds associated with the water bound to the molecule that allows the remote interconnection between the binding site to an inducer molecule and the one that binds it to the DNA.

Mechanistic studies of non-canonical amino acid mutagenesis

Over the past decade, harnessing the cellular protein synthesis machinery to incorporate non-canonical amino acids (ncAAs) into tailor-made peptides has significantly advanced many aspects of molecular science. More recently, groundbreaking progress in our ability to engineer this machinery for improved ncAA incorporation has led to significant enhancements of this powerful tool for biology and chemistry. By revealing the molecular basis for the poor or improved incorporation of ncAAs, mechanistic studies of ncAA incorporation by the protein synthesis machinery have tremendous potential for informing and directing such engineering efforts. In this chapter, we describe a set of complementary biochemical and single-molecule fluorescence assays that we have adapted for mechanistic studies of ncAA incorporation. Collectively, these assays provide data that can guide engineering of the protein synthesis machinery to expand the range of ncAAs that can be incorporated into peptides and increase the efficiency with which they can be incorporated, thereby enabling the full potential of ncAA mutagenesis technology to be realized.

Engineering molecular translation systems

Molecular translation systems provide a genetically encoded framework for protein synthesis, which is essential for all life. Engineering these systems to incorporate non-canonical amino acids (ncAAs) into peptides and proteins has opened many exciting opportunities in chemical and synthetic biology. Here, we review recent advances that are transforming our ability to engineer molecular translation systems. In cell-based systems, new processes to synthesize recoded genomes, tether ribosomal subunits, and engineer orthogonality with high-throughput workflows have emerged. In cell-free systems, adoption of flexizyme technology and cell-free ribosome synthesis and evolution platforms are expanding the limits of chemistry at the ribosome's RNA-based active site. Looking forward, innovations will deepen understanding of molecular translation and provide a path to polymers with previously unimaginable structures and functions.