The effect of Asp-His-Ser/Thr-Trp tetrad on the thermostability of WD40-repeat proteins

We recently found that Asp-His-Ser/Thr-Trp hydrogen-bonded tetrads are widely and uniquely present in the WD40-repeat proteins. WDR5 protein is a seven WD40-repeat propeller with five such tetrads. To explore the effect of the tetrad on the structure and stability of WD40-repeat proteins, the wild-type WDR5 and its seven mutants involving the substitutions of tetrad residues have been isolated. The crystal structures of the wild-type WDR5 and its three WDR5 mutants have been determined by X-ray diffraction method. The mutations of the tetrad residues are found not to change the basic structural features. The denaturing profiles of the wild type and the seven mutants with the use of denaturant guanidine hydrochloride have been studied by circular dichroism spectroscopy to determine the folding free energies of these proteins. The folding free energies of the wild type and the S62A, S146A, S188A, D192E, W330F, W330Y, and D324E mutants are measured to be about -11.6, -2.7, -3.1, -2.9, -3.6, -7.1, -7.0, and -7.5 kcal/mol, respectively. These suggest that (1) the hydrogen bonds in these hydrogen bond networks are unusually strong; (2) each hydrogen-bonded tetrad provides over 12 kcal/mol stability to the protein; thus, the removal of any single tetrad would cause unfolding of the protein; (3) since there are five tetrads, the protein must be in a highly unstable state without the tetrads, which might be related to its biological functions.     

Structures of single‐layer β‐sheet proteins evolved from β‐hairpin repeats

Free‐standing single‐layer β‐sheets are extremely rare in naturally occurring proteins, even though β‐sheet motifs are ubiquitous. Here we report the crystal structures of three homologous, single‐layer, anti‐parallel β‐sheet proteins, comprised of three or four twisted β‐hairpin repeats. The structures reveal that, in addition to the hydrogen bond network characteristic of β‐sheets, additional hydrophobic interactions mediated by small clusters of residues adjacent to the turns likely play a significant role in the structural stability and compensate for the lack of a compact hydrophobic core. These structures enabled identification of a family of secreted proteins that are broadly distributed in bacteria from the human gut microbiome and are putatively involved in the metabolism of complex carbohydrates. A conserved surface patch, rich in solvent‐exposed tyrosine residues, was identified on the concave surface of the β‐sheet. These new modular single‐layer β‐sheet proteins may serve as a new model system for studying folding and design of β‐rich proteins.     

C‐terminal deletion of leucine‐rich repeats from YopM reveals a heterogeneous distribution of stability in a cooperatively folded protein

Terminal deletions of units from α‐helical repeat proteins have provided insight into the physical origins of their cooperativity. To test if the same principles governing cooperativity apply to β‐sheet‐containing repeat proteins, we have created a series of C‐terminal deletion constructs from a large leucine‐rich repeat (LRR) protein, YopM. We have examined the structure and stability of the resulting deletion constructs by a combination of solution spectroscopy, equilibrium denaturation studies, and limited proteolysis. Surprisingly, a high degree of nonuniformity was found in the stability distribution of YopM. Unlike previously studied repeat proteins, we identified several key LRR that on deletion disrupt nearby structure, at distances as far away as up to three repeats, in YopM. This partial unfolding model is supported by limited proteolysis studies and by point substitution in repeats predicted to be disordered as a result of deletion of adjacent repeats. We show that key internal‐ and terminal‐caps must be present to maintain the structural integrity in adjacent regions (roughly four LRRs long) of decreased stability. The finding that full‐length YopM maintains a high level of cooperativity in equilibrium unfolding underscores the importance of interfacial interactions in stabilizing locally unstable regions of structure.     

Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Boosting protein stability with the computational design of β‐sheet surfaces

β‐sheets often have one face packed against the core of the protein and the other facing solvent. Mutational studies have indicated that the solvent‐facing residues can contribute significantly to protein stability, and that the preferred amino acid at each sequence position is dependent on the precise structure of the protein backbone and the identity of the neighboring amino acids. This suggests that the most advantageous methods for designing β‐sheet surfaces will be approaches that take into account the multiple energetic factors at play including side chain rotamer preferences, van der Waals forces, electrostatics, and desolvation effects. Here, we show that the protein design software Rosetta, which models these energetic factors, can be used to dramatically increase protein stability by optimizing interactions on the surfaces of small β‐sheet proteins. Two design variants of the β‐sandwich protein from tenascin were made with 7 and 14 mutations respectively on its β‐sheet surfaces. These changes raised the thermal midpoint for unfolding from 45°C to 64°C and 74°C. Additionally, we tested an empirical approach based on increasing the number of potential salt bridges on the surfaces of the β‐sheets. This was not a robust strategy for increasing stability, as three of the four variants tested were unfolded.     

Rotameric preferences of a protein spin label at edge‐strand β‐sheet sites

Protein spin labeling to yield the nitroxide‐based R1 side chain is a powerful method to measure protein dynamics and structure by electron spin resonance. However, R1 measurements are complicated by the flexibility of the side chain. While analysis approaches for solvent‐exposed α‐helical environment have been developed to partially account for flexibility, similar work in β‐sheets is lacking. The goal of this study is to provide the first essential steps for understanding the conformational preferences of R1 within edge β‐strands using X‐ray crystallography and double electron electron resonance (DEER) distance measurements. Crystal structures yielded seven rotamers for a non‐hydrogen‐bonded site and three rotamers for a hydrogen‐bonded site. The observed rotamers indicate contextual differences in R1 conformational preferences compared to other solvent‐exposed environments. For the DEER measurements, each strand site was paired with the same α‐helical site elsewhere on the protein. The most probable distance observed by DEER is rationalized based on the rotamers observed in the crystal structure. Additionally, the appropriateness of common molecular modeling methods that account for R1 conformational preferences are assessed for the β‐sheet environment. These results show that interpretation of R1 behavior in β‐sheets is difficult and indicate further development is needed for these computational methods to correctly relate DEER distances to protein structure at edge β‐strand sites.     

An amino acid code for β‐sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob‐socket model in an investigation of β‐sheet packing. Over a comprehensive set of β‐sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four‐residue packing cliques necessary to describe β‐sheet packing were characterized. Both occur between two adjacent hydrogen bonded β‐strands. First, defining the secondary structure packing within β‐sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β‐sheets, the knob‐socket XY:H+B consists of a three‐residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob‐sockets are found: side‐chain and main‐chain sockets. The amino acid composition of the pockets and knob‐sockets reveal the sequence specificity of β‐sheet packing. For β‐sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side‐chain and main‐chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β‐sheet structure and provide an intuitive topological mapping of β‐sheet packing. Proteins 2014; 82:2128–2140. © 2014 Wiley Periodicals, Inc.     

Sequence determinants of thermodynamic stability in a WW domain—An all‐β‐sheet protein

The stabilities of 66 sequence variants of the human Pin1 WW domain have been determined by equilibrium thermal denaturation experiments. All 34 residues composing the hPin1 WW three‐stranded β‐sheet structure could be replaced one at a time with at least one different natural or non‐natural amino acid residue without leading to an unfolded protein. Alanine substitutions at only four positions within the hPin1 WW domain lead to a partially or completely unfolded protein—in the absence of a physiological ligand. The side chains of these four residues form a conserved, partially solvent‐inaccessible, continuous hydrophobic minicore comprising the N‐ and C‐termini. Ala mutations at five other residues, three of which constitute the ligand binding patch on the concave side of the β‐sheet, significantly destabilize the hPin1 WW domain without leading to an unfolded protein. The remaining mutations affect protein stability only slightly, suggesting that only a small subset of side chain interactions within the hPin1 WW domain are mandatory for acquiring and maintaining a stable, cooperatively folded β‐sheet structure.     

TMEM59 defines a novel ATG16L1‐binding motif that promotes local activation of LC3

Selective autophagy underlies many of the important physiological roles that autophagy plays in multicellular organisms, but the mechanisms involved in cargo selection are poorly understood. Here we describe a molecular mechanism that can target conventional endosomes for autophagic degradation. We show that the human transmembrane protein TMEM59 contains a minimal 19‐amino‐acid peptide in its intracellular domain that promotes LC3 labelling and lysosomal targeting of its own endosomal compartment. Interestingly, this peptide defines a novel protein motif that mediates interaction with the WD‐repeat domain of ATG16L1, thus providing a mechanistic basis for the activity. The motif is represented with the same ATG16L1‐binding ability in other molecules, suggesting a more general relevance. We propose that this motif may play an important role in targeting specific membranous compartments for autophagic degradation, and therefore it may facilitate the search for adaptor proteins that promote selective autophagy by engaging ATG16L1. Endogenous TMEM59 interacts with ATG16L1 and mediates autophagy in response to Staphylococcus aureus infection.     

Capping motifs stabilize the leucine‐rich repeat protein PP32 and rigidify adjacent repeats

Capping motifs are found to flank most β‐strand‐containing repeat proteins. To better understand the roles of these capping motifs in organizing structure and stability, we carried out folding and solution NMR studies on the leucine‐rich repeat (LRR) domain of PP32, which is composed of five tandem LRR, capped by α‐helical and β‐hairpin motifs on the N‐ and C‐termini. We were able to purify PP32 constructs lacking either cap and containing destabilizing substitutions. Removing the C‐cap results in complete unfolding of PP32. Removing the N‐cap has a much less severe effect, decreasing stability but retaining much of its secondary structure. In contrast, the dynamics and tertiary structure of the first two repeats are significantly perturbed, based on ¹H‐¹⁵N relaxation studies, chemical shift perturbations, and residual dipolar couplings. However, more distal repeats (3 to C‐cap) retain their native tertiary structure. In this regard, the N‐cap drives the folding of adjacent repeats from what appears to be a molten‐globule‐like state. This interpretation is supported by extensive analysis using core packing substitutions in the full‐length and N‐cap‐truncated PP32. This work highlights the importance of caps to the stability and structural integrity of β‐strand‐containing LRR proteins, and emphasizes the different contributions of the N‐ and C‐terminal caps.     

Patterns of structural dynamics in RACK1 protein retained throughout evolution: A hydrogen‐deuterium exchange study of three orthologs

RACK1 is a member of the WD repeat family of proteins and is involved in multiple fundamental cellular processes. An intriguing feature of RACK1 is its ability to interact with at least 80 different protein partners. Thus, the structural features enabling such interactomic flexibility are of great interest. Several previous studies of the crystal structures of RACK1 orthologs described its detailed architecture and confirmed predictions that RACK1 adopts a seven‐bladed β‐propeller fold. However, this did not explain its ability to bind to multiple partners. We performed hydrogen‐deuterium (H‐D) exchange mass spectrometry on three orthologs of RACK1 (human, yeast, and plant) to obtain insights into the dynamic properties of RACK1 in solution. All three variants retained similar patterns of deuterium uptake, with some pronounced differences that can be attributed to RACK1's divergent biological functions. In all cases, the most rigid structural elements were confined to B‐C turns and, to some extent, strands B and C, while the remaining regions retained much flexibility. We also compared the average rate constants for H‐D exchange in different regions of RACK1 and found that amide protons in some regions exchanged at least 1000‐fold faster than in others. We conclude that its evolutionarily retained structural architecture might have allowed RACK1 to accommodate multiple molecular partners. This was exemplified by our additional analysis of yeast RACK1 dimer, which showed stabilization, as well as destabilization, of several interface regions upon dimer formation.     

Building a foundation for structure‐based cellulosome design for cellulosic ethanol: Insight into cohesin‐dockerin complexation from computer simulation

The organization and assembly of the cellulosome, an extracellular multienzyme complex produced by anaerobic bacteria, is mediated by the high‐affinity interaction of cohesin domains from scaffolding proteins with dockerins of cellulosomal enzymes. We have performed molecular dynamics simulations and free energy calculations on both the wild type (WT) and D39N mutant of the C. thermocellum Type I cohesin‐dockerin complex in aqueous solution. The D39N mutation has been experimentally demonstrated to disrupt cohesin‐dockerin binding. The present MD simulations indicate that the substitution triggers significant protein flexibility and causes a major change of the hydrogen‐bonding network in the recognition strips—the conserved loop regions previously proposed to be involved in binding—through electrostatic and salt‐bridge interactions between β‐strands 3 and 5 of the cohesin and α‐helix 3 of the dockerin. The mutation‐induced subtle disturbance in the local hydrogen‐bond network is accompanied by conformational rearrangements of the protein side chains and bound water molecules. Additional free energy perturbation calculations of the D39N mutation provide differences in the cohesin‐dockerin binding energy, thus offering a direct, quantitative comparison with experiments. The underlying molecular mechanism of cohesin‐dockerin complexation is further investigated through the free energy profile, that is, potential of mean force (PMF) calculations of WT cohesin‐dockerin complex. The PMF shows a high‐free energy barrier against the dissociation and reveals a stepwise pattern involving both the central β‐sheet interface and its adjacent solvent‐exposed loop/turn regions clustered at both ends of the β‐barrel structure.     

Cross‐strand disulfides in the hydrogen bonding site of antiparallel β‐sheet (aCSDhs): Forbidden disulfides that are highly strained,easily broken

Some disulfide bonds perform important structural roles in proteins, but another group has functional roles via redox reactions. Forbidden disulfides are stressed disulfides found in recognizable protein contexts, which currently constitute more than 10% of all disulfides in the PDB. They likely have functional redox roles and constitute a major subset of all redox‐active disulfides. The torsional strain of forbidden disulfides is typically higher than for structural disulfides, but not so high as to render them immediately susceptible to reduction under physionormal conditions. Previously we characterized the most abundant forbidden disulfide in the Protein Data Bank, the aCSDn: a canonical motif in which disulfide‐bonded cysteine residues are positioned directly opposite each other on adjacent anti‐parallel β‐strands such that the backbone hydrogen‐bonded moieties are directed away from each other. Here we perform a similar analysis for the aCSDh, a less common motif in which the opposed cysteine residues are backbone hydrogen bonded. Oxidation of two Cys in this context places significant strain on the protein system, with the β‐chains tilting toward each other to allow disulfide formation. Only left‐handed aCSDh conformations are compatible with the inherent right‐handed twist of β‐sheets. aCSDhs tend to be more highly strained than aCSDns, particularly when both hydrogen bonds are formed. We discuss characterized roles of aCSDh motifs in proteins of the dataset, which include catalytic disulfides in ribonucleotide reductase and ahpC peroxidase as well as a redox‐active disulfide in P1 lysozyme, involved in a major conformation change. The dataset also includes many binding proteins.     

Geometrical principles of homomeric β‐barrels and β‐helices: Application to modeling amyloid protofilaments

Examples of homomeric β‐helices and β‐barrels have recently emerged. Here we generalize the theory for the shear number in β‐barrels to encompass β‐helices and homomeric structures. We introduce the concept of the “β‐strip,” the set of parallel or antiparallel neighboring strands, from which the whole helix can be generated giving it n‐fold rotational symmetry. In this context, the shear number is interpreted as the sum around the helix of the fixed register shift between neighboring identical β‐strips. Using this approach, we have derived relationships between helical width, pitch, angle between strand direction and helical axis, mass per length, register shift, and number of strands. The validity and unifying power of the method is demonstrated with known structures including α‐hemolysin, T4 phage spike, cylindrin, and the HET‐s(218‐289) prion. From reported dimensions measured by X‐ray fiber diffraction on amyloid fibrils, the relationships can be used to predict the register shift and the number of strands within amyloid protofilaments. This was used to construct models of transthyretin and Alzheimer β(40) amyloid protofilaments that comprise a single strip of in‐register β‐strands folded into a “β‐strip helix.” Results suggest both stabilization of an individual β‐strip helix and growth by addition of further β‐strip helices can involve the same pair of sequence segments associating with β‐sheet hydrogen bonding at the same register shift. This process would be aided by a repeat sequence. Hence, understanding how the register shift (as the distance between repeat sequences) relates to helical dimensions will be useful for nanotube design.     

Identifying the Hotspots on the Top Faces of WD40-Repeat Proteins from Their Primary Sequences by β-Bulges and DHSW Tetrads

The analysis of 36 available crystal structures of WD40 repeat proteins reveals widespread existence of a beta-bulge formed at the beginning of strand a and the end of strand b, termed as WD(b-a) bulge: among a total of 259 WD40 blades, there are 243 such β-bulges. The R(1) positions in these WD(b-a) bulges have fair distributions of Arg, His, Ile, Leu, Lys, Met, Phe, Trp, Tyr and Val residues. These residues protrude on the top face of the WD40 proteins and can serve as hotspots for protein-protein interactions. An analysis of 29 protein complexes formed by 17 WD proteins reveals that these R(1) residues, along with two other residues (R(1)-2 and D-1), are indeed widely involved in protein-protein interactions. Interestingly, these WD(b-a) bulges can be easily identified by the 4-amino acid sequences of (V, L, I), R(1), R(2), (V, L, I), along with some other significant amino acids. Thus, the hotspots of WD40 proteins on the top face can be readily predicted based on the primary sequences of the proteins. The literature-reported mutagenesis studies for Met30, MDV1, Tup11, COP1 and SPA1, which crystal structures are not available, can be readily understood based on the feature-based method. Applying the method, the twelve potential hotspots on the top face of Tup11 from S. japonicas have been identified. Our ITC measurements confirm seven of them, Tyr382, Arg284, Tyr426, Tyr508, Leu559, Lys575 and Ile601, are essential for recognizing Fep1. The ITC measurements further convinced that the feature-based method provides accurate prediction of hotspots on the top face.     

Mapping of O‐linked β‐N‐acetylglucosamine modification sites in key contractile proteins of rat skeletal muscle

O‐linked β‐N‐acetylglucosamine (O‐GlcNAc) is a widespread modification of serine/threonine residues of nucleocytoplasmic proteins. Recently, several key contractile proteins in rat skeletal muscle (i.e., myosin heavy and light chains and actin) were identified as O‐GlcNAc modified. Moreover, it was demonstrated that O‐GlcNAc moieties involved in contractile protein interactions could modulate Ca²⁺ activation parameters of contraction. In order to better understand how O‐GlcNAc can modulate the contractile activity of muscle fibers, we decided to identify the sites of O‐GlcNAc modification in purified contractile protein homogenates. Using an MS‐based method that relies on mild β‐elimination followed by Michael addition of DTT (BEMAD), we determined the localization of one O‐GlcNAc site in the subdomain four of actin and four O‐GlcNAc sites in the light meromyosin region of myosin heavy chains (MHC). According to previous reports concerning the role of these regions, our data suggest that O‐GlcNAc sites might modulate the actin–tropomyosin interaction, and be involved in MHC polymerization or interactions between MHC and other contractile proteins. Thus, the results suggest that this PTM might be involved in protein–protein interactions but could also modulate the contractile properties of skeletal muscle.     

The amino terminal domain of a novel WD repeat protein from Trypanosoma cruzi contains a non-canonical mitochondrial targeting signal

WD (tryptophan/aspartic acid) repeat proteins perform a wide variety of functions in eukaryotic cells. They are characterised by the presence of a number of conserved repeat motifs that contribute to the beta-propeller structures which are the common feature of this large group of proteins. We report here the properties of the first characterised member of this family in the American trypanosome, Trypanosoma cruzi (TcBPP1). In the CL Brener clone the protein is 482 amino acids long and is predicted to contain four WD repeat motifs, flanked by amino and carboxyl terminal extensions. TcBPP1 is a single copy gene present on a 1.0/1.6 Mb pair of homologous chromosomes in a locus that is syntenic with the corresponding regions of Trypanosoma brucei and Leishmania major chromosomes. Consistent with the proposed hybrid nature of the CL Brener clone, the proteins encoded by the two different alleles share only 97% identity at the amino acid level. To determine subcellular location, we examined transfected parasites for the distribution of green fluorescent protein (GFP) fused with different regions of TcBPP1. These studies demonstrated that a 115 amino acid peptide derived from the amino terminal domain of TcBPP1 is able to target GFP to the mitochondrion. Interestingly this region lacks a typical amino terminal presequence suggesting that mitochondrial import is mediated by an alternative targeting signal.     

Optimization of a β‐sheet‐cap for long loop closure

Protein loops make up a large portion of the secondary structure in nature. But very little is known concerning loop closure dynamics and the effects of loop composition on fold stability. We have designed a small system with stable β‐sheet structures, including features that allow us to probe these questions. Using paired Trp residues that form aromatic clusters on folding, we are able to stabilize two β‐strands connected by varying loop lengths and composition (an example sequence: R W ITVTI – loop – KKIRV W E). Using NMR and CD, both fold stability and folding dynamics can be investigated for these systems. With the 16 residue loop peptide (sequence: R W ITVTI‐(GGGGKK)₂GGGG‐KKIRV W E) remaining folded (ΔG_U = 1.6 kJ/mol at 295K). To increase stability and extend the series to longer loops, we added an additional Trp/Trp pair in the loop flanking position. With this addition to the strands, the 16 residue loop (sequence: R W ITVRI W ‐(GGGGKK)₂GGGG‐ W KTIRV W E) supports a remarkably stable β‐sheet (ΔG_U = 6.3 kJ/mol at 295 K, T_m = ～55°C). Given the abundance of loops in binding motifs and between secondary structures, these constructs can be powerful tools for peptide chemists to study loop effects; with the Trp/Trp pair providing spectroscopic probes for assessing both stability and dynamics by NMR.