Engineering substrate preference in subtilisin: structural and kinetic analysis of a specificity mutant

Bacillus subtilisin has been a popular model protein for engineering altered substrate specificity. Although some studies have succeeded in increasing the specificity of subtilisin, they also demonstrate that high specificity is difficult to achieve solely by engineering selective substrate binding. In this paper, we analyze the structure and transient state kinetic behavior of Sbt160, a subtilisin engineered to strongly prefer substrates with phenylalanine or tyrosine at the P4 position. As in previous studies, we measure improvements in substrate affinity and overall specificity. Structural analysis of an inactive version of Sbt160 in complex with its cognate substrate reveals improved interactions at the S4 subsite with a P4 tyrosine. Comparison of transient state kinetic behavior against an optimal sequence (DFKAM) and a similar, but suboptimal, sequence (DVRAF) reveals the kinetic and thermodynamic basis for increased specificity, as well as the limitations of this approach. While highly selective substrate binding is achieved in Sbt160, several factors cause sequence specificity to fall short of that observed with natural processing subtilisins. First, for substrate sequences which are nearly optimal, the acylation reaction becomes faster than substrate dissociation. As a result, the level of discrimination among these substrates diminishes due to the coupling between substrate binding and the first chemical step (acylation). Second, although Sbt160 has 24-fold higher substrate affinity for the optimal substrate DFKAM than for DVRAF, the increased substrate binding energy is not translated into improved transition state stabilization of the acylation reaction. Finally, as interactions at subsites become stronger, the rate-determining step in peptide hydrolysis changes from acylation to product release. Thus, the release of the product becomes sluggish and leads to a low k(cat) for the reaction. This also leads to strong product inhibition of substrate turnover as the reaction progresses. The structural and kinetic analysis reveals that differences in the binding modes at subsites for substrates, transition states, and products are subtle and difficult to manipulate via straightforward protein engineering. These findings suggest several new strategies for engineering highly sequence selective enzymes.     

Conformational cycling in beta-phosphoglucomutase catalysis: reorientation of the beta-D-glucose 1,6-(Bis)phosphate intermediate

Activated Lactococcus lactis beta-phosphoglucomutase (betaPGM) catalyzes the conversion of beta-d-glucose 1-phosphate (betaG1P) derived from maltose to beta-d-glucose 6-phosphate (G6P). Activation requires Mg(2+) binding and phosphorylation of the active site residue Asp8. Initial velocity techniques were used to define the steady-state kinetic constants k(cat) = 177 +/- 9 s(-)(1), K(m) = 49 +/- 4 microM for the substrate betaG1P and K(m) = 6.5 +/- 0.7 microM for the activator beta-d-glucose 1,6-bisphosphate (betaG1,6bisP). The observed transient accumulation of [(14)C]betaG1,6bisP (12% at approximately 0.1 s) in the single turnover reaction carried out with excess betaPGM (40 microM) and limiting [(14)C]betaG1P (5 microM) and betaG1,6bisP (5 microM) supported the role of betaG1,6bisP as a reaction intermediate in the conversion of the betaG1P to G6P. Single turnover reactions of [(14)C]betaG1,6bisP with excess betaPGM were carried out to demonstrate that phosphoryl transfer rather than ligand binding is rate-limiting and to show that the betaG1,6bisP binds to the active site in two different orientations (one positioning the C(1)phosphoryl group for reaction with Asp8, and the other orientation positioning the C(6)phosphoryl group for reaction with Asp8) with roughly the same efficiency. Single turnover reactions carried out with betaPGM, [(14)C]betaG1P, and unlabeled betaG1,6bisP demonstrated complete exchange of label to the betaG1,6bisP during the catalytic cycle. Thus, the reorientation of the betaG1,6bisP intermediate that is required to complete the catalytic cycle occurs by diffusion into solvent followed by binding in the opposite orientation. Published X-ray structures of betaG1P suggest that the reorientation and phosphoryl transfer from betaG1,6bisP occur by conformational cycling of the enzyme between the active site open and closed forms via cap domain movement. Last, the equilibrium ratio of betaG1,6bisP to betaG1P plus G6P was examined to evidence a significant stabilization of betaPGM aspartyl phosphate.     

Interaction of nitric oxide with catalase: structural and kinetic analysis

We present the structures of bovine catalase in its native form and complexed with ammonia and nitric oxide, obtained by X-ray crystallography. Using the NO generator 1-(N,N-diethylamino)diazen-1-ium-1,2-diolate, we were able to generate sufficiently high NO concentrations within the catalase crystals that substantial occupation was observed despite a high dissociation rate. Nitric oxide seems to be slightly bent from the heme normal that may indicate some iron(II) character in the formally ferric catalase. Microspectrophotometric investigations inline with the synchrotron X-ray beam reveal photoreduction of the central heme iron. In the cases of the native and ammonia-complexed catalase, reduction is accompanied by a relaxation phase. This is likely not the case for the catalase NO complex. The kinetics of binding of NO to catalase were investigated using NO photolyzed from N,N'-bis(carboxymethyl)-N,N'-dinitroso-p-phenylenediamine using an assay that combines catalase with myoglobin binding kinetics. The off rate is 1.5 s(-1). Implications for catalase function are discussed.     

Thermodynamic and kinetic determinants of Thermotoga maritima cold shock protein stability: a structural and dynamic analysis

The cold shock protein (CSP) from hyperthermophile Thermotoga maritima (TmCSP) is only marginally stable (DeltaG(T(opt)) = 0.3 kcal/mol) at 353 K, the optimum environmental temperature (T(opt)) for T. maritima. In comparison, homologous CSPs from E. coli (DeltaG(T(opt)) = 2.2 kcal/mol) and B. subtilis (DeltaG(T(opt)) = 1.5 kcal/mol) are at least five times more stable at 310 K, the T(opt) for the mesophiles. Yet at the room temperature, TmCSP is more stable (DeltaG(T(R)) = 4.7 kcal/mol) than its homologues (DeltaG(T(R)) = 3.0 kcal/mol for E. coli CSP and DeltaG(T(R)) = 2.1 kcal/mol for B. subtilis CSP). This unique observation suggests that kinetic, rather than thermodynamic, barriers toward unfolding might help TmCSP native structure at high temperatures. Consistently, the unfolding rate of TmCSP is considerably slower than its homologues. High temperature (600 K) complete unfolding molecular dynamics (MD) simulations of TmCSP support our hypothesis and reveal an unfolding scheme unique to TmCSP. For all the studied homologues of TmCSP, the unfolding process first starts at the C-terminal region and N-terminal region unfolds in the end. But for TmCSP, both the terminals resist unfolding for consistently longer simulation times and, in the end, unfold simultaneously. In TmCSP, the C-terminal region is better fortified and has better interactions with the N-terminal region due to the charged residues, R2, E47, E49, H61, K63, and E66, being in spatial vicinity. The electrostatic interactions among these residues are unique to TmCSP. Consistently, the room temperature MD simulations show that TmCSP is more rigid at its N- and C-termini as compared to its homologues from E. coli, B. subtilis, and B. caldolyticus.     

Catalytic activity of aspartate aminotransferase in the crystal. Equilibrium and kinetic analysis.

Natural substrates and analogs rapidly diffuse through crystals of pig heart mitochondrial aspartate aminotransferase and react at the active sites causing spectral changes that can be measured by single-crystal microspectrophotometry. Dissociation constants for natural substrates and rate constants of transamination for slowly reacting substrates have been determined. A comparison between the data obtained in the crystal and in solution shows that the crystalline enzyme is catalytically competent and that events occurring in the crystal essentially parallel those occurring in solution, even though minor differences have been detected.     

Catalytic improvement and structural analysis of atrazine chlorohydrolase by site-saturation mutagenesis

To improve the catalytic activity of atrazine chlorohydrolase (AtzA), amino acid residues involved in substrate binding (Gln71) and catalytic efficiency (Val12, Ile393, and Leu395) were targeted to generate site-saturation mutagenesis libraries. Seventeen variants were obtained through Haematococcus pluvialis-based screening, and their specific activities were 1.2–5.2-fold higher than that of the wild type. For these variants, Gln71 tended to be substituted by hydrophobic amino acids, Ile393 and Leu395 by polar ones, especially arginine, and Val12 by alanine, respectively. Q71R and Q71M significantly decreased the K_m by enlarging the substrate-entry channel and affecting N-ethyl binding. Mutations at sites 393 and 395 significantly increased the k_cat/K_m, probably by improving the stability of the dual β-sheet domain and the whole enzyme, owing to hydrogen bond formation. In addition, the contradictory relationship between the substrate affinity improvement by Gln71 mutation and the catalytic efficiency improvement by the dual β-sheet domain modification was discussed.     

Catalytic cycling of human mitochondrial Lon protease

1. Download : Download high-res image (346KB)
2. Download : Download full-size image

     

Cryoenzymology: How to practice kinetic and structural studies

For a full understanding of an enzyme reaction pathway, one must identify the reaction intermediates and obtain their structures and rates of interconversion. It is impossible to obtain all this information under normal conditions. An approach is to work suboptimally, in particular at subzero temperatures. This is cryoenzymology, an approach that implies both kinetic and structural measurements on enzyme systems below 0°C. To work below 0°C one must add an antifreeze so cryoenzymology means perturbation by two agents: temperature and antifreeze, usually an organic solvent. Certain precautions are needed with these agents, which we will discuss here. In particular, we discuss the importance of choosing the right solvent: this requires extensive exploratory studies but it is the key for the successful practice of cryoenzymology. Each system has its particularities and its own ‘good’ solvent. Cryoenzymology is not only a way of reducing reaction rates, it is also a way of perturbing one's system. Thus, it is a method that allows for the accumulation of intermediates that cannot be observed under normal conditions by slowing down their kinetics of formation, by changes in rate limiting steps or by shifts in equilibria. We illustrate the usefulness of cryoenzymology by myosin and actomyosin ATPases and by creatine, arginine and 3-phosphoglycerate kinases. We also discuss recent results obtained by X-ray crystallography.     

Catalytic mechanism of biotin carboxylase: steady-state kinetic investigations

Biotin carboxylase was purified from Escherichia coli by a new procedure, and its steady-state kinetic parameters were examined. MgATP and bicarbonate add to the enzyme randomly, followed by addition of biotin. Both bicarbonate and MgATP add in rapid equilibrium. A catalytic base with a pK of 6.6 is observed in V/K profiles. Inactivation studies also revealed a sulfhydryl group in the active site that is essential for catalysis. It is proposed that the acid-base catalysts are necessary for the tautomerization of biotin, which presumably enhances its nucleophilicity toward the carboxyl group donor. A second enzymic group with a pK of 6.6, whose role is unknown, is seen in Vmax profiles. The pH profiles for the biotin carboxylase catalyzed phosphorylation of ADP by carbamoyl phosphate have the same shape as the profiles for the forward reaction, which demonstrates that the enzymic bases assume the same protonation states for catalysis of transphosphorylation in either direction. The lack of reactivity of thionucleotide analogues of ATP when Mg is used as the divalent metal ion suggests that both metal ions required for reaction coordinate to the nucleotide. The second metal ion appears to be absolutely required for reaction and not merely an activator of the reaction. Characterization of a bicabonate-dependent biotin-independent ATPase activity strongly suggests that carboxylation proceeds via a carboxyphosphate intermediate.     

Catalytic oxidation of acetaminophen by tyrosinase in the presence of L-proline: a kinetic study

A kinetic study of acetaminophen oxidation by tyrosinase in the presence of a physiological nucleophilic agent such as the amino acid L-proline is performed in the present paper. The o-quinone product of the catalytic activity, 4-acetamido-o-benzoquinone, becomes unstable through the chemical addition of L-proline, in competition with the nucleophilic addition of hydroxide ion from water. In both cases, the catechol intermediate, 3(')-hydroxyacetaminophen, is generated, as can be demonstrated by liquid chromatography. When the effect of the presence of the nucleophilic agent on the time course of the enzymatic reaction was kinetically analyzed, it was seen to decrease the duration of the lag period and increase the steady-state rate. Rate constants for the reaction of 4-acetamido-o-benzoquinone with water and L-proline were also determined. The results obtained in this paper open a new possibility to acetaminophen toxicity, that has been attributed hitherto to its corresponding p-quinone, N-acetyl-p-benzoquinone imine.     

Protein-RNA interactions: a structural analysis

A detailed computational analysis of 32 protein–RNA complexes is presented. A number of physical and chemical properties of the intermolecular interfaces are calculated and compared with those observed in protein–double-stranded DNA and protein–single-stranded DNA complexes. The interface properties of the protein–RNA complexes reveal the diverse nature of the binding sites. van der Waals contacts played a more prevalent role than hydrogen bond contacts, and preferential binding to guanine and uracil was observed. The positively charged residue, arginine, and the single aromatic residues, phenylalanine and tyrosine, all played key roles in the RNA binding sites. A comparison between protein–RNA and protein–DNA complexes showed that whilst base and backbone contacts (both hydrogen bonding and van der Waals) were observed with equal frequency in the protein–RNA complexes, backbone contacts were more dominant in the protein–DNA complexes. Although similar modes of secondary structure interactions have been observed in RNA and DNA binding proteins, the current analysis emphasises the differences that exist between the two types of nucleic acid binding protein at the atomic contact level.     

Membrane cycling and macrovacuolation under the influence of cytochalasin: kinetic and morphometric studies

Fibroblasts exposed to higher doses of cytochalasin accumulate very big discrete endoplasmic vacuoles, the membrane of which is derived by internalization of plasmalemma. Morphometry confirms that the amount of surface interiorized is equal to the difference between the original cell surface area (before CD) and the reduced surface area measurable after CD-induced rounding. Correspondingly, there is a nearly two-fold increase in the activity of the ectoenzyme 5'-nucleotidase (a marker for plasma membrane) internally within the cytoplasm, after treatment with CD. Macrovacuolation increases cell volume by approximately 30%. Surface membrane is internalized as micropinocytotic vesicles at a rate measurable by the accumulation of HRP, a marker of fluid-phase pinocytosis. Uptake of HRP is shown to be enhanced at all times during exposure to CD, and is balanced by accelerated exocytic recycling of membrane except during a phase (approximately 4-8 hr) in which pinocytic uptake exceeds exocytosis. Vesicular membrane accumulated intracellularly in this period is retained in the endoplasm, and by successive fusions forms vacuoles in close approximation to microfilament aggregates. Once established, this new macrovacuolar membrane compartment is in dynamic equilibrium with the cell surface, and its membrane is cycled like the plasma membrane, in a mutual exchange of pinosomes between the several vacuoles and the cell surface. In drug-free medium vacuole membrane apparently reverts to the surface by pinocytotic recycling, and the cells recover normal characteristics 4-6 hr after withdrawal of cytochalasin.     

Glutathione S-transferase from beef heart: structural and kinetic properties

A simple and rapid method for the purification of glutathione S-transferase is described. The physical and kinetic properties of purified enzyme are reported. The protein is constituted of two identical subunits with a total molecular weight of 46,000 daltons. The isoelectric focusing of crude cytosol or purified preparation gives a single peak of activity with a pI of 7.1. The kinetic analysis shows a relatively strict substrate specificity. Only 1-chloro-2,4-dinitrobenzene is conjugated to reduced glutathione at an appreciable rate. The peroxidase activity of the enzyme with respect to cumene hydroperoxide as substrate is negligible. Hemin and bilirubin are competitive inhibitors of catalytic activity.     

Catalytic pyrolysis of tobacco rob: kinetic study and fuel gas produced

The pyrolysis kinetics of tobacco rob (TR) was investigated using thermogravimetric analysis (TGA) under inert atmosphere, adding chemicals (dolomite and NiO) as catalysts by catalytic-mixing method. The TGA results showed that mass loss and mass loss rates were affected by catalysts. The conversion rates increased while the activation energy decreased. Moreover, the thermal decomposition behaviors of TR were studied in the fixed-bed reactor using dolomite and NiO/γ-Al₂O₃ as catalysts by catalyst-bed method. A series of experiments had been performed to explore the effects of catalysts, and reaction temperature on the composition and yield of fuel gas. The experiments demonstrated that the catalysts had a high activity of cracking tar and hydrocarbons, as well as yielding a high fuel gas production. For both methods, dolomite and NiO revealed better catalytic performance as a view of enhancing conversion rates and increasing product gas yield.     

Late events of translation initiation in bacteria: a kinetic analysis

Binding of the 50S ribosomal subunit to the 30S initiation complex and the subsequent transition from the initiation to the elongation phase up to the synthesis of the first peptide bond represent crucial steps in the translation pathway. The reactions that characterize these transitions were analyzed by quench-flow and fluorescence stopped-flow kinetic techniques. IF2-dependent GTP hydrolysis was fast (30/s) followed by slow P(i) release from the complex (1.5/s). The latter step was rate limiting for subsequent A-site binding of EF-Tu small middle dotGTP small middle dotPhe-tRNA(Phe) ternary complex. Most of the elemental rate constants of A-site binding were similar to those measured on poly(U), with the notable exception of the formation of the first peptide bond which occurred at a rate of 0.2/s. Omission of GTP or its replacement with GDP had no effect, indicating that neither the adjustment of fMet-tRNA(fMet) in the P site nor the release of IF2 from the ribosome required GTP hydrolysis.     

A substrate-induced switch in the reaction mechanism of a thermophilic esterase: kinetic evidences and structural basis

The reaction mechanism of the esterase 2 (EST2) from Alicyclobacillus acidocaldarius was studied at the kinetic and structural level to shed light on the mechanism of activity and substrate specificity increase previously observed in its double mutant M211S/R215L. In particular, the values of kinetic constants (k1, k(-1), k2, and k3) along with activation energies (E1, E(-1), E2, and E3) were measured for wild type and mutant enzyme. The previously suggested substrate-induced switch in the reaction mechanism from kcat=k3 with a short acyl chain substrate (p-nitrophenyl hexanoate) to kcat=k2 with a long acyl chain substrate (p-nitrophenyl dodecanoate) was validated. The inhibition afforded by an irreversible inhibitor (1-hexadecanesulfonyl chloride), structurally related to p-nitrophenyl dodecanoate, was studied by kinetic analysis. Moreover the three-dimensional structure of the double mutant bound to this inhibitor was determined, providing essential information on the enzyme mechanism. In fact, structural analysis explained the observed substrate-induced switch because of an inversion in the binding mode of the long acyl chain derivatives with respect to the acyl- and alcohol-binding sites.     

Non-cognate enzyme-DNA complex: structural and kinetic analysis of EcoRV endonuclease bound to the EcoRI recognition site GAATTC

The crystal structure of EcoRV endonuclease bound to non-cognate DNA at 2.0 angstroms resolution shows that very small structural adaptations are sufficient to ensure the extreme sequence specificity characteristic of restriction enzymes. EcoRV bends its specific GATATC site sharply by 50 degrees into the major groove at the center TA step, generating unusual base-base interactions along each individual DNA strand. In the symmetric non-cognate complex bound to GAATTC, the center step bend is relaxed to avoid steric hindrance caused by the different placement of the exocyclic thymine methyl groups. The decreased base-pair unstacking in turn leads to small conformational rearrangements in the sugar-phosphate backbone, sufficient to destabilize binding of crucial divalent metal ions in the active site. A second crystal structure of EcoRV bound to the base-analog GAAUTC site shows that the 50 degrees center-step bend of the DNA is restored. However, while divalent metals bind at high occupancy in this structure, one metal ion shifts away from binding at the scissile DNA phosphate to a position near the 3'-adjacent phosphate group. This may explain why the 10(4)-fold attenuated cleavage efficiency toward GAATTC is reconstituted by less than tenfold toward GAAUTC. Examination of DNA binding and bending by equilibrium and stopped-flow florescence quenching and fluorescence resonance energy transfer (FRET) methods demonstrates that the capacity of EcoRV to bend the GAATTC non-cognate site is severely limited, but that full bending of GAAUTC is achieved at only a threefold reduced rate compared with the cognate complex. Together, the structural and biochemical data demonstrate the existence of distinct mechanisms for ensuring specificity at the bending and catalytic steps, respectively. The limited conformational rearrangements observed in the EcoRV non-cognate complex provide a sharp contrast to the extensive structural changes found in a non-cognate BamHI-DNA crystal structure, thus demonstrating a diversity of mechanisms by which restriction enzymes are able to achieve specificity.     

Caught in the act: the structure of phosphorylated beta-phosphoglucomutase from Lactococcus lactis

Phosphoglucomutases catalyze the interconversion of D-glucose 1-phosphate and D-glucose 6-phosphate, a reaction central to energy metabolism in all cells and to the synthesis of cell wall polysaccharides in bacterial cells. Two classes of phosphoglucomutases (alpha-PGM and beta-PGM) are distinguished on the basis of their specificity for alpha- and beta-glucose-1-phosphate. beta-PGM is a member of the haloacid dehalogenase (HAD) superfamily, which includes the sarcoplasmic Ca(2+)-ATPase, phosphomannomutase, and phosphoserine phosphatase. beta-PGM is unusual among family members in that the common phosphoenzyme intermediate exists as a stable ground-state complex in this enzyme. Herein we report, for the first time, the three-dimensional structure of a beta-PGM and the first view of the true phosphoenzyme intermediate in the HAD superfamily. The crystal structure of the Mg(II) complex of phosphorylated beta-phosphoglucomutase (beta-PGM) from Lactococcus lactis has been determined to 2.3 A resolution by multiwavelength anomalous diffraction (MAD) phasing on selenomethionine, and refined to an R(cryst) = 0.24 and R(free) = 0.28. The active site of beta-PGM is located between the core and the cap domain and is freely solvent accessible. The residues within a 6 A radius of the phosphorylated Asp8 include Asp10, Thr16, Ser114, Lys145, Glu169, and Asp170. The cofactor Mg(2+) is liganded with octahedral coordination geometry by the carboxylate side chains of Asp8, Glu169, Asp170, and the backbone carbonyl oxygen of Asp10 along with one oxygen from the Asp8-phosphoryl group and one water ligand. The phosphate group of the phosphoaspartyl residue, Asp8, interacts with the side chains of Ser114 and Lys145. The absence of a base residue near the aspartyl phosphate group accounts for the persistence of the phosphorylated enzyme under physiological conditions. Substrate docking shows that glucose-6-P can bind to the active site of phosphorylated beta-PGM in such a way as to position the C(1)OH near the phosphoryl group of the phosphorylated Asp8 and the C(6) phosphoryl group near the carboxylate group of Asp10. This result suggests a novel two-base mechanism for phosphoryl group transfer in a phosphorylated sugar.     

Catalytic mechanism of C-C hydrolase MhpC from Escherichia coli: kinetic analysis of His263 and Ser110 site-directed mutants

C-C hydrolase MhpC (2-hydroxy-6-keto-nona-1,9-dioic acid 5,6-hydrolase) from Escherichia coli catalyses the hydrolytic C-C cleavage of the meta-ring fission product on the phenylpropionic acid catabolic pathway. The crystal structure of E. coli MhpC has revealed a number of active-site amino acid residues that may participate in catalysis. Site-directed mutants of His263, Ser110, His114, and Ser40 have been analysed using steady-state and stopped-flow kinetics. Mutants H263A, S110A and S110G show 10(4)-fold reduced catalytic efficiency, but still retain catalytic activity for C-C cleavage. Two distinct steps are observed by stopped-flow UV/Vis spectrophotometry, corresponding to ketonisation and C-C cleavage: H263A exhibits very slow ketonisation and C-C cleavage, whereas S110A and S110G exhibit fast ketonisation, an intermediate phase, and slow C-C cleavage. H114A shows only twofold-reduced catalytic efficiency, ruling out a catalytic role, but shows a fivefold-reduced K(M) for the natural substrate, and an ability to process an aryl-containing substrate, implying a role for His114 in positioning of the substrate. S40A shows only twofold-reduced catalytic efficiency, but shows a very fast (500 s(-1)) interconversion of dienol (317 nm) to dienolate (394 nm) forms of the substrate, indicating that the enzyme accepts the dienol form of the substrate. These data imply that His263 is responsible for both ketonisation of the substrate and for deprotonation of water for C-C cleavage, a novel catalytic role in a serine hydrolase. Ser110 has an important but non-essential role in catalysis, which appears not to be to act as a nucleophile. A catalytic mechanism is proposed involving stabilisation of reactive intermediates and activation of a nucleophilic water molecule by Ser110.