Hydrophobic Substrate Binding Research Articles

Structure and function of an insect α-carboxylesterase (αEsterase7) associated with insecticide resistance

Insect carboxylesterases from the αEsterase gene cluster, such as αE7 (also known as E3) from the Australian sheep blowfly Lucilia cuprina (LcαE7), play an important physiological role in lipid metabolism and are implicated in the detoxification of organophosphate (OP) insecticides. Despite the importance of OPs to agriculture and the spread of insect-borne diseases, the molecular basis for the ability of α-carboxylesterases to confer OP resistance to insects is poorly understood. In this work, we used laboratory evolution to increase the thermal stability of LcαE7, allowing its overexpression in Escherichia coli and structure determination. The crystal structure reveals a canonical α/β-hydrolase fold that is very similar to the primary target of OPs (acetylcholinesterase) and a unique N-terminal α-helix that serves as a membrane anchor. Soaking of LcαE7 crystals in OPs led to the capture of a crystallographic snapshot of LcαE7 in its phosphorylated state, which allowed comparison with acetylcholinesterase and rationalization of its ability to protect insects against the effects of OPs. Finally, inspection of the active site of LcαE7 reveals an asymmetric and hydrophobic substrate binding cavity that is well-suited to fatty acid methyl esters, which are hydrolyzed by the enzyme with specificity constants (∼10(6) M(-1) s(-1)) indicative of a natural substrate.

Large Scale Structural Rearrangement of a Serine Hydrolase from Francisella tularensis Facilitates Catalysis

Tularemia is a deadly, febrile disease caused by infection by the gram-negative bacterium, Francisella tularensis. Members of the ubiquitous serine hydrolase protein family are among current targets to treat diverse bacterial infections. Herein we present a structural and functional study of a novel bacterial carboxylesterase (FTT258) from F. tularensis, a homologue of human acyl protein thioesterase (hAPT1). The structure of FTT258 has been determined in multiple forms, and unexpectedly large conformational changes of a peripheral flexible loop occur in the presence of a mechanistic cyclobutanone ligand. The concomitant changes in this hydrophobic loop and the newly exposed hydrophobic substrate binding pocket suggest that the observed structural changes are essential to the biological function and catalytic activity of FTT258. Using diverse substrate libraries, site-directed mutagenesis, and liposome binding assays, we determined the importance of these structural changes to the catalytic activity and membrane binding activity of FTT258. Residues within the newly exposed hydrophobic binding pocket and within the peripheral flexible loop proved essential to the hydrolytic activity of FTT258, indicating that structural rearrangement is required for catalytic activity. Both FTT258 and hAPT1 also showed significant association with liposomes designed to mimic bacterial or human membranes, respectively, even though similar structural rearrangements for hAPT1 have not been reported. The necessity for acyl protein thioesterases to have maximal catalytic activity near the membrane surface suggests that these conformational changes in the protein may dually regulate catalytic activity and membrane association in bacterial and human homologues.

STRUCTURAL AND CATALYTIC ROLE OF TWO CONSERVED TYROSINES IN DELTA‐CLASS GLUTATHIONE S‐TRANSFERASE FROM Locusta migratoria

Glutathione S-transferases (GSTs) are an important family of detoxifying enzymes and play a key role in pesticide resistance in the insect. Tyrosine is essential for its detoxification function. In the present study, two conserved tyrosine residues are located at positions 108 and 116 in H-site of LmGSTD1. To elucidate how the two residues participate in the catalytic process and keeping structural stability, four mutants, Y108A, Y108E, Y116A, and Y116E, were generated. It was found that the four mutants affected the specific activity of LmGSTD1 in various degrees, depending on the types of substrate and reaction mechanism. Steady-state kinetics assay revealed that Y108E and Y116E had a significant influence on GSH-binding ability, which indicates the two tyrosine residues of H-site contribute to topology rearrangement of G-site. Both Y116A and Y116E exhibited lower CDNB-binding affinity, suggesting that Y116 takes part in hydrophobic substrate binding. The thermostability assay, intrinsic, and 8-anilino-1-naphthalenesulfonic acid (ANS) florescence results showed that the two tyrosine residues were involved in regulation of active-site conformation. Finally, homology modeling provided evidence that the two tyrosines in H-site participate in hydrophobic substrate binding. Furthermore, Y108 is closer to the S atom of S-hexylglutathione. In conclusion, the two tyrosines in LmGSTD1 are important residues in both the catalytic process and protein stability.

The structure of monoacylglycerol lipase from Bacillus sp. H257 reveals unexpected conservation of the cap architecture between bacterial and human enzymes

Monoacylglycerol lipases (MGLs) catalyse the hydrolysis of monoacylglycerol into free fatty acid and glycerol. MGLs have been identified throughout all genera of life and have adopted different substrate specificities depending on their physiological role. In humans, MGL plays an integral part in lipid metabolism affecting energy homeostasis, signalling processes and cancer cell progression. In bacteria, MGLs degrade short-chain monoacylglycerols which are otherwise toxic to the organism. We report the crystal structures of MGL from the bacterium Bacillus sp. H257 (bMGL) in its free form at 1.2Å and in complex with phenylmethylsulfonyl fluoride at 1.8Å resolution. In both structures, bMGL adopts an α/β hydrolase fold with a cap in an open conformation. Access to the active site residues, which were unambiguously identified from the protein structure, is facilitated by two different channels. The larger channel constitutes the highly hydrophobic substrate binding pocket with enough room to accommodate monoacylglycerol. The other channel is rather small and resembles the proposed glycerol exit hole in human MGL. Molecular dynamics simulation of bMGL yielded open and closed states of the entrance channel and the glycerol exit hole. Despite differences in the number of residues, secondary structure elements, and low sequence identity in the cap region, this first structure of a bacterial MGL reveals striking structural conservation of the overall cap architecture in comparison with human MGL. Thus it provides insight into the structural conservation of the cap amongst MGLs throughout evolution and provides a framework for rationalising substrate specificities in each organism.

Probing the Activator and Hydrophobic Substrate Binding Sites of PTEN by 31P NMR

PTEN (Phosphatase and Tensin Homolog) antagonizes the PI3K signaling pathway by dephosphorylating PI(3,4,5)P3 at the inositol C-3, thus suppressing cell proliferation. PI(4,5)P2, the product of dephosphorylation, activates this enzyme, possibly by binding to the N terminal region of PTEN (not visible in the crystal structure). A distinct hydrophobic pocket near the active site has also been identified (by MD simulations). We have used 31P NMR as a probe of the spatial location and functional role of these potential phospholipid binding sites. At 242.7 MHz, the linewidths of phosphodiester resonances of phosphoinositides that bind exclusively to the active site (PI, PI(4)P) are broadened consistent with intermediate exchange. Estimated linewidths for the bound activator PI(4,5)P2 are smaller. Other amphiphiles can also be affected by PTEN in a specific fashion, for example, diC7PC which is used as a matrix in assays. This ligand is likely to bind in the hydrophobic pocket identified by simulations. 31P high resolution field cycling NMR, with PTEN spin-labeled on the active site cysteine, was used to determine the distance between the unpaired electron and the phosphorus nuclei of the different bound lipids. Results of these relaxation studies indicate there are discrete sites for both substrate and activator lipids that are also separate from the region occupied by nonspecific lipids such as diC7PC (or a nonionic detergent such as Triton X-100). The distances derived are used to propose a model of PTEN regulation by the activator and matrix lipids.

Crystal structures and kinetic studies of human Kappa class glutathione transferase provide insights into the catalytic mechanism

GSTs (glutathione transferases) are a family of enzymes that primarily catalyse nucleophilic addition of the thiol of GSH (reduced glutathione) to a variety of hydrophobic electrophiles in the cellular detoxification of cytotoxic and genotoxic compounds. GSTks (Kappa class GSTs) are a distinct class because of their unique cellular localization, function and structure. In the present paper we report the crystal structures of hGSTk (human GSTk) in apo-form and in complex with GTX (S-hexylglutathione) and steady-state kinetic studies, revealing insights into the catalytic mechanism of hGSTk and other GSTks. Substrate binding induces a conformational change of the active site from an 'open' conformation in the apo-form to a 'closed' conformation in the GTX-bound complex, facilitating formations of the G site (GSH-binding site) and the H site (hydrophobic substrate-binding site). The conserved Ser(16) at the G site functions as the catalytic residue in the deprotonation of the thiol group and the conserved Asp(69), Ser(200), Asp(201) and Arg(202) form a network of interactions with γ-glutamyl carboxylate to stabilize the thiolate anion. The H site is a large hydrophobic pocket with conformational flexibility to allow the binding of different hydrophobic substrates. The kinetic mechanism of hGSTk conforms to a rapid equilibrium random sequential Bi Bi model.

A structural study of Hypocrea jecorina Cel5A

Interest in generating lignocellulosic biofuels through enzymatic hydrolysis continues to rise as nonrenewable fossil fuels are depleted. The high cost of producing cellulases, hydrolytic enzymes that cleave cellulose into fermentable sugars, currently hinders economically viable biofuel production. Here, we report the crystal structure of a prevalent endoglucanase in the biofuels industry, Cel5A from the filamentous fungus Hypocrea jecorina. The structure reveals a general fold resembling that of the closest homolog with a high-resolution structure, Cel5A from Thermoascus aurantiacus. Consistent with previously described endoglucanase structures, the H. jecorina Cel5A active site contains a primarily hydrophobic substrate binding groove and a series of hydrogen bond networks surrounding two catalytic glutamates. The reported structure, however, demonstrates stark differences between side-chain identity, loop regions, and the number of disulfides. Such structural information may aid efforts to improve the stability of this protein for industrial use while maintaining enzymatic activity through revealing nonessential and immutable regions.

Structural contributions of Delta class glutathione transferase active-site residues to catalysis

GST (glutathione transferase) is a dimeric enzyme recognized for biotransformation of xenobiotics and endogenous toxic compounds. In the present study, residues forming the hydrophobic substrate-binding site (H-site) of a Delta class enzyme were investigated in detail for the first time by site-directed mutagenesis and crystallographic studies. Enzyme kinetics reveal that Tyr111 indirectly stabilizes GSH binding, Tyr119 modulates hydrophobic substrate binding and Phe123 indirectly modulates catalysis. Mutations at Tyr111 and Phe123 also showed evidence for positive co-operativity for GSH and 1-chloro-2,4-dinitrobenzene respectively, strongly suggesting a role for these residues in manipulating subunit-subunit communication. In the present paper we report crystal structures of the wild-type enzyme, and two mutants, in complex with S-hexylglutathione. This study has identified an aromatic 'zipper' in the H-site contributing a network of aromatic pi-pi interactions. Several residues of the cluster directly interact with the hydrophobic substrate, whereas others indirectly maintain conformational stability of the dimeric structure through the C-terminal domain (domain II). The Y119E mutant structure shows major main-chain rearrangement of domain II. This reorganization is moderated through the 'zipper' that contributes to the H-site remodelling, thus illustrating a role in co-substrate binding modulation. The F123A structure shows molecular rearrangement of the H-site in one subunit, but not the other, explaining weakened hydrophobic substrate binding and kinetic co-operativity effects of Phe123 mutations. The three crystal structures provide comprehensive evidence of the aromatic 'zipper' residues having an impact upon protein stability, catalysis and specificity. Consequently, 'zipper' residues appear to modulate and co-ordinate substrate processing through permissive flexing.

Structural characterization of CYP51 from Trypanosoma cruzi and Trypanosoma brucei bound to the antifungal drugs posaconazole and fluconazole.

BackgroundChagas Disease is the leading cause of heart failure in Latin America. Current drug therapy is limited by issues of both efficacy and severe side effects. Trypansoma cruzi, the protozoan agent of Chagas Disease, is closely related to two other major global pathogens, Leishmania spp., responsible for leishmaniasis, and Trypansoma brucei, the causative agent of African Sleeping Sickness. Both T. cruzi and Leishmania parasites have an essential requirement for ergosterol, and are thus vulnerable to inhibitors of sterol 14α-demethylase (CYP51), which catalyzes the conversion of lanosterol to ergosterol. Clinically employed anti-fungal azoles inhibit ergosterol biosynthesis in fungi, and specific azoles are also effective against both Trypanosoma and Leishmania parasites. However, modification of azoles to enhance efficacy and circumvent potential drug resistance has been problematic for both parasitic and fungal infections due to the lack of structural insights into drug binding.Methodology/Principal FindingsWe have determined the crystal structures for CYP51 from T. cruzi (resolutions of 2.35 Å and 2.27 Å), and from the related pathogen T. brucei (resolutions of 2.7 Å and 2.6 Å), co-crystallized with the antifungal drugs fluconazole and posaconazole. Remarkably, both drugs adopt multiple conformations when binding the target. The fluconazole 2,4-difluorophenyl ring flips 180° depending on the H-bonding interactions with the BC-loop. The terminus of the long functional tail group of posaconazole is bound loosely in the mouth of the hydrophobic substrate binding tunnel, suggesting that the major contribution of the tail to drug efficacy is for pharmacokinetics rather than in interactions with the target.Conclusions/SignificanceThe structures provide new insights into binding of azoles to CYP51 and mechanisms of potential drug resistance. Our studies define in structural detail the CYP51 therapeutic target in T. cruzi, and offer a starting point for rationally designed anti-Chagasic drugs with improved efficacy and reduced toxicity.

Origins of catalysis by computationally designed retroaldolase enzymes

We have investigated recently reported computationally designed retroaldolase enzymes with the goal of understanding the extent and the origins of their catalytic power. Direct comparison of the designed enzymes to primary amine catalysts in solution revealed a rate acceleration of 10(5)-fold for the most active of the designed retroaldolases. Through pH-rate studies of the designed retroaldolases and evaluation of a Brønsted correlation for a series of amine catalysts, we found that lysine pK(a) values are shifted by 3-4 units in the enzymes but that the catalytic contributions from the shifted pK(a) values are estimated to be modest, about 10-fold. For the most active of the reported enzymes, we evaluated the catalytic contribution of two other design components: a motif intended to stabilize a bound water molecule and hydrophobic substrate binding interactions. Mutational analysis suggested that the bound water motif does not contribute to the rate acceleration. Comparison of the rate acceleration of the designed substrate relative to a minimal substrate suggested that hydrophobic substrate binding interactions contribute around 10(3)-fold to the enzymatic rate acceleration. Altogether, these results suggest that substrate binding interactions and shifting the pK(a) of the catalytic lysine can account for much of the enzyme's rate acceleration. Additional observations suggest that these interactions are limited in the specificity of placement of substrate and active site catalytic groups. Thus, future design efforts may benefit from a focus on achieving precision in binding interactions and placement of catalytic groups.

Structural and functional aspects of the multidrug efflux pump AcrB

The tripartite efflux system AcrA/AcrB/TolC is the main pump in Escherichia coli for the efflux of multiple antibiotics, dyes, bile salts and detergents. The inner membrane component AcrB is central to substrate recognition and energy transduction and acts as a proton/drug antiporter. Recent structural studies show that homotrimeric AcrB can adopt different monomer conformations representing consecutive states in an allosteric functional rotation transport cycle. The conformational changes create an alternate access drug transport tunnel including a hydrophobic substrate binding pocket in one of the cycle intermediates.

The crystal structure of the ring‐hydroxylating dioxygenase from Sphingomonas CHY‐1

The ring-hydroxylating dioxygenase (RHD) from Sphingomonas CHY-1 is remarkable due to its ability to initiate the oxidation of a wide range of polycyclic aromatic hydrocarbons (PAHs), including PAHs containing four- and five-fused rings, known pollutants for their toxic nature. Although the terminal oxygenase from CHY-1 exhibits limited sequence similarity with well characterized RHDs from the naphthalene dioxygenase family, the crystal structure determined to 1.85 A by molecular replacement revealed the enzyme to share the same global alpha(3)beta(3) structural pattern. The catalytic domain distinguishes itself from other bacterial non-heme Rieske iron oxygenases by a substantially larger hydrophobic substrate binding pocket, the largest ever reported for this type of enzyme. While residues in the proximal region close to the mononuclear iron atom are conserved, the central region of the catalytic pocket is shaped mainly by the side chains of three amino acids, Phe350, Phe404 and Leu356, which contribute to the rather uniform trapezoidal shape of the pocket. Two flexible loops, LI and LII, exposed to the solvent seem to control the substrate access to the catalytic pocket and control the pocket length. Compared with other naphthalene dioxygenases residues Leu223 and Leu226, on loop LI, are moved towards the solvent, thus elongating the catalytic pocket by at least 2 A. An 11 A long water channel extends from the interface between the alpha and beta subunits to the catalytic site. The comparison of these structures with other known oxygenases suggests that the broad substrate specificity presented by the CHY-1 oxygenase is primarily due to the large size and particular topology of its catalytic pocket and provided the basis for the study of its reaction mechanism.

The catalytic pocket of the ring-hydroxylating dioxygenase from Sphingomonas CHY-1

Ring-hydroxylating dioxygenases are multicomponent bacterial enzymes that catalyze the first step in the oxidative degradation of aromatic hydrocarbons. The dioxygenase from Sphingomonas CHY-1 is unique in that it can oxidize a wide range of polycyclic aromatic hydrocarbons (PAHs). With a crystal structure similar to that of the seven other known dioxygenases, its catalytic domain features the largest hydrophobic substrate binding cavity characterized so far. Molecular modeling studies indicated that the catalytic cavity is large enough to accommodate a five-ring benzo[a]pyrene molecule. The predicted positions of this and other PAHs in the substrate binding pocket are consistent with the product regio- and stereo-selectivity of the enzyme.

Structural Basis of the Suppressed Catalytic Activity of Wild-type Human Glutathione Transferase T1-1 Compared to its W234R Mutant

The crystal structures of wild-type human theta class glutathione-S-transferase (GST) T1-1 and its W234R mutant, where Trp234 was replaced by Arg, were solved both in the presence and absence of S-hexyl-glutathione. The W234R mutant was of interest due to its previously observed enhanced catalytic activity compared to the wild-type enzyme. GST T1-1 from rat and mouse naturally contain Arg in position 234, with correspondingly high catalytic efficiency. The overall structure of GST T1-1 is similar to that of GST T2-2, as expected from their 53% sequence identity at the protein level. Wild-type GST T1-1 has the side-chain of Trp234 occupying a significant portion of the active site. This bulky residue prevents efficient binding of both glutathione and hydrophobic substrates through steric hindrance. The wild-type GST T1-1 crystal structure, obtained from co-crystallization experiments with glutathione and its derivatives, showed no electron density for the glutathione ligand. However, the structure of GST T1-1 mutant W234R showed clear electron density for S-hexyl-glutathione after co-crystallization. In contrast to Trp234 in the wild-type structure, the side-chain of Arg234 in the mutant does not occupy any part of the substrate-binding site. Instead, Arg234 is pointing in a different direction and, in addition, interacts with the carboxylate group of glutathione. These findings explain our earlier observation that the W234R mutant has a markedly improved catalytic activity with most substrates tested to date compared to the wild-type enzyme. GST T1-1 catalyzes detoxication reactions as well as reactions that result in toxic products, and our findings therefore suggest that humans have gained an evolutionary advantage by a partially disabled active site.

An Electron-sharing Network Involved in the Catalytic Mechanism Is Functionally Conserved in Different Glutathione Transferase Classes

In Anopheles dirus glutathione transferase D3-3, there are electrostatic interactions between the negatively charged glutamyl alpha-carboxylate group of glutathione, the positively charged Arg-66, and the negatively charged Asp-100. This ionic interaction is stabilized by a network of hydrogen bonds from Ser-65, Thr-158, Thr-162, and a conserved water-mediated contact. This alternating ionic bridge interaction between negatively and positively charged residues stabilized by a network of hydrogen bonding we have named an electron-sharing network. We show that the electron-sharing network assists the glutamyl alpha-carboxylate of glutathione to function as a catalytic base accepting the proton from the thiol group forming an anionic glutathione, which is a crucial step in the glutathione transferase (GST) catalysis. Kinetic studies demonstrate that the mutation of electron-sharing network residues results in a decreased ability to lower the pKa of the thiol group of glutathione. Although the residues that contribute to the electron-sharing network are not conserved in the primary sequence, structural characterizations indicate that the presence of the network can be mapped to the same region in all GST classes. A structural diversification but functional conservation suggests a significant role for the electron-sharing network in catalysis as the purpose was maintained during the divergent evolution of GSTs. This network appears to be a functionally conserved motif that contributes to the "base-assisted deprotonation" model suggested to be essential for the glutathione ionization step of the catalytic mechanism.

Structure-based mutational analysis of the active site residues of d-hydantoinase

We previously proposed the hydrophobic and bulky residues of the three loops, designated stereochemistry gate loops (SGLs), to constitute a hydrophobic substrate binding pocket of d-hydantoinase from Bacillus stearothermophilus SD1. Simulation of substrate binding in the active site of d-hydantoinase and sequence alignment of various d-hydantoinases revealed the critical hydrophobic residues closely located around the exocyclic substituent of substrate. To evaluate the roles of these residues in substrate binding pocket, site-directed mutagenesis was performed specifically for Leu 65, Tyr 155, and Phe 159. When Tyr 155 was mutated to Phe and Glu, both mutants Y155F and Y155E were totally inactive for nonsubstituted hydantoin and d,l-5-hydroxyphenyl hydantoin (HPH), which indicates that Tyr 155 is involved in substrate binding via a hydrogen bond with the hydantoinic ring. Furthermore, replacement of the hydrophobic residues Leu 65 and Phe 159 with Glu, a charged amino acid, resulted in a significant decrease in activity for nonsubstituted hydantoin, but not for HPH. The K cat values of both mutants for nonsubstituted hydantoin also severely decreased, but a slight change in the K cat values was observed towards HPH. These results suggest that the hydrophobic residues in SGLs play an essential role in substrate binding, and differentially interact according to the property of the exocyclic substituent.

Solvent interaction of a Hsp70 chaperone substrate-binding domain investigated with water-NOE NMR experiments.

The interaction of solvent of the substrate binding domain of the bacterial heat shock 70 chaperone protein DnaK was studied in its apo form and with bound hydrophobic substrate peptide, using refined nuclear magnetic resonance experiments. Distinct differences between the two states of the protein were observed. According to our data, the apo form interacts more extensively with solvent than the peptide-bound form. Significantly, the open hydrophobic substrate binding cleft of DnaK in the apo form is found to contain several molecules of water which are displaced by the binding of the hydrophobic substrate, the peptide NRLLLTG. The solvent in the hydrophobic cleft has a residence time longer than 400 ps. It is predicted that the displacement of this trapped water must contribute to the binding free energy of the natural hydrophobic substrates of this class of protein-folding chaperone proteins.

Hydrophobic effects on rates and substrate selectivities in polymeric transaminase mimics.

The amination of ketoacids to amino acids by pyridoxamine is greatly accelerated when the pyridoxamine is covalently linked to polyethylenimine carrying N-methyl and N-lauryl groups. Michaelis-Menten kinetics is seen with all substrates, from which the effect of the lauryl groups and the methyl groups can be determined with respect to the strength of binding of the substrate and the rate constant k2 within the complex. The polyamine catalyzes the reaction using acid and base groups, the lauryl groups increase k2 by producing a nonpolar medium in which the reaction occurs, and the lauryl groups promote binding of hydrophobic substrates. The result is that the amination of indolepyruvic acid to produce tryptophan is accelerated by 240000-fold.

Variation in aspects of cysteine proteinase catalytic mechanism deduced by spectroscopic observation of dithioester intermediates, kinetic analysis and molecular dynamics simulations

The possibility of a slow post-acylation conformational change during catalysis by cysteine proteinases was investigated by using a new chromogenic substrate, N-acetyl-Phe-Gly methyl thionoester, four natural variants (papain, caricain, actinidin and ficin), and stopped-flow spectral analysis to monitor the pre-steady state formation of the dithioacylenzyme intermediates and their steady state hydrolysis. The predicted reversibility of acylation was demonstrated kinetically for actinidin and ficin, but not for papain or caricain. This difference between actinidin and papain was investigated by modelling using QUANTA and CHARMM. The weaker binding of hydrophobic substrates, including the new thionoester, by actinidin than by papain may not be due to the well-known difference in their S2-subsites, whereby that of actinidin in the free enzyme is shorter due to the presence of Met211. Molecular dynamics simulation suggests that during substrate binding the sidechain of Met211 moves to allow full access of a Phe sidechain to the S2-subsite. The highly anionic surface of actinidin may contribute to the specificity difference between papain and actinidin. During subsequent molecular dynamics simulations the P1 product, methanol, diffuses rapidly (over < 8ps) out of papain and caricain but ‘lingers’ around the active centre of actinidin. Uniquely in actinidin, an Asp142–Lys145 salt bridge allows formation of a cavity which appears to constrain diffusion of the methanol away from the catalytic site. The cavity then undergoes large scale movements (over 4.8 Å) in a highly correlated manner, thus controlling the motions of the methanol molecule. The changes in this cavity that release the methanol might be those deduced kinetically.

Variation in aspects of cysteine proteinase catalytic mechanism deduced by spectroscopic observation of dithioester intermediates, kinetic analysis and molecular dynamics simulations.

The possibility of a slow post-acylation conformational change during catalysis by cysteine proteinases was investigated by using a new chromogenic substrate, N-acetyl-Phe-Gly methyl thionoester, four natural variants (papain, caricain, actinidin and ficin), and stopped-flow spectral analysis to monitor the pre-steady state formation of the dithioacylenzyme intermediates and their steady state hydrolysis. The predicted reversibility of acylation was demonstrated kinetically for actinidin and ficin, but not for papain or caricain. This difference between actinidin and papain was investigated by modelling using QUANTA and CHARMM. The weaker binding of hydrophobic substrates, including the new thionoester, by actinidin than by papain may not be due to the well-known difference in their S2-subsites, whereby that of actinidin in the free enzyme is shorter due to the presence of Met211. Molecular dynamics simulation suggests that during substrate binding the sidechain of Met211 moves to allow full access of a Phe sidechain to the S2-subsite. The highly anionic surface of actinidin may contribute to the specificity difference between papain and actinidin. During subsequent molecular dynamics simulations the P1 product, methanol, diffuses rapidly (over<8 ps) out of papain and caricain but 'lingers' around the active centre of actinidin. Uniquely in actinidin, an Asp142-Lys145 salt bridge allows formation of a cavity which appears to constrain diffusion of the methanol away from the catalytic site. The cavity then undergoes large scale movements (over 4.8 A) in a highly correlated manner, thus controlling the motions of the methanol molecule. The changes in this cavity that release the methanol might be those deduced kinetically.