Structural insights into the loss of catalytic competence in pectate lyase activity at low pH.

Pectate lyase, a family 1 polysaccharide lyase, catalyses cleavage of the α-1,4 linkage of the polysaccharide homogalacturonan via an anti ß-elimination reaction. In the Michaelis complex two calcium ions bind between the C6 carboxylate of the d-galacturonate residue and enzyme aspartates at the active centre (+1 subsite), they withdraw electrons acidifying the C5 proton facilitating its abstraction by the catalytic arginine. Here we show that activity is lost at low pH because protonation of aspartates results in the loss of the two catalytic calcium-ions causing a profound failure to correctly organise the Michaelis complex.

Remeasuring HEWL pK(a) values by NMR spectroscopy: methods, analysis, accuracy, and implications for theoretical pK(a) calculations.

Site-specific pK(a) values measured by NMR spectroscopy provide essential information on protein electrostatics, the pH-dependence of protein structure, dynamics and function, and constitute an important benchmark for protein pK(a) calculation algorithms. Titration curves can be measured by tracking the NMR chemical shifts of several reporter nuclei versus sample pH. However, careful analysis of these curves is needed to extract residue-specific pK(a) values since pH-dependent chemical shift changes can arise from many sources, including through-bond inductive effects, through-space electric field effects, and conformational changes. We have re-measured titration curves for all carboxylates and His 15 in Hen Egg White Lysozyme (HEWL) by recording the pH-dependent chemical shifts of all backbone amide nitrogens and protons, Asp/Glu side chain protons and carboxyl carbons, and imidazole protonated carbons and protons in this protein. We extracted pK(a) values from the resulting titration curves using standard fitting methods, and compared these values to each other, and with those measured previously by ¹H NMR (Bartik et al., Biophys J 1994;66:1180­1184). This analysis gives insights into the true accuracy associated with experimentally measured pK(a) values. We find that apparent pK(a) values frequently differ by 0.5­1.0 units depending upon the nuclei monitored, and that larger differences occasionally can be observed. The variation in measured pK(a) values, which reflects the difficulty in fitting and assigning pH-dependent chemical shifts to specific ionization equilibria, has significant implications for the experimental procedures used for measuring protein pK(a) values, for the benchmarking of protein pK(a) calculation algorithms, and for the understanding of protein electrostatics in general.

Graphical analysis of pH-dependent properties of proteins predicted using PROPKA.

BACKGROUND: Charge states of ionizable residues in proteins determine their pH-dependent properties through their pKa values. Thus, various theoretical methods to determine ionization constants of residues in biological systems have been developed. One of the more widely used approaches for predicting pKa values in proteins is the PROPKA program, which provides convenient structural rationalization of the predicted pKa values without any additional calculations. RESULTS: The PROPKA Graphical User Interface (GUI) is a new tool for studying the pH-dependent properties of proteins such as charge and stabilization energy. It facilitates a quantitative analysis of pKa values of ionizable residues together with their structural determinants by providing a direct link between the pKa data, predicted by the PROPKA calculations, and the structure via the Visual Molecular Dynamics (VMD) program. The GUI also calculates contributions to the pH-dependent unfolding free energy at a given pH for each ionizable group in the protein. Moreover, the PROPKA-computed pKa values or energy contributions of the ionizable residues in question can be displayed interactively. The PROPKA GUI can also be used for comparing pH-dependent properties of more than one structure at the same time. CONCLUSIONS: The GUI considerably extends the analysis and validation possibilities of the PROPKA approach. The PROPKA GUI can conveniently be used to investigate ionizable groups, and their interactions, of residues with significantly perturbed pKa values or residues that contribute to the stabilization energy the most. Charge-dependent properties can be studied either for a single protein or simultaneously with other homologous structures, which makes it a helpful tool, for instance, in protein design studies or structure-based function predictions. The GUI is implemented as a Tcl/Tk plug-in for VMD, and can be obtained online at http://propka.ki.ku.dk/~luca/wiki/index.php/GUI_Web.

PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions.

In this study, we have revised the rules and parameters for one of the most commonly used empirical pKa predictors, PROPKA, based on better physical description of the desolvation and dielectric response for the protein. We have introduced a new and consistent approach to interpolate the description between the previously distinct classifications into internal and surface residues, which otherwise is found to give rise to an erratic and discontinuous behavior. Since the goal of this study is to lay out the framework and validate the concept, it focuses on Asp and Glu residues where the protein pKa values and structures are assumed to be more reliable. The new and improved implementation is evaluated and discussed; it is found to agree better with experiment than the previous implementation (in parentheses): rmsd = 0.79 (0.91) for Asp and Glu, 0.75 (0.97) for Tyr, 0.65 (0.72) for Lys, and 1.00 (1.37) for His residues. The most significant advance, however, is in reducing the number of outliers and removing unreasonable sensitivity to small structural changes that arise from classifying residues as either internal or surface.

Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of pKa Values.

The new empirical rules for protein pKa predictions implemented in the PROPKA3.0 software package (Olsson et al. J. Chem. Theory Comput.2010, 7, 525-537) have been extended to the prediction of pKa shifts of active site residues and ionizable ligand groups in protein-ligand complexes. We present new algorithms that allow pKa shifts due to inductive (i.e., covalently coupled) intraligand interactions, as well as noncovalently coupled interligand interactions in multiligand complexes, to be included in the prediction. The number of different ligand chemical groups that are automatically recognized has been increased to 18, and the general implementation has been changed so that new functional groups can be added easily by the user, aided by a new and more general protonation scheme. Except for a few cases, the new algorithms in PROPKA3.1 are found to yield results similar to or better than those obtained with PROPKA2.0 (Bas et al. Proteins: Struct., Funct., Bioinf.2008, 73, 765-783). Finally, we present a novel algorithm that identifies noncovalently coupled ionizable groups, where pKa prediction may be especially difficult. This is a general improvement to PROPKA and is applied to proteins with and without ligands.

Integrated prediction of the effect of mutations on multiple protein characteristics.

Site-directed mutagenesis is routinely used in modern biology to elucidate the functional or biophysical roles of protein residues, and plays an important role in the field of rational protein design. Over the past decade, a number of computational tools have been developed that can predict the effect of point mutations on a protein's biophysical characteristics. However, these programs usually provide predictions for only a single characteristic. Furthermore, online versions of these tools are often impractical to use for examination of large and diverse sets of mutants. We have created a new web application, (http://enzyme.ucd.ie/PEAT_SA), that can simultaneously predict the effect of mutations on stability, ligand affinity and pK(a) values. PEAT-SA also provides an expanded feature-set with respect to other online tools which includes the ability to obtain predictions for multiple mutants in one submission. As a result, researchers who use site-directed mutagenesis can access state-of-the-art protein design methods with a fraction of the effort previously required. The results of benchmarking PEAT-SA on standard test-sets demonstrate that its accuracy for all three prediction types compares well to currently available tools. We illustrate PEAT-SA's potential by using it to investigate the influence of mutations on the activity of Subtilisin BPN'. This example demonstrates how the ability to obtain a wide range of information from one source, that can be combined to obtain deeper insight into the influence of mutations, makes PEAT-SA a valuable service to both experimental and computational biologists.

Capturing, sharing and analysing biophysical data from protein engineering and protein characterization studies.

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.

Improving the analysis of NMR spectra tracking pH-induced conformational changes: removing artefacts of the electric field on the NMR chemical shift.

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.

Structural artifacts in protein-ligand X-ray structures: implications for the development of docking scoring functions.

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 .

Determination of electrostatic interaction energies and protonation state populations in enzyme active sites.

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.