Reverse engineering the (beta/alpha )8 barrel fold

Design and Selection of Heterodimerizing Helical Hairpins for Synthetic Biology

Synthetic biology applications would benefit from protein modules of reduced complexity that function orthogonally to cellular components. As many subcellular processes depend on peptide-protein or protein-protein interactions, de novo designed polypeptides that can bring together other proteins controllably are particularly useful. Thanks to established sequence-to-structure relationships, helical bundles provide good starting points for such designs. Typically, however, such designs are tested in vitro and function in cells is not guaranteed. Here, we describe the design, characterization, and application of de novo helical hairpins that heterodimerize to form 4-helix bundles in cells. Starting from a rationally designed homodimer, we construct a library of helical hairpins and identify complementary pairs using bimolecular fluorescence complementation in E. coli. We characterize some of the pairs using biophysics and X-ray crystallography to confirm heterodimeric 4-helix bundles. Finally, we demonstrate the function of an exemplar pair in regulating transcription in both E. coli and mammalian cells.

De Novo Design of a Highly Stable Ovoid TIM Barrel: Unlocking Pocket Shape towards Functional Design

The ability to finely control the structure of protein folds is an important prerequisite to functional protein design. The TIM barrel fold is an important target for these efforts as it is highly enriched for diverse functions in nature. Although a TIM barrel protein has been designed de novo, the ability to finely alter the curvature of the central beta barrel and the overall architecture of the fold remains elusive, limiting its utility for functional design. Here, we report the de novo design of a TIM barrel with ovoid (twofold) symmetry, drawing inspiration from natural beta and TIM barrels with ovoid curvature. We use an autoregressive backbone sampling strategy to implement our hypothesis for elongated barrel curvature, followed by an iterative enrichment sequence design protocol to obtain sequences which yield a high proportion of successfully folding designs. Designed sequences are highly stable and fold to the designed barrel curvature as determined by a 2.1 Å resolution crystal structure. The designs show robustness to drastic mutations, retaining high melting temperatures even when multiple charged residues are buried in the hydrophobic core or when the hydrophobic core is ablated to alanine. As a scaffold with a greater capacity for hosting diverse hydrogen bonding networks and installation of binding pockets or active sites, the ovoid TIM barrel represents a major step towards the de novo design of functional TIM barrels.

Simultaneously improving the specific activity and thermostability of α-amylase BLA by rational design

Higher activity and alkaline α-amylases are desired for textile desizing and detergent additive. Here, rational design was used to improve the specific activity and thermostability of the α-amylase BLA from Bacillus licheniformis. Seventeen mutants of BLA were designed based on sequence consensus analysis and folding free energy calculation, and characterized by measuring their respective activity and thermostability at pH 8.5. Among them, mutant Q360C exhibited nearly threefold improved activity than that of wild-type and retained a higher residual activity (75% vs 59% for wild-type) after preincubation at 70 ℃ for 30 min. The modeled structures and molecular dynamics simulations analysis demonstrated that the enhanced hydrophobic interaction near residue 360 and reduced disturbance to the conformation of catalytic residues are the possible reasons for the improved thermostability and activity of Q360C. The results suggest that 360th of BLA may act as a hotspot for engineering other enzymes in the GH13 superfamily.

Structural and Biochemical Characterization of Endo-β-1,4-glucanase from Dictyoglomus thermophilum, a Hyperthermostable and Halotolerant Cellulase

Enzymatic conversion of polysaccharides in the lignocellulosic biomass is currently the subject of intensive research and will be a key technology in future biorefineries. Using a bioinformatics approach, we previously identified a putative endo-β-1,4-glucanase (DtCel5A) from Dictyoglomus thermophilum, a chemoorganotrophic and thermophilic bacterium. Here, we structurally and functionally characterize DtCel5A and show that it is endowed with remarkable thermal and chemical stability. The structural features of DtCel5A and of its complex with cellobiose have been investigated by combining X-ray crystallography and other biophysical studies. Importantly, biochemical assays show that DtCel5A retains its activity on cellulose at high temperatures and at elevated salt concentrations. These features make DtCel5A an enzyme with interesting biotechnological applications for biomass degradation.

Triosephosphate isomerase deficiency: Effect of F240L mutation on enzyme structure

Eleven missense mutations have been describe in human triosephosphate isomerase (TPI), affecting its catalytic function. Several of these mutations generate triosephosphate isomerase deficiency, the consequences of which can in some cases be lethal. The missense F240L mutation was found in a Hungarian patient showing symptoms of chronic hemolytic anemia and neuromuscular dysfunction. In vitro studies using a recombinant version of this mutant showed that it affects kinetic parameters, thermal stability and dimeric stability. Using X-ray crystal structures, the present paper describes how this mutation affected the flexibility of catalytic residues K13 and part of the (β/α) 8-barrel fold facing the dimeric interface in the TPI.

Engineering novel S-glycosidase activity into extremo-adapted β-glucosidase by rational design

The breakdown of sulphur glycosidic bonds in thioglycosides can produce isothiocyanate, a chemoprotective agent linked to the prevention of cancers; however, only a handful of enzymes have been identified that are k0nown to catalyse this reaction. Structural studies of the myrosinase enzyme, which is capable of hydrolysing the thioglycosidic bond, have identified residues that may play important roles in sulphur bond specific activity. Using rational design, two extremo-adapted β-glycosidases from the species Thermus nonproteolyticus (TnoGH1) and Halothermothrix orenii (HorGH1) were engineered towards thioglycoside substrates. Twelve variants, six for TnoGH1and six for HorGH1, were assayed for activity. Remarkable enhancement of the specificity (k_cat/K_M) of TnoGH1 and HorGH1 towards β-thioglycoside was observed in the single mutants TnoGH1-V287R (2500 M^-1 s^-1) and HorGH1-M229R (13,260 M^-1 s^-1) which showed a 3-fold increase with no loss in turnover rate when compared with the wild-type enzymes. Thus, the role of arginine is key to induce β-thioglycosidase activity. Thorough kinetic investigation of the different mutants has shed light on the mechanism of β-glycosidases when acting on the native substrate.Key Points •Key residues were identified in the active site of Brevicoryne brassicae myrosinase. •Rationally designed mutations were introduced into two extremo-adapted β-glycosidases. •β-glycosidases mutants exhibited improved activity against thioglycosidic bonds. •The mutation to arginine in the active site yielded the best variant.

Simultaneously Improved Thermostability and Hydrolytic Pattern of Alpha-Amylase by Engineering Central Beta Strands of TIM Barrel

This study reported simultaneously improved thermostability and hydrolytic pattern of α-amylase from Bacillus subtilis CN7 by rationally engineering the mostly conserved central beta strands in TIM barrel fold. Nine single point mutations and a double mutation were introduced at the 2nd site of the β7 strand and 3rd site of the β5 strand to rationalize the weak interactions in the beta strands of the TIM barrel of α-amylase. All the five active mutants changed the compositions and percentages of maltooligosaccharides in final hydrolytic products compared to the product spectrum of the wild-type. A mutant Y204V produced only maltose, maltotriose, and maltopentaose without any glucose and maltotetraose, indicating a conversion from typical endo-amylase to novel maltooligosaccharide-producing amylase. A mutant V260I enhanced the thermal stability by 7.1 °C. To our best knowledge, this is the first report on the simultaneous improvement of thermostability and hydrolytic pattern of α-amylase by engineering central beta strands of TIM barrel and the novel "beta strands" strategy proposed here may be useful for the protein engineering of other TIM barrel proteins.

Reduced alphabet of prebiotic amino acids optimally encodes the conformational space of diverse extant protein folds

Background

There is wide agreement that only a subset of the twenty standard amino acids existed prebiotically in sufficient concentrations to form functional polypeptides. We ask how this subset, postulated as {A,D,E,G,I,L,P,S,T,V}, could have formed structures stable enough to found metabolic pathways. Inspired by alphabet reduction experiments, we undertook a computational analysis to measure the structural coding behavior of sequences simplified by reduced alphabets. We sought to discern characteristics of the prebiotic set that would endow it with unique properties relevant to structure, stability, and folding.

Results

Drawing on a large dataset of single-domain proteins, we employed an information-theoretic measure to assess how well the prebiotic amino acid set preserves fold information against all other possible ten-amino acid sets. An extensive virtual mutagenesis procedure revealed that the prebiotic set excellently preserves sequence-dependent information regarding both backbone conformation and tertiary contact matrix of proteins. We observed that information retention is fold-class dependent: the prebiotic set sufficiently encodes the structure space of α/β and α + β folds, and to a lesser extent, of all-α and all-β folds. The prebiotic set appeared insufficient to encode the small proteins. Assessing how well the prebiotic set discriminates native vs. incorrect sequence-structure matches, we found that α/β and α + β folds exhibit more pronounced energy gaps with the prebiotic set than with nearly all alternatives.

Conclusions

The prebiotic set optimally encodes local backbone structures that appear in the folded environment and near-optimally encodes the tertiary contact matrix of extant proteins. The fold-class-specific patterns observed from our structural analysis confirm the postulated timeline of fold appearance in proteogenesis derived from proteomic sequence analyses. Polypeptides arising in a prebiotic environment will likely form α/β and α + β-like folds if any at all. We infer that the progressive expansion of the alphabet allowed the increased conformational stability and functional specificity of later folds, including all-α, all-β, and small proteins. Our results suggest that prebiotic sequences are amenable to mutations that significantly lower native conformational energies and increase discrimination amidst incorrect folds. This property may have assisted the genesis of functional proto-enzymes prior to the expansion of the full amino acid alphabet.

Missense variant in TPI1 (Arg189Gln) causes neurologic deficits through structural changes in the triosephosphate isomerase catalytic site and reduced enzyme levels in vivo

Mutations in the gene triosephosphate isomerase (TPI) lead to a severe multisystem condition that is characterized by hemolytic anemia, a weakened immune system, and significant neurologic symptoms such as seizures, distal neuropathy, and intellectual disability. No effective therapy is available. Here we report a compound heterozygous patient with a novel TPI pathogenic variant (NM_000365.5:c.569G>A:p.(Arg189Gln)) in combination with the common (NM_000365.5:c.315G>C:p.(Glu104Asp)) allele. We characterized the novel variant by mutating the homologous Arg in Drosophila using a genomic engineering system, demonstrating that missense mutations at this position cause a strong loss of function. Compound heterozygote animals were generated and exhibit motor behavioural deficits and markedly reduced protein levels. Furthermore, examinations of the TPI^Arg189Gln/TPI^Glu104Asp patient fibroblasts confirmed the reduction of TPI levels, suggesting that Arg189Gln may also affect the stability of the protein. The Arg189 residue participates in two salt bridges on the backside of the TPI enzyme dimer, and we reveal that a mutation at this position alters the coordination of the substrate-binding site and important catalytic residues. Collectively, these data reveal a new human pathogenic variant associated with TPI deficiency, identify the Arg189 salt bridge as critical for organizing the catalytic site of the TPI enzyme, and demonstrates that reduced TPI levels are associated with human TPI deficiency. These findings advance our understanding of the molecular pathogenesis of the disease, and suggest new therapeutic avenues for pre-clinical trials.

Flagellin glycosylation with pseudaminic acid in Campylobacter and Helicobacter: prospects for development of novel therapeutics

Many pathogenic bacteria require flagella-mediated motility to colonise and persist in their hosts. Helicobacter pylori and Campylobacter jejuni are flagellated epsilonproteobacteria associated with several human pathologies, including gastritis, acute diarrhea, gastric carcinoma and neurological disorders. In both species, glycosylation of flagellin with an unusual sugar pseudaminic acid (Pse) plays a crucial role in the biosynthesis of functional flagella, and thereby in bacterial motility and pathogenesis. Pse is found only in pathogenic bacteria. Its biosynthesis via six consecutive enzymatic steps has been extensively studied in H. pylori and C. jejuni. This review highlights the importance of flagella glycosylation and details structural insights into the enzymes in the Pse pathway obtained via a combination of biochemical, crystallographic, and mutagenesis studies of the enzyme-substrate and -inhibitor complexes. It is anticipated that understanding the underlying structural and molecular basis of the catalytic mechanisms of the Pse-synthesising enzymes will pave the way for the development of novel antimicrobials.

A guide to the effects of a large portion of the residues of triosephosphate isomerase on catalysis, stability, druggability, and human disease

Triosephosphate isomerase (TIM) is a ubiquitous enzyme, which appeared early in evolution. TIM is responsible for obtaining net ATP from glycolysis and producing an extra pyruvate molecule for each glucose molecule, under aerobic and anaerobic conditions. It is placed in a metabolic crossroad that allows a quick balance of the triose phosphate aldolase produced by glycolysis, and is also linked to lipid metabolism through the alternation of glycerol-3-phosphate and the pentose cycle. TIM is one of the most studied enzymes with more than 199 structures deposited in the PDB. The interest for this enzyme stems from the fact that it is involved in glycolysis, but also in aging, human diseases and metabolism. TIM has been a target in the search for chemical compounds against infectious diseases and is a model to study catalytic features. Until February 2017, 62% of all residues of the protein have been studied by mutagenesis and/or using other approaches. Here, we present a detailed and comprehensive recompilation of the reported effects on TIM catalysis, stability, druggability and human disease produced by each of the amino acids studied, contributing to a better understanding of the properties of this fundamental protein. The information reviewed here shows that the role of the noncatalytic residues depend on their molecular context, the delicate balance between the short and long-range interactions in concerted action determining the properties of the protein. Each protein should be regarded as a unique entity that has evolved to be functional in the organism to which it belongs. Proteins 2017; 85:1190-1211. © 2017 Wiley Periodicals, Inc.

Pairwise contact energy statistical potentials can help to find probability of point mutations

To adopt a particular fold, a protein requires several interactions between its amino acid residues. The energetic contribution of these residue-residue interactions can be approximated by extracting statistical potentials from known high resolution structures. Several methods based on statistical potentials extracted from unrelated proteins are found to make a better prediction of probability of point mutations. We postulate that the statistical potentials extracted from known structures of similar folds with varying sequence identity can be a powerful tool to examine probability of point mutation. By keeping this in mind, we have derived pairwise residue and atomic contact energy potentials for the different functional families that adopt the (α/β)₈ TIM-Barrel fold. We carried out computational point mutations at various conserved residue positions in yeast Triose phosphate isomerase enzyme for which experimental results are already reported. We have also performed molecular dynamics simulations on a subset of point mutants to make a comparative study. The difference in pairwise residue and atomic contact energy of wildtype and various point mutations reveals probability of mutations at a particular position. Interestingly, we found that our computational prediction agrees with the experimental studies of Silverman et al. (Proc Natl Acad Sci 2001;98:3092-3097) and perform better prediction than i_Mutant and Cologne University Protein Stability Analysis Tool. The present work thus suggests deriving pairwise contact energy potentials and molecular dynamics simulations of functionally important folds could help us to predict probability of point mutations which may ultimately reduce the time and cost of mutation experiments. Proteins 2016; 85:54-64. © 2016 Wiley Periodicals, Inc.

Crystal structures of two monomeric triosephosphate isomerase variants identified via a directed-evolution protocol selecting for L-arabinose isomerase activity

The crystal structures are described of two variants of A-TIM: Ma18 (2.7 Å resolution) and Ma21 (1.55 Å resolution). A-TIM is a monomeric loop-deletion variant of triosephosphate isomerase (TIM) which has lost the TIM catalytic properties. Ma18 and Ma21 were identified after extensive directed-evolution selection experiments using an Escherichia coli L-arabinose isomerase knockout strain expressing a randomly mutated A-TIM gene. These variants facilitate better growth of the Escherichia coli selection strain in medium supplemented with 40 mM L-arabinose. Ma18 and Ma21 differ from A-TIM by four and one point mutations, respectively. Ma18 and Ma21 are more stable proteins than A-TIM, as judged from CD melting experiments. Like A-TIM, both proteins are monomeric in solution. In the Ma18 crystal structure loop 6 is open and in the Ma21 crystal structure loop 6 is closed, being stabilized by a bound glycolate molecule. The crystal structures show only small differences in the active site compared with A-TIM. In the case of Ma21 it is observed that the point mutation (Q65L) contributes to small structural rearrangements near Asn11 of loop 1, which correlate with different ligand-binding properties such as a loss of citrate binding in the active site. The Ma21 structure also shows that its Leu65 side chain is involved in van der Waals interactions with neighbouring hydrophobic side-chain moieties, correlating with its increased stability. The experimental data suggest that the increased stability and solubility properties of Ma21 and Ma18 compared with A-TIM cause better growth of the selection strain when coexpressing Ma21 and Ma18 instead of A-TIM.

Fast, cheap and out of control--Insights into thermodynamic and informatic constraints on natural protein sequences from de novo protein design

The accumulated results of thirty years of rational and computational de novo protein design have taught us important lessons about the stability, information content, and evolution of natural proteins. First, de novo protein design has complicated the assertion that biological function is equivalent to biological structure - demonstrating the capacity to abstract active sites from natural contexts and paste them into non-native topologies without loss of function. The structure-function relationship has thus been revealed to be either a generality or strictly true only in a local sense. Second, the simplification to "maquette" topologies carried out by rational protein design also has demonstrated that even sophisticated functions such as conformational switching, cooperative ligand binding, and light-activated electron transfer can be achieved with low-information design approaches. This is because for simple topologies the functional footprint in sequence space is enormous and easily exceeds the number of structures which could have possibly existed in the history of life on Earth. Finally, the pervasiveness of extraordinary stability in designed proteins challenges accepted models for the "marginal stability" of natural proteins, suggesting that there must be a selection pressure against highly stable proteins. This can be explained using recent theories which relate non-equilibrium thermodynamics and self-replication. This article is part of a Special Issue entitled Biodesign for Bioenergetics--The design and engineering of electronc transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Structure-based directed evolution of a monomeric triosephosphate isomerase: toward a pentose sugar isomerase

Through structure-based and directed evolution approaches, a new catalytic activity has been established on the (β/α)8 barrel enzyme triosephosphate isomerase (TIM). This work started from ml8bTIM, a monomeric variant of TIM, in which the phosphate-binding loop (loop-8) had been shortened. Structure analysis suggested an additional point mutation (V233A), converting ml8bTIM into A-TIM. A-TIM has no detectable TIM activity, but it binds the TIM transition state analog, 2-phosphoglycollate. In an in vivo selection approach, we aimed at transferring the activity of three sugar isomerases (L-arabinose isomerase (L-AI), D-xylose isomerase A (D-XI) and D-ribose-5-phosphate isomerase (D-RPI)) onto A-TIM. Escherichia coli knockout variants were constructed, lacking E. coli L-AI, D-XI and D-RPI activities, respectively. Through a systematic approach, new A-TIM variants were obtained only from selection experiments with the L-AI knockout strain. Selection for D-RPI activity was impossible because of an impaired strain due to the gene knockouts. The selection for D-XI activity was unsuccessful, showing the importance of the starting protein for obtaining new biocatalytic properties. The L-AI-directed evolution experiments show that A-TIM already has residual in vivo L-AI activity. Most of the mutations providing A-TIM with enhanced L-AI activity are located in the loops between β-strands and the subsequent α-helices.

Redesigning the type II' β-turn in green fluorescent protein to type I': implications for folding kinetics and stability

Both Type I' and Type II' β-turns have the same sense of the β-turn twist that is compatible with the β-sheet twist. They occur predominantly in two residue β-hairpins, but the occurrence of Type I' β-turns is two times higher than Type II' β-turns. This suggests that Type I' β-turns may be more stable than Type II' β-turns, and Type I' β-turn sequence and structure can be more favorable for protein folding than Type II' β-turns. Here, we redesigned the native Type II' β-turn in GFP to Type I' β-turn, and investigated its effect on protein folding and stability. The Type I' β-turns were designed based on the statistical analysis of residues in natural Type I' β-turns. The substitution of the native "GD" sequence of i+1 and i+2 residues with Type I' preferred "(N/D)G" sequence motif increased the folding rate by 50% and slightly improved the thermodynamic stability. Despite the enhancement of in vitro refolding kinetics and stability of the redesigned mutants, they showed poor soluble expression level compared to wild type. To overcome this problem, i and i + 3 residues of the designed Type I' β-turn were further engineered. The mutation of Thr to Lys at i + 3 could restore the in vivo soluble expression of the Type I' mutant. This study indicates that Type II' β-turns in natural β-hairpins can be further optimized by converting the sequence to Type I'.

Investigating the role of site specific synonymous variation in disease association studies

Synonymous codon changes may not always be neutral indicating their significance in disease association studies, which is almost always overlooked. Synonymous substitutions may affect protein-folding rates leading to protein misfolding and aggregation. Genome wide analysis of 2301 mitochondrial genomes is performed to evaluate the significance of synonymous codons in disease association studies. The analysis revealed usage of rare codons at several sites in mitochondrial genes with rare codon usage higher for hydrophobic amino acids. The analysis suggests that variation data in association studies should be analyzed using site-specific codon usage values to infer the potential phenotypic impact of synonymous changes.

Importance of hydrophobic cavities in allosteric regulation of formylglycinamide synthetase: insight from xenon trapping and statistical coupling analysis

Formylglycinamide ribonucleotide amidotransferase (FGAR-AT) is a 140 kDa bi-functional enzyme involved in a coupled reaction, where the glutaminase active site produces ammonia that is subsequently utilized to convert FGAR to its corresponding amidine in an ATP assisted fashion. The structure of FGAR-AT has been previously determined in an inactive state and the mechanism of activation remains largely unknown. In the current study, hydrophobic cavities were used as markers to identify regions involved in domain movements that facilitate catalytic coupling and subsequent activation of the enzyme. Three internal hydrophobic cavities were located by xenon trapping experiments on FGAR-AT crystals and further, these cavities were perturbed via site-directed mutagenesis. Biophysical characterization of the mutants demonstrated that two of these three voids are crucial for stability and function of the protein, although being ∼20 Å from the active centers. Interestingly, correlation analysis corroborated the experimental findings, and revealed that amino acids lining the functionally important cavities form correlated sets (co-evolving residues) that connect these regions to the amidotransferase active center. It was further proposed that the first cavity is transient and allows for breathing motion to occur and thereby serves as an allosteric hotspot. In contrast, the third cavity which lacks correlated residues was found to be highly plastic and accommodated steric congestion by local adjustment of the structure without affecting either stability or activity.

Structure of a simplified β-hairpin and its ATP complex

The capacity of three designed duodecamer peptides with the low diversity sequence: H1ϕ2I3K4I5D6G7K8ϕ9I10K11H12 where ϕ is His, Phe or Trp, to adopt a β-hairpin conformation was studied using NMR spectroscopy. Whereas KIAβH, the variant with His at positions two and nine, is disordered, KIAβF, the peptide with Phe at these positions, adopts a small population of β-hairpin. A high population of β-hairpin structure was detected for KIAβW, the variant with Trp. Utilizing NMR data, the structure of KIAβW was solved and it reveals a β-hairpin stabilized by hydrophobic interactions between Ile residues on one face and Trp-Trp and cation-π interactions on the opposite face. Upon adding ATP, these peptides show chemical shift changes indicative of ATP binding. The binding of ATP to KIAβW shows a KD ≈ 20 μM at pH 5, 5 °C and has a 1:1 stoichiometry. The KIAβW-ATP complex was determined using NMR spectroscopy and reveals the adenine ring sandwiched between the two Trp indole rings and that ATP binding induces important conformational changes in His1, Trp2, Lys4, Trp9 and Lys11 in the β-hairpin. The implications of these results for the hypothetic presence of β-hairpins and amyloids alongside RNAs on the prebiotic Earth are discussed.

Clusters of branched aliphatic side chains serve as cores of stability in the native state of the HisF TIM barrel protein

Imidazole-3-glycerol phosphate synthase is a heterodimeric allosteric enzyme that catalyzes consecutive reactions in imidazole biosynthesis through its HisF and HisH subunits. The unusually slow unfolding reaction of the isolated HisF TIM barrel domain from the thermophilic bacteria, Thermotoga maritima, enabled an NMR-based site-specific analysis of the main-chain hydrogen bonds that stabilize its native conformation. Very strong protection against exchange with solvent deuterium in the native state was found in a subset of buried positions in α-helices and pervasively in the underlying β-strands associated with a pair of large clusters of isoleucine, leucine and valine (ILV) side chains located in the α7(βα)8(βα)1-2 and α2(βα)3-6β7 segments of the (βα)8 barrel. The most densely packed region of the large cluster, α3(βα)4-6β7, correlates closely with the core of stability previously observed in computational, protein engineering and NMR dynamics studies, demonstrating a key role for this cluster in determining the thermodynamic and structural properties of the native state of HisF. When considered with the results of previous studies where ILV clusters were found to stabilize the hydrogen-bonded networks in folding intermediates for other TIM barrel proteins, it appears that clusters of branched aliphatic side chains can serve as cores of stability across the entire folding reaction coordinate of one of the most common motifs in biology.

Insights into the fold organization of TIM barrel from interaction energy based structure networks

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or “sequence conservation” as the basis for their understanding. Recently “interaction energy” based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the “interaction conservation” viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Proteins are polymers of amino-acids that fold into unique three-dimensional structures to perform cellular functions. This structure formation has been shown to depend on the amino-acid sequences. But examples of proteins with diverse sequences retaining a similar structural fold are quite substantial that we can no longer consider such phenomenon as exceptions. Therefore, this non-canonical relationship has been studied extensively mostly by studying the remote sequence similarities between proteins. Here we have attempted to address the above-mentioned problem by analyzing the similarities in the spatial interactions among amino-acids. Since the protein structure is a resultant of different interactions, we have considered the proteins as networks of interacting amino-acids to derive the common interactions within a popular structural fold called the TIM barrel fold. We were able to find common interactions among different families of the TIM fold and generalize the patterns of interactions by which the fold is being maintained despite sequence diversity. The results substantiate our hypothesis that interaction conservation might by a driving factor in fold formation and this new outlook can be used extensively in engineering proteins with better biophysical characteristics.

Stabilizing proteins from sequence statistics: the interplay of conservation and correlation in triosephosphate isomerase stability

Understanding the determinants of protein stability remains one of protein science's greatest challenges. There are still no computational solutions that calculate the stability effects of even point mutations with sufficient reliability for practical use. Amino acid substitutions rarely increase the stability of native proteins; hence, large libraries and high-throughput screens or selections are needed to stabilize proteins using directed evolution. Consensus mutations have proven effective for increasing stability, but these mutations are successful only about half the time. We set out to understand why some consensus mutations fail to stabilize, and what criteria might be useful to predict stabilization more accurately. Overall, consensus mutations at more conserved positions were more likely to be stabilizing in our model, triosephosphate isomerase (TIM) from Saccharomyces cerevisiae. However, positions coupled to other sites were more likely not to stabilize upon mutation. Destabilizing mutations could be removed both by removing sites with high statistical correlations to other positions and by removing nearly invariant positions at which "hidden correlations" can occur. Application of these rules resulted in identification of stabilizing mutations in 9 out of 10 positions, and amalgamation of all predicted stabilizing positions resulted in the most stable yeast TIM variant we produced (+8 °C). In contrast, a multimutant with 14 mutations each found to stabilize TIM independently was destabilized by 2 °C. Our results are a practical extension to the consensus concept of protein stabilization, and they further suggest the importance of positional independence in the mechanism of consensus stabilization.

Computational design of intermolecular stability and specificity in protein self-assembly

The ability to engineer novel proteins using the principles of molecular structure and energetics is a stringent test of our basic understanding of how proteins fold and maintain structure. The design of protein self-assembly has the potential to impact many fields of biology from molecular recognition to cell signaling to biomaterials. Most progress in computational design of protein self-assembly has focused on α-helical systems, exploring ways to concurrently optimize the stability and specificity of a target state. Applying these methods to collagen self-assembly is very challenging, due to fundamental differences in folding and structure of α- versus triple-helices. Here, we explore various computational methods for designing stable and specific oligomeric systems, with a focus on α-helix and collagen self-assembly.

Directed evolution of a thermophilic endoglucanase (Cel5A) into highly active Cel5A variants with an expanded temperature profile

Cel5A is a highly active endoglucanase from Thermoanaerobacter tengcongensis MB4, displaying an optimal temperature range between 75 and 80°C. After three rounds of error-prone PCR and screening of 4700 mutants, five variants of Cel5A with improved activities were identified by Congo Red based screening method. Compared with the wild type, the best variants 3F6 and C3-13 display 135±6% and 193±8% of the wild type specific activity for the substrate carboxymethyl cellulose (CMC), besides improvements in the relative expression level in Escherichia coli system. Remarkable are especially the improvements in activities at reduced temperatures (50% of maximum activity at 50°C and about 45°C respectively, while 65°C for the wild type). Molecular Dynamics simulations performed on the 3F6 and C3-13 variants show a decreased number of intra-Cel5A hydrogen bonds compared to the wild type, implying a more flexible protein skeleton which correlates well to the higher catalytic activity at lower temperatures. To investigate functions of each individual amino acid position site-directed (saturation) mutagenesis were generated and screened. Amino acid positions Val249 and Ile321 were found to be crucial for improving activity and residue Ile13 (encoded by rare codon AUA) yields an improved expression level in E. coli.

Comparison of the frequency of functional SH3 domains with different limited sets of amino acids using mRNA display

Although modern proteins consist of 20 different amino acids, it has been proposed that primordial proteins consisted of a small set of amino acids, and additional amino acids have gradually been recruited into the genetic code. This hypothesis has recently been supported by comparative genome sequence analysis, but no direct experimental approach has been reported. Here, we utilized a novel experimental approach to test a hypothesis that native-like globular proteins might be easily simplified by a set of putative primitive amino acids with retention of its structure and function than by a set of putative new amino acids. We performed in vitro selection of a functional SH3 domain as a model from partially randomized libraries with different sets of amino acids using mRNA display. Consequently, a library rich in putative primitive amino acids included a larger number of functional SH3 sequences than a library rich in putative new amino acids. Further, the functional SH3 sequences were enriched from the primitive library slightly earlier than from a randomized library with the full set of amino acids, while the function and structure of the selected SH3 proteins with the primitive alphabet were comparable with those from the 20 amino acid alphabet. Application of this approach to various combinations of codons in protein sequences may be useful not only for clarifying the precise order of the amino acid expansion in the early stages of protein evolution but also for efficiently creating novel functional proteins in the laboratory.

The structure of a truncated phosphoribosylanthranilate isomerase suggests a unified model for evolution of the (βα)8 barrel fold

The (βα)(8) barrel is one of the most common protein folds, and enzymes with this architecture display a remarkable range of catalytic activities. Many of these functions are associated with ancient metabolic pathways, and phylogenetic reconstructions suggest that the (βα)(8) barrel was one of the very first protein folds to emerge. Consequently, there is considerable interest in understanding the evolutionary processes that gave rise to this fold. In particular, much attention has been focused on the plausibility of (βα)(8) barrel evolution from homodimers of half barrels. However, we previously isolated a three-quarter-barrel-sized fragment of a (βα)(8) barrel, termed truncated phosphoribosylanthranilate isomerase (trPRAI), that is soluble and almost as thermostable as full-length N-(5'-phosphoribosyl)anthranilate isomerase (PRAI). Here, we report the NMR-derived structure of trPRAI. The subdomain is monomeric, is well ordered and adopts a native-like structure in solution. Side chains from strands β(1) (Glu3 and Lys5), β(2) (Tyr25) and β(6) (Lys122) of trPRAI repack to shield the hydrophobic core from the solvent. This result demonstrates that three-quarter barrels were viable intermediates in the evolution of the (βα)(8) barrel fold. We propose a unified model for (βα)(8) barrel evolution that combines our data, previously published work and plausible scenarios for the emergence of (initially error-prone) genetic systems. In this model, the earliest proto-cells contained diverse pools of part-barrel subdomains. Combinatorial assembly of these subdomains gave rise to many distinct lineages of (βα)(8) barrel proteins, that is, our model excludes the possibility that there was a single (βα)(8) barrel from which all present examples are descended.

Comparative characterization of random-sequence proteins consisting of 5, 12, and 20 kinds of amino acids

Screening of functional proteins from a random-sequence library has been used to evolve novel proteins in the field of evolutionary protein engineering. However, random-sequence proteins consisting of the 20 natural amino acids tend to aggregate, and the occurrence rate of functional proteins in a random-sequence library is low. From the viewpoint of the origin of life, it has been proposed that primordial proteins consisted of a limited set of amino acids that could have been abundantly formed early during chemical evolution. We have previously found that members of a random-sequence protein library constructed with five primitive amino acids show high solubility (Doi et al., Protein Eng Des Sel 2005;18:279-284). Although such a library is expected to be appropriate for finding functional proteins, the functionality may be limited, because they have no positively charged amino acid. Here, we constructed three libraries of 120-amino acid, random-sequence proteins using alphabets of 5, 12, and 20 amino acids by preselection using mRNA display (to eliminate sequences containing stop codons and frameshifts) and characterized and compared the structural properties of random-sequence proteins arbitrarily chosen from these libraries. We found that random-sequence proteins constructed with the 12-member alphabet (including five primitive amino acids and positively charged amino acids) have higher solubility than those constructed with the 20-member alphabet, though other biophysical properties are very similar in the two libraries. Thus, a library of moderate complexity constructed from 12 amino acids may be a more appropriate resource for functional screening than one constructed from 20 amino acids.

Structural aspects of therapeutic enzymes to treat metabolic disorders

Protein therapeutics represents a niche subset of pharmacological agents that is rapidly gaining importance in medicine. In addition to the exceptional specificity that is characteristic of protein therapeutics, several classes of proteins have also been effectively utilized for treatment of conditions that would otherwise lack effective pharmacotherapeutic options. A particularly striking class of protein therapeutics is exogenous enzymes administered for replacement therapy in patients afflicted with metabolic disorders. To date, at least 11 enzymes have either been approved for use, or are in clinical trials for the treatment of selected inherited metabolic disorders. With the recent advancement in structural biology, a significantly larger amount of structural information for several of these enzymes is now available. This article is an overview of the correlation between structural perturbations of these enzymes with the clinical presentation of the respective metabolic conditions, as well as a discussion of the relevant structural modification strategies engaged in improving these enzymes for replacement therapies.

Structural determinants of protein folding

The last several decades have seen an explosion of knowledge in the field of structural biology. With critical advances in spectroscopic techniques in examining structures of biomacromolecules, in maturation of molecular biology techniques, as well as vast improvements in computation prowess, protein structures are now being elucidated at an unprecedented rate. In spite of all the recent advances, the protein folding puzzle remains as one of the fundamental biochemical challenges. A facet to this empiric problem is the structural determinants of protein folding. What are the driving forces that pivot a polypeptide chain to a specific conformation amongst the vast conformation space? In this review, we shall discuss some of the structural determinants to protein folding that have been identified in the recent decades.

Design of thiolate rich metal binding sites within a peptidic framework

A de novo protein design strategy provides a powerful tool to elucidate how heavy metals interact with proteins.Cysteine derivatives of the TRI peptide family (Ac-G(LKALEEK)4G-NH2) have been shown to bind heavy metals in an unusual trigonal geometry. Our present objective was to design binding sites in R-helical scaffolds that are able to form higher coordination number complexes with Cd(II) and Hg(II). Herein, we evaluate the binding of Cd(II) and Hg(II) to double cysteine substituted TRI peptides lacking intervening leucines between sulfurs in the heptads. We compare a -Cysd-X-X-X-Cysa- binding motif found in TRIL12CL16C to the more common -Cysa-X-X-Cysd- sequence of native proteins found in TRIL9CL12C. Compared to TRI, these substitutions destabilize the helical aggregates,leading to mixtures of two- and three-stranded bundles. The three-stranded coiled coils are stabilized by the addition of metals. TRIL9CL12C forms distorted tetrahedral complexes with both Cd(II) and Hg(II), as supported by UV-vis,CD, 113Cd NMR, 199Hg NMR and 111mCd PAC spectroscopy. Additionally, these signatures are very similar to those found for heavy metal substituted rubredoxin. These results suggest that in terms of Hg(II) binding, TRIL9CL12Ccan be considered as a good mimic of the metallochaperone HAH1, that has previously been shown to form protein dimers. TRIL12CL16C has limited ability to generate homoleptic tetrahedral complexes (Cd(SR)42-). These type of complexes were identified only for Hg(II). However, the spectroscopic signatures suggest a different geometry around the metal ion, demonstrating that effective metal sequestration into the hydrophobic interior of the bundle requires more than simply adding two sulfur residues in adjacent layers of the peptide core. Thus, proper design of metal binding sites must also consider the orientation of cysteine sidechains in a vs d positions of the heptads.

Using diastereopeptides to control metal ion coordination in proteins

Here, we report a previously undescribed approach for controlling metal ion coordination geometry in biomolecules by reorientating amino acid side chains through substitution of L- to D-amino acids. These diastereopeptides allow us to manipulate the spatial orientation of amino acid side chains to alter the sterics of metal binding pockets. We have used this approach to design the de novo metallopeptide, Cd(TRIL12L(D)L16C)(3)(-), which is an example of Cd(II) bound to 3 L-Cys as exclusively trigonal CdS(3), as characterized by a combination of (113)Cd NMR and (111m)Cd PAC spectroscopy. We subsequently show that the physical properties of such a site, such as the high pK(a2) for Cd(II) binding of 15.1, is due to the nature of the coordination number and not the ligating group. Further more this approach allowed for the design of a construct, GRANDL12L(D)L16CL26AL30C, capable of independently binding 2 equivalents of Cd(II) to 2 very similar Cys sites as exclusively 3- and 4-, CdS(3) and CdS(3)O, respectively. Demonstrating that we are capable of controlling the Cd(II) coordination number in these 2 sites solely by varying the nature of a noncoordinating second coordination sphere amino acid, with D-leucine and L-alanine resulting in exclusively 3- and 4-coordinate structures, respectively. Cd(II) was found to selectively bind to the 4-coordinate CdS(3)O site, demonstrating that a protein can be designed that displays metal-binding selectivity based solely on coordination number control and not on the chemical identity of coordinating ligands.

Engineering the thermostability of a TIM-barrel enzyme by rational family shuffling

A possible approach to generate enzymes with an engineered temperature optimum is to create chimeras of homologous enzymes with different temperature optima. We tested this approach using two family-10 xylanases from Thermotoga maritima: the thermophilic xylanase A catalytic domain (TmxAcat, T(opt)=68 degrees C), and the hyperthermophilic xylanase B (TmxB, T(opt)=102 degrees C). Twenty-one different chimeric constructs were created by mimicking family shuffling in a rational manner. The measured temperature optima of the 16 enzymatically active chimeras do not monotonically increase with the percentage of residues coming from TmxB. Only four chimeras had a higher temperature optimum than TmxAcat, the most stable variant (T(opt)=80 degrees C) being the one in which both terminal segments came from TmxB. Further analysis suggests that the interaction between the N- and C-terminal segments has a disproportionately high contribution to the overall thermostability. The results may be generalizable to other enzymes where the N- and C-termini are in contact.

A theoretical comparison of self-assembling alpha- and beta-peptide nanostructures: toward design of beta-barrel frameworks

Self-assembling peptide-based nanotubes are among the most investigated bioactive compounds as a result of their numerous potential applications as novel biomaterials. To support rational bottom-up design of such artificial nanosystems, here we investigate structural and energetic properties of various sheet-derived nanotubes. We carried out high level quantum chemical calculations on large models, composed of up to 32 amino acids, and characterized structures from extended beta-sheets to the molecular framework of beta-barrel proteins. Surprisingly, enzyme-resistant nonnatural beta-peptides have an affinity to form nanotubes that is remarkably higher than that of natural alpha-peptides. We analyzed the stability of both systems depending on (i) parallel or antiparallel orientation, (ii) the number of peptide strands, and (iii) the formed hydrogen bond pattern. Applicability is outlined by investigating guest molecules in the tubes. It is hoped that the structural and energetic data presented here will be effectively used in the design of novel peptide nanosystems.

Genetic selection reveals the role of a buried, conserved polar residue

The burial of nonpolar surface area is known to enhance markedly the conformational stability of proteins. The contribution from the burial of polar surface area is less clear. Here, we report on the tolerance to substitution of Ser75 of bovine pancreatic ribonuclease (RNase A), a residue that has the unusual attributes of being buried, conserved, and polar. To identify variants that retain biological function, we used a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity. Cell growth at 30 degrees C, 37 degrees C, and 44 degrees C correlated with residue size, indicating that the primary attribute of Ser75 is its small size. The side-chain hydroxyl group of Ser75 forms a hydrogen bond with a main-chain nitrogen. The conformational stability of the S75A variant, which lacks this hydrogen bond, was diminished by DeltaDeltaG = 2.5 kcal/mol. Threonine, which can reinstate this hydrogen bond, provided a catalytically active RNase A variant at higher temperatures than did some smaller residues (including aspartate), indicating that a secondary attribute of Ser75 is the ability of its uncharged side chain to accept a hydrogen bond. These results provide insight on the imperatives for the conservation of a buried polar residue.

Thermodynamic properties of BPTI variants with highly simplified amino acid sequences

We report the first detailed thermodynamic analysis of simplified proteins by differential scanning calorimetry (DSC). The experiments were carried out with five simplified BPTI variants, whose structures and activities have been reported, in which several residues not essential for specifying the tertiary structure were replaced by alanine. In most aspects, the thermodynamics of simplified proteins were very similar to, if not essentially identical with, those of natural proteins. In particular, they undergo a highly cooperative two-state thermal unfolding process with a large enthalpy change, which is a thermodynamic hallmark of the native state of natural globular proteins. Furthermore, the specific enthalpy and entropy changes upon unfolding at 110 degrees C were close to values invariably observed for small natural globular proteins (55 J g(-1) and ~16 J K(-1) g(-1), respectively). On the other hand, two simplified BPTI variants, BPTI-21 and BPTI-22 (containing 21 and 22 alanine residues), were enthalpically stabilized while entropically destabilized with respect to the reference BPTI-[5,55] molecule. This peculiar type of entropy-enthalpy compensation is in sharp contrast to the usual enthalpy destabilization/entropy stabilization observed in mutational studies of natural proteins. Overall, we conclude that a thermodynamic native state can be achieved by proteins encoded with extensively simplified sequences.

Novel enzymes through design and evolution

The generation of enzymes with new catalytic activities remains a major challenge. So far, several different strategies have been developed to tackle this problem, including site-directed mutagenesis, random mutagenesis (directed evolution), antibody catalysis, computational redesign, and de novo methods. Using these techniques, a broad array of novel enzymes has been created (aldolases, decarboxylases, dehydratases, isomerases, oxidases, reductases, and others), although their low efficiencies (10 to 100 M(-1) s(-l)) compared to those of the best natural enzymes (10(6) to 10(8) M(-1) s(-1)) remains a significant concern. Whereas rational design might be the most promising and versatile approach to generating new activities, directed evolution seems to be the best way to optimize the catalytic properties of novel enzymes. Indeed, impressive successes in enzyme engineering have resulted from a combination of rational and random design.

The stability effects of protein mutations appear to be universally distributed

How the thermodynamic stability effects of protein mutations (DeltaDeltaG) are distributed is a fundamental property related to the architecture, tolerance to mutations (mutational robustness), and evolutionary history of proteins. The stability effects of mutations also dictate the rate and dynamics of protein evolution, with deleterious mutations being the main inhibitory factor. Using the FoldX algorithm that attempts to computationally predict DeltaDeltaG effects of mutations, we deduced the overall distributions of stability effects for all possible mutations in 21 different globular, single domain proteins. We found that these distributions are strikingly similar despite a range of sizes and folds, and largely follow a bi-Gaussian function: The surface residues exhibit a narrow distribution with a mildly destabilizing mean DeltaDeltaG ( approximately 0.6 kcal/mol), whereas the core residues exhibit a wider distribution with a stronger destabilizing mean ( approximately 1.4 kcal/mol). Since smaller proteins have a higher fraction of surface residues, the relative weight of these single distributions correlates with size. We also found that proteins evolved in the laboratory follow an essentially identical distribution, whereas de novo designed folds show markedly less destabilizing distributions (i.e. they seem more robust to the effects of mutations). This bi-Gaussian model provides an analytical description of the predicted distributions of mutational stability effects. It comprises a novel tool for analyzing proteins and protein models, for simulating the effect of mutations under evolutionary processes, and a quantitative description of mutational robustness.

Attempt to simplify the amino-acid sequence of photoactive yellow protein with a set of simple rules

To understand the information encoded in an amino-acid sequence, the authors have attempted to simplify the amino-acid sequence of photoactive yellow protein (PYP) with a set of simple rules. The rules are designed to reduce overlapping structural information. The simplified PYP protein, which was composed of only nine species of amino acids (Ser, Val, Asp, Lys, Phe, Met, Gly, Pro, and Cys), took a completely different structure than the native conformation. Even after the evolutionarily conserved residues were restored in the simplified protein, the PYP variant did not properly fold, indicating that the information encoded in the conserved residues is insufficient for the structure formation. Additional restorations of the substituted hydrophilic or hydrophobic residues did not lead to a variant that formed the native structure. The structural properties of these variants and the wild-type protein in aqueous solution differed. Partial simplification was successfully performed by creating chimeric proteins composed of combinations of wild-type PYP and sPYPIII. The structural characterization of each chimeric protein indicates that the important information on the structure formation is encoded in the beta-scaffold region.

Enzyme-Driven Speciation: Crystallizing Archaea via Lipid Capture

As the origin(s) of life on Earth remains an open question, detailed characteristics about the "last universal ancestor" (LUA) continue to be obscured. Here we provide arguments that strengthen the bacterial-like nature of the LUA. Our view attempts to recreate the evolution of archaeal lipids, the major components of the distinctive membrane that encapsulates these ancient prokaryotes. We show that (S)- 3-O-geranylgeranylglyceryl phosphate synthase (GGGPS), a TIM-barrel protein that performs the committed step in archaeal lipid synthesis, likely evolved from the duplication and fusion of a (betaalpha)4 half-barrel ancestor. By comparison to the well-characterized HisA and HisF TIM-barrel proteins, we propose a time line for the invention of this diagnostic archaeal biosynthetic pathway. After excluding the possibility of horizontal gene transfer, we conclude that the evolutionary history of GGGPS mirrors the emergence of Archaea from the LUA. We illustrate aspects of this "lipid capture" model that support its likelihood in recreating key evolutionary events and, as our hypothesis is built on a single initiating event, we suggest that the appearance of GGGPS represents an example of enzyme-driven speciation.

Natural polypeptide scaffolds: beta-sheets, beta-turns, and beta-hairpins

This paper provides an introduction to fundamental conformational states of polypeptides in the beta-region of phi,psi space, in which the backbone is extended near to its maximal length, and to more complex architectures in which extended segments are linked by turns and loops. There are several variants on these conformations, and they comprise versatile scaffolds for presentation of side chains and backbone amides for molecular recognition and designed catalysts. In addition, the geometry of these fundamental folds can be readily mimicked in peptidomimetics.

Intelligent design: the de novo engineering of proteins with specified functions

One of the principal successes of de novo protein design has been the creation of small, robust protein-cofactor complexes which can serve as simplified models, or maquettes, of more complicated multicofactor protein complexes commonly found in nature. Different maquettes, generated by us and others, recreate a variety of aspects, or functional elements, recognized as parts of natural enzyme function. The current challenge is to both expand the palette of functional elements and combine and/or integrate them in recreating familiar enzyme activities or generating novel catalysis in the simplest protein scaffolds.