A fast conformational search strategy for finding low energy structures of model proteins.

We describe a new computer algorithm for finding low-energy conformations of proteins. It is a chain-growth method that uses a heuristic bias function to help assemble a hydrophobic core. We call it the Core-directed chain Growth method (CG). We test the CG method on several well-known literature examples of HP lattice model proteins [in which proteins are modeled as sequences of hydrophobic (H) and polar (P) monomers], ranging from 20-64 monomers in two dimensions, and up to 88-mers in three dimensions. Previous nonexhaustive methods--Monte Carlo, a Genetic Algorithm, Hydrophobic Zippers, and Contact Interactions--have been tried on these same model sequences. CG is substantially better at finding the global optima, and avoiding local optima, and it does so in comparable or shorter times. CG finds the global minimum energy of the longest HP lattice model chain for which the global optimum is known, a 3D 88-mer that has only been reachable before by the CHCC complete search method. CG has the potential advantage that it should have nonexponential scaling with chain length. We believe this is a promising method for conformational searching in protein folding algorithms.     

A simple parameter relating sequences with folding rates of small alpha helical proteins

It is found that the helix parameter (HP), which favors clustering of non-polar residues, is linearly correlated with the logarithms of rate constants of folding of small two-state alpha-helical proteins. The definition is HP = N(H)(-1) sigma [f(i)+ (f(i-1)+f(i+1))/2], where f(i)=1 or -1, if the i'th residue is hydrophobic or hydrophilic, respectively, N(H) is the number of hydrophobic residues and the summation is taken over the hydrophobic residues.     

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N‐terminal subdomain (residues 10–41) and the autonomously folding C‐terminal subdomain (residues 42–76) which pack against each other to form a structure with a unified hydrophobic core. The X‐ray structures of the isolated C‐terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs‐ms time scale dynamics of the methyl‐bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid‐state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two‐ to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.     

Influence of substrate pattern on the adsorption of HP lattice proteins

With the highly simplified hydrophobic-polar model representation of a protein, we can study essential qualitative physics without an unnecessarily large computational overhead. Using Wang-Landau sampling in conjunction with a set of efficient Monte Carlo trial moves, we studied the adsorption of short HP lattice proteins on various simple patterned substrates and in particular for checkered patterned surfaces. A set of single-site mutated HP proteins is used to investigate the role of hydrophobicity of a protein chain and surface pattern for substrates with various pattern cell sizes relative to the protein’s native configuration. For most cases, we found that the adsorption transition occurs at a lower temperature, while the hydrophobic core formation is less affected. The flattening procedure after the HP protein is adsorbed is more sensitive to the change in surface patterns and single-site mutations. These observations stay valid for both strongly and weakly attractive surfaces.     

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.     

Contributions of CDR3 to V H H domain stability and the design of monobody scaffolds for naive antibody libraries

Camelids produce functional antibodies devoid of light chains. Autonomous heavy chain variable (V(H)H) domains in these molecules have adapted to the absence of the light chain in the following ways: bulky hydrophobic residues replace small aliphatic residues in the former light chain interface, and residues from the third complementarity-determining region (CDR3) pack against the framework and stabilize the global V(H)H domain fold. To determine the specific roles of CDR3 residues in framework stabilization, we used nai;ve phage-displayed libraries, combinatorial alanine-scanning mutagenesis and biophysical characterization of purified proteins. Our results indicate that in the most stable scaffolds, the structural residues in CDR3 reside near the boundaries of the loop and pack against the framework to form a small hydrophobic core. These results allow us to differentiate between structural CDR3 residues that should remain fixed, and CDR3 residues that are tolerant to substitution and can therefore be varied to generate functional diversity within phage-displayed libraries. These methods and insights can be applied to the rapid design of heavy chain scaffolds for the identification of novel ligands using synthetic, antibody-phage libraries. In addition, they shed light on the relationships between CDR3 sequence diversity and the structural stability of the V(H)H domain fold.     

Stabilizing and destabilizing clusters in the hydrophobic core of long two-stranded alpha-helical coiled-coils

Detailed sequence analyses of the hydrophobic core residues of two long two-stranded alpha-helical coiled-coils that differ dramatically in sequence, function, and length were performed (tropomyosin of 284 residues and the coiled-coil domain of the myosin rod of 1086 residues). Three types of regions were present in the hydrophobic core of both proteins: stabilizing clusters and destabilizing clusters, defined as three or more consecutive core residues of either stabilizing (Leu, Ile, Val, Met, Phe, and Tyr) or destabilizing (Gly, Ala, Cys, Ser, Thr, Asn, Gln, Asp, Glu, His, Arg, Lys, and Trp) residues, and intervening regions that consist of both stabilizing and destabilizing residues in the hydrophobic core but no clusters. Subsequently, we designed a series of two-stranded coiled-coils to determine what defines a destabilizing cluster and varied the length of the destabilizing cluster from 3 to 7 residues to determine the length effect of the destabilizing cluster on protein stability. The results showed a dramatic destabilization, caused by a single Leu to Ala substitution, on formation of a 3-residue destabilizing cluster (DeltaT(m) of 17-21 degrees C) regardless of the stability of the coiled-coil. Any further substitution of Leu to Ala that increased the size of the destabilizing cluster to 5 or 7 hydrophobic core residues in length had little effect on stability (DeltaT(m) of 1.4-2.8 degrees C). These results suggested that the contribution of Leu to protein stability is context-dependent on whether the hydrophobe is in a stabilizing cluster or its proximity to neighboring destabilizing and stabilizing clusters.     

Defining the minimum size of a hydrophobic cluster in two-stranded alpha-helical coiled-coils: effects on protein stability

The alpha-helical coiled-coil motif is characterized by a heptad repeat pattern (abcdefg)(n) in which residues a and d form the hydrophobic core. Long coiled-coils (e.g., tropomyosin, 284 residues per polypeptide chain) typically do not have a continuous hydrophobic core of stabilizing residues, but rather one that consists of alternating clusters of stabilizing and destabilizing residues. We have arbitrarily defined a cluster as a minimum of three consecutive stabilizing or destabilizing residues in the hydrophobic core. We report here on a series of two-stranded, disulfide-bridged parallel alpha-helical coiled-coils that contain a central cassette of three consecutive hydrophobic core positions (d, a, and d) with a destabilizing cluster of three consecutive Ala residues in the hydrophobic core on each side of the cassette. The effect of adding one to three stabilizing hydrophobes in these positions (Leu or Ile; denoted as [see text]) was investigated. Alanine residues (denoted as [see text]) are used to represent destabilizing residues. The peptide with three Ala residues in the d a d cassette positions ([see text]) was among the least stable coiled-coil (T(m) = 39.3 degrees C and Urea(1/2) = 1.9 M). Surprisingly, the addition of one stabilizing hydrophobe (Leu) to the cassette or two stabilizing hydrophobes (Leu), still interspersed by an Ala in the cassette ([see text]), also did not lead to any gain in stability. However, peptides with two adjacent hydrophobes in the cassette ([see text])([see text]) did show a gain in stability of 0.9 kcal/mole over the peptide with two interspersed hydrophobes ([see text]). Because the latter three peptides have the same inherent hydrophobicity, the juxtaposition of stabilizing hydrophobes leads to a synergistic effect, and thus a clustering effect. The addition of a third stabilizing hydrophobe to the cassette ([see text]) resulted in a further synergistic gain in stability of 1.7 kcal/mole (T(m) = 54.1 degrees C and Urea(1/2) = 3.3M). Therefore, the role of hydrophobicity in the hydrophobic core of coiled-coils is extremely context dependent and clustering is an important aspect of protein folding and stability.     

Effect of solvent conditions upon refolding pathways and intermediates for a simple lattice protein

Folding pathways and intermediates for a two-dimensional lattice protein have been investigated via computer simulation at various denaturant concentrations. The protein is represented as a chain of 8 hydrophobic (H) and 12 polar (P) beads on a square lattice sequenced in such a way that the native state is a compact hydrophobic core surrounded by a shell of polar beads. Two nonbonded H beads are said to attract each other with a potential of mean force of strength ϵ. Increasing |ϵ/kT| mimics decreasing the denaturant concentration in the solution. Dynamic Monte Carlo simulations have been performed in order to investigate the folding transition and the folding pathways. Sharp folding—unfolding transitions are observed and the folding process proceeds along well-defined pathways that are populated by partially folded intermediates. The folding pathways as well as the populations of the intermediates are strongly dependent upon the denaturant concentration. Generally, intermediates containing long open stretches of H beads are more populated at high denaturant concentration, whereas compact intermediates containing a substantial number of hydrophobic contacts are more populated at low denaturant concentrations. The folding process is also observed to be cooperative in nature in that the chain does not start folding until a key fold in the middle section of the chain is formed correctly. © 1997 John Wiley & Sons, Inc. Biopoly 42: 399–409, 1997     

Histones H3/H4 form a tight complex with the inner nuclear membrane protein LBR and heterochromatin protein 1

We have recently shown that heterochromatin protein 1 (HP1) interacts with the nuclear envelope in an acetylation-dependent manner. Using purified components and in vitro assays, we now demonstrate that HP1 forms a quaternary complex with the inner nuclear membrane protein LBR and a sub-set of core histones. This complex involves histone H3/H4 oligomers, which mediate binding of LBR to HP1 and cross-link these two proteins that do not interact directly with each other. Consistent with previous observations, HP1 and LBR binding to core histones is strongly inhibited when H3/H4 are modified by recombinant CREB-binding protein, revealing a new mechanism for anchoring domains of under-acetylated chromatin to the inner nuclear membrane.     

Clustering of large hydrophobes in the hydrophobic core of two-stranded alpha-helical coiled-coils controls protein folding and stability

The de novo design and biophysical characterization of two 60-residue peptides that dimerize to fold as parallel coiled-coils with different hydrophobic core clustering is described. Our goal was to investigate whether designing coiled-coils with identical hydrophobicity but with different hydrophobic clustering of non-polar core residues (each contained 6 Leu, 3 Ile, and 7 Ala residues in the hydrophobic core) would affect helical content and protein stability. The disulfide-bridged P3 and P2 differed dramatically in alpha-helical structure in benign conditions. P3 with three hydrophobic clusters was 98% alpha-helical, whereas P2 was only 65% alpha-helical. The stability profiles of these two analogs were compared, and the enthalpy and heat capacity changes upon denaturation were determined by measuring the temperature dependence by circular dichroism spectroscopy and confirmed by differential scanning calorimetry. The results showed that P3 assembled into a stable alpha-helical two-stranded coiled-coil and exhibited a native protein-like cooperative two-state transition in thermal melting, chemical denaturation, and calorimetry experiments. Although both peptides have identical inherent hydrophobicity (the hydrophobic burial of identical non-polar residues in equivalent heptad coiled-coil positions), we found that the context dependence of an additional hydrophobic cluster dramatically increased stability of P3 (Delta Tm approximately equal to 18 degrees C and Delta[urea](1/2) approximately equal to 1.5 M) as compared with P2. These results suggested that hydrophobic clustering significantly stabilized the coiled-coil structure and may explain how long fibrous proteins like tropomyosin maintain chain integrity while accommodating polar or charged residues in regions of the protein hydrophobic core.     

Wang–Landau sampling of the interplay between surface adsorption and folding of HP lattice proteins

Generic features associated with the adsorption of proteins on solid surfaces are reviewed within the framework of the hydrophobic-polar (HP) lattice protein model. The thermodynamic behaviour and structural properties of various HP protein sequences interacting with attractive surfaces have been studied using extensive Wang–Landau sampling with different types of surfaces, each of which attracts either: all monomers, only hydrophobic (H) monomers or only polar (P) monomers, respectively. Consequently, different types of folding behaviour occur for varied surface strengths. Analysis of the combined patterns of various structural observables, e.g. the derivatives of the number of interaction contacts, together with the specific heat, leads to the identification of fundamental categories of folding and transition hierarchies. We also inferred a connection between the transition categories and the relative surface strengths, i.e. the ratios of the surface attractive strengths to the intra-chain attraction among H monomers. Thus, we believe that the folding hierarchies and identification scheme are generic for different HP sequences interacting with attractive surfaces, regardless of the chain length, sequence or surface attraction.     

Detection of secondary structure elements in proteins by hydrophobic cluster analysis.

Hydrophobic cluster analysis (HCA) is a protein sequence comparison method based on alpha-helical representations of the sequences where the size, shape and orientation of the clusters of hydrophobic residues are primarily compared. The effectiveness of HCA has been suggested to originate from its potential ability to focus on the residues forming the hydrophobic core of globular proteins. We have addressed the robustness of the bidimensional representation used for HCA in its ability to detect the regular secondary structure elements of proteins. Various parameters have been studied such as those governing cluster size and limits, the hydrophobic residues constituting the clusters as well as the potential shift of the cluster positions with respect to the position of the regular secondary structure elements. The following results have been found to support the alpha-helical bidimensional representation used in HCA: (i) there is a positive correlation (clearly above background noise) between the hydrophobic clusters and the regular secondary structure elements in proteins; (ii) the hydrophobic clusters are centred on the regular secondary structure elements; (iii) the pitch of the helical representation which gives the best correspondence is that of an alpha-helix. The correspondence between hydrophobic clusters and regular secondary structure elements suggests a way to implement variable gap penalties during the automatic alignment of protein sequences.     

Probing phenylalanine environments in oligomeric structures with pentafluorophenylalanine and cyclohexylalanine

Stabilization of protein structures and protein-protein interactions are critical in the engineering of industrially useful enzymes and in the design of pharmaceutically valuable ligands. Hydrophobic interactions involving phenylalanine residues play crucial roles in protein stability and protein-protein/peptide interactions. To establish an effective method to explore the hydrophobic environments of phenylalanine residues, we present a strategy that uses pentafluorophenylalanine (F5Phe) and cyclohexylalanine (Cha). In this study, substitution of F5Phe or Cha for three Phe residues at positions 328, 338, and 341 in the tetramerization domain of the tumor suppressor protein p53 was performed. These residues are located at the interfaces of p53-p53 interactions and are important in the stabilization of the tetrameric structure. The stability of the p53 tetrameric structure did not change significantly when F5Phe-containing peptides at positions Phe328 or Phe338 were used. In contrast, the substitution of Cha for Phe341 in the hydrophobic core enhanced the stability of the tetrameric structure with a T(m) value of 100 degrees C. Phe328 and Phe338 interact with each other through pi-interactions, whereas Phe341 is buried in the surrounding alkyl side-chains of the hydrophobic core of the p53 tetramerization domain. Furthermore, high pressure-assisted denaturation analysis indicated improvement in the occupancy of the hydrophobic core. Considerable stabilization of the p53 tetramer was achieved by filling the identified cavity in the hydrophobic core of the p53 tetramer. The results indicate the status of the Phe residues, indicating that the "pair substitution" of Cha and F5Phe is highly suitable for probing the environments of Phe residues.     

Features of a Spatially Constrained Cystine Loop in the p10 FAST Protein Ectodomain Define a New Class of Viral Fusion Peptides

The reovirus fusion-associated small transmembrane (FAST) proteins are the smallest known viral membrane fusion proteins. With ectodomains of only ∼20–40 residues, it is unclear how such diminutive fusion proteins can mediate cell-cell fusion and syncytium formation. Contained within the 40-residue ectodomain of the p10 FAST protein resides an 11-residue sequence of moderately apolar residues, termed the hydrophobic patch (HP). Previous studies indicate the p10 HP shares operational features with the fusion peptide motifs found within the enveloped virus membrane fusion proteins. Using biotinylation assays, we now report that two highly conserved cysteine residues flanking the p10 HP form an essential intramolecular disulfide bond to create a cystine loop. Mutagenic analyses revealed that both formation of the cystine loop and p10 membrane fusion activity are highly sensitive to changes in the size and spatial arrangement of amino acids within the loop. The p10 cystine loop may therefore function as a cystine noose, where fusion peptide activity is dependent on structural constraints within the noose that force solvent exposure of key hydrophobic residues. Moreover, inhibitors of cell surface thioreductase activity indicate that disruption of the disulfide bridge is important for p10-mediated membrane fusion. This is the first example of a viral fusion peptide composed of a small, spatially constrained cystine loop whose function is dependent on altered loop formation, and it suggests the p10 cystine loop represents a new class of viral fusion peptides.     

Structural determinants of antiproliferative activity of heparin on pulmonary artery smooth muscle cells

In addition to its anticoagulant properties, heparin (HP), a complex polysaccharide covalently linked to a protein core, inhibits proliferation of several cell types including pulmonary artery smooth muscle cells (PASMCs). Commercial lots of HP exhibit varying degrees of antiproliferative activity on PASMCs that may due to structural differences in the lots. Fractionation of a potent antiproliferative HP preparation into high and low molecular weight components does not alter the antiproliferative effect on PASMCs, suggesting that the size of HP is not the major determinant of this biological activity. The protein core of HP obtained by cleaving the carbohydrate-protein linkage has no growth inhibition on PASMCs, demonstrating that the antiproliferative activity resides in the glycosaminoglycan component. Basic sugar residues of glucosamine can be replaced with another basic sugar, i.e., galactosamine, without affecting growth inhibition of PASMCs. N-sulfonate groups on these sugar residues of HP are not essential for growth inhibition. However, O-sulfonate groups on both sugar residues are essential for the antiproliferative activity on PASMCs. In whole HP, in contrast to an earlier finding based on a synthetic pentasaccharide of HP, 3-O-sulfonation is not critical for the antiproliferative activity against PASMCs. The amounts and distribution of sulfonate groups on both sugar residues of the glycosaminoglycan chain are the major determinant of antiproliferative activity.     

Comparative modeling and protein-like features of hydrophobic-polar models on a two-dimensional lattice

Lattice models of proteins have been extensively used to study protein thermodynamics, folding dynamics, and evolution. Our study considers two different hydrophobic-polar (HP) models on the 2D square lattice: the purely HP model and a model where a compactness-favoring term is added. We exhaustively enumerate all the possible structures in our models and perform the study of their corresponding folds, HP arrangements in space and shapes. The two models considered differ greatly in their numbers of structures, folds, arrangements, and shapes. Despite their differences, both lattice models have distinctive protein-like features: (1) Shapes are compact in both models, especially when a compactness-favoring energy term is added. (2) The residue composition is independent of the chain length and is very close to 50% hydrophobic in both models, as we observe in real proteins. (3) Comparative modeling works well in both models, particularly in the more compact one. The fact that our models show protein-like features suggests that lattice models incorporate the fundamental physical principles of proteins. Our study supports the use of lattice models to study questions about proteins that require exactness and extensive calculations, such as protein design and evolution, which are often too complex and computationally demanding to be addressed with more detailed models.     

Protein‐spanning water networks and implications for prediction of protein–protein interactions mediated through hydrophobic effects

Hydrophobic effects, often conflated with hydrophobic forces, are implicated as major determinants in biological association and self‐assembly processes. Protein–protein interactions involved in signaling pathways in living systems are a prime example where hydrophobic effects have profound implications. In the context of protein–protein interactions, a priori knowledge of relevant binding interfaces (i.e., clusters of residues involved directly with binding interactions) is difficult. In the case of hydrophobically mediated interactions, use of hydropathy‐based methods relying on single residue hydrophobicity properties are routinely and widely used to predict propensities for such residues to be present in hydrophobic interfaces. However, recent studies suggest that consideration of hydrophobicity for single residues on a protein surface require accounting of the local environment dictated by neighboring residues and local water. In this study, we use a method derived from percolation theory to evaluate spanning water networks in the first hydration shells of a series of small proteins. We use residue‐based water density and single‐linkage clustering methods to predict hydrophobic regions of proteins; these regions are putatively involved in binding interactions. We find that this simple method is able to predict with sufficient accuracy and coverage the binding interface residues of a series of proteins. The approach is competitive with automated servers. The results of this study highlight the importance of accounting of local environment in determining the hydrophobic nature of individual residues on protein surfaces. Proteins 2014; 82:3312–3326. © 2014 Wiley Periodicals, Inc.     

Sequential mechanism of refolding of carbonic anhydrase B

The kinetics of refolding of bovine carbonic anhydrase B was studied by a variety of methods over a wide range of times (from milliseconds to hours). It has been shown that protein refolding proceeds through three stages. At the first stage (t1/2 approximately equal to 0.03 s) hydrophobic clusters and a compact state of the chain are formed. At the second stage (t1/2 approximately equal to 140 s) hydrophobic clusters are desolvated and the rigid native-like hydrophobic core is formed. At the third stage (t1/2 approximately equal to 600 s) the native active protein is formed.     

Automatic recognition of hydrophobic clusters and their correlation with protein folding units.

A method is described to objectively identify hydrophobic clusters in proteins of known structure. Clusters are found by examining a protein for compact groupings of side chains. Compact clusters contain seven or more residues, have an average of 65% hydrophobic residues, and usually occur in protein interiors. Although smaller clusters contain only side-chain moieties, larger clusters enclose significant portions of the peptide backbone in regular secondary structure. These clusters agree well with hydrophobic regions assigned by more intuitive methods and many larger clusters correlate with protein domains. These results are in striking contrast with the clustering algorithm of J. Heringa and P. Argos (1991, J Mol Biol 220:151-171). That method finds that clusters located on a protein's surface are not especially hydrophobic and average only 3-4 residues in size. Hydrophobic clusters can be correlated with experimental evidence on early folding intermediates. This correlation is optimized when clusters with less than nine hydrophobic residues are removed from the data set. This suggests that hydrophobic clusters are important in the folding process only if they have enough hydrophobic residues.