Linderstrom-Lang-Schellman's model for protein stabilization revisited

The fact that cleavage of single peptide linkages in proteins often leads to extensive conformational alteration, including regions far removed from the cleavage site is not fully understood. We propose, based on the work of Linderstrom-Lang and Schellman, that disruption primarily occurs within protein structural domains that are stabilized by cooperative interactions and that cleavage of single peptide linkages of the domain perturbs the entire cooperative interaction. For this model we review experimental observations: on fragment complexation (ribonuclease A, staphylococcal nuclease and cytochrome c), destabilized N-terminal large fragments (ribonuclease A and nuclease), cooperative folding and stabilization of proteins (ribonuclease A, nuclease and cytochrome c), the close relationship of the three-dimensional structure between fragment complexes and the original protein (ribonuclease A and nuclease), ligand induced stabilization (nuclease), 3D domain swapping, circular permutation (dihydrofolate reductase), evolutionary conservation (cytochrome c fold). Based on analysis of these observations, we conclude that the cooperative interactions of domains are important for the mechanism of 3D domain swapping as well as for stabilization and thereby, determination of the ground state of native proteins. Furthermore, analysis of the observations reveals that domains generally contain a hydrophobic core. Further, based on studies of cytochrome c and the Tsao, Evans and Wennerstrom model of electrostatic interactions between two hydrophobic monolayers, we propose the model that the hydrophobic core of a domain is polarizable and responds to the surface charges through its polarizability to stabilize the domain, explaining in part the nature of the cooperative interactions.     

Domain swapping proceeds via complete unfolding: a 19F- and 1H-NMR study of the Cyanovirin-N protein

Domain swapping creates protein oligomers by exchange of structural units between identical monomers. At present, no unifying molecular mechanism of domain swapping has emerged. Here we used the protein Cyanovirin-N (CV-N) and (19)F-NMR to investigate the process of domain swapping. CV-N is an HIV inactivating protein that can exist as a monomer or a domain-swapped dimer. We measured thermodynamic and kinetic parameters of the conversion process and determined the size of the energy barrier between the two species. The barrier is very large and of similar magnitude to that for equilibrium unfolding of the protein. Therefore, for CV-N, overall unfolding of the polypeptide is required for domain swapping.     

Domains in bovine seminal ribonuclease

Bovine seminal ribonuclease is the only pancreatic-type ribonuclease to possess a dimeric structure: the two identical subunits are covalently linked by two disulfide bridges. Actually, the protein exists in two different dimeric structures owing to the possibility of swapping the N-terminal α-helical segments: the swapped MxM dimer, and the non-swapped M=M dimer. The thermal denaturation of the two separated forms is investigated by differential scanning calorimetry. The process is reversible and can be represented by two sequential two-state transitions, indicating the presence of two domains in BS-RNase, regardless of the swapping phenomenon. Inspection of the structural models leads to the tentative identification of an external domain and a core domain, the latter more stable.     

Why not consider a spherical protein? Implications of backbone hydrogen bonding for protein structure and function

The intrinsic ability of protein structures to exhibit the geometric features required for molecular function in the absence of evolution is examined in the context of three systems: the reference set of real, single domain protein structures, a library of computationally generated, compact homopolypeptides, artificial structures with protein-like secondary structural elements, and quasi-spherical random proteins packed at the same density as proteins but lacking backbone secondary structure and hydrogen bonding. Without any evolutionary selection, the library of artificial structures has similar backbone hydrogen bonding, global shape, surface to volume ratio and statistically significant structural matches to real protein global structures. Moreover, these artificial structures have native like ligand binding cavities, and a tiny subset has interfacial geometries consistent with native-like protein-protein interactions and DNA binding. In contrast, the quasi-spherical random proteins, being devoid of secondary structure, have a lower surface to volume ratio and lack ligand binding pockets and intermolecular interaction interfaces. Surprisingly, these quasi-spherical random proteins exhibit protein like distributions of virtual bond angles and almost all have a statistically significant structural match to real protein structures. This implies that it is local chain stiffness, even without backbone hydrogen bonding, and compactness that give rise to the likely completeness of the library solved single domain protein structures. These studies also suggest that the packing of secondary structural elements generates the requisite geometry for intermolecular binding. Thus, backbone hydrogen bonding plays an important role not only in protein structure but also in protein function. Such ability to bind biological molecules is an inherent feature of protein structure; if combined with appropriate protein sequences, it could provide the non-zero background probability for low-level function that evolution requires for selection to occur.     

Noncovalent cross-links in context with other structural and functional elements of proteins

Proteins are heteropolymers with evolutionary selected native sequences of residues. These native sequences code for unique and stable 3D structures indispensable for biochemical activity and for proteolysis resistance, the latter which guarantees an appropriate lifetime for the protein in the protease rich cellular environment. Cross-links between residues close in space but far in the primary structure are required to maintain the folded structure of proteins. Some of these cross-links are covalent, most frequently disulfide bonds, but the majority of the cross-links are sets of cooperative noncovalent long-range interactions. In this paper we focus on special clusters of noncovalent long-range interactions: the Stabilization Centers (SCs). The relation between the SCs and secondary structural elements as well as the relation between SCs and functionally important regions of proteins are presented to show a detailed picture of these clusters, which are believed to be primarily responsible for major aspects of protein stability.     

Identification of Regions in Apomyoglobin that Form Intermolecular Interactions in Amyloid Aggregates Using High-Performance Mass Spectrometry

The formation of amyloid aggregates in human organs and tissues causes the development of incurable diseases. However, experimental studies of the mechanism of amyloid formation by proteins and the structural characteristics of amyloids are complicated because of the heterogeneity and high molecular weight of the aggregates. We used limited proteolysis and mass spectrometry for the identification of regions in the apomyoglobin polypeptide chain, which give rise to intermolecular interactions in amyloid structures. Tandem mass spectroscopy enabled the identification of regions in the myoglobin polypeptide chain, which form the core of amyloid structures. It was shown that the main structural elements for the formation of the core of amyloid fibrils in myoglobin were regions from 60 through 90 and from 97 through 124 amino acid residues. These regions coincide well with those theoretically predicted. This approach yielded important data on the structure of protein molecules in aggregates and on conformational rearrangements of apomyoglobin upon amyloid formation.     

Interaction of Adjacent Amino Acids

Ramachandran plots display the dihedral angles of a single protein residue. We propose a crossed torsion angle plot called SSY‐plot between two neighboring amino acids and demonstrate that a special coherence motion can exist between some very special amino acid pairs leading to spontaneous unusual structures. A 6mer was extracted from a BBA polypeptide chain which in this plot shows two diagonal domains for the Ser‐Arg pair after some induction time. Other amino acid pairs in general do not show this kind of split domain. This shows that a special pair is required for stabilizing two distinct native structures in protein folding. We suggest that the existence of these two domains corresponds to a bifurcation between two different protein structures and that the special pair is the key to producing these two structures. These two different structures are produced spontaneously without an external agent.     

The structural determinants that lead to the formation of particular oligomeric structures in the pancreatic-type ribonuclease family

Pancreatic-type ribonucleases are a family of RNA degrading enzymes that share different degrees of sequence identity but a very similar 3D-structure. The prototype of this family is bovine pancreatic ribonuclease or ribonuclease A. This enzyme has been the object of landmark work on the folding, stability, protein chemistry, catalysis, enzyme-substrate interaction and molecular evolution. In the recent years, the interest in the study of pancreatic-type ribonucleases has increased due to the involvement of some members of this family in special biological functions. In addition, dimeric and also higher oligomeric structures can be attained by the members of this family. The oligomers described structurally to date are mainly formed by 3D-domain swapping, a process which consists of the exchange of identical domains (i.e. identical structural elements, usually the N- and C-termini) between the subunits and is considered to be a mechanism for amyloid-type aggregate formation. This review compares the dimeric and oligomeric structures of different members of the pancreatic-type ribonuclease family which are able to acquire these structures, namely, bovine seminal ribonuclease, ribonuclease A and its human counterpart, human pancreatic ribonuclease. A specific focus is placed on what is known about the structural determinants that lead to the acquisition of a particular oligomeric structure and on the proposed mechanism of 3D-swapping.     

Synthetic non-peptide mimetics of alpha-helices

Proteins in nature fold into native conformations in which combinations of peripherally projected aliphatic, aromatic and ionic functionalities direct a wide range of properties. Alpha-helices, one of the most common protein secondary structures, serve as important recognition regions on protein surfaces for numerous protein-protein, protein-DNA and protein-RNA interactions. These interactions are characterized by conserved structural features within the alpha-helical domain. Rational design of structural mimetics of these domains with synthetic small molecules has proven an effective means to modulate such protein functions. In this tutorial review we discuss strategies that utilize synthetic small-molecule antagonists to selectively target essential protein-protein interactions involved in certain diseases. We also evaluate some of the protein-protein interactions that have been or are potential targets for alpha-helix mimetics.     

Trigger factor slows co-translational folding through kinetic trapping while sterically protecting the nascent chain from aberrant cytosolic interactions

The E. coli chaperone trigger factor (TF) interacts directly with nascent polypeptide chains as they emerge from the ribosome exit tunnel. Small protein domains can fold under the cradle created by TF, but the co-translational folding of larger proteins is slowed down by its presence. Because of the great experimental challenges in achieving high spatial and time resolution, it is not yet known whether or not TF alters the folding properties of small proteins and if the reduced rate of folding of larger proteins is the result of kinetic or thermodynamic effects. We show, by molecular simulations employing a coarse-grained model of a series of ribosome nascent-chain complexes, that TF does not alter significantly the co-translational folding process of a small protein G domain but delays that of a large β-galactosidase domain as a result of kinetic trapping of its unfolded ensemble. We demonstrate that this trapping occurs through a combination of three distinct mechanisms: a decrease in the rate of structural rearrangements within the nascent chain, an increase in the effective exit tunnel length due to folding outside the cradle, and entanglement of the nascent chain with TF. We present evidence that this TF-induced trapping represents a trade-off between promoting co-translational folding and sterically shielding the nascent chain from aberrant cytosolic interactions that could lead to its aggregation or degradation.     

Targeting metastable coiled-coil domains by computational design

Approximately 30% of eukaryotic genomes are predicted to encode partially unfolded proteins. Many of these unstructured domains contact multiple partners in short-lived interactions critical for cellular homeostasis. Understanding the functional implications of these transient binding events is a current challenge that could be addressed with designed peptide inhibitors. Most current protein design methodologies, however, target only structurally well-defined, stable structures. To address this limitation, we implemented a computational design strategy that alternates between a fixed backbone sequence search for binding specificity and structural optimization of the designed interfaces. We applied this method to create specific peptide inhibitors of the C-terminal metastable coiled-coil domain of the essential yeast septin Cdc12p. Specific binding of the designed sequences was demonstrated by circular dichroism and equilibrium ultracentrifugation. Our results validate computational methods to design specific peptide ligands to protein domains lacking intrinsic structural stability and set the stage for functional analysis of Cdc12p coiled coil function in vivo.     

Mechanisms of molecular recognition in the pikromycin polyketide synthase

BACKGROUND: Modular polyketide synthases (PKSs) produce a wide range of medically significant compounds. In the case of the pikromycin PKS of Streptomyces venezuelae, four separate polypeptides (PikAI-PikAIV), comprising a total of one loading domain and six extension modules, generate the 14-membered ring macrolactone narbonolide. The polypeptide PikAIV contains a thioesterase (TE) domain and is responsible for catalyzing both the last elongation step with methylmalonyl CoA, and subsequent release of the final polyketide chain elongation intermediate from the PKS. Under certain growth conditions this polypeptide is synthesized from an alternative translational start site, giving rise to an N-terminal truncated form of PikAIV, containing only half of the ketosynthase (KS(6)) domain. The truncated form of PikAIV is unable to catalyze the final elongation step, but is able to cleave a polyketide chain from the preceding module on PikAIII (ACP(5)), giving rise to the 12-membered ring product 10-deoxymethynolide. RESULTS: S. venezuelae mutants expressing hybrid PikAIV polypeptides containing acyl carrier protein (ACP) and malonyl CoA specific acyltransferase (AT) domains from the rapamycin PKS were unable to catalyze production of 12- or 14-membered ring macrolactone products. Plasmid-based expression of a hybrid PikAIV containing the native KS(6) and TE domains, however, restored production of both narbonolide and 10-deoxymethynolide in the S. venezuelae AX912 mutant that generates a TE-deleted form of PikAIV. Use of alternative KS domains or deletion of the KS(6) domain within the hybrid PikAIV resulted in loss of both products. Plasmid-based expression of a TE domain of PikAIV as a separate polypeptide in the AX912 mutant resulted in greater than 50% restoration of 10-deoxymethynolide, but not in mutants expressing a hybrid PikAIV bearing an unnatural AT domain. Mutants expressing hybrid PikAIV polypeptides containing the natural AT(6) domains and different ACP domains efficiently produced polyketide products, but with a significantly higher 10-deoxymethynolide/narbonolide ratio than observed with native PikAIV. CONCLUSIONS: Dimerization of KS(6) modules allows in vivo formation of a PKS heterodimer using PikAIV polypeptides containing different AT and ACP domains. In such heterodimers, the TE domain and the AT(6) domain responsible for formation of the narbonolide product are located on different polypeptide chains. The AT(6) domain of PikAIV plays an important role in facilitating TE-catalyzed chain termination (10-deoxymethynolide formation) at the proceeding module in PikAIII. The pikromycin PKS can tolerate the presence of multiple forms (active and inactive) of PikAIV, and decreased efficiency of elongation by PikAIV can result in increased levels of 10-deoxymethynolide. These results provide new insight into functional molecular interactions and interdomain recognition in modular PKSs.     

Engineering of an "unnatural" natural product by swapping polyketide synthase domains in Aspergillus nidulans

An StcA-AfoE hybrid polyketide synthase (PKS), generated by swapping the AfoE (asperfuranone biosynthesis) SAT domain with the StcA (sterigmatocystin biosynthesis) SAT domian, produced a major new metabolite with the same chain length as the native AfoE product. Structure elucidation allowed us to propose a likely pathway, and feeding studies supported the hypothesis that the chain length of PKS metabolites may be under precise control of KS and PT domains.     

Simultaneous prediction of binding free energy and specificity for PDZ domain–peptide interactions

Interactions between protein domains and linear peptides underlie many biological processes. Among these interactions, the recognition of C-terminal peptides by PDZ domains is one of the most ubiquitous. In this work, we present a mathematical model for PDZ domain–peptide interactions capable of predicting both affinity and specificity of binding based on X-ray crystal structures and comparative modeling with Rosetta. We developed our mathematical model using a large phage display dataset describing binding specificity for a wild type PDZ domain and 91 single mutants, as well as binding affinity data for a wild type PDZ domain binding to 28 different peptides. Structural refinement was carried out through several Rosetta protocols, the most accurate of which included flexible peptide docking and several iterations of side chain repacking and backbone minimization. Our findings emphasize the importance of backbone flexibility and the energetic contributions of side chain-side chain hydrogen bonds in accurately predicting interactions. We also determined that predicting PDZ domain–peptide interactions became increasingly challenging as the length of the peptide increased in the N-terminal direction. In the training dataset, predicted binding energies correlated with those derived through calorimetry and specificity switches introduced through single mutations at interface positions were recapitulated. In independent tests, our best performing protocol was capable of predicting dissociation constants well within one order of magnitude of the experimental values and specificity profiles at the level of accuracy of previous studies. To our knowledge, this approach represents the first integrated protocol for predicting both affinity and specificity for PDZ domain–peptide interactions.     

Progress in protein structural class prediction and its impact to bioinformatics and proteomics

The structural class is an important attribute used to characterize the overall folding type of a protein or its domain. Since the concept of protein structural class was developed about 3 decades ago based on a visual inspection of polypeptide chain topologies in a dataset of only 31 gloular proteins, the number of structure-known proteins has been increased rapidly. For example, as of 12-July-2005, the entries deposited into RCSB PDB Protein Data Bank for proteins, peptides, and viruses whose 3-dimensional structures were determined by X-ray and NMR techniques have been increased to 28,920. To properly cover more and more structure-known proteins, some modification and expansion from the original structural classification scheme have been developed. Meanwhile, many different approaches have been proposed for predicting the structural class of proteins. In this review, the new classification schemes are briefly introduced. The attention is focused on the progress in structural class prediction and its impact in stimulating the development of identifying the other attributes of proteins. It is interesting to point out that the development of the latter has actually in turn greatly enriched the power of the former. Also, some promising approaches for the further development of protein structural class prediction are also addressed.     

Mechanical properties of the domains of titin in a Go-like model

Comparison of properties of three domains of titin, I1, I27, and I28, in a simple geometry-based model shows that despite a high structural homology between their native states different domains show similar but distinguishable mechanical properties. Folding properties of the separate domains are predicted to be diversified which reflects sensitivity of the kinetics to the details of native structures. The Go-like model corresponding to the experimentally resolved native structure of the I1 domain is found to provide the biggest thermodynamic and mechanical stability compared to the other domains studied here. We analyze elastic, thermodynamic, and kinetic properties of several structures corresponding to the I28 domain as obtained through homology-based modeling. We discuss the ability of the models of the I28 domain to reproduce experimental results qualitatively. A strengthening of contacts that involve hydrophobic amino acids does not affect theoretical comparisons of the domains. Tandem linkages of up to five identical or different domains unravel in a serial fashion at low temperatures. We study the nature of the intermediate state that arises in the early stages of the serial unraveling and find it to qualitatively agree with the results of Marszalek et al.     

Differential ordering of the protein backbone and side chains during protein folding revealed by site-specific recombinant infrared probes

The time scale for ordering of the polypeptide backbone relative to the side chains is a critical issue in protein folding. The interplay between ordering of the backbone and ordering of the side chains is particularly important for the formation of β-sheet structures, as the polypeptide chain searches for the native stabilizing cross-strand interactions. We have studied these issues in the N-terminal domain of protein L9 (NTL9), a model protein with mixed α/β structure. We have developed a general approach for introducing site-specific IR probes for the side chains (azide) and backbone ((13)C═(18)O) using recombinant protein expression. Temperature-jump time-resolved IR spectroscopy combined with site-specific labeling enables independent measurement of the respective backbone and side-chain dynamics with single residue resolution. We have found that side-chain ordering in a key region of the β-sheet structure occurs on a slower time scale than ordering of the backbone during the folding of NTL9, likely as a result of the transient formation of non-native side-chain interactions.     

Redesign of protein domains using one-bead-one-compound combinatorial chemistry

A novel combinatorial strategy for the redesign of proteins based on the strength and specificity of intra- and interprotein interactions is described. The strategy has been used to redesign the hydrophobic core of the B domain of protein A. Using one-bead-one-compound combinatorial chemistry, 300 analogues of the C-terminal alpha-helix of the B domain, H3x, have been synthesized using a biocompatible resin and the HMFS linker, allowing the screening to occur while the peptides were bound to the resin. The screening was based on the ability of the H3x analogues to interact with the N-terminal helices of the B domain, H1-H2, and retain the native B domain activity, that is binding to IgG. Eight different analogues containing some nonconservative mutations were obtained from the library, the two most frequent of which, H3P1 and H3P2, were studied in detail. CD analysis revealed that the active analogues interact with H1-H2. To validate the redesign strategy the covalent modified domains H1-H2-H3P1 and H1-H2-H3P2 were synthesized and characterized. CD and NMR analysis revealed that they had a unique, stable, and well-defined three-dimensional structure similar to that for the wild-type B domain. This combinatorial strategy allows us to select for redesigned proteins with the desired activity or the desired physicochemical properties provided the right screening test is used. Furthermore, it is rich in potential for the chemical modification of proteins overcoming the drawbacks associated with the total synthesis of large protein domains.     

Improved greedy algorithm for protein structure reconstruction

This article concerns the development of an improved greedy algorithm for protein structure reconstruction. Our stochastic greedy algorithm, which attempts to locate the ground state of an approximate energy function, exploits the fact that protein structures consist of overlapping structural building blocks that are not independent. Application of this approach to a series of 16 proteins with 50-250 amino acids leads to predicted models deviating from the experimental structures by 0.5 A RMSD using an RMSD-based energy function and within 1.5 to 4.8 A RMSD using a Go-based energy function. The Go-based results are significant because they illustrate the strength of combining structural fragments and stochastic greedy algorithms in capturing the native structures of proteins stabilized by long-range interactions separated by more than 30 amino acids. These results clearly open the door to less computationally demanding solutions to predict structures from sequences.