Long-range Tertiary Interactions Research Articles

Loop-Sequence Features and Stability Determinants in Antibody Variable Domains by High-Throughput Experiments

Protein loops are frequently considered as critical determinants in protein structure and function. Recent advances in high-throughput methods for DNA sequencing and thermal stability measurement have enabled effective exploration of sequence-structure-function relationships in local protein regions. Using these data-intensive technologies, we investigated the sequence-structure-function relationships of six complementarity-determining regions (CDRs) and ten non-CDR loops in the variable domains of a model vascular endothelial growth factor (VEGF)-binding single-chain antibody variable fragment (scFv) whose sequence had been optimized via a consensus-sequence approach. The results show that only a handful of residues involving long-range tertiary interactions distant from the antigen-binding site are strongly coupled with antigen binding. This implies that the loops are passive regions in protein folding; the essential sequences of these regions are dictated by conserved tertiary interactions and the consensus local loop-sequence features contribute little to protein stability and function.

An RNA folding motif: GNRA tetraloop–receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Average Conformations Determined from PRE Data Provide High-Resolution Maps of Transient Tertiary Interactions in Disordered Proteins

In the last decade it has become evident that disordered states of proteins play important physiological and pathological roles and that the transient tertiary interactions often present in these systems can play a role in their biological activity. The structural characterization of such states has so far largely relied on ensemble representations, which in principle account for both their local and global structural features. However, these approaches are inherently of low resolution due to the large number of degrees of freedom of conformational ensembles and to the sparse nature of the experimental data used to determine them. Here, we overcome these limitations by showing that tertiary interactions in disordered states can be mapped at high resolution by fitting paramagnetic relaxation enhancement data to a small number of conformations, which can be as low as one. This result opens up the possibility of determining the topology of cooperatively collapsed and hidden folded states when these are present in the vast conformational landscape accessible to disordered states of proteins. As a first application, we study the long-range tertiary interactions of acid-unfolded apomyoglobin from experimentally measured paramagnetic relaxation enhancement data.

Probing Retroviral and Retrotransposon Genome Structures: The “SHAPE” of Things to Come

Understanding the nuances of RNA structure as they pertain to biological function remains a formidable challenge for retrovirus research and development of RNA-based therapeutics, an area of particular importance with respect to combating HIV infection. Although a variety of chemical and enzymatic RNA probing techniques have been successfully employed for more than 30 years, they primarily interrogate small (100–500 nt) RNAs that have been removed from their biological context, potentially eliminating long-range tertiary interactions (such as kissing loops and pseudoknots) that may play a critical regulatory role. Selective 2′ hydroxyl acylation analyzed by primer extension (SHAPE), pioneered recently by Merino and colleagues, represents a facile, user-friendly technology capable of interrogating RNA structure with a single reagent and, combined with automated capillary electrophoresis, can analyze an entire 10,000-nucleotide RNA genome in a matter of weeks. Despite these obvious advantages, SHAPE essentially provides a nucleotide “connectivity map,” conversion of which into a 3-D structure requires a variety of complementary approaches. This paper summarizes contributions from SHAPE towards our understanding of the structure of retroviral genomes, modifications to which technology that have been developed to address some of its limitations, and future challenges.

Transient Structure and SH3 Interaction Sites in an Intrinsically Disordered Fragment of the Hepatitis C Virus Protein NS5A

Understanding the molecular mechanisms involved in virus replication and particle assembly is of primary fundamental and biomedical importance. Intrinsic conformational disorder plays a prominent role in viral proteins and their interaction with other viral and host cell proteins via transiently populated structural elements. Here, we report on the results of an investigation of an intrinsically disordered 188-residue fragment of the hepatitis C virus non-structural protein 5A (NS5A), which contains a classical poly-proline Src homology 3 (SH3) binding motif, using sensitivity- and resolution-optimized multidimensional NMR methods, complemented by small-angle X-ray scattering data. Our study provides detailed atomic-resolution information on transient local and long-range structure, as well as fast time scale dynamics in this NS5A fragment. In addition, we could characterize two distinct interaction modes with the SH3 domain of Bin1 (bridging integrator protein 1), a pro-apoptotic tumor suppressor. Despite being largely disordered, the protein contains three regions that transiently adopt α-helical structures, partly stabilized by long-range tertiary interactions. Two of these transient α-helices form a noncanonical SH3-binding motif, which allows low-affinity SH3 binding. Our results contribute to a better understanding of the role of the NS5A protein during hepatitis C virus infection. The present work also highlights the power of NMR spectroscopy to characterize multiple binding events including short-lived transient interactions between globular and highly disordered proteins.

Structural architecture of an RNA that competitively inhibits RNase L

Activation of RNase L endonuclease activity is part of the mammalian innate immune response to viral infection. The poliovirus RNA genome contains a sequence in its protein-coding region that can act as a competitive inhibitor of RNase L. Mutation, sequence, and functional analysis of this competitive inhibitor RNA (ciRNA) revealed that its activity depends on specific sequences, showed that a loop-loop hairpin interaction forms in the ciRNA, and suggested the presence of a loop E motif. These features lead to the hypothesis that the ciRNA's function is conferred in part by a specific three-dimensional folded RNA architecture. By using a combination of biophysical, mutational, and functional studies, we have mapped features of the three-dimensional architecture of the ciRNA in its unbound form. We show that the loop-loop interaction forms in the free ciRNA and affects the overall structure, perhaps forming long-range tertiary interactions with the loop E motif. Local tight RNA-RNA backbone packing occurs in parts of the structure, but the fold appears to be less stable than many other tightly packed RNAs. This feature may allow the ciRNA to accommodate the translocation of ribosomes and polymerase across this multifunctional region of the viral RNA but also to function as an RNase L inhibitor.

Long-range pseudoknot interactions dictate the regulatory response in the tetrahydrofolate riboswitch

Tetrahydrofolate (THF), a biologically active form of the vitamin folate (B(9)), is an essential cofactor in one-carbon transfer reactions. In bacteria, expression of folate-related genes is controlled by feedback modulation in response to specific binding of THF and related compounds to a riboswitch. Here, we present the X-ray structures of the THF-sensing domain from the Eubacterium siraeum riboswitch in the ligand-bound and unbound states. The structure reveals an "inverted" three-way junctional architecture, most unusual for riboswitches, with the junction located far from the regulatory helix P1 and not directly participating in helix P1 formation. Instead, the three-way junction, stabilized by binding to the ligand, aligns the riboswitch stems for long-range tertiary pseudoknot interactions that contribute to the organization of helix P1 and therefore stipulate the regulatory response of the riboswitch. The pterin moiety of the ligand docks in a semiopen pocket adjacent to the junction, where it forms specific hydrogen bonds with two moderately conserved pyrimidines. The aminobenzoate moiety stacks on a guanine base, whereas the glutamate moiety does not appear to make strong interactions with the RNA. In contrast to other riboswitches, these findings demonstrate that the THF riboswitch uses a limited number of available determinants for ligand recognition. Given that modern antibiotics target folate metabolism, the THF riboswitch structure provides insights on mechanistic aspects of riboswitch function and may help in manipulating THF levels in pathogenic bacteria.

The Role of Electrostatic Relaxation on the Folding Kinetics of a Bacterial Ribozyme

The self-assembly of catalytic RNA ribozymes into compact, native structures is critical for their functions in the cell. The first step in forming RNA tertiary structure is the neutralization by cations of the negative charges of the phosphates. This electrostatic stabilization enables dynamical exploration of more compact conformations, and the search for long-range tertiary interactions. Our previous time-resolved Small Angle X-ray Scattering (SAXS) studies showed that the Azoarcus ribozyme exhibits triphasic folding kinetics: up to 90% folds in less than 10 ms, which we attribute to specific collapse.

Folding of the lysine riboswitch: importance of peripheral elements for transcriptional regulation

The Bacillus subtilis lysC lysine riboswitch modulates its own gene expression upon lysine binding through a transcription attenuation mechanism. The riboswitch aptamer is organized around a single five-way junction that provides the scaffold for two long-range tertiary interactions (loop L2–loop L3 and helix P2–loop L4)—all of this for the creation of a specific lysine binding site. We have determined that the interaction P2–L4 is particularly important for the organization of the ligand-binding site and for the riboswitch transcription attenuation control. Moreover, we have observed that a folding synergy between L2–L3 and P2–L4 allows both interactions to fold at lower magnesium ion concentrations. The P2–L4 interaction is also critical for the close juxtaposition involving stems P1 and P5. This is facilitated by the presence of lysine, suggesting an active role of the ligand in the folding transition. We also show that a previously uncharacterized stem–loop located in the expression platform is highly important for the riboswitch activity. Thus, folding elements located in the aptamer and the expression platform both influence the lysine riboswitch gene regulation.

Long-range tertiary interactions in single hammerhead ribozymes bias motional sampling toward catalytically active conformations

Enzymes generally are thought to derive their functional activity from conformational motions. The limited chemical variation in RNA suggests that such structural dynamics may play a particularly important role in RNA function. Minimal hammerhead ribozymes are known to cleave efficiently only in ∼ 10-fold higher than physiologic concentrations of Mg(2+) ions. Extended versions containing native loop-loop interactions, however, show greatly enhanced catalytic activity at physiologically relevant Mg(2+) concentrations, for reasons that are still ill-understood. Here, we use Mg(2+) titrations, activity assays, ensemble, and single molecule fluorescence resonance energy transfer (FRET) approaches, combined with molecular dynamics (MD) simulations, to ask what influence the spatially distant tertiary loop-loop interactions of an extended hammerhead ribozyme have on its structural dynamics. By comparing hammerhead variants with wild-type, partially disrupted, and fully disrupted loop-loop interaction sequences we find that the tertiary interactions lead to a dynamic motional sampling that increasingly populates catalytically active conformations. At the global level the wild-type tertiary interactions lead to more frequent, if transient, encounters of the loop-carrying stems, whereas at the local level they lead to an enrichment in favorable in-line attack angles at the cleavage site. These results invoke a linkage between RNA structural dynamics and function and suggest that loop-loop interactions in extended hammerhead ribozymes-and Mg(2+) ions that bind to minimal ribozymes-may generally allow more frequent access to a catalytically relevant conformation(s), rather than simply locking the ribozyme into a single active state.

Nonhierarchical Ribonucleoprotein Assembly Suggests a Strain-Propagation Model for Protein-Facilitated RNA Folding

Proteins play diverse and critical roles in cellular ribonucleoproteins (RNPs) including promoting formation of and stabilizing active RNA conformations. Yet, the conformational changes required to convert large RNAs into active RNPs have proven difficult to characterize fully. Here we use high-resolution approaches to monitor both local nucleotide flexibility and solvent accessibility for nearly all nucleotides in the bI3 group I intron RNP in four assembly states: the free RNA, maturase-bound RNA, Mrs1-bound RNA, and the complete six-component holocomplex. The free RNA is misfolded relative to the secondary structure required for splicing. The maturase and Mrs1 proteins each stabilized long-range tertiary interactions, but neither protein alone induced folding into the functional secondary structure. In contrast, simultaneous binding by both proteins results in large secondary structure rearrangements in the RNA and yielded the catalytically active group I intron structure. Secondary and tertiary folding of the RNA component of the bI3 RNP are thus not independent: RNA folding is strongly nonhierarchical. These results emphasize that protein-mediated stabilization of RNA tertiary interactions functions to pull the secondary structure into an energetically disfavored, but functional, conformation and emphasize a new role for facilitator proteins in RNP assembly.

NMR characterization of hydrophobic collapses in amyloidogenic unfolded states and their implications for amyloid formation

NMR spectroscopy was used to characterize hydrophobic clusters in amyloidogenic unfolded states of a protein and their implications for amyloid formation. Three local hydrophobic clusters were observed in the amyloidogenic state of the phosphatidylinositol 3-kinase (PI3K) SH3 domain. Our NMR studies showed that residues with high average area buried upon folding (AABUF) parameter collapsed to form the clusters. Interestingly, the hydrophobic collapses were not stabilized by long-range tertiary interactions among the clusters that were typically observed in non-amyloidogenic unfolded states of various proteins. The lack of the long-range interactions may be a critical property of the amyloidogenic unfolded state. The SH3 domain was also engineered to disrupt one of the clusters by a single-point mutagenesis (W55G), which allowed us to investigate the effect of the clustering on folding and misfolding. The mutant form of the SH3 domain was not able to fold under folding conditions of the wild type protein (pH 3.6–4.0), supporting the cooperative folding hypothesis. However, aggregation properties of the mutant form were not influenced by the mutation, suggesting the SH3 domain forms amyloid via a non-cooperative process.

Dissecting the influence of Mg 2+ on 3D architecture and ligand-binding of the guanine-sensing riboswitch aptamer domain

Long-range tertiary interactions determine the three-dimensional structure of a number of metabolite-binding riboswitch RNA elements and were found to be important for their regulatory function. For the guanine-sensing riboswitch of the Bacillus subtilis xpt-pbuX operon, our previous NMR-spectroscopic studies indicated pre-formation of long-range tertiary contacts in the ligand-free state of its aptamer domain. Loss of the structural pre-organization in a mutant of this RNA (G37A/C61U) resulted in the requirement of Mg2+ for ligand binding. Here, we investigate structural and stability aspects of the wild-type aptamer domain (Gsw) and the G37A/C61U-mutant (Gswloop) of the guanine-sensing riboswitch and their Mg2+-induced folding characteristics to dissect the role of long-range tertiary interactions, the link between pre-formation of structural elements and ligand-binding properties and the functional stability. Destabilization of the long-range interactions as a result of the introduced mutations for Gswloop or the increase in temperature for both Gsw and Gswloop involves pronounced alterations of the conformational ensemble characteristics of the ligand-free state of the riboswitch. The increased flexibility of the conformational ensemble can, however, be compensated by Mg2+. We propose that reduction of conformational dynamics in remote regions of the riboswitch aptamer domain is the minimal pre-requisite to pre-organize the core region for specific ligand binding.

The Mrs1 Splicing Factor Binds the bI3 Group I Intron at Each of Two Tetraloop-Receptor Motifs

Most large ribozymes require protein cofactors in order to function efficiently. The yeast mitochondrial bI3 group I intron requires two proteins for efficient splicing, Mrs1 and the bI3 maturase. Mrs1 has evolved from DNA junction resolvases to function as an RNA cofactor for at least two group I introns; however, the RNA binding site and the mechanism by which Mrs1 facilitates splicing were unknown. Here we use high-throughput RNA structure analysis to show that Mrs1 binds a ubiquitous RNA tertiary structure motif, the GNRA tetraloop-receptor interaction, at two sites in the bI3 RNA. Mrs1 also interacts at similar tetraloop-receptor elements, as well as other structures, in the self-folding Azoarcus group I intron and in the RNase P enzyme. Thus, Mrs1 recognizes general features found in the tetraloop-receptor motif. Identification of the two Mrs1 binding sites now makes it possible to create a model of the complete six-component bI3 ribonucleoprotein. All protein cofactors bind at the periphery of the RNA such that every long-range RNA tertiary interaction is stabilized by protein binding, involving either Mrs1 or the bI3 maturase. This work emphasizes the strong evolutionary pressure to bolster RNA tertiary structure with RNA-binding interactions as seen in the ribosome, spliceosome, and other large RNA machines.

The Role of Bulges and Hinges in RNA Folding

Long-range tertiary interactions govern RNA structural assembly, which is a critical step toward RNA biological functionality. Thereby, a universal strategy has emerged for global conformational change; flexible junctions enable unpaired nucleotides to act as beacons between helical regions. Bulges, for example, are versatile secondary structural elements implicated in helix recognition and packing. The P4-P6 domain of the Tetrahymena ribozyme utilizes this folding strategy by hinging to form two inter-helical tertiary contacts, the adenine-rich (A-rich) bulge and the tetraloop-tetraloop receptor interactions. To explore the kinetic and thermodynamic properties of tertiary contact formation, we probe the P4-P6 domain hinging and ribose zippering that forms the A-rich bulge interaction using single-molecule FRET methods. We obtain the docking and undocking rates of the A-rich bulge and P4 helix as function of cation concentration and temperature. Docking is accelerated and undocking decelerated by Mg2+. In spite of rapid docking at high [Mg2+] (kdock = 20 ± 2 s−1), the A-rich bulge interaction is only marginally stable (Kdock = 1.2 ± 0.1). These results support that the role of the A-rich bulge is to kinetically direct P4-P6 domain folding while thermodynamic stability is added through the tetraloop-receptor interaction. Formation of the A-rich bulge contact shows specificity for divalent cations, with a preference for Mg2+ as anticipated from the Mg2+ coordination observed in structural data. A significant kinetic heterogeneity is characterized; only 50% of the molecules exhibit efficient folding at high [Mg2+]. Mutations of the A-rich bulge construct reveal a crucial role of the P4-P6 secondary architecture in enabling the A-rich bulge contact.

Structural insights into amino acid binding and gene control by a lysine riboswitch

In bacteria, the intracellular concentration of several amino acids is controlled by riboswitches. One of the important regulatory circuits involves lysine-specific riboswitches, which direct the biosynthesis and transport of lysine and precursors common for lysine and other amino acids. To understand the molecular basis of amino acid recognition by riboswitches, here we present the crystal structure of the 174-nucleotide sensing domain of the Thermotoga maritima lysine riboswitch in the lysine-bound (1.9 ångström (A)) and free (3.1 A) states. The riboswitch features an unusual and intricate architecture, involving three-helical and two-helical bundles connected by a compact five-helical junction and stabilized by various long-range tertiary interactions. Lysine interacts with the junctional core of the riboswitch and is specifically recognized through shape-complementarity within the elongated binding pocket and through several direct and K(+)-mediated hydrogen bonds to its charged ends. Our structural and biochemical studies indicate preformation of the riboswitch scaffold and identify conformational changes associated with the formation of a stable lysine-bound state, which prevents alternative folding of the riboswitch and facilitates formation of downstream regulatory elements. We have also determined several structures of the riboswitch bound to different lysine analogues, including antibiotics, in an effort to understand the ligand-binding capabilities of the lysine riboswitch and understand the nature of antibiotic resistance. Our results provide insights into a mechanism of lysine-riboswitch-dependent gene control at the molecular level, thereby contributing to continuing efforts at exploration of the pharmaceutical and biotechnological potential of riboswitches.

MS3D structural elucidation of the HIV-1 packaging signal

The structure of HIV-1 Psi-RNA has been elucidated by a concerted approach combining structural probes with mass spectrometric detection (MS3D), which is not affected by the size and crystallization properties of target biomolecules. Distance constraints from bifunctional cross-linkers provided the information required for assembling an all-atom model from the high-resolution coordinates of separate domains by triangulating their reciprocal placement in 3D space. The resulting structure revealed a compact cloverleaf morphology stabilized by a long-range tertiary interaction between the GNRA tetraloop of stemloop 4 (SL4) and the upper stem of stemloop 1 (SL1). The preservation of discrete stemloop structures ruled out the possibility that major rearrangements might produce a putative supersite with enhanced affinity for the nucleocapsid (NC) domain of the viral Gag polyprotein, which would drive genome recognition and packaging. The steric situation of single-stranded regions exposed on the cloverleaf structure offered a valid explanation for the stoichiometry exhibited by full-length Psi-RNA in the presence of NC. The participation of SL4 in a putative GNRA loop-receptor interaction provided further indications of the plasticity of this region of genomic RNA, which can also anneal with upstream sequences to stabilize alternative conformations of the 5' untranslated region (5'-UTR). Considering the ability to sustain specific NC binding, the multifaceted activities supported by the SL4 sequence suggest a mechanism by which Gag could actively participate in regulating the vital functions mediated by 5'-UTR. Substantiated by the 3D structure of Psi-RNA, the central role played by SL4 in specific RNA-RNA and protein-RNA interactions advances this domain as a primary target for possible therapeutic intervention.

Self-similarity of biopolymer backbones in the ribosome

Self-similar properties of the biopolymer backbones in the ribosome are investigated in terms of the fractal dimension. We especially estimate the chain fractal and capacity dimensions of the ribosomal RNAs and proteins, which are constituents of the ribosome. The fractal dimensions of both biopolymers are compared with that of the self-avoiding walk, which is a typical model of a polymer without interaction between monomers. We demonstrate that the fractality found in the ribosomal RNAs is pertinent to explain their structural characteristics: local helix formation and long-range tertiary interaction forming three-dimensional structures. The fractal dimension of the ribosomal protein supports the existence of the long and extended domain, which is hardly seen in the globular protein. The self-similarity also upholds the fact that the ribosomal proteins function primarily to stabilize the structure of the ribosome by both the long-extended domain of the protein penetrating into the inside of the RNA, and the globular domain interacting with the RNA on the exterior of it. These results partially, if not whole, unravel the structural characteristics of the biopolymers in the ribosome.

Structural basis of a ribozyme's thermostability: P1–L9 interdomain interaction in RNase P RNA

For stability, many catalytic RNAs rely on long-range tertiary interactions, the precise role of each often being unclear. Here we demonstrate that one of the three interdomain architectural struts of RNase P RNA (P RNA) is the key to activity at higher temperatures: disrupting the P1-L9 helix-tetraloop interaction in P RNA of the thermophile Thermus thermophilus decreased activity at high temperatures in the RNA-alone reaction and at low Mg2+ concentrations in the holoenzyme reaction. Conversely, implanting the P1-P9 module of T. thermophilus in the P RNA from the mesophile Escherichia coli converted the latter RNA into a thermostable one. Moreover, replacing the E. coli P1-P9 elements with a pseudoknot module that mediates the homologous interaction in Mycoplasma P RNAs not only conferred thermostability upon E. coli P RNA but also increased its maximum turnover rate at 55 degrees C to the highest yet described for a P RNA ribozyme.

How does the coupling of secondary and tertiary interactions control the folding of helical macromolecules?

The authors study how the simultaneous presence of short-range secondary and long-range tertiary interactions controls the folding and collapse behavior of a helical macromolecule. The secondary interactions stabilize the helical conformation of the chain, while the tertiary interactions govern its overall three-dimensional shape. The authors have carried out Monte Carlo simulations to study the effect of chain length on the folding and collapse behavior of the chain. They have calculated state diagrams for four chain lengths and found that the physics is very rich with a plethora of stable conformational states. In addition to the helix-coil and coil-globule transitions, their model describes the coupling between them which takes place at low temperatures. Under these conditions, their model predicts a cascade of continuous, conformational transitions between states with an increase in the strength of the tertiary interactions. During each transition the chain shrinks, i.e., collapses, in a rapid and specific manner. In addition, the number of the transitions increases with increasing chain length. They have also found that the low-temperature regions of the state diagram between the transition lines cannot be associated with specific structures of the chain, but rather, with ensembles of various configurations of the chain with similar characteristics. Based on these results the authors propose a mechanism for the folding and collapse of helical macromolecules which is further supported by the analysis of configurational, configurational, and thermodynamic properties of the chain.