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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

We have used small angle X-ray scattering (SAXS) to monitor changes in the overall size and shape of the Tetrahymena ribozyme as it folds. The native ribozyme, formed in the presence of Mg2+, is much more compact and globular than the ensemble of unfolded conformations. Time-resolved measurements show that most of the compaction occurs at least 20-fold faster than the overall folding to the native state, suggesting that a compact intermediate or family of intermediates is formed early and then rearranges in the slow steps that limit the overall folding rate. These results lead to a kinetic folding model in which an initial 'electrostatic collapse' of the RNA is followed by slower rearrangements of elements that are initially mispositioned.

Prothymosin alpha has previously been shown to be unfolded at neutral pH, thus belonging to a growing family of "natively unfolded" proteins. The structural properties and conformational stability of recombinant human prothymosin alpha were characterized at neutral and acidic pH by gel filtration, SAXS, circular dichroism, ANS fluorescence, (1)H NMR, and resistance to urea-induced unfolding. Interestingly, prothymosin alpha underwent a cooperative transition from the unfolded state into a partially folded conformation on lowering the pH. This conformation of prothymosin alpha is a compact denatured state, with structural properties different from those of the molten globule. The formation of alpha-helical structure by the glutamic acid-rich elements of the protein accompanied by the partial hydrophobic collapse is expected at lower pH due to the neutralization of the negatively charged residues. It is possible that such conformational changes may be associated with the protein function.

Solution X-ray scattering was used to study the equilibrium unfolding of cytochrome c as a function of guanidine hydrochloride concentration at neutral pH. The radius of gyration (Rg) shows a cooperative transition with increasing denaturant with a similar Cm to that observed with circular dichroism. However, the lack of an isoscattering point in the X-ray scattering patterns suggests the equilibrium unfolding is not simply a two-state process. Singular value decomposition (SVD) analysis was applied to the scattering patterns to determine the number of distinct scattering species. SVD analysis reveals the existence of three components, suggesting that at least three equilibrium states of the protein exist. A model was employed to determine the thermodynamic parameters and the scattering profiles of the three equilibrium states. These scattering profiles show that one state is native (N). The other two states (U1, U2) are unfolded, with U2 being fully unfolded and U1 having some residual structure. Using the thermodynamic parameters to calculate fractional populations, U1 is maximally populated at intermediate denaturant concentrations while U2 is maximally populated at high denaturant concentrations. It is likely that there is a multiplicity of denatured states with U1 and U2 representing an average of the denatured states populated at intermediate and high denaturant concentrations, respectively.

We present an algorithm for simulating the long time scale dynamics of proteins and other macromolecules. Our method applies the concept of multiple time step integration to the diffusive Langevin equation, in which short time scale dynamics are replaced by friction and noise. The macromolecular force field is represented at atomic resolution. Slow motions are modeled by constrained Langevin dynamics with very large time steps, while faster degrees of freedom are kept in local thermal equilibrium. In the limit of a sufficiently large molecule, our algorithm is shown to reduce the CPU time required by two orders of magnitude. We test the algorithm on two systems, alanine dipeptide and bovine pancreatic trypsin inhibitor (BPTI), and find that it accurately calculates a variety of equilibrium and dynamical properties. In the case of BPTI, the CPU time required is reduced by nearly a factor of 60 compared to a conventional, unconstrained Langevin simulation using the same force field. Proteins 30:215–227, 1998. © 1998 Wiley-Liss, Inc.

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.

Micelle-forming detergents provide an amphipathic environment that can mimic lipid bilayers and are important tools for solubilizing membrane proteins for functional and structural investigations in vitro. However, the formation of a soluble protein-detergent complex (PDC) currently relies on empirical screening of detergents, and a stable and functional PDC is often not obtained. To provide a foundation for systematic comparisons between the properties of the detergent micelle and the resulting PDC, a comprehensive set of detergents commonly used for membrane protein studies are systematically investigated. Using small-angle X-ray scattering (SAXS), micelle shapes and sizes are determined for phosphocholines with 10, 12, and 14 alkyl carbons, glucosides with 8, 9, and 10 alkyl carbons, maltosides with 8, 10, and 12 alkyl carbons, and lysophosphatidyl glycerols with 14 and 16 alkyl carbons. The SAXS profiles are well described by two-component ellipsoid models, with an electron rich outer shell corresponding to the detergent head groups and a less electron dense hydrophobic core composed of the alkyl chains. The minor axis of the elliptical micelle core from these models is constrained by the length of the alkyl chain, and increases by 1.2–1.5 A° per carbon addition to the alkyl chain. The major elliptical axis also increases with chain length; however, the ellipticity remains approximately constant for each detergent series. In addition, the aggregation number of these detergents increases by ,16 monomers per micelle for each alkyl carbon added. The data provide a comprehensive view of the determinants of micelle shape and size and provide a baseline for correlating micelle properties with protein-detergent interactions.

Freshly-prepared blue membranes from Halobacterium halobium, previously reported to be disordered, are shown to have a distinct crystal lattice structure, slightly different from the native form. The lattice of the blue form is disrupted irreversibly when dehydrated. The disorder process was observed using time-resolved small-angle X-ray diffraction and analyzed by radial autocorrelation functions. The diffraction peaks of the in-plane lattice first sharpen and increase due to improved membrane orientation, then the trimer lattice becomes disordered and the unit cell dimension decreases by 1.8 A. In contrast, dehydration of purple membranes does not disorder the lattice, and the unit cell dimension shrinks by only 1.0 A. Comparisons of radial autocorrelation functions for the blue membrane during drying show drastic loss of inter-trimer, long-range correlation while the intra-trimer, short-range correlations remain more or less unchanged. This suggests that the deionized protein trimers can maintain their overall structure during the dehydration, even though the lattice dimension decreases appreciably and the two-dimensional crystallinity is disrupted.

The nitrogenase enzyme complex, consisting of the molybdenum-iron protein and the iron protein, plays a critical role in the biological reduction of dinitrogen to ammonia (nitrogen fixation). The nitrogen-fixing site within the molybdenum-iron protein is an iron-molybdenum-sulfur cofactor (FeMoco) of roughly 1000-2000 Dalton mass. Structural aspects of FeMoco have been determined by spectroscopic and more recently by crystallographic studies. In order to determine the radius of gyration (Rg) of isolated FeMoco, we have performed small-angle x-ray scattering studies of FeMoco in N-methylformamide solution, in the absence of the molybdenum-iron protein. Model compounds of known structure have also been examined in similar solvents, N,N-dimethylformamide and acetonitrile, as controls and for calibration purposes. The Rg values obtained for the models are in good agreement with calculations based upon their respective crystal structures. However, the Rg obtained for FeMoco clearly indicates that the cofactor is not monomeric in solution, but rather aggregated and possibly polydisperse. Further, Rg values were also measured after addition of thiol, dithionite, and thiol and dithionite, to the FeMoco samples. The results indicate, surprisingly, that oxidation state and putative thiol coordination have no detectable effect on the aggregation behavior of FeMoco in solution, as determined by these measurements.

A Monte Carlo algorithm that searches for the optimal docking configuration of hen egg white lysozyme to an antibody is de- veloped. Both the lysozyme and the antibody are kept rigid. Unlike the work of other au- thors, our algorithm does not attempt to explic- itly maximize surface contact, but minimizes the energy computed using coarse-grained pair potentials. The final refinement of our best so- lutions using all-atom OPLS potentials (Jor- gensen and Tirado-Rives') consistently yields the native conformation as the preferred solu- tion for three different antibodies. We find that the use of an exponential distance-dependent dielectric function is an improvement over the more commonly used linear form. 0 1993 Wiley-Liss, Inc.

We have directly characterized the extent of chain collapse early in the folding of protein L using time-resolved small angle X-ray scattering. We find that, immediately after the initiation of refolding, the protein exhibits dimensions indistinguishable from those observed under highly denaturing, equilibrium conditions and that this expanded initial state collapses with the same rate as that of the overall folding reaction. The observation that chain compaction need not significantly precede the rate-limiting step of folding demonstrates that rapid chain collapse is not an obligatory feature of protein folding reactions.

The ion atmosphere around nucleic acids critically affects biological and physical processes such as chromosome packing, RNA folding, and molecular recognition. However, the dynamic nature of the ion atmosphere renders it difficult to characterize. The basic thermodynamic description of this atmosphere, a full accounting of the type and number of associated ions, has remained elusive. Here we provide the first complete accounting of the ion atmosphere, using buffer equilibration and atomic emission spectroscopy (BE-AES) to accurately quantitate the cation association and anion depletion. We have examined the influence of ion size and charge on ion occupancy around simple, well-defined DNA molecules. The relative affinity of monovalent and divalent cations correlates inversely with their size. Divalent cations associate preferentially over monovalent cations; e.g., with Na+ in four-fold excess of Mg2+ (20 vs. 5 mM), the ion atmosphere nevertheless has three-fold more Mg2+ than Na+. Further, the dicationic polyamine putrescine2+ does not compete effectively for association relative to divalent metal ions, presumably because of its lower charge density. These and other BE-AES results can be used to evaluate and guide the improvement of electrostatic treatments. As a first step, we compare the BE-AES results to predictions from the widely-used nonlinear Poisson Boltzmann (NLPB) theory and assess the applicability and precision of this theory. In the future, BE-AES in conjunction with improved theoretical models, can be applied to complex binding and folding equilibria of nucleic acids and their complexes, to parse the electrostatic contribution from the overall thermodynamics of important biological processes.

Hepatitis C virus NS3 helicase is an enzyme that unwinds double-stranded polynu- cleotides in an ATP-dependent reaction. It pro- vides a promising target for small molecule thera- peutic agents against hepatitis C. Design of such drugs requires a thorough understanding of the dynamical nature of the mechanochemical func- tioning of the helicase. Despite recent progress, the detailed mechanism of the coupling between ATPase activity and helicase activity remains un- clear. Based on an elastic network model (ENM), we apply two computational analysis tools to probe the dynamical mechanism underlying the allosteric coupling between ATP binding and polynucleotide binding in this enzyme. The corre- lation analysis identifies a network of hot-spot residues that dynamically couple the ATP-bind- ing site and the polynucleotide-binding site. Sev- eral of these key residues have been found by mutational experiments as functionally impor- tant, while our analysis also reveals previously unexplored hot-spot residues that are potential targets for future mutational studies. The confor- mational changes between different crystal struc- tures of NS3 helicase are found to be dominated by the lowest frequency mode solved from the ENM. This mode corresponds to a hinge motion of the highly flexible domain 2. This motion simulta- neously modulates the opening/closing of the domains 1-2 cleft where ATP binds, and the domains 2-3 cleft where the polynucleotide binds. Additionally, a small twisting motion of domain 1, observed in both mode 1a nd the computed ATP binding induced conformational change, fine- tunes the binding affinity of the domains 1-3 interface for the polynucleotide. The combination of these motions facilitates the translocation of a single-stranded polynucleotide in an inchworm- like manner. Proteins 2007;67:886-896. V C 2007 Wiley-Liss, Inc. y

We performed a systematic exploration of the use of structural information derived from small angle X-ray scattering (SAXS) measurements to improve fold recognition. SAXS data provide the Fourier transform of the histogram of atomic pair distances (pair distribution function) for a given protein and hence can serve as a structural constraint on methods used to determine the native conformational fold of the protein. Here we used it to construct a similarity-based fitness score with which to evaluate candidate structures generated by a threading procedure. In order to combine the SAXS scores with the standard energy scores and other 1D profile-based scores used in threading, we made use both of a linear regression method and of a neural network-based technique to obtain optimal combined fitness scores and applied them to the ranking of candidate structures. Our results show that the use of SAXS data with gapless threading significantly improves the performance of fold recognition.

Substantial evidence suggests that the aggregation of the presynaptic protein alpha-synuclein is a key step in the etiology of Parkinson's disease (PD). Although the molecular mechanisms underlying alpha-synuclein aggregation remain unknown, oxidative stress has been implicated in the pathogenesis of PD. Here, we report the effects of tyrosine nitration on the propensity of human recombinant alpha-synuclein to fibrillate in vitro. The properties of nitrated alpha-synuclein were investigated using a variety of biophysical and biochemical techniques, which revealed that nitration led to formation of a partially folded conformation with increased secondary structure relative to the intrinsically disordered structure of the monomer, and to oligomerization at neutral pH. The degree of self-association was concentration-dependent, but at 1 mg/mL, nitrated alpha-synuclein was predominantly an octamer. At low pH, small-angle X-ray scattering data indicated that the nitrated protein was monomeric. alpha-Synuclein fibrillation at neutral pH was completely inhibited by nitrotyrosination and is attributed to the formation of stable soluble oligomers. The presence of heparin or metals did not overcome the inhibition; however, the inhibitory effect was eliminated at low pH. The addition of nitrated alpha-synuclein inhibited fibrillation of non-modified alpha-synuclein at neutral pH. Potential implications of these findings to the etiology of Parkinson's disease are discussed.

The structure of C-terminal do- main (CaD136, C-terminal residues 636 -771) of chicken gizzard caldesmon has been analyzed by a variety of physico-chemical methods. We are show- ing here that CaD136 does not have globular struc- ture, has low secondary structure content, is essen- tially noncompact, as it follows from high Rg and RS values, and is characterized by the absence of dis- tinct heat absorption peaks, i.e. it belongs to the family of natively unfolded (or intrinsically unstruc- tured) proteins. Surprisingly, effective binding of single calmodulin molecule (Kd 1.4 0.2 M) leads only to a very moderate folding of this protein and CaD136 remains substantially unfolded within its tight complex with calmodulin. The biological sig- nificance of these observations is discussed. Proteins

Light chain amyloidosis involves the systemic deposition of fibrils in patients overproducing monoclonal immunoglobulin light chains. The kinetics of fibril formation of LEN, a benign light chain variable domain, were investigated at physiological pH in the presence of urea. Despite the lack of in vivofibril formation, LEN readily forms fibrils in vitro under mildly destabilizing conditions. The effect of low to moderate concentrations of urea on the conformation, association state, stability, and kinetics of fibrillation of LEN were investigated. The conformation of LEN was only slightly affected by the addition of up to 4 m urea. The fibrillation kinetics were highly dependent on protein and urea concentrations, becoming faster with decreasing protein concentration and increasing urea concentration. Changes in spectral probes were concomitant to fibril formation throughout the protein and urea concentration ranges, indicating the absence of off-pathway oligomeric species or amorphous aggregates prior to fibril formation. Reducing the amount of dimers initially present in solution by either decreasing the protein concentration or adding urea resulted in faster fibril formation. Thus, increasing concentrations of urea, by triggering dissociation of dimeric LEN, lead to increased rates of fibrillation.

An important element of protein folding theory has been the identification of equilibrium parameters that might uniquely distinguish rapidly folding polypeptide sequences from those that fold slowly. One such parameter, termed sigma, is a dimensionless, equilibrium measure of the coincidence of chain compaction and folding that is predicted to be an important determinant of relative folding kinetics. To test this prediction and improve our understanding of the putative relationship between nonspecific compaction of the unfolded state and protein folding kinetics, we have used small-angle X-ray scattering and circular dichroism spectroscopy to measure the sigma of five well-characterized proteins. Consistent with theoretical predictions, we find that near-perfect coincidence of the unfolded state contraction and folding (sigma approximately 0) is associated with the rapid kinetics of these naturally occurring proteins. We do not, however, observe any significant correlation between sigma and either the relative folding rates of these proteins or the presence or absence of well-populated kinetic intermediates. Thus, while sigma approximately 0 may be a necessary condition to ensure rapid folding, differences in sigma do not account for the wide range of rates and mechanisms with which naturally occurring proteins fold.

In 1977, Z. Kam theorized that correlations of scattering patterns, measured by exposing a solution of randomly oriented identical particles to x-ray radiation, could yield detailed information on the internal structure of the individual particles [1]. During a single exposure (whose duration should be short compared to the particle rotational diffusion time), physical correlations arise whenever multiple photons scatter from the same particle into different directions. By averaging correlations from many exposures, we have demonstrated that one can extract this correlated signal from a background of uncorrelated single-direction scattering events from different particles [2]. This additional information can be used to place constraints on model structures of the particles under investigation, providing a method of structure refinement to atomic resolution. We recently observed correlated scattering from solutions of ~10^9 silver nanoparticles exposed to synchrotron radiation at a microfocus beamline at SSRL [2]. By auto- and cross- correlating the Bragg rings 111 and 200, five correlation peaks were resolved corresponding to the structure and symmetry of silver's reciprocal lattice. To transition from nanoparticles to biomolecule studies, we have performed several experiments at x-ray free electron laser centers (SLAC and SPring-8), and are working to refine analysis techniques.

Human recombinant prothymosin alpha (ProTalpha) is known to have coil-like conformation at neutral pH; i.e., it belongs to the class of "natively unfolded" proteins. By means of circular dichroism, SAXS, and ANS fluorescence, we have investigated the effect of several divalent cations on the structure of this protein. Results of these studies are consistent with the conclusion that ProTalpha conformation is unaffected by large excess of Ca(2+), Mg(2+), Mn(2+), Cu(2+), and Ni(2+). However, Zn(2+) induces compaction and considerable rearrangement of the protein structure. This means that ProTalpha can specifically interact with Zn(2+) (K(D) approximately 10(-3) M), and such interactions induce folding of the natively unfolded protein into a compact partially folded (premolten globule-like) conformation. It is possible that these structural changes may be important for the function of this protein.

The equilibrium unfolding and the kinetic refolding of cytochrome c (Cyt c) in the presence of imidazole were studied with small-angle X-ray scattering (SAXS). The equilibrium unfolding experiments showed the radius of gyration, R(g), of native Cyt c to swell approximately 1 A with the addition of imidazole. The thermodynamic parameter m also reflects an expansion of the protein as its lower value demonstrates an increase in solvent-accessible surface area over that of native Cyt c in the absence of imidazole. Refolding was studied in the presence of imidazole as it prevents misligated intermediate states from forming during the refolding process, simplifying the kinetics, and making them easier to resolve. Time-resolved decreases in the forward scattering amplitude, I(0), demonstrated the transient formation of an aggregated intermediate. Final protein and denaturant concentrations were varied in the refolding kinetics, and the singular value decomposition (SVD) method was employed to characterize the associated state. This state was determined to be a dimer, with properties consistent with a molten globule.

Small angle x-ray scattering experiments have been carried out on the purified iron proteins of nitrogenase from wild-type Azotobacter vinelandii and from a Nif- mutant strain, A. vinelandii UW91 (which has an A157S mutation). This study was designed to investigate the influence of MgATP and MgADP binding on the protein structure in solution. For the wild-type protein, the binding of MgATP induces a significant conformational change that is observed as a decrease of about 2.0 A in the radius of gyration. In contrast, the binding of MgADP to the wild-type iron protein does not detectably affect the radius of gyration. In the absence of nucleotides, the radius of gyration for the UW91 mutant is indistinguishable from that of the wild-type. However, unlike for the wild-type protein, the radius of gyration of the UW91 iron protein is unaffected by the addition of MgATP. We have previously shown that the UW91 iron protein has a normal [4Fe-4S] cluster and MgATP binding ability but that it is completely blocked for electron transfer and MgATP hydrolysis (Gavini, N., and Burgess, B. K. (1992) J. Biol. Chem. 267, 21179-21186). These x-ray scattering measurements suggest that a conformation different from that of the native state is therefore required for the iron protein to perform electron transfer to the MoFe protein. These results also support the hypothesis that Ala-157 is crucial for the iron protein to establish the electron-transfer-favored conformation induced by MgATP binding.

The current work describes a simplified representation of protein structure with uses in the simulation of protein folding. The model assumes that a protein can be represented by a freely rotating rigid chain with a single atom approximating the effect of each side chain. Potentials describing the attraction or repulsion between different types of amino acids are determined directly from the distribution of amino acids in the database of known protein structures. The optimization technique of simulated annealing has been used to dynamically sample the conformations available to this simple model, allowing the protein to evolve from an extended, random coil into a compact globular structure. Many characteristics expected of true proteins, such as the sequence-dependent formation of secondary structure, the partitioning of hydrophobic residues, and specific disulfide pairing, are reproduced by the simulation, suggesting the model may accurately simulate the folding process.

X-ray crystallographic studies of troponin C (Herzberg, O., and James, M.N.G. (1985) Nature 313, 653-659; Sundaralingam, M., Bergstrom, R., Strasburg, G., Rao, S.T., and Roychowdhury, P. (1985a) Science 227, 945-948) have revealed a novel protein structure consisting of two globular domains, each containing two Ca2+-binding sites, connected via a nine-turn alpha-helix, three turns of which are fully exposed to solvent. Since the crystals were grown at pH approximately 5, it is of interest to determine whether this structure is applicable to the protein in solution under physiological conditions. We have used small-angle x-ray scattering to examine the solution structure of troponin C at pH 6.8 and the effect of Ca2+ on the structure. The scattering data are consistent with an elongated structure in solution with a radius of gyration of approximately 23.0 A, which is quite comparable to that computed for the crystal structure. The experimental scattering profile and the scattering profile computed from the crystal structure coordinates do, however, exhibit differences at the 40-A level. A weak Ca2+-facilitated dimerization of troponin C was observed. The data rule out large Ca2+-induced structural changes, indicating rather that the molecule with Ca2+ bound is only slightly more compact than the Ca2+-free molecule.

Linearly polarized synchrotron radiation has been used to investigate the orientation dependence of the K absorption edge of a single crystal of ammonium dithiomolybdate, the crystal structure of which is also reported herein. Following the collection of spectra with the polarization vector parallel to the oxygen-oxygen interatomic vector, the sulfur-sulfur vector, and the twofold rotation axis, self-consistent field X..cap alpha.. multiple scattered wave (SCF X..cap alpha.. MSW) calculations were done to interpret the observed orientation dependence. It is shown that this orientation dependence can be used to highlight the edge features which derive from molecular orbitals contributed to principally by the oxygen of the sulfur ligands. Anisotropy is also apparent in the extended X-ray absorption fine structure (EXAFS) spectra, which were recorded for orientations along the S-S and O-O interatomic vectors. It is thus demonstrated that polarized single-crystal X-ray absorption spectroscopy can provide orientationally selective information about the absorbing atom. It could, in principle, be used to deduce certain types of geometric information about the x-ray absorbing species.

One major obstacle to membrane protein structure determination is the selection of a detergent micelle that mimics the native lipid bilayer. Currently, detergents are selected by exhaustive screening because the effects of protein-detergent interactions on protein structure are poorly understood. In this study, the structure and dynamics of an integral membrane protein in different detergents is investigated by nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) spectroscopy, and small angle X-ray scattering (SAXS). The results suggest that matching of the micelle dimensions to the protein’s hydrophobic surface avoids exchange processes that reduce the completeness of the NMR observations. Based on these dimensions, several mixed micelles were designed that improved the completeness of NMR observations. These findings provide a basis for the rational design of mixed micelles that may advance membrane protein structure determination by NMR.