Beyond consensus: statistical free energies reveal hidden interactions in the design of a TPR motif

Statistical Analysis and Tokenization of Epitopes to Construct Artificial Neoepitope Libraries

Epitopes are specific regions on an antigen's surface that the immune system recognizes. Epitopes are usually protein regions on foreign immune-stimulating entities such as viruses and bacteria, and in some cases, endogenous proteins may act as antigens. Identifying epitopes is crucial for accelerating the development of vaccines and immunotherapies. However, mapping epitopes in pathogen proteomes is challenging using conventional methods. Screening artificial neoepitope libraries against antibodies can overcome this issue. Here, we applied conventional sequence analysis and methods inspired in natural language processing to reveal specific sequence patterns in the linear epitopes deposited in the Immune Epitope Database (www.iedb.org) that can serve as building blocks for the design of universal epitope libraries. Our results reveal that amino acid frequency in annotated linear epitopes differs from that in the human proteome. Aromatic residues are overrepresented, while the presence of cysteines is practically null in epitopes. Byte pair encoding tokenization shows high frequencies of tryptophan in tokens of 5, 6, and 7 amino acids, corroborating the findings of the conventional sequence analysis. These results can be applied to reduce the diversity of linear epitope libraries by orders of magnitude.

Engineered Protein-Driven Synthesis of Tunable Iron Oxide Nanoparticles as T1 and T2 Magnetic Resonance Imaging Contrast Agents

Iron oxide nanoparticles (IONPs) have become one of the most promising nanomaterials for biomedical applications because of their biocompatibility and physicochemical properties. This study demonstrates the use of protein engineering as a novel approach to design scaffolds for the tunable synthesis of ultrasmall IONPs. Rationally designed proteins, containing different number of metal-coordination sites, were evaluated to control the size and the physicochemical and magnetic properties of a set of protein-stabilized IONPs (Prot-IONPs). Prot-IONPs, synthesized through an optimized coprecipitation approach, presented good T1 and T2 relaxivity values, stability, and biocompatibility, showing potential for magnetic resonance imaging (MRI) applications.

Unraveling the Mechanics of a Repeat-Protein Nanospring: From Folding of Individual Repeats to Fluctuations of the Superhelix

Tandem-repeat proteins comprise small secondary structure motifs that stack to form one-dimensional arrays with distinctive mechanical properties that are proposed to direct their cellular functions. Here, we use single-molecule optical tweezers to study the folding of consensus-designed tetratricopeptide repeats (CTPRs), superhelical arrays of short helix-turn-helix motifs. We find that CTPRs display a spring-like mechanical response in which individual repeats undergo rapid equilibrium fluctuations between partially folded and unfolded conformations. We rationalize the force response using Ising models and dissect the folding pathway of CTPRs under mechanical load, revealing how the repeat arrays form from the center toward both termini simultaneously. Most strikingly, we also directly observe the protein's superhelical tertiary structure in the force signal. Using protein engineering, crystallography, and single-molecule experiments, we show that the superhelical geometry can be altered by carefully placed amino acid substitutions, and we examine how these sequence changes affect intrinsic repeat stability and inter-repeat coupling. Our findings provide the means to dissect and modulate repeat-protein stability and dynamics, which will be essential for researchers to understand the function of natural repeat proteins and to exploit artificial repeats proteins in nanotechnology and biomedical applications.

Unraveling the Mechanics of a Repeat-Protein Nanospring: From Folding of Individual Repeats to Fluctuations of the Superhelix

Tandem-repeat proteins comprise small secondary structure motifs that stack to form one-dimensional arrays with distinctive mechanical properties that are proposed to direct their cellular functions. Here, we use single-molecule optical tweezers to study the folding of consensus-designed tetratricopeptide repeats (CTPRs), superhelical arrays of short helix-turn-helix motifs. We find that CTPRs display a spring-like mechanical response in which individual repeats undergo rapid equilibrium fluctuations between partially folded and unfolded conformations. We rationalize the force response using Ising models and dissect the folding pathway of CTPRs under mechanical load, revealing how the repeat arrays form from the center toward both termini simultaneously. Most strikingly, we also directly observe the protein's superhelical tertiary structure in the force signal. Using protein engineering, crystallography, and single-molecule experiments, we show that the superhelical geometry can be altered by carefully placed amino acid substitutions, and we examine how these sequence changes affect intrinsic repeat stability and inter-repeat coupling. Our findings provide the means to dissect and modulate repeat-protein stability and dynamics, which will be essential for researchers to understand the function of natural repeat proteins and to exploit artificial repeats proteins in nanotechnology and biomedical applications.

Antigenic evolution of human influenza H3N2 neuraminidase is constrained by charge balancing

As one of the main influenza antigens, neuraminidase (NA) in H3N2 virus has evolved extensively for more than 50 years due to continuous immune pressure. While NA has recently emerged as an effective vaccine target, biophysical constraints on the antigenic evolution of NA remain largely elusive. Here, we apply combinatorial mutagenesis and next-generation sequencing to characterize the local fitness landscape in an antigenic region of NA in six different human H3N2 strains that were isolated around 10 years apart. The local fitness landscape correlates well among strains and the pairwise epistasis is highly conserved. Our analysis further demonstrates that local net charge governs the pairwise epistasis in this antigenic region. In addition, we show that residue coevolution in this antigenic region is correlated with the pairwise epistasis between charge states. Overall, this study demonstrates the importance of quantifying epistasis and the underlying biophysical constraint for building a model of influenza evolution.

Partial Consensus Design and Enhancement of Protein Function by Secondary-Structure-Guided Consensus Mutations

Consensus design (CD) is a representative sequence-based protein design method that enables the design of highly functional proteins by analyzing vast amounts of protein sequence data. This study proposes a partial consensus design (PCD) of a protein as a derivative approach of CD. The method replaces the target protein sequence with a consensus sequence in a secondary-structure-dependent manner (i.e., regionally dependent and divided into α-helix, β-sheet, and loop regions). In this study, we generated several artificial partial consensus l-threonine 3-dehydrogenases (PcTDHs) by PCD using the TDH from Cupriavidus necator (CnTDH) as a target protein. Structural and functional analysis of PcTDHs suggested that thermostability would be independently improved when consensus mutations are introduced into the loop region of TDHs. On the other hand, enzyme kinetic parameters (k_cat/K_m) and average productivity would be synergistically enhanced by changing the combination of the mutations-replacement of one region of CnTDH with a consensus sequence provided only negative effects, but the negative effects were nullified when the two regions were replaced simultaneously. Taken together, we propose the hypothesis that there are protein regions that encode individual protein properties, such as thermostability and activity, and that the introduction of consensus mutations into these regions could additively or synergistically modify their functions.

The Effect of Mutations in the TPR and Ankyrin Families of Alpha Solenoid Repeat Proteins

Protein repeats are short, highly similar peptide motifs that occur several times within a single protein, for example the TPR and Ankyrin repeats. Understanding the role of mutation in these proteins is complicated by the competing facts that 1) the repeats are much more restricted to a set sequence than non-repeat proteins, so mutations should be harmful much more often because there are more residues that are heavily restricted due to the need of the sequence to repeat and 2) the symmetry of the repeats in allows the distribution of functional contributions over a number of residues so that sometimes no specific site is singularly responsible for function (unlike enzymatic active site catalytic residues). To address this issue, we review the effects of mutations in a number of natural repeat proteins from the tetratricopeptide and Ankyrin repeat families. We find that mutations are context dependent. Some mutations are indeed highly disruptive to the function of the protein repeats while mutations in identical positions in other repeats in the same protein have little to no effect on structure or function.

A cluster of Ankyrin and Ankyrin-TPR repeat genes is associated with panicle branching diversity in rice

The number of grains per panicle is an important yield-related trait in cereals which depends in part on panicle branching complexity. One component of this complexity is the number of secondary branches per panicle. Previously, a GWAS site associated with secondary branch and spikelet numbers per panicle in rice was identified. Here we combined gene capture, bi-parental genetic population analysis, expression profiling and transgenic approaches in order to investigate the functional significance of a cluster of 6 ANK and ANK-TPR genes within the QTL. Four of the ANK and ANK-TPR genes present a differential expression associated with panicle secondary branch number in contrasted accessions. These differential expression patterns correlate in the different alleles of these genes with specific deletions of potential cis-regulatory sequences in their promoters. Two of these genes were confirmed through functional analysis as playing a role in the control of panicle architecture. Our findings indicate that secondary branching diversity in the rice panicle is governed in part by differentially expressed genes within this cluster encoding ANK and ANK-TPR domain proteins that may act as positive or negative regulators of panicle meristem’s identity transition from indeterminate to determinate state.

Grain yield is one of the most important indexes in rice breeding, which is controlled in part by panicle branching complexity. A new QTL with co-location of spikelet number (SpN) and secondary branch number (SBN) traits was identified by genome-wide association study in a Vietnamese rice landrace panel. A set of four Ankyrin and Tetratricopeptide repeat domain-encoding genes was identified from this QTL based on their difference of expression levels between two contrasted haplotypes for the SpN and SBN traits. The differential expression is correlated with deletions in the promoter regions of these genes. Two of the genes act as negative regulators of the panicle meristem’s identity transition from indeterminate to determinate state while the other two act as positive regulators of this meristem fate transition. Based on the different phenotypes between overexpressed and mutant plants, two of these genes were confirmed as playing a role in the control of panicle architecture. These findings can be directly used to assist selection for grain yield improvement.

An Analytical Review of the Structural Features of Pentatricopeptide Repeats: Strategic Amino Acids, Repeat Arrangements and Superhelical Architecture

Tricopeptide repeats are common in natural proteins, and are exemplified by 34- and 35-residue repeats, known respectively as tetratricopeptide repeats (TPRs) and pentatricopeptide repeats (PPRs). In both classes, each repeat unit forms an antiparallel bihelical structure, so that multiple such units in a polypeptide are arranged in a parallel fashion. The primary structures of the motifs are nonidentical, but amino acids of similar properties occur in strategic positions. The focus of the present work was on PPR, but TPR, its better-studied cousin, is often included for comparison. The analyses revealed that critical amino acids, namely Gly, Pro, Ala and Trp, were placed at distinct locations in the higher order structure of PPR domains. While most TPRs occur in repeats of three, the PPRs exhibited a much greater diversity in repeat numbers, from 1 to 30 or more, separated by spacers of various sequences and lengths. Studies of PPR strings in proteins showed that the majority of PPR units are single, and that the longer tandems (i.e., without space in between) occurred in decreasing order. The multi-PPR domains also formed superhelical vortices, likely governed by interhelical angles rather than the spacers. These findings should be useful in designing and understanding the PPR domains.

Improved thermostability of creatinase from Alcaligenes Faecalis through non-biased phylogenetic consensus-guided mutagenesis

Background

Enzymatic quantification of creatinine has become an essential method for clinical evaluation of renal function. Although creatinase (CR) is frequently used for this purpose, its poor thermostability severely limits industrial applications. Herein, we report a novel creatinase from Alcaligenes faecalis (afCR) with higher catalytic activity and lower K_M value, than currently used creatinases. Furthermore, we developed a non-biased phylogenetic consensus method to improve the thermostability of afCR.

Results

We applied a non-biased phylogenetic consensus method to identify 59 candidate consensus residues from 24 creatinase family homologs for screening afCR mutants with improved thermostability. Twenty-one amino acids of afCR were selected to mutagenesis and 11 of them exhibited improved thermostability compared to the parent enzyme (afCR-M0). Combination of single-site mutations in sequential screens resulted in a quadruple mutant D17V/T199S/L6P/T251C (M4-2) which showed ~ 1700-fold enhanced half-life at 57 °C and a 4.2 °C higher T₅₀¹⁵ than that of afCR-M0. The mutant retained catalytic activity equivalent to afCR-M0, and thus showed strong promise for application in creatinine detection. Structural homology modeling revealed a wide range of potential molecular interactions associated with individual mutations that contributed to improving afCR thermostability.

Conclusions

Results of this study clearly demonstrated that the non-biased-phylogenetic consensus design for improvement of thermostability in afCR is effective and promising in improving the thermostability of more enzymes.

Protein Sequence Selection Method That Enables Full Consensus Design of Artificial l-Threonine 3-Dehydrogenases with Unique Enzymatic Properties

Exponentially increasing protein sequence data enables artificial enzyme design using sequence-based protein design methods, including full-consensus protein design (FCD). The success of artificial enzyme design is strongly dependent on the nature of the sequences used. Hence, sequences must be selected from databases and curated libraries prepared to enable a successful design by FCD. In this study, we proposed a selection approach regarding several key residues as sequence motifs. We used l-threonine 3-dehydrogenase (TDH) as a model to test the validity of this approach. In the classification, four residues (143, 174, 188, and 214) were used as key residues. We classified thousands of TDH homologous sequences into five groups containing hundreds of sequences. Utilizing sequences in the libraries, we designed five artificial TDHs by FCD. Among the five, we successfully expressed four in soluble form. Biochemical analysis of artificial TDHs indicated that their enzymatic properties vary; half of the maximum measured enzyme activity (t_1/2) and activation energies were distributed from 53 to 65 °C and from 38 to 125 kJ/mol, respectively. The artificial TDHs had unique kinetic parameters, distinct from one another. Structural analysis indicates that consensus mutations are mainly introduced in the secondary or outer shell. The functional diversity of the artificial TDHs is due to the accumulation of mutations that affect their physicochemical properties. Taken together, our findings indicate that our proposed approach can help generate artificial enzymes with unique enzymatic properties.

Consensus protein engineering on the thermostable histone-like bacterial protein HUs significantly improves stability and DNA binding affinity

Consensus-based protein engineering strategy has been applied to various proteins and it can lead to the design of proteins with enhanced biological performance. Histone-like HUs comprise a protein family with sequence variety within a highly conserved 3D-fold. HU function includes compacting and regulating bacterial DNA in a wide range of biological conditions in bacteria. To explore the possible impact of consensus-based design in the thermodynamic stability of HU proteins, the approach was applied using a dataset of sequences derived from a group of 40 mesostable, thermostable, and hyperthermostable HUs. The consensus-derived HU protein was named HUBest, since it is expected to perform best. The synthetic HU gene was overexpressed in E. coli and the recombinant protein was purified. Subsequently, HUBest was characterized concerning its correct folding and thermodynamic stability, as well as its ability to interact with plasmid DNA. A substantial increase in HUBest stability at high temperatures is observed. HUBest has significantly improved biological performance at ambience temperature, presenting very low K_d values for binding plasmid DNA as indicated from the Gibbs energy profile of HUBest. This K_d may be associated to conformational changes leading to decreased thermodynamic stability and, therefore, higher flexibility at ambient temperature.

Competition of individual domain folding with inter-domain interaction in WW domain engineered repeat proteins

Engineered repeat proteins have proven to be a fertile ground for studying the competition between folding, misfolding and transient aggregation of tethered protein domains. We examine the interplay between folding and inter-domain interactions of engineered FiP35 WW domain repeat proteins with n = 1 through 5 repeats. We characterize protein expression, thermal and guanidium melts, as well as laser T-jump kinetics. All experimental data is fitted by a global fitting model with two states per domain (U, N), plus a third state M to account for non-native states due to domain interactions present in all but the monomer. A detailed structural model is provided by coarse-grained simulated annealing using the AWSEM Hamiltonian. Tethered FiP35 WW domains with n = 2 and 3 domains are just slightly less stable than the monomer. The n = 4 oligomer is yet less stable, its expression yield is much lower than the monomer's, and depends on the purification tag used. The n = 5 plasmid did not express at all, indicating the sudden onset of aggregation past n = 4. Thus, tethered FiP35 has a critical nucleus size for inter-domain aggregation of n ≈ 4. According to our simulations, misfolded structures become increasingly prevalent as one proceeds from monomer to pentamer, with extended inter-domain beta sheets appearing first, then multi-sheet 'intramolecular amyloid' structures, and finally novel motifs containing alpha helices. We discuss the implications of our results for oligomeric aggregate formation and structure, transient aggregation of proteins whilst folding, as well as for protein evolution that starts with repeat proteins.

Phylogenetic spread of sequence data affects fitness of consensus enzymes: Insights from triosephosphate isomerase

The concept of consensus in multiple sequence alignments (MSAs) has been used to design and engineer proteins previously with some success. However, consensus design implicitly assumes that all amino acid positions function independently, whereas in reality, the amino acids in a protein interact with each other and work cooperatively to produce the optimum structure required for its function. Correlation analysis is a tool that can capture the effect of such interactions. In a previously published study, we made consensus variants of the triosephosphate isomerase (TIM) protein using MSAs that included sequences form both prokaryotic and eukaryotic organisms. These variants were not completely native-like and were also surprisingly different from each other in terms of oligomeric state, structural dynamics, and activity. Extensive correlation analysis of the TIM database has revealed some clues about factors leading to the unusual behavior of the previously constructed consensus proteins. Among other things, we have found that the more ill-behaved consensus mutant had more broken correlations than the better-behaved consensus variant. Moreover, we report three correlation and phylogeny-based consensus variants of TIM. These variants were more native-like than the previous consensus mutants and considerably more stable than a wild-type TIM from a mesophilic organism. This study highlights the importance of choosing the appropriate diversity of MSA for consensus analysis and provides information that can be used to engineer stable enzymes.

Benchmark Analysis of Native and Artificial NAD+-Dependent Enzymes Generated by a Sequence-Based Design Method with or without Phylogenetic Data

The expansion of protein sequence databases has enabled us to design artificial proteins by sequence-based design methods, such as full-consensus design (FCD) and ancestral-sequence reconstruction (ASR). Artificial proteins with enhanced activity levels compared with native ones can potentially be generated by such methods, but successful design is rare because preparing a sequence library by curating the database and selecting a method is difficult. Utilizing a curated library prepared by reducing conservation energies, we successfully designed two artificial l-threonine 3-dehydrogenases (SDR-TDH) with higher activity levels than native SDR-TDH, FcTDH-N1, and AncTDH, using FCD and ASR, respectively. The artificial SDR-TDHs had excellent thermal stability and NAD⁺ recognition compared to native SDR-TDH from Cupriavidus necator (CnTDH); the melting temperatures of FcTDH-N1 and AncTDH were about 10 and 5 °C higher than that of CnTDH, respectively, and the dissociation constants toward NAD⁺ of FcTDH-N1 and AncTDH were 2- and 7-fold lower than that of CnTDH, respectively. Enzymatic efficiency of the artificial SDR-TDHs were comparable to that of CnTDH. Crystal structures of FcTDH-N1 and AncTDH were determined at 2.8 and 2.1 Å resolution, respectively. Structural and MD simulation analysis of the SDR-TDHs indicated that only the flexibility at specific regions was changed, suggesting that multiple mutations introduced in the artificial SDR-TDHs altered their flexibility and thereby affected their enzymatic properties. Benchmark analysis of the SDR-TDHs indicated that both FCD and ASR can generate highly functional proteins if a curated library is prepared appropriately.

Structural Insight into Anaphase Promoting Complex 3 Structure and Docking with a Natural Inhibitory Compound

BACKGROUND

Anaphase promoting complex (APC) is the biggest Cullin-RING E3 ligase and is very important in cell cycle control; many anti-cancer agents target this. APC controls the onset of chromosome separation and mitotic exit through securin and cyclin B degradation, respectively. Its APC3 subunit identifies the APC activators-Cdh1 and Cdc20.

MATERIALS AND METHODS

The structural model of the APC3 subunit of APC was developed by means of computational techniques; the binding of a natural inhibitory compound to APC3 was also investigated.

RESULTS

It was found that APC3 structure consists of numerous helices organized in anti-parallel and the overall model is superhelical of tetratrico-peptide repeat (TPR) domains. Furthermore, binding pocket of the natural inhibitory compound as APC3 inhibitor was shown.

CONCLUSION

The findings are beneficial to understand the mechanism of the APC activation and design inhibitory compounds.

Intensification: A Resource for Amplifying Population-Genetic Signals with Protein Repeats

Large-scale genome sequencing holds great promise for the interpretation of protein structures through the discovery of many, rare functional variants in the human population. However, because protein-coding regions are under high selective constraints, these variants occur at low frequencies, such that there is often insufficient statistics for downstream calculations. To address this problem, we develop the Intensification approach, which uses the modular structure of repeat protein domains to amplify signals of selection from population genetics and traditional interspecies conservation. In particular, we are able to aggregate variants at the codon level to identify important positions in repeat domains that show strong conservation signals. This allows us to compare conservation over different evolutionary timescales. It also enables us to visualize population-genetic measures on protein structures. We make available the Intensification results as an online resource (http://intensification.gersteinlab.org) and illustrate the approach through a case study on the tetratricopeptide repeat.

Origin of a folded repeat protein from an intrinsically disordered ancestor

Repetitive proteins are thought to have arisen through the amplification of subdomain-sized peptides. Many of these originated in a non-repetitive context as cofactors of RNA-based replication and catalysis, and required the RNA to assume their active conformation. In search of the origins of one of the most widespread repeat protein families, the tetratricopeptide repeat (TPR), we identified several potential homologs of its repeated helical hairpin in non-repetitive proteins, including the putatively ancient ribosomal protein S20 (RPS20), which only becomes structured in the context of the ribosome. We evaluated the ability of the RPS20 hairpin to form a TPR fold by amplification and obtained structures identical to natural TPRs for variants with 2-5 point mutations per repeat. The mutations were neutral in the parent organism, suggesting that they could have been sampled in the course of evolution. TPRs could thus have plausibly arisen by amplification from an ancestral helical hairpin.

Combining Design and Selection to Create Novel Protein-Peptide Interactions

The ability to design new protein-protein interactions (PPIs) has many applications in biotechnology and medicine. The goal of designed PPIs is to achieve both high affinity and specificity for the target protein. A great challenge in protein design is to identify such proteins from an enormous number of potential sequences. Many computational and experimental methods have been developed to contend with this challenge. Here we describe one particularly powerful approach-semirational design-that combines design and selection. This approach has been applied to generate new PPIs for many applications, including novel affinity reagents for protein detection/purification and bioorthogonal modules for synthetic biology (Jackrel, Valverde, & Regan, 2009; Sawyer et al., 2014; Speltz, Brown, Hajare, Schlieker, & Regan, 2015; Speltz, Nathan, & Regan, 2015).

Protein stability: computation, sequence statistics, and new experimental methods

Calculating protein stability and predicting stabilizing mutations remain exceedingly difficult tasks, largely due to the inadequacy of potential functions, the difficulty of modeling entropy and the unfolded state, and challenges of sampling, particularly of backbone conformations. Yet, computational design has produced some remarkably stable proteins in recent years, apparently owing to near ideality in structure and sequence features. With caveats, computational prediction of stability can be used to guide mutation, and mutations derived from consensus sequence analysis, especially improved by recent co-variation filters, are very likely to stabilize without sacrificing function. The combination of computational and statistical approaches with library approaches, including new technologies such as deep sequencing and high throughput stability measurements, point to a very exciting near term future for stability engineering, even with difficult computational issues remaining.

Biomolecular templating of functional hybrid nanostructures using repeat protein scaffolds

The precise synthesis of materials and devices with tailored complex structures and properties is a requisite for the development of the next generation of products based on nanotechnology. Nowadays, the technology for the generation of this type of devices lacks the precision to determine their properties and is accomplished mostly by 'trial and error' experimental approaches. The use of bottom-up approaches that rely on highly specific biomolecular interactions of small and simple components is an attractive approach for the templating of nanoscale elements. In nature, protein assemblies define complex structures and functions. Engineering novel bio-inspired assemblies by exploiting the same rules and interactions that encode the natural diversity is an emerging field that opens the door to create nanostructures with numerous potential applications in synthetic biology and nanotechnology. Self-assembly of biological molecules into defined functional structures has a tremendous potential in nano-patterning and the design of novel materials and functional devices. Molecular self-assembly is a process by which complex 3D structures with specified functions are constructed from simple molecular building blocks. Here we discuss the basis of biomolecular templating, the great potential of repeat proteins as building blocks for biomolecular templating and nano-patterning. In particular, we focus on the designed consensus tetratricopeptide repeats (CTPRs), the control on the assembly of these proteins into higher order structures and their potential as building blocks in order to generate functional nanostructures and materials.

Detection and sequence/structure mapping of biophysical constraints to protein variation in saturated mutational libraries and protein sequence alignments with a dedicated server

Background

Protein variability can now be studied by measuring high-resolution tolerance-to-substitution maps and fitness landscapes in saturated mutational libraries. But these rich and expensive datasets are typically interpreted coarsely, restricting detailed analyses to positions of extremely high or low variability or dubbed important beforehand based on existing knowledge about active sites, interaction surfaces, (de)stabilizing mutations, etc.

Results

Our new webserver PsychoProt (freely available without registration at http://psychoprot.epfl.ch or at http://lucianoabriata.altervista.org/psychoprot/index.html) helps to detect, quantify, and sequence/structure map the biophysical and biochemical traits that shape amino acid preferences throughout a protein as determined by deep-sequencing of saturated mutational libraries or from large alignments of naturally occurring variants.

Discussion

We exemplify how PsychoProt helps to (i) unveil protein structure-function relationships from experiments and from alignments that are consistent with structures according to coevolution analysis, (ii) recall global information about structural and functional features and identify hitherto unknown constraints to variation in alignments, and (iii) point at different sources of variation among related experimental datasets or between experimental and alignment-based data. Remarkably, metabolic costs of the amino acids pose strong constraints to variability at protein surfaces in nature but not in the laboratory. This and other differences call for caution when extrapolating results from in vitro experiments to natural scenarios in, for example, studies of protein evolution.

Conclusion

We show through examples how PsychoProt can be a useful tool for the broad communities of structural biology and molecular evolution, particularly for studies about protein modeling, evolution and design.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1124-4) contains supplementary material, which is available to authorized users.

Consensus protein design

A popular and successful strategy in semi-rational design of protein stability is the use of evolutionary information encapsulated in homologous protein sequences. Consensus design is based on the hypothesis that at a given position, the respective consensus amino acid contributes more than average to the stability of the protein than non-conserved amino acids. Here, we review the consensus design approach, its theoretical underpinnings, successes, limitations and challenges, as well as providing a detailed guide to its application in protein engineering.

Design of Protein-Peptide Interaction Modules for Assembling Supramolecular Structures in Vivo and in Vitro

Synthetic biology and protein origami both require protein building blocks that behave in a reliable, predictable fashion. In particular, we require protein interaction modules with known specificity and affinity. Here, we describe three designed TRAP (Tetratricopeptide Repeat Affinity Protein)-peptide interaction pairs that are functional in vivo. We show that each TRAP binds to its cognate peptide and exhibits low cross-reactivity with the peptides bound by the other TRAPs. In addition, we demonstrate that the TRAP-peptide interactions are functional in many cellular contexts. In extensions of these designs, we show that the binding affinity of a TRAP-peptide pair can be systematically varied. The TRAP-peptide pairs we present thus represent a powerful set of new building blocks that are suitable for a variety of applications.

Protein evolution analysis of S-hydroxynitrile lyase by complete sequence design utilizing the INTMSAlign software

Development of software and methods for design of complete sequences of functional proteins could contribute to studies of protein engineering and protein evolution. To this end, we developed the INTMSAlign software, and used it to design functional proteins and evaluate their usefulness. The software could assign both consensus and correlation residues of target proteins. We generated three protein sequences with S-selective hydroxynitrile lyase (S-HNL) activity, which we call designed S-HNLs; these proteins folded as efficiently as the native S-HNL. Sequence and biochemical analysis of the designed S-HNLs suggested that accumulation of neutral mutations occurs during the process of S-HNLs evolution from a low-activity form to a high-activity (native) form. Taken together, our results demonstrate that our software and the associated methods could be applied not only to design of complete sequences, but also to predictions of protein evolution, especially within families such as esterases and S-HNLs.

A general computational approach for repeat protein design

Repeat proteins have considerable potential for use as modular binding reagents or biomaterials in biomedical and nanotechnology applications. Here we describe a general computational method for building idealized repeats that integrates available family sequences and structural information with Rosetta de novo protein design calculations. Idealized designs from six different repeat families were generated and experimentally characterized; 80% of the proteins were expressed and soluble and more than 40% were folded and monomeric with high thermal stability. Crystal structures determined for members of three families are within 1Å root-mean-square deviation to the design models. The method provides a general approach for fast and reliable generation of stable modular repeat protein scaffolds.

Crystallization and preliminary X-ray crystallographic analysis of an artificial molten-globular-like triosephosphate isomerase protein of mixed phylogenetic origin

A bioinformatics-based protein-engineering approach called consensus design led to the construction of a chimeric triosephosphate isomerase (TIM) protein called ccTIM (curated consensus TIM) which is as active as Saccharomyces cerevisiae TIM despite sharing only 58% sequence identity with it. The amino-acid sequence of this novel protein is as identical to native sequences from eukaryotes as to those from prokaryotes and shares some biophysical traits with a molten globular protein. Solving its crystal structure would help in understanding the physical implications of its bioinformatics-based sequence. In this report, the ccTIM protein was successfully crystallized using the microbatch-under-oil method and a full X-ray diffraction data set was collected to 2.2 Å resolution using a synchrotron-radiation source. The crystals belonged to space group C2221, with unit-cell parameters a=107.97, b=187.21, c=288.22 Å. Matthews coefficient calculations indicated the presence of six dimers in the asymmetric unit, with an approximate solvent content of 46.2%.

Consensus design of a NOD receptor leucine rich repeat domain with binding affinity for a muramyl dipeptide, a bacterial cell wall fragment

Repeat proteins have recently emerged as especially well-suited alternative binding scaffolds due to their modular architecture and biophysical properties. Here we present the design of a scaffold based on the consensus sequence of the leucine rich repeat (LRR) domain of the NOD family of cytoplasmic innate immune system receptors. Consensus sequence design has emerged as a protein design tool to create de novo proteins that capture sequence-structure relationships and interactions present in nature. The multiple sequence alignment of 311 individual LRRs, which are the putative ligand-recognition domain in NOD proteins, resulted in a consensus sequence protein containing two internal and N- and C-capping repeats named CLRR2. CLRR2 protein is a stable, monomeric, and cysteine free scaffold that without any affinity maturation displays micromolar binding to muramyl dipeptide, a bacterial cell wall fragment. To our knowledge, this is the first report of direct interaction of a NOD LRR with a physiologically relevant ligand.

Functional insight from the tetratricopeptide repeat-like motifs of the type III secretion chaperone SicA in Salmonella enterica serovar Typhimurium

SicA functions both as a class II chaperone for SipB and SipC of the type III secretion system (T3SS)-1 and as a transcriptional cofactor for the AraC-type transcription factor InvF in Salmonella enterica subsp. enterica serovar Typhimurium. Bioinformatic analysis has predicted that SicA possesses three tetratricopeptide repeat (TPR)-like motifs, which are important for protein-protein interactions and serve as multiprotein complex mediators. To investigate whether the TPR-like motifs in SicA are critical for its transcriptional cofactor function, the canonical residues in these motifs were mutated to glutamate (SicAA44E , SicAA78E , and SicAG112E ). None of these mutants except SicAA44E were able to activate the expression of the sipB and sigD genes. SicAA44E still has a capacity to interact with InvF in vitro, and despite its instability in cell, it could activate the sigDE operon. This suggests that TPR motifs are important for the transcriptional cofactor function of the SicA chaperone.

Highly covarying residues have a functional role in antibody constant domains

The ability to generate and design antibodies recognizing specific targets has revolutionized the pharmaceutical industry and medical imaging. Engineering antibody therapeutics in some cases requires modifying their constant domains to enable new and altered interactions. Engineering novel specificities into antibody constant domains has proved challenging due to the complexity of inter-domain interactions. Covarying networks of residues that tend to cluster on the protein surface and near binding sites have been identified in some proteins. However, the underlying role these networks play in the protein resulting in their conservation remains unclear in most cases. Resolving their role is crucial, because residues in these networks are not viable design targets if their role is to maintain the fold of the protein. Conversely, these networks of residues are ideal candidates for manipulating specificity if they are primarily involved in binding, such as the myriad interdomain interactions maintained within antibodies. Here, we identify networks of evolutionarily-related residues in C-class antibody domains by evaluating covariation, a measure of propensity with which residue pairs vary dependently during evolution. We computationally test whether mutation of residues in these networks affects stability of the folded antibody domain, determining their viability as design candidates. We find that members of covarying networks cluster at domain-domain interfaces, and that mutations to these residues are diverse and frequent during evolution, precluding their importance to domain stability. These results indicate that networks of covarying residues exist in antibody domains for functional reasons unrelated to thermodynamic stability, making them ideal targets for antibody design.

All repeats are not equal: a module-based approach to guide repeat protein design

Repeat proteins composed of tandem arrays of a short structural motif often mediate protein-protein interactions. Past efforts to design repeat protein-based molecular recognition tools have focused on the creation of templates from the consensus of individual repeats, regardless of their natural context. Such an approach assumes that all repeats are essentially equivalent. In this study, we present the results of a "module-based" approach in which modules composed of tandem repeats are aligned to identify repeat-specific features. Using this approach to analyze tetratricopeptide repeat modules that contain three tandem repeats (3TPRs), we identify two classes of 3TPR modules with distinct structural signatures that are correlated with different sets of functional residues. Our analyses also reveal a high degree of correlation between positions across the entire ligand-binding surface, indicative of a coordinated, coevolving binding surface. Extension of our analyses to different repeat protein modules reveals more examples of repeat-specific features, especially in armadillo repeat modules. In summary, the module-based analyses that we present effectively capture key repeat-specific features that will be important to include in future repeat protein design templates.

Protein engineering and stabilization from sequence statistics: variation and covariation analysis

The concepts of consensus and correlation in multiple sequence alignments (MSAs) have been used in the past to understand and engineer proteins. However, there are multiple ways of acquiring MSA databases and also numerous mathematical metrics that can be applied to calculate each of the parameters. This chapter describes an overall methodology that we have chosen to employ for acquiring and statistically analyzing MSAs. We have provided a step-by-step protocol for calculating relative entropy and mutual information metrics and describe how they can be used to predict mutations that have a high probability of stabilizing a protein. This protocol allows for flexibility for modification of formulae and parameters without using anything more complicated than Microsoft Excel. We have also demonstrated various aspects of data analysis by carrying out a sample analysis on the BPTI-Kunitz family of proteins and identified mutations that would be predicted to stabilize this protein based on consensus and correlation values.

The complex interactions of Chs5p, the ChAPs, and the cargo Chs3p

The exomer complex, consisting of ChAPs and Chs5p, exports specialized cargoes from the TGN. ChAPs bind to Chs5p through TPR repeats, whereas cargo specificity of the ChAPs is outside these interaction modules. Chs3p and Chs6p may require a complex interaction to form a complex.

The exomer complex is a putative vesicle coat required for the direct transport of a subset of cargoes from the trans-Golgi network (TGN) to the plasma membrane. Exomer comprises Chs5p and the ChAPs family of proteins (Chs6p, Bud7p, Bch1p, and Bch2p), which are believed to act as cargo receptors. In particular, Chs6p is required for the transport of the chitin synthase Chs3p to the bud neck. However, how the ChAPs associate with Chs5p and recognize cargo is not well understood. Using domain-switch chimeras of Chs6p and Bch2p, we show that four tetratricopeptide repeats (TPRs) are involved in interaction with Chs5p. Because these roles are conserved among the ChAPs, the TPRs are interchangeable among different ChAP proteins. In contrast, the N-terminal and the central parts of the ChAPs contribute to cargo specificity. Although the entire N-terminal domain of Chs6p is required for Chs3p export at all cell cycle stages, the central part seems to predominantly favor Chs3p export in small-budded cells. The cargo Chs3p probably also uses a complex motif for the interaction with Chs6, as the C-terminus of Chs3p interacts with Chs6p and is necessary, but not sufficient, for TGN export.

Exploring the molecular linkage of protein stability traits for enzyme optimization by iterative truncation and evolution

The stability of proteins is paramount for their therapeutic and industrial use and, thus, is a major task for protein engineering. Several types of chemical and physical stabilities are desired, and discussion revolves around whether each stability trait needs to be addressed separately and how specific and compatible stabilizing mutations act. We demonstrate a stepwise perturbation-compensation strategy, which identifies mutations rescuing the activity of a truncated TEM β-lactamase. Analyses relating structural stress with the external stresses of heat, denaturants, and proteases reveal our second-site suppressors as general stability centers that also improve the full-length enzyme. A library of lactamase variants truncated by 15 N-terminal and three C-terminal residues (Bla-NΔ15CΔ3) was subjected to activity selection and DNA shuffling. The resulting clone with the best in vivo performance harbored eight mutations, surpassed the full-length wild-type protein by 5.3 °C in T(m), displayed significantly higher catalytic activity at elevated temperatures, and showed delayed guanidine-induced denaturation. The crystal structure of this mutant was determined and provided insights into its stability determinants. Stepwise reconstitution of the N- and C-termini increased its thermal, denaturant, and proteolytic resistance successively, leading to a full-length enzyme with a T(m) increased by 15.3 °C and a half-denaturation concentration shifted from 0.53 to 1.75 M guanidinium relative to that of the wild type. These improvements demonstrate that iterative truncation-optimization cycles can exploit stability-trait linkages in proteins and are exceptionally suited for the creation of progressively stabilized variants and/or downsized proteins without the need for detailed structural or mechanistic information.

From artificial antibodies to nanosprings: the biophysical properties of repeat proteins

In this chapter we review recent studies of repeat proteins, a class of proteins consisting of tandem arrays of small structural motifs that stack approximately linearly to produce elongated structures. We discuss the observation that, despite lacking the long-range tertiary interactions that are thought to be the hallmark of globular protein stability, repeat proteins can be as stable and as co-orperatively folded as their globular counterparts. The symmetry inherent in the structures of repeat arrays, however, means there can be many partly folded species (whether it be intermediates or transition states) that have similar stabilities. Consequently they do have distinct folding properties compared with globular proteins and these are manifest in their behaviour both at equilibrium and under kinetic conditions. Thus, when studying repeat proteins one appears to be probing a moving target: a relatively small perturbation, by mutation for example, can result in a shift to a different intermediate or transition state. The growing literature on these proteins illustrates how their modular architecture can be adapted to a remarkable array of biological and physical roles, both in vivo and in vitro. Further, their simple architecture makes them uniquely amenable to redesign-of their stability, folding and function-promising exciting possibilities for future research.

Structure of minimal tetratricopeptide repeat domain protein Tah1 reveals mechanism of its interaction with Pih1 and Hsp90

Tah1 and Pih1 are novel Hsp90 interactors. Tah1 acts as a cofactor of Hsp90 to stabilize Pih1. In yeast, Hsp90, Tah1, and Pih1 were found to form a complex that is required for ribosomal RNA processing through their effect on box C/D small nucleolar ribonucleoprotein assembly. Tah1 is a minimal tetratricopeptide repeat protein of 111 amino acid residues that binds to the C terminus of the Hsp90 molecular chaperone, whereas Pih1 consists of 344 residues of unknown fold. The NMR structure of Tah1 has been solved, and this structure shows the presence of two tetratricopeptide repeat motifs followed by a C helix and an unstructured region. The binding of Tah1 to Hsp90 is mediated by the EEVD C-terminal residues of Hsp90, which bind to a positively charged channel formed by Tah1. Five highly conserved residues, which form a two-carboxylate clamp that tightly interacts with the ultimate Asp-0 residue of the bound peptide, are also present in Tah1. Tah1 was found to bind to the C terminus of Pih1 through the C helix and the unstructured region. The C terminus of Pih1 destabilizes the protein in vitro and in vivo, whereas the binding of Tah1 to Pih1 allows for the formation of a stable complex. Based on our data, a model for an Hsp90-Tah1-Pih1 ternary complex is proposed.

Modulating repeat protein stability: the effect of individual helix stability on the collective behavior of the ensemble

Repeat proteins are tandem arrays of a small structural motif, in which tertiary structure is stabilized by interactions within a repeat and between neighboring repeats. Several studies have shown that this modular structure is manifest in modular thermodynamic properties. Specifically, the global stability of a repeat protein can be described by simple linear models, considering only two parameters: the stability of the individual repeated units (H) and the coupling interaction between the units (J). If the repeat units are identical, single values of H and J, together with the number of repeated units, is sufficient to completely describe the thermodynamic behavior of any protein within a series. In this work, we demonstrate how the global stability of a repeat protein can be changed, in a predictable fashion, by modifying only the H parameter. Taking a previously characterized series of consensus tetratricopeptide repeats (TPR) (CTPRa) proteins, we introduced mutations into the basic repeating unit, such that the stability of the individual repeat unit was increased, but its interaction with neighboring units was unchanged. In other words, we increased H but kept J constant. We demonstrated that the denaturation curves for a series of such repeat proteins can be fit and additional curves can be predicted by the one-dimensional Ising model in which only H has changed from the original fit for the CTPRa series. Our results show that we can significantly increase the stability of a repeat protein by rationally increasing the stability of the units (H), whereas the interaction between repeats (J) remains unchanged.

Protein stability by number: high-throughput and statistical approaches to one of protein science's most difficult problems

Most proteins are only barely stable, which impedes research, complicates therapeutic applications, and makes proteins susceptible to pathologically destabilizing mutations. Our ability to predict the thermodynamic consequences of even single point mutations is still surprisingly limited, and established methods of measuring stability are slow. Recent advances are bringing protein stability studies into the high-throughput realm. Some methods are based on inferential read-outs such as activity, proteolytic resistance or split-protein fragment reassembly. Other methods use miniaturization of direct measurements, such as intrinsic fluorescence, H/D exchange, cysteine reactivity, aggregation and hydrophobic dye binding (DSF). Protein engineering based on statistical analysis (consensus and correlated occurrences of amino acids) is promising, but much work remains to understand and implement these methods.

Designed proteins to modulate cellular networks

A major challenge of protein design is to create useful new proteins that interact specifically with biological targets in living cells. Such binding modules have many potential applications, including the targeted perturbation of protein networks. As a general approach to create such modules, we designed a library with approximately 10(9) different binding specificities based on a small 3-tetratricopeptide repeat (TPR) motif framework. We employed a novel strategy, based on split GFP reassembly, to screen the library for modules with the desired binding specificity. Using this approach, we identified modules that bind tightly and specifically to Dss1, a small human protein that interacts with the tumor suppressor protein BRCA2. We showed that these modules also bind the yeast homologue of Dss1, Sem1. Furthermore, we demonstrated that these modules inhibit Sem1 activity in yeast. This strategy will be generally applicable to make novel genetically encoded tools for systems/synthetic biology applications.

Screening libraries to identify proteins with desired binding activities using a split-GFP reassembly assay

Designer protein modules, which bind specifically to a desired target, have numerous potential applications. One approach to creating such proteins is to construct and screen libraries. Here we present a detailed description of using a split-GFP reassembly assay to screen libraries and identify proteins with novel binding properties. Attractive features of the split-GFP based screen are the absence of false positives and the simplicity, robustness, and ease of automation of the screen. Here, we describe both the construction of a naive protein library, and screening of the library using the split-GFP assay to identify proteins that bind specifically to chosen peptide sequences.

Creating novel proteins by combining design and selection

We present the results of combining design and selection to remodel a protein-peptide binding interface, using the peptide PTIEEVD and the TPR1 module interaction as our test case. We initially used the program Rosetta to interrogate possible TPR1 sequences compatible with binding the peptide PTIEEVD. Based on these results, we screened a small library of TPR1 variants, using a split GFP fluorescent assay to identify proteins that are able to bind to the PTIEEVD peptide. We discuss the similarities and differences between the modeling and selection results at each position. We show that a new 'consensus' TPR1, created based on the results of the sequences identified in the screen, indeed binds to the PTIEEVD peptide. These results demonstrate the utility of combining design and selection in a synergistic fashion to remodel protein recognition interfaces.

Calcium-induced folding of a beta roll motif requires C-terminal entropic stabilization

Beta roll motifs are associated with several proteins secreted by the type 1 secretion system (T1SS). Located just upstream of the C-terminal T1SS secretion signal, they are believed to act as calcium-induced switches that prevent folding before secretion. Bordetella pertussis adenylate cyclase (CyaA) toxin has five blocks of beta roll motifs (or repeats-in-toxin motifs) separated by linkers. The block V motif on its own has been reported to be non-responsive to calcium. Only when the N- and C-terminal linkers, or flanking groups, were fused did the motif bind calcium and fold. In an effort to understand the requirements for beta roll folding, we have truncated the N- and C-terminal flanks at several locations to determine the minimal essential sequences. Calcium-responsive beta roll folding occurred even in the absence of the natural N-terminal flank. The natural C-terminal flank could not be truncated without decreased calcium affinity and only partially truncated before losing calcium-responsiveness. Globular protein fusion at the C-terminus likewise enabled calcium-induced folding but fusions solely at the N-terminus failed. This demonstrates that calcium-induced folding is an inherent property of the beta roll motif rather than the flanking groups. Given the disparate nature of the observed functional flanking groups, C-terminal fusions appear to confer calcium-responsiveness to the beta roll motif via a non-specific mechanism, suggesting that entropic stabilization of the unstructured C-terminus can enable beta roll folding. Increased calcium affinity was observed when the natural C-terminal flank was used to enable calcium-induced folding, pointing to its cooperative participation in beta roll formation. This work indicates that a general principle of C-terminal entropic stabilization can enable stimulus-responsive repeat protein folding, while the C-terminal flank has a specific role in tuning calcium-responsive beta roll formation. These observations are in stark contrast to what has been reported for other repeat proteins.

Consensus protein design without phylogenetic bias

Consensus design is an appealing strategy for the stabilization of proteins. It exploits amino acid conservation in sets of homologous proteins to identify likely beneficial mutations. Nevertheless, its success depends on the phylogenetic diversity of the sequence set available. Here, we show that randomization of a single protein represents a reliable alternative source of sequence diversity that is essentially free of phylogenetic bias. A small number of functional protein sequences selected from binary-patterned libraries suffice as input for the consensus design of active enzymes that are easier to produce and substantially more stable than individual members of the starting data set. Although catalytic activity correlates less consistently with sequence conservation in these extensively randomized proteins, less extreme mutagenesis strategies might be adopted in practice to augment stability while maintaining function.

Distribution of activator of G-protein signaling 3 within the aggresomal pathway: role of specific residues in the tetratricopeptide repeat domain and differential regulation by the AGS3 binding partners Gi(alpha) and mammalian inscuteable

AGS3, a receptor-independent activator of G-protein signaling, is involved in unexpected functional diversity for G-protein signaling systems. AGS3 has seven tetratricopeptide (TPR) motifs upstream of four G-protein regulatory (GPR) motifs that serve as docking sites for Gialpha-GDP. The positioning of AGS3 within the cell and the intramolecular dynamics between different domains of the proteins are likely key determinants of their ability to influence G-protein signaling. We report that AGS3 enters into the aggresome pathway and that distribution of the protein is regulated by the AGS3 binding partners Gialpha and mammalian Inscuteable (mInsc). Gialpha rescues AGS3 from the aggresome, whereas mInsc augments the aggresome-like distribution of AGS3. The distribution of AGS3 to the aggresome is dependent upon the TPR domain, and it is accelerated by disruption of the TPR organizational structure or introduction of a nonsynonymous single-nucleotide polymorphism. These data present AGS3, G-proteins, and mInsc as candidate proteins involved in regulating cellular stress associated with protein-processing pathologies.