Crystal structure of the complex of carboxypeptidase A with a strongly bound phosphonate in a new crystalline form: comparison with structures of other complexes

Phosphinic Peptides as Tool Compounds for the Study of Pharmacologically Relevant Zn-Metalloproteases

Phosphinic peptides constitute an important class of bioactive compounds that have found a wide range of applications in the field of biology and pharmacology of Zn-metalloproteases, the largest family of proteases in humans. They are designed to mimic the structure of natural substrates during their proteolysis, thus acting as mechanism-based, transition state analogue inhibitors. A combination of electrostatic interactions between the phosphinic acid group and the Zn cation as well as optimal noncovalent enzyme-ligand interactions can result in both high binding affinity for the desired target and selectivity against other proteases. Due to these unique properties, phosphinic peptides have been mainly employed as tool compounds for (a) the purposes of rational drug design by serving as ligands in X-ray crystal structures of target enzymes and allowing the identification of crucial interactions that govern optimal molecular recognition, and (b) the delineation of biological pathways where Zn-metalloproteases are key regulators. For the latter objective, inhibitors of the phosphinopeptidic type have been used either unmodified or after being transformed to probes of various types, thus expanding the arsenal of functional tools available to researchers. The aim of this review is to summarize all recent research achievements in which phosphinic peptides have played a central role as tool compounds in the understanding of the mechanism and biological functions of Zn-metalloproteases in both health and disease.

Low barrier hydrogen bonds in protein structure and function

Low-barrier hydrogen bonds (LBHBs) are a special type of short hydrogen bond (HB) that is characterized by the equal sharing of a hydrogen atom. The existence and catalytic role of LBHBs in proteins has been intensely contested. Advancements in X-ray and neutron diffraction methods has revealed delocalized hydrogen atoms involved in potential LBHBs in a number of proteins, while also demonstrating that short HBs are not necessarily LBHBs. More importantly, a series of experiments on ketosteroid isomerase (KSI) have suggested that LBHBs are significantly stronger than standard HBs in the protein microenvironment in terms of enthalpy, but not free energy. The discrepancy between the enthalpy and free energy of LBHBs offers clues to the challenges, and potential solutions, of the LBHB debate, where the unique strength of LBHBs plays a special role in the kinetic processes of enzyme function and structure, together with other molecular forces in a pre-organized environment.

A new benchmark illustrates that integration of geometric constraints inferred from enzyme reaction chemistry can increase enzyme active site modeling accuracy

Enzymes play a critical role in a wide array of industrial, medical, and research applications and with the recent explosion of genomic sequencing, we now have sequences for millions of enzymes for which there is no known structure. In order to utilize modern computational design tools for constructing inhibitors or engineering novel catalysts, the ability to accurately model enzymes is critical. A popular approach for modeling enzymes are comparative modeling techniques which can often accurately predict the global structural features. However, achieving atomic accuracy of an active site remains a challenge and is an issue when trying to utilize the molecular details for designing inhibitors or enhanced catalysts. Here we explore integrating knowledge about the required geometric orientation of conserved catalytic residues into the comparative modeling process in order to improve modeling accuracy. In order to investigate the utility of adding this information, we first carefully construct a benchmark set of reference structures to use. Consistent with previous findings, our benchmark demonstrates that the geometry between catalytic residues across an enzyme family is conserved and does not tend to deviate by more than 0.5Å. We then find that by integrating these geometric constraints during modeling, we can double the number of atomic level accuracy models (<1Å RMSD to the crystal structure ligand) within our benchmarking dataset, even for targets with templates as low as 20-30% sequence identity. Catalytic residues within an enzyme family are highly conserved and can often be readily identified through comparative sequence analysis to a known structure within the enzyme family. Therefore utilizing this readily available information has the potential to significantly improve drug design and enzyme engineering efforts for which there is no known structure for the enzyme of interest.

Predictive methods for computational metalloenzyme redesign - a test case with carboxypeptidase A

Computational metalloenzyme design is a multi-scale problem. It requires treating the metal coordination quantum mechanically, extensive sampling of the protein backbone, and additionally accounting for the polarization of the active site by both the metal cation and the surrounding protein (a phenomenon called electrostatic preorganization). We bring together a combination of theoretical methods that jointly offer these desired qualities: QM/DMD for mixed quantum-classical dynamic sampling, quantum theory of atoms in molecules (QTAIM) for the assessment of electrostatic preorganization, and Density Functional Theory (DFT) for mechanistic studies. Within this suite of principally different methods, there are both complementarity of capabilities and cross-validation. Using these methods, predictions can be made regarding the relative activities of related enzymes, as we show on the native Zn²⁺-dependent carboxypeptidase A (CPA), and its mutant proteins, which are hypothesized to hydrolyze modified substrates. For the native CPA, we replicated the catalytic mechanism and the rate in close agreement with the experiment, giving validity to the QM/DMD predicted structure, the DFT mechanism, and the QTAIM assessment of catalytic activity. For most sequences of the modified substrate and tried CPA mutants, substantially worsened activity is predicted. However, for the substrate mutant that contains Asp instead of Phe at the C-terminus, one CPA mutant exhibits a reasonable activity, as predicted across the theoretical methods. CPA is a well-studied system, and here it serves as a testing ground for the offered methods.

Relationships of human α/β hydrolase fold proteins and other organophosphate-interacting proteins

Organophosphates (OPs) are either found in nature or synthetized for use as pesticides, flame retardants, neurotoxic warfare agents or drugs (cholinergic enhancers in Alzheimer's disease and myasthenia gravis, or inhibitors of lipases in metabolic diseases). Because of the central role of acetylcholinesterase cholinergic neurotransmission in humans, one of the main purposes for using OPs is inactivation of the enzyme by phosphorylation of the nucleophilic serine residue in the active center. However, hundreds of serine hydrolases are expressed in the human proteome, and many of them are potential targets for OP adduction. In this review, we first situate the α/β hydrolase fold proteins among the distinctively folded proteins known to interact with OPs, in particular the different lipases, peptidases, and enzymes hydrolyzing OPs. Second, we compile the human α/β hydrolases and review those that have been experimentally shown to interact with OPs. Among the 120 human α/β hydrolase fold proteins, 102 have a serine in the consensus GXSXG pentapeptide compatible with an active site, 6 have an aspartate or a cysteine as the active site nucleophile residue, and 12 evidently lack an active site. 76 of the 120 have been experimentally shown to bind an OP.

High-resolution structure of the M14-type cytosolic carboxypeptidase from Burkholderia cenocepacia refined exploiting PDB_REDO strategies

A potential cytosolic metallocarboxypeptidase from Burkholderia cenocepacia has been crystallized and a synchrotron-radiation microfocus beamline allowed the acquisition of diffraction data to 1.9 Å resolution. The asymmetric unit comprises a tetramer containing over 1500 amino acids, and the high-throughput automated protocols embedded in PDB_REDO were coupled with model-map inspections in refinement. This approach has highlighted the value of such protocols for efficient analyses. The subunit is constructed from two domains. The N-terminal domain has previously only been observed in cytosolic carboxypeptidase (CCP) proteins. The C-terminal domain, which carries the Zn2+-containing active site, serves to classify this protein as a member of the M14D subfamily of carboxypeptidases. Although eukaryotic CCPs possess deglutamylase activity and are implicated in processing modified tubulin, the function and substrates of the bacterial family members remain unknown. The B. cenocepacia protein did not display deglutamylase activity towards a furylacryloyl glutamate derivative, a potential substrate. Residues previously shown to coordinate the divalent cation and that contribute to peptide-bond cleavage in related enzymes such as bovine carboxypeptidase are conserved. The location of a conserved basic patch in the active site adjacent to the catalytic Zn2+, where an acetate ion is identified, suggests recognition of the carboxy-terminus in a similar fashion to other carboxypeptidases. However, there are significant differences that indicate the recognition of substrates with different properties. Of note is the presence of a lysine in the S1' recognition subsite that suggests specificity towards an acidic substrate.

The novel structure of a cytosolic M14 metallocarboxypeptidase (CCP) from Pseudomonas aeruginosa: a model for mammalian CCPs

PaCCP is a metallocarboxypeptidase (MCP) of the M14 family from Pseudomonas aeruginosa, which belongs to a bacterial clade of carboxypeptidases that are homologous to the recently discovered M14D subfamily of human nonsecretory cytosolic carboxypeptidases (CCPs). CCPs are intracellular peptidases involved, among other roles, in the post-translational modifications of tubulin. Here we report the crystal structure of PaCCP at high resolution (1.6 Å). Its 375 residues are folded in a novel β-sandwich N-terminal domain followed by the classical carboxypeptidase α/β-hydrolase domain, this one in a shorter and more compact form. The former is unique in the whole family and does not have sequential or structural homology with other domains that are usually flanking the latter, like the prodomain of the M14A subfamily or the C-terminal transthyretin/prealbumin-like domains of the M14B subfamily. PaCCP does not display activity against small carboxypeptidase substrates, so in this form it might constitute an inactive precursor of the protease. Structural results derived from cocrystallization with well-known inhibitors of MCPs indicate that the enzyme might only possess C-terminal hydrolase activity against cellular substrates of particular specificity and/or when undergoes structural rearrangements. The derived PaCCP structure allows a first structural insight into the more complex and largely unknown mammalian CCP subfamily.

Crystal structure of bacteriophage ϕNIT1 zinc peptidase PghP that hydrolyzes γ-glutamyl linkage of bacterial poly-γ-glutamate

Poly-γ-glutamate hydrolase P (PghP) of Bacillus subtilis bacteriophage ΦNIT1 hydrolyzes the γ-glutamyl peptide linkage of extracellular poly-γ-glutamate produced by bacilli, which facilitates infection and propagation of phage progenies. Crystal structure of PghP was determined at a resolution of 1.9 Å. Structure of PghP was elucidated as a globular protein with an open α/β mixed core structure and a seven-stranded parallel/anti-parallel β-sheet. The β-sheet contained a core four-stranded parallel β-sheet. A zinc-binding motif, His-Glu-His, was identified at the C-terminal end of the β-sheet. Structure analysis demonstrated that PghP, which had not been previously classified into any peptidase/protease family due to lack of amino acid sequence similarity with known enzymes, had a catalytic center containing a zinc ion and an overall topology resembling mammalian carboxypeptidase A and related enzymes. Structural comparisons indicated important amino acid residues of PghP for catalysis and recognition of the γ-peptide bond of poly-γ-glutamate, which was confirmed by site-directed mutagenesis of PghP.

pH-Dependent reactivity for glycyl-L-tyrosine in carboxypeptidase-A-catalyzed hydrolysis

The dipeptide glycyl-L-tyrosine (GY) can be either a substrate for carboxypeptidase A (CPA) or an inhibitor, depending on pH. In this work, we investigate the pH-dependent reactivity of this dipeptide in CPA-catalyzed hydrolysis using a combined quantum mechanical and molecular mechanical method. It is shown that the monoionic form of the dipeptide, prevalent at high pH, chelates the active site zinc ion, rendering the enzyme inactive. This inhibitory form is consistent with an earlier X-ray structure of the CPA-GY complex. On the other hand, the prevailing di-ionic form of the dipeptide at low pH was found to undergo hydrolysis via a nucleophilic mechanism, leading to an acyl-enzyme complex. The stability of this reaction intermediate is consistent with previous low-temperature solid-state NMR results. The calculated overall free-energy barrier of 20.1 kcal/mol is in excellent agreement with the experimental value of 19.9 kcal/mol.

Catalysis of carboxypeptidase A: promoted-water versus nucleophilic pathways

The catalytic mechanism of carboxypeptidase A (CPA) for the hydrolysis of ester substrates is investigated using hybrid quantum mechanical/molecular mechanical (QM/MM) methods and high-level density functional theory. The prevailing mechanism was found to utilize an active-site water molecule assisted by Glu270, and this so-called promoted-water pathway is similar to that in the CPA catalyzed proteolytic reaction (D. Xu and H. Guo, J. Am. Chem. Soc. 2009, 131, 9780). On the other hand, our simulations indicated the existence of an alternative pathway due to direct nucleophilic attack of Glu270 on the scissile carbonyl carbon. This so-called nucleophilic pathway, which is not viable in proteolytic reactions, leads to a stable acyl-enzyme complex. However, the nucleophilic pathway is nonproductive as it is blocked by a high barrier in the deacylation step. On the basis of results reported here and in our earlier publication, a unified model is proposed to account for nearly all experimental observations concerning the catalysis of CPA.

Progress in metallocarboxypeptidases and their small molecular weight inhibitors

In what corresponds to a life span, metallocarboxypeptidases (MCPs) have jumped from being mere contaminants in animal pancreas powders (in depression year 1929) to be key players in cellular and molecular processes (in yet-another-depression years 2009-2010). MCPs are unique zinc-dependent enzymes that catalyze the breakdown of the amide bond at the C-terminus of peptide and protein substrates and participate in the recovery of dietary amino acids, tissue organogenesis, neurohormone and cytokine maturation and other important physiological processes. More than 26 genes code for MCPs in the human genome, many of them still waiting to be fully understood in terms of physiological function. A variety of MCPs have been linked to diseases in man: acute pancreatitis and pancreas cancer, type 2 diabetes, Alzheimer's Disease, various types of cancer, and fibrinolysis and inflammation. Many of these discoveries have been made possible thanks to recent advances, as exemplified by plasma carboxypeptidases N and B, known for fifty and twenty years, respectively, which have had their structures released only very recently. Plasma carboxypeptidase B is a biological target for therapy because of its involvement in the coagulation/fibrinolysis processes. Besides, the widespread use of carboxypeptidase A as a benchmark metalloprotease since the early days of Biochemistry has allowed the identification and design of an increasingly vast repertory of small molecular weight inhibitors. With these two examples we wish to emphasize that MCPs have become part of the drug discovery portfolio of pharmaceutical companies and academic research laboratories. This paper will review key developments in the discovery and design of MCP small molecular weight inhibitors, with an emphasis on the discovery of chemically diverse entities. Although encouraging advances have been achieved in the last few years, the specificity and oral bioavailability of the new chemotherapeutic agents seem to pose a challenge to medicinal chemists.

Quantum mechanical/molecular mechanical and density functional theory studies of a prototypical zinc peptidase (carboxypeptidase A) suggest a general acid-general base mechanism

Carboxypeptidase A is a zinc-containing enzyme that cleaves the C-terminal residue in a polypeptide substrate. Despite much experimental work, there is still a significant controversy concerning its catalytic mechanism. In this study, the carboxypeptidase A-catalyzed hydrolysis of the hippuryl-L-Phe molecule (k(cat) = 17.7 +/- 0.7 s(-1)) is investigated using both density functional theory and a hybrid quantum mechanical/molecular mechanical approach. The enzymatic reaction was found to proceed via a promoted-water pathway with Glu270 serving as the general base and general acid. Free-energy calculations indicate that the first nucleophilic addition step is rate-limiting, with a barrier of 17.9 kcal/mol. Besides activating the zinc-bound water nucleophile, the zinc cofactor also serves as an electrophilic catalyst that stabilizes the substrate carbonyl oxygen during the formation of the tetrahedral intermediate. In the Michaelis complex, Arg127, rather than Zn(II), is responsible for the polarization of the substrate carbonyl and it also serves as the oxyanion hole. As a result, its mutation leads to a higher free-energy barrier, in agreement with experimental observations.

APIF: a new interaction fingerprint based on atom pairs and its application to virtual screening

A new interaction fingerprint (IF) called APIF (atom-pairs-based interaction fingerprint) has been developed for postprocessing protein-ligand docking results. Unlike other existing fingerprints which employ absolute locations of individual interactions, APIF considers the relative positions of pairs of interacting atoms. Docking-based virtual screening was performed with GOLD using the crystal structures of trypsin, rhinovirus, HIV protease, carboxypeptidase, and estrogen receptor-alpha as targets. A score derived from the similarity of the bit strings for each docking solution to that of a known reference binding mode was obtained. Comparisons between APIF, GoldScore function, and standard interaction fingerprint (CHIF) scores were performed using enrichment plots. Superior recovery rates were observed in the IF score cases. Comparable results were achieved by using either of the two interaction fingerprints, substantially improving GoldScore function enrichment factors. Binding mode analyses were also carried out in order to study the best method for selecting conformations with a binding mode similar to that of the reference crystallized complex. These showed that the first conformations retrieved by interaction fingerprint scores had a more similar binding mode to the reference complex than those retrieved by the GoldScore function.

Characterization of alpha-nitromethyl ketone as a new zinc-binding group based on structural analysis of its complex with carboxypeptidase A

Zinc-binding groups (ZBGs) are exhaustively applied in the development of the new inhibitors against a wide variety of physiologically and pathologically important zinc proteases. Here the alpha-nitro ketone was presented as a new ZBG, which is a transition-state analog featured by the unique bifurcated hydrogen bonds at the active site of carboxypeptidase A based on the structural analysis. Introduction of a nitro group at the alpha-position of the ketone could provide more non-covalent interactions without loss of the abilities to form a tetrahedral transition-state analog.

How similar are enzyme active site geometries derived from quantum mechanical theozymes to crystal structures of enzyme-inhibitor complexes? Implications for enzyme design

Quantum mechanical optimizations of theoretical enzymes (theozymes), which are predicted catalytic arrays of biological functionalities stabilizing a transition state, have been carried out for a set of nine diverse enzyme active sites. For each enzyme, the theozyme for the rate-determining transition state plus the catalytic groups modeled by side-chain mimics was optimized using B3LYP/6-31G(d) or, in one case, HF/3-21G(d) quantum mechanical calculations. To determine if the theozyme can reproduce the natural evolutionary catalytic geometry, the positions of optimized catalytic atoms, i.e., covalent, partial covalent, or stabilizing interactions with transition state atoms, are compared to the positions of the atoms in the X-ray crystal structure with a bound inhibitor. These structure comparisons are contrasted to computed substrate-active site structures surrounded by the same theozyme residues. The theozyme/transition structure is shown to predict geometries of active sites with an average RMSD of 0.64 A from the crystal structure, while the RMSD for the bound intermediate complexes are significantly higher at 1.42 A. The implications for computational enzyme design are discussed.

Catalytic antibodies: contributions from engineering and expression in Escherichia coli

Antibodies have been raised against the transition state of many reactions and shown to catalyse the relevant reaction. Their moderate catalytic efficiencies can be increased by protein engineering, if ways can be found to express the engineered antibody. We have developed a system by which fully functional Fv and Fab fragments can be expressed in Escherichia coli. The Fv fragment dissociates at low concentrations; we therefore devised methods to stabilize the fragment. We showed that the Fv fragment of the antibody McPC603, a phosphorylcholine-binding immunoglobulin A, binds the antigen with the same affinity as does the intact antibody isolated from mouse ascites. Phosphorylcholine is an analogue of the transition state for the hydrolysis of choline carboxylate ester. The Fv fragment of McPC603 catalysed this hydrolysis. Mutational analysis of the residues in the binding site of the antibody has shown which are essential for binding and for catalysis, and the importance of charged residues in certain positions. The E. coli expression system combined with protein engineering and screening methods will facilitate understanding of enzyme catalysis and the development of new catalytic antibodies.

Mutational analysis of aspartoacylase: implications for Canavan disease

Mutations that result in near undetectable activity of aspartoacylase, which catalyzes the deacetylation of N-acetyl-l-aspartate, correlate with Canavan Disease, a neurodegenerative disorder usually fatal during childhood. The underlying biochemical mechanisms of how these mutations ablate activity are poorly understood. Therefore, we developed and tested a three-dimensional homology model of aspartoacylase based on zinc dependent carboxypeptidase A. Mutations of the putative zinc-binding residues (H21G, E24D/G, and H116G), the general proton donor (E178A), and mutants designed to switch the order of the zinc-binding residues (H21E/E24H and E24H/H116E) yielded wild-type aspartoacylase protein levels and undetectable ASPA activity. Mutations that affect substrate carboxyl binding (R71N) and transition state stabilization (R63N) also yielded wild-type aspartoacylase protein levels and undetectable aspartoacylase activity. Alanine substitutions of Cys124 and Cys152, residues indicated by homology modeling to be in close proximity and in the proper orientation for disulfide bonding, yielded reduced ASPA protein and activity levels. Finally, expression of several previously tested (E24G, D68A, C152W, E214X, D249V, E285A, and A305E) and untested (H21P, A57T, I143T, P183H, M195R, K213E/G274R, G274R, and F295S) Canavan Disease mutations resulted in undetectable enzyme activity, and only E285A and P183H showed wild-type aspartoacylase protein levels. These results show that aspartoacylase is a member of the caboxypeptidase A family and offer novel explanations for most loss-of-function aspartoacylase mutations associated with Canavan Disease.

Nna1-like proteins are active metallocarboxypeptidases of a new and diverse M14 subfamily

Nna1 has some sequence similarity to metallocarboxypeptidases, but the biochemical characterization of Nna1 has not previously been reported. In this work we performed a detailed genomic scan and found >100 Nna1 homologues in bacteria, Protista, and Animalia, including several paralogs in most eukaryotic species. Phylogenetic analysis of the Nna1-like sequences demonstrates a major divergence between Nna1-like peptidases and the previously known metallocarboxypeptidases subfamilies: M14A, M14B, and M14C. Conformational modeling of representative Nna1-like proteins from a variety of species indicates an unusually open active site, a property that might facilitate its action on a wide variety of peptide and protein substrates. To test this, we expressed a recombinant form of one of the Nna1-like peptidases from Caenorhabditis elegans and demonstrated that this protein is a fully functional metallocarboxypeptidase that cleaves a range of C-terminal amino acids from synthetic peptides. The enzymatic activity is activated by ATP/ADP and salt-inactivated, and is preferentially inhibited by Z-Glu-Tyr dipeptide, which is without precedent in metallocarboxypeptidases and resembles tubulin carboxypeptidase functioning; this hypothesis is strongly reinforced by the results depicted in Kalinina et al. published as accompanying paper in this journal. Our findings demonstrate that the M14 family of metallocarboxypeptidases is more complex and diverse than expected, and that Nna1-like peptidases are functional variants of such enzymes, representing a novel subfamily (we propose the name M14D) that contributes substantially to such diversity.

Structure of aspartoacylase, the brain enzyme impaired in Canavan disease

Aspartoacylase catalyzes hydrolysis of N-acetyl-l-aspartate to aspartate and acetate in the vertebrate brain. Deficiency in this activity leads to spongiform degeneration of the white matter of the brain and is the established cause of Canavan disease, a fatal progressive leukodystrophy affecting young children. We present crystal structures of recombinant human and rat aspartoacylase refined to 2.8- and 1.8-A resolution, respectively. The structures revealed that the N-terminal domain of aspartoacylase adopts a protein fold similar to that of zinc-dependent hydrolases related to carboxypeptidases A. The catalytic site of aspartoacylase shows close structural similarity to those of carboxypeptidases despite only 10-13% sequence identity between these proteins. About 100 C-terminal residues of aspartoacylase form a globular domain with a two-stranded beta-sheet linker that wraps around the N-terminal domain. The long channel leading to the active site is formed by the interface of the N- and C-terminal domains. The C-terminal domain is positioned in a way that prevents productive binding of polypeptides in the active site. The structures revealed that residues 158-164 may undergo a conformational change that results in opening and partial closing of the channel entrance. We hypothesize that the catalytic mechanism of aspartoacylase is closely analogous to that of carboxypeptidases. We identify residues involved in zinc coordination, and propose which residues may be involved in substrate binding and catalysis. The structures also provide a structural framework necessary for understanding the deleterious effects of many missense mutations of human aspartoacylase.

The crystal structure of the bacteriophage PSA endolysin reveals a unique fold responsible for specific recognition of Listeria cell walls

Bacteriophage murein hydrolases exhibit high specificity towards the cell walls of their host bacteria. This specificity is mostly provided by a structurally well defined cell wall-binding domain that attaches the enzyme to its solid substrate. To gain deeper insight into this mechanism we have crystallized the complete 314 amino acid endolysin from the temperate Listeria monocytogenes phage PSA. The crystal structure of PlyPSA was determined by single wavelength anomalous dispersion methods and refined to 1.8 A resolution. The two functional domains of the polypeptide, providing cell wall-binding and enzymatic activities, can be clearly distinguished and are connected via a linker segment of six amino acid residues. The core of the N-acetylmuramoyl-L-alanine amidase moiety is formed by a twisted, six-stranded beta-sheet flanked by six helices. Although the catalytic domain is unique among the known Listeria phage endolysins, its structure is highly similar to known phosphorylase/hydrolase-like alpha/beta-proteins, including an autolysin amidase from Paenibacillus polymyxa. In contrast, the C-terminal domain of PlyPSA features a novel fold, comprising two copies of a beta-barrel-like motif, which are held together by means of swapped beta-strands. The architecture of the enzyme with its two separate domains explains its unique substrate recognition properties and also provides insight into the lytic mechanisms of related Listeria phage endolysins, a class of enzymes that bear biotechnological potential.

Mapping the energetics of water-protein and water-ligand interactions with the "natural" HINT forcefield: predictive tools for characterizing the roles of water in biomolecules

The energetics and hydrogen bonding pattern of water molecules bound to proteins were mapped by analyzing structural data (resolution better than 2.3A) for sets of uncomplexed and ligand-complexed proteins. Water-protein and water-ligand interactions were evaluated using hydropatic interactions (HINT), a non-Newtonian forcefield based on experimentally determined logP(octanol/water) values. Potential water hydrogen bonding ability was assessed by a new Rank algorithm. The HINT-derived binding energies and Ranks for second shell water molecules were -0.04 kcal mol(-1) and 0.0, respectively, for first shell water molecules -0.38 kcal mol(-1) and 1.6, for active site water molecules -0.45 kcal mol(-1) and 2.3, for cavity water molecules -0.55 kcal mol(-1) and 3.3, and for buried water molecules -0.56 kcal mol(-1) and 4.4. For the last four classes, similar energies indicate that internal and external water molecules interact with protein almost equally, despite different degrees of hydrogen bonding. The binding energies and Ranks for water molecules bridging ligand-protein were -1.13 kcal mol(-1) and 4.5, respectively. This energetic contribution is shared equally between protein and ligand, whereas Rank favors the protein. Lastly, by comparing the uncomplexed and complexed forms of proteins, guidelines were developed for prediction of the roles played by active site water molecules in ligand binding. A water molecule with high Rank and HINT score is unlikely to make further interactions with the ligand and is largely irrelevant to the binding process, while a water molecule with moderate Rank and high HINT score is available for ligand interaction. Water molecule displaced for steric reasons were characterized by lower Rank and HINT score. These guidelines, tested by calculating HINT score and Rank for 50 water molecules bound in the active site of four uncomplexed proteins (for which the structures of the liganded forms were also available), correctly predicted the ultimate roles (in the complex) for 76% of water molecules. Some failures were likely due to ambiguities in the structural data.

Structural basis of the resistance of an insect carboxypeptidase to plant protease inhibitors

Corn earworm (Helicoverpa zea), also called tomato fruitworm, is a common pest of many Solanaceous plants. This insect is known to adapt to the ingestion of plant serine protease inhibitors by using digestive proteases that are insensitive to inhibition. We have now identified a B-type carboxypeptidase of H. zea (CPBHz) insensitive to potato carboxypeptidase inhibitor (PCI) in corn earworm. To elucidate the structural features leading to the adaptation of the insect enzyme, the crystal structure of the recombinant CPBHz protein was determined by x-ray diffraction. CPBHz is a member of the A/B subfamily of metallocarboxypeptidases, which displays the characteristic metallocarboxypeptidase alpha/beta-hydrolase fold, and does not differ essentially from the previously described Helicoverpa armigera CPA, which is very sensitive to PCI. The data provide structural insight into several functional properties of CPBHz. The high selectivity shown by CPBHz for C-terminal lysine residues is due to residue changes in the S1' substrate specificity pocket that render it unable to accommodate the side chain of an arginine. The insensitivity of CPBHz to plant inhibitors is explained by the exceptional positioning of two of the main regions that stabilize other carboxypeptidase-PCI complexes, the beta8-alpha9 loop, and alpha7 together with the alpha7-alpha8 loop. The rearrangement of these two regions leads to a displacement of the active-site entrance that impairs the proper interaction with PCI. This report explains a crystal structure of an insect protease and its adaptation to defensive plant protease inhibitors.

Identification of molting fluid carboxypeptidase A (MF-CPA) in Bombyx mori

Using microarray analyses, we identified carboxypeptidase A (MF-CPA), which was induced during pupal ecdysis in the wing discs of Bombyx mori. Here, we report the functional characterization of MF-CPA. MF-CPA has amino acid sequence similarities with the proteins in the carboxypeptidase A/B subfamily, from human to nematode. The MF-CPA gene is expressed during the molting periods in the epithelial tissues. MF-CPA is detected in the molting fluid, which fills the space between the old and new cuticle during molting. By Western blot analysis, we show that MF-CPA is secreted as a zymogen and processed in the molting fluid. Recombinant MF-CPA expressed in the insect cells has carboxypeptidase A activity. We propose that MF-CPA degrades the proteins from the old cuticle during the molting periods and contributes to recycling of the amino acids.

The role of water in protein-DNA recognition

Is it by design or by default that water molecules are observed at the interfaces of some protein-DNA complexes? Both experimental and theoretical studies on the thermodynamics of protein-DNA binding overwhelmingly support the extended hydrophobic view that water release from interfaces favors binding. Structural and energy analyses indicate that the waters that remain at the interfaces of protein-DNA complexes ensure liquid-state packing densities, screen the electrostatic repulsions between like charges (which seems to be by design), and in a few cases act as linkers between complementary charges on the biomolecules (which may well be by default). This review presents a survey of the current literature on water in protein-DNA complexes and a critique of various interpretations of the data in the context of the role of water in protein-DNA binding and principles of protein-DNA recognition in general.

A practical approach to docking of zinc metalloproteinase inhibitors

Forty zinc-dependent metalloproteinase/ligand complexes with known crystal structures were re-docked using five docking/scoring approaches (DOCK, FlexX, DrugScore, GOLD, and AutoDock). Correct geometry of the coordination bonds between the ligand's zinc binding group (ZBG) and the catalytic zinc is important for docking accuracy and scoring reliability. More than 75% of docked poses with RMSD less than 2A were found to have appropriate ZBG binding, but for poor ZBG binding, about 95% of poses failed to dock correctly. Elimination of poses with inappropriate zinc binding resulted in better binding energy predictions that were further improved by dividing the ligands into subsets according to the ZBG (carboxylates, hydroxamates, and phosphorus containing groups). After a subset re-scoring using the regression functions obtained for individual subsets, DrugScore was able to explain 77% and the consensus scoring scheme X-CSCORE even 88% of variance in binding energies. The approach combining ZBG-based pose selection and subset re-scoring improved the hit rate in virtual screening for metalloproteinase inhibitors for all tested methods by 4-16%.

Single-strand-specific nucleases

Single-strand-specific nucleases are multifunctional enzymes and widespread in distribution. Their ability to act selectively on single-stranded nucleic acids and single-stranded regions in double-stranded nucleic acids has led to their extensive application as probes for the structural determination of nucleic acids. Intracellularly, they have been implicated in recombination, repair and replication, whereas extracellular enzymes have a role in nutrition. Although more than 30 single-strand-specific nucleases from various sources have been isolated till now, only a few enzymes (S1 nuclease from Aspergillus oryzae, P1 nuclease from Penicillium citrinum and nucleases from Alteromonas espejiana, Neurospora crassa, Ustilago maydis and mung bean) have been characterized to a significant extent. Recently, some of these enzymes have been cloned, their crystal structures solved and their interactions with different substrates have been established. The detection, purification, characteristics, structure-function correlations, biological role and applications of single-strand-specific nucleases are reviewed.

Loopy proteins appear conserved in evolution

Over the last decade, structural biologists have unravelled many proteins that appear natively disordered. Common assumptions are that many of these proteins adopt structure through binding and that the structural flexibility enables them to adopt different functions. Here, we investigated regions of more than 70 sequence-consecutive residues that have no regular secondary structure (NORS). Analysing 31 entirely sequenced organisms, we predicted five times as many proteins with NORS regions (loopy proteins) in eukaryotes (20%) than in prokaryotes and archaeas (4%). Thousands of these NORS regions were over 150 residues long. The amino acid composition of NORS regions differed from that of loops in PDB. Although NORS proteins had significantly more residues in low-complexity regions than other proteins, simple cut-off thresholds for sequence bias missed most NORS regions. On average, NORS regions were evolutionarily at least as conserved as their flanking regions. Furthermore, yeast proteins with NORS regions had more protein-protein interaction partners than other proteins. Regulatory and transcription-related functions were over-represented in loopy proteins, biosynthesis and energy metabolism were under-represented. Overall, our analysis confirmed that proteins with non-regular structures appear to play important functional roles, and they may adopt as yet unknown types of protein structures.

Do water molecules mediate protein-DNA recognition?

A comprehensive analysis of interfacial water molecules in the structures of 109 unique protein-DNA complexes is presented together with a new view on their role in protein-DNA recognition. Location of interfacial water molecules as reported in the crystal structures and as emerging from a series of molecular dynamics studies on protein-DNA complexes with explicit solvent and counterions, was analyzed based on their acceptor, donor hydrogen bond relationships with the atoms and residues of the macromolecules, electrostatic field calculations and packing density considerations. Water molecules for the purpose of this study have been categorized into four classes: viz. (I) those that contact both the protein and the DNA simultaneously and thus mediate recognition directly; (II) those that contact either the protein or the DNA exclusively via hydrogen bonds solvating each solute separately; (III) those that contact the hydrophobic groups in either the protein or the DNA; and, lastly (IV) those that contact another water molecule. Of the 17,963 crystallographic water molecules under examination, about 6% belong to class I and 76% belong to class II. About three-fourths of class I and class II water molecules are exclusively associated with hydrogen bond acceptor atoms of both protein and DNA. Noting that DNA is polyanionic, it is significant that a majority of the crystallographically observed water molecules as well as those from molecular dynamics simulations should be involved in facilitating binding by screening unfavorable electrostatics. Less than 2% of the reported water molecules occur between hydrogen bond donor atoms of protein and acceptor atoms of DNA. These represent cases where protein atoms cannot reach out to DNA to make favorable hydrogen bond interactions due to packing/structural restrictions and interfacial water molecules provide an extension to side-chains to accomplish hydrogen bonding.

Crystal structure of a novel mid-gut procarboxypeptidase from the cotton pest Helicoverpa armigera

The cotton bollworm Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) is one of the most serious insect pests in Australia, India and China. The larva causes substantial economical losses to legume, fibre, cereal oilseed and vegetable crops. This pest has proven to be difficult to control by conventional means, mainly due to the development of pesticide resistance. We present here the 2.5 A crystal structure from the novel procarboxypeptidase (PCPAHa) found in the gut extracts from H. armigera larvae, the first one reported for an insect. This metalloprotease is synthesized as a zymogen of 46.6 kDa which, upon in vitro activation with Lys-C endoproteinase, yields a pro-segment of 91 residues and an active carboxypeptidase moiety of 318 residues. Both regions show a three-dimensional structure quite similar to the corresponding structures in mammalian digestive carboxypeptidases, the most relevant structural differences being located in the loops between conserved secondary structure elements, including the primary activation site. This activation site contains the motif (Ala)(5)Lys at the C terminus of the helix connecting the pro- and the carboxypeptidase domains. A remarkable feature of PCPAHa is the occurrence of the same (Ala)(6)Lys near the C terminus of the active enzyme. The presence of Ser255 in PCPAHa instead of Ile and Asp found in the pancreatic A and B forms, respectively, enlarges the S1' specificity pocket and influences the substrate preferences of the enzyme. The C-terminal tail of the leech carboxypeptidase inhibitor has been modelled into the PCPAHa active site to explore the substrate preferences and the enzymatic mechanism of this enzyme.

Crystal structure of the stromelysin-3 (MMP-11) catalytic domain complexed with a phosphinic inhibitor mimicking the transition-state

Stromelysin-3 (ST3) is a matrix metalloproteinase (MMP-11) whose proteolytic activity plays an important role in tumorigenicity enhancement. In breast cancer, ST3 is a bad prognosis marker: its expression is associated with a poor clinical outcome. This enzyme therefore represents an attractive therapeutic target. The topology of matrix metalloproteinases (MMPs) is remarkably well conserved, making the design of highly specific inhibitors difficult. The major difference between MMPs lies in the S(1)' subsite, a well-defined hydrophobic pocket of variable depth. The present crystal structure, the first 3D-structure of the ST3 catalytic domain in interaction with a phosphinic inhibitor mimicking a (d, l) peptide, clearly demonstrates that its S(1)' pocket corresponds to a tunnel running through the enzyme. This open channel is filled by the inhibitor P(1)' group which adopts a constrained conformation to fit this pocket, together with two water molecules interacting with the ST3-specific residue Gln215. These observations provide clues for the design of more specific inhibitors and show how ST3 can accommodate a phosphinic inhibitor mimicking a (d, l) peptide. The presence of a water molecule interacting with one oxygen atom of the inhibitor phosphinyl group and the proline residue of the Met-turn suggests how the intermediate formed during proteolysis may be stabilized. Furthermore, the hydrogen bond distance observed between the methyl of the phosphinic group and the carbonyl group of Ala182 mimics the interaction between this carbonyl group and the amide group of the cleaved peptidic bond. Our crystal structure provides a good model to study the MMPs mechanism of proteolysis.

Thrombin-activatable fibrinolysis inhibitor (TAFI, plasma procarboxypeptidase B, procarboxypeptidase R, procarboxypeptidase U)

Recently, a new inhibitor of fibrinolysis was described. This inhibitor downregulated fibrinolysis after it was activated by thrombin, and was therefore named TAFI (thrombin-activatable fibrinolysis inhibitor; EC 3.4.17.20). TAFI turned out to be identical to previously described proteins, procarboxypeptidase U, procarboxypeptidase R, and plasma procarboxypeptidase B. In this overview, the protein will be referred to as TAFI. TAFI is a procarboxypeptidase and a member of the family of metallocarboxypeptidases. These enzymes are circulating in plasma and are present in several tissues such as pancreas. In this review, we will describe the properties of basic carboxypeptidases with the emphasis on the role of TAFI in coagulation and fibrinolysis. It cannot be ruled out, however, that TAFI has other, yet undefined, functions in biology.

Novel bifunctional hyperthermostable carboxypeptidase/aminoacylase from Pyrococcus horikoshii OT3

Genome sequencing of the thermophilic archaeon Pyrococcus horikoshii OT3 revealed a gene which had high sequence similarity to the gene encoding the carboxypeptidase of Sulfolobus solfataricus and also to that encoding the aminoacylase from Bacillus stearothermophilus. The gene from P. horikoshii comprises an open reading frame of 1,164 bp with an ATG initiation codon and a TGA termination codon, encoding a 43,058-Da protein of 387 amino acid residues. However, some of the proposed active-site residues for carboxypeptidase were not found in this gene. The gene was overexpressed in Escherichia coli with the pET vector system, and the expressed enzyme had high hydrolytic activity for both carboxypeptidase and aminoacylase at high temperatures. The enzyme was stable at 90 degrees C, with the highest activity above 95 degrees C. The enzyme contained one bound zinc ion per one molecule that was essential for the activity. The results of site-directed mutagenesis of Glu367, which corresponds to the essential Glu270 in bovine carboxypeptidase A and the essential Glu in other known carboxypeptidases, revealed that Glu367 was not essential for this enzyme. The results of chemical modification of the SH group and site-directed mutagenesis of Cys102 indicated that Cys102 was located at the active site and was related to the activity. From these findings, it was proven that this enzyme is a hyperthermostable, bifunctional, new zinc-dependent metalloenzyme which is structurally similar to carboxypeptidase but whose hydrolytic mechanism is similar to that of aminoacylase. Some characteristics of this enzyme suggested that carboxypeptidase and aminoacylase might have evolved from a common origin.

Alternate modes of binding in two crystal structures of alkaline phosphatase-inhibitor complexes

Two high resolution crystal structures of Escherichia coli alkaline phosphatase (AP) in the presence of phosphonate inhibitors are reported. The phosphonate compounds, phosphonoacetic acid (PAA) and mercaptomethylphosphonic acid (MMP), bind competitively to AP with dissociation constants of 5.5 and 0.6 mM, respectively. The structures of the complexes of AP with PAA and MMP were refined at high resolution to crystallographic R-values of 19.0 and 17.5%, respectively. Refinement of the AP-inhibitor complexes was carried out using X-PLOR. The final round of refinement was done using SHELXL-97. Crystallographic analyses of the inhibitor complexes reveal different binding modes for the two phosphonate compounds. The significant difference in binding constants can be attributed to these alternative binding modes observed in the high resolution X-ray structures. The phosphinyl group of PAA coordinates to the active site zinc ions in a manner similar to the competitive inhibitor and product inorganic phosphate. In contrast, MMP binds with its phosphonate moiety directed toward solvent. Both enzyme-inhibitor complexes exhibit close contacts, one of which has the chemical and geometrical potential to be considered an unconventional hydrogen bond of the type C-H...X.