Mutations in alpha-helical solvent-exposed sites of eglin c have long-range effects: evidence from molecular dynamics simulations

Eglin c is a small protease inhibitor whose structural and thermodynamic properties have been well studied. Previous thermodynamic measurements on mutants at solvent-accessible positions in the protein's helix have shown the unexpected result that the data could be best fit by the inclusion of residue- and position-specific parameters to the model. To explore the origins of this surprising result, long molecular dynamics simulations in explicit solvent have been performed. These simulations indicate specific long-range interactions between the solvent-exposed residues in the eglin c alpha-helix and binding loop, an unexpected observation for such a small protein. The residues involved in the interaction are on opposite sides of the protein, about 25 A apart. Simulations of alanine substitutions at the solvent-exposed helix positions, arginine 22, glutamic acid 23, threonine 26, and leucine 27, show both small and large perturbations of eglin c dynamics. Two mutations exhibit large impacts on the long-range helix-loop interactions. Previous stability measurements (Yi et al., Biochemistry 2003;42:7594-7603) had indicated that an alanine substitution at position 27 was less stabilizing than at other solvent-exposed positions in the helix. The L27A mutation effects observed in these simulations suggest that the position-dependent loss of stability measured in wet bench experiments is derived from changes in dynamics that involve long-range interactions; thus, these simulations support the hypothesis that solvent-exposed positions in helices are not always equivalent.     

Dynamic coupling and allosteric behavior in a nonallosteric protein

Long-range intraprotein interactions give rise to many important protein behaviors. Understanding how energy is transduced through protein structures to either transmit a signal or elicit conformational changes is therefore a current challenge in structural biology. In an effort to understand such linkages, multiple V --> A mutations were made in the small globular protein eglin c. The physical responses, as mapped by NMR spin relaxation, residual dipolar couplings (RDCs), and scalar couplings, illustrate that the interior of this nonallosteric protein forms a dynamic network and that local perturbations are transmitted as dynamic and structural changes to distal sites as far as 16 A away. Two basic types of propagation responses were observed: contiguous pathways of enhanced (attenuated) dynamics with no change in structure; and dispersed (noncontiguous) changes in methyl rotation rates that appear to result from subtle deformation of backbone structure. In addition, energy transmission is found to be unidirectional. In one mutant, an allosteric conformational change of a side chain is seen in the context of a pathway of propagated changes in picosecond to nanosecond dynamics. The observation of so many long-range interactions in a small, rigid system lends experimental weight to the idea that all well-folded proteins inherently possess allosteric features [Gunasekaran et al. (2004) Proteins 57, 433-443], and that dynamics are a rich source of information for mapping and gaining mechanistic insight into communication pathways in individual proteins.     

Monitoring aromatic picosecond to nanosecond dynamics in proteins via 13C relaxation: expanding perturbation mapping of the rigidifying core mutation, V54A, in eglin c

Long-range effects, such as allostery, have evolved in proteins as a means of regulating function via communication between distal sites. An NMR-based perturbation mapping approach was used to more completely probe the dynamic response of the core mutation V54A in the protein eglin c by monitoring changes in picosecond to nanosecond aromatic side-chain dynamics and H/D exchange stabilities. Previous side-chain dynamics studies on this mutant were limited to methyl-bearing residues, most of which were found to rigidify on the picosecond to nanosecond time scale in the form of a contiguous "network". Here, high precision (13)C relaxation data from 13 aromatic side chains were acquired by applying canonical relaxation experiments to a newly developed carbon labeling scheme [Teilum et al. (2006) J. Am. Chem. Soc. 128, 2506-2507]. The fitting of model-free parameters yielded S (2) variability which is intermediate with respect to backbone and methyl-bearing side-chain variability and tau e values that are approximately 1 ns. Inclusion of the aromatic dynamic response results in an expanded network of dynamically coupled residues, with some aromatics showing increases in flexibility, which partially offsets the rigidification in methyl side chains. Using amide hydrogen exchange, dynamic propagation on a slower time scale was probed in response to the V54A perturbation. Surprisingly, regional stabilization (slowed exchange) 10-12 A from the site of mutation was observed despite a global destabilization of 1.5 kcal x mol (-1). Furthermore, this unlikely pocket of stabilized residues colocalizes with increases in aromatic flexibility on the faster time scale. Because the converse is also true (destabilized residues colocalize with rigidification on the fast time scale), a plausible entropy-driven mechanism is discussed for relating colocalization of opposing dynamic trends on vastly different time scales.     

Crystal and molecular structure of the bovine alpha-chymotrypsin-eglin c complex at 2.0 A resolution.

The crystal structure of the complex between bovine alpha-chymotrypsin and the leech (Hirudo medicinalis) protein proteinase inhibitor eglin c has been refined at 2.0 A resolution to a crystallographic R-factor of 0.167. The structure of the complex includes 2290 protein and 143 solvent atoms. Eglin c is bound to the cognate enzyme through interactions involving 11 residues of the inhibitor (sites P5-P4' in the reactive site loop, P10' and P23') and 17 residues from chymotrypsin. Binding of eglin c to the enzyme causes a contained hinge-bending movement around residues P4 and P4' of the inhibitor. The tertiary structure of chymotrypsin is little affected, with the exception of the 10-13 region, where an ordered structure for the polypeptide chain is observed. The overall binding mode is consistent with those found in other serine proteinase-protein-inhibitor complexes, including those from different inhibition families. Contained, but significant differences are observed in the establishment of intramolecular hydrogen bonds and polar interactions stabilizing the structure of the intact inhibitor, if the structure of eglin c in its complex with chymotrypsin is compared with that of other eglin c-serine proteinase complexes.     

Hydration-coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin

The influence of hydration on the internal dynamics of a typical EF-hand calciprotein, parvalbumin, was investigated by incoherent quasi-elastic neutron scattering (IQNS) and solid-state 13C-NMR spectroscopy using the powdered protein at different hydration levels. Both approaches establish an increase in protein dynamics upon progressive hydration above a threshold that only corresponds to partial coverage of the protein surface by the water molecules. Selective motions are apparent by NMR in the 10-ns time scale at the level of the polar lysyl side chains (externally located), as well as of more internally located side chains (from Ala and Ile), whereas IQNS monitors diffusive motions of hydrogen atoms in the protein at time scales up to 20 ps. Hydration-induced dynamics at the level of the abundant lysyl residues mainly involve the ammonium extremity of the side chain, as shown by NMR. The combined results suggest that peripheral water-protein interactions influence the protein dynamics in a global manner. There is a progressive induction of mobility at increasing hydration from the periphery toward the protein interior. This study gives a microscopic view of the structural and dynamic events following the hydration of a globular protein.     

Hydrophobic core mutations in CI2 globally perturb fast side-chain dynamics similarly without regard to position

Protein dynamics is currently an area of intense research because of its importance as complementary information to the huge quantity of available data relating protein structure and function. Because it is usually the influence of dynamics on function that is studied, the physical determinants of the distribution of flexibility in proteins have not been explored as thoroughly. In the present NMR study, an expanded suite of five (2)H relaxation experiments was used to characterize the picosecond-to-nanosecond side-chain dynamics of chymotrypsin inhibitor 2 (CI2) and five hydrophobic core mutants, some of which are members of the folding nucleus. Because CI2 is a homologue of the serine protease inhibitor eglin c, which has already been extensively characterized in terms of its dynamics, it was possible to compare not only side-chain dynamics but also the responses of these dynamics to analogous mutations. Remarkably, each of the five core mutations in CI2 led to similar and reproducible increases in side-chain flexibility throughout the entire structure. Although the expanded suite of (2)H relaxation experiments does not affect model selection for the vast majority of residues, it did enable the detection of increasing levels of nanosecond-scale motions in CI2's reactive site binding loop as the L68 side chain was progressively shortened by mutation. Collectively, we observed that the CI2 mutants are more dynamically similar to each other than to the more rigid wild-type CI2, from which we propose that wild-type CI2 has been optimized to a specific level of rigidity which may aid in its function as a serine protease inhibitor. We also observed that the pattern of side-chain dynamics of CI2 is quantitatively similar to eglin c, but that this similarity is lost upon mutating both proteins at an equivalent position. Finally, (15)N relaxation was used to characterize the backbone dynamics of wild-type and mutant CI2. Interestingly, mutation at folding nucleus positions led to widespread increases in backbone flexibility, whereas non-folding-nucleus positions led to increases in flexibility in the C-terminal half of the protein only.     

Missense mutations and evolutionary conservation of amino acids: evidence that many of the amino acids in factor IX function as "spacer" elements.

We report 31 point mutations in the factor IX gene and explore the relationship between the level of evolutionary conservation of an amino acid and the probability of a mutation causing hemophilia B. From our total sample of 125 hemophiliacs and from those reported by others, we identify 95 independent missense mutations, 94 of which occur at amino acids that are evolutionarily conserved in the available mammalian factor IX sequences. The likelihood of a missense mutation causing hemophilia B depends on whether the residue is also conserved in the factor IX-related proteases: factor VII, factor X, and protein C. Most of the possible missense mutations in generically conserved residues (i.e., those conserved in factor IX and in all the related proteases) should cause disease. In contrast, missense mutations in factor IX-specific residues (i.e., those conserved in human, cow, dog, and mouse factor IX but not in the related proteases) are sixfold less likely to cause disease. Missense mutations at nonconserved residues are 33-fold less likely to cause disease. At least three models are compatible with these observations. A comparison of sequence alignments from four and nine species of factor IX and an examination of the missense mutations occurring at CpG residues suggest a model in which most residues fall on opposite ends of a spectrum. In about 40% of residues, virtually any missense mutation in a minority of the residues will cause disease, while virtually no missense mutations will cause disease in most of the remaining residues. Thus, many of the residues in factor IX are spacers; that is, the main chains are presumably necessary to keep other amino acid interactions in register, but the nature of the side chain is unimportant.     

Backbone and side chain dynamics of mutant calmodulin-peptide complexes

The mechanism of long-range coupling of allosteric sites in calcium-saturated calmodulin (CaM) has been explored by characterizing structural and dynamics effects of mutants of calmodulin in complex with a peptide corresponding to the smooth muscle myosin light chain kinase calmodulin-binding domain (smMLCKp). Four CaM mutants were examined: D95N and D58N, located in Ca2+-binding loops; and M124L and E84K, located in the target domain-binding site of CaM. Three of these mutants have altered allosteric coupling either between Ca2+-binding sites (D58N and D95N) or between the target- and Ca2+-binding sites (E84K). The structure and dynamics of the mutant calmodulins in complex with smMLCKp were characterized using solution NMR. Analysis of chemical shift perturbations was employed to detect largely structural perturbations. 15N and 2H relaxation was employed to detect perturbations of the dynamics of the backbone and methyl-bearing side chains of calmodulin. The least median squares method was found to be robust in the detection of perturbed sites. The main chain dynamics of calmodulin are found to be largely unresponsive to the mutations. Three mutants show significantly perturbed dynamics of methyl-bearing side chains. Despite the pseudosymmetric location of Ca2+-binding loop mutations D58N and D95N, the dynamic response of CaM is asymmetric, producing long-range perturbation in D58N and almost none in D95N. The mutations located at the target domain-binding site have quite different effects. For M124L, a local perturbation of the methyl dynamics is observed, while the E84K mutation produces a long-range propagation of dynamic perturbations along the target domain-binding site.     

Determinants of the factor IX mutational spectrum in haemophilia B: an analysis of missense mutations using a multi-domain molecular model of the activated protein

A multi-domain molecular model of factor IXa was constructed by comparative methods. The quaternary structure of the protein was assembled by docking individual domains through consideration of their shape complementarity, polaric properties and the location of cross-reacting material positive/negative (CRM^+/–) variants on domain surfaces. Some 217 different missense mutations in the factor IX (F9) gene were then selected for study. Using maximum likelihood analysis, missense mutations affecting highly conserved amino acid residues of factor IX were shown to be 15–20 times more likely to result in haemophilia B than those affecting non-conserved residues. However, about one quarter of this increase in likelihood of clinical observation could be attributed to the magnitude of the amino acid exchange. Missense mutations in structurally conserved residues were found to be 2.1-fold more likely to come to clinical attention than those in structurally variable residues. Missense mutations in residues whose side chains were inwardly pointing were 3.6-fold more likely to be observed than those in surface residues. These observations imply a complex hierarchy of sequence/structure conservation in the protein. The severity of the clinical phenotype correlated with both the extent of the evolutionary sequence conservation of the residue at the site of mutation and the magnitude of the amino acid exchange. Further, the substitution of residues exhibiting minimal side chain solvent accessibility was associated disproportionately with severe haemophilia compared with that of surface residues. Clusters of CRM⁺ mutations were observed at factor IX-specific residues on the surface of the molecule. These clusters may reflect factor IX-specific docking interactions. The likelihood that a given factor IX mutation will come to clinical attention is therefore a complex function of the sequence characteristics of the F9 gene, the nature of the amino acid substitution, its precise location and immediate environment within the protein molecule, and its resulting effects on the structure and function of the protein.This paper is dedicated to the memory of Andrew Wacey     

The role of packing interactions in stabilizing folded proteins.

In order to investigate the role of nonpolar side chains in determining protein stability, we have carried out a molecular dynamics simulation study of the thermodynamics of interconverting isoleucine and valine side chains in the core of ribonuclease T1. The free energy change in the unfolded state, which we take to be fully solvated, was small and agrees qualitatively with experimental studies of alkane solvation. In the two Ile----Val mutations studied, the protein was able to relax around the smaller side chains, while in the case of the two Val----Ile mutations, the ability of the core to accommodate the extra methylene group depended on where the mutation took place. We argue that the experimentally observed decrease in stability for mutating isoleucine into valine results from a loss of favorable packing interactions of the side chain in the folded form of the protein. This supports the view that packing interactions in the folded state are an important contributor to the overall stability of the folded protein and that the core of the native protein is packed efficiently and almost completely.     

Structural assessment of glycyl mutations in invariantly conserved motifs

Motifs that are evolutionarily conserved in proteins are crucial to their structure and function. In one of our earlier studies, we demonstrated that the conserved motifs occurring invariantly across several organisms could act as structural determinants of the proteins. We observed the abundance of glycyl residues in these invariantly conserved motifs. The role of glycyl residues in highly conserved motifs has not been studied extensively. Thus, it would be interesting to examine the structural perturbations induced by mutation in these conserved glycyl sites. In this work, we selected a representative set of invariant signature (IS) peptides for which both the PDB structure and mutation information was available. We thoroughly analyzed the conformational features of the glycyl sites and their local interactions with the surrounding residues. Using Ramachandran angles, we showed that the glycyl residues occurring in these IS peptides, which have undergone mutation, occurred more often in the L-disallowed as compared with the L-allowed region of the Ramachandran plot. Short range contacts around the mutation site were analyzed to study the steric effects. With the results obtained from our analysis, we hypothesize that any change of activity arising because of such mutations must be attributed to the long-range interaction(s) of the new residue if the glycyl residue in the IS peptide occurred in the L-allowed region of the Ramachandran plot. However, the mutation of those conserved glycyl residues that occurred in the L-disallowed region of the Ramachandran plot might lead to an altered activity of the protein as a result of an altered conformation of the backbone in the immediate vicinity of the glycyl residue, in addition to long range effects arising from the long side chains of the new residue. Thus, the loss of activity because of mutation in the conserved glycyl site might either relate to long range interactions or to local perturbations around the site depending upon the conformational preference of the glycyl residue.     

X-ray crystal structure of the serine proteinase inhibitor eglin c at 1.95 A resolution.

The crystal structure of eglin c, naturally occurring in the leech Hirudo medicinalis, is known from its complexes with various serine proteinases, but the crystallization of free eglin c has not yet been reported. A method is described for growing well-diffracting crystals of free eglin c from highly concentrated protein solutions (approximately 200 mg/ml). The space group of the orthorhombic crystals was determined to be P2(1)2(1)2(1) with unit cell parameters a = 32.6, b = 42.0, c = 44.1 A. The structure of free eglin c was resolved at 1.95 A resolution by Patterson search methods. The final model contains all 70 amino acids of eglin c and 125 water molecules. In comparison to the eglin structure known from its complexes with proteinases, only small differences have been observed in free eglin c. However, the reactive site-binding loop and a few residues on the surface of eglin have been found in different conformations due to crystal contacts. In contrast to the complex structures, the first seven amino acids of the highly flexible amino terminus can be located. Crystallographic refinement comprised molecular dynamics refinement, classical restrained least-squares refinement and individual isotropic atomic temperature refinement. The final R-factor is 15.8%.     

Probing subtle conformational changes induced by phosphorylation and point mutations in the TIR domains of TLR2 and TLR3

Extensive research performed on Toll‐like receptor (TLR) signaling has identified residues in the Toll/interleukin‐1 receptor (TIR) domains that are essential for its proper functioning. Among these residues, those in BB loop are particularly significant as single amino acid mutations in this region can cause drastic changes in downstream signaling. However, while the effect of these mutations on the function is well studied (like the P681H mutation in TLR2, the A795P mutation in TLR3, and the P714H mutation in TLR4), their influence on the dynamics and inter‐residue networks is not well understood. The effects of local perturbations induced by these mutations could propagate throughout the TIR domain, influencing interactions with other TIR domain‐containing proteins. The identification of these subtle changes in inter‐residue interactions can provide new insights and structural rationale for how single‐point mutations cause drastic changes in TIR–TIR interactions. We employed molecular dynamics simulations and protein structure network (PSN) analyses to investigate the structural transitions with special emphasis on TLR2 and TLR3. Our results reveal that phosphorylation of the Tyr 759 residue in the TIR domain of TLR3 introduces rigidity to its BB loop. Subtle differences in the intra BB loop hydrogen bonding network between TLR3 and TLR2 are also observed. The PSN analyses indicate that the TIR domain is highly connected and pinpoints key differences in the inter‐residue interactions between the wild‐type and mutant TIR domains, suggesting that TIR domain structure is prone to allosteric effects, consistent with the current view of the influence of allostery on TLR signaling.     

Conservative mutation Met8 --> Leu affects the folding process and structural stability of squash trypsin inhibitor CMTI-I

Protein molecules can accommodate a large number of mutations without noticeable effects on their stability and folding kinetics. On the other hand, some mutations can have quite strong effects on protein conformational properties. Such mutations either destabilize secondary structures, e.g., alpha-helices, are incompatible with close packing of protein hydrophobic cores, or lead to disruption of some specific interactions such as disulfide cross links, salt bridges, hydrogen bonds, or aromatic-aromatic contacts. The Met8 --> Leu mutation in CMTI-I results in significant destabilization of the protein structure. This effect could hardly be expected since the mutation is highly conservative, and the side chain of residue 8 is situated on the protein surface. We show that the protein destabilization is caused by rearrangement of a hydrophobic cluster formed by side chains of residues 8, Ile6, and Leu17 that leads to partial breaking of a hydrogen bond formed by the amide group of Leu17 with water and to a reduction of a hydrophobic surface buried within the cluster. The mutation perturbs also the protein folding. In aerobic conditions the reduced wild-type protein folds effectively into its native structure, whereas more then 75% of the mutant molecules are trapped in various misfolded species. The main conclusion of this work is that conservative mutations of hydrophobic residues can destabilize a protein structure even if these residues are situated on the protein surface and partially accessible to water. Structural rearrangement of small hydrophobic clusters formed by such residues can lead to local changes in protein hydration, and consequently, can affect considerably protein stability and folding process.     

Free energy simulation to investigate the effect of amino acid sequence environment on the severity of osteogenesis imperfecta by glycine mutations in collagen

Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Gly --> Ser mutations causing two different osteogenesis imperfecta phenotypes. These simulations were performed to investigate the impact of local amino acid sequence environment adjacent to a mutation site on the stability of the collagen. The average free energy differences for a Gly --> Ser mutant relative to a wild type are 3.4 kcal/mol and 8.2 kcal/mol for a nonlethal site and a lethal site, respectively. The free energy change differences of mutant containing two Ser residues relative to the wild type at the nonlethal and lethal mutation sites are 4.6 and 9.8 kcal/mol, respectively. Although electrostatic interactions stabilize mutants containing one or two Ser residues at both mutation sites, van der Waals interactions are of sufficient magnitude to cause a net destabilization. The presence of Gln and Arg near the mutation site, which contain large and polar side chains, provide more destabilization than amino acids containing small and nonpolar side chains.     

Consequences of single-site mutations in the intestinal fatty acid binding protein

The intestinal fatty acid binding protein (IFABP) is a small (15 kDa) protein consisting mostly of 10 antiparallel beta-strands (A-J) and a small helical region that serves as a portal for the ligand. Two beta-sheet structures (strands A-E and F-J) surround a cavity into which the ligand binds. In this work, we investigated how changes in the side chains of specific residues are propagated through the structure. To determine what these changes were and how they relate to changes in stability, (15)N chemical shift perturbations were measured and compared to those of the wild-type protein. Seven mutations, five of which change either valine or leucine to glycine, have been examined. All these mutants were less stable than wild-type IFABP, suggesting some structural changes. For five of the mutants, the data suggest that destabilization of a small region of the protein propagates throughout the structure, resulting in an overall decrease in stability. In two (Leu38Gly and Leu89Gly), the loss of cooperativity in the equilibrium denaturation curves suggests that the destabilization of one region may not be transmitted to other regions in a cooperative manner. It is shown that the effect of mutating hydrophobic residues is much greater than that observed upon mutation of a solvent-exposed polar residue.     

Severity of osteogenesis imperfecta and structure of a collagen-like peptide modeling a lethal mutation site

We show that there are correlations between the severities of osteogenesis imperfecta (OI) phenotypes and changes in the residues near the mutation site. Our results show the correlations between the severity of various forms of the inherited disease OI and alteration of residues near the site of OI causing mutations. Among our many observed correlations are particularly striking ones between the presence of nearby proline residues and lethal mutations, and the presence of nearby alanines residues and nonlethal mutations. We investigated the possibility that these correlations have a structural basis using molecular dynamics simulations of collagen-like molecules designed to mimic the site of a lethal OI mutation in collagen type I. Our significant finding is that interchain hydrogen bonding is greatly affected by variations in residue type. We found that the strength of hydrogen bond networks between backbone atoms on different chains depends on the local residue sequence and is weaker in proline-rich regions of the molecule. We also found that an alanine at a site near an OI mutation causes less structural disruption than a proline, and that residue side chains also form interchain hydrogen bonds with frequencies that are dependent on residue type. For example, arginine side chains form strong hydrogen bonds with the backbone of the subsequent peptide chain, while lysine and glutamine less frequently form similar hydrogen bonds. This decrease in the observed hydrogen bond frequency correlates with a decrease in the experimentally determined thermal stability. We contrasted general structural properties of model collagen peptides with and without the mutation to examine the effect of the single-point mutation on the surrounding residues.     

Nucleosome Interaction Surface of Linker Histone H1c Is Distinct from That of H10

The fully organized structure of the eukaryotic nucleosome remains unsolved, in part due to limited information regarding the binding site of the H1 or linker histone. The central globular domain of H1 is believed to interact with the nucleosome core at or near the dyad and to bind at least two strands of DNA. We utilized site-directed mutagenesis and in vivo photobleaching to identify residues that contribute to the binding of the globular domain of the somatic H1 subtype H1c to the nucleosome. As was previously observed for the H1⁰ subtype, the binding residues for H1c are clustered on the surface of one face of the domain. Despite considerable structural conservation between the globular domains of these two subtypes, the locations of the binding sites identified for H1c are distinct from those of H1⁰. We suggest that the globular domains of these two linker histone subtypes will bind to the nucleosome with distinct orientations that may contribute to higher order chromatin structure heterogeneity or to differences in dynamic interactions with other DNA or chromatin-binding proteins.     

Polar residues in membrane domains of proteins: molecular basis for helix-helix association in a mutant CFTR transmembrane segment

Polar side chains constitute over 20% of residues in the transmembrane (TM) helices of membrane proteins, where they may serve as hydrogen bond interaction sites for phenotypic polar mutations that arise in membrane protein-related diseases. To systematically explore the structural consequences of H-bonds between TM helices, we focused on TM4 of the cystic fibrosis conductance regulator (CFTR) and its cystic fibrosis- (CF-) phenotypic mutation, V232D, as a model system. Synthetic peptides corresponding to wild-type (TM4-wt) (residues 219-242: LQASAFCGLGFLIVLALFQAGLGR) and mutant (TM4-V232D) sequences both adopt helical structures in SDS micelles and display dimer bands on SDS-PAGE arising from disulfide bond formation via wild-type residue Cys-225. However, the TM4-V232D peptide additionally forms a ladder of noncovalent oligomers, including tetramers, hexamers, and octamers, mediated by a hydrogen bond network involving Asp-Gln side chain-side chain interactions. Ala-scanning mutagenesis of the TM4 sequence indicated that ladder formation minimally required the simultaneous presence of the Cys-225, Asp-232, and Gln-237 residues. As random hydrophobic sequences containing these three residues at TM4 equivalent positions did not oligomerize, specific van der Waals packing interactions between helix side chains were also shown to play a crucial role. Overall, the results suggest that polar mutations in membrane domains, in conjunction with critically positioned polar partner residues, potentially constitute a source of aberrant helix interactions that could contribute to loss of function when they arise in protein transmembrane domains.     

Contribution of hydrophobic residues to ice binding by fish type III antifreeze protein

Type III antifreeze protein (AFP) is a 7-kDa globular protein with a flat ice-binding face centered on Ala 16. Neighboring hydrophilic residues Gln 9, Asn 14, Thr 15, Thr 18 and Gln 44 have been implicated by site-directed mutagenesis in binding to ice. These residues have the potential to form hydrogen bonds with ice, but the tight packing of side chains on the ice-binding face limits the number and strength of possible hydrogen bond interactions. Recent work with alpha-helical AFPs has emphasized the hydrophobicity of their ice-binding sites and suggests that hydrophobic interactions are important for antifreeze activity. To investigate the contribution of hydrophobic interactions between type III AFP and ice, Leu, Ile and Val residues on the rim of the ice-binding face were changed to alanine. Mutant AFPs with single alanine substitutions, L19A, V20A, and V41A, showed a 20% loss in activity. Doubly substituted mutants, L19A/V41A and L10A/I13A, had less than 50% of the activity of the wild type. Thus, side chain substitutions that leave a cavity or undercut the contact surface are almost as deleterious to antifreeze activity as those that lengthen the side chain. These mutations emphasize the importance of maintaining a specific surface contour on the ice-binding face for docking to ice.