Search

Background: In all organisms, the ribosome performs the unique and essential function of translating mRNA into protein. Since deletion of any ribosomal protein (RP) is embryonic lethal in complex eukaryotes, the ribosome is believed to be structurally uniform throughout the organism, with all RPs essential for cell viability. Results: We tested the structural homogeneity of RP use in ribosomes for normal and cancer tissues and cell lines using RP mRNA, ribosome profiling, and protein data. Analysis of RP mRNA data from 11,688 normal samples for 53 human tissues from 714 subjects and 10,363 tumor samples for 33 human cancers, normalized by total RP mRNA level per sample, showed that both normal (non-diseased) and tumor samples cluster by tissue type in humans. Matrix factorization showed that at least 3 RP mRNA signatures are necessary to describe normal blood and brain tissues, and a minimum of 16 RP mRNA signatures are necessary to describe data from 53 different normal tissues. A pan-cancer analysis of copy number variation (CNV) in 10,845 tumor samples for 33 human cancers showed that loss of one or both copies of RP genes was prevalent in every cancer type. Furthermore, there was no association between the number of double-deleted RP genes in tumors and patient survival, showing that RP loss does not reduce tumor fitness. CRISPR-Cas9 deletion of RP genes shows that many RPs are not essential in many cell lines. In several cancers, multiple RP-subtypes exist, with significant survival and genomic differences among them. These RP-subtypes often map to known molecular subtypes, and various RP genes are often deleted in one or both subtypes. This suggests that when genomic landscapes are modified in tumors, genes coding RPs are often lost, but these losses do not affect tumor viability. Analysis of mRNA data and ribosome profiling data of RP genes for cells and tissues from human, mouse and rat showed that these are highly correlated, showing that transcripts encoding ribosomal proteins are being translated into ribosomal proteins at rates proportional to their mRNA levels. Consistently, both mRNA data and ribosome profiling data of RP genes, normalized by total level per sample, showed tissue specific and development-stage specific clusters. Finally, analysis of RP protein levels in human adult and fetal tissues, standardized per sample, showed both development-stage and tissue specificity, showing that there is both tissue and development-stage specific heterogeneity of RP protein usage. Conclusions: These results suggest that there are multiple ribosome types in complex eukaryotes, with different RP composition which are regulated in a tissue and development-stage specific manner by some novel, yet unknown mechanism. Citation Format: Anshuman Panda, Anupama Yadav, Huwate Yeerna, Amartya Singh, Michael Biehl, Markus Lux, Alexander Schulz, Tyler Klecha, Sebastian Doniach, Hossein Khiabanian, Shridar Ganesan, Pablo Tamayo, Gyan Bhanot. The composition of the human ribosome varies significantly in different normal and malignant tissues [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 5865.

We give results from a detailed analysis of human Ribosomal Protein (RP) levels in normal and cancer samples and cell lines from large mRNA, copy number variation and ribosome profiling datasets. After normalizing total RP mRNA levels per sample, we find highly consistent tissue specific RP mRNA signatures in normal and tumor samples. Multiple RP mRNA-subtypes exist in several cancers, with significant survival and genomic differences. Some RP mRNA variations among subtypes correlate with copy number loss of RP genes. In kidney cancer, RP subtypes map to molecular subtypes related to cell-of-origin. Pan-cancer analysis of TCGA data showed widespread single/double copy loss of RP genes, without significantly affecting survival. In several cancer cell lines, CRISPR-Cas9 knockout of RP genes did not affect cell viability. Matched RP ribosome profiling and mRNA data in humans and rodents stratified by tissue and development stage and were strongly correlated, showing that RP translation rates were proportional to mRNA levels. In a small dataset of human adult and fetal tissues, RP protein levels showed development stage and tissue specific heterogeneity of RP levels. Our results suggest that heterogeneous RP levels play a significant functional role in cellular physiology, in both normal and disease states. © The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.

An atomic twinning structure is observed by averaging intensity correlations from many snapshots of gold nanoparticles in solution., During X-ray exposure of a molecular solution, photons scattered from the same molecule are correlated. If molecular motion is insignificant during exposure, then differences in momentum transfer between correlated photons are direct measurements of the molecular structure. In conventional small- and wide-angle solution scattering, photon correlations are ignored. This report presents advances in a new biomolecular structural analysis technique, correlated X-ray scattering (CXS), which uses angular intensity correlations to recover hidden structural details from molecules in solution. Due to its intense rapid pulses, an X-ray free electron laser (XFEL) is an excellent tool for CXS experiments. A protocol is outlined for analysis of a CXS data set comprising a total of half a million X-ray exposures of solutions of small gold nanoparticles recorded at the Spring-8 Ångström Compact XFEL facility (SACLA). From the scattered intensities and their correlations, two populations of nanoparticle domains within the solution are distinguished: small twinned, and large probably non-twinned domains. It is shown analytically how, in a solution measurement, twinning information is only accessible via intensity correlations, demonstrating how CXS reveals atomic-level information from a disordered solution of like molecules.

Gold labels provide a sequence–to–low-resolution structure map in SAXS measurements of biological macromolecules., Small-angle x-ray scattering (SAXS) is a powerful technique to probe the structure of biological macromolecules and their complexes under virtually arbitrary solution conditions, without the need for crystallization. While it is possible to reconstruct molecular shapes from SAXS data ab initio, the resulting electron density maps have a resolution of ~1 nm and are often insufficient to reliably assign secondary structure elements or domains. We show that SAXS data of gold-labeled samples significantly enhance the information content of SAXS measurements, allowing the unambiguous assignment of macromolecular sequence motifs to specific locations within a SAXS structure. We first demonstrate our approach for site-specifically internally and end-labeled DNA and an RNA motif. In addition, we present a protocol for highly uniform and site-specific labeling of proteins with small (~1.4 nm diameter) gold particles and apply our method to the signaling protein calmodulin. In all cases, the position of the small gold probes can be reliably identified in low-resolution electron density maps. Enhancing low-resolution measurements by site-selective gold labeling provides an attractive approach to aid modeling of a large range of macromolecular systems.

In this chapter we focus on scattering from non-crystalline solutions of molecules or nanoparticles in which the scattering objects are rotationally disordered.

Using mRNA-Seq and clinical data for 469 clear cell Renal Cell Carcinoma (ccRCC) samples from The Cancer Genome Atlas (TCGA), we develop a protocol to identify patients likely to have early recurrence of their disease. We first split the data into two sets, with 380 samples in the training set and 89 samples in the test set. Using the training set, we identify genes whose outlier status (high or low mRNA expression) is predictive of recurrence, based on Kaplan-Meier recurrence free survival log-rank p-value. We find a significant overlap among genes identified as predictive biomarkers in Reads per Kilobase Million (RPKM) normalized data and Raw Reads mRNA-Seq data. Using 80 consensus genes predictive in both RPKM and Raw Reads data, we define an outlier-based risk score R to stratify patients into two groups, a high-risk (early recurrence) group (R 2). The KM recurrence curve using this stratification shows excellent separation in training and test sets. Restricting the analysis to patients who had recurrence within two years (109 cases) and those who had no recurrence in five years (107 cases) we find that the risk predictor achieves ca. 80 percent sensitivity and specificity. The 80 genes identified by the outlier analysis were used to develop a more intuitive classifier based on Generalized Matrix Learning Vector Quantization (GMLVQ). This method stratifies samples into risk classes based on defining prototypes in feature space and an appropriate distance metric. GMLVQ identified a subset of 12 genes that have high accuracy in predicting recurrence, which suggests that an assay with a small number of genes might be able to predict recurrence in ccRCC.

Accurate determination of molecular distances is fundamental to understanding the structure, dynamics, and conformational ensembles of biological macromolecules. Here we present a method to determine the full distance distribution between small (∼7 A radius) gold labels attached to macromolecules with very high-precision (≤1 A) and on an absolute distance scale. Our method uses anomalous small-angle X-ray scattering close to a gold absorption edge to separate the gold–gold interference pattern from other scattering contributions. Results for 10–30 bp DNA constructs achieve excellent signal-to-noise and are in good agreement with previous results obtained by single-energy SAXS measurements without requiring the preparation and measurement of single labeled and unlabeled samples. The use of small gold labels in combination with ASAXS read out provides an attractive approach to determining molecular distance distributions that will be applicable to a broad range of macromolecular systems.

In vivo chromosomal behavior is dictated by the organization of genomic DNA at length scales ranging from nanometers to microns. At these disparate scales, the DNA conformation is influenced by a range of proteins that package, twist and disentangle the DNA double helix, leading to a complex hierarchical structure that remains undetermined. Thus, there is a critical need for methods of structural characterization of DNA that can accommodate complex environmental conditions over biologically relevant length scales. Based on multiscale molecular simulations, we report on the possibility of measuring supercoiling in complex environments using angular correlations of scattered X-rays resulting from X-ray free electron laser (xFEL) experiments. We recently demonstrated the observation of structural detail for solutions of randomly oriented metallic nanoparticles [D. Mendez et al., Philos. Trans. R. Soc. B 360 (2014) 20130315]. Here, we argue, based on simulations, that correlated X-ray scattering (CXS) has the potential for measuring the distribution of DNA folds in complex environments, on the scale of a few persistence lengths.

RNA folding is enabled by interactions between the nucleic acid and its ion atmosphere, the mobile sheath of aqueous ions that surrounds and stabilizes it. Understanding the ion atmosphere requires the interplay of experiment and theory. However, even an apparently simple experiment to probe the ion atmosphere, measuring the dependence of DNA duplex stability upon ion concentration and identity, suffers from substantial complexity, because the unfolded ensemble contains many conformational states that are difficult to treat accurately with theory. To minimize this limitation, we measured the unfolding equilibrium of a DNA hairpin using a single-molecule optical trapping assay, in which the unfolded state is constrained to a limited set of elongated conformations. The unfolding free energy increased linearly with the logarithm of monovalent cation concentration for several cations, such that smaller cations tended to favor the folded state. Mg(2+) stabilized the hairpin much more effectively at low concentrations than did any of the monovalent cations. Poisson-Boltzmann theory captured trends in hairpin stability measured for the monovalent cation titrations with reasonable accuracy, but failed to do so for the Mg(2+) titrations. This finding is consistent with previous work, suggesting that Poisson-Boltzmann and other mean-field theories fail for higher valency cations where ion-ion correlation effects may become significant. The high-resolution data herein, because of the straightforward nature of both the folded and the unfolded states, should serve as benchmarks for the development of more accurate electrostatic theories that will be needed for a more quantitative and predictive understanding of nucleic acid folding.

International audience; Small Angle X-ray Scattering (SAXS) techniques are becoming more and more useful for structural biologists and biochemists, thanks to better access to dedicated synchrotron beamlines, better detectors and the relative easiness of sample preparation. The ability to compute the theoretical SAXS profile of a given structural model, and to compare this profile with the measured scattering intensity, yields crucial structural informations about the macromolecule under study and/or its complexes in solution. An important contribution to the profile, besides the macromolecule itself and its solvent-excluded volume, is the excess density due to the hydration layer. AquaSAXS takes advantage of recently developed methods, such as AquaSol, that give the equilibrium solvent density map around macromolecules, to compute an accurate SAXS/WAXS profile of a given structure and to compare it to the experimental one. Here, we describe the interface architecture and capabilities of the AquaSAXS web server (http://lorentz.dynstr.pasteur.fr/aquasaxs.php).

Membrane proteins impose enormous challenges to cellular protein homeostasis during their post-translational targeting, and they require chaperones to keep them soluble and translocation competent. Here we show that a novel targeting factor in the chloroplast signal recognition particle (cpSRP), cpSRP43, is a highly specific molecular chaperone that efficiently reverses the aggregation of its substrate proteins. In contrast to 'ATPases associated with various cellular activities' (AAA(+)) chaperones, cpSRP43 uses specific binding interactions with its substrate to mediate its 'disaggregase' activity. This disaggregase capability can allow targeting machineries to more effectively capture their protein substrates and emphasizes a close connection between protein folding and trafficking processes. Moreover, cpSRP43 provides the first example to our knowledge of an ATP-independent disaggregase and shows that efficient reversal of protein aggregation can be attained by specific binding interactions between a chaperone and its substrate.

Riboswitches are gene-regulating RNAs that are usually found in the 5′-untranslated regions of messenger RNA. As the sugar-phosphate backbone of RNA is highly negatively charged, the folding and ligand-binding interactions of riboswitches are strongly dependent on the presence of cations. Using small angle X-ray scattering (SAXS) and hydroxyl radical footprinting, we examined the cation dependence of the different folding stages of the glycine-binding riboswitch from Vibrio cholerae. We found that the partial folding of the tandem aptamer of this riboswitch in the absence of glycine is supported by all tested mono- and divalent ions, suggesting that this transition is mediated by nonspecific electrostatic screening. Poisson–Boltzmann calculations using SAXS-derived low-resolution structural models allowed us to perform an energetic dissection of this process. The results showed that a model with a constant favorable contribution to folding that is opposed by an unfavorable electrostatic term that varies with ion concentration and valency provides a reasonable quantitative description of the observed folding behavior. Glycine binding, on the other hand, requires specific divalent ions binding based on the observation that Mg2+, Ca2+, and Mn2+ facilitated glycine binding, whereas other divalent cations did not. The results provide a case study of how ion-dependent electrostatic relaxation, specific ion binding, and ligand binding can be coupled to shape the energetic landscape of a riboswitch and can begin to be quantitatively dissected.

Structured RNAs must fold into their native structures and discriminate against a large number of alternative ones, an especially difficult task given the limited information content of RNA's nucleotide alphabet. The simplest motifs within structured RNAs are two helices joined by nonhelical junctions. To uncover the fundamental behavior of these motifs and to elucidate the underlying physical forces and challenges faced by structured RNAs, we computationally and experimentally studied a tethered duplex model system composed of two helices joined by flexible single- or double-stranded polyethylene glycol tethers, whose lengths correspond to those typically observed in junctions from structured RNAs. To dissect the thermodynamic properties of these simple motifs, we computationally probed how junction topology, electrostatics, and tertiary contact location influenced folding stability. Small-angle X-ray scattering was used to assess our predictions. Single- or double-stranded junctions, independent of sequence, greatly reduce the space of allowed helical conformations and influencing the preferred location and orientation of their adjoining helices. A double-stranded junction guides the helices along a hinge-like pathway. In contrast, a single-stranded junction samples a broader set of conformations and has different preferences than the double-stranded junction. In turn, these preferences determine the stability and distinct specificities of tertiary structure formation. These sequence-independent effects suggest that properties as simple as a junction's topology can generally define the accessible conformational space, thereby stabilizing desired structures and assisting in discriminating against misfolded structures. Thus, junction topology provides a fundamental strategy for transcending the limitations imposed by the low information content of RNA primary sequence.

Using small-angle x-ray diffraction from centrifugally oriented acetylcholine receptor (AChR) enriched membranes coupled with anomalous scattering from terbium ions (Tb3+) titrated into presumed Ca2+ binding sites, we have mapped the distribution of Tb3+ perpendicular to the membrane plane using a heavy atom refinement algorithm. We have compared the distribution of Tb3+ in the closed resting state with that in the carbamylcholine-desensitized state. In the closed resting state we find 45 Tb3+ ions distributed in 10 narrow peaks perpendicular to the membrane plane. Applying the same refinement procedure to the data from carbamylcholine desensitized AChR we find 18 fewer Tb3+ ions in eight peaks, and slight rearrangements of Tb3+ density in the peaks near the ends of the AChR ion channel pore. These agonist dependent changes in the Tb3+ stoichiometry and distribution suggest a likely role for multivalent cations in stabilizing the different functional states of the AChR, and the changes in the Tb3+ distribution at the two ends of the pore suggest a potential role for multivalent cations in the gating of the ion channel.

Staphylococcal nuclease forms three different partially-folded intermediates at low pH in the presence of low to moderate concentration of anions, differing in the amount of secondary structure, globularity, stability, and compactness. Although these intermediates are monomeric at low protein concentration (< or =0.25 mg/mL), increasing concentrations of protein result in the formation of dimers and soluble oligomers, ultimately leading to larger insoluble aggregates. Unexpectedly, increasing protein concentration not only led to association, but also to increased structure of the intermediates. The secondary structure, stability, and globularity of the two less-ordered partially-folded intermediates (A1 and A2) were substantially increased upon association, suggesting that aggregation induces structure. An excellent correlation was found between degree of association and amount of structure measured by different techniques, including circular dichroism, fluorescence, Fourier transform infrared spectroscopy (FTIR), and small-angle X-ray scattering. The associated states were also substantially more stable toward urea denaturation than the monomeric forms. A mechanism is proposed, in which the observed association of monomeric intermediates involves intermolecular interactions which correspond to those found intramolecularly in normal folding to the native state.

The large electrostatic repulsion arising from the negatively-charged backbone of RNA molecules presents a large barrier to folding. Solution counterions assist in the folding process by screening this electrostatic repulsion. While early research interpreted the effect of these counterions in terms of an empirical ligand-binding model, theories based on physical models have supplanted them and revised our view of the roles that ions play in folding. Instead of specific ion-binding sites, most ions in solution interact inside an “ion atmosphere” – a fluctuating cloud of non-specifically associated ions surrounding any charged molecule. Recent advances in experiments have begun the task of characterizing the ion atmosphere, yielding valuable data that has revealed deficiencies in Poisson-Boltzmann theory, the most widely-used theory of the ion atmosphere. The continued development of experiments will help guide the development of improved theories, with the ultimate goal of understanding RNA folding and function and nucleic acid/protein interactions from a quantitative perspective.

Protein alpha-helices are ubiquitous secondary structural elements, seldom considered to be stable without tertiary contacts. However, amino acid sequences in proteins that are based on alternating repeats of four glutamic acid (E) residues and four positively charged residues, a combination of arginine (R) and lysine (K), have been shown to form stable alpha-helices in a few proteins, in the absence of tertiary interactions. Here, we find that this ER/K motif is more prevalent than previously reported, being represented in proteins of diverse function from archaea to humans. By using molecular dynamics (MD) simulations, we characterize a dynamic pattern of side-chain interactions that extends along the backbone of ER/K alpha-helices. A simplified model predicts that side-chain interactions alone contribute substantial bending rigidity (0.5 pN/nm) to ER/K alpha-helices. Results of small-angle x-ray scattering (SAXS) and single-molecule optical-trap analyses are consistent with the high bending rigidity predicted by our model. Thus, the ER/K alpha-helix is an isolated secondary structural element that can efficiently span long distances in proteins, making it a promising tool in designing synthetic proteins. We propose that the significant rigidity of the ER/K alpha-helix can help regulate protein function, as a force transducer between protein subdomains.

We have used small-angle X-ray solution scattering to obtain ab initio shape reconstructions of the complete VS ribozyme. The ribozyme occupies an electron density envelope with an irregular shape, into which helical sections have been fitted. The ribozyme is built around a core comprising a near-coaxial stack of three helices, organized by two three-way helical junctions. An additional three-way junction formed by an auxiliary helix directs the substrate stem-loop, juxtaposing the cleavage site with an internal loop to create the active complex. This is consistent with the current view of the probable mechanism of transesterification in which adenine and guanine nucleobases contributed by the interacting loops combine together in general acid-base catalysis.

Electrostatic forces, acting between helices and modulated by the presence of the ion atmosphere, are key determinants in the energetic balance that governs RNA folding. Previous studies have employed Poisson-Boltzmann (PB) theory to compute the energetic contribution of these forces in RNA folding. However, the complex interaction of these electrostatic forces with RNA features such as tertiary contact formation, specific ion-binding, and complex interhelical junctions present in prior studies precluded a rigorous evaluation of PB theory, especially in physiologically important Mg(2+) solutions. To critically assess PB theory, we developed a model system that isolates these electrostatic forces. The model system, composed of two DNA duplexes tethered by a polyethylene glycol junction, is an analog for the unfolded state of canonical helix-junction-helix motifs found in virtually all structured RNAs. This model system lacks the complicating features that have precluded a critical assessment of PB in prior studies, ensuring that interhelical electrostatic forces dominate the behavior of the system. The system's simplicity allows PB predictions to be directly compared with small-angle X-ray scattering experiments over a range of monovalent and divalent ion concentrations. These comparisons indicate that PB is a reasonable description of the underlying electrostatic energies for monovalent ions, but large deviations are observed for divalent ions. The validation of PB for monovalent solutions allows analysis of the change in the conformational ensemble of this simple motif as salt concentration is changed. Addition of ions allows the motif to sample more compact microstates, increasing its conformational entropy. The increase of conformational entropy presents an additional barrier to folding by stabilizing the unfolded state. Neglecting this effect will adversely impact the accuracy of folding analyses and models.

Myosin VI has challenged the lever arm hypothesis of myosin movement because of its ability to take approximately 36-nm steps along actin with a canonical lever arm that seems to be too short to allow such large steps. Here we demonstrate that the large step of dimeric myosin VI is primarily made possible by a medial tail in each monomer that forms a rare single alpha-helix of approximately 10 nm, which is anchored to the calmodulin-bound IQ domain by a globular proximal tail. With the medial tail contributing to the approximately 36-nm step, rather than dimerizing as previously proposed, we show that the cargo binding domain is the dimerization interface. Furthermore, the cargo binding domain seems to be folded back in the presence of the catalytic head, constituting a potential regulatory mechanism that inhibits dimerization.

We present a systematic analysis of the aggregation number and shape of micelles formed by nine detergents commonly used in the study of membrane proteins. Small-angle X-ray scattering measurements are reported for glucosides with 8 and 9 alkyl carbons (OG/NG), maltosides and phosphocholines with 10 and 12 alkyl carbons (DM/DDM and FC-10/FC-12), 1,2-dihexanoyl-sn-glycero-phosphocholine (DHPC), 1-palmitoyl-2-hydroxy-sn-glycero-3-[phospho-rac-(1-glycerol)] (LPPG), and 3-[(3-cholamidopropyl)dimethylammonio]-1-propane sulfonate (CHAPS). The SAXS intensities are well described by two-component ellipsoid models, with a dense outer shell corresponding to the detergent head groups and a less electron dense hydrophobic core. These models provide an intermediate resolution view of micelle size and shape. In addition, we show that Guinier analysis of the forward scattering intensity can be used to obtain an independent and model-free measurement of the micelle aggregation number and radius of gyration. This approach has the advantage of being easily generalizable to protein-detergent complexes, where simple geometric models are inapplicable. Furthermore, we have discovered that the position of the second maximum in the scattering intensity provides a direct measurement of the characteristic head group-head group spacing across the micelle core. Our results for the micellar aggregation numbers and dimensions agree favorably with literature values as far as they are available. We de novo determine the shape of FC-10, FC-12, DM, LPPG, and CHAPS micelles and the aggregation numbers of FC-10 and OG to be ca. 50 and 250, respectively. Combined, these data provide a comprehensive view of the determinants of micelle formation and serve as a starting point to correlate detergent properties with detergent-protein interactions.

Small-angle X-ray scattering can be a valuable tool in the structural characterization of membrane protein–detergent complexes (PDCs). However, a major challenge is to separate the PDC scattering signal from that of the `empty' detergent micelle in a protein–detergent mixture. We briefly review an approach that allows approximate determination of the PDC scattering signal at low momentum transfer and present a novel approach that employs a singular value decomposition (SVD) and fitting of scattering data collected at different protein–detergent stoichiometries. The SVD approach allows the scattering profile for the PDC over the entire measured momentum transfer range to be obtained, it is applicable to strongly scattering detergents and can take into account interparticle interference. The two approaches are contrasted and an application to the membrane protein TM0026 from Thermotoga maritima is presented.

Several algorithms are available to reconstruct low-resolution electron density maps of biological macromolecules from small-angle solution scattering data. These algorithms have been extensively applied to proteins and protein complexes. Here, we demonstrate their applicability to nucleic acids by reconstructing a set of RNA and DNA molecules of known three-dimensional structure from their small-angle X-ray scattering profiles. The overall size and shape of the molecules get reproduced well in all tested cases. Furthermore, we show that the generated bead models can be used as inputs for electrostatic calculations. The number of ions bound under different solution conditions computed from numerical solutions of the Poisson–Boltzmann equation for bead models agrees very well with results of calculations on all atom models derived from crystallography. The predictions from Poisson–Boltzmann theory also agree generally well with experimentally determined ion binding numbers.

Intense femtosecond x-ray pulses from free-electron laser sources allow the imaging of individual particles in a single shot. Early experiments at the Linac Coherent Light Source (LCLS) have led to rapid progress in the field and, so far, coherent diffractive images have been recorded from biological specimens, aerosols, and quantum systems with a few-tens-of-nanometers resolution. In March 2014, LCLS held a workshop to discuss the scientific and technical challenges for reaching the ultimate goal of atomic resolution with single-shot coherent diffractive imaging. This paper summarizes the workshop findings and presents the roadmap toward reaching atomic resolution, 3D imaging at free-electron laser sources.

We recently used computational protein design to create an extremely stable, globular protein, Top7, with a sequence and fold not observed previously in nature. Since Top7 was created in the absence of genetic selection, it provides a rare opportunity to investigate aspects of the cellular protein production and surveillance machinery that are subject to natural selection. Here we show that a portion of the Top7 protein corresponding to the final 49 C-terminal residues is efficiently mis-translated and accumulates at high levels in Escherichia coli. We used circular dichroism, size-exclusion chromatography, small-angle X-ray scattering, analytical ultra-centrifugation, and NMR spectroscopy to show that the resulting C-terminal fragment (CFr) protein adopts a compact, extremely stable, homo-dimeric structure. Based on the solution structure, we engineered an even more stable variant of CFr by disulfide-induced covalent circularisation that should be an excellent platform for design of novel functions. The accumulation of high levels of CFr exposes the high error rate of the protein translation machinery. The rarity of correspondingly stable fragments in natural proteins coupled with the observation that high quality ribosome binding sites are found to occur within E. coli protein-coding regions significantly less often than expected by random chance implies a stringent evolutionary pressure against protein sub-fragments that can independently fold into stable structures. The symmetric self-association between two identical mis-translated CFr sub-domains to generate an extremely stable structure parallels a mechanism for natural protein-fold evolution by modular recombination of protein sub-structures.

Structural studies of integral membrane proteins typically rely upon detergent micelles as faithful mimics of the native lipid bilayer. Therefore, membrane protein structure determination would be greatly facilitated by biophysical techniques that are capable of evaluating and assessing the fold and oligomeric state of these proteins solubilized in detergent micelles. In this study, an approach to the characterization of detergent-solubilized integral membrane proteins is presented. Eight Thermotoga maritima membrane proteins were screened for solubility in 11 detergents, and the resulting soluble protein–detergent complexes were characterized with small angle X-ray scattering (SAXS), nuclear magnetic resonance (NMR) spectroscopy, circular dichroism (CD) spectroscopy, and chemical cross-linking to evaluate the homogeneity, oligomeric state, radius of gyration, and overall fold. A new application of SAXS is presented, which does not require density matching, and NMR methods, typically used to evaluate soluble proteins, are successfully applied to detergent-solubilized membrane proteins. Although detergents with longer alkyl chains solubilized the most proteins, further characterization indicates that some of these protein–detergent complexes are not well suited for NMR structure determination due to conformational exchange and protein oligomerization. These results emphasize the need to screen several different detergents and to characterize the protein–detergent complex in order to pursue structural studies. Finally, the physical characterization of the protein–detergent complexes indicates optimal solution conditions for further structural studies for three of the eight overexpressed membrane proteins.

Polyproline type II (PPII) helix has emerged recently as the dominant paradigm for describing the conformation of unfolded polypeptides. However, most experimental observables used to characterize unfolded proteins typically provide only short-range, sequence-local structural information that is both time- and ensemble-averaged, giving limited detail about the long-range structure of the chain. Here, we report a study of a long-range property: the radius of gyration of an alanine-based peptide, Ace-(diaminobutyric acid) 2 -(Ala) 7 -(ornithine) 2 -NH2. This molecule has previously been studied as a model for the unfolded state of proteins under folding conditions and is believed to adopt a PPII fold based on short-range techniques such as NMR and CD. By using synchrotron radiation and small-angle x-ray scattering, we have determined the radius of gyration of this peptide to be 7.4 ± 0.5 Å, which is significantly less than the value expected from an ideal PPII helix in solution (13.1 Å). To further study this contradiction, we have used molecular dynamics simulations using six variants of the AMBER force field and the GROMOS 53A6 force field. However, in all cases, the simulated ensembles underestimate the PPII content while overestimating the experimental radius of gyration. The conformational model that we propose, based on our small angle x-ray scattering results and what is known about this molecule from before, is that of a very flexible, fluctuating structure that on the level of individual residues explores a wide basin around the ideal PPII geometry but is never, or only rarely, in the ideal extended PPII helical conformation.

We study the early steps of amyloid formation of the seven residue peptide GNNQQNY from yeast prion-like protein Sup35 by simulating the random coil to beta-sheet and alpha-helix to beta-sheet transition both in the absence and presence of a cross-beta amyloid nucleus. The simulation method at atomic resolution employs a new implementation of a Langevin dynamics "reaction path annealing" algorithm. The results indicate that the presence of amyloid-like cross-beta-sheet strands both facilitates the transition into the cross-beta conformation and substantially lowers the free energy barrier for this transition. This model systems allows us to investigate the energetic and kinetic details of this transition, which is consistent with an auto-catalyzed, nucleation-like mechanism for the formation of beta-amyloid. In particular, we find that electrostatic interactions of peptide backbone dipoles contribute significantly to the stability of the beta-amyloid state. Furthermore, we find water exclusion and interactions of polar side-chains to be driving forces of amyloid formation: the cross-beta conformation is stabilized by burial of polar side-chains and inter-residue hydrogen bonds in the presence of an amyloid-like "seed". These findings are in support of a "dry, polar zipper model" of amyloid formation.

Summary The open/closed transition in polymerases is a crucial event in DNA replication and transcription. We hypothesize that the residues that transmit the signal for the open/closed transition are also strongly conserved. To identify the dynamically relevant residues, we use an elastic network model of polymerases and probe the residue-specific response to a local perturbation. In a variety of DNA/RNA polymerases, a network of residues spanning the fingers and palm domains is involved in the open/closed transition. The similarity in the network of residues responsible for large-scale domain movements supports the notion of a common induced-fit mechanism in the polymerase families for the formation of a closed ternary complex. Multiple sequence alignment shows that many of these residues are also strongly conserved. Residues with the largest sensitivity to local perturbations include those that are not so obviously involved in the polymerase catalysis. Our results suggest that mutations of the mechanical "hot spots" can compromise the efficiency of the enzyme.

Spectroscopic studies have identified a number of proteins that appear to retain significant residual structure under even strongly denaturing conditions. Intrinsic viscosity, hydrodynamic radii, and small-angle x-ray scattering studies, in contrast, indicate that the dimensions of most chemically denatured proteins scale with polypeptide length by means of the power-law relationship expected for random-coil behavior. Here we further explore this discrepancy by expanding the length range of characterized denatured-state radii of gyration ( R G ) and by reexamining proteins that reportedly do not fit the expected dimensional scaling. We find that only 2 of 28 crosslink-free, prosthetic-group-free, chemically denatured polypeptides deviate significantly from a power-law relationship with polymer length. The R G of the remaining 26 polypeptides, which range from 16 to 549 residues, are well fitted ( r 2 = 0.988) by a power-law relationship with a best-fit exponent, 0.598 ± 0.028, coinciding closely with the 0.588 predicted for an excluded volume random coil. Therefore, it appears that the mean dimensions of the large majority of chemically denatured proteins are effectively indistinguishable from the mean dimensions of a random-coil ensemble.

Fibrillar deposits of insulin cause serious problems in implantable insulin pumps, commercial production of insulin, and for some diabetics. We performed a systematic investigation of the effect of urea-induced structural perturbations on the mechanism of fibrillation of insulin. The addition of as little as 0.5 m urea to zinc-bound hexameric insulin led to dissociation into dimers. Moderate concentrations of urea led to accumulation of a partially unfolded dimer state, which dissociates into an expanded, partially folded monomeric state. Very high concentrations of urea resulted in an unfolded monomer with some residual structure. The addition of even very low concentrations of urea resulted in increased fibrillation. Accelerated fibrillation correlated with population of the partially folded intermediates, which existed at up to 8 m urea, accounting for the formation of substantial amounts of fibrils under such conditions. Under monomeric conditions the addition of low concentrations of urea slowed down the rate of fibrillation, e.g. 5-fold at 0.75 m urea. The decreased fibrillation of the monomer was due to an induced non-native conformation with significantly increased α-helical content compared with the native conformation. The data indicate a close-knit relationship between insulin conformation and propensity to fibrillate. The correlation between fibrillation and the partially unfolded monomer indicates that the latter is a critical amyloidogenic intermediate in insulin fibrillation.

Native zinc-bound insulin exists as a hexamer at neutral pH. Under destabilizing conditions, the hexamer dissociates, and is very prone to forming fibrils. Insulin fibrils exhibit the typical properties of amyloid fibrils, and pose a problem in the purification, storage, and delivery of therapeutic insulin solutions. We have carried out a systematic investigation of the effect of guanidine hydrochloride (Gdn.HCl)-induced structural perturbations on the mechanism of fibrillation of insulin. At pH 7.4, the addition of as little as 0.25 M Gdn.HCl leads to dissociation of insulin hexamers into dimers. Moderate concentrations of Gdn.HCl lead to formation of a novel partially unfolded dimer state, which dissociates into a partially unfolded monomer state. High concentrations of Gdn.HCl resulted in unfolded monomers with some residual structure. The addition of even very low concentrations of Gdn.HCl resulted in substantially accelerated fibrillation, although the yield of fibrils decreased at high concentrations. Accelerated fibrillation correlated with the population of the expanded (partially folded) monomer, which existed up to6 M Gdn.HCl, accounting for the formation of substantial amounts of fibrils under such conditions. In the presence of 20% acetic acid, where insulin exists as the monomer, fibrillation was also accelerated by Gdn.HCl. The enhanced fibrillation of the monomer was due to the increased ionic strength at low denaturant concentrations, and due to the presence of the partially unfolded, expanded conformation at Gdn.HCl concentrations above 1 M. The data suggest that under physiological conditions, the fibrillation of insulin involves both changes in the association state (with rate-limiting hexamer dissociation) and conformational changes, leading to formation of the amyloidogenic expanded monomer intermediate.

Large RNAs can collapse into compact conformations well before the stable formation of the tertiary contacts that define their final folds. This study identifies likely physical mechanisms driving these early compaction events in RNA folding. We have employed time-resolved small-angle X-ray scattering to monitor the fastest global shape changes of the Tetrahymena ribozyme under different ionic conditions and with RNA mutations that remove long-range tertiary contacts. A partial collapse in each of the folding time-courses occurs within tens of milliseconds with either monovalent or divalent cations. Combined with comparison to predictions from structural models, this observation suggests a relaxation of the RNA to a more compact but denatured conformational ensemble in response to enhanced electrostatic screening at higher ionic concentrations. Further, the results provide evidence against counterion-correlation-mediated attraction between RNA double helices, a recently proposed model for early collapse. A previous study revealed a second 100 ms phase of collapse to a globular state. Surprisingly, we find that progression to this second early folding intermediate requires RNA sequence motifs that eventually mediate native long-range tertiary interactions, even though these regions of the RNA were observed to be solvent-accessible in previous footprinting studies under similar conditions. These results help delineate an analogy between the early conformational changes in RNA folding and the "burst phase" changes and molten globule formation in protein folding.

The aggregation of alpha-synuclein is believed to play an important role in the pathogenesis of Parkinson's disease as well as other neurodegenerative disorders ("synucleinopathies"). However, the function of alpha-synuclein under physiologic and pathological conditions is unknown, and the mechanism of alpha-synuclein aggregation is not well understood. Here we show that alpha-synuclein forms a tight 2:1 complex with histones and that the fibrillation rate of alpha-synuclein is dramatically accelerated in the presence of histones in vitro. We also describe the presence of alpha-synuclein and its co-localization with histones in the nuclei of nigral neurons from mice exposed to a toxic insult (i.e., injections of the herbicide paraquat). These observations indicate that translocation into the nucleus and binding with histones represent potential mechanisms underlying alpha-synuclein pathophysiology.

Chaperonins use ATPase cycling to promote conformational changes leading to protein folding. The prokaryotic chaperonin GroEL requires a cofactor, GroES, which serves as a “lid” enclosing substrates in the central cavity and confers an asymmetry on GroEL required for cooperative transitions driving the reaction. The eukaryotic chaperonin TRiC/CCT does not have such a cofactor but appears to have a “built-in” lid. Whether this seemingly symmetric chaperonin also operates through an asymmetric cycle is unclear. We show that unlike GroEL, TRiC does not close its lid upon nucleotide binding, but instead responds to the trigonal-bipyramidal transition state of ATP hydrolysis. Further, nucleotide analogs inducing this transition state confer an asymmetric conformation on TRiC. Similar to GroEL, lid closure in TRiC confines the substrates in the cavity and is essential for folding. Understanding the distinct mechanisms governing eukaryotic and bacterial chaperonin function may reveal how TRiC has evolved to fold specific eukaryotic proteins.

Human insulin exists in different association states, from monomer to hexamer, depending on the conditions. In the presence of zinc the “normal” state is a hexamer. The structural properties of 20 variants of human insulin were studied by near-UV circular dichroism, fluorescence spectroscopy, and small-angle X-ray scattering (SAXS). The mutants showed different degrees of association (monomer, dimers, tetramers, and hexamers) at neutral pH. A correlation was shown between the accessibility of tyrosines to acrylamide quenching and the degree of association of the insulin mutants. The near-UV CD spectra of the insulins were affected by protein association and by mutation-induced structural perturbations. However, the shape and intensity of difference CD spectra, obtained by subtraction of the spectra measured in 20% acetic acid (where all insulin species were monomeric) from the corresponding spectra measured at neutral pH, correlate well with the degree of insulin association. In fact, the near-UV CD difference spectra for monomeric, dimeric, tetrameric, and hexameric insulin are very distinctive, both in terms of intensity and shape. The results show that the spectral properties of the insulins reflect their state of association, and can be used to predict their oligomeric state. © 2003 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 92:847–858, 2003

The pathological hallmark of Parkinson's disease is the presence of intracellular inclusions, Lewy bodies, and Lewy neurites, in the dopaminergic neurons of the substantia nigra and several other brain regions. Filamentous α-synuclein is the major component of these deposits and its aggregation is believed to play an important role in Parkinson's disease and several other neurodegenerative diseases. Two homologous proteins, β- and γ-synucleins, are also abundant in the brain. The synucleins are natively unfolded proteins. β-Synuclein, which lacks 11 central hydrophobic residues compared with its homologs, exhibited the properties of a random coil, whereas α- and γ-synucleins were slightly more compact and structured. γ-Synuclein, unlike its homologs, formed a soluble oligomer at relatively low concentrations, which appears to be an off-fibrillation pathway species. Here we show that, although they have similar biophysical properties to α-synuclein, β- And γ-synucleins inhibit α-synuclein fibril formation. Complete inhibition of α-synuclein fibrillation was observed at 4:1 molar excess of β- and γ-synucleins. No significant incorporation of β-synuclein into the fibrils was detected. The lack of fibrils formed by β-synuclein is most readily explained by the absence of a stretch of hydrophobic residues from the middle region of the protein. A model for the inhibition is proposed.

Light chain amyloidosis involves the systemic pathologic deposition of monoclonal light chain variable domains of immunoglobulins as insoluble fibrils. The variable domain LEN was obtained from a patient who had no overt amyloidosis; however, LEN forms fibrils in vitro, under mildly destabilizing conditions. The in vitro kinetics of fibrillation were investigated using a wide variety of probes. The rate of fibril formation was highly dependent on the initial protein concentration. In contrast to most amyloid systems, the kinetics became slower with increasing LEN concentrations. At high protein concentrations a significant lag in time was observed between the conformational changes and the formation of fibrils, consistent with the formation of soluble off-pathway oligomeric species and a branched pathway. The presence of off-pathway species was confirmed by small angle x-ray scattering. At low protein concentrations the structural rearrangements were concurrent with fibril formation, indicating the absence of formation of the off-pathway species. The data are consistent with a model for fibrillation in which a dimeric form of LEN (at high protein concentration) inhibits fibril formation by interaction with an intermediate on the fibrillation pathway and leads to formation of the off-pathway intermediate.

Detergent micelles are used in many areas of research and technology, in particular, as mimics of the cellular membranes in the purification and biochemical and structural characterization of membrane proteins. Applications of detergent micelles are often hindered by the limited set of properties of commercially available detergents. Mixtures of micelle-forming detergents provide a means to systematically obtain additional micellar properties and expand the repertoire of micelle features available; however, our understanding of the properties of detergent mixtures is still limited. In this study, the shape and size of binary mixtures of seven different detergents commonly used in molecular host-guest systems and membrane protein research were investigated. The data suggests that the detergents form ideally mixed micelles with sizes and shapes different from those of pure individual micelles. For most measurements of size, the mixtures varied linearly with detergent mole fraction and therefore can be calculated from the values of the pure detergents. We propose that properties such as the geometry, size, and surface charge can be systematically and predictably tuned for specific applications.

Ions surround nucleic acids in what is referred to as an ion atmosphere. As a result, the folding and dynamics of RNA and DNA and their complexes with proteins and with each other cannot be understood without a reasonably sophisticated appreciation of these ions' electrostatic interactions. However, the underlying behavior of the ion atmosphere follows physical rules that are distinct from the rules of site binding that biochemists are most familiar and comfortable with. The main goal of this review is to familiarize nucleic acid experimentalists with the physical concepts that underlie nucleic acid–ion interactions. Throughout, we provide practical strategies for interpreting and analyzing nucleic acid experiments that avoid pitfalls from oversimplified or incorrect models. We briefly review the status of theories that predict or simulate nucleic acid–ion interactions and experiments that test these theories. Finally, we describe opportunities for going beyond phenomenological fits to a next-generation, truly predictive understanding of nucleic acid–ion interactions.