Improving small-angle X-ray scattering data for structural analyses of the RNA world

Defining the shape, conformation, or assembly state of an RNA in solution often requires multiple investigative tools ranging from nucleotide analog interference mapping to X-ray crystallography. A key addition to this toolbox is small-angle X-ray scattering (SAXS). SAXS provides direct structural information regarding the size, shape, and flexibility of the particle in solution and has proven powerful for analyses of RNA structures with minimal requirements for sample concentration and volumes. In principle, SAXS can provide reliable data on small and large RNA molecules. In practice, SAXS investigations of RNA samples can show inconsistencies that suggest limitations in the SAXS experimental analyses or problems with the samples. Here, we show through investigations on the SAM-I riboswitch, the Group I intron P4-P6 domain, 30S ribosomal subunit from Sulfolobus solfataricus (30S), brome mosaic virus tRNA-like structure (BMV TLS), Thermotoga maritima asd lysine riboswitch, the recombinant tRNA^val, and yeast tRNA^phe that many problems with SAXS experiments on RNA samples derive from heterogeneity of the folded RNA. Furthermore, we propose and test a general approach to reducing these sample limitations for accurate SAXS analyses of RNA. Together our method and results show that SAXS with synchrotron radiation has great potential to provide accurate RNA shapes, conformations, and assembly states in solution that inform RNA biological functions in fundamental ways.     

Cooperative reorganization of a 420 subunit virus capsid

The complex protein capsids of many viruses exhibit dramatic reorganizations at critical stages in their life-cycle. Here, time-resolved solution X-ray scattering was used to study a dynamic, large-scale conformational maturation of the 420 subunit, 13 MDa, icosahedrally symmetric HK97 bacteriophage capsid. Isoscattering points in the time-resolved scattering patterns and singular value decomposition revealed that the expansion occurs as a cooperative, two-state reaction. The analysis demonstrates that the population shift from Prohead-II to Expansion Intermediate I, EI-I (60 A larger than Prohead-II) occurs in minutes, but does not reveal the time required for individual transitions that occur stochastically. Any intermediate forms that may be traversed during this conversion are unstable and do not constitute an appreciable population of the ensemble of particles. In an energetic landscape view, particles must undergo an energy barrier-crossing event in order to successfully convert from Prohead-II to EI-I. This implies that the particles "hop" over the energy barrier stochastically as they individually attain an expansion-active state. Interestingly, systematic deviations from single-exponential kinetics were observed for the population shift. This may indicate that in undergoing the irreversible conversion from Prohead-II to EI-I, particles are subject to a complex energy landscape that links the initial and final particle forms.     

Comparison of the structure of ribosomal 5S RNA from E. coli and from rat liver using X-ray scattering and dynamic light scattering

The structures of eukaryotic ribosomal 5S RNA from rat liver and of prokaryotic 5S RNA from E. coli (A-conformer) have been investigated by scattering methods. For both molecules, a molar mass of 44,500±4,000 was determined from small angle X-ray scattering as well as from dynamic light scattering. The shape parameters of the two rRNAs, volume V _c, surface O _c, radius of gyration R _s, maximum dimension of the molecule L, thickness D, and cross section radius of gyration R _sq, agree within the experimental error limits. The mean values are V _c=57±3 nm³, O _c=165±10 nm², R _s=3.37±0.05 nm, L=10.8±0.7 nm, D=1.57±0.07 nm, R _sa=0.92±0.01 nm.Identical structures for the E. coli 5S rRNA and the rat liver 5S rRNA at a resolution of 1 nm can be deduced from this agreement and from the comparison of experimental X-ray scattering curves and of experimental electron distance distribution function. The flat shape model derived for prokaryotic and eukaryotic 5S rRNA shows a compact region and two protruding arms. Double helical stems are eleven-fold helices with a mean base pair distance of 0.28 nm. Combining the shape information obtained from X-ray scattering with the information about the frictional behaviour of the molecules, deduced from the diffusion coefficients D _20,w ⁰ =(5.9±0.2)·10^-7 cm²s^-1 and (6.2±0.2)·10^-7 cm²s^-1 for rat liver 5S rRNA and E. coli 5S rRNA, respectively, a solvation shell of about 0.3 nm thickness around both molecules is determined. This structural similarity and the consensus secondary structure pattern derived from comparative sequence analyses suggest that all 5S rRNAs may indeed have conserved essentially the same type of folding of their polynucleotide strands during evolution, despite having very different sequences.     

The evidence of large-scale DNA-induced compaction in the mycobacterial chromosomal ParB

The bacterial chromosome trafficking apparatus or the segrosome participates in the mitotic-like segregation of the chromosomes prior to cell division in several bacteria. ParB, which is the parS DNA-binding component of the segrosome, polymerizes on the parS-adjacent chromosome to form a nucleoprotein filament of unknown nature for the segregation function. We combined static light scattering, circular dichroism and small-angle X-ray scattering to present evidence that the apo form of the mycobacterial ParB forms an elongated dimer with intrinsically disordered regions as well as folded domains in solution. A comparison of the solution scattering of the apo and the parS-bound ParBs indicates a rather drastic compaction of the protein upon DNA binding. We propose that this binding-induced conformational transition is priming the ParB for polymerization on the DNA template.     

Changes in the molecular packing of fibrillin microfibrils during extension indicate intrafibrillar and interfibrillar reorganization in elastic response

Fibrillin-rich microfibrils are the major structural components of the extracellular matrix that provide elasticity in a majority of connective tissues. The basis of elastic properties lies in the organization of fibrillin molecules, which, unfortunately, is still poorly understood. An X-ray diffraction study of hydrated fibrillin-rich microfibrils from zonular filaments has been conducted to give an insight into the molecular structure of microfibrils in intact tissue. A series of measurements was taken during controlled tissue extension to observe alterations in the lateral packing of microfibrils. Computer-generated simulated patterns were used to fit the experimental X-ray scattering data and to obtain the fibril diameter and lateral distance between the fibrils. The results suggest a nonlinear correlation between external strain and decrease in fibril diameter and lateral spacing. This was accompanied by a nonlinear increase in axial periodicity and a structure with a 160-nm periodicity, which is reported here for the first time using X-ray diffraction. These changes may reflect the unraveling of fibrillin from the complex folded arrangement into a linear structure. This finding supports a pleating model where fibrillin molecules are highly folded within the microfibrils; more importantly, the connection is made between the interaction of individual microfibrils and the change in their suprafibrillar coherent organization during extension. We suggest that the intermediate states observed in our study reflect sequential unfolding of fibrillin and can explain the process of its reversible unraveling.     

The Ligand-Free State of the TPP Riboswitch: A Partially Folded RNA Structure

Riboswitches are elements of mRNA that regulate gene expression by undergoing structural changes upon binding of small ligands. Although the structures of several riboswitches have been solved with their ligands bound, the ligand-free states of only a few riboswitches have been characterized. The ligand-free state is as important for the functionality of the riboswitch as the ligand-bound form, but the ligand-free state is often a partially folded structure of the RNA, with conformational heterogeneity that makes it particularly challenging to study. Here, we present models of the ligand-free state of a thiamine pyrophosphate riboswitch that are derived from a combination of complementary experimental and computational modeling approaches. We obtain a global picture of the molecule using small-angle X-ray scattering data and use an RNA structure modeling software, MC-Sym, to fit local structural details to these data on an atomic scale. We have used two different approaches to obtaining these models. Our first approach develops a model of the RNA from the structures of its constituent junction fragments in isolation. The second approach treats the RNA as a single entity, without bias from the structure of its individual constituents. We find that both approaches give similar models for the ligand-free form, but the ligand-bound models differ for the two approaches, and only the models from the second approach agree with the ligand-bound structure known previously from X-ray crystallography. Our models provide a picture of the conformational changes that may occur in the riboswitch upon binding of its ligand. Our results also demonstrate the power of combining experimental small-angle X-ray scattering data with theoretical structure prediction tools in the determination of RNA structures beyond riboswitches.     

Hinge stiffness is a barrier to RNA folding

Cation-mediated RNA folding from extended to compact, biologically active conformations relies on a temporal balance of forces. The Mg^2 +-mediated folding of the Tetrahymena thermophila ribozyme is characterized by rapid nonspecific collapse followed by tertiary-contact-induced compaction. This article focuses on an autonomously folding portion of the Tetrahymena ribozyme, its P4-P6 domain, in order to probe one facet of the rapid collapse: chain flexibility. The time evolution of P4-P6 folding was followed by global and local measures as a function of Mg^2 + concentration. While all concentrations of Mg^2 + studied are sufficient to screen the charge on the helices, the rates of compaction and tertiary contact formation diverge as the concentration of Mg^2 + increases; collapse is greatly accelerated by Mg^2 +, while tertiary contact formation is not. These studies highlight the importance of chain stiffness to RNA folding; at 10 mM Mg^2 +, a stiff hinge limits the rate of P4-P6 folding. At higher magnesium concentrations, the rate-limiting step shifts from hinge bending to tertiary contact formation.     

Structural transitions and thermodynamics of a glycine-dependent riboswitch from Vibrio cholerae

Riboswitches are complex folded RNA domains found in noncoding regions of mRNA that regulate gene expression upon small molecule binding. Recently, Breaker and coworkers reported a tandem aptamer riboswitch (VCI-II) that binds glycine cooperatively. Here, we use hydroxyl radical footprinting and small-angle X-ray scattering (SAXS) to study the conformations of this tandem aptamer as a function of Mg(2+) and glycine concentration. We fit a simple three-state thermodynamic model that describes the energetic coupling between magnesium-induced folding and glycine binding. Furthermore, we characterize the structural conformations of each of the three states: In low salt with no magnesium present, the VCI-II construct has an extended overall conformation, presumably representing unfolded structures. Addition of millimolar concentrations of Mg(2+) in the absence of glycine leads to a significant compaction and partial folding as judged by hydroxyl radical protections. In the presence of millimolar Mg(2+) concentrations, the tandem aptamer binds glycine cooperatively. The glycine binding transition involves a further compaction, additional tertiary packing interactions and further uptake of magnesium ions relative to the state in high Mg(2+) but no glycine. Employing density reconstruction algorithms, we obtain low resolution 3-D structures for all three states from the SAXS measurements. These data provide a first glimpse into the structural conformations of the VCI-II aptamer, establish rigorous constraints for further modeling, and provide a framework for future mechanistic studies.     

Comparison of the structures of the native form of rat liver 5S rRNA and yeast tRNAphe: small-angle and wide-angle X-ray scattering study

The small-angle and wide-angle X-ray scattering of tRNA^phe (yeast) and ribosomal 5S RNA (rat liver) in solution have been analysed and compared. tRNA^phe in solution is folded into a compact L-shaped structure similar to its structure in crystals. The geometry of the secondary structure of the double helical regions is also equivalent to the A-form in the crystalline state. Despite differences between the molar mosses of 5S rRNA (40 000 g mol^?1) and tRNA^phe (25 000 g mol^?1), and the fact that the 5S rRNA molecule is more anisometric than the tRNA^phe molecule, there are many structural similarities. The geometrical parameters of the secondary structure of double helical regions in both RNA molecules are almost identical; the mean rise per base pair is about 0.253–0.28 nm and the mean turn angle is about 32.5–33.5. Identical cross-sectional radii of gyration, R_sq,1 ≈ 1.16 nm and R_sq,2 = 0.92 nm, identical molar mass per unit length, , and a mean thickness of the molecules D ≈ 1.65 nm suggest a similar, nearly coplanar organization of isolated, double helical arms. Furthermore, there are compact regions in the central parts of both molecules, which are the sites of tertiary interactions in the tRNA^phe molecule and are a potential site of tertiary interactions in the SS rRNA molecule for stabilization of the complicated L-shape of the two molecules. Both molecules have a pseudo-twofold axis,w hich may play a role in recognition for binding of specific proteins.     

水、甲醇、乙醇对含水脑磷脂液晶结构影响的X射线衍射研究

用小角X射线散射法(SAXS)对含水脑磷脂分别用水、甲醇和乙醇实验所形成的液晶体系结构进行了研究。实验结果证明:在脑磷脂-水体系中,随着水含量增加,脑磷脂和水形成的双分子层液晶体系的层间距变大;在脑磷脂-甲醇体系中,随着甲醇含量增加,它们的层间距变小。在脑磷脂-乙醇体系中,随着乙醇含量的增加,它们的层间距先由小到大,继而又由大变小,然后液晶相逐渐消失,最后变成液态。水和甲醇、乙醇相比,水有使层间距变大的趋势,醇类有使层间距变小的趋势,随着醇中碳链的增长,层间距减小的趋势增大。     

SAXSDom: Modeling multidomain protein structures using small-angle X-ray scattering data

Many proteins are composed of several domains that pack together into a complex tertiary structure. Multidomain proteins can be challenging for protein structure modeling, particularly those for which templates can be found for individual domains but not for the entire sequence. In such cases, homology modeling can generate high quality models of the domains but not for the orientations between domains. Small-angle X-ray scattering (SAXS) reports the structural properties of entire proteins and has the potential for guiding homology modeling of multidomain proteins. In this article, we describe a novel multidomain protein assembly modeling method, SAXSDom that integrates experimental knowledge from SAXS with probabilistic Input-Output Hidden Markov model to assemble the structures of individual domains together. Four SAXS-based scoring functions were developed and tested, and the method was evaluated on multidomain proteins from two public datasets. Incorporation of SAXS information improved the accuracy of domain assembly for 40 out of 46 critical assessment of protein structure prediction multidomain protein targets and 45 out of 73 multidomain protein targets from the ab initio domain assembly dataset. The results demonstrate that SAXS data can provide useful information to improve the accuracy of domain-domain assembly. The source code and tool packages are available at https://github.com/jianlin-cheng/SAXSDom .     

Dissection of the ATP-Dependent Conformational Change Cycle of a Group II Chaperonin

Group II chaperonin captures an unfolded protein while in its open conformation and then mediates the folding of the protein during ATP-driven conformational change cycle. In this study, we performed kinetic analyses of the group II chaperonin from a hyperthermophilic archaeon, Thermococcus sp. KS-1 (TKS1-Cpn), by stopped-flow fluorometry and stopped-flow small-angle X-ray scattering to reveal the reaction cycle. Two TKS1-Cpn variants containing a Trp residue at position 265 or position 56 exhibit nearly the same fluorescence kinetics induced by rapid mixing with ATP. Fluorescence started to increase immediately after the start of mixing and reached a maximum at 1–2 s after mixing. Only in the presence of K⁺ that a gradual decrease in fluorescence was observed after the initial peak. Similar results were obtained by stopped-flow small-angle X-ray scattering. A rapid fluorescence increase, which reflects nucleotide binding, was observed for the mutant containing a Trp residue near the ATP binding site (K485W), irrespective of the presence or absence of K⁺. Without K⁺, a small, rapid fluorescence decrease followed the initial increase, and then a gradual decrease was observed. In contrast, with K⁺, a large, rapid fluorescence decrease occurred just after the initial increase, and then the fluorescence gradually increased. Finally, we observed ATP binding signal and also subtle conformational change in an ATPase-deficient mutant with K485W mutation. Based on these results, we propose a reaction cycle model for group II chaperonins.     

1964: The first model for the shape of a transfer RNA molecule. An account of an unpublished small-angle X-ray scattering study

The shape of non-fractionated Escherichia coli transfer RNA molecules in solution was investigated using small-angle X-ray scattering during the years 1960-1962 at the Centre de Recherche sur les Macromolécules in Strasbourg. The innermost region of the scattering curve yielded the average molecular weight (Mr) and the radius of gyration (Rg) of the particles, whereas the experimental data at large angles could be approximated at best by the scattering curve of a kinked rod-shaped molecule. The simplest model that was compatible with Mr, Rg, and the mass per unit length of the rod was a boomerang-shaped particle made of two double helical stems connected by a sharp kink. This model that eventually proved similar to the high-resolution L-shaped structure, was presented in my Ph.D. dissertation (J. Witz, Etude de la structure de quelques polynucléotides en solution par diffusion centrale des rayons X, Ph.D. dissertation, University of Strasbourg, France, 1964) but has never been published in detail. It is the purpose of this note to recall this story.     

Nanoscale uniformity of pore architecture in diatomaceous silica: a combined small and wide angle x-ray scattering study

Combined small and wide angle X‐ray scattering (SAXS and WAXS) analysis was applied to purified biogenic silica of cultured diatom frustules and of natural populations sampled on marine tidal flats. The overall WAXS patterns did not reveal crystalline phases (WAXS domain between 0.07 to 0.5 nm) in this biogenic silica, which is in line with previous reports on the amorphous character of the SiO₂ matrix of diatom frustules. One exception was the silica of the pennate species Cylindrotheca fusiformis Reimann et Lewin, which revealed wide peaks in the WAXS spectra. These peaks either indicate the presence of a yet unknown crystalline phase with a repetitive distance (d‐value ≈0.06 nm) or are caused by the ordering of the fibrous silica fragments; numerous girdle bands. The SAXS spectra revealed the size range of pores (diameter d between 3.0 and 65 nm), the presence of distinct pores (slope transitions), and structure factors (oscillation of the spectra). All slopes varied in the range of ?4.0 to ?2.5, with two clear common regions among species: d < 10 nm (slopes –4, denoted as region I and also called the Porod region), and 10.0 < d < 40.0 nm (slopes ?2.9 to ?3.8, denoted as region II). The existence of these common regions suggests the presence of comparable form (region I) and structure (region II) factors, respectively the shape of the primary building units of the silica and the geometry of the pores. Contrast variation experiments using dibromomethane to fill pores in the SiO₂ matrix showed that scattering was caused by pores rather than silica particles. Electron microscopic analysis confirmed the presence of circular, elliptical, and rectangular pores ranging in size from 3 to 65 nm, determining the structure factor. The fine architecture (length/width ratio of pore diameters) and distribution of the pores, however, seemed to be influenced by environmental factors, such as the salinity of and additions of AlCl₃ to the growth medium. The results indicate that diatoms deposit silica with pores <50 nm in size and are highly homologous with respect to geometry. Consequently, it is suggested that in diatoms, whether pennate or centric, the formation of silica at a nanoscale level is a uniform process.     

Characterization of protein fold by wide-angle X-ray solution scattering

Wide-angle X-ray solution scattering (WAXS) patterns contain substantial information about the three-dimensional structure of a protein. Although WAXS data have far less information than is required for determination of a full three-dimensional structure, the actual amount of information contained in a WAXS pattern has not been carefully quantified. Here we carry out an analysis of the amount of information that can be extracted from a WAXS pattern and demonstrate that it is adequate to estimate the secondary-structure content of a protein and to strongly limit its possible tertiary structures. WAXS patterns computed from the atomic coordinates of a set of 498 protein domains representing all of known fold space were used as the basis for constructing a multidimensional space of all corresponding WAXS patterns (‘WAXS space’). Within WAXS space, each scattering pattern is represented by a single vector. A principal components analysis was carried out to identify those directions in WAXS space that provide the greatest discrimination among patterns. The number of dimensions that provide significant discrimination among protein folds agrees well with the number of independent parameters estimated from a naïve Shannon sampling theorem approach. Estimates of the relative abundances of secondary structures were made using training/test sets derived from this data set. The average error in the estimate of α-helical content was 11%, and of β-sheet content was 9%. The distribution of proteins that are members of the four structure classes, α, β, α/β and α+β, are well separated in WAXS space when data extending to a spacing of 2.2 Å are used. Quantification of the information embedded within a WAXS pattern indicates that these data can be used as a powerful constraint in homology modeling of protein structures.     

Concordant exploration of the kinetics of RNA folding from global and local perspectives

Time-resolved small-angle X-ray scattering (SAXS) with millisecond time-resolution reveals two discrete phases of global compaction upon Mg2+-mediated folding of the Tetrahymena thermophila ribozyme. Electrostatic relaxation of the RNA occurs rapidly and dominates the first phase of compaction during which the observed radius of gyration (R(g)) decreases from 75 angstroms to 55 angstroms. A further decrease in R(g) to 45 angstroms occurs in a well-defined second phase. An analysis of mutant ribozymes shows that the latter phase depends upon the formation of long-range tertiary contacts within the P4-P6 domain of the ribozyme; disruption of the three remaining long-range contacts linking the peripheral helices has no effect on the 55-45 angstroms compaction transition. A better understanding of the role of specific tertiary contacts in compaction was obtained by concordant time-resolved hydroxyl radical (OH) analyses that report local changes in the solvent accessibility of the RNA backbone. Comparison of the global and local measures of folding shows that formation of a subset of native tertiary contacts (i.e. those defining the ribozyme core) can occur within a highly compact ensemble whose R(g) is close to that of the fully folded ribozyme. Analyses of additional ribozyme mutants and reaction conditions establish the generality of the rapid formation of a partially collapsed state with little to no detectable tertiary structure. These studies directly link global RNA compaction with formation of tertiary structure as the molecule acquires its biologically active structure, and underscore the strong dependence on salt of both local and global measures of folding kinetics.     

Evolution of quaternary structure in a homotetrameric enzyme

Dihydrodipicolinate synthase (DHDPS) is an essential enzyme in (S)-lysine biosynthesis and an important antibiotic target. All X-ray crystal structures solved to date reveal a homotetrameric enzyme. In order to explore the role of this quaternary structure, dimeric variants of Escherichia coli DHDPS were engineered and their properties were compared to those of the wild-type tetrameric form. X-ray crystallography reveals that the active site is not disturbed when the quaternary structure is disrupted. However, the activity of the dimeric enzymes in solution is substantially reduced, and a tetrahedral adduct of a substrate analogue is observed to be trapped at the active site in the crystal form. Remarkably, heating the dimeric enzymes increases activity. We propose that the homotetrameric structure of DHDPS reduces dynamic fluctuations present in the dimeric forms and increases specificity for the first substrate, pyruvate. By restricting motion in a key catalytic motif, a competing, non-productive reaction with a substrate analogue is avoided. Small-angle X-ray scattering and mutagenesis data, together with a B-factor analysis of the crystal structures, support this hypothesis and lead to the suggestion that in at least some cases, the evolution of quaternary enzyme structures might serve to optimise the dynamic properties of the protein subunits.     

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

1. Download : Download high-res image (195KB)
2. Download : Download full-size image

     

Structural Rearrangements Linked to Global Folding Pathways of the Azoarcus Group I Ribozyme

Stable RNAs must fold into specific three-dimensional structures to be biologically active, yet many RNAs form metastable structures that compete with the native state. Our previous time-resolved footprinting experiments showed that Azoarcus group I ribozyme forms its tertiary structure rapidly (τ < 30 ms) without becoming significantly trapped in kinetic intermediates. Here, we use stopped-flow fluorescence spectroscopy to probe the global folding kinetics of a ribozyme containing 2-aminopurine in the loop of P9. The modified ribozyme was catalytically active and exhibited two equilibrium folding transitions centered at 0.3 and 1.6 mM Mg²⁺, consistent with previous results. Stopped-flow fluorescence revealed four kinetic folding transitions with observed rate constants of 100, 34, 1, and 0.1 s^− 1 at 37 °C. From comparison with time-resolved Fe(II)-ethylenediaminetetraacetic acid footprinting of the modified ribozyme under the same conditions, these folding transitions were assigned to formation of the I_C intermediate, tertiary folding and docking of the nicked P9 tetraloop, reorganization of the P3 pseudoknot, and refolding of nonnative conformers, respectively. The footprinting results show that 50-60% of the modified ribozyme folds in less than 30 ms, while the rest of the RNA population undergoes slow structural rearrangements that control the global folding rate. The results show how small perturbations to the structure of the RNA, such as a nick in P9, populate kinetic folding intermediates that are not observed in the natural ribozyme.