Lessons from the lysozyme of phage T4

An overview is presented of some of the major insights that have come from studies of the structure, stability, and folding of T4 phage lysozyme. A major purpose of this review is to provide the reader with a complete tabulation of all of the variants that have been characterized, including melting temperatures, crystallographic data, Protein Data Bank access codes, and references to the original literature. The greatest increase in melting temperature (T_m) for any point mutant is 5.1°C for the mutant Ser 117 → Val. This is achieved in part not only by hydrophobic stabilization but also by eliminating an unusually short hydrogen bond of 2.48 Å that apparently has an unfavorable van der Waals contact. Increases in T_m of more than 3–4°C for point mutants are rare, whereas several different types of destabilizing substitutions decrease T_m by 20°C or thereabouts. The energetic cost of cavity creation and its relation to the hydrophobic effect, derived from early studies of “large‐to‐small” mutants in the core of T4 lysozyme, has recently been strongly supported by related studies of the intrinsic membrane protein bacteriorhodopsin. The L99A cavity in the C‐terminal domain of the protein, which readily binds benzene and many other ligands, has been the subject of extensive study. Crystallographic evidence, together with recent NMR analysis, suggest that these ligands are admitted by a conformational change involving Helix F and its neighbors. A total of 43 nonisomorphous crystal forms of different monomeric lysozyme mutants were obtained plus three more for synthetically‐engineered dimers. Among the 43 space groups, P2₁2₁2₁ and P2₁ were observed most frequently, consistent with the prediction of Wukovitz and Yeates.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Folding of the human protein fkbp. lattice monte-carlo simulations

Monte-Carlo simulations of folding of the human protein FKBP are presented. The protein is confined in a simple cubic lattice and only nearest-neighbour interactions are considered. The evolution of protein structure, energy and diameter is followed over time. Starting from different extended conformations, compact globular forms with a hydrophobic core are reached above a critical temperature T_c, while below T_c the protein ‘freezes’ into high-energy, non-compact states. In the temperature range of folding, all the recorded intermediate states belong to two structural groups, where the process spends most of its time, separated by relatively fast transitions. During folding, the protein is successively composed of three and two compact fragments, whose separation occurs at loop positions. From comparisons performed on a domain of the family sharing 24 % identity with FKBP, it appears that the number of fragments, and therefore their location, are sequence dependent.     

Three‐dimensional change in temperature sensitivity of northern vegetation phenology

Understanding how the temperature sensitivity of phenology changes with three spatial dimensions (altitude, latitude, and longitude) is critical for the prediction of future phenological synchronization. Here we investigate the spatial pattern of temperature sensitivity of spring and autumn phenology with altitude, latitude, and longitude during 1982–2016 across mid‐ and high‐latitude Northern Hemisphere (north of 30°N). We find distinct spatial patterns of temperature sensitivity of spring phenology (hereafter “spring S_T”) among altitudinal, latitudinal, and longitudinal gradient. Spring S_T decreased with altitude mostly over eastern Europe, whereas the opposite occurs in eastern North America and the north China plain. Spring S_T decreased with latitude mainly in the boreal regions of North America, temperate Eurasia, and the arid/semi‐arid regions of Central Asia. This distribution may be related to the increased temperature variance, decreased precipitation, and radiation with latitude. Compared to spring S_T, the spatial pattern of temperature sensitivity of autumn phenology (hereafter “autumn S_T”) is more heterogeneous, only showing a clear spatial pattern of autumn S_T along the latitudinal gradient. Our results highlight the three‐dimensional view to understand the phenological response to climate change and provide new metrics for evaluating phenological models. Accordingly, establishing a dense, high‐quality three‐dimensional observation system of phenology data is necessary for enhancing our ability to both predict phenological changes under changing climatic conditions and to facilitate sustainable management of ecosystems.     

Uncover the conserved property underlying sequence‐distant and structure‐similar proteins

It is widely accepted that a protein's sequence determines its structure. The surprising finding that proteins of distant sequence can adopt similar 3D structures has raised interesting questions regarding underlying conserved properties that are essential for protein folding and stability. Uncovering the conserved properties may shed light on the folding mechanism of proteins and help with the development of computational tools for protein structure prediction. We compiled and analyzed a structure pair dataset of 66 high‐resolution and low sequence identity (16–38%) soluble proteins. Structure deviation for each pair was confirmed by calculating its C_α SiMax value and comparing its potential energy per residue. Analysis of favorable inter‐residue interactions for each structure pair indicated that the average number of inter‐residue interactions within each structure represents a conserved feature of homologous structures of distant sequence. Detailed comparison of individual types of interactions showed that the average number of either hydrophobic or hydrogen bonding interactions remains unchanged for each structure pair. These findings should be of help to improving the quality of homology models based on templates of low sequence identity, thus broadening the application of homology modeling techniques for protein studies. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 340–347, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

Visualization of the nature of protein folding by a study of a distance constraint approach in two-dimensional models

A distance constraint approach is applied to two-dimensional models of proteins in order to visualize the nature of protein folding and to examine the relative roles of different ranges of interaction. Three different native structures (I, II, and III) are considered; they have two different kinds of residues, viz., hydrophobic and hydrophilic, and different sequences of these residues. We examine how the distance constraint approach functions in the prediction of protein folding when we know the sequence of the residues, the (fixed) bond lengths, the mean distances between residues i and i + 2, and i and i + 3, and the mean distances for hydrophobic–hydrophobic, hydrophobic–hydrophilic, and hydrophilic–hydrophilic contacts between residues i and i + j, where j ≥ 4. This approach involves optimization of an object function with respect to 98 variables and is not free of the multiple-minimum problem. The optimization is always terminated if the chain is entangled and/or the segments (residues) are packed too compactly to move. In order to escape from such situations and to take the excluded-volume effect into account, a Monte Carlo method is used after the optimization is trapped in local minima. Success in the prediction of folding is found to depend on the starting conformations and on the native conformations. Fair success is obtained in predicting the helix-like structure in protein I and the overall structure of protein III, but not the β-like structures of proteins I and II. Insofar as the prediction of the structure of protein III is reasonable, it appears that some sequences of residues produce greater constraints on their conformations than others, if one considers only the hydrophobic and hydrophilic nature of the residues. These results imply that, in the folding of real proteins in three dimensions, the competition for hydrophobic (and hydrophilic) residues for inside (outside) positions in the molecule probably constitutes a necessary but not a sufficient condition to form and stabilize the native structure. The failure to predict the structure of protein II, and part of that of protein I, suggests that there are two types of long-range interactions. One (which we considered here) is nonspecific (i.e., is defined only in terms of contacts between residues of the same or different polarity) and acts at any stage of protein folding; the other (which we did not consider here) is a specific interaction between residues in pairs and contributes only when the residues in the specific pair take on the native conformation. Presumably, incorporation of such specific long-range interactions, together with the nonspecific ones, is necessary for successful protein folding, using the distance constraint approach.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of β‐hairpin structure

We previously studied a 16‐amino acid‐residue fragment of the C‐terminal β‐hairpin of the B3 domain (residues 46–61), [IG(46–61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46–61) by systematically shortening the peptide by one residue at a time from both the C‐ and the N‐terminus. To determine the structure and stability of two resulting 12‐ and 14‐amino acid‐residue peptides, IG(48–59) and IG(47–60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48–59) possesses organized three‐dimensional structure stabilized by hydrophobic interactions (Tyr50–Phe57 and Trp48–Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (T_m = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47–60) determined by DSC is T_m = 330 K and its structure is similar to that of the native β‐hairpin at all (lower) temperatures examined (283–313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a β‐hairpin structural element. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Fluorescence of proteins in aqueous neutral salt solutions. II. Influence of monovalent cation chlorides, particularly cesium chloride.

The effects of varying concentrations of monovalent cation chlorides on the fluorescence of nine proteins were studied. These are discussed in terms of “direct” or “indirect” interactions with the aromatic amino acid residues. Cs⁺ is the only cation that quenches fluorescence of proteins due to “direct” interaction with aromatic amino acid residues. Quenching is due to collisional processes. An agreement with the Stern-Volmer relationship is shown and the values of [(K_Q)eff] and [(f_a)eff] are calculated. These values confirm that the fraction of fluorescence accessible to Cs⁺ belongs to the “exposed” fluorophors. The mechanism of quenching by Cs⁺ is due to the heavy-atom effect because phosphorescence enhancement is also seen at the same time. The chlorides of Na⁺, K⁺, Rb⁺, NH₄⁺, and Li⁺ do not have a similar effect on the fluorescence of all proteins. For a given protein a gradation of the same effect (i.e., quenching or dequenching) is seen. Interactions with factors that “inderectly” affect fluorescence of any protein are involved and the structural features of the protein are responsible for such “indirect” effects. The results indicate that neutral salts can act in more than one manner. The changes in fluorescence are indicative of electrostatic and lyotropic effects of ions. Only electrostatic interactions which occur in the vicinity of tryptophan in proteins are reflected. Li⁺ shows strong interactions with proteins. In 4 M LiCl, BSA, papain, and trypsin show fluorescence changes that are indicative of changes in protein structure.     

Refolding of cytochromeb 562 and its structural stabilization by introducing a disulfide bond

The packing mechanism of the secondary structures (4-α-helices and 3₁₀-helix) of cytochromeb ₅₆₂ is simulated by the “island model,” where the formation of protein structure is accomplished by the growth-type mechanism with the driving force of packing of the long-range and specific hydrophobic interactions. Packing proceeds through the formation of the structure at the nonhelical part, where a lot of hydrophobic pairs are distributed. Consequently, conformation, nearly similar to the native one, is successfully obtained. With the help of this result, the theoretical prediction of the possibility of forming this disulfide mutant (N22C/G82C) ofb ₅₆₂ can be performed prior to the experiments by our geometrical criterion (“lampshade”). This criterion is expected to be a significant principle for introducing possible disulfide bonds into a protein to be engineered.     

Copper(II)-DNA denaturation. II. The model of DNA denaturation

In a continuing study of the denaturation of DNA as brought abought about by Cu(II) ions, results are presented for the dependence of T_m and τ (the terminal relaxation time) on ionic strength, pH, reactant concentrations, and temperature. Maximum stability of the double helix, as reflected by the longest relaxation times and highest T_m values, was observed between pH 5.3 and 6.2. Outside this range both T_m and τ decreased sharply. A second, faster relaxation time was deduced from the kinetic cureves. The apparent activation energies of the rapid and slow (“terminal”) relaxations were found to be 12 and 55 kcal/mole, respectively. Several lines of evidence led to the conclusions that (1) the rate-determining step in DNA denaturation, when occurring in the transition region, is determined by chemical events and (2) the interactions which are disrupted kinetically in the rate-determining step are those which account for the major portion of the thermal (T_m) stability of helical DNA.     

The “Pre-Molten Globule,” a New Intermediate in Protein Folding

In vitro folding studies of several proteins revealed the formation, within 2–4 msec, of transient intermediates with a large far-UV ellipticity but no amide proton protection. To solve the contradiction between the secondary structure contents estimated by these two methods, we characterized the isolated C-terminal fragment F2 of the tryptophan synthase β₂ subunit. In β₂, F2 forms its tertiary interactions with the F1 N-terminal region. Hence, in the absence of F1, isolated F2 should remain at an early folding stage with no long-range interactions. We shall show that isolated F2 folds into, and remains in, a “state” called the pre-molten globule, that indeed corresponds to a 2- to 4-msec intermediate. This condensed, but not compact, “state” corresponds to an array of conformations in rapid equilibrium comprising native as well as nonnative secondary structures. It fits the “new view” on the folding process.     

Amino acid size, charge, hydropathy indices and matrices for protein structure analysis

Background  

Prediction of protein folding and specific interactions from only the sequence (ab initio) is a major challenge in bioinformatics. It is believed that such prediction will prove possible if Anfinsen's thermodynamic principle is correct for all kinds of proteins, and all the information necessary to form a concrete 3D structure is indeed present in the sequence.     

The stability of Taq DNA polymerase results from a reduced entropic folding penalty; identification of other thermophilic proteins with similar folding thermodynamics

The thermal stability of Taq DNA polymerase is well known, and is the basis for its use in PCR. A comparative thermodynamic characterization of the large fragment domains of Taq (Klentaq) and E. coli (Klenow) DNA polymerases has been performed by obtaining full Gibbs‐Helmholtz stability curves of the free energy of folding (ΔG) versus temperature. This analysis provides the temperature dependencies of the folding enthalpy and entropy (ΔH and ΔS), and the heat capacity (ΔC_p) of folding. If increased or enhanced non‐covalent bonding in the native state is responsible for enhanced thermal stabilization of a protein, as is often proposed, then an enhanced favourable folding enthalpy should, in general, be observed for thermophilic proteins. However, for the Klenow–Klentaq homologous pair, the folding enthalpy (ΔH_fold) of Klentaq is considerably less favorable than that of Klenow at all temperatures. In contrast, it is found that Klentaq's extreme free energy of folding (ΔG_fold) originates from a significantly reduced entropic penalty of folding (ΔS_fold). Furthermore, the heat capacity changes upon folding are similar for Klenow and Klentaq. Along with this new data, comparable extended analysis of available thermodynamic data for 17 other mesophilic–thermophilic protein pairs (where enough applicable thermodynamic data exists) shows a similar pattern in seven of the 18 total systems. When analyzed with this approach, the more familiar “reduced ΔC_p mechanism” for protein thermal stabilization (observed in a different six of the 18 systems) frequently manifests as a temperature dependent shift from enthalpy driven stabilization to a reduced‐entropic‐penalty model. Proteins 2014; 82:785–793. © 2013 Wiley Periodicals, Inc.     

Dielectric study of the cooperative order–disorder transition in aqueous solutions of schizophyllan,a triple-helical polysaccharide

Schizophyllan exists in aqueous solution as a triple helix, which is intact at room temperature. Its aqueous solution forms some ordered structure at low temperatures but undergoes a sharp transition to a disordered structure as the temperature is raised. The transition temperature T_c is about 7 and 18°C for H₂O and D₂O solutions, respectively. This transition was followed by time-domain reflectometry to investigate dynamic aspects of the transition. In addition to a major peak around 10 GHz, the dielectric dispersion curve of a 20 wt % schizophyllan in D₂O exhibited a small peak around 100 MHz below T_c and around 10 MHz above T_c. The major peak is due to bulk water, whereas the 100 MH_z peak is assigned to “bound” or “structured” water, and that around 10 MHz to side-chain glucose residues. However, unlike usual bound water reported for biopolymer solutions, this “structured” water disappears abruptly when the temperature becomes close to T_c without accompanying a conformational transition of the main chain. The above assignment is consistent with the structure of the ordered phase derived from previous static data that it consists of side-chain glucose residues along with nearby water molecules surrounding the helix core that are interacting with each other loosely through hydrogen bonds, and spreads radially only a layer of one or two water molecules but a long distance along the helix axis. © 1995 John Wiley & Sons, Inc.     

A Stoichiometry Driven Universal Spatial Organization of Backbones of Folded Proteins: Are there Chargaff's Rules for Protein Folding?

Abstract

Protein folding is at least a six decade old problem, since the times of Pauling and Anfinsen. However, rules of protein folding remain elusive till date. In this work, rigorous analyses of several thousand crystal structures of folded proteins reveal a surprisingly simple unifying principle of backbone organization in protein folding. We find that protein folding is a direct consequence of a narrow band of stoichiometric occurrences of amino-acids in primary sequences, regardless of the size and the fold of a protein. We observe that “preferential interactions” between amino-acids do not drive protein folding, contrary to all prevalent views. We dedicate our discovery to the seminal contribution of Chargaff which was one of the major keys to elucidation of the stoichiometry-driven spatially organized double helical structure of DNA.     

An empirical energy function for threading protein sequence through the folding motif

In this paper we present a new residue contact potantial derived by statistical analysis of protein crystal structures. This gives mean hydrophobic and pairwise contact energies as a function of residue type and distance interval. To test the accuracy of this potential we generate model structures by “threading” different sequences through backbone folding motifs found in the structural data base. We find that conformational energies calculated by summing contact potentials show perfect specificity in matching the correct sequences with each globular folding motif in a 161-protcin data set. They also identify correct models with the core folding motifs of heme-rythrin and immunoglobulin McPC603 V₁-do- main, among millions of alternatives possible when we align subsequences with α-helices and β-strands, and allow for variation in the lengths of intervening loops. We suggest that contact potentials reflect important constraints on nonbonded interaction in native proteins, and that “threading” may be useful for structure prediction by recognition of folding motif. © 1993 Wiley-Liss, Inc.     

Folding of S6 structures with divergent amino acid composition: pathway flexibility within partly overlapping foldons

Studies of circular permutants have demonstrated that the folding reaction of S6 from Thermus thermophilus (S6_T) is malleable and responds in an ordered manner to changes of the sequence separation between interacting residues: the S6_T permutants retain a common nucleation pattern in the form of a two-strand-helix motif that can be recruited from different parts of the structure. To further test the robustness of the two-strand-helix nucleus we have here determined the crystal structure and folding reaction of an evolutionary divergent S6 protein from the hyperthermophilic bacterium Aquifex aeolicus (S6_A). Although the overall topology of S6_A is very similar to that of S6_T the architecture of the hydrophobic core is radically different by containing a large proportion of stacked Phe side-chains. Despite this disparate core composition, the folding rate constant and the kinetic m values of S6_A are identical to those of S6_T. The folding nucleus of S6_A is also found to retain the characteristic two-strand-helix motif of the S6_T permutants, but with a new structural emphasis. The results suggest that the protein folding reaction is linked to topology only in the sense that the native-state topology determines the repertoire of accessible nucleation motifs. If the native structure allows several equivalent ways of recruiting a productive nucleus the folding reaction is free to redistribute within these topological constraints.     

Structural energetics of the molten globule state

Certain partly ordered protein conformations, commonly called “moltenglobule states,” are widely believed to represent protein folding intermediates. Recentstructural studies of molten globule states ofdifferent proteins have revealed features whichappear to be general in scope. The emergingconsensus is that these partly ordered forms exhibit a high content of secondary structure, considerable compactness, nonspecific tertiary structure, and significant structural flexibility. These characteristics may be used to define ageneral state of protein folding called “the molten globule state,” which is structurally andthermodynamically distinct from both the native state and the denatured state. Despite exaatensive knowledge of structural features of afew molten globule states, a cogent thermodynamic argument for their stability has not yetbeen advanced. The prevailing opinion of thelast decade was that there is little or no enthalpy difference or heat capacity differencebetween the molten globule state and the unfolded state. This view, however, appears to beat variance with the existing database of protein structural energetics and with recent estimates of the energetics of denaturation of α-lactalbumin, cytochrome c, apomyoglobin, and T4 lysozyme. We discuss these four proteins at length. The results of structural studies, together with the existing thermodynamic values for fundamental interactions in proteins, provide the foundation for a structural thermodynamic framework which can account for the observed behavior of molten globule states. Within this framework, we analyze the physical basis for both the high stability of several molten globule states and the low probability of other protential folding intermediates. Additionally, we consider, in terms of reduced enthalpy changes and disrupted cooperative interactions, the thermodynamic basis for the apparent absence of a thermally induced, cooperative unfolding transition for some molten globule states. © 1993 Wiley-Liss, Inc.