The blind watchmaker and rational protein engineering

Interpreting conformational effects in protein nano-ESI-MS spectra

Nano-electrospray-ionization mass spectrometry (nano-ESI-MS) is employed here to describe equilibrium protein conformational transitions and to analyze the influence of instrumental settings, pH, and solvent surface tension on the charge-state distributions (CSD). A first set of experiments shows that high flow rates of N(2) as curtain gas can induce unfolding of cytochrome c (cyt c) and myoglobin (Mb), under conditions in which the stability of the native protein structure has already been reduced by acidification. However, it is possible to identify conditions under which the instrumental settings are not limiting factors for the conformational stability of the protein inside ESI droplets. Under such conditions, equilibrium unfolding transitions described by ESI-MS are comparable with those obtained by other established biophysical methods. Experiments with the very stable proteins ubiquitin (Ubq) and lysozyme (Lyz) enable testing of the influence of extreme pH changes on the ESI process, uncoupled from acid-induced unfolding. When HCl is used for acidification, Ubq and Lyz mass spectra do not change between pH~7 and pH 2.2, indicating that the CSD is highly characteristic of a given protein conformation and not directly affected by even large pH changes. Use of formic or acetic acid for acidification of Ubq solutions results in major spectral changes that can be interpreted in terms of protein unfolding as a result of the increased hydrophobicity of the solvent. On the other hand, Lyz, cyt c, and Mb enable direct comparison of protein CSD (corresponding to either the folded or the unfolded protein) in HCl or acetic acid solutions at low pH. The values of surface tension for these solutions differ significantly. Confirming indications already present in the literature, we observe very similar CSD under these solvent conditions for several proteins in either compact or disordered conformations. The same is true for comparison between water and water-acetic acid for folded cyt c and Lyz. Thus, protein CSD from water-acetic solutions do not seem to be limited by the low surface tension of acetic acid as previously suggested. This result could reflect a general lack of dependence of protein CSD on the surface tension of the solvent. However, it is also possible that the effect of acetic acid on the precursor ESI droplets is smaller than generally assumed.

Engineering the pH-optimum of a triglyceride lipase: from predictions based on electrostatic computations to experimental results

The optimisation of enzymes for particular purposes or conditions remains an important target in virtually all protein engineering endeavours. Here, we present a successful strategy for altering the pH-optimum of the triglyceride lipase cutinase from Fusarium solani pisi. The computed electrostatic pH-dependent potentials in the active site environment are correlated with the experimentally observed enzymatic activities. At pH-optimum a distinct negative potential is present in all the lipases and esterases that we studied so far. This has prompted us to propose the "The Electrostatic Catapult Model" as a model for product release after cleavage of the ester bond. The origin of the negative potential is associated with the titration status of specific residues in the vicinity of the active site cleft. In the case of cutinase, the role of Glu44 was systematically investigated by mutations into Ala and Lys. Also, the neighbouring Thr45 was mutated into Proline, with the aim of shifting the spatial location of Glu44. All the charge mutants displayed altered titration behaviour of active site electrostatic potentials. Typically, the substitution of the residue Glu44 pushes the onset of the active site negative potential towards more alkaline conditions. We, therefore, predicted more alkaline pH optima, and this was indeed the experimentally observed. Finally, it was found that the pH-dependent computed Coulombic energy displayed a strong correlation with the observed melting temperatures of native cutinase.

How do lipases and esterases work: the electrostatic contribution

This work explores the role of one of the factors explaining lipase/esterase activity: the contribution of electrostatic interactions to lipase/esterase activity. The electrostatic potential distribution on the molecular surface of an enzyme as a function of pH determines, to a large extent, the enzyme's pH activity profile. Other important factors include the presence and distribution of polar and hydrophobic residues in the active cleft. We have mapped the electrostatic potential distribution as a function of pH on the molecular surface of nine lipases/esterases for which the 3D structure is experimentally known. A comparison of these potential maps at different pH values with the corresponding pH-activity profile, pH optimum or pH range where the activity displayed by the enzyme is maximum, has revealed a considerable correlation. A negative potential in the active site appears correlated with maximum activity towards triglycerides, which has prompted us to propose a model for product release ('The electrostatic catapult model') after cleavage of an ester bond. At the same time as the bottom of the active site cleft becomes negatively charged, other nearby regions also titrate and become negatively charged when pH becomes more alkaline, for some of the studied lipases. If such lipases also show phospholipase activity (such as guinea pig lipase-related proteins 2 chimera) we raise the hypothesis that such other titratable regions after becoming negatively charged might stabilise the positive charge present in the polar head of phospholipids, such as phosphatidylcholine, phosphatidylethanolamine and phosphatidylserine. The distribution of polar, weak polar and non-polar residues on the molecular surface of each studied lipase, in particular the active site region, was compared for all the lipases studied. The combination of graphical visualisation of the electrostatic potential maps and the polarity maps combined with knowledge about the location of key residues on the protein surface allows us to envision atomic models for lipolytic activity.

Amino acid neighbours and detailed conformational analysis of cysteines in proteins

Here we present an investigation of the contacts that cysteines make with residues in their three-dimensional environment and a comprehensive analysis of the conformational features of 351 disulphide bridges in 131 non-homologous single-chain protein structures. Upstream half-cystines preferentially have downstream neighbours, whereas downstream half-cystines have mainly upstream neighbours. Non-disulphide bridged cysteines (free cysteines) have no preference for upstream or downstream neighbours. Free cysteines have more contacts to non-polar residues and fewer contacts to polar/charged residues than half-cystines, which correlates with our observation that free cysteines are more buried than half-cystines. Free cysteines prefer to be located in alpha-helices while no clear preference is observed for half-cystines. Histidine and methionine are preferentially seen nearby free cysteines. Tryptophan is found preferentially nearby half-cystines. We have merged sequential and spatial information, and highly interesting novel patterns have been discovered. The number of cysteines per protein is typically an even number, peaking at four. The number of residues separating two half-cystines is preferentially 11 and 16. Left-handed and right-handed disulphide bridges display different conformational parameters. Here we present side chain torsion angle information based on a 5-12 times larger number of disulphide bridges than has previously been published. Considering the importance of cysteines for maintaining the 3D-structural scaffold of proteins, it is essential to have as accurate information as possible concerning the packing and conformational preferences. The present work may provide key information for engineering the protein environment around cysteines.

Protein engineering the surface of enzymes

The protein surface is the interface through which a protein senses the external world. Its composition of charged, polar and hydrophobic residues is crucial for the stability and activity of the protein. The charge state of seven of the twenty naturally occurring amino acids is pH dependent. A total of 95% of all titratable residues are located on the surface of soluble proteins. In evolutionary related families of proteins such residues are particularly prone to substitutions, insertions and deletions. We present here an analysis of the residue composition of 4038 proteins, selected from 125 protein families with < 25% identity between core members of each family. Whereas only 16.8% of the residues were truly buried, 40.7% were > 30% exposed on the surface and the remainder were < 30% exposed. The individual residue types show distinct differences. The data presented provides an important new approach to protein engineering of protein surfaces. Guidelines for the optimization of solvent exposure for a given residue are given. The cutinase family of enzymes has been investigated. The stability of native cutinase has been studied as a function of pH, and has been compared with the cutinase activity towards tributyrin. Whereas the onset of enzymatic activity is linked with the deprotonation of the active site HIS188, destabilization of the 3D structure as determined by differential scanning calorimetry is coupled with the loss of activity at very basic pH values. A modeling investigation of the pH dependence of the electrostatic potentials reveals that the activity range is accompanied by the development of a highly significant negative potential in the active site cleft. The 3D structures of three mutants of the Fusarium solani pisi cutinase have been solved to high resolution using X-ray diffraction analysis. Preliminary X-ray data are presented.

Lipases and esterases: a review of their sequences, structure and evolution

This chapter aims to provide a brief review on the enzyme family of lipases and esterases. The sequences, 3D structures and pH dependent electrostatic signatures are presented and analyzed. Since the family comprises more than 100 sequences, we have tried to focus on the most interesting features from our perspective, which translates into finding similarities and differences between members of this family, in particular in and around the active sites, and to identify residues that are partially or totally conserved. Such residues we believe are either important for maintaining the structural scaf-fold of the protein or to maintain activity or specificity. The structure function relationship for these proteins is therefore of central interest. Can we uniquely identify a protein from this large family of sequences--and if so, what is the identifier? The protein family displays some highly complex features: many of the proteins are interfacially activated, i.e. they need to be in physical contact with the aggregated substrate. Access to the active site is blocked with either a loop fragment or an alpha-helical fragment in the absence of interfacial contact. Although the number of known, relevant protein 3D structures is growing steadily, we are nevertheless faced with a virtual explosion in the number of known or deduced amino acid sequences. It is therefore unrealistic to expect that all protein sequences within the foreseeable future will have their 3D structure determined by X-ray diffractional analysis or through other methods. When feasible the gene and/or the amino acid sequences will be analyzed from an evolutionary perspective. As the 3D folds are often remarkably similar, both among the triglyceride lipases as well as among the esterases, the functional diversities (e.g. specificity) must originate in differences in surface residue utilization, in particular of charged residues. The pH variations in the isopotential surfaces of some of the most interesting lipases are presented and a qualitative interpretation proposed. Finally we illustrate that NMR has potential for becoming an important tool in the study of lipases, esterases and their kinetics.