Molecular dynamics simulations of the acyl-enzyme and the tetrahedral intermediate in the deacylation step of serine proteases

Despite the availability of many experimental data and some modeling studies, questions remain as to the precise mechanism of the serine proteases. Here we report molecular dynamics simulations on the acyl-enzyme complex and the tetrahedral intermediate during the deacylation step in elastase catalyzed hydrolysis of a simple peptide. The models are based on recent crystallographic data for an acyl-enzyme intermediate at pH 5 and a time-resolved study on the deacylation step. Simulations were carried out on the acyl enzyme complex with His-57 in protonated (as for the pH 5 crystallographic work) and deprotonated forms. In both cases, a water molecule that could provide the nucleophilic hydroxide ion to attack the ester carbonyl was located between the imidazole ring of His-57 and the carbonyl carbon, close to the hydrolytic position assigned in the crystal structure. In the "neutral pH" simulations of the acyl-enzyme complex, the hydrolytic water oxygen was hydrogen bonded to the imidazole ring and the side chain of Arg-61. Alternative stable locations for water in the active site were also observed. Movement of the His-57 side-chain from that observed in the crystal structure allowed more solvent waters to enter the active site, suggesting that an alternative hydrolytic process directly involving two water molecules may be possible. At the acyl-enzyme stage, the ester carbonyl was found to flip easily in and out of the oxyanion hole. In contrast, simulations on the tetrahedral intermediate showed no significant movement of His-57 and the ester carbonyl was constantly located in the oxyanion hole. A comparison between the simulated tetrahedral intermediate and a time-resolved crystallographic structure assigned as predominantly reflecting the tetrahedral intermediate suggests that the experimental structure may not precisely represent an optimal arrangement for catalysis in solution. Movement of loop residues 216-223 and P3 residue, seen both in the tetrahedral simulation and the experimental analysis, could be related to product release. Furthermore, an analysis of the geometric data obtained from the simulations and the pH 5 crystal structure of the acyl-enzyme suggests that since His-57 is protonated, in some aspects, this crystal structure resembles the tetrahedral intermediate.     

Alkyl isocyanates as active site-specific reagents for serine proteases. Location of alkyl binding site in chymotrypsin by X-ray diffraction.

The structure of octylcarbamoyl-alpha-chymotrypsin to a resolution of 3.0 A is described. The n-octyl side chain of the active site directed irreversible inactivator octyl isocyanate is bound exclusively in the hydrophobic substrate binding pocket. The n-octyl isocyanate forms a planar urethane bond with the Ser-195 Ogamma and extends approximately 1 A deeper into the hydrophobic pocket than the indolyl group of indoleacryloyl-alpha-chymotrypsin (Henderson, R. (1970), J. Mol. Biol. 54, 341). All the structural changes are essentially identical with those observed in indoleacryloyl-alpha-chymotrypsin including the observation of a hydrogen bonded water molecule between the carbonyl oxygen of the octylcarbamoyl group and the imidazole group of His-57. The observed mode of n-octyl alkyl binding to chymotrypsin is consistent with the hypothesis proposed earlier (Brown, W. E. and Wold, F. (1973), Biochemistry 12, 828).     

Structure of the complex formed by bovine trypsin and bovine pancreatic trypsin inhibitor. Crystal structure determination and stereochemistry of the contact region

The structure of the complex of bovine trypsin and bovine pancreatic trypsin inhibitor has been determined by crystal structure analysis at 2.8 Å resolution. The structure is closely similar to the model predicted from the structures of the components. The complex is a tetrahedral adduct with a covalent bond between the carbonyl carbon of Lys-15I of the inhibitor and the γ-oxygen of Ser-195 of the enzyme. The imidazole of His-57 is hydrogen-bonded to Asp-102 and the bound seryl γ-oxygen in accord with the histidine being charged. The negatively charged carbonyl oxygen of Lys-15I forms two hydrogen bonds with the amide nitrogens of Gly-193 and Ser-195. Protonation of the leaving group N-H of Ala-16I to form an acyl-complex requires a conformational change of the imidazole of His-57. The tetrahedral adduct is further stabilized by hydrogen bonds between groups at the leaving group side and inhibitor and enzyme, which would be weakened in the acyl-enzyme. The kinetic data of inhibitor-enzyme interaction are reconciled with the structural model, and relations between enzyme-inhibitor interaction and productive enzyme-substrate interaction are proposed.     

Stereochemical aspects of peptide hydrolysis catalyzed by serine proteases of the chymotrypsin type

The steric course of peptide hydrolysis catalyzed by serine proteases has been studied on the basis of the available, extensive structural data and taking into account the stereoelectronic theory of Deslongchamps (Heterocycles, 7, 1271 (1977)). These studies allowed elucidation of the structure of intermediates, in particular of the tetrahedral intermediate, and of the main structural events taking place during catalysis. They reveal a difficulty inherent in the generally accepted mechanism of peptide hydrolysis: protonation of the leaving nitrogen in the configuration arising from nucleophilic attack of Ser-195 on the carbonyl carbon cannot take place internally from His-57. Two alternative mechanisms are discussed which are compatible with all implications of the stereoelectronic theory. The main features of the more probable mechanism are: (i) a conformational change allowing the imidazole ring of His-57 to occupy two distinct positions; in one position a proton is abstracted from O^γ of Ser-195, and in the other this proton is donated to the leaving nitrogen; (ii) a configurational change (inversion) of the pyramidal leaving nitrogen reorienting the lone-pair orbital developed during nucleophilic attack; in one orientation CO bond breaking, and in the other CN bond breaking, is allowed. This inversion process confers on the nitrogen the property of a switch controlling the breakdown of the tetrahedral intermediate.     

Crystal structure of a mini-intein reveals a conserved catalytic module involved in side chain cyclization of asparagine during protein splicing

We have determined the crystal structure of a 154-residue intein derived from the dnaB gene of Synechocystis sp. strain PCC6803 and refined it to a 2.0-A resolution. The x-ray structure suggests that this intein possesses two catalytic sites that appear to be separately responsible for splicing and cleavage of the N- and C-terminal scissile bonds. The conserved intein block F residues are the important components of a catalytic site for side chain cyclization of the last intein residue, Asn-154. The data suggest that the imidazole ring of His-143 is involved in the activation of the side chain Ndelta atom of Asn-154, leading to a nucleophilic attack on the carbonyl carbon of Asn-154. Substitution of His-143 with Ala or Gln resulted in the inhibition of C-terminal cleavage. His-153, Asp-136, and a water molecule appear to constitute an oxyanion binding site by contacting the carbonyl oxygen of Asn-154 to stabilize the transition state. The structure and mutagenesis data also support that the close contact between the hydroxyl groups of Thr-138 and Ser-155, whose side chain participates in an S --> O acyl shift, plays an important role in the nucleophile orientation. Our structural modeling suggests that this catalytic module is conserved in the C-terminal subdomains of inteins from diverse organisms.     

Crystal structure at 1.63 A resolution of the native form of porcine beta-trypsin: revealing an acetate ion binding site and functional water network

The active center of a serine protease is the catalytic triad composed of His-57, Ser-195 and Asp-102. The existing crystal structure data on serine proteases have not fully answered a number of fundamental questions relating to the catalytic activity of serine proteases. The new high resolution native porcine beta-trypsin (BPT) structure is aimed at extending the knowledge on the conformation of the active site and the ordered water structure within and around the active site. The crystal structure of BPT has been determined at 1.63 A resolution. An acetate ion bound at the active site of a trypsin molecule by both classical hydrogen bonds and C-HellipsisO hydrogen bonds has been identified for the first time. A large network of water molecules extending from the recognition amino acid Asp-184 to the entry of the active site has been observed in the BPT structure. A detailed comparison with inhibitor complexes and autolysates indicates that the sulfate ion and the acetate ion bind at the same site of the trypsin molecule. The Ser-195 Cbeta-Ogamma-His-57 Nepsilon angle in the catalytic triad of BPT is intermediate between the corresponding values of the complex and native structure due to acetate ion binding. The network of waters from the recognition amino acid to the active site entry is probably the first ever complete picture of functional waters around the active site. Structural comparisons show that the functional waters involved in the binding of small molecule inhibitors and protease inhibitors are distinctly different.     

Structural diversity within the mononuclear and binuclear active sites of N-acetyl-D-glucosamine-6-phosphate deacetylase

NagA catalyzes the hydrolysis of N-acetyl-d-glucosamine-6-phosphate to d-glucosamine-6-phosphate and acetate. X-ray crystal structures of NagA from Escherichia coli were determined to establish the number and ligation scheme for the binding of zinc to the active site and to elucidate the molecular interactions between the protein and substrate. The three-dimensional structures of the apo-NagA, Zn-NagA, and the D273N mutant enzyme in the presence of a tight-binding N-methylhydroxyphosphinyl-d-glucosamine-6-phosphate inhibitor were determined. The structure of the Zn-NagA confirms that this enzyme binds a single divalent cation at the beta-position in the active site via ligation to Glu-131, His-195, and His-216. A water molecule completes the ligation shell, which is also in position to be hydrogen bonded to Asp-273. In the structure of NagA bound to the tight binding inhibitor that mimics the tetrahedral intermediate, the methyl phosphonate moiety has displaced the hydrolytic water molecule and is directly coordinated to the zinc within the active site. The side chain of Asp-273 is positioned to activate the hydrolytic water molecule via general base catalysis and to deliver this proton to the amino group upon cleavage of the amide bond of the substrate. His-143 is positioned to help polarize the carbonyl group of the substrate in conjunction with Lewis acid catalysis by the bound zinc. The inhibitor is bound in the alpha-configuration at the anomeric carbon through a hydrogen bonding interaction of the hydroxyl group at C-1 with the side chain of His-251. The phosphate group of the inhibitor attached to the hydroxyl at C-6 is ion paired with Arg-227 from the adjacent subunit. NagA from Thermotoga maritima was shown to require a single divalent cation for full catalytic activity.     

Stabilization of the imidazole ring of His-195 at the active site of chloramphenicol acetyltransferase

The imidazole of His-195 plays an essential role in the proposed general base mechanism of chloramphenicol acetyltransferase (CAT). The structure of the binary complex of CATIII and chloramphenicol suggests that two unusual interactions might determine the conformation of the side chain of His-195: (i) an intraresidue hydrogen bond between its main chain carbonyl and the protonated N delta 1 of the imidazole ring and (ii) face-to-face van der Waals contact between the His-195 imidazole group and the aromatic side chain of Tyr-25. Tyr-25 also makes a hydrogen bond, via its phenolic hydroxyl, to the carbonyl oxygen of the substrate chloramphenicol. Replacement of Tyr-25 of CATIII by phenylalanine results in a modest increase in the Km for chloramphenicol (from 11.6 to 14.6 microM) and a 2-fold fall in kcat (599 to 258 s-1), indicative of a free energy contribution to transition state binding of 0.6 kcal mol-1 for the hydrogen bond between Tyr-25 and chloramphenicol. In contrast, substitution of Tyr-25 by alanine yields an enzyme that is dramatically impaired in its ability to bind chloramphenicol (Km = 173 microM). As kcat for Ala-25 CAT is also reduced (130 s-1), the loss of the aryl group results in a 69-fold decrease in kcat/Km, corresponding to a free energy contribution to binding and catalysis of 2.5 kcal mol-1. In addition to the loss of the hydrogen bond between Tyr-25 and chloramphenicol, the loss of substrate affinity in Ala-25 CAT may be a direct consequence of reduced hydrophobicity of the chloramphenicol-binding site and/or the loss of critical constraints on the precise conformation of the catalytic imidazole. However, as with wild type CAT, inactivation of Ala-25 CAT by the affinity reagent 3-(bromoacetyl) chloramphenicol is accompanied by modification solely at N epsilon 2 of His-195. Hence, the results demonstrate that tautomeric stabilization of the imidazole ring persists in the absence of van der Waals interactions with the side chain of Tyr-25, probably as a consequence of hydrogen bonding between the protonated N delta 1 and the carbonyl oxygen of His-195.     

On the mechanism of action of methyl chymotrypsin

The properties of a-chymotrypsin methylated at histidine-57 were examined to explain the mechanism of this enzyme which is about 10⁵ times less active than chymotrypsin. Studies on the protein showed (i) an alteration in the acyl and leaving group specificity, (ii) decreased binding of some protein protease inhibitors by methyl chymotrypsin, (iii) lack of dimerization of methyl chymotrypsin at low pH, (iv) decreased stability of methyl chymotrypsin in urea, (v) a larger solvent deuterium isotope effect with methyl chymotrypsin, and (vi) decreased binding of a tetrahedral intermediate analog to methyl chymotrypsin. These properties suggest that while only subtle alterations occur in the active site upon methylation of His-57, the transition state and the tetrahedral intermediate are destabilized but not to the same extent. General base catalysis remains an integral feature of the hydrolytic mechanism of the modified chymotrypsin, and the base appears to be the methylated nitrogen of the imidazole moiety of His-57.     

X-ray crystallographic and site-directed mutagenesis analysis of the mechanism of Schiff-base formation in phosphonoacetaldehyde hydrolase catalysis

Phosphonoacetaldehyde hydrolase (phosphonatase) catalyzes the hydrolytic P-C bond cleavage of phosphonoacetaldehyde (Pald) to form orthophosphate and acetaldehyde. The reaction proceeds via a Schiff-base intermediate formed between Lys-53 and the Pald carbonyl. The x-ray crystal structures of the wild-type phosphonatase complexed with Mg(II) alone or with Mg(II) plus vinylsulfonate (a phosphonoethylenamine analog) were determined to 2.8 and 2.4 A, respectively. These structures were used to determine the identity and positions of active site residues surrounding the Lys-53 ammonium group and the Pald carbonyl. These include Cys-22, His-56, Tyr-128, and Met-49. Site-directed mutagenesis was then employed to determine whether or not these groups participate in catalysis. Based on rate contributions, Tyr-128 and Cys-22 were eliminated as potential catalytic groups. The Lys-53 epsilon-amino group, positioned for reaction with the Pald carbonyl, forms a hydrogen bond with water 120. Water 120 is also within hydrogen bond distance of an imidazole nitrogen of His-56 and the sulfur atom of Met-49. Kinetic constants for mutants indicated that His-56 (1000-fold reduction in k(cat)/K(m) upon Ala substitution) and Met-49 (17,000-fold reduction in k(cat)/K(m) upon Leu substitution) function in catalysis of Schiff-base formation. Based on these results, it is proposed that a network of hydrogen bonds among Lys-53, water 120, His-56, and Met-49 facilitate proton transfer from Lys-53 to the carbinolamine intermediate. Comparison of the vinylsulfonate complex versus unliganded structures indicated that association of the cap and core domains is essential for the positioning of the Lys-53 for attack at the Pald carbonyl and that substrate binding at the core domain stabilizes cap domain binding.     

Two crystal structures of the leupeptin-trypsin complex.

Three-dimensional structures of trypsin with the reversible inhibitor leupeptin have been determined in two different crystal forms. The first structure was determined at 1.7 A resolution with R-factor = 17.7% in the trigonal crystal space group P3(1)21, with unit cell dimensions of a = b = 55.62 A, c = 110.51 A. The second structure was determined at a resolution of 1.8 A with R-factor = 17.5% in the orthorhombic space group P2(1)2(1)2(1), with unit cell dimensions of a = 63.69 A, b = 69.37 A, c = 63.01 A. The overall protein structure is very similar in both crystal forms, with RMS difference for main-chain atoms of 0.27 A. The leupeptin backbone forms four hydrogen bonds with trypsin and a fifth hydrogen bond interaction is mediated by a water molecule. The aldehyde carbonyl of leupeptin forms a covalent bond of 1.42 A length with side-chain oxygen of Ser-195 in the active site. The reaction of trypsin with leupeptin proceeds through the formation of stable tetrahedral complex in which the hemiacetal oxygen atom is pointing out of the oxyanion hole and forming a hydrogen bond with His-57.     

Cleavage of beta-lactone ring by serine protease. Mechanistic implications

Both enantiomers of 3-benzyl-2-oxetanone (1) were found to be slowly hydrolyzed substrates of alpha-chymotrypsin having k(cat) values of 0.134+/-0.008 and 0.105+/-0.004 min(-1) for (R)-1 and (S)-1, respectively, revealing that alpha-CT is virtually unable to differentiate the enantiomers in the hydrolysis of 1. The initial step to form the acyl-enzyme intermediate by the attack of Ser-195 hydroxyl on the beta-lactone ring at the 2-position in the hydrolysis reaction may not be enzymatically driven, but the relief of high ring strain energy of beta-lactone may constitute a major driving force. The deacylation step is also attenuated, which is possibly due to the hydrogen bond that would be formed between the imidazole nitrogen of His-57 and the hydroxyl group generated during the acylation in the case of (R)-1, but in the alpha-CT catalyzed hydrolysis of (S)-1 the imidazole nitrogen may form a hydrogen bond with the ester carbonyl oxygen.     

Identification of important residues in diketoreductase from Acinetobacter baylyi by molecular modeling and site-directed mutagenesis

Diketoreductase (DKR) from Acinetobacter baylyi exhibits a unique property of double reduction of a β, δ-diketo ester with excellent stereoselectivity, which can serve as an efficient biocatalyst for the preparation of an important chiral intermediate for cholesterol lowering statin drugs. Taken the advantage of high homology between DKR and human heart 3-hydroxyacyl-CoA dehydrogenase (HAD), a molecular model was created to compare the tertiary structures of DKR and HAD. In addition to the possible participation of His-143 in the enzyme catalysis by pH profile, three key amino acid residues, Ser-122, His-143 and Glu-155, were identified and mutated to explore the possibility of involving in the catalytic process. The catalytic activities for mutants S122A/C, H143A/K and E155Q were below detectable level, while their binding affinities to the diketo ester substrate and cofactor NADH did not change obviously. The experimental results were further supported by molecular docking, suggesting that Ser-122 and His-143 were essential for the proton transfer to the carbonyl functional groups of the substrate. Moreover, Glu-155 was crucial for maintaining the proper orientation and protonation of the imidazole ring of His-143 for efficient catalysis.     

Quantum chemistry analysis of the mechanism of action of proteolytic enzymes. III. Proton transport in serine proteases

The model system for the proton transfer on the amide atom of the substrate leaving group based on the existence of "charge relay system" in the serine type proteases was analysed by the CNDO/2 method. The unfitness of this model to explain the action mechanism of serine proteases was shown. The model system for proton transfer with the water molecule as the intermediate acceptor of the Ser-195 proton was suggested and analysed by the same method. The acylation activation barrier of this system was shown to localize on the stage of synchronous transfer of the Ser-195 alcoholic proton and the water molecule proton hydrogen bound to the His-57 N epsilon 2-atom on the water molecule oxygen atom and the N epsilon 2-atom, respectively. The protonation of substrate in the case of the model system with the water molecule as the intermediate acceptor of proton was demonstrated to begin before the completion of the tetrahedral intermediate substance and the protonated from of the tetrahedral intermediate was shown to form only. A hypothesis considering the role of this water molecule as the nucleophilic reagent on the deacylation stage is presented.     

Catalytic roles for two water bridged residues (Asp-98 and His-255) in the active site of copper-containing nitrite reductase

Two active site residues, Asp-98 and His-255, of copper-containing nitrite reductase (NIR) from Alcaligenes faecalis have been mutated to probe the catalytic mechanism. Three mutations at these two sites (D98N, H255D, and H255N) result in large reductions in activity relative to native NIR, suggesting that both residues are involved intimately in the reaction mechanism. Crystal structures of these mutants have been determined using data collected to better than 1. 9-A resolution. In the native structure, His-255 Nepsilon2 forms a hydrogen bond through a bridging water molecule to the side chain of Asp-98, which also forms a hydrogen bond to a water or nitrite oxygen ligated to the active site copper. In the D98N mutant, reorientation of the Asn-98 side chain results in the loss of the hydrogen bond to the copper ligand water, consistent with a negatively charged Asp-98 directing the binding and protonation of nitrite in the native enzyme. An additional solvent molecule is situated between residues 255 and the bridging water in the H255N and H255D mutants and likely inhibits nitrite binding. The interaction of His-255 with the bridging water appears to be necessary for catalysis and may donate a proton to reaction intermediates in addition to Asp-98.     

Reaction of porcine pancreatic elastase with 7-substituted 3-alkoxy-4-chloroisocoumarins: design of potent inhibitors using the crystal structure of the complex formed with 4-chloro-3-ethoxy-7-guanidinoisocoumarin

The crystal structure of the acyl enzyme formed upon inhibition of porcine pancreatic elastase (PPE) by 4-chloro-3-ethoxy-7-guanidinoisocoumarin has been determined at a 1.85-A effective resolution. The chlorine atom is still present in this acyl enzyme, in contrast to the previously reported structure of the 7-amino-4-chloro-3-methoxyisocoumarin-PPE complex where the chlorine atom has been replaced by an acetoxy group. The guanidino group forms hydrogen bonds with the carbonyl group and side-chain hydroxyl group of Thr-41, and the acyl carbonyl group has been twisted out of the oxyanion hole. Molecular modeling indicates that the orientation of the initial Michaelis enzyme-inhibitor complex is quite different from that of the acyl enzyme since simple reconstruction of the isocoumarin ring would result in unfavorable interactions with Ser-195 and His-57. Molecular models were used to design a series of new 7-(alkylureido)- and 7-(alkylthioureido)-substituted derivatives of 3-alkoxy-7-amino-4-chloroisocoumarin as PPE inhibitors. All the 3-ethoxyisocoumarins were better inhibitors than those in the 3-methoxy series due to better interactions with the S1 pocket of PPE. The best ureido inhibitor also contained a tert-butylureido group at the 7-position of the isocoumarin. Due to a predicted interaction with a small hydrophobic pocket on the surface of PPE, this isocoumarin and a related phenylthioureido derivative are among the best irreversible inhibitors thus far reported for PPE (kobs/[I] = 8100 M-1 s-1 and 12,000 M-1 s-1). Kinetic studies of the stability of enzyme-inhibitor complexes suggest that many isocoumarins are alkylating the active site histidine at pH 7.5 via a quinone imine methide intermediate, while at pH 5.0, the predominant pathway appears to be simple formation of a stable acyl enzyme derivative.     

Structures of trihydroxynaphthalene reductase-fungicide complexes: implications for structure-based design and catalysis

BACKGROUND: Trihydroxynaphthalene reductase catalyzes two intermediate steps in the fungal melanin biosynthetic pathway. The enzyme, a typical short-chain dehydrogenase, is the biochemical target of three commercial fungicides. The fungicides bind preferentially to the NADPH form of the enzyme. RESULTS: Three X-ray structures of the Magnaporthe grisea enzyme complexed with NADPH and two commercial and one experimental fungicide were determined at 1.7 A (pyroquilon), 2.0 A (2,3-dihydro-4-nitro-1H-inden-1-one, 1), and 2.1 A (phthalide) resolutions. The chemically distinct inhibitors occupy similar space within the enzyme's active site. The three inhibitors share hydrogen bonds with the side chain hydroxyls of Ser-164 and Tyr-178 via a carbonyl oxygen (pyroquilon and 1) or via a carbonyl oxygen and a ring oxygen (phthalide). Active site residues occupy similar positions among the three structures. A buried water molecule that is hydrogen bonded to the NZ nitrogen of Lys-182 in each of the three structures likely serves to stabilize the cationic form of the residue for participation in catalysis. CONCLUSIONS: The pro S hydrogen of NADPH (which is transferred as a hydride to the enzyme's naphthol substrates) is directed toward the carbonyl carbon of the inhibitors that mimic an intermediate along the reaction coordinate. Modeling tetrahydroxynaphthalene and trihydroxynaphthalene in the active site shows steric and electrostatic repulsion between the extra hydroxyl oxygen of the former substrate and the sulfur atom of Met-283 (the C-terminal residue), which accounts, in part, for the 4-fold greater substrate specificity for trihydroxynaphthalene over tetrahydroxynaphthalene.     

Metal-activated histidine carbon donor hydrogen bonds contribute to metalloprotein folding and function

Carbon donor hydrogen bonds are typically weak interactions that contribute less than 2 kcal/mol, and provide only modest stabilization in proteins. One exception is the class of hydrogen bonds donated by heterocyclic side chain carbons. Histidine is capable of particularly strong interactions through the Cε¹ and Cδ² carbons when the imidazole is protonated or bound to metal. Given the frequent occurrence of metal-bound histidines in metalloproteins, we characterized the energies of these interactions through DFT calculations on model compounds. Imidazole-water hydrogen bonding could vary from −11.0 to −17.0 kcal/mol, depending on the metal identity and oxidation state. A geometric search of metalloprotein structures in the PDB identified a number of candidate His C-H···O hydrogen bonds which may be important for folding or function. DFT calculations on model complexes of superoxide reductase show a carbon donor hydrogen bond positioning a water molecule above the active site.     

Crystal structures of two engineered thiol trypsins

We have determined the three-dimensional structures of engineered rat trypsins which mimic the active sites of two classes of cysteine proteases. The catalytic serine was replaced with cysteine (S195C) to test the ability of sulfur to function as a nucleophile in a serine protease environment. This variant mimics the cysteine trypsin class of thiol proteases. An additional mutation of the active site aspartate to an asparagine (D102N) created the catalytic triad of the papain-type cysteine proteases. Rat trypsins S195C and D102N,S195C were solved to 2.5 and 2.0 A, respectively. The refined structures were analyzed to determine the structural basis for the 10(6)-fold loss of activity of trypsin S195C and the 10(8)-fold loss of activity of trypsin D102N,S195C, relative to rat trypsin. The active site thiols were found in a reduced state in contrast to the oxidized thiols found in previous thiol protease structures. These are the first reported structures of serine proteases with the catalytic centers of sulfhydryl proteases. Structure analysis revealed only subtle global changes in enzyme conformation. The substrate binding pocket is unaltered, and active site amino acid 102 forms hydrogen bonds to H57 and S214 as well as to the backbone amides of A56 and H57. In trypsin S195C, D102 is a hydrogen-bond acceptor for H57 which allows the other imidazole nitrogen to function as a base during catalysis. In trypsin D102N,S195C, the asparagine at position 102 is a hydrogen-bond donor to H57 which places a proton on the imidazole nitrogen proximal to the nucleophile. This tautomer of H57 is unable to act as a base in catalysis.(ABSTRACT TRUNCATED AT 250 WORDS)     

Influence of a lipid interface on protein dynamics in a fungal lipase.

Lipases catalyze lipolytic reactions and for optimal activity they require a lipid interface. To study the effect of a lipid aggregate on the behavior of the enzyme at the interfacial plane and how the aggregate influences an attached substrate or product molecule in time and space, we have performed molecular dynamics simulations. The simulations were performed over 1 to 2 ns using explicit SPC water. The interaction energies between protein and lipid are mainly due to van der Waals contributions reflecting the hydrophobic nature of the lipid molecules. Estimations of the protonation state of titratable residues indicated that the negative charge on the fatty acid is stabilized by interactions with the titratable residues Tyr-28, His-143, and His-257. In the presence of a lipid patch, the active site lid opens wider than observed in the corresponding simulations in an aqueous environment. In that lid conformation, the hydrophobic residues Ile-85, Ile-89, and Leu-92 are embedded in the lipid patch. The behavior of the substrate or product molecule is sensitive to the environment. Entering and leaving of substrate molecules could be observed in presence of the lipid patch, whereas the product forms strong hydrogen bonds with Ser-82, Ser-144, and Trp-88, suggesting that the formation of hydrogen bonds may be an important contribution to the mechanism by which product inhibition might take place.