An amino acid code for protein folding

Structural Biology in the AlphaFold Era: How Far Is Artificial Intelligence from Deciphering the Protein Folding Code?

Proteins are biomolecules characterized by uncommon chemical and physicochemical complexities coupled with extreme responsiveness to even minor chemical modifications or environmental variations. Since the shape that proteins assume is fundamental for their function, understanding the chemical and structural bases that drive their three-dimensional structures represents the central problem for an atomic-level interpretation of biology. Not surprisingly, this question has progressively become the Holy Grail of structural biology (the folding problem). From this perspective, we initially describe and discuss the different formulations of the folding problem. In the present manuscript, the folding problem is framed from a historical perspective, effectively highlighting the progress made in the last lustrum. We chronologically summarize the major contributions that traditional methodologies provide in approaching this multifaceted problem. We then describe the recent advent and evolution of predictive approaches based on machine learning techniques that are revolutionizing the field by pointing out the potentialities and limitations of this approach. In the final part of the perspective, we illustrate the contribution that computational approaches will make in current structural biology to overcome the limitations of the reductionist approach of studying individual molecules to afford the atomic-level characterization of entire cellular compartments.

Global Patterns of Subgenome Evolution in Organelle-Targeted Genes of Six Allotetraploid Angiosperms

Whole-genome duplications (WGDs) are a prominent process of diversification in eukaryotes. The genetic and evolutionary forces that WGD imposes on cytoplasmic genomes are not well understood, despite the central role that cytonuclear interactions play in eukaryotic function and fitness. Cellular respiration and photosynthesis depend on successful interaction between the 3,000+ nuclear-encoded proteins destined for the mitochondria or plastids and the gene products of cytoplasmic genomes in multi-subunit complexes such as OXPHOS, organellar ribosomes, Photosystems I and II, and Rubisco. Allopolyploids are thus faced with the critical task of coordinating interactions between the nuclear and cytoplasmic genes that were inherited from different species. Because the cytoplasmic genomes share a more recent history of common descent with the maternal nuclear subgenome than the paternal subgenome, evolutionary "mismatches" between the paternal subgenome and the cytoplasmic genomes in allopolyploids might lead to the accelerated rates of evolution in the paternal homoeologs of allopolyploids, either through relaxed purifying selection or strong directional selection to rectify these mismatches. We report evidence from six independently formed allotetraploids that the subgenomes exhibit unequal rates of protein-sequence evolution, but we found no evidence that cytonuclear incompatibilities result in altered evolutionary trajectories of the paternal homoeologs of organelle-targeted genes. The analyses of gene content revealed mixed evidence for whether the organelle-targeted genes are lost more rapidly than the non-organelle-targeted genes. Together, these global analyses provide insights into the complex evolutionary dynamics of allopolyploids, showing that the allopolyploid subgenomes have separate evolutionary trajectories despite sharing the same nucleus, generation time, and ecological context.

Cytochrome c: Using Biological Insight toward Engineering an Optimized Anticancer Biodrug

The heme protein cytochrome c (Cyt c) plays pivotal roles in cellular life and death processes. In the respiratory chain of mitochondria, it serves as an electron transfer protein, contributing to the proliferation of healthy cells. In the cell cytoplasm, it activates intrinsic apoptosis to terminate damaged cells. Insight into these mechanisms and the associated physicochemical properties and biomolecular interactions of Cyt c informs on the anticancer therapeutic potential of the protein, especially in its ability to subvert the current limitations of small molecule-based chemotherapy. In this review, we explore the development of Cyt c as an anticancer drug by identifying cancer types that would be receptive to the cytotoxicity of the protein and factors that can be finetuned to enhance its apoptotic potency. To this end, some information is obtained by characterizing known drugs that operate, in part, by triggering Cyt c induced apoptosis. The application of different smart drug delivery systems is surveyed to highlight important features for maintaining Cyt c stability and activity and improving its specificity for cancer cells and high drug payload release while recognizing the continuing limitations. This work serves to elucidate on the optimization of the strategies to translate Cyt c to the clinical market.

Protein folding - seeing is deceiving

This Perspective is intended to raise questions about the conventional interpretation of protein folding. According to the conventional interpretation, developed over many decades, a protein population can visit a vast number of conformations under unfolding conditions, but a single dominant native population emerges under folding conditions. Accordingly, folding comes with a substantial loss of conformational entropy. How is this price paid? The conventional answer is that favorable interactions between and among the side chains can compensate for entropy loss, and moreover, these interactions are responsible for the structural particulars of the native conformation. Challenging this interpretation, the Perspective introduces a proposal that high energy (i.e., unfavorable) excluding interactions winnow the accessible population substantially under physical-chemical conditions that favor folding. Both steric clash and unsatisfied hydrogen bond donors and acceptors are classified as excluding interactions, so called because conformers with such disfavored interactions will be largely excluded from the thermodynamic population. Both excluding interactions and solvent factors that induce compactness are somewhat nonspecific, yet together they promote substantial chain organization. Moreover, proteins are built on a backbone scaffold consisting of α-helices and strands of β-sheet, where the number of hydrogen bond donors and acceptors is exactly balanced. These repetitive secondary structural elements are the only two conformers that can be both completely hydrogen-bond satisfied and extended indefinitely without encountering a steric clash. Consequently, the number of fundamental folds is limited to no more than ~10,000 for a protein domain. Once excluding interactions are taken into account, the issue of "frustration" is largely eliminated and the Levinthal paradox is resolved. Putting the "bottom line" at the top: it is likely that hydrogen-bond satisfaction represents a largely under-appreciated parameter in protein folding models.

Looking Beyond the Core: The Role of Flanking Regions in the Aggregation of Amyloidogenic Peptides and Proteins

Amyloid proteins are involved in many neurodegenerative disorders such as Alzheimer's disease [Tau, Amyloid β (Aβ)], Parkinson's disease [alpha-synuclein (αSyn)], and amyotrophic lateral sclerosis (TDP-43). Driven by the early observation of the presence of ordered structure within amyloid fibrils and the potential to develop inhibitors of their formation, a major goal of the amyloid field has been to elucidate the structure of the amyloid fold at atomic resolution. This has now been achieved for a wide variety of sequences using solid-state NMR, microcrystallography, X-ray fiber diffraction and cryo-electron microscopy. These studies, together with in silico methods able to predict aggregation-prone regions (APRs) in protein sequences, have provided a wealth of information about the ordered fibril cores that comprise the amyloid fold. Structural and kinetic analyses have also shown that amyloidogenic proteins often contain less well-ordered sequences outside of the amyloid core (termed here as flanking regions) that modulate function, toxicity and/or aggregation rates. These flanking regions, which often form a dynamically disordered "fuzzy coat" around the fibril core, have been shown to play key parts in the physiological roles of functional amyloids, including the binding of RNA and in phase separation. They are also the mediators of chaperone binding and membrane binding/disruption in toxic amyloid assemblies. Here, we review the role of flanking regions in different proteins spanning both functional amyloid and amyloid in disease, in the context of their role in aggregation, toxicity and cellular (dys)function. Understanding the properties of these regions could provide new opportunities to target disease-related aggregation without disturbing critical biological functions.

Exploring the candidates for a new protein folding - cross-α amyloid - in available protein databases

Amyloid fibrils are formed from the assembly of soluble proteins and are responsible for many diseases. They are known to have a cross-β structure, where the fibril runs perpendicular to the β-sheets. A new type of tertiary structure formed by the aggregation of peptides in their α-helical form, in naturally occurring as well as synthetic peptides, termed cross-α amyloid has been reported recently. We have studied the interactions responsible for the formation of these cross-α amyloids and proposed a model to determine the peptides that could form these structures. Eight such peptides obtained using the model have been shown to form a cross-α structure using molecular dynamics simulations. The formation of a cross-α structure from eight copies of a randomly chosen peptide and its stability over a microsecond simulation have been demonstrated. A software named Cross-Alpha-Det has been developed that can determine whether a protein can form a cross-α structure from its secondary structure.

Paramagnetic relaxation enhancement-assisted structural characterization of a partially disordered conformation of ubiquitin

Nuclear magnetic resonance (NMR) is a powerful tool to study three-dimensional structures as well as protein conformational fluctuations in solution, but it is compromised by increases in peak widths and missing signals. We previously reported that ubiquitin has two folded conformations, N₁ and N₂ and plus another folded conformation, I, in which some amide group signals of residues 33-41 almost disappeared above 3 kbar at pH 4.5 and 273 K. Thus, well-converged structural models could not be obtained for this region owing to the absence of distance restraints. Here, we reexamine the problem using the ubiquitin Q41N variant as a model for this locally disordered conformation, I. We demonstrate that the variant shows pressure-induced loss of backbone amide group signals at residues 28, 33, 36, and 39-41 like the wild-type, with a similar but smaller effect on CαH and CβH signals. In order to characterize this I structure, we measured paramagnetic relaxation enhancement (PRE) under high pressure to obtain distance restraints, and calculated the structure assisted by Bayesian inference. We conclude that the more disordered I conformation observed at pH 4.0, 278 K, and 2.5 kbar largely retained the N₂ conformation, although the amide groups at residues 33-41 have more heterogeneous conformations and more contact with water, which differ from the N₁ and N₂ states. The PRE-assisted strategy has the potential to improve structural characterization of proteins that lack NMR signals, especially for relatively more open and hydrated protein conformations.

Autonomous Sequences in Myoglobin Emerging from X-ray Structure of Holomyoglobin

The proposed continuous folding structure units are fundamental to analyze protein structure. Here, we could elucidate for the first time two types of hydrophobic core networks in apomyoglobin using continuous folding structure units. In myoglobin, two autonomous sequences emerged clearly. We could thus characterize the autonomous sequences using well-defined hydrophobic core networks within respective semifolds. A hydrophobic core is defined as a pair of topology-local hydrophobic amino acids in different folding structures. Hydrophobic core formation is indispensable to stabilize the different folding structures via an efficient hydrophobic interaction. Autonomous sequences in myoglobin encode tertiary structure information for semifolds. These sequences fold autonomously into small sets of continuous folding structure units to grow separate semifolds on each separate framework. The autonomous sequence can be defined as the local sequence assigned to the small set of continuous folding structure units. They create the discrete hydrophobic region in a semifold by assembly of their hydrophobic regions. Semifolds were characterized by discrete hydrophobic regions stabilized by respective type I hydrophobic core networks, which were present within each semifold. The discrete hydrophobic region of a semifold propagated itself with that of a different semifold by hydrophobic interactions in type II hydrophobic core network, which was present between different semifolds, as observed by the X-ray structures of semifolds. The most significant feature of semifolds in apomyoglobin was that they could be verified by the X-ray structure of holomyoglobin regardless of the instability of folds characteristic to autonomous sequence fragments. This work presents the first description of autonomous sequences.

Radical amino acid mutations persist longer in the absence of sex

Harmful mutations are ubiquitous and inevitable, and the rate at which these mutations are removed from populations is a critical determinant of evolutionary fate. Closely related sexual and asexual taxa provide a particularly powerful setting to study deleterious mutation elimination because sexual reproduction should facilitate mutational clearance by reducing selective interference between sites and by allowing the production of offspring with different mutational complements than their parents. Here, we compared the rate of removal of conservative (i.e., similar biochemical properties) and radical (i.e., distinct biochemical properties) nonsynonymous mutations from mitochondrial genomes of sexual versus asexual Potamopyrgus antipodarum, a New Zealand freshwater snail characterized by coexisting and ecologically similar sexual and asexual lineages. Our analyses revealed that radical nonsynonymous mutations are cleared at higher rates than conservative changes and that sexual lineages eliminate radical changes more rapidly than asexual counterparts. These results are consistent with reduced efficacy of purifying selection in asexual lineages allowing harmful mutations to remain polymorphic longer than in sexual lineages. Together, these data illuminate some of the population-level processes contributing to mitochondrial mutation accumulation and suggest that mutation accumulation could influence the outcome of competition between sexual and asexual lineages.

Protein Folding Structures: Formation of Folding Structures Based on Probability Theory

To the best of our knowledge, this is the first study that shows that the X-ray structures of proteins can be dissected into their continuous folding structure units. Each folding structure unit was designed such that both the terminal di- or tri-peptide sequences shared common sequences with the two adjacent folding structure units. To encode the folding structure information of proteins into their amino acid sequences, we proposed 44 kinds of folding elements, which covered all of the amino acids in the protein chains, and defined all folding structure units. The folding element was defined to mean a minimum structural piece, which covered the frame of the main chain of each amino acid in a protein chain. A folding structure unit of a local sequence could be fully characterized by the sequential combination of individual folding elements assigned to each amino acid. The folding structure information showed amino acid preferences in various positions in folding structure units. Folding structure formation proceeded on the basis of probability theory. Strikingly, relative formation ability analysis clearly indicated that we can decode the types and the chain length of folding structure units from the amino acid sequence of a protein.

The role of the backbone torsion in protein folding

BACKGROUND

The set of forces and sequence of events that govern the transition from an unfolded polypeptide chain to a functional protein with correct spatial structure remain incompletely known, despite the importance of the problem and decades of theory development, computer simulations, and laboratory experiments. Information about the correctly folded state of most proteins is likely to be present in their sequences, and yet many proteins fail to attain native structure after overexpression in a non-native environment or upon experimental denaturation and refolding.

PRESENTATION OF THE HYPOTHESIS

We hypothesize that correct protein folding in vivo is an active, energy-dependent process that most likely applies torque force co-translationally to all proteins and possibly also post-translationally to many proteins in every cell. When a site on an unfolded polypeptide is rotationally constrained, torsion applied at another site would induce twisting of the main chain, which would initiate the formation of a local secondary structure, such as an alpha-helical turn or a beta-turn/beta-hairpin. The nucleation of structural elements is a rate-limiting, energetically unfavorable step in the process of protein folding, and energy-dependent chain torsion is likely to help overcome this barrier in vivo. Several molecular machines in a cell, primarily ribosomes, but also possibly signal recognition particles and chaperone systems, may play a role in applying torque to an unfolded protein chain, using the energy of GTP or ATP hydrolysis. Lack of such force in the in vitro systems may be the main reason of the failure of many longer proteins to attain the correct functional conformation.

TESTING THE HYPOTHESIS

The hypothesis can be tested using single-molecule approaches, by measuring directly the forces applied to polypeptide chains under controlled conditions in vitro, and in bulk, by assessing folding rates and extent of misfolding in proteins that are engineered to experience transient spatial constraint during their synthesis.

IMPLICATIONS OF THE HYPOTHESIS

Learning about the role of main chain torsion in protein folding will improve our understanding of folding mechanisms and may lead to bioengineering solutions that would enhance the yields of correctly folded proteins in heterologous expression systems.

REVIEWERS

This article was reviewed by Frank Eisenhaber, Igor Berezovsky and Michael Gromiha.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Multiple Docking for Protein Structure Prediction

Protein structure prediction and protein docking prediction are two related problems in molecular biology. We suggest the use of multiple docking in the process of protein structure prediction. Once reliable structural models are predicted to disjoint fragments of the protein target sequence, a combinatorial assembly may be used to predict their native arrangement. Here, we present CombDock, a combinatorial docking algorithm for the structural units assembly problem. We have tested the algorithm on various examples using both domains and domain substructures as input. Inaccurate models of the structural units were also used, to test the robustness of the algorithm. The algorithm was able to predict a near-native arrangement of the input structural units in almost all of the cases, showing that the combinatorial approach succeeds in overcoming the inexact shape complementarity caused by the inaccuracy of the models.

UniCon3D: de novo protein structure prediction using united-residue conformational search via stepwise, probabilistic sampling

MOTIVATION

Recent experimental studies have suggested that proteins fold via stepwise assembly of structural units named 'foldons' through the process of sequential stabilization. Alongside, latest developments on computational side based on probabilistic modeling have shown promising direction to perform de novo protein conformational sampling from continuous space. However, existing computational approaches for de novo protein structure prediction often randomly sample protein conformational space as opposed to experimentally suggested stepwise sampling.

RESULTS

Here, we develop a novel generative, probabilistic model that simultaneously captures local structural preferences of backbone and side chain conformational space of polypeptide chains in a united-residue representation and performs experimentally motivated conditional conformational sampling via stepwise synthesis and assembly of foldon units that minimizes a composite physics and knowledge-based energy function for de novo protein structure prediction. The proposed method, UniCon3D, has been found to (i) sample lower energy conformations with higher accuracy than traditional random sampling in a small benchmark of 6 proteins; (ii) perform comparably with the top five automated methods on 30 difficult target domains from the 11th Critical Assessment of Protein Structure Prediction (CASP) experiment and on 15 difficult target domains from the 10th CASP experiment; and (iii) outperform two state-of-the-art approaches and a baseline counterpart of UniCon3D that performs traditional random sampling for protein modeling aided by predicted residue-residue contacts on 45 targets from the 10th edition of CASP.

AVAILABILITY AND IMPLEMENTATION

Source code, executable versions, manuals and example data of UniCon3D for Linux and OSX are freely available to non-commercial users at http://sysbio.rnet.missouri.edu/UniCon3D/ CONTACT: chengji@missouri.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Network representation of conformational transitions between hidden intermediates of Rd-apocytochrome b562

The folding kinetics of Rd-apocytochrome b562 is two-state, but native-state hydrogen exchange experiments show that there are discrete partially unfolded (PUF) structures in equilibrium with the native state. These PUF structures are called hidden intermediates because they are not detected in kinetic experiments and they exist after the rate-limiting step. Structures of the mimics of hidden intermediates of Rd-apocytochrome b562 are resolved by NMR. Based upon their relative stability and structural features, the folding mechanism was proposed to follow a specific pathway (unfolded → rate-limiting transition state → PUF1 → PUF2 → native). Investigating the roles of equilibrium PUF structures in folding kinetics and their interrelationship not only deepens our understanding of the details of folding mechanism but also provides guides in protein design and prevention of misfolding. We performed molecular dynamics simulations starting from a hidden intermediate and the native state of Rd-apocytochrome b562 in explicit solvent, for a total of 37.18 μs mainly with Anton. We validated our simulations by detailed comparison with experimental data and other computations. We have verified that we sampled the post rate-limiting transition state region only. Markov state model was used to analyze the simulation results. We replace the specific pathway model with a network model. Transition-path theory was employed to calculate the net effective flux from the most unfolded state towards the most folded state in the network. The proposed sequential folding pathway via PUF1 then more stable, more native-like PUF2 is one of the routes in our network, but it is not dominant. The dominant path visits PUF2 without going through PUF1. There is also a route from PUF1 directly to the most folded state in the network without visiting PUF2. Our results indicate that the PUF states are not necessarily sequential in the folding. The major routes predicted in our network are testable by future experiments such as single molecule experiment.

General mechanism of two-state protein folding kinetics

We describe here a general model of the kinetic mechanism of protein folding. In the Foldon Funnel Model, proteins fold in units of secondary structures, which form sequentially along the folding pathway, stabilized by tertiary interactions. The model predicts that the free energy landscape has a volcano shape, rather than a simple funnel, that folding is two-state (single-exponential) when secondary structures are intrinsically unstable, and that each structure along the folding path is a transition state for the previous structure. It shows how sequential pathways are consistent with multiple stochastic routes on funnel landscapes, and it gives good agreement with the 9 order of magnitude dependence of folding rates on protein size for a set of 93 proteins, at the same time it is consistent with the near independence of folding equilibrium constant on size. This model gives estimates of folding rates of proteomes, leading to a median folding time in Escherichia coli of about 5 s.

Measurement of histidine pKa values and tautomer populations in invisible protein states

The histidine imidazole side chain plays a critical role in protein function and stability. Its importance for catalysis is underscored by the fact that histidines are localized to active sites in ∼ 50% of all enzymes. NMR spectroscopy has become an important tool for studies of histidine side chains through the measurement of site-specific pK(a)s and tautomer populations. To date, such studies have been confined to observable protein ground states; however, a complete understanding of the role of histidine electrostatics in protein function and stability requires that similar investigations be extended to rare, transiently formed conformers that populate the energy landscape, yet are often "invisible" in standard NMR spectra. Here we present NMR experiments and a simple strategy for studies of such conformationally excited states based on measurement of histidine (13)Cγ, (13)Cδ2 chemical shifts and (1)Hε-(13)Cε one-bond scalar couplings. The methodology is first validated and then used to obtain pKa values and tautomer distributions for histidine residues of an invisible on-pathway folding intermediate of the colicin E7 immunity protein. Our results imply that the side chains of H40 and H47 are exposed in the intermediate state and undergo significant conformational rearrangements during folding to the native structure. Further, the pKa values explain the pH-dependent stability differences between native and intermediate states over the pH range 5.5-6.5 and they suggest that imidazole deprotonation is not a barrier to the folding of this protein.

Stepwise protein folding at near amino acid resolution by hydrogen exchange and mass spectrometry

The kinetic folding of ribonuclease H was studied by hydrogen exchange (HX) pulse labeling with analysis by an advanced fragment separation mass spectrometry technology. The results show that folding proceeds through distinct intermediates in a stepwise pathway that sequentially incorporates cooperative native-like structural elements to build the native protein. Each step is seen as a concerted transition of one or more segments from an HX-unprotected to an HX-protected state. Deconvolution of the data to near amino acid resolution shows that each step corresponds to the folding of a secondary structural element of the native protein, termed a "foldon." Each folded segment is retained through subsequent steps of foldon addition, revealing a stepwise buildup of the native structure via a single dominant pathway. Analysis of the pertinent literature suggests that this model is consistent with experimental results for many proteins and some current theoretical results. Two biophysical principles appear to dictate this behavior. The principle of cooperativity determines the central role of native-like foldon units. An interaction principle termed "sequential stabilization" based on native-like interfoldon interactions orders the pathway.

Mass spectrometry combined with oxidative labeling for exploring protein structure and folding

This review discusses various mass spectrometry (MS)-based approaches for exploring structural aspects of proteins in solution. Electrospray ionization (ESI)-MS, in particular, has found fascinating applications in this area. For example, when used in conjunction with solution-phase hydrogen/deuterium exchange (HDX), ESI-MS is a highly sensitive tool for probing conformational dynamics. The main focus of this article is a technique that is complementary to HDX, that is, the covalent labeling of proteins by hydroxyl radicals. The reactivity of individual amino acid side chains with *OH is strongly affected by their degree of solvent exposure. Thus, analysis of the oxidative labeling pattern by peptide mapping and tandem mass spectrometry provides detailed structural information. A convenient method for *OH production is the photolysis of H(2)O(2) by a pulsed UV laser, resulting in oxidative labeling on the microsecond time scale. Selected examples demonstrate the use of this technique for structural studies on membrane proteins, and the combination with rapid mixing devices for characterizing the properties of short-lived protein (un)folding intermediates.

Atomic interaction networks in the core of protein domains and their native folds

Vastly divergent sequences populate a majority of protein folds. In the quest to identify features that are conserved within protein domains belonging to the same fold, we set out to examine the entire protein universe on a fold-by-fold basis. We report that the atomic interaction network in the solvent-unexposed core of protein domains are fold-conserved, extraordinary sequence divergence notwithstanding. Further, we find that this feature, termed protein core atomic interaction network (or PCAIN) is significantly distinguishable across different folds, thus appearing to be “signature” of a domain's native fold. As part of this study, we computed the PCAINs for 8698 representative protein domains from families across the 1018 known protein folds to construct our seed database and an automated framework was developed for PCAIN-based characterization of the protein fold universe. A test set of randomly selected domains that are not in the seed database was classified with over 97% accuracy, independent of sequence divergence. As an application of this novel fold signature, a PCAIN-based scoring scheme was developed for comparative (homology-based) structure prediction, with 1–2 angstroms (mean 1.61A) C_α RMSD generally observed between computed structures and reference crystal structures. Our results are consistent across the full spectrum of test domains including those from recent CASP experiments and most notably in the ‘twilight’ and ‘midnight’ zones wherein <30% and <10% target-template sequence identity prevails (mean twilight RMSD of 1.69A). We further demonstrate the utility of the PCAIN protocol to derive biological insight into protein structure-function relationships, by modeling the structure of the YopM effector novel E3 ligase (NEL) domain from plague-causative bacterium Yersinia Pestis and discussing its implications for host adaptive and innate immune modulation by the pathogen. Considering the several high-throughput, sequence-identity-independent applications demonstrated in this work, we suggest that the PCAIN is a fundamental fold feature that could be a valuable addition to the arsenal of protein modeling and analysis tools.

NMR investigations on residue level unfolding thermodynamics in DLC8 dimer by temperature dependent native state hydrogen exchange

Understanding protein stability at residue level detail in the native state ensemble of a protein is crucial to understanding its biological function. At the same time, deriving thermodynamic parameters using conventional spectroscopic and calorimetric techniques remains a major challenge for some proteins due to protein aggregation and irreversibility of denaturation at higher temperature values. In this regard, we describe here the NMR investigations on the conformational stabilities and related thermodynamic parameters such as local unfolding enthalpies, heat capacities and transition midpoints in DLC8 dimer, by using temperature dependent native state hydrogen exchange; this protein aggregates at high (>65 degrees C) temperatures. The stability (free energy) of the native state was found to vary substantially with temperature at every residue. Significant differences were found in the thermodynamic parameters at individual residue sites indicating that the local environments in the protein structure would respond differently to external perturbations; this reflects on plasticity differences in different regions of the protein. Further, comparison of this data with similar data obtained from GdnHCl dependent native state hydrogen exchange indicated many similarities at residue level, suggesting that local unfolding transitions may be similar in both the cases. This has implications for the folding/unfolding mechanisms of the protein.

Lack of cytochrome c in mouse fibroblasts disrupts assembly/stability of respiratory complexes I and IV

Cytochrome c (cyt c) is a heme-containing protein that participates in electron transport in the respiratory chain and as a signaling molecule in the apoptotic cascade. Here we addressed the effect of removing mammalian cyt c on the integrity of the respiratory complexes in mammalian cells. Mitochondria from cyt c knockout mouse cells lacked fully assembled complexes I and IV and had reduced levels of complex III. A redox-deficient mutant of cyt c was unable to rescue the levels of complexes I and IV. We found that cyt c is associated with both complex IV and respiratory supercomplexes, providing a potential mechanism for the requirement for cyt c in the assembly/stability of complex IV.

Residue-wise conformational stability of DLC8 dimer from native-state hydrogen exchange

Dynein light chain (DLC8) is the smallest subunit of the dynein motor complex, which is known to act as a cargo adaptor in intracellular trafficking. The protein exists as a pure dimer at physiological pH and a completely folded monomer below pH 4. Here, we have determined the energy landscape of the dimeric protein using a combination of optical techniques and native-state hydrogen exchange of amide groups, the former giving the global features and the latter yielding the residue level details. The data indicated the presence of intermediates along the equilibrium unfolding transition. The hydrogen exchange data suggested that the molecule has differential stability in its various segments. We deduce from the free energy data that the antiparallel beta-sheets (beta4 and beta5) that form the hydrophobic core of the protein and the alpha2 helix, all of which are highly protected with regard to hydrogen exchange, contribute significantly to the initial step of the protein folding mechanism. Denaturant-dependent hydrogen exchange indicated further that some amides exchange via local fluctuations, whereas there are others which exchange via global unfolding events. Implications of these to cargo adaptability of the dimer are discussed.

Structure and energetics of the hydrogen-bonded backbone in protein folding

We seek to understand the link between protein thermodynamics and protein structure in molecular detail. A classical approach to this problem involves assessing changes in protein stability resulting from added cosolvents. Under any given conditions, protein molecules in aqueous buffer are in equilibrium between unfolded and folded states, U(nfolded) <==> N(ative). Addition of organic osmolytes, small uncharged compounds found throughout nature, shift this equilibrium. Urea, a denaturing osmolyte, shifts the equilibrium toward U; trimethylamine N-oxide (TMAO), a protecting osmolyte, shifts the equilibrium toward N. Using the Tanford Transfer Model, the thermodynamic response to many such osmolytes has been dissected into groupwise free energy contributions. It is found that the energetics involving backbone hydrogen bonding controls these shifts in protein stability almost entirely, with osmolyte cosolvents simply dialing between solvent-backbone versus backbone-backbone hydrogen bonds, as a function of solvent quality. This reciprocal relationship establishes the essential link between protein thermodynamics and the protein's hydrogen-bonded backbone structure.

Protein folding: independent unrelated pathways or predetermined pathway with optional errors

The observation of heterogeneous protein folding kinetics has been widely interpreted in terms of multiple independent unrelated pathways (IUP model), both experimentally and in theoretical calculations. However, direct structural information on folding intermediates and their properties now indicates that all of a protein population folds through essentially the same stepwise pathway, determined by cooperative native-like foldon units and the way that the foldons fit together in the native protein. It is essential to decide between these fundamentally different folding mechanisms. This article shows, contrary to previous supposition, that the heterogeneous folding kinetics observed for the staphylococcal nuclease protein (SNase) does not require alternative parallel pathways. SNase folding kinetics can be fit equally well by a single predetermined pathway that allows for optional misfolding errors, which are known to occur ubiquitously in protein folding. Structural, kinetic, and thermodynamic information for the folding intermediates and pathways of many proteins is consistent with the predetermined pathway-optional error (PPOE) model but contrary to the properties implied in IUP models.

NMR analysis of partially folded states and persistent structure in the alpha subunit of tryptophan synthase: implications for the equilibrium folding mechanism of a 29-kDa TIM barrel protein

Structural insights into the equilibrium folding mechanism of the alpha subunit of tryptophan synthase (alpha TS) from Escherichia coli, a (beta alpha)(8) TIM barrel protein, were obtained with a pair of complementary nuclear magnetic resonance (NMR) spectroscopic techniques. The secondary structures of rare high-energy partially folded states were probed by native-state hydrogen-exchange NMR analysis of main-chain amide hydrogens. 2D heteronuclear single quantum coherence NMR analysis of several (15)N-labeled nonpolar amino acids was used to probe the side chains involved in stabilizing a highly denatured intermediate that is devoid of secondary structure. The dynamic broadening of a subset of isoleucine and leucine side chains and the absence of protection against exchange showed that the highest energy folded state on the free-energy landscape is stabilized by a hydrophobic cluster lacking stable secondary structure. The core of this cluster, centered near the N-terminus of alpha TS, serves as a nucleus for the stabilization of what appears to be nonnative secondary structure in a marginally stable intermediate. The progressive decrease in protection against exchange from this nucleus toward both termini and from the N-termini to the C-termini of several beta-strands is best described by an ensemble of weakly coupled conformers. Comparison with previous data strongly suggests that this ensemble corresponds to a marginally stable off-pathway intermediate that arises in the first few milliseconds of folding and persists under equilibrium conditions. A second, more stable intermediate, which has an intact beta-barrel and a frayed alpha-helical shell, coexists with this marginally stable species. The conversion of the more stable intermediate to the native state of alpha TS entails the formation of a stable helical shell and completes the acquisition of the tertiary structure.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.

Importance of metastable states in the free energy landscapes of polypeptide chains

We show that the interplay between excluded volume effects, hydrophobicity, and hydrogen bonding in a tubelike representation of a polypeptide chain gives rise to free energy landscapes that, in addition to a clear global minimum, are characterized by the general presence of a small number of metastable minima, which correspond to common structural motifs observed in proteins. The complexity of the landscape increases only moderately with the length of the chain. Analysis of the temperature dependence of these landscapes reveals that the stability of specific metastable states is maximal at a temperature close to the midpoint of folding. These mestastable states are therefore likely to be of particular significance in determining the generic tendency of proteins to aggregate into potentially pathogenic agents.

Nanoscopic and redox characterization of engineered horse cytochrome C chemisorbed on a bare gold electrode

In this paper, we exploit the potential offered by site-directed mutagenesis to achieve direct adsorption of horse cyt c on a bare gold electrode surface. To this issue, the side chain T102 has been replaced by a cysteine. T102 is close to the surface exposed C-terminal residue (E104), therefore the T102C mutation is expected to generate an exposed cysteine side chain able to facilitate protein binding to the electrode via the sulphur atom (analogously to what observed for yeast iso-1-cyt c). Scanning Tunnelling and Tapping Mode Atomic Force Microscopy measurements show that the T102C mutant stably adsorbs on an Au(111) surface and retains the morphological characteristics of the native form. Cyclic voltammetry reveals that the adsorbed variant is electroactive; however, the heterogeneous electron transfer with the electrode surface is slower than that observed for yeast iso-1-cyt c. We ascribe it to differences in the tertiary architecture of the two proteins, characterized by different flexibility and stability. In particular, the region where the N- and C-terminal helices get in contact (and where the mutation occurs) is analyzed in detail, since the interactions between these two helices are considered crucial for the stability of the overall protein fold.

Deletion mutations in N-terminal alpha1 helix render heat labile enterotoxin B subunit susceptible to degradation

Heat-labile enterotoxin (LT) from enterotoxigenic Escherichia coli is a heterohexameric protein consisting of an enzymatically active A subunit, LTA, and a carrier pentameric B subunit, LTB. It is clear from the crystal structure of LTB that the N-terminal alpha1 helix lies outside the core structure. However, the function of the N-terminal alpha1 helix of LTB is unknown. The present work was carried out to investigate the effect of site-directed mutagenesis of the alpha1 helix on LTB synthesis. Six amino acids (PQSITE) located at positions 2-7 from the N terminus, including 4 aa from the alpha1 helix, were deleted by site-directed mutagenesis. The deletion resulted in complete inhibition of LTB expression in E. coli when expressed along with its signal sequence. A single amino acid deletion within the alpha1 helix also resulted in loss of expression. However, a single amino acid deletion outside the alpha1 helix did not affect LTB synthesis. Mutant proteins, whose synthesis was not detected in vivo, could be successfully translated in vitro by using the coupled transcription-translation system. Immunoblot analysis, Northern blot analysis, and in vitro transcription-translation data collectively indicate that the lack of synthesis of the mutant proteins is caused by the immediate degradation of the expressed product by cellular proteases rather than by faulty translation of mutant LTB mRNA. Coexpression of the LTA could not rescue the degradation of LTB mutants.

Branching in the sequential folding pathway of cytochrome c

Previous results indicate that the folding pathways of cytochrome c and other proteins progressively build the target native protein in a predetermined stepwise manner by the sequential formation and association of native-like foldon units. The present work used native state hydrogen exchange methods to investigate a structural anomaly in cytochrome c results that suggested the concerted folding of two segments that have little structural relationship in the native protein. The results show that the two segments, an 18-residue omega loop and a 10-residue helix, are able to unfold and refold independently, which allows a branch point in the folding pathway. The pathway that emerges assembles native-like foldon units in a linear sequential manner when prior native-like structure can template a single subsequent foldon, and optional pathway branching is seen when prior structure is able to support the folding of two different foldons.

Evidence of a two-stage thermal denaturation process in lysozyme: a Raman scattering and differential scanning calorimetry investigation

Raman spectroscopy (in the low-frequency range and the amide I band region) and modulated differential scanning calorimetry investigations have been used to analyze temperature-induced structural changes in lysozyme dissolved in 1H2O and 2H2O in the thermal denaturation process. Low-frequency Raman data reveal a change in tertiary structure without concomitant unfolding of the secondary structure. Calorimetric data show that this structural change is responsible for the configurational entropy change associated with the strong-to-fragile liquid transition and correspond to about 1/3 of the native-denaturated transition enthalpy. This is the first stage of the thermal denaturation which is a precursor of the secondary structure change and is determined to be strongly dependent on the stability of the hydrogen-bond network in water. Low-frequency Raman spectroscopy provides information on the flexibility of the tertiary structure (in the native state and the transient folding state) in relation to the fragility of the mixture. The unfolding of the secondary structure appears as a consequence of the change in the tertiary structure and independent of the solvent. Protein conformational stability is directly dependent on the stability of the native tertiary structure. The structural transformation of tertiary structure can be detected through the enhanced 1H/2H exchange inhibited in native proteins. Taking into account similar features reported in the literature observed for different proteins it can be considered that the two-stage transformation observed in lysozyme dissolved in water is a general mechanism for the thermal denaturation of proteins.

Plasticity of the protein folding landscape: switching between on- and off-pathway intermediates

Proteins may fold via parallel routes partitioned by the relative effect of solvent conditions on the relevant transition states. Thus, intermediates may or may not necessarily be obligatory species accumulated during the folding process, but rather kinetic traps due to the ruggedness of the folding landscape. Implicit in this view is the notion of plasticity of the folding pathway: proteins can be rerouted through the energy landscape by mutational, topological or solvent perturbations. Our work was specifically aimed to the experimental identification of a switch in the folding mechanism of a c-type cytochrome from the thermophilic bacterium Hydrogenobacter thermophilus (HT cyt c(552)) induced by acidic conditions. We present evidence that, by destabilizing the relevant transition state, the native state of HT cyt c(552) can be reached along alternative folding routes, which may involve an off-pathway intermediate.

Effects of serine-to-cysteine mutations on beta-lactamase folding

B. licheniformis exo-small beta-lactamase (ESBL) has two nonsequential domains and a complex architecture. We replaced ESBL serine residues 126 and 265 with cysteine to probe the conformation of buried regions in each domain. Spectroscopic, hydrodynamic, and chemical methods revealed that the mutations do not alter the native fold but distinctly change stability (S-126C > wild-type > S-126/265C > S-265C ESBL) and the features of partially folded states. The observed wild-type ESBL equilibrium intermediate has decreased fluorescence but full secondary structure. S-126C ESBL intermediate has the fluorescence of the unfolded state, no thiol reactivity, and partial secondary structure. S-265C and S-126/265C ESBL populate intermediate states unfolded by fluorescence and thiol reactivity but with full secondary structure. Mass analysis of S-126/265C ESBL in the partially folded state proved that both thiol groups become exposed simultaneously. None of the intermediates is compatible with sequential domain unfolding. Molecular dynamics simulation suggests that the stabilizing effect of the S-126C substitution is due to optimization of van der Waals interactions and packing. On the other hand, destabilization induced by the S-265C mutation results from alteration of the hydrogen-bond network. The results illustrate the large impact that seemingly conservative serine-to-cysteine changes can have on the energy landscape of proteins.

Denaturation of jack-bean urease by sodium n-dodecyl sulphate: A kinetic study below the critical micelle concentration

Kinetics of urease denaturation by anionic surfactant (sodium n-dodecyl sulphate, SDS) at concentrations below the critical micelle concentration (CMC) is investigated spectrophotometrically at neutral pH and the corresponding two-phase kinetic parameters of the process are estimated from a three-state reversible process using a binomial exponential relation based on the relaxation time method as: Using a prepared computer program, the experimental data are properly fitted into a binomial exponential relation, considering a two-phase denaturation pathway including a kinetically stable folded intermediate formed at SDS concentration of 1.1 mM. Forward and backward rate constants are estimated as: k(1)=0.2141+/-4.5 x 10(-3), k(2)=5.173 x 10(-3)+/-8.3 x 10(-5), k(-1)=0.09432+/-3.6 x 10(-4) and k(-2)=2.079 x 10(-3)+/-5.6 x 10(-5)s(-1) for the proposed mechanism. The rate-limiting step as well as the reaction coordinates in the denaturation mechanism are established. The mechanism involves formation of a kinetically stable folded native like intermediate through the electrostatic interactions. The intermediate was found to be more stable even than the native form (by about 9 kJmol(-1)) and still hexamer, because no loss of amplitude was observed. Electrophoresis experiments on the native and surfactant/urease complexes indicated a higher mobility for the kinetically folded native like intermediate.

Partially unfolded forms and non-two-state folding of a beta-sandwich: FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

FHA domains adopt a beta-sandwich fold with 11 strands. The first evidence of partially unfolded forms of a beta-sandwich is derived from native-state hydrogen exchange (NHX) of the forkhead-associated (FHA) domain from kinase-associated protein phosphatase from Arabidopsis. The folding kinetics of this FHA domain indicate that EX2 behavior prevails at pH 6.3. In the chevron plot, rollover in the folding arm and bends in the unfolding arm suggest folding intermediates. NHX of this FHA domain suggests a core of six most stable beta-strands and two loops, characterized by rare global unfolding events. Flanking this stable core are beta-strands and recognition loops with less stability, termed subglobal motifs. These suggest partially unfolded forms (near-native intermediates) with two levels of stability. The spatial separation of the subglobal motifs on the flanks suggests possible parallelism in their folding as additional beta-strands align with the stable core of six strands. Intermediates may contribute to differences in stabilities and m-values suggested by NHX or kinetics relative to chemical denaturation. Residual structure in the unfolded regime is suggested by superprotection of beta-strand 6 and by GdmCl-dependence of adjustments in amide NMR spectra and residual optical signal. The global folding stability depends strongly on pH, with at least 3 kcal/mol more stability at pH 7.3 than at pH 6.3. This FHA domain is hypothesized to fold progressively with initial hydrophobic collapse of its stable six-stranded core followed by addition of less stable flanking beta-strands and ordering of recognition loops.

Heme-peptide/protein ions and phosphorous ligands: search for site-specific addition reactions

High-resolution Fourier transform ion cyclotron resonance mass spectrometry is employed to gain thorough kinetics and thermodynamics information on the reaction of free and ligated heme-type ions with selected ligands, with the aim of obtaining an insight into the coordination environment of the prosthetic group in a variety of biomolecular ions. Adopting a stepwise approach towards systems of increasing complexity, we examined the reactivity of free gaseous iron(III) protoporphyrin IX ions, Fe(III)-heme(+), of the charged species from microperoxidase-11 (MP11) (covalently peptide bound heme), and of the multiply charged ions from heme proteins, namely, cytochrome c (cyt c) and myoglobin (examples of noncovalently protein bound hemes). Among an array of test compounds allowed to react with Fe(III)-heme(+), OP(OMe)(3) and P(OMe)(3) proved to be similarly efficient ligands in the first addition step, yet displayed markedly distinct reactivity towards heme iron already engaged in axial coordination. The ease with which P(OMe)(3) acts as a second axial ligand is exploited to probe structural and conformational features of biomolecular ions. In this way, circumstantial evidence is gained of a folded conformation of +2 charge state ions from MP11 and an elongated one for the +3 charge state ions. Similarly, both the general reaction pattern and detailed kinetics and thermodynamics data point to a regiospecific addition reaction of P(OMe)(3) directed at the heme iron within multiply charged ions from cyt c. This unprecedented example of ion-molecule reaction which specifically involves a prosthetic group belonging to protein ions stands in contrast to the multiple, nonspecific interactions established by OP(OMe)(3) molecules with the protonated sites of multiply charged cyt c and apomyoglobin ions. This finding may develop and provide sensitive probes of the structure and bonding features of protein ions in the gas phase.

Conformational Transition between Four and Five-stranded Phenylalanine Zippers Determined by a Local Packing Interaction

Alpha-helical coiled coils play a crucial role in mediating specific protein-protein interactions. However, the rules and mechanisms that govern helix-helix association in coiled coils remain incompletely understood. Here we have engineered a seven heptad "Phe-zipper" protein (Phe-14) with phenylalanine residues at all 14 hydrophobic a and d positions, and generated a further variant (Phe-14(M)) in which a single core Phe residue is substituted with Met. Phe-14 forms a discrete alpha-helical pentamer in aqueous solution, while Phe-14(M) folds into a tetrameric helical structure. X-ray crystal structures reveal that in both the tetramer and the pentamer the a and d side-chains interlock in a classical knobs-into-holes packing to produce parallel coiled-coil structures enclosing large tubular cavities. However, the presence of the Met residue in the apolar interface of the tetramer markedly alters its local coiled-coil conformation and superhelical geometry. Thus, short-range interactions involving the Met side-chain serve to preferentially select for tetramer formation, either by inhibiting a nucleation step essential for pentamer folding or by abrogating an intermediate required to form the pentamer. Although specific trigger sequences have not been clearly identified in dimeric coiled coils, higher-order coiled coils, as well as other oligomeric multi-protein complexes, may require such sequences to nucleate and direct their assembly.

The folding nucleus of the insulin superfamily: a flexible peptide model foreshadows the native state

Oxidative folding of insulin-like growth factor I (IGF-I) and single-chain insulin analogs proceeds via one- and two-disulfide intermediates. A predominant one-disulfide intermediate in each case contains the canonical A20-B19 disulfide bridge (cystines 18-61 in IGF-I and 19-85 in human proinsulin). Here, we describe a disulfide-linked peptide model of this on-pathway intermediate. One peptide fragment (19 amino acids) spans IGF-I residues 7-25 (canonical positions B8-B26 in the insulin superfamily); the other (18 amino acids) spans IGF-I residues 53-70 (positions A12-A21 and D1-D8). Containing only half of the IGF-I sequence, the disulfide-linked polypeptide (designated IGF-p) is not well ordered. Nascent helical elements corresponding to native alpha-helices are nonetheless observed at 4 degrees C. Furthermore, (13)C-edited nuclear Overhauser effects establish transient formation of a native-like partial core; no non-native nuclear Overhauser effects are observed. Together, these observations suggest that early events in the folding of insulin-related polypeptides are nucleated by a native-like molten subdomain containing Cys(A20) and Cys(B19). We propose that nascent interactions within this subdomain orient the A20 and B19 thiolates for disulfide bond formation and stabilize the one-disulfide intermediate once formed. Substitutions in the corresponding region of insulin are associated with inefficient chain combination and impaired biosynthetic expression. The intrinsic conformational propensities of a flexible disulfide-linked peptide thus define a folding nucleus, foreshadowing the structure of the native state.

Order of steps in the cytochrome C folding pathway: evidence for a sequential stabilization mechanism

Previous work used hydrogen exchange (HX) experiments in kinetic and equilibrium modes to study the reversible unfolding and refolding of cytochrome c (Cyt c) under native conditions. Accumulated results now show that Cyt c is composed of five individually cooperative folding units, called foldons, which unfold and refold as concerted units in a stepwise pathway sequence. The first three steps of the folding pathway are linear and sequential. The ordering of the last two steps has been unclear because the fast HX of the amino acid residues in these foldons has made measurement difficult. New HX experiments done under slower exchange conditions show that the final two foldons do not unfold and refold in an obligatory sequence. They unfold separately and neither unfolding obligately contains the other, as indicated by their similar unfolding surface exposure and the specific effects of destabilizing and stabilizing mutations, pH change, and oxidation state. These results taken together support a sequential stabilization mechanism in which folding occurs in the native context with prior native-like structure serving to template the stepwise formation of subsequent native-like foldon units. Where the native structure of Cyt c requires sequential folding, in the first three steps, this is found. Where structural determination is ambiguous, in the final two steps, alternative parallel folding is found.

Determination of an ensemble of structures representing the intermediate state of the bacterial immunity protein Im7

We present a detailed structural characterization of the intermediate state populated during the folding and unfolding of the bacterial immunity protein Im7. We achieve this result by incorporating a variety of experimental data available for this species in molecular dynamics simulations. First, we define the structure of the exchange-competent intermediate state of Im7 by using equilibrium hydrogen-exchange protection factors. Second, we use this ensemble to predict Phi-values and compare the results with the experimentally determined Phi-values of the kinetic refolding intermediate. Third, we predict chemical-shift measurements and compare them with the measured chemical shifts of a mutational variant of Im7 for which the kinetic folding intermediate is the most stable state populated at equilibrium. Remarkably, we found that the properties of the latter two species are predicted with high accuracy from the exchange-competent intermediate that we determined, suggesting that these three states are characterized by a similar architecture in which helices I, II, and IV are aligned in a native-like, but reorganized, manner. Furthermore, the structural ensemble that we obtained enabled us to rationalize the results of tryptophan fluorescence experiments in the WT protein and a series of mutational variants. The results show that the integration of diverse sets of experimental data at relatively low structural resolution is a powerful approach that can provide insights into the structural organization of this conformationally heterogeneous three-helix intermediate with unprecedented detail and highlight the importance of both native and non-native interactions in stabilizing its structure.

'Protein Peeling': an approach for splitting a 3D protein structure into compact fragments

MOTIVATION

The object of this study is to propose a new method to identify small compact units that compose protein three-dimensional structures. These fragments, called 'protein units (PU)', are a new level of description to well understand and analyze the organization of protein structures. The method only works from the contact probability matrix, i.e. the inter Calpha-distances translated into probabilities. It uses the principle of conventional hierarchical clustering, leading to a series of nested partitions of the 3D structure. Every step aims at dividing optimally a unit into 2 or 3 subunits according to a criterion called 'partition index' assessing the structural independence of the subunits newly defined. Moreover, an entropy-derived squared correlation R is used for assessing globally the protein structure dissection. The method is compared to other splitting algorithms and shows relevant performance.

AVAILABILITY

An Internet server with dedicated tools is available at http://www.ebgm.jussieu.fr/~gelly/