Locating the stabilizing residues in (alpha/beta)8 barrel proteins based on hydrophobicity, long-range interactions, and sequence conservation

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.

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.

The rs61742690 (S783N) single nucleotide polymorphism is a suitable target for disrupting BCL11A-mediated foetal-to-adult globin switching

Background

B-cell lymphoma/leukaemia 11A (BCL11A) is a C2H2-type zinc-finger transcription factor protein that is a critical modulator of haemoglobin switching and suppresses the production of foetal haemoglobin. Variation in the BCL11A gene ameliorates the severity of sickle cell disease (SCD) and β-thalassemia (β-thal). The BCL11A gene is located on chromosome 2p16.1 and encodes an 835-amino acid protein.

Method

Using state-of-the-art in silico tools, this study examined the most pathogenic non-synonymous single nucleotide polymorphisms (nsSNPs) that disrupt the BCL11A protein and mediate foetal-to-adult globin switching. A total of 11,463 SNPs were retrieved from the Single Nucleotide Polymorphism database (dbSNP). These included 799 in the 5′ untranslated region (UTR), 486 in the 3′ UTR, and 266 non-synonymous, 189 coding synonymous, six nonsense, and six stop-gained SNPs.

Results and discussion

In silico tools (SIFT, SNAP, PolyPhen-2, PANTHER, I-Mutant, PROVEAN, SNPs&GO, mCSM, and PhD-SNP) predicted the five most-deleterious nsSNPs: rs61742690, rs62142605, rs17028351, rs115666026, and rs74987258. Molecular dynamic simulation and homology modelling of the mutated proteins (S783N, D643N, G451S, K670R, and M313L) of the most deleterious nsSNPs revealed their functional and structural impact. nsSNP rs61742690 was predicted to be the most deleterious, as supported by eight of the nine in silico tools.

Conclusions

Complete failure in the protein–protein interactions with functional partners (KLF1 and others) and significant changes (±100% variation) in the interface energy revealed that rs61742690 (S783N) in the zinc-finger domain is a suitable target for disrupting BCL11A-mediated foetal-to-adult globin switching.

Identification and Analysis of Key Residues in Protein-RNA Complexes

Protein-RNA complexes play important roles in various biological processes. The functions of protein-RNA complexes are dictated by their interactions, binding, stability, and affinity. In this work, we have identified the key residues (KRs), which are involved in both stability and binding. We found that 42 percent of considered proteins share common binding and stabilizing residues, whereas these residues are distinct in 58 percent of the proteins. Overall, 5 percent of stabilizing and 3 percent of binding residues serve as key residues. These residues are enriched with the combination of polar, charged, aliphatic, and aromatic residues. Analysis on subclasses of protein-RNA complexes based on protein structural class, function and RNA type showed that regulatory proteins, and complexes with single stranded RNA and rRNA have appreciable number of key residues. Specifically, Arg, Tyr, and Thr are preferred in most of the subclasses of protein-RNA complexes. In addition, residues with similar chemical behavior have different preferences to be KRs, such that Arg, Tyr, Val, and Thr are preferred over Lys, Trp, Ile, and Ser, respectively. Atomic level contacts revealed that charged and polar-nonpolar contacts are dominant in enzymes, polar in structural, and nonpolar in regulatory proteins. On the other hand, polar-nonpolar contacts are enriched in all these classes of protein-RNA complexes. Further, the influence of sequence and structural features such as conservation score, surrounding hydrophobicity, solvent accessibility, secondary structure, and long-range order in key residues are also discussed. We envisage that the present study provides insights to understand the structural and functional aspects of protein-RNA complexes.

Prediction of structural consequences for disease causing variants in C21orf2 protein using computational approaches

Amyotrophic lateral sclerosis (ALS), a progressive motor-neurone disease, affects individuals usually aged between 50 and 70 years. C21orf2, recently identified as the new ALS susceptibility gene, harbours rare missense mutations that cause this fatal disease. We used bioinformatics and molecular modelling approaches to study specific ALS-associated mutations in C21orf2. Both native and mutant structures of the protein obtained from homology modelling were analysed in detail to gain insights into the potential impact of these mutations on the protein structure and its function. Our analyses reveal that more than 75% of the mutations are likely to be deleterious. These effects seem to carry through to mouse C21orf2 as well, indicating that mouse would make a viable animal model to study this ALS gene in detail.

Correlation of fitness landscapes from three orthologous TIM barrels originates from sequence and structure constraints

Sequence divergence of orthologous proteins enables adaptation to environmental stresses and promotes evolution of novel functions. Limits on evolution imposed by constraints on sequence and structure were explored using a model TIM barrel protein, indole-3-glycerol phosphate synthase (IGPS). Fitness effects of point mutations in three phylogenetically divergent IGPS proteins during adaptation to temperature stress were probed by auxotrophic complementation of yeast with prokaryotic, thermophilic IGPS. Analysis of beneficial mutations pointed to an unexpected, long-range allosteric pathway towards the active site of the protein. Significant correlations between the fitness landscapes of distant orthologues implicate both sequence and structure as primary forces in defining the TIM barrel fitness landscape and suggest that fitness landscapes can be translocated in sequence space. Exploration of fitness landscapes in the context of a protein fold provides a strategy for elucidating the sequence-structure-fitness relationships in other common motifs.

The TIM barrel fold is an evolutionarily conserved motif found in proteins with a variety of enzymatic functions. Here the authors explore the fitness landscape of the TIM barrel protein IGPS and uncover evolutionary constraints on both sequence and structure, accompanied by long range allosteric interactions.

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.

A multifaceted computational report on the variants effect on KIR2DL3 and IFNL3 candidate gene in HCV clearance

HCV infection causes acute and chronic liver diseases including, cirrhosis and hepatocellular carcinoma. Following HCV infection, spontaneous clearance occurs in approximately 20 % of the population dependant upon HCV genotype. In this study, functional and non-functional variant analysis was executed for the classical and the latest HCV clearance candidate genes namely, KIR2DL3 and IFNL3. Initially, the functional effects of non-synonymous SNPs were assigned on exposing to homology based tools, SIFT, PolyPhen-2 and PROVEAN. Further, UTR and splice sites variants were scanned for the gene expression and regulation changes. Subsequently, the haplotype and CNV were also identified. The mutation H77Y of KIR2DL3 and R157Q, H156Y, S63L, R157W, F179V, H128R, T101M, R180C, and F176I of IFNL3 results in conservation, RMSD, total energy, stability, and secondary structures revealed a negative impact on the structural fitness. UTRscan and the splice site result indicate functional change, which may affect gene regulation and expression. The graphical display of selected population shows alleles like rs270779, rs2296370, rs10423751, rs12982559, rs9797797, and rs35987710 of KIR2DL3 and rs12972991, rs12980275, rs4803217, rs8109886, and rs8099917 of IFNL3 are in high LD with a measure of [Formula: see text] broadcasting its protective effect in HCV clearance. Similarly, CNV report suggests major DNA fragment loss that could have a profound impact on the gene expression affecting the overall phenotype. This roundup report specifies the effect of NK cell receptor, KIR2DL3 and IFNL3 variants that can have a better prospect in GWAS and immunogenetic studies leading to better understanding of HCV clearance and progression.

In silico transcriptional regulation and functional analysis of dengue shock syndrome associated SNPs in PLCE1 and MICB genes

Single nucleotide polymorphisms (SNPs) in PLCE1 and MICB genes increase risk for the development of dengue shock syndrome (DSS). We used Bioinformatics tools to predict alterations at the transcriptional and posttranslational levels driven by PLCE1 and MICB SNPs associated with DSS. Functional and phenotypic analysis conducted to determine deleterious SNPs and impact of amino acid substitution on the structure and function of proteins identified rs2274223 (H1619R) as deleterious to protein coding as it induces structural change in the C2 domain of PLCε, with the mutant residue more positively charged than the wild-type residue (RMSD score, 1.75 Å). Moreover, rs2274223 condenses the chromatin-repressing PLCε expression in DSS. Briefly, this study presents the impact of a single nucleotide transition at SNPs associated with DSS on differential protein binding patterns with PLCE1 and MICB genes and on protein structure modification and their possible role in the pathogenesis of DSS.

Elucidating the Key Determinants of Structure, Folding, and Stability for the ( 4β+ α ) Conformation of the B1 Domain of Protein G Using Bioinformatics Approaches

The B1 domain of protein G (GB1) is a small, 56 amino acid bacterial immunoglobulin-binding protein with a 4β+ α fold. Architecturally, it is composed of a two-layer sandwich consisting of a four-stranded β -sheet that packs against an α -helix. Using several bioinformatics approaches, we investigated which residues may be key determinants of this fold. We identified nine structurally conserved amino acids using a conservation analysis and propose they are critical to forming and stabilizing the fold. The nine conserved residues form a predominantly hydrophobic nucleus within the core of GB1. A network analysis of all the long-range interactions in the structure of GB1 in concert with a betweenness centrality analysis revealed the relative significance of each conserved amino acid residue based on the number and location of the interactions. This bioinformatics analysis provides an important foundation for the design and interpretation of both computational and experimental work which may be helpful in solving the protein folding problem.

Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 β-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the β-barrel core is highly rigid, while the α-helices that flank the β-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.

In-Silico Computing of the Most Deleterious nsSNPs in HBA1 Gene

Background

α-Thalassemia (α-thal) is a genetic disorder caused by the substitution of single amino acid or large deletions in the HBA1 and/or HBA2 genes.

Method

Using modern bioinformatics tools as a systematic in-silico approach to predict the deleterious SNPs in the HBA1 gene and its significant pathogenic impact on the functions and structure of HBA1 protein was predicted.

Results and Discussion

A total of 389 SNPs in HBA1 were retrieved from dbSNP database, which includes: 201 non-coding synonymous (nsSNPs), 43 human active SNPs, 16 intronic SNPs, 11 mRNA 3′ UTR SNPs, 9 coding synonymous SNPs, 9 5′ UTR SNPs and other types. Structural homology-based method (PolyPhen) and sequence homology-based tool (SIFT), SNPs&Go, PROVEAN and PANTHER revealed that 2.4% of the nsSNPs are pathogenic.

Conclusions

A total of 5 nsSNPs (G60V, K17M, K17T, L92F and W15R) were predicted to be responsible for the structural and functional modifications of HBA1 protein. It is evident from the deep comprehensive in-silico analysis that, two nsSNPs such as G60Vand W15R in HBA1 are highly deleterious. These “2 pathogenic nsSNPs” can be considered for wet-lab confirmatory analysis.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

Modulation of the thermostability and substrate specificity of Candida rugosa lipase1 by altering the acyl-binding residue Gly414 at the α-helix-connecting bend

Candida rugosa Lipase1 (LIP1) is widely used in industrial applications. Optimizing its catalytic performance is still a challenging goal for protein engineers. Mutagenesis of key residues in the active site of the enzyme may provide an effective strategy for enhancing stability and altering substrate specificity. In this study, multiple sequence alignment and structural analysis revealed that the acyl-binding residue, Gly414, of LIP1, which is located at a bend connecting α-helixes, was the non-conserved residue in five other isoenzymes. Using saturation mutagenesis, four mutants with improved stability (G414A, G414M, G414H and G414W) were obtained. Compared to the wild type, the best mutant (G414W) exhibited a remarkable 6.5-fold enhancement in half-life at 60 °C and a 14 °C higher T50(15). Its optimum temperature was increased by 15 °C. Simultaneously, G414W displayed a shift in substrate preference from medium-chain to short-chain pNP-ester. Modeling analysis showed that the multiple interactions formed by hydrophobic clusters and hydrogen bonds in the acyl-binding tunnel might lead to the observed thermostability improvement. Additionally, the bulky tryptophan substitution formed a strong steric hindrance to the accommodation of long-chain substrates in the tunnel. These results indicate that the key acyl-binding residues at the α-helix-connecting bend could mediate enzyme stability and catalytic substrate spectra.

Structure-based network analysis of activation mechanisms in the ErbB family of receptor tyrosine kinases: the regulatory spine residues are global mediators of structural stability and allosteric interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.

Radical SAM enzyme QueE defines a new minimal core fold and metal-dependent mechanism

7-carboxy-7-deazaguanine synthase (QueE) catalyzes a key S-adenosyl-L-methionine (AdoMet)- and Mg(2+)-dependent radical-mediated ring contraction step, which is common to the biosynthetic pathways of all deazapurine-containing compounds. QueE is a member of the AdoMet radical superfamily, which employs the 5'-deoxyadenosyl radical from reductive cleavage of AdoMet to initiate chemistry. To provide a mechanistic rationale for this elaborate transformation, we present the crystal structure of a QueE along with structures of pre- and post-turnover states. We find that substrate binds perpendicular to the [4Fe-4S]-bound AdoMet, exposing its C6 hydrogen atom for abstraction and generating the binding site for Mg(2+), which coordinates directly to the substrate. The Burkholderia multivorans structure reported here varies from all other previously characterized members of the AdoMet radical superfamily in that it contains a hypermodified (β6/α3) protein core and an expanded cluster-binding motif, CX14CX2C.

Computational analysis of C-reactive protein for assessment of molecular dynamics and interaction properties

Serum C-reactive protein (CRP) is used as a marker of inflammation in several diseases including autoimmune disease and cardiovascular disease. CRP, a member of the pentraxin family, is comprised of five identical subunits. CRP has diverse ligand-binding properties which depend upon different structural states of CRP. However, little is known about the molecular dynamics and interaction properties of CRP. In this study, we used SAPS, SCRATCH protein predictor, PDBsum, ConSurf, ProtScale, Drawhca, ASAView, SCide and SRide server and performed comprehensive analyses of molecular dynamics, protein-protein and residue-residue interactions of CRP. We used 1GNH.pdb file for the crystal structure of human CRP which generated two pentamers (ABCDE and FGHIJ). The number of residues involved in residue-residue interactions between A-B, B-C, C-D, D-E, F-G, G-H, H-I, I-J, A-E and F-J subunits were 12, 11, 10, 11, 12, 11, 10, 11, 10 and 10, respectively. Fifteen antiparallel β sheets were involved in β-sheet topology, and five β hairpins were involved in forming the secondary structure. Analysis of hydrophobic segment distribution revealed deviations in surface hydrophobicity at different cavities present in CRP. Approximately 33 % of all residues were involved in the stabilization centers. We show that the bioinformatics tools can provide a rapid method to predict molecular dynamics and interaction properties of CRP. Our prediction of molecular dynamics and interaction properties of CRP combined with the modeling data based on the known 3D structure of CRP is helpful in designing stable forms of CRP mutants for structure-function studies of CRP and may facilitate in silico drug design for therapeutic targeting of CRP.

Octarellin VI: using rosetta to design a putative artificial (β/α)8 protein

The computational protein design protocol Rosetta has been applied successfully to a wide variety of protein engineering problems. Here the aim was to test its ability to design de novo a protein adopting the TIM-barrel fold, whose formation requires about twice as many residues as in the largest proteins successfully designed de novo to date. The designed protein, Octarellin VI, contains 216 residues. Its amino acid composition is similar to that of natural TIM-barrel proteins. When produced and purified, it showed a far-UV circular dichroism spectrum characteristic of folded proteins, with α-helical and β-sheet secondary structure. Its stable tertiary structure was confirmed by both tryptophan fluorescence and circular dichroism in the near UV. It proved heat stable up to 70°C. Dynamic light scattering experiments revealed a unique population of particles averaging 4 nm in diameter, in good agreement with our model. Although these data suggest the successful creation of an artificial α/β protein of more than 200 amino acids, Octarellin VI shows an apparent noncooperative chemical unfolding and low solubility.

Investigations on the role of π-π interactions and π-π networks in eNOS and nNOS proteins

π-π Interactions play an important role in the stability of protein structures. In the present study, we have analyzed the influence of π-π interactions in eNOS and nNOS proteins. The contribution of these π-π interacting residues in sequential separation, secondary structure involvement, solvent accessibility and stabilization centers has been evaluated. π-π interactions stabilize the core regions within eNOS and nNOS proteins. π-π interacting residues are evolutionary conserved. There is a significant number of π-π interactions in spite of the lesser natural occurrences of π-residues in eNOS and nNOS proteins. In addition to π-π interactions, π residues also form π-π networks in both eNOS and nNOS proteins which might play an important role in the structural stability of these protein structures.

Influence of C-H...O interactions on the structural stability of β-lactamases

β-Lactamases produced by pathogenic bacteria cleave β-lactam antibiotics and render them ineffective. Understanding the principles that govern the structural stability of β-lactamases requires elucidation of the nature of the interactions that are involved in stabilization. In the present study, we systematically analyze the influence of CH...O interactions on determining the specificity and stability of β-lactamases in relation to environmental preferences. It is interesting to note that all the residues located in the active site of β-lactamases are involved in CH...O interactions. A significant percentage of CH...O interactions have a higher conservation score and short-range interactions are the predominant type of interactions in β-lactamases. These results will be useful in understanding the stability patterns of β-lactamases.

Emerging role of N- and C-terminal interactions in stabilizing (β/α)8 fold with special emphasis on Family 10 xylanases

Xylanases belong to an important class of industrial enzymes. Various xylanases have been purified and characterized from a plethora of organisms including bacteria, marine algae, plants, protozoans, insects, snails and crustaceans. Depending on the source, the enzymatic activity of xylanases varies considerably under various physico-chemical conditions such as temperature, pH, high salt and in the presence of proteases. Family 10 or glycosyl hydrolase 10 (GH10) xylanases are one of the well characterized and thoroughly studied classes of industrial enzymes. The TIM-barrel fold structure which is ubiquitous in nature is one of the characteristics of family 10 xylanases. Family 10 xylanases have been used as a “model system” due to their TIM-barrel fold to dissect and understand protein stability under various conditions. A better understanding of structure-stability-function relationships of family 10 xylanases allows one to apply these governing molecular rules to engineer other TIM-barrel fold proteins to improve their stability and retain function(s) under adverse conditions. In this review, we discuss the implications of N-and C-terminal interactions, observed in family 10 xylanases on protein stability under extreme conditions. The role of metal binding and aromatic clusters in protein stability is also discussed. Studying and understanding family 10 xylanase structure and function, can contribute to our protein engineering knowledge.

Computer-Aided Protein Directed Evolution: a Review of Web Servers, Databases and other Computational Tools for Protein Engineering

The combination of computational and directed evolution methods has proven a winning strategy for protein engineering. We refer to this approach as computer-aided protein directed evolution (CAPDE) and the review summarizes the recent developments in this rapidly growing field. We will restrict ourselves to overview the availability, usability and limitations of web servers, databases and other computational tools proposed in the last five years. The goal of this review is to provide concise information about currently available computational resources to assist the design of directed evolution based protein engineering experiment.

Improving the thermostability of alpha-amylase by combinatorial coevolving-site saturation mutagenesis

Background

The generation of focused mutant libraries at hotspot residues is an important strategy in directed protein evolution. Existing methods, such as combinatorial active site testing and residual coupling analysis, depend primarily on the evolutionary conserved information to find the hotspot residues. Hardly any attention has been paid to another important functional and structural determinants, the functionally correlated variation information--coevolution.

Results

In this paper, we suggest a new method, named combinatorial coevolving-site saturation mutagenesis (CCSM), in which the functionally correlated variation sites of proteins are chosen as the hotspot sites to construct focused mutant libraries. The CCSM approach was used to improve the thermal stability of α-amylase from Bacillus subtilis CN7 (Amy7C). The results indicate that the CCSM can identify novel beneficial mutation sites, and enhance the thermal stability of wild-type Amy7C by 8°C ( T5030), which could not be achieved with the ordinarily rational introduction of single or a double point mutation.

Conclusions

Our method is able to produce more thermostable mutant α-amylases with novel beneficial mutations at new sites. It is also verified that the coevolving sites can be used as the hotspots to construct focused mutant libraries in protein engineering. This study throws new light on the active researches of the molecular coevolution.

A compact review on the comparison of conventional and non-conventional interactions on the structural stability of therapeutic proteins

Therapeutic proteins carry out the most difficult tasks in living cells. They do so by interacting specifically with other molecules. This requires that they fold to a unique and more stable conformation. A prerequisite for comprehending the folding processes in their immense complexity entails a thorough understanding of many weak interactions. The purpose of this review is to systematically study the role of weak interactions such as cation-π, C-H......π, N-H......π and O-H......π, in the set of 49 therapeutic proteins. The importance of many of these interactions (for example, cationic residues interacting with π system) is revealed by the higher degree of conservation observed for them in protein structures. These interactions are mainly formed by long-range contacts and significant percentage of cation-π, C-H......π, N-H......π and O-H......π interacting residues had one or more stabilization centers. Further, a comparison of conventional and nonconventional interactions in the present data set unambiguously highlights the significance of these weak interactions in the structural stability of therapeutic proteins. We propose that the incorporation of the entirety of these interactions leads to a more complete description of the problem, and that this could provide new perspectives and new possible answers.

Salt bridges: geometrically specific, designable interactions

Salt bridges occur frequently in proteins, providing conformational specificity and contributing to molecular recognition and catalysis. We present a comprehensive analysis of these interactions in protein structures by surveying a large database of protein structures. Salt bridges between Asp or Glu and His, Arg, or Lys display extremely well-defined geometric preferences. Several previously observed preferences are confirmed, and others that were previously unrecognized are discovered. Salt bridges are explored for their preferences for different separations in sequence and in space, geometric preferences within proteins and at protein-protein interfaces, co-operativity in networked salt bridges, inclusion within metal-binding sites, preference for acidic electrons, apparent conformational side chain entropy reduction on formation, and degree of burial. Salt bridges occur far more frequently between residues at close than distant sequence separations, but, at close distances, there remain strong preferences for salt bridges at specific separations. Specific types of complex salt bridges, involving three or more members, are also discovered. As we observe a strong relationship between the propensity to form a salt bridge and the placement of salt-bridging residues in protein sequences, we discuss the role that salt bridges might play in kinetically influencing protein folding and thermodynamically stabilizing the native conformation. We also develop a quantitative method to select appropriate crystal structure resolution and B-factor cutoffs. Detailed knowledge of these geometric and sequence dependences should aid de novo design and prediction algorithms.

In silico characterization of thermostable lipases

Thermostable lipases are of high priority for industrial applications as they are endowed with the capability of carrying out diversified reactions at elevated temperatures. Extremophiles are their potential source. Sequence and structure annotation of thermostable lipases can elucidate evolution of lipases from their mesophilic counterparts with enhanced thermostability hence better industrial potential. Sequence analysis highlighted the conserved residues in bacterial and fungal thermostable lipases. Higher frequency of AXXXA motif and poly Ala residues in lid domain of thermostable Bacillus lipases were distinguishing characteristics. Comparison of amino acid composition among thermostable and mesostable lipases brought into light the role of neutral, charged and aromatic amino acid residues in enhancement of thermostability. Structural annotation of thermostable lipases with that of mesostable lipases revealed some striking features which are increment of gamma turns in thermostable lipases; being first time reported in our paper, longer beta strands, lesser beta-branched residues in helices, increase in charged-neutral hydrogen bonding pair, hydrophobic-hydrophobic contact and differences in the N-cap and C-cap residues of the α helices. Conclusively, it can be stated that subtle changes in the arrangement of amino acid residues in the tertiary structure of lipases contributes to enhanced thermostability.

Nature of cation-pi interactions and their role in structural stability of immunoglobulin proteins

Cation-pi interactions are known to be important contributors to protein stability and ligand-protein interactions. In this study, we have analyzed the influence of cation-pi interactions in single chain immunoglobulin proteins. We observed 87 cation-pi interactions in a data set of 33 proteins. These interactions are mainly formed by long-range contacts, and there is preference of Arg over Lys in these interactions. Arg-Tyr interactions are predominant among the various pairs analyzed. Despite the scarcity of interactions involving Trp, the average energy for Trp-cation interactions is quite high. This information suggests that the cation-pi interactions involving Trp might be of high relevance to the proteins. Secondary structure analysis reveals that cation-pi interactions are formed preferably between residues in which at least one is in beta-strand. Proteins having beta-strand regions have the highest number of cation-pi interaction-forming residues.

The critical role of N- and C-terminal contact in protein stability and folding of a family 10 xylanase under extreme conditions

Background

Stabilization strategies adopted by proteins under extreme conditions are very complex and involve various kinds of interactions. Recent studies have shown that a large proportion of proteins have their N- and C-terminal elements in close contact and suggested they play a role in protein folding and stability. However, the biological significance of this contact remains elusive.

Methodology

In the present study, we investigate the role of N- and C-terminal residue interaction using a family 10 xylanase (BSX) with a TIM-barrel structure that shows stability under high temperature, alkali pH, and protease and SDS treatment. Based on crystal structure, an aromatic cluster was identified that involves Phe4, Trp6 and Tyr343 holding the N- and C-terminus together; this is a unique and important feature of this protein that might be crucial for folding and stability under poly-extreme conditions.

Conclusion

A series of mutants was created to disrupt this aromatic cluster formation and study the loss of stability and function under given conditions. While the deletions of Phe4 resulted in loss of stability, removal of Trp6 and Tyr343 affected in vivo folding and activity. Alanine substitution with Phe4, Trp6 and Tyr343 drastically decreased stability under all parameters studied. Importantly, substitution of Phe4 with Trp increased stability in SDS treatment. Mass spectrometry results of limited proteolysis further demonstrated that the Arg344 residue is highly susceptible to trypsin digestion in sensitive mutants such as ΔF4, W6A and Y343A, suggesting again that disruption of the Phe4-Trp6-Tyr343 (F-W-Y) cluster destabilizes the N- and C-terminal interaction. Our results underscore the importance of N- and C-terminal contact through aromatic interactions in protein folding and stability under extreme conditions, and these results may be useful to improve the stability of other proteins under suboptimal conditions.

Thermodynamic database for proteins: features and applications

We have developed a thermodynamic database for proteins and mutants, ProTherm, which is a collection of a large number of thermodynamic data on protein stability along with the sequence and structure information, experimental methods and conditions, and literature information. This is a valuable resource for understanding/predicting the stability of proteins, and it can be accessible at http://www.gibk26.bse.kyutech.ac.jp/jouhou/Protherm/protherm.html . ProTherm has several features including various search, display, and sorting options and visualization tools. We have analyzed the data in ProTherm to examine the relationship among thermodynamics, structure, and function of proteins. We describe the progress on the development of methods for understanding/predicting protein stability, such as (i) relationship between the stability of protein mutants and amino acid properties, (ii) average assignment method, (iii) empirical energy functions, (iv) torsion, distance, and contact potentials, and (v) machine learning techniques. The list of online resources for predicting protein stability has also been provided.

Influence of cation-pi interactions to the structural stability of prokaryotic and eukaryotic translation elongation factors

We have investigated the role of cation-pi interactions on translation elongation factors. In our investigation, an average of four significant cation-pi interactions were found, that is, an average of one cation-pi interaction per 44 residues in the ten elongation factors were observed. The analysis on the influence of short (< + or - 4), medium (> + or - 4 to < + or - 20) and long (>20) range contacts showed that cation-pi interactions are mainly formed by medium and long-range contacts. Arg-Tyr pair was found largest in number but energetic contribution of Arg-Trp pair was found most. Preferred secondary structural conformation analysis of the residues involved in cation-pi interaction indicates that the cationic Arg prefers to be in helix and Lys having equal probability for helix and strand, whereas the aromatic Phe and Trp were found mostly in helix while Tyr in strand regions. The cation-pi interaction residues involved in these proteins were found highly conserved with 48.86% residues having conservation score of > or = 6. Analysis of secondary structure preference of the energetically significant cation-pi residues in different solvent accessible range indicates that most of the pi residues are found buried or partially buried whereas cationic residues were found mostly at the protein surface. The results presented in this study will be useful for structural stability studies in translation elongation factors.

Betaalpha-hairpin clamps brace betaalphabeta modules and can make substantive contributions to the stability of TIM barrel proteins

Non-local hydrogen bonding interactions between main chain amide hydrogen atoms and polar side chain acceptors that bracket consecutive βα or αβ elements of secondary structure in αTS from E. coli, a TIM barrel protein, have previously been found to contribute 4–6 kcal mol⁻¹ to the stability of the native conformation. Experimental analysis of similar βα-hairpin clamps in a homologous pair of TIM barrel proteins of low sequence identity, IGPS from S. solfataricus and E. coli, reveals that this dramatic enhancement of stability is not unique to αTS. A survey of 71 TIM barrel proteins demonstrates a 4-fold symmetry for the placement of βα-hairpin clamps, bracing the fundamental βαβ building block and defining its register in the (βα)₈ motif. The preferred sequences and locations of βα-hairpin clamps will enhance structure prediction algorithms and provide a strategy for engineering stability in TIM barrel proteins.

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.

Investigations on unconventional hydrogen bonds in RNA binding proteins: The role of CH⋯OC interactions

We have investigated the roles played by C-H...O=C interactions in RNA binding proteins. There was an average of 78 CH...O=C interactions per protein and also there was an average of one significant CH...O=C interactions for every 6 residues in the 59 RNA binding proteins studied. Main chain-Main chain (MM) CH...O=C interactions are the predominant type of interactions in RNA binding proteins. The donor atom contribution to CH...O=C interactions was mainly from aliphatic residues. The acceptor atom contribution for MM CH...O=C interactions was mainly from Val, Phe, Leu, Ile, Arg and Ala. The secondary structure preference analysis of CH...O=C interacting residues showed that, Arg, Gln, Glu and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Ile, Leu, Lys, Met, Phe, Trp and Val preferred to be in strand conformation. Most of the CH...O=C interacting polar amino acid residues were solvent exposed while, majority of the CH...O=C interacting non polar residues were excluded from the solvent. Long and medium-range CH...O=C interactions are the predominant type of interactions in RNA binding proteins. More than 50% of CH...O=C interacting residues had a higher conservation score. Significant percentage of CH...O=C interacting residues had one or more stabilization centers. Sixty-six percent of the theoretically predicted stabilizing residues were also involved in CH...O=C interactions and hence these residues may also contribute additional stability to RNA binding proteins.

Thermal-unfolding Reaction of Triosephosphate Isomerase from Trypanosoma cruzi

Thermal denaturation of triosephosphate isomerase from Trypanosoma cruzi was studied by circular dicrhoism and fluorescence spectroscopies. The unfolding transition was found to be highly irreversible even at the very early stages of the reaction. Kinetic studies, allowed us to identify consecutive reactions. Firstly, only the tryptophan environment is altered. Next, changes on the secondary structure and hydrophobic surface exposure measured by 1-anilino-8-naphthalenesulfonate (ANS) binding were observed. Further conformational changes imply additional modifications on the secondary and tertiary structures and release of the hydrophobic dye leading to the formation of the unfolded state that is prone to aggregate.

Investigations on C–H⋯π interactions in RNA binding proteins

We have investigated the roles played by C-Hcdots, three dots, centeredpi interactions in RNA binding proteins. There was an average of 55 C-Hcdots, three dots, centeredpi interactions per protein and also there was an average of one significant C-Hcdots, three dots, centeredpi interaction for every nine residues in the 59 RNA binding proteins studied. Main-chain to side-chain C-Hcdots, three dots, centeredpi interactions is the predominant type of interactions in RNA binding proteins. The donor atom contribution to C-Hcdots, three dots, centeredpi interactions was mainly from Phe, Tyr, Trp, Pro, Gly, Lys, His and Ala residues. The acceptor atom contribution to main-chain to side-chain C-Hcdots, three dots, centeredpi and side-chain to side-chain C-Hcdots, three dots, centeredpi interactions was mainly from Phe and Tyr residues. On the contrary, the acceptor atoms of Trp residues contributed to all the four types of C-Hcdots, three dots, centeredpi interactions. Also, Trp contributed both donor and acceptor atoms in main-chain to side-chain, main-chain to side-chain five-member aromatic ring and side-chain to side-chain C-Hcdots, three dots, centeredpi interactions. The secondary structure preference analysis of C-Hcdots, three dots, centeredpi interacting residues showed that, Arg, Gln, Glu, His, Ile, Leu, Lys, Met, Phe and Tyr preferred to be in helix, while Ala, Asp, Cys, Gly, Trp and Val preferred to be in strand conformation. Long-range C-Hcdots, three dots, centeredpi interactions are the predominant type of interactions in RNA binding proteins. More than 50% of C-Hcdots, three dots, centeredpi interacting residues had a higher conservation score. Significant percentage of C-Hcdots, three dots, centeredpi interacting residues had one or more stabilization centers. Seven percent of the theoretically predicted stabilizing residues were also involved in C-Hcdots, three dots, centeredpi interactions and hence these residues may also contribute additional stability to RNA binding proteins.

Exploring the role of cation–π interactions in glycoproteins lipid-binding proteins and RNA-binding proteins

We have analyzed and compared the influence of cation-pi interactions in glycoproteins (GPs), lipid-binding proteins (LBPs) and RNA-binding proteins (RBPs) in this study. We observed that all the proteins included in the study had profound cation-pi interactions. There is an average of one energetically significant cation-pi interaction for every 71 residues in GPs, for every 58 residues in LBPs and for every 64 residues in RBPs. Long-range contacts are predominant in all the three types of proteins studied. The pair-wise cation-pi interaction energy between the positively charged and aromatic residues shows that Arg-Trp pair energy was the strongest among all six possible pairs in all the three types of proteins studied. There were considerable differences in the preference of cation-pi interacting residues to different secondary structure elements and ASA and these might contribute to differences in biochemical functions of GPs, LBPs and RBPs. It was interesting to note that all the five residues involved in cation-pi interactions were found to have stabilization centers in GPs, LBPs and RBPs. Majority of the cation-pi interacting residues investigated in the present study had a conservation score of 6, the cutoff value used to identify the stabilizing residues. A small percentage of cation-pi interacting residues were also present as stabilizing residues. The cation-pi interaction-forming residues play an important role in the structural stability of in GPs, LBPs and RBPs. The results obtained in this study will be helpful in further understanding the stability, specificity and differences in the biochemical functions of GPs, LBPs and RBPs.

Mapping the structure of folding cores in TIM barrel proteins by hydrogen exchange mass spectrometry: the roles of motif and sequence for the indole-3-glycerol phosphate synthase from Sulfolobus solfataricus

To test the roles of motif and amino acid sequence in the folding mechanisms of TIM barrel proteins, hydrogen-deuterium exchange was used to explore the structure of the stable folding intermediates for the of indole-3-glycerol phosphate synthase from Sulfolobus solfataricus (sIGPS). Previous studies of the urea denaturation of sIGPS revealed the presence of an intermediate that is highly populated at approximately 4.5 M urea and contains approximately 50% of the secondary structure of the native (N) state. Kinetic studies showed that this apparent equilibrium intermediate is actually comprised of two thermodynamically distinct species, I(a) and I(b). To probe the location of the secondary structure in this pair of stable on-pathway intermediates, the equilibrium unfolding process of sIGPS was monitored by hydrogen-deuterium exchange mass spectrometry. The intact protein and pepsin-digested fragments were studied at various concentrations of urea by electrospray and matrix-assisted laser desorption ionization time-of-flight mass spectrometry, respectively. Intact sIGPS strongly protects at least 54 amide protons from hydrogen-deuterium exchange in the intermediate states, demonstrating the presence of stable folded cores. When the protection patterns and the exchange mechanisms for the peptides are considered with the proposed folding mechanism, the results can be interpreted to define the structural boundaries of I(a) and I(b). Comparison of these results with previous hydrogen-deuterium exchange studies on another TIM barrel protein of low sequence identify, alpha-tryptophan synthase (alphaTS), indicates that the thermodynamic states corresponding to the folding intermediates are better conserved than their structures. Although the TIM barrel motif appears to define the basic features of the folding free energy surface, the structures of the partially folded states that appear during the folding reaction depend on the amino acid sequence. Markedly, the good correlation between the hydrogen-deuterium exchange patterns of sIGPS and alphaTS with the locations of hydrophobic clusters defined by isoleucine, leucine, and valine residues suggests that branch aliphatic side-chains play a critical role in defining the structures of the equilibrium intermediates.

Bioinformatic analysis of protein structure-function relationships: case study of leukocyte elastase (ELA2) missense mutations

Cyclic and congenital neutropenia are caused by mutations in the human neutrophil elastase (HNE) gene (ELA2), leading to an immunodeficiency characterized by decreased or oscillating levels of neutrophils in the blood. The HNE mutations presumably cause loss of enzyme activity, consequently leading to compromised immune system function. To understand the structural basis for the disease, we implemented methods from bioinformatics to analyze all the known HNE missense mutations at both the sequence and structural level. Our results demonstrate that the 32 different mutations have diverse effects on HNE structure and function, affecting structural disorder and aggregation tendencies, stability maintaining contacts, and electrostatic properties. A large proportion of the mutations are located at conserved amino acids, which are usually essential in determining protein structure and function. The majority of the disease-causing HNE missense mutations lead to major structural changes and loss of stability in the protein. A few mutations also affect functional residues, leading into decreased catalytic activity or altered ligand binding. Our analysis reveals the putative effects of all known missense mutations in HNE, thus allowing the structural basis of cyclic and congenital neutropenia to be elucidated. We have employed and analyzed a set of some 30 different methods for predicting the effects of amino acid substitutions. We present results and experience from the analysis of the applicability of these methods in the analysis of numerous genes, proteins, and diseases to reveal protein structure-function relationships and disease genotype-phenotype correlations.