Perturbations in the free energy and enthalpy of ionization of histidine-159 at the active site of papain as determined by fluorescence spectroscopy

Thiosemicarbazone derivatives: Evaluation as cruzipain inhibitors and molecular modeling study of complexes with cruzain

The development of cruzipain inhibitors represents one of the most attractive challenges in the search for drugs for the treatment of Chagas disease. A recombinant form of this enzyme, cruzain, has been crystallized with numerous inhibitors, excluding thiosemicarbazones. These compounds have been established as potent inhibitors of cruzain, although there is very little data in the literature of thiosemicarbazones tested on cruzipain. In this work, we present the results of the evaluation of eleven thiosemicarbazones on cruzipain, isolated from T. cruzi epimastigotes, six of them previously evaluated on cruzain. For these latter, we studied through computational methods, the mode of interaction with the active site of cruzain and the contribution of geometric parameters to the possible mechanism of action involved in the observed inhibition. Finally, from some geometric parameters analyzed on modeled TSC-cruzain complexes, a semi-quantitative relationship was established that could explain the inhibitory activity of thiosemicarbazones on cruzipain, the enzyme actually present in the parasite.

Mass spectral characterization of organophosphate-labeled, tyrosine-containing peptides: characteristic mass fragments and a new binding motif for organophosphates

We have identified organophosphorus agent (OP)-tyrosine adducts on 12 different proteins labeled with six different OP. Labeling was achieved by treating pure proteins with up to 40-fold molar excess of OP at pH 8-8.6. OP-treated proteins were digested with trypsin, and peptides were separated by HPLC. Fragmentation patterns for 100 OP-peptides labeled on tyrosine were determined in the mass spectrometer. The goals of the present work were (1) to determine the common features of the OP-reactive tyrosines, and (2) to describe non-sequence MSMS fragments characteristic of OP-tyrosine peptides. Characteristic ions at 272 and 244 amu for tyrosine-OP immonium ions were nearly always present in the MSMS spectrum of peptides labeled on tyrosine by chlorpyrifos-oxon. Characteristic fragments also appeared from the parent ions that had been labeled with diisopropylfluorophosphate (216 amu), sarin (214 amu), soman (214 amu) or FP-biotin (227, 312, 329, 691 and 708 amu). In contrast to OP-reactive serines, which lie in the consensus sequence GXSXG, the OP-reactive tyrosines have no consensus sequence. Their common feature is the presence of nearby positively charged residues that activate the phenolic hydroxyl group. The significance of these findings is the recognition of a new binding motif for OP to proteins that have no active site serine. Modified peptides are difficult to find when the OP bears no radiolabel and no tag. The characteristic MSMS fragment ions are valuable because they are identifiers for OP-tyrosine, independent of the peptide.

Differing views of the role of selenium in thioredoxin reductase

This review covers three different chemical explanations that could account for the requirement of selenium in the form of selenocysteine in the active site of mammalian thioredoxin reductase. These views are the following: (1) the traditional view of selenocysteine as a superior nucleophile relative to cysteine, (2) the superior leaving group ability of a selenol relative to a thiol due to its significantly lower pK (a) and, (3) the superior ability of selenium to accept electrons (electrophilicity) relative to sulfur. We term these chemical explanations as the "chemico-enzymatic" function of selenium in an enzyme. We formally define the chemico-enzymatic function of selenium as its specific chemical property that allows a selenoenzyme to catalyze its individual reaction. However we, and others, question whether selenocysteine is chemically necessary to catalyze an enzymatic reaction since cysteine-homologs of selenocysteine-containing enzymes catalyze their specific enzymatic reactions with high catalytic efficiency. There must be a unique chemical reason for the presence of selenocysteine in enzymes that explains the biological pressure on the genome to maintain the complex selenocysteine-insertion machinery. We term this biological pressure the "chemico-biological" function of selenocysteine. We discuss evidence that this chemico-biological function is the ability of selenoenzymes to resist inactivation by irreversible oxidation. The way in which selenocysteine confers resistance to oxidation could be due to the superior ability of the oxidized form of selenocysteine (Sec-SeO(2)(-), seleninic acid) to be recycled back to its parent form (Sec-SeH, selenocysteine) in comparison to the same cycling of cysteine-sulfinic acid to cysteine (Cys-SO(2)(-) to Cys-SH).

Challenging a paradigm: theoretical calculations of the protonation state of the Cys25-His159 catalytic diad in free papain

A central mechanistic paradigm of cysteine proteases is that the His-Cys catalytic diad forms an ion-pair NH(+)/S(-) already in the catalytically active free enzyme. Most molecular modeling studies of cysteine proteases refer to this paradigm as their starting point. Nevertheless, several recent kinetics and X-ray crystallography studies of viral and bacterial cysteine proteases depart from the ion-pair mechanism, suggesting general base catalysis. We challenge the postulate of the ion-pair formation in free papain. Applying our QM/SCRF(VS) molecular modeling approach, we analyzed all protonation states of the catalytic diad in free papain and its SMe derivative, comparing the predicted and experimental pK(a) data. We conclude that the His-Cys catalytic diad in free papain is fully protonated, NH(+)/SH. The experimental pK(a) = 8.62 of His159 imidazole in free papain, obtained by NMR-controlled titration and originally interpreted as the NH(+)/S(-) <==> N/S(-) NH(+)/S(-) <==> N/S(-) equilibrium, is now assigned to the NH(+)/SH <==> N/SH NH(+)/SH <==> N/SH equilibrium.

Covalent binding of the organophosphorus agent FP-biotin to tyrosine in eight proteins that have no active site serine

Organophosphorus (OP) esters are known to bind covalently to the active site serine of enzymes in the serine hydrolase family. It was a surprise to find that proteins with no active site serine are also covalently modified by OP. The binding site in albumin, transferrin, and tubulin was identified as tyrosine. The goal of the present work was to determine whether binding to tyrosine is a general phenomenon. Fourteen proteins were treated with a biotin-tagged organophosphorus agent called FP-biotin. The proteins were digested with trypsin and the labeled peptides enriched by binding to monomeric avidin. Peptides were purified by HPLC and fragmented by collision induced dissociation in a tandem ion trap mass spectrometer. Eight proteins were labeled and six were not. Tyrosine was labeled in human alpha-2-glycoprotein 1 zinc-binding protein (Tyr 138, Tyr 174 and Tyr 181), human kinesin 3C motor domain (Tyr 145), human keratin 1 (Tyr 230), bovine actin (Tyr 55 and Tyr 200), murine ATP synthase beta (Tyr 431), murine adenine nucleotide translocase 1 (Tyr 81), bovine chymotrypsinogen (Tyr 201) and porcine pepsin (Tyr 310). Only 1-3 tyrosines per protein were modified, suggesting that the reactive tyrosine was activated by nearby residues that facilitated ionization of the hydroxyl group of tyrosine. These results suggest that OP binding to tyrosine is a general phenomenon. It is concluded that organophosphorus-reactive proteins include not only enzymes in the serine hydrolase family, but also proteins that have no active site serine. The recognition of a new OP-binding motif to tyrosine suggests new directions to search for mechanisms of long-term effects of OP exposure. Another application is in the search for biomarkers of organophosphorus agent exposure. Previous searches have been limited to serine hydrolases. Now proteins such as albumin and keratin can be considered.

Nonequilibrium capillary electrophoresis of equilibrium mixtures--a single experiment reveals equilibrium and kinetic parameters of protein-DNA interactions

We introduce a novel electrophoretic method, nonequilibrium capillary electrophoresis of equilibrium mixtures (NECEEM), and demonstrate its use for studying protein-DNA interactions. The equilibrium mixture of protein and DNA contains three components: free protein, free DNA, and the protein-DNA complex. A short plug of such a mixture is injected into the capillary, and the three components are separated under nonequilibrium conditions. The resulting electropherograms are composed of characteristic peaks and exponential curves. An easy nonnumerical analysis of a single electropherogram reveals two parameters: the equilibrium binding constant and the monomolecular rate constant of complex decay. The bimolecular rate constant of complex formation can then be calculated as the product of the two experimentally determined constants. NECEEM was applied to study the interaction between single-stranded DNA binding protein and a fluorescently labeled 15-mer oligonucleotide. It allowed us to measure for the first time the rate constant of complex decay for this important protein-DNA pair, k-1 = 0.03 s-1. The value of the equilibrium binding constant, Kb = 3.6 x 10-6 M-1, was in good agreement with those measured by other methods. As low as 10-18 mol of the protein was sufficient for the measurements. Thus, the new method is simple, informative, and highly sensitive. Moreover, it can be equally applied to other noncovalent protein-ligand complexes. These features of NECEEM make this method an indispensable tool in studies of macromolecular interactions. They also emphasize the potential role of NECEEM in the development of extremely sensitive protein assays using nucleotide aptamers.

Computational studies on nonenzymatic and enzymatic pyridoxal phosphate catalyzed decarboxylations of 2-aminoisobutyrate

A computational study of nonenzymatic and enzymatic pyridoxal phosphate-catalyzed decarboxylation of 2-aminoisobutyrate (AIB) is presented. Four prototropic isomers of a model aldimine between AIB and 5'-deoxypyridoxal, with acetate interacting with the pyridine nitrogen, were employed in calculations of both gas phase and water model (PM3 and PM3-SM3) decarboxylation reaction paths. Calculations employing the transition state structures obtained for the four isomers allow the demonstration of stereoelectronic effects in transition state stabilization as well as a separation of the contributions of the Schiff base and pyridine ring moieties to this stabilization. The unprotonated Schiff base contribution (approximately 16 kcal/mol) is larger than that of the pyridine ring even when it is protonated (approximately 10 kcal/mol), providing an explanation of the catalytic power of pyruvoyl-dependent amino acid decarboxylases. An active site model of dialkylglycine decarboxylase was constructed and validated, and enzymatic decarboxylation reaction paths were calculated. The reaction coordinate is shown to be complex, with proton transfer from Lys272 to the coenzyme C4' likely simultaneous with C alpha--CO(2)(-) bond cleavage. The proposed concerted decarboxylation/proton-transfer mechanism provides a simple explanation for the observed specificity of this enzyme toward oxidative decarboxylation.

Crystal structure of a catalytic antibody with a serine protease active site

The three-dimensional structure of an unusually active hydrolytic antibody with a phosphonate transition state analog (hapten) bound to the active site has been solved to 2.5 A resolution. The antibody (17E8) catalyzes the hydrolysis of norleucine and methionine phenyl esters and is selective for amino acid esters that have the natural alpha-carbon L configuration. A plot of the pH-dependence of the antibody-catalyzed reaction is bell-shaped with an activity maximum at pH 9.5; experiments on mechanism lend support to the formation of a covalent acyl-antibody intermediate. The structural and kinetic data are complementary and support a hydrolytic mechanism for the antibody that is remarkably similar to that of the serine proteases. The antibody active site contains a Ser-His dyad structure proximal to the phosphorous atom of the bound hapten that resembles two of the three components of the Ser-His-Asp catalytic triad of serine proteases. The antibody active site also contains a Lys residue to stabilize oxyanion formation, and a hydrophobic binding pocket for specific substrate recognition of norleucine and methionine side chains. The structure identifies active site residues that mediate catalysis and suggests specific mutations that may improve the catalytic efficiency of the antibody. This high resolution structure of a catalytic antibody-hapten complex shows that antibodies can converge on active site structures that have arisen through natural enzyme evolution.

Reversible covalent binding of peptide nitriles to papain

The dissociation constants for reversible covalent binding of twelve peptide nitrile inhibitors to the active site of papain have been measured by means of fluorescence titration. The binding constants generally parallel the kinetic specificity constants (kcat/Km) for related papain substrates, supporting earlier suggestions that peptide nitriles behave as transition state analog inhibitors of papain. In ten cases the temperature dependence of binding was analyzed to determine the enthalpic and entropic contributions to the binding energy. A compensation plot of delta H vs. T delta S resulted in two parallel lines, one for 'specific' nitriles (i.e., N-Ac-L-aa-NHCH2CN; aa = Phe, Leu, Met) and the other for 'non-specific' nitriles (e.g., N-Ac-D-Phe-NHCH2CN, PhCH2CH2CONHCH2CN hippurylnitrile, etc.). For both specific and nonspecific nitriles representing an 1800-fold range of Kd values (0.27 microM-490 microM), the solvent deuterium isotope effect on binding (Kd(H2O)/Kd(D2O) = DKd) was very close to 2.0. This isotope effect could be accounted for entirely by the simple protonic change which occurs upon the reversible addition of the active site sulfhydryl of papain to the nitrile group of the peptide derivative to form a covalent thioimidate linkage. In contrast, six closely related non-nitrile ligands containing identical peptide side chains but having C-terminal groups incapable of binding covalently to papain had unmeasureably high dissociation constants. Collectively, these results indicate that strong binding of peptide nitrile substrate analogs to papain requires a combination of (1) hydrophobic interaction (especially at the P2 position), (2) specific intermolecular hydrogen bonding and (3) covalent interaction of the nitrile with the active site sulfhydryl group.

Active-site ionizations of papain. An evaluation of the potentiometric difference titration method

The underlying assumption of the potentiometric difference titration method, as applied to the evaluation of the sulfhydryl-dependent ionizations in the active site of papain, is that the pKa of His-159 is independent of whether the neighboring sulfhydryl (Cys-25) is protonated or methylthiolated. That this idealized assumption may not strictly apply is indirectly indicated by the larger pKa of His-159 in S-methylpapain versus that in S-methylthiopapain, as determined from fluorometric titrations (delta pKa = 0.32 +/- 0.05, 25 degrees C). On the basis of the Wegscheider principle of the equivalence of protons and methyl groups, the potentiometric difference titration method will underestimate the concentration of thiolate-imidazolium ion pair in the active site versus that of the thiol-imidazole tautomer, provided that there is no significant H-bonding interaction in the latter species.

Electrostatic control of enzyme reactions: the mechanism of inhibition of glucose oxidase by putrescine

The interaction of putrescine dihydrochloride with glucose oxidase is reported. At pH 7.65 glucose oxidase is strongly anionic (Z = -80). The pKa of an essential acidic group on the reduced form of the enzyme is extremely sensitive to ionic strength, as predicted by simple electrostatic theory [J. G. Voet, J. Coe, J. Epstein, V. Matossian, and T. Shipley (1981) Biochemistry 20, 7182-7185]. Putrescine dihydrochloride was found to inhibit glucose oxidase at pH 7.65 at a constant ionic strength of 0.05. The kinetics do not obey simple competitive inhibition, however. The data can best be explained by a model in which change in the electrostatic potential of the enzyme on putrescine binding changes the observed pKa of the essential acidic group. The pH dependence of putrescine inhibition supports this interpretation. At I = 0.05, 5 mM putrescine was found to change the pKa of the essential acidic group from 7.6 to 7.1. The shift in the pKa as a function of putrescine concentration at pH 7.7 and I = 0.05 also supports the model presented. The Ka for putrescine to the active form of the enzyme was calculated to be 4.2 mM.

Calculations of electrostatic interactions in biological systems and in solutions

Correlating the structure and action of biological molecules requires knowledge of the corresponding relation between structure and energy. Probably the most important factors in such a structure– energy correlation are associated with electrostatic interactions. Thus the key requirement for quantative understanding of the action of biological molecules is the ability to correlate electrostatic interactions with structural information. To appreciate this point it is useful to compare the electrostatic energy of a charged amino acid in a polar solvent to the corresponding van der Waals energy. The electrostatic free energy, ΔGel, can be approximated (as will be shown in Section II) by the Born formula (ΔGel= –(166Q2/ā) (I – I/E)). Where ΔGelis given in kcal/mol,Qis the charge of the given group, in units of electron charge,āis the effective radius of the group, andEis the dielectric constant of the solvent. With an effective radius of charged amino acids of approximately 2 Å, Born's formula gives about – 80 kcal/mol for their energy in polar solvents whereEis larger than 10. This energy is two orders of magnitude larger than the van der Waals interaction of such groups and their surroundings.

Electrostatic control of enzyme reactions: effect of ionic strength on the pKa of an essential acidic group on glucose oxidase

The dissociation constant of an essential acidic group on the reduced form of glucose oxidase from Aspergillus niger (K4) has been found to be extremely sensitive to ionic strength. Increasing the ionic strength from 0.025 to 0.225 causes a decrease in pK4,obsd of 0.9 pH unit, from 8.2 to 7.3. Analysis of the ionic strength dependence of pK4,obsd, making the assumption that the enzyme is a homogeneously charged impenetrable sphere [Edsall, J. T., & Wyman, J. (1958) Biophysical Chemistry, Vol. 1, pp 282-289, 512-514, Academic Press, New York], predicts that the intrinsic pKa of the acidic group is 6.7 and that the charge on the protein is -78. The enzyme was titrated from its isoelectric point (pH 4.05) to pH 7.7, the pH at which the ionic strength dependence was determined. It was found to have an actual charge at that pH of -77, in remarkable agreement with the theoretical prediction. Thus, glucose oxidase exerts electrostatic control on pK4,obsd as though it were a uniformly charged sphere. The group responsible for pK4,obsd has not been identified. However, its measured delta H degrees obsd of 8.0 kcal mol-1 and delta S degrees obsd of -6.1 cal mol-1 K-1, together with its pKa of 6.7, are consistent with the group being a histidine residue.

Kinetic solvent isotope effects on the deacylation of specific acyl-papains. Proton inventory studies on the papain-catalysed hydrolyses of specific ester substrates: analysis of possible transition state structures

1. The hydrolyses of the p-nitrophenyl esters of N-benzyloxycarbonylglycine, alpha-N-benzyloxycarbonyl-L-lysine and N-methoxycarbonyl-L-phenylalanylglycine catalysed by papain (EC 3.4.22.2) have been studied in solvents having a variable composition of 2H2O and H2O. 2. kcat., which represents deacylation in the papain-catalysed hydrolysis of reactive esters, is some 2.3-fold less in 2H2O compared with H2O. The magnitude of kcat. has been determined as a function of the 2H atom fraction of the solvent. 3. Both linear and non-linear methods of least-square regression analysis have been applied to the data in order to obtain best-fit parameter values for several three-parameter models which express kcat. in terms of the 2H atom fraction of the solvent. These models represent some possible modes of restructuring of the active site protonic configuration consequent upon transition state formation. 4. The results of curve fitting reveal an essentially linear dependence of kcat. upon the 2H atom fraction, and it may therefore be concluded that the isotope effect originates from a single proton which is in the process of transfer in the transition state. 5. It is postulated on the basis of this and other evidence that the mobile proton is transferred from an attacking water molecule to the imidazole side chain of His-159 during tetrahedral intermediate formation. This has the effect of stabilizing the transition state and promoting catalysis. The role of His-159 in deacylation is therefore to provide general base catalysis. 6. Models that involve two or more protons, such as a two-proton relay system analogous to that proposed for the serine proteinases, or a multiproton 'medium' effect, are considered unlikely on the basis of the data reported in this paper. 7. A more detailed examination of possible transition state structures reveals that the only structure compatible with available experimental data and consistent with certain theoretical predictions is one in which the proton translocated in concern with reorganization of the heavy atom framework. In addition, the transition state vibrations of the mobile proton are strongly coupled to those of the heavy atoms. These properties of the transition state are also manifest in the transition state for the deacylation of serine proteinases.