Crystallographic analysis of transition-state mimics bound to penicillopepsin: phosphorus-containing peptide analogues

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.

iMOLSDOCK: Induced-fit docking using mutually orthogonal Latin squares (MOLS)

We have earlier reported the MOLSDOCK technique to perform rigid receptor/flexible ligand docking. The method uses the MOLS method, developed in our laboratory. In this paper we report iMOLSDOCK, the 'flexible receptor' extension we have carried out to the algorithm MOLSDOCK. iMOLSDOCK uses mutually orthogonal Latin squares (MOLS) to sample the conformation and the docking pose of the ligand and also the flexible residues of the receptor protein. The method then uses a variant of the mean field technique to analyze the sample to arrive at the optimum. We have benchmarked and validated iMOLSDOCK with a dataset of 44 peptide-protein complexes with peptides. We have also compared iMOLSDOCK with other flexible receptor docking tools GOLD v5.2.1 and AutoDock Vina. The results obtained show that the method works better than these two algorithms, though it consumes more computer time.

Extensive expansion of A1 family aspartic proteinases in fungi revealed by evolutionary analyses of 107 complete eukaryotic proteomes

The A1 family of eukaryotic aspartic proteinases (APs) forms one of the 16 AP families. Although one of the best characterized families, the recent increase in genome sequence data has revealed many fungal AP homologs with novel sequence characteristics. This study was performed to explore the fungal AP sequence space and to obtain an in-depth understanding of fungal AP evolution. Using a comprehensive phylogeny of approximately 700 AP sequences from the complete proteomes of 87 fungi and 20 nonfungal eukaryotes, 11 major clades of APs were defined of which clade I largely corresponds to the A1A subfamily of pepsin-archetype APs. Clade II largely corresponds to the A1B subfamily of nepenthesin-archetype APs. Remarkably, the nine other clades contain only fungal APs, thus indicating that fungal APs have undergone a large sequence diversification. The topology of the tree indicates that fungal APs have been subject to both “birth and death” evolution and “functional redundancy and diversification.” This is substantiated by coclustering of certain functional sequence characteristics. A meta-analysis toward the identification of Cluster Determining Positions (CDPs) was performed in order to investigate the structural and biochemical basis for diversification. Seven CDPs contribute to the secondary structure of the enzyme. Three other CDPs are found in the vicinity of the substrate binding cleft. Tree topology, the large sequence variation among fungal APs, and the apparent functional diversification suggest that an amendment to update the current A1 AP classification based on a comprehensive phylogenetic clustering might contribute to refinement of the classification in the MEROPS peptidase database.

PeptiSite: a structural database of peptide binding sites in 4D

We developed PeptiSite, a comprehensive and reliable database of biologically and structurally characterized peptide-binding sites, in which each site is represented by an ensemble of its complexes with protein, peptide and small molecule partners. The unique features of the database include: (1) the ensemble site representation that provides a fourth dimension to the otherwise three dimensional data, (2) comprehensive characterization of the binding site architecture that may consist of a multimeric protein assembly with cofactors and metal ions and (3) analysis of consensus interaction motifs within the ensembles and identification of conserved determinants of these interactions. Currently the database contains 585 proteins with 650 peptide-binding sites. http://peptisite.ucsd.edu/ link allows searching for the sites of interest and interactive visualization of the ensembles using the ActiveICM web-browser plugin. This structural database for protein-peptide interactions enables understanding of structural principles of these interactions and may assist the development of an efficient peptide docking benchmark.

Structure-activity relationships of the phosphonate antibiotic dehydrophos

Synthetic derivatives of the phosphonate antibiotic dehydrophos were tested for antimicrobial activity. Both the phosphonate monomethyl ester and the vinyl phosphonate moiety proved to be important for bacteriocidal activity of the natural product.

Synthesis of α-carboxyphosphinopeptides derived from norleucine

In the present study, we describe in detail the synthesis of a relatively rare class of phosphorus compounds, α-carboxyphosphinopeptides. We prepared several norleucine-derived α-carboxyphosphinic pseudopeptides of the general formula Nle-Ψ[PO(OH)]-Gly. These compounds could have important applications as transition state-mimicking inhibitors for methionine or leucine aminopeptidases or other enzymes. For the preparation of the key α-carboxyphosphinate protected precursors, we investigated, compared and improved two different synthetic methods described in literature: the Arbuzov reaction of a silylated N-protected phosphinic acid with a bromoacetate ester and the nucleophilic addition of a mixed O-methyl S-phenyl N-protected phosphonic acid or a methyl N-protected phosphonochloridate with tert-butyl lithioacetate. We also prepared two N-Fmoc protected synthons, Fmoc-Nle-Ψ[PO(OH)]-Gly-COOH and Fmoc-Nle-Ψ[PO(OAd)]-Gly-COOH, and demonstrated that these precursors are suitable building blocks for the solid-phase synthesis of α-carboxyphosphinopeptides.

Grassystatins A-C from marine cyanobacteria, potent cathepsin E inhibitors that reduce antigen presentation

In our efforts to explore marine cyanobacteria as a source of novel bioactive compounds, we discovered a statine unit-containing linear decadepsipeptide, grassystatin A (1), which we screened against a diverse set of 59 proteases. We describe the structure determination of 1 and two natural analogues, grassystatins B (2) and C (3), using NMR, MS, and chiral HPLC techniques. Compound 1 selectively inhibited cathepsins D and E with IC(50)s of 26.5 nM and 886 pM, respectively. Compound 2 showed similar potency and selectivity against cathepsins D and E (IC(50)s of 7.27 nM and 354 pM, respectively), whereas the truncated peptide analogue grassystatin C (3), which consists of two fewer residues than 1 and 2, was less potent against both but still selective for cathepsin E. The selectivity of compounds 1-3 for cathepsin E over D (20-38-fold) suggests that these natural products may be useful tools to probe cathepsin E function. We investigated the structural basis of this selectivity using molecular docking. We also show that 1 can reduce antigen presentation by dendritic cells, a process thought to rely on cathepsin E.

Slow-onset inhibition of fumarylacetoacetate hydrolase by phosphinate mimics of the tetrahedral intermediate: kinetics, crystal structure and pharmacokinetics

FAH (fumarylacetoacetate hydrolase) catalyses the final step of tyrosine catabolism to produce fumarate and acetoacetate. HT1 (hereditary tyrosinaemia type 1) results from deficiency of this enzyme. Previously, we prepared a partial mimic of the putative tetrahedral intermediate in the reaction catalysed by FAH co-crystallized with the enzyme to reveal details of the mechanism [Bateman, Bhanumoorthy, Witte, McClard, Grompe and Timm (2001) J. Biol. Chem. 276, 15284-15291]. We have now successfully synthesized complete mimics CEHPOBA {4-[(2-carboxyethyl)-hydroxyphosphinyl]-3-oxobutyrate} and COPHPAA {3-[(3-carboxy-2-oxopropyl)hydroxyphosphinyl]acrylate}, which inhibit FAH in slow-onset tight-binding mode with K(i) values of 41 and 12 nM respectively. A high-resolution (1.35 A; 1 A=0.1 nm) crystal structure of the FAH.CEHPOBA complex was solved to reveal the affinity determinants for these compounds and to provide further insight into the mechanism of FAH catalysis. These compounds are active in vivo, and CEHPOBA demonstrated a notable dose-dependent increase in SA (succinylacetone; a metabolite seen in patients with HT1) in mouse serum after repeated injections, and, following a single injection (1 mumol/g; intraperitoneal), only a modest regain of FAH enzyme activity was detected in liver protein isolates after 24 h. These potent inhibitors provide a means to chemically phenocopy the metabolic defects of either HT1 or FAH knockout mice and promise future pharmacological utility for hepatocyte transplantation.

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.

Kinetic studies on beta-site amyloid precursor protein-cleaving enzyme (BACE). Confirmation of an iso mechanism

The steady-state kinetic mechanism of beta-amyloid precursor protein-cleaving enzyme (BACE)-catalyzed proteolytic cleavage was evaluated using product and statine- (Stat(V)) or hydroxyethylene-containing (OM99-2) peptide inhibition data, solvent kinetic isotope effects, and proton NMR spectroscopy. The noncompetitive inhibition pattern observed for both cleavage products, together with the independence of Stat(V) inhibition on substrate concentration, suggests a uni-bi-iso kinetic mechanism. According to this mechanism, the enzyme undergoes multiple conformation changes during the catalytic cycle. If any of these steps are rate-limiting to turnover, an enzyme form preceding the rate-limiting conformational change should accumulate. An insignificant solvent kinetic isotope effect (SKIE) on k(cat)/K(m), a large inverse solvent kinetic isotope effect on k(cat), and the absence of any SKIE on the inhibition onset by Stat(V) during catalysis together indicate that the rate-limiting iso-step occurs after formation of a tetrahedral intermediate. A moderately short and strong hydrogen bond (at delta 13.0 ppm and phi of 0.6) has been observed by NMR spectroscopy in the enzyme-hydroxyethylene peptide (OM99-2) complex that presumably mimics the tetrahedral intermediate of catalysis. Collapse of this intermediate, involving multiple steps and interconversion of enzyme forms, has been suggested to impose a rate limitation, which is manifested in a significant SKIE on k(cat). Multiple enzyme forms and their distribution during catalysis were evaluated by measuring the SKIE on the noncompetitive (mixed) inhibition constants for the C-terminal reaction product. Large, normal SKIE values were observed for these inhibition constants, suggesting that both kinetic and thermodynamic components contribute to the K(ii) and K(is) expressions, as has been suggested for other iso-mechanism featuring enzymes. We propose that a conformational change related to the reprotonation of aspartates during or after the bond-breaking event is the rate-limiting segment in the catalytic reaction of beta-amyloid precursor protein-cleaving enzyme, and ligands binding to other than the ground-state forms of the enzyme might provide inhibitors of greater pharmacological relevance.

Five atomic resolution structures of endothiapepsin inhibitor complexes: implications for the aspartic proteinase mechanism

Endothiapepsin is derived from the fungus Endothia parasitica and is a member of the aspartic proteinase class of enzymes. This class of enzyme is comprised of two structurally similar lobes, each lobe contributing an aspartic acid residue to form a catalytic dyad that acts to cleave the substrate peptide bond. The three-dimensional structures of endothiapepsin bound to five transition state analogue inhibitors (H189, H256, CP-80,794, PD-129,541 and PD-130,328) have been solved at atomic resolution allowing full anisotropic modelling of each complex. The active sites of the five structures have been studied with a view to studying the catalytic mechanism of the aspartic proteinases by locating the active site protons by carboxyl bond length differences and electron density analysis. In the CP-80,794 structure there is excellent electron density for the hydrogen on the inhibitory statine hydroxyl group which forms a hydrogen bond with the inner oxygen of Asp32. The location of this proton has implications for the catalytic mechanism of the aspartic proteinases as it is consistent with the proposed mechanism in which Asp32 is the negatively charged aspartate. A number of short hydrogen bonds (approximately 2.6 A) with ESD values of around 0.01 A that may have a role in catalysis have been identified within the active site of each structure; the lengths of these bonds have been confirmed using NMR techniques. The possibility and implications of low barrier hydrogen bonds in the active site are considered.

Analysis of crystal structures of aspartic proteinases: on the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes

To elucidate the role of amino acid residues adjacent to the catalytic site of pepsin-like enzymes, we analyzed and compared the crystal structures of these enzymes, their complexes with inhibitors, and zymogens in the active site area (a total of 82 structures). In addition to the water molecule (W1) located between the active carboxyls and playing a role of the nucleophile during catalytic reaction, another water molecule (W2) at the vicinity of the active groups was found to be completely conserved. This water molecule plays an essential role in formation of a chain of hydrogen-bonded residues between the active site flap and the active carboxyls on ligand binding. These data suggest a new approach to understanding the role of residues around the catalytic site, which can assist the development of the catalytic reaction. The influence of groups adjacent to the active carboxyls is manifested by pepsin activity at pH 1.0. Some features of pepsin-like enzymes and their mutants are discussed in the framework of the approach.

Mechanistic inferences from the crystal structure of fumarylacetoacetate hydrolase with a bound phosphorus-based inhibitor

Fumarylacetoacetate hydrolase (FAH) catalyzes the hydrolytic cleavage of a carbon-carbon bond in fumarylacetoacetate to yield fumarate and acetoacetate as the final step of Phe and Tyr degradation. This unusual reaction is an essential human metabolic function, with loss of FAH activity causing the fatal metabolic disease hereditary tyrosinemia type I (HT1). An enzymatic mechanism involving a catalytic metal ion, a Glu/His catalytic dyad, and a charged oxyanion hole was previously proposed based on recently determined FAH crystal structures. Here we report the development and characterization of an FAH inhibitor, 4-(hydroxymethylphosphinoyl)-3-oxo-butanoic acid (HMPOBA), that competes with the physiological substrate with a K(i) of 85 microM. The crystal structure of FAH complexed with HMPOBA refined at 1.3-A resolution reveals the molecular basis for the competitive inhibition, supports the proposed formation of a tetrahedral alkoxy transition state intermediate during the FAH catalyzed reaction, and reveals a Mg(2+) bound in the enzyme's active site. The analysis of FAH structures corresponding to different catalytic states reveals significant active site side-chain motions that may also be related to catalytic function. Thus, these results advance the understanding of an essential catabolic reaction associated with a fatal metabolic disease and provide insight into the structure-based development of FAH inhibitors.

A novel serine protease inhibition motif involving a multi-centered short hydrogen bonding network at the active site

We describe a new serine protease inhibition motif in which binding is mediated by a cluster of very short hydrogen bonds (<2.3 A) at the active site. This protease-inhibitor binding paradigm is observed at high resolution in a large set of crystal structures of trypsin, thrombin, and urokinase-type plasminogen activator (uPA) bound with a series of small molecule inhibitors (2-(2-phenol)indoles and 2-(2-phenol)benzimidazoles). In each complex there are eight enzyme-inhibitor or enzyme-water-inhibitor hydrogen bonds at the active site, three of which are very short. These short hydrogen bonds connect a triangle of oxygen atoms comprising O(gamma)(Ser195), a water molecule co-bound in the oxyanion hole (H(2)O(oxy)), and the phenolate oxygen atom of the inhibitor (O6'). Two of the other hydrogen bonds between the inhibitor and active site of the trypsin and uPA complexes become short in the thrombin counterparts, extending the three-centered short hydrogen-bonding array into a tetrahedral array of atoms (three oxygen and one nitrogen) involved in short hydrogen bonds. In the uPA complexes, the extensive hydrogen-bonding interactions at the active site prevent the inhibitor S1 amidine from forming direct hydrogen bonds with Asp189 because the S1 site is deeper in uPA than in trypsin or thrombin. Ionization equilibria at the active site associated with inhibitor binding are probed through determination and comparison of structures over a wide range of pH (3.5 to 11.4) of thrombin complexes and of trypsin complexes in three different crystal forms. The high-pH trypsin-inhibitor structures suggest that His57 is protonated at pH values as high as 9.5. The pH-dependent inhibition of trypsin, thrombin, uPA and factor Xa by 2-(2-phenol)benzimidazole analogs in which the pK(a) of the phenol group is modulated is shown to be consistent with a binding process involving ionization of both the inhibitor and the enzyme. These data further suggest that the pK(a) of His57 of each protease in the unbound state in solution is about the same, approximately 6.8. By comparing inhibition constants (K(i) values), inhibitor solubilities, inhibitor conformational energies and corresponding structures of short and normal hydrogen bond-mediated complexes, we have estimated the contribution of the short hydrogen bond networks to inhibitor affinity ( approximately 1.7 kcal/mol). The structures and K(i) values associated with the short hydrogen-bonding motif are compared with those corresponding to an alternate, Zn(2+)-mediated inhibition motif at the active site. Structural differences among apo-enzymes, enzyme-inhibitor and enzyme-inhibitor-Zn(2+) complexes are discussed in the context of affinity determinants, selectivity development, and structure-based inhibitor design.

Contribution of peptide bonds to inhibitor-protease binding: crystal structures of the turkey ovomucoid third domain backbone variants OMTKY3-Pro18I and OMTKY3-psi[COO]-Leu18I in complex with Streptomyces griseus proteinase B (SGPB) and the structure of the free inhibitor, OMTKY-3-psi[CH2NH2+]-Asp19I

X-ray crystallography has been used to determine the 3D structures of two complexes between Streptomyces griseus proteinase B (SGPB), a bacterial serine proteinase, and backbone variants of turkey ovomucoid third domain (OMTKY3). The natural P1 residue (Leu18I) has been substituted by a proline residue (OMTKY3-Pro18I) and in the second variant, the peptide bond between Thr17I and Leu18I was replaced by an ester bond (OMTKY3-psi[COO]-Leu18I). Both variants lack the P1 NH group that donates a bifurcated hydrogen bond to the carbonyl O of Ser214 and O(gamma) of the catalytic Ser195, one of the common interactions between serine proteinases and their canonical inhibitors. The SGPB:OMTKY3-Pro18I complex has many structural differences in the vicinity of the S1 pocket when compared with the previously determined structure of SGPB:OMTKY3-Leu18I. The result is a huge difference in the DeltaG degrees of binding (8.3 kcal/mol), only part of which can be attributed to the missing hydrogen bond. In contrast, very little structural difference exists between the complexes of SGPB:OMTKY3-psi[COO]-Leu18I and SGPB:OMTKY3-Leu18I, aside from an ester O replacing the P1 NH group. Therefore, the difference in DeltaG degrees, 1.5 kcal/mol as calculated from the measured equilibrium association constants, can be attributed to the contribution of the P1 NH hydrogen bond toward binding. A crystal structure of OMTKY3 having a reduced peptide bond between P1 Leu18I and P'1 Asp19I, (OMTKY3-psi[CH2NH2+]-Asp19I) has also been determined by X-ray crystallography. This variant has very weak association equilibrium constants with SGPB and with chymotrypsin. The structure of the free inhibitor suggests that the reduced peptide bond has not introduced any major structural changes in the inhibitor. Therefore, its poor ability to inhibit serine proteinases is likely due to the disruptions of the canonical interactions at the oxyanion hole.

Carbonic anhydrase inhibitors; phosphoryl-sulfonamides--a new class of high affinity inhibitors of isozymes I and II

A series of phosphorylated aromatic/heterocyclic sulfonamides with the general formula ArSO2NHPO3H2 have been prepared by condensing ArSO2NH2 with phosphorus pentachloride, followed by controlled hydrolysis in the presence of formic acid. The new derivatives generally act as stronger inhibitors of two carbonic anhydrase (CA) isozymes, CA I and CA II, as compared to the parent unsubstituted sulfonamides from which they were obtained. The inhibition mechanism by this new class of CA inhibitors, as well as structure activity correlations for the series of investigated derivatives, are also discussed.

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.

Reduced amide pseudopeptide analogues of a malaria peptide possess secondary structural elements responsible for induction of functional antibodies which react with native proteins expressed in Plasmodium falciparum erythrocyte stages

A psi[CH2NH] isoster bond was introduced by replacing one peptide bond at a time within the 1513 malaria peptide KEKMV motif to obtain a set of five pseudopeptides. The motif belongs to a Plasmodium falciparum malarial peptide coded 1513, derived from the MSP-1 protein. This high-binding motif included in the 1513 peptide is involved in the attachment of the malarial parasite to human erythrocytes. The novel malaria 1513 psi[CH2NH] surrogates were analyzed using RP-HPLC and MALDI-TOF mass spectrometry techniques. Nuclear magnetic resonance experiments allowed definition of the five pseudopeptide analogues' secondary structural features. Such structures are present in only a very few molecules in the 1513 parent peptide. A molecular model demonstrating the solution of the three-dimensional structure of the 1 513 peptide Pse-437 analogue was constructed on the basis of 1H-NMR spectral parameters. Monoclonal antibodies were generated to the five 1513 malaria peptide pseudopeptide analogues. These antibodies not only recognize the native MSP-1 (195 kDa) and its 83 kDa and 42 kDa proteolytic processing proteins but also different SPf(66)n malaria vaccine batches containing the native sequence. In addition, the mAbs were able to modify the kinetics of Plasmodium falciparum parasites' intraerythrocytic development and their ability to invade new RBCs. The presented evidence suggests that peptide bond-modified peptides could reproduce a transient state in 1513's native sequence and represent useful candidates in the development of a second generation of effective malarial vaccines.

"Strong" hydrogen bonds in chemistry and biology

Hydrogen bonds are a key feature of chemical structure and reactivity. Recently there has been much interest in a special class of hydrogen bonds called "strong" or "low-barrier" and characterized by great strength, short distances, a low or vanishing barrier to hydrogen transfer, and distinctive features in the NMR spectrum. Although the energy of an ordinary hydrogen bond is ca 5 kcal mol-1, the strength of these hydrogen bonds may be > or = 10 kcal mol-1. The properties of these hydrogen bonds have been investigated by many experimental techniques, as well as by calculation and by correlations among those properties. Although it has been proposed that strong, short, low-barrier hydrogen bonds are important in enzymatic reactions, it is concluded that the evidence for them in small molecules and in biomolecules is inconclusive.

Similarities of hydrolytic antibodies revealed by their X-ray structures: a review

Numerous antibodies have been programmed to catalyse the hydrolysis of esters as well as other acyl transfer reactions. They were raised against stable analogues that model the structure of the tetrahedral transition states of these reactions. The three-dimensional structures of four hydrolytic antibodies complexed to their respective phosphonate transition state analogues (TSAs) reveal a similar orientation of hapten relative to the antibody. Analysis of the four combining sites suggests that residues binding the phosphonate TSA stabilise the oxyanion intermediate of the reaction and play a preponderant role in catalysis. Comparison of catalytic antibodies selected from the same hybridoma fusion indicates a high similarity of the motifs that catalyse the hydrolysis of a given substrate.

Empirical scoring functions: I. The development of a fast empirical scoring function to estimate the binding affinity of ligands in receptor complexes

This paper describes the development of a simple empirical scoring function designed to estimate the free energy of binding for a protein-ligand complex when the 3D structure of the complex is known or can be approximated. The function uses simple contact terms to estimate lipophilic and metal-ligand binding contributions, a simple explicit form for hydrogen bonds and a term which penalises flexibility. The coefficients of each term are obtained using a regression based on 82 ligand-receptor complexes for which the binding affinity is known. The function reproduces the binding affinity of the complexes with a cross-validated error of 8.68 kJ/mol. Tests on internal consistency indicate that the coefficients obtained are stable to changes in the composition of the training set. The function is also tested on two test sets containing a further 20 and 10 complexes, respectively. The deficiencies of this type of function are discussed and it is compared to approaches by other workers.

A low energy short hydrogen bond in very high resolution structures of protein receptor--phosphate complexes

A very short hydrogen bond between an Asp and a phosphate is established in two high resolution structures (0.98 and 1.05 A). A mutant complex that changes the Asp to an Asn, which forms a normal hydrogen bond, has a similar free energy of binding to the wild type complex, suggesting that the contribution of the short hydrogen bond is not extraordinarily strong.

The three-dimensional structure at 2.4 A resolution of glycosylated proteinase A from the lysosome-like vacuole of Saccharomyces cerevisiae

The crystal structures of glycosylated native proteinase A, an aspartic proteinase found in the vacuole of Saccharomyces cerevisiae, and its complex with a difluorostatone-containing tripeptide have been determined by molecular replacement to 3.5 A and 2.4 A resolutions, respectively. Superposition of the bound and native forms gave an r.m.s. difference of 0.6 A largely reflecting the poor resolution of the native crystal structure. The secondary and tertiary structures are highly similar to those found in porcine pepsin and lysosomal cathepsin D; superposition of the structure of proteinase A bound to the difluorostatone inhibitor on those of pepsin and cathepsin D gave pairwise r.m.s. differences for C(alpha) atoms of 1.36 A and 0.88 A. Most differences occur in loop regions. Comparison of the structure of the proteinase A-difluorostatone complex with that of endothiapepsin bound with the same inhibitor shows that the conformation and hydrogen bond interactions of the inhibitor in the active site are very similar, even though the enzymes have only 27% sequence identity. Electron density for the crystal structure of the proteinase A complex reveals five residues of the oligosaccharide structure attached to Asn67: Man-(1 --> 2)-alpha-Man-(1 --> 3)-beta-Man-(1 --> 4)-beta-GlcNAc-(1 --> 4)-beta-GlcNAc-Asn-67. The first three residues of the oligosaccharide cover the same region of the protein surface as those of the oligosaccharide attached to the equivalent position in cathepsin D. The second carbohydrate attachment site is disordered beyond the first carbohydrate residue in both enzymes.

Expanding the 43C9 class of catalytic antibodies using a chain-shuffling approach

We employed a chain-shuffling technique to determine if the light chain of the catalytic antibody, 43C9, provides the best partner for the 43C9 heavy chain. Previously, we reported construction and screening of a 43C9 HC CROSS library, where the 43C9 heavy-chain gene was crossed with a library of light-chain genes in a lambda bacteriophage system. The library contained a high frequency of reconstituted antibodies recognizing the transition-state analogue. Here, we report the isolation and characterization of four of these clones. Recovered light-chain proteins share 92-96% sequence identity to the 43C9 light-chain protein. Somatic mutations of these light chains occur randomly at positions distant from the active site. Residues required for binding and catalysis were conserved. Mutations affected the topology of the binding site. Nevertheless, catalysis was not affected. Isolation of these light chains suggests the best partner for the 43C9 heavy chain is the original light chain. These clones attempt to broaden a class of 43C9-like antibodies, where the catalytic residues, His91 and Arg96, have been reproducibly selected. Similar catalytic properties between the 43C9-like antibodies suggests binding has been optimized, thus further maturation of the light chain would not lead to a better catalyst. To improve catalysis, other approaches must be considered.

Molecular docking using surface complementarity

A method is described to dock a ligand into a binding site in a protein on the basis of the complementarity of the intermolecular atomic contacts. Docking is performed by maximization of a complementarity function that is dependent on atomic contact surface area and the chemical properties of the contacting atoms. The generality and simplicity of the complementarity function ensure that a wide range of chemical structures can be handled. The ligand and the protein are treated as rigid bodies, but displacement of a small number of residues lining the ligand binding site can be taken into account. The method can assist in the design of improved ligands by indicating what changes in complementarity may occur as a result of the substitution of an atom in the ligand. The capabilities of the method are demonstrated by application to 14 protein-ligand complexes of known crystal structure.

Refined 1.7 A X-ray crystallographic structure of P-30 protein, an amphibian ribonuclease with anti-tumor activity

The X-ray crystallographic structure of P-30 protein (Onconase) has been solved by multiple isomorphous replacement and the structure has been refined at 1.7 A resolution to a conventional R-factor of 0.178. The molecular model comprises all 826 non-hydrogen protein atoms, 96 solvent molecules and a sulfate anion that is bound at the active site. The molecular structure is similar to that of ribonuclease A. The active site cleft is located at the junction of two three-stranded beta-sheets and the N-terminal helix. A sulfate anion is non-covalently bound by Lys9, His10, His97, Phe98 and an intermolecular contact involving Lys55' from a neighboring molecule. The N-terminal pyroglutamyl (Pyr) residue is part of the active site and its O epsilon 1 atom forms a hydrogen bond with the Lys9 N zeta. The previously constructed comparative molecular model of P-30 based on ribonuclease A correctly predicted the overall fold of P-30 and the conformation of its active site residues. The model failed to predict the conformation of Pyr1 and the conformation of the two loops following helix alpha 3 and strand beta 3.

Use of steady state kinetic methods to elucidate the kinetic and chemical mechanisms of retroviral proteases

Despite the current plethora of structural data of HIV-1 protease and the availability of potent inhibitors, whose structures are based in part on the presumed mechanism of action of this enzyme, our actual understanding of its chemical mechanism has been until now based largely on the precedents of the mammalian and fungal aspartic proteases and static three-dimensional data. The available steady state kinetic data of the protease, as reviewed here, constitute a first step in a detailed description of the mechanism of the enzyme to complement the structural data.

Use of nitrogen-15 kinetic isotope effects to elucidate details of the chemical mechanism of human immunodeficiency virus 1 protease

We have used 15N kinetic isotope effects of the HIV-1 protease-catalyzed peptidolysis of Ac-Ser-Gln-Asn-Tyr-Pro-Val-Val-NH2 to characterize the chemical mechanism of this enzyme. In addition, the multiple isotope effects have been determined by measuring the 15N kinetic isotope effects in both H2O and D2O. The isotope effects, measured on values of V/K, were determined by the incorporation of a radiolabel (tritium and 14C in peptides bearing the heavy and light isotopes, respectively) at a position remote from the isotopically labeled scissile peptide bond, such that the isotope effect was determined by measurement of the change in the 14C/3H ratio in recovered substrates at various fractions of reaction. At pH = 6.0 (37 degrees C), the nitrogen isotope effects were slightly, but significantly, inverse in both solvents: 15(V/K)H2O = 0.995 +/- 0.002, and 15(V/K)D2O = 0.992 +/- 0.003. The observation of an inverse nitrogen kinetic isotope effect implies that bonding to the nitrogen atom is becoming stiffened in a reaction transition state, and since this inverse isotope effect is enhanced in D2O, this isotope effect likely arises from protonation of the proline nitrogen atom.(ABSTRACT TRUNCATED AT 250 WORDS)