Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

Inherent protein structural flexibility at the RNA-binding interface of L30e

The Saccharomyces cerevisiae ribosomal protein L30 autoregulates its own expression by binding to a purine-rich internal loop in its pre-mRNA and mRNA. NMR studies of L30 and its RNA complex showed that both the internal loop of the RNA as well as a region of the protein become substantially more ordered upon binding. A crystal structure of a maltose binding protein (MBP)-L30 fusion protein with two copies in the asymmetric unit has been determined. The flexible RNA-binding region in the L30 copies has two distinct conformations, one resembles the RNA bound form solved by NMR and the other is unique. Structure prediction algorithms also had difficulty accurately predicting this region, which is consistent with conformational flexibility seen in the NMR and X-ray crystallography studies. Inherent conformational flexibility may be a hallmark of regions involved in intermolecular interactions.     

Identification of a cDNA clone encoding DIP1-binding protein in Drosophila melanogaster

The Drosophila melanogaster L27a gene encodes a ribosomal protein which is a member of the L15 family of ribosomal proteins. D.m. L27a is closely related to the mammalian protein that has been found differentially expressed in lung cancer tissues and therefore could be involved in the control of cell proliferation such as the ribosomal protein S6. Our work elucidates the role of DIP1 which is a novel protein that we found in Drosophila. We performed a two-hybrid system assay and identified the L27a protein as an interactor of DIP1. The interaction was then validated by in vitro binding assays. DIP1, similar to other nuclear proteins in eukaryotes, is localized to the nuclear periphery and chromatin domain in all nuclei, but disappears at the metaphase. It is possible that in D.m. L27a protein, via interaction with DIP1, could be involved in protein synthesis as well as in cell cycle regulation.     

New insights into the interaction of ribosomal protein L1 with RNA

The RNA-binding ability of ribosomal protein L1 is of profound interest, since L1 has a dual function as a ribosomal structural protein that binds rRNA and as a translational repressor that binds its own mRNA. Here, we report the crystal structure at 2.6 A resolution of ribosomal protein L1 from the bacterium Thermus thermophilus in complex with a 38 nt fragment of L1 mRNA from Methanoccocus vannielii. The conformation of RNA-bound T.thermophilus L1 differs dramatically from that of the isolated protein. Analysis of four copies of the L1-mRNA complex in the crystal has shown that domain II of the protein does not contribute to mRNA-specific binding. A detailed comparison of the protein-RNA interactions in the L1-mRNA and L1-rRNA complexes identified amino acid residues of L1 crucial for recognition of its specific targets on the both RNAs. Incorporation of the structure of bacterial L1 into a model of the Escherichia coli ribosome revealed two additional contact regions for L1 on the 23S rRNA that were not identified in previous ribosome models.     

Structural effects of Parkinson's disease linked DJ-1 mutations

Mutations in the protein DJ-1 are associated with familial forms of Parkinson's disease, indicating that DJ-1 may be involved in pathways related to the etiology of this disorder. Here we have used solution state NMR and circular dichroism spectroscopies to evaluate the extent of structural perturbations associated with five different Parkinson's disease linked DJ-1mutations: L166P, E64D, M26I, A104T, and D149A. Comparison of the data with those obtained for the wild-type protein shows that the L166P mutation leads to severe and global destabilization and unfolding of the protein structure, while the structure of the E64D mutation, as expected, is nearly unperturbed. Interestingly, the remaining three mutants all show different degrees of structural perturbation, which are accompanied by a reduction in the thermodynamic stability of the protein. The observed structural and thermodynamic differences are likely to underlie any functional variations between these mutants and the wild type, which in turn are likely responsible for the pathogenicity of these mutations.     

Accurate computer-based design of a new backbone conformation in the second turn of protein L.

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.     

Structural compliance: A new metric for protein flexibility

Proteins are the active players in performing essential molecular activities throughout biology, and their dynamics has been broadly demonstrated to relate to their mechanisms. The intrinsic fluctuations have often been used to represent their dynamics and then compared to the experimental B-factors. However, proteins do not move in a vacuum and their motions are modulated by solvent that can impose forces on the structure. In this paper, we introduce a new structural concept, which has been called the structural compliance, for the evaluation of the global and local deformability of the protein structure in response to intramolecular and solvent forces. Based on the application of pairwise pulling forces to a protein elastic network, this structural quantity has been computed and sometimes is even found to yield an improved correlation with the experimental B-factors, meaning that it may serve as a better metric for protein flexibility. The inverse of structural compliance, namely the structural stiffness, has also been defined, which shows a clear anticorrelation with the experimental data. Although the present applications are made to proteins, this approach can also be applied to other biomolecular structures such as RNA. This present study considers only elastic network models, but the approach could be applied further to conventional atomic molecular dynamics. Compliance is found to have a slightly better agreement with the experimental B-factors, perhaps reflecting its bias toward the effects of local perturbations, in contrast to mean square fluctuations. The code for calculating protein compliance and stiffness is freely accessible at https://jerniganlab.github.io/Software/PACKMAN/Tutorials/compliance .     

A database of protein structure families with common folding motifs.

The availability of fast and robust algorithms for protein structure comparison provides an opportunity to produce a database of three-dimensional comparisons, called families of structurally similar proteins (FSSP). The database currently contains an extended structural family for each of 154 representative (below 30% sequence identity) protein chains. Each data set contains: the search structure; all its relatives with 70-30% sequence identity, aligned structurally; and all other proteins from the representative set that contain substructures significantly similar to the search structure. Very close relatives (above 70% sequence identity) rarely have significant structural differences and are excluded. The alignments of remote relatives are the result of pairwise all-against-all structural comparisons in the set of 154 representative protein chains. The comparisons were carried out with each of three novel automatic algorithms that cover different aspects of protein structure similarity. The user of the database has the choice between strict rigid-body comparisons and comparisons that take into account interdomain motion or geometrical distortions; and, between comparisons that require strictly sequential ordering of segments and comparisons, which allow altered topology of loop connections or chain reversals. The data sets report the structurally equivalent residues in the form of a multiple alignment and as a list of matching fragments to facilitate inspection by three-dimensional graphics. If substructures are ignored, the result is a database of structure alignments of full-length proteins, including those in the twilight zone of sequence similarity.(ABSTRACT TRUNCATED AT 250 WORDS)     

Analysis of the Structure of the PsbO Protein and its Implications

The PsbO protein is a ubiquitous extrinsic subunit of Photosystem II (PS II), the water splitting enzyme of photosynthesis. A recently determined 3D X-ray structure of a cyanobacterial protein bound to PS II has given an opportunity to conduct complete analyses of its sequence and structural characteristics using bioinformatic methods. Multiple sequence alignments for the PsbO family are constructed and correlated with the cyanobacterial structure. We identify the most conserved regions of PsbO and the mapping of their positions within the structure indicates their functional roles especially in relation to interactions of this protein with the lumenal surface of PS II. Homologous models for eukaryotic PsbO were built in order to compare with the prokaryotic protein. We also explore structural homology between PsbO and other proteins for which 3D structures are known and determine its structural classification. These analyses contribute to the understanding of the function and evolutionary origin of the PS II manganese stabilising protein.     

A ribosomal protein from Thermus thermophilus is homologous to a general shock protein

The gene encoding the ribosomal protein from Thermus thermophilus, TL5, which binds to the 5S rRNA, has been cloned and sequenced. The codon usage shows a clear preference for G/C rich codons that is characteristic for many genes in thermophilic bacteria. The deduced amino acid sequence consists of 206 residues. The sequence of TL5 shows a strong similarity to a general shock protein from Bacillus subtilis, named CTC. The protein CTC is homologous in its N-terminal part to the 5S rRNA binding protein, L25, from E coli. An alignment of the TL5, CTC and L25 sequences displays a number of residues that are totally conserved. No clear sequence similarity was found between TL5 and other proteins which are known to bind to 5S rRNA. The evolutionary relationship of a heat shock protein in mesophiles and a ribosomal protein in thermophilic bacteria as well as a possible role of TL5 in the ribosome are discussed.     

A molecular dynamics study of the C-terminal fragment of the L7/L12 ribosomal protein

The crystallographic dimer of the C-terminal fragment (CTF) of the L7/L12 ribosomal protein has been subjected to molecular dynamics (MD) simulations. A 90 picosecond (ps) trajectory for the protein dimer, 19 water molecules and two counter ions has been calculated at constant temperature. Effects of intermolecular interactions on the structure and dynamics have been studied. The exact crystallographic symmetry is lost and the atomic fluctuations differ from one monomer to the other. The average MD structure is more stable than the X-ray one, as judged by accessible surface area and energy calculations. Crystal (non-dimeric) interactions have been simulated in another 40 ps trajectory by using harmonic restraints to represent intermolecular hydrogen bonds. The conformational changes with respect ot the X-ray structure are then virtually suppressed.The unrestrained dimer trajectory has been scanned for cooperative motions involving secondary structure elements. The intrinsic collective motions of the monomer are transmitted via intermolecular contacts to the dimer structure.The existence of a stable dimeric form of CTF, resembling the crystallographic one, has been documented. At the cost of fairly small energy expenditure the dimer has considerable conformational flexibility. This flexibility may endow the dimer with some functional potential as an energy transducer.     

Conformational analysis of peptides corresponding to all the secondary structure elements of protein L B1 domain: secondary structure propensities are not conserved in proteins with the same fold.

The solution conformation of three peptides corresponding to the two beta-hairpins and the alpha-helix of the protein L B1 domain have been analyzed by circular dichroism (CD) and nuclear magnetic resonance spectroscopy (NMR). In aqueous solution, the three peptides show low populations of native and non-native locally folded structures, but no well-defined hairpin or helix structures are formed. In 30% aqueous trifluoroethanol (TFE), the peptide corresponding to the alpha-helix adopts a high populated helical conformation three residues longer than in the protein. The hairpin peptides aggregate in TFE, and no significant conformational change occurs in the NMR observable fraction of molecules. These results indicate that the helical peptide has a significant intrinsic tendency to adopt its native structure and that the hairpin sequences seem to be selected as non-helical. This suggests that these sequences favor the structure finally attained in the protein, but the contribution of the local interactions alone is not enough to drive the formation of a detectable population of native secondary structures. This pattern of secondary structure tendencies is different to those observed in two structurally related proteins: ubiquitin and the protein G B1 domain. The only common feature is a certain propensity of the helical segments to form the native structure. These results indicate that for a protein to fold, there is no need for large native-like secondary structure propensities, although a minimum tendency to avoid non-native structures and to favor native ones could be required.     

A study of the membrane-water interface region of membrane proteins

The most conspicuous structural characteristic of the alpha-helical membrane proteins is their long transmembrane alpha-helices. However, other structural elements, as yet largely ignored in statistical studies of membrane protein structure, are found in those parts of the protein that are located in the membrane-water interface region. Here, we show that this region is enriched in irregular structure and in interfacial helices running roughly parallel with the membrane surface, while beta-strands are extremely rare. The average amino acid composition is different between the interfacial helices, the parts of the transmembrane helices located in the interface region, and the irregular structures. In this region, hydrophobic and aromatic residues tend to point toward the membrane and charged/polar residues tend to point away from the membrane. The interface region thus imposes different constraints on protein structure than do the central hydrocarbon core of the membrane and the surrounding aqueous phase.     

A rice (Oryza sativa L.) cDNA encodes a protein sequence homologous to the eukaryotic ribosomal 5S RNA-binding protein

A rice (Oryza sativa L.) cDNA clone coding for the cytoplasmic ribosomal protein L5, which associates with 5 S rRNA for ribosome assembly, was cloned and its nucleotide sequence was determined. The primary structure of rice L5, deduced from the nucleotide sequence, contains 294 amino acids and has intriguing features some of which are also conserved in other eucaryotic homologues. These include: four clusters of basic amino acids, one of which may serve as a nucleolar localization signal; three repeated amino acid sequences; the conservation of glycine residues. This protein was identified as the nuclear-encoded cytoplasmic ribosomal protein L5 of rice by sequence similarity to other eucaryotic ribosomal 5 S RNA-binding proteins of rat, chicken, Xenopus laevis, and Saccharomyces cerevisiae. Rice L5 shares 51 to 62% amino acid sequence identity with the homologues. A group of ribosomal proteins from archaebacteria including Methanococcus vanniellii L18 and Halobacterium cutirubrum L13, which are known to be associated with 5 S rRNA, also related to rice L5 and the other eucaryotic counterparts, suggesting an evolutionary relationship in these ribosomal 5 S RNA-binding proteins.     

The solution structure of ribosomal protein L18 from Bacillus stearothermophilus

A medium resolution solution structure has been obtained for L18 from Bacillus stearothermophilus (BstL18), a ribosomal protein that stabilizes the tertiary structure of 5S rRNA and mediates its interaction with the rest of the large subunit. The N-terminal 22 amino acid residues of BstL18 are unstructured in solution. Its remaining 98 residues form a globular domain that has the same topology as the globular domains of other L18s, but the orientation of helices is different. This conformational peculiarity should not prevent BstL18 from functioning in the ribosome the same way as other L18s.     

Intermediates and the folding of proteins L and G

We use a minimalist protein model, in combination with a sequence design strategy, to determine differences in primary structure for proteins L and G, which are responsible for the two proteins folding through distinctly different folding mechanisms. We find that the folding of proteins L and G are consistent with a nucleation-condensation mechanism, each of which is described as helix-assisted beta-1 and beta-2 hairpin formation, respectively. We determine that the model for protein G exhibits an early intermediate that precedes the rate-limiting barrier of folding, and which draws together misaligned secondary structure elements that are stabilized by hydrophobic core contacts involving the third beta-strand, and presages the later transition state in which the correct strand alignment of these same secondary structure elements is restored. Finally, the validity of the targeted intermediate ensemble for protein G was analyzed by fitting the kinetic data to a two-step first-order reversible reaction, proving that protein G folding involves an on-pathway early intermediate, and should be populated and therefore observable by experiment.     

Structural comparison of yeast snoRNP and spliceosomal protein Snu13p with its homologs

Snu13p is a bifunctional yeast protein involved in both messenger RNA splicing as well as ribosomal RNA maturation. Snu13p initiates assembly of ribonucleoprotein particles by interacting with a conserved RNA motif called kink turn. Unlike its archaeal homolog, L7Ae, Snu13p displays differential specificity for functionally distinct kink turns. Thus, the structures of Snu13p at different functional states, including those alone and bound with RNAs, are required to understand how the protein differentially interacts with kink turns. Although the structure of the human homolog of Snu13p bound with a spliceosomal RNA is known, there has not been a report of a structure of free Snu13p. This has hindered our ability to understand the structural basis for Snu13p's substrate specificity. We report a crystal structure of free Snu13p at 1.9A and a detailed structural comparison with its homologs. We show that free Snu13p has nearly an identical conformation as that of its human homolog bound with RNA. Interestingly, both eukaryotic proteins exhibit notable structural differences in their central beta-sheets as compared to their archaeal homolog, L7Ae. The observed structural differences offer a possible explanation to the observed difference in RNA specificity between Snu13p and L7Ae.