Insight Derived from Molecular Dynamics Simulation into Substrate-Induced Changes in Protein Motions of Proteinase K

Abstract

Because of the significant industrial, agricultural and biotechnological importance of serine protease proteinase K, it has been extensively investigated using experimental approaches such as X-ray crystallography, site-directed mutagenesis and kinetic measurement. However, detailed aspects of enzymatic mechanism such as substrate binding, release and relevant regulation remain unstudied. Molecular dynamics (MD) simulations of the proteinase K alone and in complex with the peptide substrate AAPA were performed to investigate the effect of substrate binding on the dynamics/molecular motions of proteinase K. The results indicate that during simulations the substrate-complexed proteinase K adopt a more compact and stable conformation than the substrate-free form. Further essential dynamics (ED) analysis reveals that the major internal motions are confined within a subspace of very small dimension. Upon substrate binding, the overall flexibility of the protease is reduced; and the noticeable displacements are observed not only in substrate-binding regions but also in regions opposite the substrate-binding groove/pockets. The dynamic pockets caused by the large concerted motions are proposed to be linked to the substrate recognition, binding, orientation and product release; and the significant displacements in regions opposite the binding groove/pockets are considered to play a role in modulating the dynamics of enzyme-substrate interaction. Our simulation results complement the biochemical and structural studies, highlighting the dynamic mechanism of the functional properties of proteinase K.     

MD simulations reveal alternate conformations of the oxyanion hole in the Zika virus NS2B/NS3 protease

Recent crystallography studies have shown that the binding site oxyanion hole plays an important role in inhibitor binding, but can exist in two conformations (active/inactive). We have undertaken molecular dynamics (MD) calculations to better understand oxyanion hole dynamics and thermodynamics. We find that the Zika virus (ZIKV) NS2B/NS3 protease maintains a stable closed conformation over multiple 100-ns conventional MD simulations in both the presence and absence of inhibitors. The S1, S2, and S3 pockets are stable as well. However, in two of eight simulations, the A132-G133 peptide bond in the binding pocket of S1' spontaneously flips to form a 3₁₀-helix that corresponds to the inactive conformation of the oxyanion hole, and then maintains this conformation until the end of the 100-ns conventional MD simulations without inversion of the flip. This conformational change affects the S1' pocket in ZIKV NS2B/NS3 protease active site, which is important for small molecule binding. The simulation results provide evidence at the atomic level that the inactive conformation of the oxyanion hole is more favored energetically when no specific interactions are formed between substrate/inhibitor and oxyanion hole residues. Interestingly, however, transition between the active and inactive conformation of the oxyanion hole can be observed by boosting the valley potential in accelerated MD simulations. This supports a proposed induced-fit mechanism of ZIKV NS2B/NS3 protease from computational methods and provides useful direction to enhance inhibitor binding predictions in structure-based drug design.     

Molecular basis for substrate selectivity of a mono- and diacylglycerol lipase from Malassezia globosa

The lipase from Malassezia globosa (SMG1) was identified to be strictly specific for mono- and diacylglycerol but not triacylglycerol. The crystal structures of SMG1 were solved in the closed conformation, but they failed to provide direct evidence of factors responsible for this unique selectivity. To address this problem, we constructed a structure in the open, active conformation and modeled a diacylglycerol analogue into the active site. Molecular dynamics simulations were performed on this enzyme-analogue complex to relax steric clashes. This bound diacylglycerol analogue unambiguously identified the position of two pockets which accommodated two alkyl chains of substrate. The structure of SMG1-analogue complex revealed that Leu103 and Phe278 divided the catalytic pocket into two separated moieties, an exposed groove and a narrow tunnel. Analysis of the binding model suggested that the unique selectivity of this lipase mainly resulted from the shape and size of this narrow tunnel, in which there was no space for the settlement of the third chain of triacylglycerol. These results expand our understanding on the mechanism underlying substrate selectivity of enzyme, and could pave the way for site-directed mutagenesis experiments to improve the enzyme for application.     

Regioselectivity of CYP2B6: homology modeling,molecular dynamics simulation,docking

Human cytochrome P450 (CYP) 2B6 activates the anticancer prodrug cyclophosphamide (CPA) by 4-hydroxylation. In contrast, the same enzyme catalyzes N-deethylation of a structural isomer, the prodrug ifosfamide (IFA), thus causing severe adverse drug effects. To model the molecular interactions leading to a switch in regioselectivity, the structure of CYP2B6 was modeled based on the structure of rabbit CYP2C5. We modeled the missing 22-residue loop in CYP2C5 between helices F and G (the F-G loop), which is not resolved in the X-ray structure, by molecular dynamics (MD) simulations using a simulated annealing protocol. The modeled conformation of the loop was validated by unconstrained MD simulations of the complete enzymes (CYP2C5 and CYP2B6) in water for 70 and 120 ps, respectively. The simulations were stable and led to a backbone r.m.s. deviation of 1.7 A between the two CYPs.The shape of the substrate binding site of CYP2B6 was further analyzed. It consists of three well-defined hydrophobic binding pockets adjacent to the catalytic heme. Size, shape and hydrophobicity of these pockets were compared to the shapes of the two structurally isomeric substrates. In their preferred orientation in the binding site, both substrates fill all three binding pockets without repulsive interactions. The distance to the heme iron is short enough for 4-hydroxylation and N-deethylation to occur for CPA and IFA, respectively. However, if the substrates are docked in the non-preferred orientation (such that 4-hydroxylation and N-deethylation would occur for IFA and CPA, respectively), one pocket is left empty, and clashes were observed between the substrates and the enzyme.     

Mutation of a critical arginine in microsomal prostaglandin E synthase-1 shifts the isomerase activity to a reductase activity that converts prostaglandin H2 into prostaglandin F2alpha

Microsomal prostaglandin E synthase type 1 (mPGES-1) converts prostaglandin endoperoxides, generated from arachidonic acid by cyclooxygenases, into prostaglandin E2. This enzyme belongs to the membrane-associated proteins in eicosanoid and glutathione metabolism (MAPEG) family of integral membrane proteins, and because of its link to inflammatory conditions and preferential coupling to cyclooxygenase 2, it has received considerable attention as a drug target. Based on the high resolution crystal structure of human leukotriene C4 synthase, a model of mPGES-1 has been constructed in which the tripeptide co-substrate glutathione is bound in a horseshoe-shaped conformation with its thiol group positioned in close proximity to Arg-126. Mutation of Arg-126 into an Ala or Gln strongly reduces the enzyme's prostaglandin E synthase activity (85-95%), whereas mutation of a neighboring Arg-122 does not have any significant effect. Interestingly, R126A and R126Q mPGES-1 exhibit a novel, glutathione-dependent, reductase activity, which allows conversion of prostaglandin H2 into prostaglandin F2alpha. Our data show that Arg-126 is a catalytic residue in mPGES-1 and suggest that MAPEG enzymes share significant structural components of their active sites.     

Molecular dynamics simulations of arachidonic acid complexes with COX-1 and COX-2: insights into equilibrium behavior

The cyclooxygenase (COX) enzymes are responsible for the committed step in prostaglandin biosynthesis, the generation of prostaglandin H(2). As a result, these enzymes are pharmacologically important targets for nonsteroidal antiinflammatory drugs, such as aspirin and newer COX-2 selective inhibitors. The cyclooxygenases are functional homodimers, and each subunit contains both a cyclooxygenase and a peroxidase active site. These enzymes are quite interesting mechanistically, as the conversion of arachidonic acid to prostaglandin H(2) requires two oxygenation and two cyclization reactions, resulting in the formation of five new chiral centers with nearly absolute regio- and stereochemical fidelity. We have used molecular dynamics (MD) simulations to investigate the equilibrium behavior of both COX-1 and COX-2 enzyme isoforms with bound arachidonate. These simulations were compared with reference simulations of arachidonate in solution to explore the effect of enzyme on substrate conformation and positioning in the active site. The simulations suggest that the substrate has greater conformational freedom in the COX-2 active site, consistent with the larger COX-2 active site volume observed in X-ray crystal structures. The simulations reveal different conformational behavior for arachidonate in each subunit over the course of extended equilibrium MD simulations. The simulations also provide detailed information for several protein channels that might be important for oxygen and water transport to or from active sites or for intermediate trafficking between the cyclooxygenase and peroxidase active sites. The detailed comparisons for COX-1 versus COX-2 active site structural fluctuations may also provide useful information for design of new isozyme-selective inhibitors.     

Structural and functional characterization of human microsomal prostaglandin E synthase-1 by computational modeling and site-directed mutagenesis

Microsomal prostaglandin (PG) E synthase-1 (mPGES-1) has recently been recognized as a novel, promising drug target for inflammation-related diseases. Functional and pathological studies on this enzyme further stimulate to understand its structure and the structure-function relationships. Using an approach of the combined structure prediction, molecular docking, site-directed mutagenesis, and enzymatic activity assay, we have developed the first three-dimensional (3D) model of the substrate-binding domain (SBD) of mPGES-1 and its binding with substrates prostaglandin H2 (PGH2) and glutathione (GSH). In light of the 3D model, key amino acid residues have been identified for the substrate binding and the obtained experimental activity data have confirmed the computationally determined substrate-enzyme binding mode. Both the computational and experimental results show that Y130 plays a vital role in the binding with PGH2 and, probably, in the catalytic reaction process. R110 and T114 interact intensively with the carboxyl tail of PGH2, whereas Q36 and Q134 only enhance the PGH2-binding affinity. The modeled binding structure indicates that substrate PGH2 interacts with GSH through hydrogen binding between the peroxy group of PGH2 and the -SH group of GSH. The -SH group of GSH is expected to attack the peroxy group of PGH2, initializing the catalytic reaction transforming PGH2 to prostaglandin E2 (PGE2). The overall agreement between the calculated and experimental results demonstrates that the predicted 3D model could be valuable in future rational design of potent inhibitors of mPGES-1 as the next-generation inflammation-related therapeutic.     

A short-chain carbonyl reductase mutant is an efficient catalyst in the production of (R)-1,3-butanediol

R-1,3-butanediol (R-1,3-BDO) is an important chiral intermediate of penem and carbapenem synthesis. Among the different synthesis methods to obtain pure enantiomer R-1,3-BDO, oxidation–reduction cascades catalysed by enzymes are promising strategies for its production. Dehydrogenases have been used for the reduction step, but the enantio-selectivity is not high enough for further organic synthesis efforts. Here, a short-chain carbonyl reductase (LnRCR) was evaluated for the reduction step and developed via protein engineering. After docking result analysis with the substrate 4-hydroxy-2-butanone (4H2B), residues were selected for virtual mutagenesis, their substrate-binding energies were compared, and four sites were selected for saturation mutagenesis. High-throughput screening helped identify a Ser154Lys mutant which increased the catalytic efficiency by 115% compared to the parent enzyme. Computer-aided simulations indicated that after single residue replacement, movements in two flexible areas (VTDPAF and SVGFANK) facilitated the volumetric compression of the 4H2B-binding pocket. The number of hydrogen bonds between the stabilized 4H2B-binding pocket of the mutant enzyme and substrate was higher (from four to six) than the wild-type enzyme, while the substrate-binding energy was decreased (from −17.0 kJ/mol to −29.1 kJ/mol). Consequently, the catalytic efficiency increased by approximately 115% and enantio-selectivity increased from 95% to 99%. Our findings indicate that compact and stable substrate-binding pockets are critical for enzyme catalysis. Lastly, the utilization of a microbe expressing the Ser154Lys mutant enzyme was proven to be a robust process to conduct the oxidation–reduction cascade at larger scales.     

Phage display and crystallographic analysis reveals potential substrate/binding site interactions in the protein secretion chaperone CsaA from Agrobacterium tumefaciens

The protein CsaA has been proposed to function as a protein secretion chaperone in bacteria that lack the Sec-dependent protein-targeting chaperone SecB. CsaA is a homodimer with two putative substrate-binding pockets, one in each monomer. To test the hypothesis that these cavities are indeed substrate-binding sites able to interact with other polypeptide chains, we selected a peptide that bound to CsaA from a random peptide library displayed on phage. Presented here is the structure of CsaA from Agrobacterium tumefaciens (AtCsaA) solved in the presence and absence of the selected peptide. To promote co-crystallization, the sequence for this peptide was genetically fused to the amino-terminus of AtCsaA. The resulting 1.65 Å resolution crystal structure reveals that the tethered peptide from one AtCsaA molecule binds to the proposed substrate-binding pocket of a symmetry-related molecule possibly mimicking the interaction between a pre-protein substrate and CsaA. The structure shows that the peptide lies in an extended conformation with alanine, proline and glutamine side chains pointing into the binding pocket. The peptide interacts with the atoms of the AtCsaA-binding pocket via seven direct hydrogen bonds. The side chain of a conserved pocket residue, Arg76, has an “up” conformation when the CsaA-binding site is empty and a “down” conformation when the CsaA-binding site is occupied, suggesting that this residue may function to stabilize the peptide in the binding cavity. The presented aggregation assays, phage-display analysis and structural analysis are consistent with AtCsaA being a general chaperone. The properties of the proposed CsaA-binding pocket/peptide interactions are compared to those from other structurally characterized molecular chaperones.     

A molecular dynamics investigation of mono and dimeric states of the outer membrane enzyme OMPLA

OMPLA is a phospholipase found in the outer membranes of many Gram-negative bacteria. Enzyme activation requires calcium-induced dimerisation plus bilayer perturbation. As the conformation of OMPLA in the different crystal forms (monomer versus dimer; with/without bound Ca(2+)) is remarkably similar we have used multi-nanosecond molecular dynamics (MD) simulations to probe possible differences in conformational dynamics that may be related to enzyme activation. Simulations of calcium-free monomeric OMPLA, of the Ca(2+)-bound dimer, and of the Ca(2+)-bound dimer with a substrate analogue covalently linked to the active site serine have been performed, all with the protein embedded in a phospholipid (POPC) bilayer. All simulations were stable, but differences in the dynamic behaviour of the protein between the various states were observed. In particular, the stability of the active site and the hydrophobic substrate-binding cleft varied. Dimeric OMPLA is less flexible than monomeric OMPLA, especially around the active site. In the absence of bound substrate analogue, the hydrophobic substrate-binding cleft of dimeric OMPLA collapses. A model is proposed whereby the increased stability of the active site in dimeric OMPLA is a consequence of the local ordering of water around the nearby calcium ion. The observed collapse of the substrate-binding cleft may explain the experimentally observed occurrence of multiple dimer conformations of OMPLA, one of which is fully active while the other shows significantly reduced activity.     

HTRF-based assay for microsomal prostaglandin E2 synthase-1 activity

Microsomal prostaglandin E(2) synthase-1 (mPGES-1) catalyzes the formation of prostaglandin E(2) (PGE(2)) from the endoperoxide prostaglandin H( 2) (PGH(2)). Expression of this enzyme is induced during the inflammatory response, and mouse knockout experiments suggest it may be an attractive target for antiarthritic therapies. Assaying the activity of this enzyme in vitro is challenging because of the unstable nature of the PGH( 2) substrate. Here, the authors present an mPGES-1 activity assay suitable for characterization of enzyme preparations and for determining the potency of inhibitor compounds. This plate-based competition assay uses homogenous time-resolved fluorescence to measure PGE(2) produced by the enzyme. The assay is insensitive to DMSO concentration up to 10% and does not require extensive washes after the initial enzyme reaction is concluded, making it a simple and convenient way to assess mPGES-1 inhibition.     

Crystal structure of the terminal oxygenase component of biphenyl dioxygenase derived from Rhodococcus sp. strain RHA1

Biphenyl dioxygenase is the enzyme that catalyzes the stereospecific dioxygenation of the aromatic ring. This enzyme has attracted the attention of researchers due to its ability to oxidize polychlorinated biphenyls, which is one of the serious environmental contaminants. We determined the crystal structure of the terminal oxygenase component of the biphenyl dioxygenase (BphA1A2) derived from Rhodococcus strain sp. RHA1 in substrate-free and complex forms. These crystal structures revealed that the substrate-binding pocket makes significant conformational changes upon substrate binding to accommodate the substrate into the pocket. Our analysis of the crystal structures suggested that the residues in the substrate-binding pocket can be classified into three groups, which, respectively, seem to be responsible for the catalytic reaction, the orientation/conformation of the substrate, and the conformational changes of the substrate-binding pocket. The cooperative actions of residues in the three groups seem to determine the substrate specificity of the enzyme.     

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.     

Cyclooxygenase-2, Not Microsomal Prostaglandin E Synthase-1, Is the Mechanism for Interleukin-1β-induced Prostaglandin E2 Production and Inhibition of Insulin Secretion in Pancreatic Islets

Arachidonic acid is converted to prostaglandin E(2) (PGE(2)) by a sequential enzymatic reaction performed by two isoenzyme groups, cyclooxygenases (COX-1 and COX-2) and terminal prostaglandin E synthases (cPGES, mPGES-1, and mPGES-2). mPGES-1 is widely considered to be the final enzyme regulating COX-2-dependent PGE(2) synthesis. These generalizations have been based in most part on experiments utilizing gene expression analyses of cell lines and tumor tissue. To assess the relevance of these generalizations to a native mammalian tissue, we used isolated human and rodent pancreatic islets to examine interleukin (IL)-1β-induced PGE(2) production, because PGE(2) has been shown to mediate IL-1β inhibition of islet function. Rat islets constitutively expressed mRNAs of COX-1, COX-2, cPGES, and mPGES-1. As expected, IL-1β increased mRNA levels for COX-2 and mPGES-1, but not for COX-1 or cPGES. Basal protein levels of COX-1, cPGES, and mPGES-2 were readily detected in whole cell extracts but were not regulated by IL-1β. IL-1β increased protein levels of COX-2, but unexpectedly mPGES-1 protein levels were low and unaffected. In microsomal extracts, mPGES-1 protein was barely detectable in rat islets but clearly present in human islets; however, in neither case did IL-1β increase mPGES-1 protein levels. To further assess the importance of mPGES-1 to IL-1β regulation of an islet physiologic response, glucose-stimulated insulin secretion was examined in isolated islets of WT and mPGES-1-deficient mice. IL-1β inhibited glucose-stimulated insulin secretion equally in both WT and mPGES-1(-/-) islets, indicating that COX-2, not mPGES-1, mediates IL-1β-induced PGE(2) production and subsequent inhibition of insulin secretion.     

Biarylimidazoles as inhibitors of microsomal prostaglandin E2 synthase-1

Microsomal prostaglandin E(2) synthase (mPGES-1) represents a potential target for novel analgesic and anti-inflammatory agents. High-throughput screening identified several leads of mPGES-1 inhibitors which were further optimized for potency and selectivity. A series of inhibitors bearing a biaryl imidazole scaffold exhibits excellent inhibition of PGE(2) production in enzymatic and cell-based assays. The synthesis of these molecules and their activities will be discussed.     

Molecular dynamics simulations of the Apo-, Holo-, and acyl-forms of Escherichia coli acyl carrier protein

Acyl carrier protein (ACP) is an essential co-factor protein in fatty acid biosynthesis that shuttles covalently bound fatty acyl intermediates in its hydrophobic pocket to various enzyme partners. To characterize acyl chain-ACP interactions and their influence on enzyme interactions, we performed 19 molecular dynamics (MD) simulations of Escherichia coli apo-, holo-, and acyl-ACPs. The simulations were started with the acyl chain in either a solvent-exposed or a buried conformation. All four short-chain (< or = C10) and one long-chain (C16) unbiased acyl-ACP MD simulation show the transition of the solvent-exposed acyl chain into the hydrophobic pocket of ACP, revealing its pathway of acyl chain binding. Although the acyl chain resides inside the pocket, Thr-39 and Glu-60 at the entrance stabilize the phosphopantetheine linker through hydrogen bonding. Comparisons of the different ACP forms indicate that the loop region between helices II and III and the prosthetic linker may aid in substrate recognition by enzymes of fatty acid synthase systems. The MD simulations consistently show that the hydrophobic binding pocket of ACP is best suited to accommodate an octanoyl group and is capable of adjusting in size to accommodate chain lengths as long as decanoic acid. The simulations also reveal a second, novel binding mode of the acyl chains inside the hydrophobic binding pocket directed toward helix I. This study provides a detailed dynamic picture of acyl-ACPs that is in excellent agreement with available experimental data and, thereby, provides a new understanding of enzyme-ACP interactions.     

New structural motifs on the chymotrypsin fold and their potential roles in complement factor B

Factor B and C2 are two central enzymes for complement activation. They are multidomain serine proteases and require cofactor binding for full expression of proteolytic activities. We present a 2.1 A crystal structure of the serine protease domain of factor B. It shows a number of structural motifs novel to the chymotrypsin fold, which by sequence homology are probably present in C2 as well. These motifs distribute characteristically on the protein surface. Six loops surround the active site, four of which shape substrate-binding pockets. Three loops next to the oxyanion hole, which typically mediate zymogen activation, are much shorter or absent. Three insertions including the linker to the preceding domain bulge from the side opposite to the active site. The catalytic triad and non-specific substrate-binding site display active conformations, but the oxyanion hole displays a zymogen-like conformation. The bottom of the S1 pocket has a negative charge at residue 226 instead of the typical 189 position. These unique structural features may play different roles in domain-domain interaction, cofactor binding and substrate binding.     

Effects of protonation state of Asp181 and position of active site water molecules on the conformation of PTP1B

In protein tyrosine phosphatase 1B (PTP1B), the flexible WPD loop adopts a closed conformation (WPD_closed) in the active state of PTP1B, bringing the catalytic Asp181 close to the active site pocket, while WPD loop is in an open conformation (WPD_open) in the inactive state. Previous studies showed that Asp181 may be protonated at physiological pH, and ordered water molecules exist in the active site. In the current study, molecular dynamics simulations are employed at different Asp181 protonation states and initial positions of active site water molecules, and compared with the existing crystallographic data of PTP1B. In WPD_closed conformation, the active site is found to maintain its conformation only in the protonated state of Asp181 in both free and liganded states, while Asp181 is likely to be deprotonated in WPD_open conformation. When the active site water molecule network that is a part of the free WPD_closed crystal structure is disrupted, intermediate WPD loop conformations, similar to that in the PTPRR crystal structure, are sampled in the MD simulations. In liganded PTP1B, one active site water molecule is found to be important for facilitating the orientation of Cys215 and the phosphate ion, thus may play a role in the reaction. In conclusion, conformational stability of WPD loop, and possibly catalytic activity of PTP1B, is significantly affected by the protonation state of Asp181 and position of active site water molecules, showing that these aspects should be taken into consideration both in MD simulations and inhibitor design. © Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Conformational dynamics of threonine 195 and the S1 subsite in functional trypsin variants

Replacing the catalytic serine in trypsin with threonine (S195T variant) leads to a nearly complete loss of catalytic activity, which can be partially restored by eliminating the C42-C58 disulfide bond. The 0.69 μs of combined explicit solvent molecular dynamics (MD) simulations revealed continuous rearrangement of T195 with different conformational preferences between five trypsin variants tested. Among three conformational families observed for the T195 residue, one showed the T195 hydroxyl in a conformation analogous to that of the serine residue in wild-type trypsin, positioning the hydroxyl oxygen atom for attack on the carbonyl carbon of the peptide substrate. MD simulations demonstrated that this conformation was more populated for the C42A/C58V/S195T and C42A/C58A/S195T triple variants than for the catalytically inactive S195T variant and correlated with restored enzymatic activities for triple variants. In addition, observation of the increased motion of the S214-G219 segment in the S195T substituted variants suggested an existence of open and closed conformations for the substrate binding pocket. The closed conformation precludes access to the S1 binding site and could further reduce enzymatic activities for triple variants. Double variants with intact serine residues (C42A/C58A/S195 and C42A/C58V/S195) also showed interchange between closed and open conformations for the S214-G219 segment, but to a lesser extent than the triple variants. The increased conformational flexibility of the S1 subsite, which was not observed for the wild-type, correlated with reduced enzymatic activities and suggested a possible mode of substrate regulation for the trypsin variants tested.