A thermodynamic and kinetic analysis of the folding pathway of an SH3 domain entropically stabilised by a redesigned hydrophobic core

The folding thermodynamics and kinetics of the alpha-spectrin SH3 domain with a redesigned hydrophobic core have been studied. The introduction of five replacements, A11V, V23L, M25V, V44I and V58L, resulted in an increase of 16% in the overall volume of the side-chains forming the hydrophobic core but caused no remarkable changes to the positions of the backbone atoms. Judging by the scanning calorimetry data, the increased stability of the folded structure of the new SH3-variant is caused by entropic factors, since the changes in heat capacity and enthalpy upon the unfolding of the wild-type and mutant proteins were identical at 298 K. It appears that the design process resulted in an increase in burying both the hydrophobic and hydrophilic surfaces, which resulted in a compensatory effect upon the changes in heat capacity and enthalpy. Kinetic analysis shows that both the folding and unfolding rate constants are higher for the new variant, suggesting that its transition state becomes more stable compared to the folded and unfolded states. The phi(double dagger-U) values found for a number of side-chains are slightly lower than those of the wild-type protein, indicating that although the transition state ensemble (TSE) did not change overall, it has moved towards a more denatured conformation, in accordance with Hammond's postulate. Thus, the acceleration of the folding-unfolding reactions is caused mainly by an improvement in the specific and/or non-specific hydrophobic interactions within the TSE rather than by changes in the contact order. Experimental evidence showing that the TSE changes globally according to its hydrophobic content suggests that hydrophobicity may modulate the kinetic behaviour and also the folding pathway of a protein.     

Analysis of anisotropic side-chain packing in proteins and application to high-resolution structure prediction

pi-pi, Cation-pi, and hydrophobic packing interactions contribute specificity to protein folding and stability to the native state. As a step towards developing improved models of these interactions in proteins, we compare the side-chain packing arrangements in native proteins to those found in compact decoys produced by the Rosetta de novo structure prediction method. We find enrichments in the native distributions for T-shaped and parallel offset arrangements of aromatic residue pairs, in parallel stacked arrangements of cation-aromatic pairs, in parallel stacked pairs involving proline residues, and in parallel offset arrangements for aliphatic residue pairs. We then investigate the extent to which the distinctive features of native packing can be explained using Lennard-Jones and electrostatics models. Finally, we derive orientation-dependent pi-pi, cation-pi and hydrophobic interaction potentials based on the differences between the native and compact decoy distributions and investigate their efficacy for high-resolution protein structure prediction. Surprisingly, the orientation-dependent potential derived from the packing arrangements of aliphatic side-chain pairs distinguishes the native structure from compact decoys better than the orientation-dependent potentials describing pi-pi and cation-pi interactions.     

The relationship between conservation, thermodynamic stability, and function in the SH3 domain hydrophobic core

To investigate the relationships between sequence conservation, protein stability, and protein function, we have measured the thermodynamic stability, folding kinetics, and in vitro peptide-binding activity of a large number of single-site substitutions in the hydrophobic core of the Fyn SH3 domain. Comparison of these data to that derived from an analysis of a large alignment of SH3 domain sequences revealed a very good correlation between the distinct pattern of conservation observed at each core position and the thermodynamic stability of mutants. Conservation was also found to correlate well with the unfolding rates of mutants, but not to the folding rates, suggesting that evolution selects more strongly for optimal native state packing interactions than for maximal folding rates. Structural analysis suggests that residue-residue core packing interactions are very similar in all SH3 domains, which provides an explanation for the correlation between conservation and mutant stability effects studied in a single SH3 domain. We also demonstrate a correlation between stability and the in vivo activity of mutants, and between conservation and activity. However, the relationship between conservation and activity was very strong only for the three most conserved hydrophobic core positions. The weaker correlation between activity and conservation seen at the other seven core positions indicates that maintenance of protein stability is the dominant selective pressure at these positions. In general, the pattern of conservation at hydrophobic core positions appears to arise from conserved packing constraints, and can be effectively utilized to predict the destabilizing effects of amino acid substitutions.     

Appendant structure plays an important role in amyloidogenic cystatin dimerization prior to domain swapping

It has been hypothesized that prior to protein domain swapping, unfolding occurs in regions important for the stability of the native monomeric structure, which probably increases the possibility of intermolecular interaction. In order to explore the detailed information of the important unfolding regions in cystatin prior to domain swapping, 20?ns molecular dynamic simulations were performed at atomic level with typical amyloidogenic chicken cystatin (cC) mutant I66Q monomer under conditions that enable forming amyloid fibrils in biological experiments. Our results showed that I66Q mutant exhibited relatively large secondary structure changes and obvious expanding tendency of hydrophobic core compared to wild-type cC. More importantly, the appendant structure (AS) showed a large displacement and distortion towards the hydrophobic core in amyloidogenic cystatin. The structural analysis on cystatin monomer suggested that structural changes of the AS might make the hydrophobic core expand more easily. In addition, analysis on docking dimer has shown that the distorted AS was favor to intermolecular interactions between two cystatin monomers. Data from an independent theoretical derived algorithm as well as biological experiments also support this hypothesis.     

Crystal structure of a cAMP-dependent protein kinase mutant at 1.26A: new insights into the catalytic mechanism

The catalytic subunit of cAMP-dependent protein kinase has served as a paradigm for the entire kinase family. In the course of studying the structure-function relationship of the P+1 loop (Leu198-Leu205) of the kinase, we have solved the crystal structure of the Tyr204 to Ala mutant in complexes with Mg.ATP and an inhibitory peptide at 1.26A, with overall structure very similar to that of the wild-type protein. However, at the nucleotide binding site, ATP was found largely hydrolyzed, with the products ADP-PO(4) retained in the structure. High-resolution refinement suggests that 26% of the molecules contain the intact ATP, whereas 74% have the hydrolyzed products. The observation of the substrate and product states in the same structure adds significant information to our understanding of the phosphoryl transfer process. Structural examination of the mutation site substantiates and extends the emerging concept that the hydrophobic core in the large lobe of the kinase might serve as a stable platform for anchoring key segments involved in catalysis. We propose that Tyr204 is critical for anchoring the P+1 loop to the core. Further analysis has highlighted two major connections between the P+1 loop and the catalytic loop (Arg165-Asn171). One emphasizes the hydrophobic packing of Tyr204 and Leu167 mediated through residues from the alphaF-helix, recently recognized as a signal integration motif, which together with the alphaE-helix forms the center of the hydrophobic core network. The other connection is mediated by the hydrogen bond interaction between Thr201 and Asp166, in a substrate-dependent manner. We speculate that the latter interaction may be important for the kinase to sense the presence of substrate and prepare itself for the catalytic reaction. Thus, the P+1 loop is not merely involved in substrate binding; it mediates the communication between substrate and catalytic residues.     

Solution structure of switch Arc,a mutant with 3(10) helices replacing a wild-type beta-ribbon

Adjacent N11L and L12N mutations in the antiparallel beta-ribbon of Arc repressor result in dramatic changes in local structure in which each beta-strand is replaced by a right-handed helix. The full solution structure of this "switch" Arc mutant shows that irregular 3(10) helices compose the new secondary structure. This structural metamorphosis conserves the number of main-chain and side-chain to main-chain hydrogen bonds and the number of fully buried core residues. Apart from a slight widening of the interhelical angle between alpha-helices A and B and changes in side-chain conformation of a few core residues in Arc, no large-scale structural adjustments in the remainder of the protein are necessary to accommodate the ribbon-to-helix change. Nevertheless, some changes in hydrogen-exchange rates are observed, even in regions that have very similar structures in the two proteins. The surface of switch Arc is packed poorly compared to wild-type, leading to approximately 1000A(2) of additional solvent-accessible surface area, and the N termini of the 3(10) helices make unfavorable head-to-head electrostatic interactions. These structural features account for the positive m value and salt dependence of the ribbon-to-helix transition in Arc-N11L, a variant that can adopt either the mutant or wild-type structures. The tertiary fold is capped in different ways in switch and wild-type Arc, showing how stepwise evolutionary transformations can arise through small changes in amino acid sequence.     

The dimerization mechanism of LIS1 and its implication for proteins containing the LisH motif

Miller-Dieker lissencephaly, or "smooth-brain" is a debilitating genetic developmental syndrome of the cerebral cortex, and is linked to mutations in the Lis1 gene. The LIS1 protein contains a so-called LisH motif at the N terminus, followed by a coiled-coil region and a seven WD-40 repeat forming beta-propeller structure. In vivo and in vitro, LIS1 is a dimer, and the dimerization is mediated by the N-terminal fragment and is essential for the protein's biological function. The recently determined crystal structure of the murine LIS1 N-terminal fragment encompassing residues 1-86 (N-LIS1) revealed that the LisH motif forms a tightly associated homodimer with a four-helix antiparallel bundle core, while the parallel coiled-coil situated downstream is stabilized by three canonical heptad repeats. This homodimer is uniquely asymmetric because of a distinct kink in one of the helices. Because the LisH motif is widespread among many proteins, some of which are implicated in human diseases, we investigated in detail the mechanism of N-LIS1 dimerization. We found that dimerization is dependent on both the LisH motif and the residues downstream of it, including the first few turns of the helix. We also have found that the coiled-coil does not contribute to dimerization, but instead is very labile and can adopt both supercoiled and helical conformations. These observations suggest that the presence of the LisH motif alone is not sufficient for high-affinity homodimerization and that other structural elements are likely to play an important role in this large family of proteins. The observed lability of the coiled-coil fragment in LIS1 is most likely of functional importance.     

The structure of the N-terminal domain of riboflavin synthase in complex with riboflavin at 2.6A resolution

Riboflavin synthase of Escherichia coli is a homotrimer with a molecular mass of 70 kDa. The enzyme catalyzes the dismutation of 6,7-dimethyl-8-(1'-D-ribityl)-lumazine, affording riboflavin and 5-amino-6-ribitylamino-2,4(1H,3H)-pyrimidinedione. The N-terminal segment (residues 1-87) and the C-terminal segment (residues 98-187) form beta-barrels with similar fold and a high degree of sequence similarity. A recombinant peptide comprising amino acid residues 1-97 forms a dimer, which binds riboflavin with high affinity. Here, we report the structure of this construct in complex with riboflavin at 2.6A resolution. It is demonstrated that the complex can serve as a model for ligand-binding in the native enzyme. The structure and riboflavin-binding mode is in excellent agreement with structural information obtained from the native enzyme from Escherichia coli and riboflavin synthase from Schizosaccharomyces pombe. The implications for the binding specificity and the regiospecificity of the catalyzed reaction are discussed.     

The jigsaw puzzle model: search for conformational specificity in protein interiors

The jigsaw puzzle model postulates that the predominant factor relating primary sequence to three-dimensional fold lies in the stereospecific packing of interdigitating side-chains within densely packed protein interiors. An attempt has been made to check the validity of the model by means of a surface complementarity function. Out of a database of 100 highly resolved protein structures the contacts between buried hydrophobic residues (Leu, Ile, Val, Phe) and their neighbours have been categorized in terms of the extent of side-chain surface area involved in a contact (overlap) and their steric fit (Sm). The results show that the majority of contacts between a buried residue and its immediate neighbours (side-chains) are of high steric fit and in the case of extended overlap at least one of the angular parameters characterizing interresidue geometry to have pronounced deviation from a random distribution, estimated by chi(2). The calculations thus tend to support the "jigsaw puzzle" model in that 75-85% of the contacts involving hydrophobic residues are of high surface complementarity, which, coupled to high overlap, exercise fairly stringent constraints over the possible geometrical orientations between interacting residues. These constraints manifest in simple patterns in the distributions of orientational angles. Approximately 60-80% of the buried side-chain surface packs against neighbouring side-chains, the rest interacting with main-chain atoms. The latter partition of the surface maintains an equally high steric fit (relative to side-chain contacts) emphasizing a non-trivial though secondary role played by main-chain atoms in interior packing. The majority of this class of contacts, though of high complementarity, is of reduced overlap. All residues whether hydrophobic or polar/charged show similar surface complementarity measures upon burial, indicating comparable competence of all amino acids in packing effectively with their atomic environments. The specificity thus appears to be distributed over the entire network of contacts within proteins. The study concludes with a proposal to classify contacts as specific and non-specific (based on overlap and fit), with the former perhaps contributing more to the specificity between sequence and fold than the latter.