Additivity in both thermodynamic stability and thermal transition temperature for rubredoxin chimeras via hybrid native partitioning

Given any operational definition of pairwise interaction, the set of residues that differ between two structurally homologous proteins can be uniquely partitioned into subsets of clusters for which no such interactions occur between clusters. Although hybrid protein sequences that preserve such clustering are consistent with tertiary structures composed of only parental native-like interactions, the stability of such predicted structures will depend upon the physical robustness of the assumed interaction potential. A simple distance cutoff criterion was applied to the most thermostable protein known to predict such a seven-residue cluster in the metal binding site region of Pyrococcus furiosus rubredoxin and a mesophile homolog. Both conformational stability and thermal transition temperature measurements demonstrate that 39% of the differential stability arises from these seven residues.     

Residue cluster additivity of thermodynamic stability in the hydrophobic core of mesophile vs. hyperthermophile rubredoxins

The branched sidechain residues 24 and 33 in the hydrophobic core of rubredoxin differ between the Clostridium pasteurianum (Cp) and Pyrococcus furiosus (Pf) sequences. Their X-ray structures indicate that these two sidechains are in van der Waals contact with each other, while neither appears to significantly interact with the other nonconserved residues. The simultaneous interchange of residues 24 and 33 between the Cp and Pf rubredoxin sequences yield a complementary pair of hybrid proteins for which the sum of their thermodynamic stabilities equals that of the parental rubredoxins. The 1.2 kcal/mol change arising from this two residues interchange accounts for 21% of the differential thermodynamic stability between the mesophile and hyperthermophile proteins. The additional interchange of the sole nonconserved aromatic residue in the hydrophobic core yields a 0.78 kcal/mol deviation from thermodynamic additivity.     

Long- and short-range interactions in native protein structures are consistent/minimally frustrated in sequence space

We show that long- and short-range interactions in almost all protein native structures are actually consistent with each other for coarse-grained energy scales; specifically we mean the long-range inter-residue contact energies and the short-range secondary structure energies based on peptide dihedral angles, which are potentials of mean force evaluated from residue distributions observed in protein native structures. This consistency is observed at equilibrium in sequence space rather than in conformational space. Statistical ensembles of sequences are generated by exchanging residues for each of 797 protein native structures with the Metropolis method. It is shown that adding the other category of interaction to either the short- or long-range interactions decreases the means and variances of those energies for essentially all protein native structures, indicating that both interactions consistently work by more-or-less restricting sequence spaces available to one of the interactions. In addition to this consistency, independence by these interaction classes is also indicated by the fact that there are almost no correlations between them when equilibrated using both interactions and significant but small, positive correlations at equilibrium using only one of the interactions. Evidence is provided that protein native sequences can be regarded approximately as samples from the statistical ensembles of sequences with these energy scales and that all proteins have the same effective conformational temperature. Designing protein structures and sequences to be consistent and minimally frustrated among the various interactions is a most effective way to increase protein stability and foldability.     

Sequence-structure signals of 3D domain swapping in proteins

Three-dimensional domain swapping occurs when two or more identical proteins exchange identical parts of their structure to generate an oligomeric unit. It affects proteins with diverse sequences and structures, and is expected to play important roles in evolution, functional regulation and even conformational diseases. Here, we search for traces of domain swapping in the protein sequence, by means of algorithms that predict the structure and stability of proteins using database-derived potentials. Regions whose sequences are not optimal with regard to the stability of the native structure, or showing marked intrinsic preferences for non-native conformations in absence of tertiary interactions are detected in most domain-swapping proteins. These regions are often located in areas crucial in the swapping process and are likely to influence it on a kinetic or thermodynamic level. In addition, cation-pi interactions are frequently observed to zip up the edges of the interface between intertwined chains or to involve hinge loop residues, thereby modulating stability. We end by proposing a set of mutations altering the swapping propensities, whose experimental characterization would contribute to refine our in silico derived hypotheses.     

Design of combinatorial protein libraries of optimal size

In this article we introduce a computational procedure, OPTCOMB (Optimal Pattern of Tiling for COMBinatorial library design), for designing protein hybrid libraries that optimally balance library size with quality. The proposed procedure is directly applicable to oligonucleotide ligation-based protocols such as GeneReassembly, DHR, SISDC, and many more. Given a set of parental sequences and the size ranges of the parental sequence fragments, OPTCOMB determines the optimal junction points (i.e., crossover positions) and the fragment contributing parental sequences at each one of the junction points. By rationally selecting the junction points and the contributing parental sequences, the number of clashes (i.e., unfavorable interactions) in the library is systematically minimized with the aim of improving the overall library quality. Using OPTCOMB, hybrid libraries containing fragments from three different dihydrofolate reductase sequences (Escherichia coli, Bacillus subtilis, and Lactobacillus casei) are computationally designed. Notably, we find that there exists an optimal library size when both the number of clashes between the fragments composing the library and the average number of clashes per hybrid in the library are minimized. Results reveal that the best library designs typically involve complex tiling patterns of parental segments of unequal size hard to infer without relying on computational means.     

Interaction of DnaK with native proteins and membrane proteins correlates with their accessible hydrophobicity

Molecular chaperones are involved in protein folding, protein targeting to membranes, and protein renaturation after stress. They interact specifically with hydrophobic sequences that are exposed in unfolded proteins, and buried in native proteins. We have studied the interaction of DnaK with native water-soluble proteins and membrane proteins. DnaK–native protein interactions are characterized by dissociation constants between 1 and 50 μM (compared with 0.01–1 μM for unfolded proteins). This affinity is within the range of most intracellular protein concentrations, suggesting that DnaK interacts with a greater number of native proteins than previously suspected. We found a correlation between the affinity of native proteins for DnaK and their affinity for hydrophobic-interaction chromatography adsorbents, suggesting that DnaK interacts with exposed hydrophobic groups in native proteins. The interaction between DnaK and membrane proteins is characterized by DnaK's high affinity for detergent-solubilized membrane proteins, and its lower affinity for membrane proteins inserted in lipid bilayers, suggesting that the chaperone can interact with the hydrophobic sequences of the former, while it cannot penetrate the hydrophobic core of lipid bilayers. Thus, the specificity of DnaK for hydrophobic sequences is involved in its interaction with not only unfolded proteins, but also native water-soluble proteins and membrane proteins. All proteins interact with DnaK according to their exposed hydrophobicity.     

How strong are side chain interactions in the folding intermediate?

Influence of 12 nonpolar amino acids residues from the hydrophobic core of apomyoglobin on stability of its native state and folding intermediate was studied. Six of the selected residues are from the A, G and H helices; these are conserved in structure of the globin family, although nonfunctional, that is, not involved in heme binding. The rest are nonconserved hydrophobic residues that belong to the B, C, D, and E helices. Each residue was substituted by alanine, and equilibrium pH‐induced transitions in apomyoglobin and its mutants were studied by circular dichroism and fluorescent spectroscopy. The obtained results allowed estimating changes in their free energy during formation of the intermediate state. It was first shown that the strength of side chain interactions in the apomyoglobin intermediate state amounts to 15–50% of that in its native state for conserved residues, and practically to 0% for nonconserved residues. These results allow a better understanding of interactions occurring in the intermediate state and shed light on involvement of certain residues in protein folding at different stages.     

Assembly of protein tertiary structures from secondary structures using optimized potentials

We present a simulated annealing-based method for the prediction of the tertiary structures of proteins given knowledge of the secondary structure associated with each amino acid in the sequence. The backbone is represented in a detailed fashion whereas the sidechains and pairwise interactions are modeled in a simplified way, following the LINUS model of Srinivasan and Rose. A perceptron-based technique is used to optimize the interaction potentials for a training set of three proteins. For these proteins, the procedure is able to reproduce the tertiary structures to below 3 A in root mean square deviation (rmsd) from the PDB targets. We present the results of tests on twelve other proteins. For half of these, the lowest energy decoy has a rmsd from the native state below 6 A and, in 9 out of 12 cases, we obtain decoys whose rmsd from the native states are also well below 5 A.     

Failures of inverse folding and threading with gapped alignment

To calculate the tertiary structure of a protein from its amino acid sequence, the thermodynamic approach requires a potential function of sequence and conformation that has its global minimum at the native conformation for many different proteins. Here we study the behavior of such functions for the simplest model system that still has some of the features of the protein folding problem, namely two-dimensional square lattice chain configurations involving two residue types. First we show that even the given contact potential, which by definition is used to identify the folding sequences and their unique native conformations, cannot always correctly select which sequences will fold to a given structure. Second, we demonstrate that the given contact potential is not always able to favor the native alignment of a native sequence on its own native conformation over other gapped alignments of different folding sequences onto that same conformation. Because of these shortcomings, even in this simple model system in which all conformations and all native sequences are known and determined directly by the given potential, we must reexamine our expectations for empirical potentials used for inverse folding and gapped alignment on more realistic representations of proteins. © 1996 Wiley-Liss, Inc.     

Electrostatic interactions in the reconstitution of an SH2 domain from constituent peptide fragments

Fragment complementation has been used to delineate the essential recognition elements for stable folding in Src homology 2 (SH2) domains by using NMR spectroscopy, alanine scanning, and surface plasmon resonance. The unfolded 9-kD and 5-kD peptide fragments formed by limited proteolytic digestion of the N-terminal SH2 domain from the p85alpha subunit of phosphatidylinositol 3'-kinase fold into an active native-like structure on interaction with one another. The corresponding 5-kD fragment of the homologous Src protein, however, was not capable of structurally complementing the p85 9-kD fragment, indicating that fragment complementation among these SH2 domains is sensitive to the sequence differences between the Src and p85 domains. Partial complementation and folding activity could be recovered with hybrid sequences of these SH2 domains. Complete alanine scanning of the 5-kD p85 fragment was used to identify the sequence recognition elements required for complex formation. The alanine substitutions in the p85 5-kD fragment that abolished binding affinity with the cognate 9-kD fragment correlate well with highly conserved residues among SH2 domains that are either integrally involved in core packing or found at the interface between fragments. Surprisingly, however, mutation of a nonconserved surface-exposed aspartic acid to alanine was found to have a significant effect on complementation. A single additional mutation of arginine to aspartic acid allowed for recovery of native structure and increased the thermal stability of the designed Src-p85 chimera by 18 degrees C. This modification appears to relieve an unfavorable surface electrostatic interaction, demonstrating the importance of surface charge interactions in protein stability.     

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

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

Electric charge balance mechanism of extended soluble proteins

Extended proteins such as calmodulin and troponin C have two globular terminal domains linked by a central region that is exposed to water and often acts as a function-regulating element. The mechanisms that stabilize the tertiary structure of extended proteins appear to differ greatly from those of globular proteins. Identifying such differences in physical properties of amino acid sequences between extended proteins and globular proteins can provide clues useful for identification of extended proteins from complete genomes including orphan sequences. In the present study, we examined the structure and amino acid sequence of extended proteins. We found that extended proteins have a large net electric charge, high charge density, and an even balance of charge between the terminal domains, indicating that electrostatic interaction is a dominant factor in stabilization of extended proteins. Additionally, the central domain exposed to water contained many amphiphilic residues. Extended proteins can be identified from these physical properties of the tertiary structure, which can be deduced from the amino acid sequence. Analysis of physical properties of amino acid sequences can provide clues to the mechanism of protein folding. Also, structural changes in extended proteins may be caused by formation of molecular complexes. Long-range effects of electrostatic interactions also appear to play important roles in structural changes of extended proteins.     

The heat capacity of proteins

The heat capacity plays a major role in the determination of the energetics of protein folding and molecular recognition. As such, a better understanding of this thermodynamic parameter and its structural origin will provide new insights for the development of better molecular design strategies. In this paper we have analyzed the absolute heat capacity of proteins in different conformations. The results of these studies indicate that three major terms account for the absolute heat capacity of a protein: (1) one term that depends only on the primary or covalent structure of a protein and contains contributions from vibrational frequencies arising from the stretching and bending modes of each valence bond and internal rotations; (2) a term that contains the contributions of noncovalent interactions arising from secondary and tertiary structure; and (3) a term that contains the contributions of hydration. For a typical globular protein in solution the bulk of the heat capacity at 25°C is given by the covalent structure term (close to 85% of the total). The hydration term contributes about 15 and 40% to the total heat capacity of the native and unfolded states, respectively. The contribution of non-covalent structure to the total heat capacity of the native state is positive but very small and does not amount to more than 3% at 25°C. The change in heat capacity upon unfolding is primarily given by the increase in the hydration term (about 95%) and to a much lesser extent by the loss of noncovalent interactions (up to ～5%). It is demonstrated that a single universal mathematical function can be used to represent the partial molar heat capacity of the native and unfolded states of proteins in solution. This function can be experimentally written in terms of the molecular weight, the polar and apolar solvent accessible surface areas, and the total area buried from the solvent. This unique function accurately predicts the different magnitude and temperature dependences of the heat capacity of both the native and unfolded states, and therefore of the heat capacity changes associated with folding/unfolding transitions. © 1995 Wiley-Liss, Inc.     

Roles of physical interactions in determining protein-folding mechanisms: molecular simulation of protein G and alpha spectrin SH3

Protein-folding mechanisms of two small globular proteins, IgG binding domain of protein G and alpha spectrin SH3 domain are investigated via Brownian dynamics simulations with a model made of coarse-grained physical energy functions responsible for sequence-specific interactions and weak Gō-like energies. The folding pathways of alpha spectrin SH3 are known to be mainly controlled by the native topology, while protein G folding is anticipated to be more sensitive to the sequence-specific effects than native topology. We found in the folding of protein G that the C terminal beta hairpin is formed earlier and is rigid, once ordered, in the presence of an intact C terminal turn. The alpha helix is found to exhibit repeated partial formations/deformations during folding and to be stabilized via the tertiary contact with preformed beta sheets. This predicted scenario is fully consistent with experimental phi value data. Moreover, we found that the folding route is critically affected when the hydrophobic interaction is excluded from physical energy terms, suggesting that the hydrophobicity critically contributes to the folding propensity of protein G. For the folding of alpha spectrin SH3, we found that the distal beta hairpin and diverging turn are parts formed early, fully in harmony with previous results of simple Gō-like and experimental analysis, supporting that the folding route of SH3 domain is robust and coded by the native topology. The hybrid method provides useful tools for analyzing roles of physical interactions in determining folding mechanisms.     

Heterogeneity in non‐epitope loop sequence and outer membrane protein complexes alters antibody binding to the major porin protein PorB in serogroup B Neisseria meningitidis

PorB is a well‐characterized outer membrane protein that is common among Neisseria species and is required for survival. A vaccine candidate, PorB induces antibody responses that are directed against six variable surface‐exposed loops that differ in sequence depending on serotype. Although Neisseria meningitidis is naturally competent and porB genetic mosaicism provides evidence for strong positive selection, the sequences of PorB serotypes commonly associated with invasive disease are often conserved, calling into question the interaction of specific PorB loop sequences in immune engagement. In this report, we provide evidence that antibody binding to a PorB epitope can be altered by sequence mutations in non‐epitope loops. Through the construction of hybrid PorB types and PorB molecular dynamics simulations, we demonstrate that loops both adjacent and non‐adjacent to the epitope loop can enhance or diminish antibody binding, a phenotype that correlates with serum bactericidal activity. We further examine the interaction of PorB with outer membrane‐associated proteins, including PorA and RmpM. Deletion of these proteins alters the composition of PorB‐containing native complexes and reduces antibody binding and serum killing relative to the parental strain, suggesting that both intramolecular and intermolecular PorB interactions contribute to host adaptive immune evasion.     

Computational design of proteins stereochemically optimized in size, stability, and folding speed

Artificial proteins potentially barrier-free in the folding kinetics are approached computationally under the guidance of protein-folding theories. The smallest and fastest folding globular protein triple-helix-bundle (THB) is so modified as to minimize or eliminate its presumed barriers in folding speed. As the barriers may reside in the ordering of either secondary or tertiary structure, the elements of both secondary and tertiary structure in the protein are targeted for prenucleation with suitable stereochemically constrained amino acid residues. The required elements of topology and sequence for the THB are optimized independently; first the topology is optimized with simulated annealing in polypeptides of highly simplified alphabet; next, the sequence in side chains is optimized using the standard inverse design methods. The resultant three best-adapted THBs, variable in topology and distinctive in sequences, are assessed by comparing them with a few benchmark proteins. The results of mainly molecular dynamics (MD) comparisons, undertaken in explicit water at different temperatures, show that the designed sequences are favorably placed against the chosen benchmarks as THB proteins potentially thermostable in the native folds. Folding simulation experiments with MD establish that the designed sequences are rapid in the folding of individual helices, but not in the evolution of tertiary structure; energetic cum topological frustrations remain but could be the artifacts of the starting conformations that were chosen in the THBs in the folding simulations. Overall, a practical high-throughput approach for de novo protein design has been developed that may have fruitful application for any type of tertiary structure.     

Improved Gō-like models demonstrate the robustness of protein folding mechanisms towards non-native interactions

The use of simple theoretical models has provided a considerable contribution to our present understanding of the means by which proteins adopt their native fold from the plethora of available unfolded states. A common assumption in building computationally tractable models has been the neglect of stabilizing non-native interactions in the class of models described as "Gō-like." The focus of this study is the characterization of the folding of a number of proteins via a Gō-like model, which aims to map a maximal amount of information reflecting the protein sequence onto a "minimalist" skeleton. This model is shown to contain sufficient information to reproduce the folding transition states of a number of proteins, including topologically analogous proteins that fold via different transition states. Remarkably, these models also demonstrate consistency with the general features of folding transition states thought to be stabilized by non-native interactions. This suggests that native interactions are the primary determinant of most protein folding transition states, and that non-native interactions lead only to local structural perturbations. A prediction is also included for an asymmetrical folding transition state of bacteriophage lambda protein W, which has yet to be subjected to experimental characterization.     

Charge centers and formation of the protein folding core

Electrostatic interactions are important for protein folding. At low resolution, the electrostatic field of the whole molecule can be described in terms of charge center(s). To study electrostatic effects, the centers of positive and negative charge were calculated for 20 small proteins of known structure, for which hydrogen exchange cores had been determined experimentally. Two observations seem to be important. First, in all 20 proteins studied 30-100% of the residues forming hydrogen exchange core(s) were clustered around the charge centers. Moreover, in each protein more than half of the core sequences are located near the centers of charge. Second, the general architecture of all-alpha proteins from the set seems to be stabilized by interactions of residues surrounding the charge centers. In most of the alpha-beta proteins, either or both of the centers are located near a pair of consecutive strands, and this is even more characteristic for alpha/Beta and all-beta structures. Consecutive strands are very probable sites of early folding events. These two points lead to the conclusion that charge centers, defined solely from the structure of the folded protein may indicate the location of a protein's hydrogen exchange/folding core. In addition, almost all the proteins contain well-conserved continuous hydrophobic sequences of three or more residues located in the vicinity of the charge centers. These hydrophobic sequences may be primary nucleation sites for protein folding. The results suggest the following scheme for the order of events in folding: local hydrophobic nucleation, electrostatic collapse of the core, global hydrophobic collapse, and slow annealing to the native state. This analysis emphasizes the importance of treating electrostatics during protein-folding simulations.     

Protein structure and import into the peroxisomal matrix

Proteins destined for the peroxisomal matrix are synthesized in the cytosol, and imported post-translationally. It has been previously demonstrated that stably folded proteins are substrates for peroxisomal import. Mammalian peroxisomes do not contain endogenous chaperone molecules. Therefore, it is possible that proteins are required to fold into their stable, tertiary conformation in order to be imported into the peroxisome. These investigations were undertaken to determine whether proteins rendered incapable of folding were also substrates for import into peroxisomes. Reduction of albumin resulted in a less compact tertiary structure as measured by analytical centrifugation. Microinjection of unfolded albumin molecules bearing the PTS1 targeting signal resulted in their import into peroxisomes. Kinetic analysis indicated that native and unfolded molecules were imported into peroxisomes at comparable rates. While import was unaffected by treatment with cycloheximide, hsc70 molecules were observed to be imported along with the unfolded albumin molecules. These results indicate that proteins, which are incapable of assuming their native conformation, are substrates for peroxisomal import. When combined with previous observations demonstrating the import of stably folded proteins, these results support the model that tertiary structure has no effect on protein import into the peroxisomal matrix .