Secondary structure prediction and folding of globular protein: refolding of ferredoxin

The physicochemical mechanism of protein folding has been elucidated by the island model, describing a growth type of folding. The folding pathway is closely related with nucleation on the polypeptide chain and thus the formation of small local structures or secondary structures at the earliest stage of folding is essential to all following steps. The island model is applicable to any protein, but a high precision of secondary structure prediction is indispensable to folding simulation. The secondary structures formed at the earliest stage of folding are supposed to be of standard form, but they are usually deformed during the folding process, especially at the last stage, although the degree of deformation is different for each protein. Ferredoxin is an example of a protein having this property. According to X-ray investigation (1FDX), ferredoxin is not supposed to have secondary structures. However, if we assumed that in ferredoxin all the residues are in a coil state, we could not attain the correct structure similar to the native one. Further, we found that some parts of the chain are not flexible, suggesting the presence of secondary structures, in agreement with the recent PDB data (1DUR). Assuming standard secondary structures (-helices and -strands) at the nonflexible parts at the early stage of folding, and deforming these at the final stage, a structure similar to the native one was obtained. Another peculiarity of ferredoxin is the absence of disulfide bonds, in spite of its having eight cysteines. The reason cysteines do not form disulfide bonds became clear by applying the lampshade criterion, but more importantly, the two groups of cysteines are ready to make iron complexes, respectively, at a rather later stage of folding. The reason for poor prediction accuracy of secondary structure with conventional methods is discussed.     

Physical principles of protein structure and protein folding

Physical principles determining the protein structure and protein folding are reviewed: (i) the molecular theory of protein secondary structure and the method of its prediction based on this theory; (ii) the existence of a limited set of thermodynamically favourable folding patterns of α- and β-regions in a compact globule which does not depend on the details of the amino acid sequence; (iii) the moderns approaches to the prediction of the folding patterns of α- and β-regions in concrete proteins; (iv) experimental approaches to the mechanism of protein folding. The review reflects theoretical and experimental works of the author and his collaborators as well as those of other groups.     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

A heteropolymer model study for the mechanism of protein folding

A heteropolymer model of randomly self-interacting chains in two dimensions is studied with numerical simulations in order to elucidate the folding mechanism of protein. We find that the model occasionally shows folding propensity depending on the sequence of random numbers given to the chain. We study the thermodynamic and kinematic roles in the folding mechanism by grouping the local energy minima found in the simulations into clusters according to the similarity of their conformations. It is suggested that the local minima to which some heteropolymers show a folding tendency are always the lowest energy states of the energy spectrum within a cluster, though which cluster is selected depends on the sequence. For the eight random sequences we study, we find that the energy gap between the ground state and excited states is little correlated with folding or nonfolding. We rather find that folding propensities are correlated with the global structure of the average energy surface, implying a dominant kinetic role in the folding mechanism, although thermal factors cannot be ignored as the mechanism of choosing the ground state within a cluster of states connected by small deformations. We suggest that a hierarchical cluster structure plays an important role in selecting a unique folded state out of the huge number of local minima of heteropolymers. © 1997 John Wiley & Sons, Inc.     

What determines protein folding type? An investigation of intrinsic structural properties and its implications for understanding folding mechanisms

Protein folding experiments demonstrate that the folding behaviors of many proteins can be roughly classified into two types: two-state kinetics and multi-state kinetics. Although the two types of protein folding kinetics have been observed for a long time, what determines the folding type of a protein is still largely unclear. The present work performed a comparative study based on a dataset of 43 two-state and 42 multi-state folders at different levels of proteins' intrinsic properties from the simplest sequence length to native structure topology. The results show that protein's amino acids composition and the long-range interaction-based topological complexity rather than secondary structure contents are the major determinants of protein folding type. Furthermore, a sequence-based folding type prediction achieved an accuracy of more than 80%. These findings implicate that there is no clear boundary between secondary and tertiary structure formation during the protein folding process and support the existence of a continuum of folding mechanism between the two ends of hierarchic and nucleation folding scenarios.     

Characterization of a quaternary-structured folding intermediate of an antibody Fab-fragment.

Antibody folding is a complex process comprising folding and association reactions. Although it is usually difficult to characterize kinetic folding intermediates, in the case of the antibody Fab fragment, domain-domain interactions lead to a rate-limiting step of folding, thus accumulating folding intermediates at a late step of folding. Here, we analyzed a late folding intermediate of the Fab fragment of the monoclonal antibody MAK 33 from mouse (kappa/IgG1). As a strategy for accumulation of this intermediate we used partial denaturation of the native Fab by guanidinium chloride. This denaturation intermediate, which can be populated to about 90%, is indistinguishable from a late-folding intermediate with respect to denaturation and renaturation kinetics. The spectroscopic analysis reveals a native-like secondary structure of this intermediate with aromatic side chains only slightly more solvent exposed than in the native state. The respective partner domains are weekly associated. From these data we conclude that the intramolecular association of the two chains during folding, with all domains in a native-like structure, follows a two-step mechanism. In this mechanism, presumably hydrophobic interactions are followed by rearrangements leading to the exact complementarity of the contact sites of the respective domains.     

Secondary and tertiary structure formation of the beta-barrel membrane protein OmpA is synchronized and depends on membrane thickness

The mechanism of membrane insertion and folding of a beta-barrel membrane protein has been studied using the outer membrane protein A (OmpA) as an example. OmpA forms an eight-stranded beta-barrel that functions as a structural protein and perhaps as an ion channel in the outer membrane of Escherichia coli. OmpA folds spontaneously from a urea-denatured state into lipid bilayers of small unilamellar vesicles. We have used fluorescence spectroscopy, circular dichroism spectroscopy, and gel electrophoresis to investigate basic mechanistic principles of structure formation in OmpA. Folding kinetics followed a second-order rate law and is strongly depended on the hydrophobic thickness of the lipid bilayer. When OmpA was refolded into model membranes of dilaurylphosphatidylcholine, fluorescence kinetics were characterized by a rate constant that was about fivefold higher than the rate constants of formation of secondary and tertiary structure, which were determined by circular dichroism spectroscopy and gel electrophoresis, respectively. The formation of beta-sheet secondary structure and closure of the beta-barrel of OmpA were correlated with the same rate constant and coupled to the insertion of the protein into the lipid bilayer. OmpA, and presumably other beta-barrel membrane proteins therefore do not follow a mechanism according to the two-stage model that has been proposed for the folding of alpha-helical bundle membrane proteins. These different folding mechanisms are likely a consequence of the very different intramolecular hydrogen bonding and hydrophobicity patterns in these two classes of membrane proteins.     

Simulating the Folding Pathway of RNA Secondary Structure Using the Modified Ant Colony Algorithm

<正> A new method for simulating the folding pathway of RNA secondary structure using the modified ant colony algorithmis proposed.For a given RNA sequence,the set of all possible stems is obtained and the energy of each stem iscalculated and stored at the initial stage.Furthermore,a more realistic formula is used to compute the energy ofmulti-branch loop in the following iteration.Then a folding pathway is simulated,including such processes as constructionof the heuristic information,the rule of initializing the pheromone,the mechanism of choosing the initial andnext stem and the strategy of updating the pheromone between two different stems.Finally by testing RNA sequences withknown secondary structures from the public databases,we analyze the experimental data to select appropriate values forparameters.The measure indexes show that our procedure is more consistent with phylogenetically proven structures thansoftware RNAstructure sometimes and more effective than the standard Genetic Algorithm.     

同义密码子使用模式对蛋白产物表达及构象形成的影响*

鉴于遗传密码子的简并性能够将基因遗传信息的容量提升,同义密码子使用偏嗜性得以在生物体的基因组中广泛存在。虽然同义密码子之间碱基的变化并不能导致氨基酸种类的改变,在研究mRNA半衰期、编码多肽翻译效率及肽链空间构象正确折叠的准确性和翻译等这一系列过程中发现,同义密码子使用的偏嗜性在某种程度上通过精微调控翻译机制体现其遗传学功能。同义密码子指导tRNA在翻译过程中识别核糖体的速率变化是由氨基酸的特定顺序决定,并且在新生多肽链合成时,蛋白质共翻译转运机制同时调节其空间构象的正确折叠从而保证蛋白的正常生物学功能。某些同义密码子使用偏嗜性与特定蛋白结构的形成具有显著相关性,密码子使用偏嗜性一旦改变将可能导致新生多肽空间构象出现错误折叠。结合近些年来国内外在此领域的研究成果,阐述同义密码子使用偏嗜性如何发挥精微调控翻译的生物学功能与作用。     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Malleability of the folding mechanism of the outer membrane protein PagP: parallel pathways and the effect of membrane elasticity

Understanding the interactions between membrane proteins and the lipid bilayer is key to increasing our ability to predict and tailor the folding mechanism, structure and stability of membrane proteins. Here, we have investigated the effects of changing the membrane composition and the relative concentrations of protein and lipid on the folding mechanism of the bacterial outer membrane protein PagP. The folding pathway, monitored by tryptophan fluorescence, was found to be characterized by a burst phase, representing PagP adsorption to the liposome surface, followed by a time course that reflects the folding and insertion of the protein into the membrane. In 1,2-dilauroyl-sn-glycero-3-phosphocholine (diC(12:0)PC) liposomes, the post-adsorption time course fits well to a single exponential at high lipid-to-protein ratios (LPRs), but at low LPRs, a second exponential phase with a slower folding rate constant is observed. Interrupted refolding assays demonstrated that the two exponential phases reflect the presence of parallel folding pathways. Partitioning between these pathways was found to be modulated by the elastic properties of the membrane. Folding into mixed 1,2-dilauroyl-sn-glycero-3-phosphoethanolamine:diC(12:0)PC liposomes resulted in a decrease in PagP adsorption to the liposomes and a switch to the slower folding pathway. By contrast, inclusion of 1,2-dilauroyl-sn-glycero-3-phosphoserine into diC(12:0)PC liposomes resulted in a decrease in the folding rate of the fast pathway. The results highlight the effect of lipid composition in tailoring the folding mechanism of a membrane protein, revealing that membrane proteins have access to multiple, competing folding routes to a unique native structure.     

Pathway of oxidative folding of bovine alpha-interferon: predominance of native disulfide-bonded folding intermediates

Bovine alpha-interferon (BoINF-alpha) is a single polypeptide protein containing 166 amino acids, two disulfide bonds (Cys1-Cys99 and Cys29-Cys138), and five stretches of alpha-helical structure. The pathway of oxidative folding of BoINF-alpha has been investigated here. Of the eight possible one- and two-disulfide isomers, only two nativelike one-disulfide isomers, BoINF-alpha (Cys1-Cys99) and BoINF-alpha (Cys29-Cys138), predominate as intermediates along the folding pathway. More strikingly, alpha-helical structures formed almost quantitatively before any detectable formation of a disulfide bond. This is demonstrated by the observation that fully reduced BoINF-alpha (starting material of oxidative folding) and reduced carboxymethylated BoINF-alpha both exhibit alpha-helical structure content indistinguishable form that of native BoINF-alpha. The folding mechanism of BoINF-alpha appears to be compatible with the framework model, in which secondary structures fold first, followed by docking (compaction) of preformed secondary structural elements yielding the native structure.     

Topology-based modeling of intrinsically disordered proteins: balancing intrinsic folding and intermolecular interactions

Coupled binding and folding is frequently involved in specific recognition of so-called intrinsically disordered proteins (IDPs), a newly recognized class of proteins that rely on a lack of stable tertiary fold for function. Here, we exploit topology-based Gō-like modeling as an effective tool for the mechanism of IDP recognition within the theoretical framework of minimally frustrated energy landscape. Importantly, substantial differences exist between IDPs and globular proteins in both amino acid sequence and binding interface characteristics. We demonstrate that established Gō-like models designed for folded proteins tend to over-estimate the level of residual structures in unbound IDPs, whereas under-estimating the strength of intermolecular interactions. Such systematic biases have important consequences in the predicted mechanism of interaction. A strategy is proposed to recalibrate topology-derived models to balance intrinsic folding propensities and intermolecular interactions, based on experimental knowledge of the overall residual structure level and binding affinity. Applied to pKID/KIX, the calibrated Gō-like model predicts a dominant multistep sequential pathway for binding-induced folding of pKID that is initiated by KIX binding via the C-terminus in disordered conformations, followed by binding and folding of the rest of C-terminal helix and finally the N-terminal helix. This novel mechanism is consistent with key observations derived from a recent NMR titration and relaxation dispersion study and provides a molecular-level interpretation of kinetic rates derived from dispersion curve analysis. These case studies provide important insight into the applicability and potential pitfalls of topology-based modeling for studying IDP folding and interaction in general.     

The thermodynamic stability of insulin disulfides is not affected by the C-domain of insulin-like growth factor 1

Both Insulin and insulin-like growth factor 1 are members of insulin superfamily. They share homologous primary and tertiary structure as well as weakly overlapping biological activity. However, their folding behavior is different: insulin and its recombinant precursor (PIP) fold into one unique tertiary structure, while IGF-1 folds into two disulfides isomers with similar thermody-namic stability. To elucidate the molecular mechanism of their different folding behavior, we prepared a single-chain hybrid of insulin and IGF-1, [B10Glu]lns/IGF-1(C), and studied its folding behavior compared with that of PIP and IGF-1. We also separated a major non-native disulfides iso-mer of the hybrid and studied its refolding. The data showed that the C-domain of IGF-1 did not affect the folding thermodynamics of insulin, that is, the primary structure of the hybrid encoded only one thermodynamically stable disulfides linkage. However, the folding kinetics of insulin was affected by the C-domain of IGF-1.     

The thermodynamic stability of insulin disulfides is not affected by the C-domain of insulin-like growth factor 1

Both Insulin and insulin-like growth factor 1 are members of insulin superfamily. They share homologous primary and tertiary structure as well as weakly overlapping biological activity. However, their folding behavior is different: insulin and its recombinant precursor (PIP) fold into one unique tertiary structure, while IGF-1 folds into two disulfides isomers with similar thermodynamic stability. To elucidate the molecular mechanism of their different folding behavior, we prepared a singlechain hybrid of insulin and IGF-1, [B10Glu]Ins/IGF-1(C), and studied its folding behavior compared with that of PIP and IGF-1. We also separated a major non-native disulfides isomer of the hybrid and studied its refolding. The data showed that the C-domain of IGF-1 did not affect the folding thermodynamics of insulin, that is, the primary structure of the hybrid encoded only one thermodynamically stable disulfides linkage. However, the folding kinetics of insulin was affected by the C-domain of IGF-1.     

How the folding rates of two‐ and multistate proteins depend on the amino acid properties

Proteins fold by either two‐state or multistate kinetic mechanism. We observe that amino acids play different roles in different mechanism. Many residues that are easy to form regular secondary structures (α helices, β sheets and turns) can promote the two‐state folding reactions of small proteins. Most of hydrophilic residues can speed up the multistate folding reactions of large proteins. Folding rates of large proteins are equally responsive to the flexibility of partial amino acids. Other properties of amino acids (including volume, polarity, accessible surface, exposure degree, isoelectric point, and phase transfer energy) have contributed little to folding kinetics of the proteins. Cysteine is a special residue, it triggers two‐state folding reaction and but inhibits multistate folding reaction. These findings not only provide a new insight into protein structure prediction, but also could be used to direct the point mutations that can change folding rate. Proteins 2014; 82:2375–2382. © 2014 Wiley Periodicals, Inc.     

The folding pathway of ubiquitin from all-atom molecular dynamics simulations

The folding (unfolding) pathway of ubiquitin is probed using all-atom molecular dynamics simulations. We dissect the folding pathway using two techniques: first, we probe the folding pathway of ubiquitin by calculating the evolution of structural properties over time and second, we identify the rate determining transition state for folding. The structural properties that we look at are hydrophobic solvent accessible surface area (SASA) and Calpha-root-mean-square deviation (rmsd). These properties on their own tell us relatively little about the folding pathway of ubiquitin; however, when plotted against each other, they become powerful tools for dissecting ubiquitin's folding mechanism. Plots of Calpha-rmsd against SASA serve as a phase space trajectories for the folding of ubiquitin. In this study, these plots show that ubiquitin folds to the native state via the population of an intermediate state. This is shown by an initial hydrophobic collapse phase followed by a second phase of secondary structure arrangement. Analysis of the structure of the intermediate state shows that it is a collapsed species with very little secondary structure. In reconciling these observations with recent experimental data, the transition that we observe in our simulations from the unfolded state (U) to the intermediate state (I) most likely occurs in the dead-time of the stopped flow instrument. The folding pathway of ubiquitin is probed further by identification of the rate-determining transition state for folding. The method used for this is essential dynamics, which utilizes a principal component analysis (PCA) on the atomic fluctuations throughout the simulation. The five transition state structures identified in silico are in good agreement with the experimentally determined transition state. The calculation of phi-values from the structures generated in the simulations is also carried out and it shows a good correlation with the experimentally measured values. An initial analysis of the denatured state shows that it is compact with fluctuating regions of nonnative secondary structure. It is found that the compactness in the denatured state is due to the burial of some hydrophobic residues. We conclude by looking at a correlation between folding kinetics and residual structure in the denatured state. A hierarchical folding mechanism is then proposed for ubiquitin.     

Contact patterns between helices and strands of sheet define protein folding patterns

Comparing and classifying protein folding patterns allows organizing the known structures and enumerating possible protein structural patterns including those not yet observed. We capture the essence of protein folding patterns in a concise tableau representation based on the order and contact patterns of secondary structures: helices and strands of sheet. The tableaux are intelligible to both humans and computers. They provide a database, derived from the Protein Data Bank, mineable in studies of protein architecture. Using this database, we have: (i) determined statistical properties of secondary structure contacts in an unbiased set of protein domains from ASTRAL, (ii) observed that in 98% of cases, the tableau is a faithful representation of the folding pattern as classified in SCOP, (iii) demonstrated that to a large extent the local structure of proteins indicates their complete folding topology, and (iv) studied the use of the representation for fold identification.     

Peptide folding driven by Van der Waals interactions

Contrary to the widespread view that hydrogen bonding and its entropy effect play a dominant role in protein folding, folding into helical and hairpin-like structures is observed in molecular dynamics (MD) simulations without hydrogen bonding in the peptide-solvent system. In the widely used point charge model, hydrogen bonding is calculated as part of the interaction between atomic partial charges. It is removed from these simulations by setting atomic charges of the peptide and water to zero. Because of the structural difference between the peptide and water, van der Waals (VDW) interactions favor peptide intramolecular interactions and are a major contributing factor to the structural compactness. These compact structures are amino acid sequence dependent and closely resemble standard secondary structures, as a consequence of VDW interactions and covalent bonding constraints. Hydrogen bonding is a short range interaction and it locks the approximate structure into the specific secondary structure when it is included in the simulation. In contrast to standard molecular simulations where the total energy is dominated by charge-charge interactions, these simulation results will give us a new view of the folding mechanism.