Peptides as transmembrane segments: decrypting the determinants for helix-helix interactions in membrane proteins

Although the structural analysis of membrane proteins is advancing, an understanding of the basic principles that underlie their folding and assembly remains limited because of the high insolubility intrinsic to these molecules and concomitant challenges in obtaining crystals. Fortunately, from an experimental standpoint, membrane protein folding can be approximated as the rigid-body docking of pre-formed alpha-helical transmembrane segments one with another to form the final functional protein structure. Peptides derived from the sequences of native alpha-helical transmembrane segments and those that mimic their properties are therefore valuable in the experimental evaluation of protein folding within the membrane. Here we present an overview of the progress made in our laboratory and elsewhere in using peptide models toward defining the sequence requirements and forces stabilizing membrane protein folds.     

Physicochemical bases for protein folding,dynamics, and protein-ligand binding

Proteins are essential parts of living organisms and participate in virtually every process within cells.As the genomic sequences for increasing number of organisms are completed,research into how proteins can perform such a variety of functions has become much more intensive because the value of the genomic sequences relies on the accuracy of understanding the encoded gene products.Although the static three-dimensional structures of many proteins are known,the functions of proteins are ultimately governed by their dynamic characteristics,including the folding process,conformational fluctuations,molecular motions,and protein-ligand interactions.In this review,the physicochemical principles underlying these dynamic processes are discussed in depth based on the free energy landscape(FEL)theory.Questions of why and how proteins fold into their native conformational states,why proteins are inherently dynamic,and how their dynamic personalities govern protein functions are answered.This paper will contribute to the understanding of structure-function relationship of proteins in the post-genome era of life science research.     

Simulating protein evolution in sequence and structure space

Naturally occurring proteins comprise a special subset of all plausible sequences and structures selected through evolution. Simulating protein evolution with simplified and all-atom models has shed light on the evolutionary dynamics of protein populations, the nature of evolved sequences and structures, and the extent to which today's proteins are shaped by selection pressures on folding, structure and function. Extensive mapping of the native structure, stability and folding rate in sequence space using lattice proteins has revealed organizational principles of the sequence/structure map important for evolutionary dynamics. Evolutionary simulations with lattice proteins have highlighted the importance of fitness landscapes, evolutionary mechanisms, population dynamics and sequence space entropy in shaping the generic properties of proteins. Finally, evolutionary-like simulations with all-atom models, in particular computational protein design, have helped identify the dominant selection pressures on naturally occurring protein sequences and structures.     

Class-specific correlations between protein folding rate, structure-derived, and sequence-derived descriptors

Small single-domain proteins that fold by simple two-state kinetics have been shown to exhibit a wide variation in their folding rates. It has been proposed that folding mechanisms in these proteins are largely determined by the native-state topology, and a significant correlation between folding rate and measures of the average topological complexity, such as relative contact order (RCO), has been reported. We perform a statistical analysis of folding rate and RCO in all three major structural classes (alpha, beta, and alpha/beta) of small two-state proteins and of RCO in groups of analogous and homologous small single-domain proteins with the same topology. We also study correlation between folding rate and the average physicochemical properties of amino acid sequences in two-state proteins. Our results indicate that 1) helical proteins have statistically distinguishable, class-specific folding rates; 2) RCO accounts for essentially all the variation of folding rate in helical proteins, but for only a part of the variation in beta-sheet-containing proteins; and 3) only a small fraction of the protein topologies studied show a topology-specific RCO. We also report a highly significant correlation between the folding rate and average intrinsic structural propensities of protein sequences. These results suggest that intrinsic structural propensities may be an important determinant of the rate of folding in small two-state proteins.     

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.     

Prediction of the aggregation propensity of proteins from the primary sequence: aggregation properties of proteomes

In the cell, protein folding into stable globular conformations is in competition with aggregation into non-functional and usually toxic structures, since the biophysical properties that promote folding also tend to favor intermolecular contacts, leading to the formation of β-sheet-enriched insoluble assemblies. The formation of protein deposits is linked to at least 20 different human disorders, ranging from dementia to diabetes. Furthermore, protein deposition inside cells represents a major obstacle for the biotechnological production of polypeptides. Importantly, the aggregation behavior of polypeptides appears to be strongly influenced by the intrinsic properties encoded in their sequences and specifically by the presence of selective short regions with high aggregation propensity. This allows computational methods to be used to analyze the aggregation properties of proteins without the previous requirement for structural information. Applications range from the identification of individual amyloidogenic regions in disease-linked polypeptides to the analysis of the aggregation properties of complete proteomes. Herein, we review these theoretical approaches and illustrate how they have become important and useful tools in understanding the molecular mechanisms underlying protein aggregation.     

Extending the concept of template-assembled synthetic proteins.

The creation of native-like macromolecules in copying nature's way represents a fascinating challenge in protein chemistry today. In the absence of a detailed knowledge of the complex folding pathway the ultimate goal in protein de novo design, the construction of artificial proteins with predetermined three-dimensional structure and tailor-made functions based on a defined, generally valid set of rules, appears to be still out of reach. With progress in synthesis strategies and biostructural characterization methods, topological templates have become a versatile tool for inducing and stabilizing secondary and tertiary structures, such as protein loops, beta-turns, alpha-helices, beta-sheets and a variety of folding motifs. In this article, we extend the concept of template-assembled synthetic proteins for the construction of protein-like topologies with multiply bridged, oligocyclic chain architectures termed locked-in tertiary folds that exhibit unique physicochemical and folding properties because of the highly confined conformational space. Furthermore, we show that some fundamental questions in protein assembly can be approached applying the template concept. Using covalent template trapping of self-associated peptide assemblies in aqueous solution the structural and physical forces guiding protein folding, supramolecular assembly and molecular recognition processes can be studied on a molecular level.     

A protein sequence that can encode native structure by disfavoring alternate conformations

The linear sequence of amino acids contains all the necessary information for a protein to fold into its unique three-dimensional structure. Native protein sequences are known to accomplish this by promoting the formation of stable, kinetically accessible structures. Here we describe a Pro residue in the center of the third transmembrane helix of the cystic fibrosis transmembrane conductance regulator that promotes folding by a distinct mechanism: disfavoring the formation of a misfolded structure. The generality of this mechanism is supported by genome-wide transmembrane sequence analyses. Furthermore, the results provide an explanation for the increased frequency of Pro residues in transmembrane alpha-helices. Incorporation by nature of such 'negative folding determinants', aimed at preventing the formation of off-pathway structures, represents an additional mechanism by which folding information is encoded within the evolved sequences of proteins.     

Intrinsically unstructured proteins and their functions

Many gene sequences in eukaryotic genomes encode entire proteins or large segments of proteins that lack a well-structured three-dimensional fold. Disordered regions can be highly conserved between species in both composition and sequence and, contrary to the traditional view that protein function equates with a stable three-dimensional structure, disordered regions are often functional, in ways that we are only beginning to discover. Many disordered segments fold on binding to their biological targets (coupled folding and binding), whereas others constitute flexible linkers that have a role in the assembly of macromolecular arrays.     

Photosynthesis research: advances through molecular biology – the beginnings, 1975–1980s and on...

Restriction endonuclease recognition sites and genes for rRNAs were first mapped on chloroplast chromosomes in 1975–1976. This marked the beginning of the application of molecular biology tools to photosynthesis research. In the first phase, knowledge about proteins involved in photosynthesis was used to identify plastid and nuclear genes encoding these proteins on cloned segments of DNA. Soon afterwards the DNA sequences of the cloned genes revealed the full primary sequences of the proteins. Knowledge of the primary amino acid sequences provided deeper understanding of the functioning of the protein and interactions among proteins of the photosynthetic apparatus. Later, as chloroplast DNA sequencing proceeded, genes were discovered that encoded proteins that had not been known to be part of the photosynthetic apparatus. This more complete knowledge of the composition of reaction centers and of the primary amino acid sequences of individual proteins comprising the reaction centers opened the way to determining the three-dimensional structures of reaction centers. At present, the availability of cloned genes, knowledge of the gene sequences and systems developed to genetically manipulate photosynthetic organisms is permitting experimental inquiries to be made into crucial details of the photosynthetic process. This revised version was published online in August 2006 with corrections to the Cover Date.     

Prediction of Protein Structural Features from Sequence Data Based on Shannon Entropy and Kolmogorov Complexity

While the genome for a given organism stores the information necessary for the organism to function and flourish it is the proteins that are encoded by the genome that perhaps more than anything else characterize the phenotype for that organism. It is therefore not surprising that one of the many approaches to understanding and predicting protein folding and properties has come from genomics and more specifically from multiple sequence alignments. In this work I explore ways in which data derived from sequence alignment data can be used to investigate in a predictive way three different aspects of protein structure: secondary structures, inter-residue contacts and the dynamics of switching between different states of the protein. In particular the use of Kolmogorov complexity has identified a novel pathway towards achieving these goals.     

Amino acid sequence predicts folding rate for middle-size two-state proteins

The significant correlation between protein folding rates and the sequence-predicted secondary structure suggests that folding rates are largely determined by the amino acid sequence. Here, we present a method for predicting the folding rates of proteins from sequences using the intrinsic properties of amino acids, which does not require any information on secondary structure prediction and structural topology. The contribution of residue to the folding rate is expressed by the residue's Omega value. For a given residue, its Omega depends on the amino acid properties (amino acid rigidity and dislike of amino acid for secondary structures). Our investigation achieves 82% correlation with folding rates determined experimentally for simple, two-state proteins studied until the present, suggesting that the amino acid sequence of a protein is an important determinant of the protein-folding rate and mechanism.     

枯草芽孢杆菌同义密码子使用偏性对蛋白质折叠速率的影响

目前,有关同义密码子使用偏性对蛋白质折叠的影响研究中,样本蛋白均来源于不同的物种。考虑到同义密码子使用偏性的物种差异性,选取枯草杆菌的核蛋白为研究对象。首先,将每条核蛋白按二级结构截取为α螺旋片段、β折叠片段和无规卷曲（α-β混合）片段,并计算其蛋白质折叠速率。然后,整理每个片段相应的核酸序列信息,计算其同义密码子使用度。在此基础上,分析枯草芽孢杆菌核蛋白的同义密码子使用偏性与蛋白质折叠速率的相关性。发现对于不同二级结构的肽链片段,都有部分密码子的使用偏性与其对应的肽链折叠速率显著相关。进一步分析发现,与肽链片段折叠速率显著相关的密码子绝大部分为枯草杆菌全序列或核蛋白序列的每一组同义密码子中使用度最高的密码子。结果表明,在蛋白质的折叠过程中,枯草芽孢杆菌的同义密码子使用偏性起着重要作用。     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

Automated search of natively folded protein fragments for high-throughput structure determination in structural genomics

Structural genomic projects envision almost routine protein structure determinations, which are currently imaginable only for small proteins with molecular weights below 25,000 Da. For larger proteins, structural insight can be obtained by breaking them into small segments of amino acid sequences that can fold into native structures, even when isolated from the rest of the protein. Such segments are autonomously folding units (AFU) and have sizes suitable for fast structural analyses. Here, we propose to expand an intuitive procedure often employed for identifying biologically important domains to an automatic method for detecting putative folded protein fragments. The procedure is based on the recognition that large proteins can be regarded as a combination of independent domains conserved among diverse organisms. We thus have developed a program that reorganizes the output of BLAST searches and detects regions with a large number of similar sequences. To automate the detection process, it is reduced to a simple geometrical problem of recognizing rectangular shaped elevations in a graph that plots the number of similar sequences at each residue of a query sequence. We used our program to quantitatively corroborate the premise that segments with conserved sequences correspond to domains that fold into native structures. We applied our program to a test data set composed of 99 amino acid sequences containing 150 segments with structures listed in the Protein Data Bank, and thus known to fold into native structures. Overall, the fragments identified by our program have an almost 50% probability of forming a native structure, and comparable results are observed with sequences containing domain linkers classified in SCOP. Furthermore, we verified that our program identifies AFU in libraries from various organisms, and we found a significant number of AFU candidates for structural analysis, covering an estimated 5 to 20% of the genomic databases. Altogether, these results argue that methods based on sequence similarity can be useful for dissecting large proteins into small autonomously folding domains, and such methods may provide an efficient support to structural genomics projects.     

Identification of structured peptide segments in folding proteins.

Prediction of the structures of long polypeptides and small proteins has been carried out using conformational energy calculations. These calculations can be applied to large proteins if structured regions of their sequences can be identified. Three different approaches to identifying such sequences are presented. First, sequences of five or more contiguous hydrophobic residues tend to nucleate alpha-helices. Second, peptide sequences from parent proteins that have the same biological activities as the parent proteins are highly structured. Third, structured synthetic peptide segments from proteins inhibit the folding of the parent proteins by competing with the corresponding segment of the protein chain for associating with complementary regions. Examples of each of these approaches are presented.     

PROFALIGN algorithm identifies the regions containing folding determinants by scoring pairs of hydrophobic profiles of remotely related proteins.

Profile comparison methods have been shown to be very powerful in creating accurate alignments of protein sequences, especially in the case of remotely related proteins (RRP). These methods take advantage of the observation that hydrophobic profiles are more conserved than the corresponding amino acid sequences. Here, we present the PROFALIGN algorithm, which allows one to perform a detailed comparative analysis, at both local and global levels of two protein sequence profiles. The user can either choose among four different hydrophobic scales (Miyazawa-Jernigan, Eisenberg, Engelman-Steiz, and Kyte-Doolittle) or can add a personal scale. The interface is designed for a wide range of users, including those who are not involved in protein research. It allows one to vary the alignment parameters (such as gap penalties, embedding, and profile smoothness). Secondary structure propensity is added as an optional alignment filter. Similar segments of two proteins are singled out on the basis of score. We have tested the algorithm with different Src homology 3 (SH3) domain fragments sharing low sequence homology but very similar three-dimensional (3D) structures. By using the Miyazawa-Jernigan hydrophobic scale, PROFALIGN was able to detect the strong correlation between the regions that are known to be crucial for SH3 transition state topology. PROFALIGN seems able to identify most of the mutual alignment of structures on the basis of their hydrophobic profiles, delimiting the regions containing the key determinants of folding. Therefore, the present methodology may be useful for the detection of the most structurally relevant positions inside remote related proteins.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

植物凝集素的超级家族

凝集素是一类专一、可逆地和糖类结合的蛋白质,迄今已经分离纯化并测定了氨基酸序列的凝集素已有不少,一些凝集素以及它们与配体糖相互结合的复合物的高级结构也已经给出,许多工作已深入到基因水平．就目前已有的知识,说明植物凝集素是一个庞大的蛋白质家族．