Abstract: Large amounts of data are being generated annually on the connection between the sequence, structure and function of proteins using site-directed mutagenesis, protein design and directed evolution techniques. These data provide the fundamental building blocks for our understanding of protein function, molecular biology and living organisms in general. However, much experimental data are never deposited in databases and is thus 'lost' in journal publications or in PhD theses. At the same time theoretical scientists are in need of large amounts of experimental data for benchmarking and calibrating novel predictive algorithms, and theoretical progress is therefore often hampered by the lack of suitable data to validate or disprove a theoretical assumption. We present PEAT (Protein Engineering Analysis Tool), an application that integrates data deposition, storage and analysis for researchers carrying out protein engineering projects or biophysical characterization of proteins. PEAT contains modules for DNA sequence manipulation, primer design, fitting of biophysical characterization data (enzyme kinetics, circular dichroism spectroscopy, NMR titration data, etc.), and facilitates sharing of experimental data and analyses for a typical university-based research group. PEAT is freely available to academic researchers at http://enzyme.ucd.ie/PEAT.
Abstract: pH-induced chemical shift perturbations (CSPs) can be used to study pH-dependent conformational transitions in proteins. Recently, an elegant principal component analysis (PCA) algorithm was developed and used to study the pH-dependent structural transitions in bovine beta-lactoglobulin (betaLG) by analyzing its NMR pH-titration spectra. Here, we augment this analysis method by filtering out changes in the NMR chemical shift that stem from effects that are electrostatic in nature. Specifically, we examine how many CSPs can be explained by purely electrostatic effects arising from titrational events in betaLG. The results show that around 20% of the amide nuclei CSPs in betaLG originate exclusively from "through-space" electric field effects. A PCA of NMR data where electric field artefacts have been removed gives a different picture of the pH-dependent structural transitions in betaLG. The method implemented here is well suited to be applied on a whole range of proteins, which experience at least one pH-dependent conformational change. Proteins 2010. (c) 2009 Wiley-Liss, Inc.
Abstract: NMR-monitored pH titration experiments are routinely used to measure site-specific protein pKa values. Accurate experimental pKa values are essential in dissecting enzyme catalysis, in studying the pH-dependence of protein stability and ligand binding, in benchmarking pKa prediction algorithms, and ultimately in understanding electrostatic effects in proteins. However, due to the complex ways in which pH-dependent electrostatic and structural changes manifest themselves in NMR spectra, reported apparent pKa values are often dependent on the way that NMR pH-titration curves are analyzed. It is therefore important to retain the raw NMR spectroscopic data to allow for documentation and possible re-interpretation. We have constructed a database of primary NMR pH-titration data, which is accessible via a web interface. Here, we report statistics of the database contents and analyze the data with a global perspective to provide guidelines on best practice for fitting NMR titration curves. Titration_DB is available at http://enzyme.ucd.ie/Titration_DB. Proteins 2010. (c) 2009 Wiley-Liss, Inc.
Abstract: The glycosaminoglycan heparin is known to possess antimetastatic activity in experimental models and preclinical studies, but there is still uncertainty over its mechanism of action in this respect. As an anticoagulant, heparin enhances inhibition of thrombin by the serpin antithrombin III, but a similar cofactor role has not been previously investigated for proteases linked to metastasis. The squamous cell carcinoma antigens (serpins B3 and B4) are tumor-associated proteins that can inhibit papain-like cysteine proteases, including cathepsins L, K, and S. In this study, we show that SCCA-1 (B3) and SCCA-2 (B4) can bind heparin as demonstrated by affinity chromatography, native PAGE gel shifts, and intrinsic fluorescence quenching. Binding was specific for heparin and heparan sulfate but not other glycosaminoglycans. The presence of heparin accelerated inhibition of cathepsin L by both serpins, and in the case of SCCA-1, heparin increased the second order inhibition rate constant from 5.4 x 10(5) to >10(8), indicating a rate enhancement of at least 180-fold. A templating mechanism was shown, consistent with ternary complex formation. Furthermore, SCCA-1 inhibition of cathepsin L-like proteolytic activity secreted from breast and melanoma cancer cell lines was significantly enhanced by heparin. This is the first example of glycosaminoglycan enhancement of B-clade serpin activity and the first report of heparin acting as a cofactor in serpin cross-class inhibition of cysteine proteases. Most importantly, this finding raises the possibility that the anticancer properties of heparin may be due, at least partly, to enhanced inhibition of prometastatic proteases.
Abstract: The development of docking scoring functions requires high-resolution 3D structures of protein-ligand complexes for which the binding affinity of the ligand has been measured experimentally. Protein-ligand binding affinities are measured in solution experiments, and high resolution protein-ligand structures can be determined only by X-ray crystallography. Protein-ligand scoring functions must therefore reproduce solution binding energies using analyses of proteins in a crystal environment. We present an analysis of the prevalence of crystal-induced artifacts and water-mediated contacts in protein-ligand complexes and demonstrate the effect that these can have on the performance of protein-ligand scoring functions. We find 36% of ligands in the PDBBind 2007 refined data set to be influenced by crystal contacts and find the performance of a scoring function to be affected by these. A Web server for detecting crystal contacts in protein-ligand complexes is available at http://enzyme.ucd.ie/LIGCRYST .
Abstract: Hyaluronidases from diverse species and sources have different pH optima. Distinct mechanisms with regard to dynamic structural changes, which control hyaluronidase activity at varying pH, are unknown. Human serum hyaluronidase 1 (HYAL1) is active solely below pH 5.1. Here we report the design of a HYAL1 variant that degrades hyaluronan up to pH 5.9. Besides highly conserved residues in close proximity of the active site of most hyaluronidases, we identified a bulky loop formation located at the end of the substrate binding crevice of HYAL1 to be crucial for substrate hydrolysis. The stretch between cysteine residues 207 and 221, which normally contains 13 amino acids, could be replaced by a tetrapeptide sequence of alternating glycine serine residues, thereby yielding an active enzyme with an extended binding cleft. This variant exhibited hyaluronan degradation at elevated pH. This is indicative for appropriate substrate binding and proper positioning being decisively affected by sites far off from the active center.
Abstract: The pH dependence of protein biophysical characteristics is often analyzed to gain an improved understanding of protein stability, enzyme activity, and protein-ligand-binding processes. Indeed, much of our understanding of the catalytic mechanisms of enzymes derives from studies of the pH dependence of catalytic activity, and the ability to redesign the pH-dependent properties of enzymes continues to be of high relevance for both industrial and medical applications of proteins. This chapter discusses current theoretical methods for calculating protein pK(a) values and illustrates how one can analyze protein pK(a) calculation results to study calculation accuracy, pH stability profiles, and enzymatic pH activity profiles. A description of how one can analyze the importance of individual titratable groups is presented along with details on methods for redesigning protein pK(a) values and enzymatic pH activity profiles. Finally, I discuss novel methods for fitting experimental nuclear magnetic resonance titration curves and enzymatic pH activity profiles that can be used to derive information on electrostatic interaction energies in proteins.
Abstract: The pH-dependence of the NMR chemical shift for titratable groups in proteins often deviate from a standard Henderson-Hasselbalch (HH) titration curve. A non-HH dependence of the chemical shift for a given residue can arise from a single-site, non-HH titrational event for that residue, or if the chemical shift of the group is influenced by additional titrational events occurring in other residues. We show that simultaneous fits of several non-HH NMR titration curves of interacting protein residues to a statistical mechanical model can be used to distinguish between these two cases. From fitting of non-HH titrations, we can extract electrostatic interaction energies between protein residues. Furthermore, by performing simultaneous fits of NMR titration curves and enzymatic pH-activity profiles, we can gain information on the identity and populations of the catalytically competent protonation states in enzymes. We apply the global fitting of titrational events (GloFTE) method to experimental data on five enzyme systems and on a single non-enzyme system, and show that the extracted electrostatic interaction energies and effective dielectric constants for a subset of these systems agree excellently with experimentally determined values as well as with theoretical calculations. In the case of reduced Escherichia coli thioredoxin we use GloFTE analysis to distinguish between two possible interpretations of the NMR titration curves of the active site residues. We also show that for the strongly coupled system of titratable groups in the active site of the Bacillus circulans xylanase (BCX) N35D mutant, GloFTE fits of a single titration curve and an enzymatic pH-activity profile can give a full description of the energetics of the titrational events in the enzyme's active site. Using only the X-ray crystallographic structure of the enzyme and the electrostatic interaction energies extracted from such a GloFTE fit, we can uniquely identify the three catalytic groups in this system. This raises the prospect of completely characterising active site titrational events from a single unassigned NMR titration curve and an enzymatic pH-activity profile.
Abstract: Real-world observable physical and chemical characteristics are increasingly being calculated from the 3D structures of biomolecules. Methods for calculating pK(a) values, binding constants of ligands, and changes in protein stability are readily available, but often the limiting step in computational biology is the conversion of PDB structures into formats ready for use with biomolecular simulation software. The continued sophistication and integration of biomolecular simulation methods for systems- and genome-wide studies requires a fast, robust, physically realistic and standardized protocol for preparing macromolecular structures for biophysical algorithms. As described previously, the PDB2PQR web server addresses this need for electrostatic field calculations (Dolinsky et al., Nucleic Acids Research, 32, W665-W667, 2004). Here we report the significantly expanded PDB2PQR that includes the following features: robust standalone command line support, improved pK(a) estimation via the PROPKA framework, ligand parameterization via PEOE_PB charge methodology, expanded set of force fields and easily incorporated user-defined parameters via XML input files, and improvement of atom addition and optimization code. These features are available through a new web interface (http://pdb2pqr.sourceforge.net/), which offers users a wide range of options for PDB file conversion, modification and parameterization.
Abstract: We have characterized by NMR spectroscopy the three active site (His80, His85, and His205) and two non-active site (His107 and His114) histidines in the 34 kDa catalytic domain of Cellulomonas fimi xylanase Cex in its apo, noncovalently aza-sugar-inhibited, and trapped glycosyl-enzyme intermediate states. Due to protection from hydrogen exchange, the level of which increased upon inhibition, the labile 1Hdelta1 and 1H epsilon1 atoms of four histidines (t1/2 approximately 0.1-300 s at 30 degrees C and pH approximately 7), as well as the nitrogen-bonded protons in the xylobio-imidazole and -isofagomine inhibitors, could be observed with chemical shifts between 10.2 and 17.6 ppm. The histidine pKa values and neutral tautomeric forms were determined from their pH-dependent 13C epsilon1-1H epsilon1 chemical shifts, combined with multiple-bond 1H delta2/epsilon1-15N delta1/epsilon2 scalar coupling patterns. Remarkably, these pKa values span more than 8 log units such that at the pH optimum of approximately 6 for Cex activity, His107 and His205 are positively charged (pKa > 10.4), His85 is neutral (pKa < 2.8), and both His80 (pKa = 7.9) and His114 (pKa = 8.1) are titrating between charged and neutral states. Furthermore, upon formation of the glycosyl-enzyme intermediate, the pKa value of His80 drops from 7.9 to <2.8, becoming neutral and accepting a hydrogen bond from an exocyclic oxygen of the bound sugar moiety. Changes in the pH-dependent activity of Cex due to mutation of His80 to an alanine confirm the importance of this interaction. The diverse ionization behaviors of the histidine residues are discussed in terms of their structural and functional roles in this model glycoside hydrolase.
Abstract: In this paper we present a method for the multi-resolution comparison of biomolecular electrostatic potentials without the need for global structural alignment of the biomolecules. The underlying computational geometry algorithm uses multi-resolution attributed contour trees (MACTs) to compare the topological features of volumetric scalar fields. We apply the MACTs to compute electrostatic similarity metrics for a large set of protein chains with varying degrees of sequence, structure, and function similarity. For calibration, we also compute similarity metrics for these chains by a more traditional approach based upon 3D structural alignment and analysis of Carbo similarity indices. Moreover, because the MACT approach does not rely upon pairwise structural alignment, its accuracy and efficiency promises to perform well on future large-scale classification efforts across groups of structurally-diverse proteins. The MACT method discriminates between protein chains at a level comparable to the Carbo similarity index method; i.e., it is able to accurately cluster proteins into functionally-relevant groups which demonstrate strong dependence on ligand binding sites. The results of the analyses are available from the linked web databases http://ccvweb.cres.utexas.edu/MolSignature/ and http://agave.wustl.edu/similarity/. The MACT analysis tools are available as part of the public domain library of the Topological Analysis and Quantitative Tools (TAQT) from the Center of Computational Visualization, at the University of Texas at Austin (http://ccvweb.csres.utexas.edu/software). The Carbo software is available for download with the open-source APBS software package at http://apbs.sf.net/.
Abstract: pK(a) calculations for macromolecules are normally performed by solving the Poisson-Boltzmann equation, accounting for the different dielectric constants of solvent and solute, as well as the ionic strength. Despite the large number of successful applications, there are some situations where the current algorithms are not suitable: (1) large scale, high-throughput analysis which requires calculations to be completed within a fraction of a second, e.g. when permanently monitoring pK(a) shifts during a molecular dynamics simulation; (2) prediction of pK(a)s in periodic boundaries, e.g. when reconstructing entire protein crystal unit cells from PDB files, including the correct protonation patterns at experimental pH. Such in silico crystals are needed by 'self-parameterizing' molecular dynamics force fields like YASARA YAMBER, that optimize their parameters while energy-minimizing high-resolution protein crystals. To address both problems, we define an empirical equation that expresses the pK(a) as a function of electrostatic potential, hydrogen bonds and accessible surface area. The electrostatic potential is evaluated by Ewald summation, which captures periodic crystal environments and the uncertainty in atom positions using Gaussian charge densities. The empirical proportionality constants are derived from 217 experimentally determined pK(a)s, and despite its simplicity, this pK(a) calculation method reaches a high overall jack-knifed accuracy, and is fast enough to be used during a molecular dynamics simulation. A reliable null-model to judge pK(a) prediction accuracies is also presented.
Abstract: Continuum solvation models, such as Poisson-Boltzmann and Generalized Born methods, have become increasingly popular tools for investigating the influence of electrostatics on biomolecular structure, energetics and dynamics. However, the use of such methods requires accurate and complete structural data as well as force field parameters such as atomic charges and radii. Unfortunately, the limiting step in continuum electrostatics calculations is often the addition of missing atomic coordinates to molecular structures from the Protein Data Bank and the assignment of parameters to biomolecular structures. To address this problem, we have developed the PDB2PQR web service (http://agave.wustl.edu/pdb2pqr/). This server automates many of the common tasks of preparing structures for continuum electrostatics calculations, including adding a limited number of missing heavy atoms to biomolecular structures, estimating titration states and protonating biomolecules in a manner consistent with favorable hydrogen bonding, assigning charge and radius parameters from a variety of force fields, and finally generating 'PQR' output compatible with several popular computational biology packages. This service is intended to facilitate the setup and execution of electrostatics calculations for both experts and non-experts and thereby broaden the accessibility to the biological community of continuum electrostatics analyses of biomolecular systems.
Abstract: Glycoside hydrolysis by retaining family 18 chitinases involves a catalytic acid (Glu) which is part of a conserved DXDXE sequence motif that spans strand four of a (betaalpha)8 barrel (TIM barrel) structure. These glycoside hydrolases are unusual in that the positive charge emerging on the anomeric carbon after departure of the leaving group is stabilized by the substrate itself (the N-acetyl group of the distorted -1 sugar), rather than by a carboxylate group on the enzyme. We have studied seven conserved residues in the catalytic center of chitinase B from Serratia marcescens. Putative roles for these residues are proposed on the basis of the observed mutational effects, the pH-dependency of these effects, pKa calculations and available structural information. The results indicate that the pKa of the catalytic acid (Glu144) is 'cycled' during catalysis as a consequence of substrate-binding and release and, possibly, by a back and forth movement of Asp142 between Asp140 and Glu144. Rotation of Asp142 towards Glu144 also contributes to an essential distortion of the N-acetyl group of the -1 sugar. Two other conserved residues (Tyr10 and Ser93) are important because they stabilize the charge on Asp140 while Asp142 points towards Glu144. Asp215, lying opposite Glu144 on the other side of the scissile glycosidic bond, contributes to catalysis by promoting distortion of the -1 sugar and by increasing the pKa of the catalytic acid. The hydroxyl group of Tyr214 makes a major contribution to the positioning of the N-acetyl group of the -1 sugar. Taken together, the results show that catalysis in family 18 chitinases depends on a relatively large number of (partly mobile) residues that interact with each other and the substrate.
Abstract: The ionization properties of the active-site residues in enzymes are of considerable interest in the study of the catalytic mechanisms of enzymes. Knowledge of these ionization constants (pKa values) often allows the researcher to identify the proton donor and the catalytic nucleophile in the reaction mechanism of the enzyme. Estimates of protein residue pKa values can be obtained by applying pKa calculation algorithms to protein X-ray structures. We show that pKa values accurate enough for identifying the proton donor in an enzyme active site can be calculated by considering in detail only the active-site residues and their immediate electrostatic interaction partners, thus allowing for a large decrease in calculation time. More specifically we omit the calculation of site-site interaction energies, and the calculation of desolvation and background interaction energies for a large number of pairs of titratable groups. The method presented here is well suited to be applied on a genomic scale, and can be implemented in most pKa calculation algorithms to give significant reductions in calculation time with little or no impact on the accuracy of the results. The work presented here has implications for the understanding of enzymes in general and for the design of novel biocatalysts.
Abstract: The calculation of the physical properties of a protein from its X-ray structure is of importance in virtually every aspect of modern biology. Although computational algorithms have been developed for calculating everything from the dynamics of a protein to its binding specificity, only limited information is available on the ability of these methods to give accurate results when used with a particular X-ray structure. We examine the ability of a pKa calculation algorithm to predict the proton-donating residue in the catalytic mechanism of hen egg white lysozyme. We examine the correlation between the ability of the pKa calculation method to obtain the correct result and the overall characteristics of 41 X-ray structures such as crystallization conditions, resolution, and the output of structure validation software. We furthermore examine the ability of energy minimizations (EM), molecular dynamics (MD) simulations, and structure-perturbation methods to optimize the X-ray structures such that these give correct results with the pKa calculation algorithm. We propose a set of criteria for identifying the proton donor in a catalytic mechanism, and demonstrate that the application of these criteria give highly accurate prediction results when using unmodified X-ray structures. More specifically, we are able to successfully identify the proton donor in 85% of the X-ray structures when excluding structures with crystal contacts near the active site. Neither the use of the overall characteristics of the X-ray structures nor the optimization of the structure by EM, MD, or other methods improves the results of the pKa calculation algorithm. We discuss these results and their implications for the design of structure-based energy calculation algorithms in general.
Abstract: We have developed a computer algorithm, FOLDEF (for FOLD-X energy function), to provide a fast and quantitative estimation of the importance of the interactions contributing to the stability of proteins and protein complexes. The predictive power of FOLDEF was tested on a very large set of point mutants (1088 mutants) spanning most of the structural environments found in proteins. FOLDEF uses a full atomic description of the structure of the proteins. The different energy terms taken into account in FOLDEF have been weighted using empirical data obtained from protein engineering experiments. First, we considered a training database of 339 mutants in nine different proteins and optimised the set of parameters and weighting factors that best accounted for the changes in stability of the mutants. The predictive power of the method was then tested using a blind test mutant database of 667 mutants, as well as a database of 82 protein-protein complex mutants. The global correlation obtained for 95 % of the entire mutant database (1030 mutants) is 0.83 with a standard deviation of 0.81 kcal mol(-1) and a slope of 0.76. The present energy function uses a minimum of computational resources and can therefore easily be used in protein design algorithms, and in the field of protein structure and folding pathways prediction where one requires a fast and accurate energy function. FOLDEF is available via a web-interface at http://fold-x.embl-heidelberg.de
Abstract: The impact of long range electrostatic interactions on catalysis in the thermolysin-like protease from Bacillus stearothermophilus was studied by analyzing the effects of inserting or removing charges on the protein surface. Various mutations were introduced at six different positions, and double-mutant cycle analysis was used to study the extent to which mutational effects were interdependent. The effects of single point mutations on the k(cat)/K(m) were non-additive, even in cases where the point mutations were located 10 A or more from the active site Zn(2+) and separated from each other by up to 25 A. This shows that catalysis is affected by large electrostatic networks that involve major parts of the enzyme. The interdependence of mutations at positions as much as 25 A apart in space also indicates that other effects, such as active site dynamics, play an important role in determining active site electrostatics. Several mutations yielded a significant increase in the activity, the most active (quadruple) mutant being almost four times as active as the wild type. In some cases the shape of the pH-activity profile was changed significantly. Remarkably, large changes in the pH-optimum were not observed.
Abstract: The glycosyl hydrolases present a large family of enzymes that are of great significance for industry. Consequently, there is considerable interest in engineering the enzymes in this family for optimal performance under a range of very diverse conditions. Until recently, tailoring glycosyl hydrolases for specific industrial processes mainly involved stability engineering, but lately there has also been considerable interest in engineering their pH-activity profiles. We mutated four neutral residues (N190, F290, N326 and Q360) in the chimeric Bacillus Ba2 alpha-amylase to both charged and neutral amino acids. The results show that the pH-activity profile of the Ba2 alpha-amylase can be changed by inserting charged residues close to the active site. The changes in the pH-activity profile for these neutral --> charged mutations do not, however, correlate with the predictions from calculations of the p K(a) values of the active site residues. More surprisingly, the neutral --> neutral mutations change the pH-activity profile as much as the neutral --> charged mutations. From these results, it is concluded that factors other than electrostatics, presumably the dynamic aspects of the active site, are important for the shape of the pH-activity profiles of the alpha-amylases.
Abstract: pK(a) calculation methods that are based on finite difference solutions to the Poisson-Boltzmann equation (FDPB) require that energy calculations be performed for a large number of different protonation states of the protein. Normally, the differences between these protonation states are modeled by changing the charges on a few atoms, sometimes the differences are modeled by adding or removing hydrogens, and in a few cases the positions of these hydrogens are optimized locally. We present an FDPB-based pK(a) calculation method in which the hydrogen-bond network is globally optimized for every single protonation state used. This global optimization gives a significant improvement in the accuracy of calculated pK(a) values, especially for buried residues. It is also shown that large errors in calculated pK(a) values are often due to structural artifacts induced by crystal packing. Optimization of the force fields and parameters used in pK(a) calculations should therefore be performed with X-ray structures that are corrected for crystal artifacts.
Abstract: Previous studies of the low molecular mass family 11 xylanase from Bacillus circulans show that the ionization state of the nucleophile (Glu78, pK(a) 4.6) and the acid/base catalyst (Glu172, pK(a) 6.7) gives rise to its pH-dependent activity profile. Inspection of the crystal structure of BCX reveals that Glu78 and Glu172 are in very similar environments and are surrounded by several chemically equivalent and highly conserved active site residues. Hence, there are no obvious reasons why their apparent pK(a) values are different. To address this question, a mutagenic approach was implemented to determine what features establish the pK(a) values (measured directly by (13)C NMR and indirectly by pH-dependent activity profiles) of these two catalytic carboxylic acids. Analysis of several BCX variants indicates that the ionized form of Glu78 is preferentially stabilized over that of Glu172 in part by stronger hydrogen bonds contributed by two well-ordered residues, namely, Tyr69 and Gln127. In addition, theoretical pK(a) calculations show that Glu78 has a lower pK(a) value than Glu172 due to a smaller desolvation energy and more favorable background interactions with permanent partial charges and ionizable groups within the protein. The pK(a) value of Glu172 is in turn elevated due to electrostatic repulsion from the negatively charged glutamate at position 78. The results also indicate that all of the conserved active site residues act concertedly in establishing the pK(a) values of Glu78 and Glu172, with no particular residue being singly more important than any of the others. In general, residues that contribute positive charges and hydrogen bonds serve to lower the pK(a) values of Glu78 and Glu172. The degree to which a hydrogen bond lowers a pK(a) value is largely dependent on the length of the hydrogen bond (shorter bonds lower pK(a) values more) and the chemical nature of the donor (COOH > OH > CONH(2)). In contrast, neighboring carboxyl groups can either lower or raise the pK(a) values of the catalytic glutamic acids depending upon the electrostatic linkage of the ionization constants of the residues involved in the interaction. While the pH optimum of BCX can be shifted from -1.1 to +0.6 pH units by mutating neighboring residues within the active site, activity is usually compromised due to the loss of important ground and/or transition state interactions. These results suggest that the pH optima of an enzyme might be best engineered by making strategic amino acid substitutions, at positions outside of the "core" active site, that electrostatically influence catalytic residues without perturbing their immediate structural environment.
Abstract: alpha-Amylases constitute a very diverse family of glycosyl hydrolases that cleave alpha1-->4 linkages in amylose and related polymers. Recent structural and mutagenic studies of archeael, mammalian and bacterial alpha-amylases have resulted in a wealth of information on the catalytic mechanism and on the structural features of this enzyme class. Because of their high thermo-stability, the Bacillus alpha-amylases have found widespread use in industrial processes, and much attention has been devoted to optimising these enzymes for the very harsh conditions encountered there. Stability has been a major area of focus in this respect, and several remarkably stable bacterial alpha-amylases have been produced by bioengineering techniques. Protein engineering studies of pH-activity profiles and of substrate specificities have also been initiated, although without much success. In the coming years it is likely, however, that the focus of alpha-amylase engineering will shift from engineering stability to these new areas.
Abstract: The amino acid sequence of Leishmania mexicana triose phosphate isomerase is unique in having at position 65 a glutamic acid instead of a glutamine. The stability properties of LmTIM and the E65Q mutant were investigated by pH and guanidinium chloride-induced unfolding. The crystal structure of E65Q was determined. Three important observations were made: (a) there are no structural rearrangements as the result of the substitution; (b) the mutant is more stable than the wild-type; and (c) the stability of the wild-type enzyme shows strong pH dependence, which can be attributed to the ionization of Glu65. Burying of the Glu65 side chain in the uncharged environment of the dimer interface results in a shift in pKa of more than 3 units. The pH-dependent decrease in overall stability is due to weakening of the monomer-monomer interactions (in the dimer). The E65Q substitution causes an increase in stability as the result of the formation of an additional hydrogen bond in each subunit (DeltaDeltaG degrees of 2 kcal.mol-1 per monomer) and the elimination of a charged group in the dimer interface (DeltaDeltaG degrees of at least 9 kcal.mol-1 per dimer). The computated shift in pKa and the stability of the dimer calculated from the charge distribution in the protein structure agree closely with the experimental results. The guanidinium chloride dependence of the unfolding constant was smaller than expected from studies involving monomeric model proteins. No intermediates could be identified in the unfolding equilibrium by combining fluorescence and CD measurements. Study of a stable monomeric triose phosphate isomerase variant confirmed that the phenomenon persists in the monomer.
Abstract: The importance of electrostatics in catalysis has been emphasized in the literature for a large number of enzymes. We examined this hypothesis for the Bacillus licheniformis alpha-amylase by constructing site-directed mutants that were predicted to change the pKa values of the catalytic residues and thus change the pH-activity profile of the enzyme. To change the pKa of the catalytic residues in the active site, we constructed mutations that altered the hydrogen bonding network, mutations that changed the solvent accessibility, and mutations that altered the net charge of the molecule. The results show that changing the hydrogen bonding network near an active site residue or changing the solvent accessibility of an active site residue will very likely result in an enzyme with drastically reduced activity. The differences in the pH-activity profiles for these mutants were modest. pH-activity profiles of mutants which change the net charge on the molecule were significantly different from the wild-type pH-activity profile. The differences were, however, difficult to correlate with the electrostatic field changes calculated. In several cases we observed that pH-activity profiles shifted in the opposite direction compared to the shift predicted from electrostatic calculations. This strongly suggests that electrostatic effects cannot be solely responsible for the pH-activity profile of the B. licheniformis alpha-amylase.
Abstract: Electrostatic interactions play a key role in many aspects of protein engineering. Consequently, much effort has been put into the design of software for calculating electrostatic fields around macromolecules. We show that optimization of hydrogen bonding networks can improve both the results of pK(a) calculations and the results of electrostatic calculations performed by commonly used programs such as DelPhi. Further optimization can often be achieved by flipping the side chains of asparagine, histidine and glutamine around their chi2, chi2 and chi3 torsion angles, respectively, when this improves the local hydrogen bonding network. These optimizations are applied to some well characterized proteins: BPTI, hen egg white lysozyme and superoxide dismutase. A search for flipped residues in the PDB revealed that significant improvements in electrostatic calculations in or near the active site of enzymes can be expected for about one quarter of all enzymes in the PDB.
