The QM/MM molecular dynamics and free energy simulations of the acylation reaction catalyzed by the serine-carboxyl peptidase kumamolisin-As

Lowering pH optimum of activity of SshEstI, a slightly alkaliphilic archaeal esterase of the hormone-sensitive lipase family

SshEstI, a carboxylesterase from the thermoacidophilic archaeon Saccharolobus shibatae, is a member of the hormone-sensitive lipase family that displays slightly alkaliphilic activity with an optimum activity at pH 8.0. In this study, three distinct strategies were explored to confer acidophilic properties to SshEstI. The first strategy involved engineering the oxyanion hole by replacing Gly81 with serine or aspartic acid. The G81S mutant showed optimum activity at pH 7.0, whereas the aspartic acid mutant (G81D) rendered the enzyme slightly acidophilic with optimum activity observed at pH 6.0; however, kcat and kcat/Km values were reduced by these substitutions. The second strategy involved examining the effects of surfactant additives on the pH-activity profiles of SshEstI. The results showed that cetyltrimethylammonium bromide (CTAB) enhanced wild-type enzyme (WT) activity at acidic pH values. In the presence of 0.1 mM CTAB, G81S and G81D were acidophilic enzymes with optimum activity at pH 6.0 and 4.0, respectively, although their enzyme activities were low. The third strategy involved engineering the active site to resemble that of kumamolisin-As (kuma-As), an acidophilic peptidase of the sedolisin family. The catalytic triad of kuma-As was exchanged into SshEstI using site-directed mutagenesis. X-ray crystallographic analysis of the mutants (H274D and H274E) revealed that the potential hydrogen donor-acceptor distances around the active site of WT were fully maintained in these mutants. However, these mutants were inactive at pH 4-8.

Serine-Carboxyl Peptidases, Sedolisins: From Discovery to Evolution

Sedolisin is a proteolytic enzyme, listed in the peptidase database MEROPS as a founding member of clan SB, family S53. This enzyme, although active at low pH, was originally shown not to be inhibited by an aspartic peptidase specific inhibitor, S-PI (pepstatin Ac). In this Perspective, the S53 family is described from the moment of original identification to evolution. The representative enzymes of the family are sedolisin, kumamolisin, and TPP-1. They exhibit the following unique features. (1) The fold of the molecule is similar to that of subtilisin, but the catalytic residues consist of a triad, Ser/Glu/Asp, that is unlike the Ser/His/Asp triad of subtilisin. (2) The molecule is expressed as a pro-form composed of the amino-terminal prosegment and the active domain. Additionally, some members of this family have an additional, carboxy-terminal prosegment. (3) Their optimum pH for activity is in the acidic region, not in the neutral to alkaline region where subtilisin is active. (4) Their distribution in nature is very broad across the three kingdoms of life. (5) Some of these enzymes from fungi and bacteria are pathogens to plants. (6) Some of them have significant potential applications for industry. (7) The lack of a TPP-1 gene in human brain is the cause of incurable juvenile neuronal ceroid lipofuscinosis (Batten's disease).

Mechanisms of Proteolytic Enzymes and Their Inhibition in QM/MM Studies

Experimental evidence for enzymatic mechanisms is often scarce, and in many cases inadvertently biased by the employed methods. Thus, apparently contradictory model mechanisms can result in decade long discussions about the correct interpretation of data and the true theory behind it. However, often such opposing views turn out to be special cases of a more comprehensive and superior concept. Molecular dynamics (MD) and the more advanced molecular mechanical and quantum mechanical approach (QM/MM) provide a relatively consistent framework to treat enzymatic mechanisms, in particular, the activity of proteolytic enzymes. In line with this, computational chemistry based on experimental structures came up with studies on all major protease classes in recent years; examples of aspartic, metallo-, cysteine, serine, and threonine protease mechanisms are well founded on corresponding standards. In addition, experimental evidence from enzyme kinetics, structural research, and various other methods supports the described calculated mechanisms. One step beyond is the application of this information to the design of new and powerful inhibitors of disease-related enzymes, such as the HIV protease. In this overview, a few examples demonstrate the high potential of the QM/MM approach for sophisticated pharmaceutical compound design and supporting functions in the analysis of biomolecular structures.

Quantum Mechanics/Molecular Mechanics Simulations Identify the Ring-Opening Mechanism of Creatininase

Creatininase catalyzes the conversion of creatinine (a biosensor for kidney function) to creatine via a two-step mechanism: water addition followed by ring opening. Water addition is common to other known cyclic amidohydrolases, but the precise mechanism for ring opening is still under debate. The proton donor in this step is either His178 or a water molecule bound to one of the metal ions, and the roles of His178 and Glu122 are unclear. Here, the two possible reaction pathways have been fully examined by means of combined quantum mechanics/molecular mechanics simulations at the SCC-DFTB/CHARMM22 level of theory. The results indicate that His178 is the main catalytic residue for the whole reaction and explain its role as proton shuttle during the ring-opening step. In the first step, His178 provides electrostatic stabilization to the gem-diolate tetrahedral intermediate. In the second step, His178 abstracts the hydroxyl proton of the intermediate and delivers it to the cyclic amide nitrogen, leading to ring opening. The latter is the rate-limiting step with a free energy barrier of 18.5 kcal/mol, in agreement with the experiment. We find that Glu122 must be protonated during the enzyme reaction, so that it can form a stable hydrogen bond with its neighboring water molecule. Simulations of the E122Q mutant showed that this replacement disrupts the H-bond network formed by three conserved residues (Glu34, Ser78, and Glu122) and water, increasing the energy barrier. Our computational studies provide a comprehensive explanation for previous structural and kinetic observations, including why the H178A mutation causes a complete loss of activity but the E122Q mutation does not.

QM/MM Investigation of Substrate and Product Specificities of Suv4-20h2: How Does This Enzyme Generate Dimethylated H4K20 from Monomethylated Substrate?

Protein lysine methyltransferases (PKMTs) catalyze the methylation of lysine residues on histone proteins in the regulation of chromatin structure and gene expression. In contrast to many other PKMTs for which unmodified lysine is the methylation target, the enzymes in the Suv4-20 family are able to generate dimethylated product (H4K20me2) based exclusively on the monomethylated H4K20 substrate (H4K20me1). The origin of such substrate/product specificity is still not clear. Here, molecular dynamics (MD) and free energy (potential of mean force) simulations are undertaken using quantum mechanical/molecular mechanical (QM/MM) potentials to understand the substrate/product specificities of Suv4-20h2, a member of the Suv4-20 family. The free energy barriers for mono-, di-, and trimethylation in Suv4-20h2 obtained from the simulations are found to be well correlated with the specificities observed experimentally with the allowed dimethylation based on the H4K20me1 substrate and prohibited monomethylation and trimethylation based on H4K20 and H4K20me2, respectively. It is demonstrated that the reason for the relatively efficient dimethylation is an effective transition state (TS) stabilization through strengthening the CH···O interactions as well as the presence of a cation-π interaction at the transition state. The simulations also show that the failures of Suv4-20h2 to catalyze monomethylation and trimethylation are due, respectively, to a less effective TS stabilization and inability of the reactant complex containing H4K20me2 to adopt a reactive (near attack) configuration for methyl transfer. The results suggest that care must be exercised in the prediction of the substrate specificity based only on the existence of near attack configurations in substrate complexes.

Analysis of catalytic properties of tripeptidyl peptidase I (TTP-I), a serine carboxyl lysosomal protease, and its detection in tissue extracts using selective FRET peptide substrate

Tripeptidyl peptidase I (TPP-I), also named ceroid lipofuscinosis 2 protease (CLN2p), is a serine carboxyl lysosomal protease involved in neurodegenerative diseases, and has both tripeptidyl amino- and endo- peptidase activities under different pH conditions. We developed fluorescence resonance energy transfer (FRET) peptides using tryptophan (W) as the fluorophore to study TPP-I hydrolytic properties based on previous detailed substrate specificity study (Tian Y. et al., J. Biol. Chem. 2006, 281:6559-72). Tripeptidyl amino peptidase activity is enhanced by the presence of amino acids in the prime side and the peptide NH2-RWFFIQ-EDDnp is so far the best substrate described for TPP-I. The hydrolytic parameters of this peptide and its analogues indicated that the S4 subsite of TPP-I is occluded and there is an electrostatic interaction of the positively charged substrate N-terminus amino group and a negative locus in the region of the enzyme active site. KCl activated TPP-I in contrast to the inhibition by Ca(2+) and NaCl. Solvent kinetic isotope effects (SKIEs) show the importance of the free N-terminus amino group of the substrates, whose absence results in a more complex solvent-dependent enzyme: substrate interaction and catalytic process. Like pure TPP-I, rat spleen and kidney homogenates cleaved NH2-RWFFIQ-EDDnp only at F-F bond and is not inhibited by pepstatin, E-64, EDTA or PMSF. The selectivity of NH2-RWFFIQ-EDDnp to TPP-I was also demonstrated by the 400 times higher k(cat)/K(M) compared to generally used substrate, NH2-AAF-MCA and by its resistance to hydrolysis by cathepsin D that is present in high levels in kidneys.

QM/MM MD and Free Energy Simulation Study of Methyl Transfer Processes Catalyzed by PKMTs and PRMTs

Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here, we review the results of quantum mechanics/molecular mechanics molecular dynamics and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.

Diversity, Structures, and Collagen-Degrading Mechanisms of Bacterial Collagenolytic Proteases

Bacterial collagenolytic proteases are important because of their essential role in global collagen degradation and because of their virulence in some human bacterial infections. Bacterial collagenolytic proteases include some metalloproteases of the M9 family from Clostridium or Vibrio strains, some serine proteases distributed in the S1, S8, and S53 families, and members of the U32 family. In recent years, there has been remarkable progress in discovering new bacterial collagenolytic proteases and in investigating the collagen-degrading mechanisms of bacterial collagenolytic proteases. This review provides comprehensive insight into bacterial collagenolytic proteases, especially focusing on the structures and collagen-degrading mechanisms of representative bacterial collagenolytic proteases in each family. The roles of bacterial collagenolytic proteases in human diseases and global nitrogen cycling, together with the biotechnological and medical applications for these proteases, are also briefly discussed.

QM/MM MD and free energy simulation study of methyl transfer processes catalyzed by PKMTs and PRMTs

Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here we review the results of quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.

QM/MM MD and free energy simulations of the methylation reactions catalyzed by protein arginine methyltransferase PRMT3

Protein arginine N-methyltransferases (PRMTs) catalyze the transfer of methyl group(s) from S-adenosyl-l-methionine (AdoMet) to the guanidine group of arginine residue in abundant eukaryotic proteins. Two major types of PRMTs have been identified in mammalian cells. Type I PRMTs catalyze the formation of asymmetric ω-NG, NG-dimethylarginine (ADMA), while Type II PRMTs catalyze the formation of symmetric ω-NG, N′G-dimethylarginine (SDMA). The two different methylation products (ADMA or SDMA) of the substrate could lead to different biological consequences. Although PRMTs have been the subject of extensive experimental investigations, the origin of the product specificity remains unclear. In this study, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed to study the reaction mechanism for one of Type I PRMTs, PRMT3, and to gain insights into the energetic origin of its product specificity (ADMA). Our simulations have identified some important interactions and proton transfers involving the active site residues. These interactions and proton transfers seem to be responsible, at least in part, in making the Nη2 atom of the substrate arginine the target of the both 1st and 2nd methylations, leading to the asymmetric dimethylation product. The simulations also suggest that the methyl transfer and proton transfer appear to be somehow concerted processes and that Glu326 is likely to function as the general base during the catalysis.

Mechanism of C-terminal intein cleavage in protein splicing from QM/MM molecular dynamics simulations

Protein splicing is a post-translational process in which a biologically inactive protein is activated by the release of a segment denoted as an intein. The process involves four steps. In the third, the scission of the intein takes place after the cyclization of the last amino acid of the segment, an asparagine. Little is known about the chemical reaction necessary for this cyclization. Experiments demonstrate that two histidines (the penultimate amino acid of the intein, and a histidine located 10 amino acids upstream) are relevant in the cyclization of the asparagine. We have investigated the mechanism and determinants of reaction in the GyrA intein focusing on the requirements for asparagine activation for its cyclization. First, the influence that the protonation states of these two histidines have on the orientation of the asparagine side chain is investigated by means of molecular dynamics simulation. Molecular dynamics simulations using the CHARMM27 force field were carried out on the three possible protonation states for each of these two histidines. The results indicate that the only protonation state in which the conformation of the system is suitable for cyclization is when the penultimate histidine is fully protonated (positively charged), and the upstream histidine is in the His(ε) neutral tautomeric form. The free energy profile for the reaction in which the asparagine is activated by a proton transfer to the upstream histidine is presented, computed by hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics at the SCCDFTB/CHARMM27 level of theory. The calculated free energy barrier for the reaction is 19.0 kcal mol(-1). B3LYP/6-31+G(d) QM/MM single-point calculations give a qualitatively a similar energy profile, although with somewhat higher energy barriers, in good agreement with the value derived from experiment of 25 kcal mol(-1) at 60 °C. QM/MM molecular dynamics simulations of the reactant, activated reactant and intermediate states highlight the importance of the Arg181-Val182-Asp183 segment in catalysing the reaction. Overall, the results indicate that nucleophilic activation of the asparagine for its cyclization by the upstream histidine acting as the base is a plausible mechanism for the C-terminal cleavage in protein splicing.

New families of carboxyl peptidases: serine-carboxyl peptidases and glutamic peptidases

Peptidases or proteinases are now classified into seven families based on the nature of the catalytic residues [MEROPS-the peptidase database (http://merops.sanger.ac.uk/)]. They are aspartic- (first described in 1993), cysteine- (1993), serine- (1993) metallo- (1993), threonine- (1997), glutamic- (2004) and asparagine-peptidase (2010). By using an S-PI (pepstatin Ac) as a probe, a new subfamily of serine peptidase, serine-carboxyl peptidase (sedolisin) was discovered in 2001. In addition, the sixth family of peptidase, glutamic peptidase (eqolisin) was also discovered in 2004. The former peptidase is widely distributed in nature from archea to mammals, including humans. One of these enzymes is related to a human fatal hereditable disease, Batten disease. In contrast, the distribution of the latter peptidases is limited, with most of them found in human or plant pathogenic fungi. One such enzyme was isolated from a fungal infection in an HIV-infected patient. In this review, the background of the findings, and crystal structures, catalytic mechanisms, substrates specificities and distribution of the new peptidase families are described.

Clarification of the mechanism of acylation reaction and origin of substrate specificity of the serine-carboxyl peptidase sedolisin through QM/MM free energy simulations

Quantum mechanical/molecular mechanical (QM/MM) free energy simulations are applied for understanding the mechanism of the acylation reaction catalyzed by sedolisin, a representative serine-carboxyl peptidase, leading to the acyl-enzyme (AE) and first product from the enzyme-catalyzed reaction. One of the interesting questions to be addressed in this work is the origin of the substrate specificity of sedolisin that shows a relatively high activity on the substrates with Glu at P(1) site. It is shown that the bond making and breaking events of the acylation reaction involving a peptide substrate (LLE*FL) seem to be accompanied by local conformational changes, proton transfers as well as the formation of alternative hydrogen bonds. The results of the simulations indicate that the conformational change of Glu at P(1) site and its formation of a low barrier hydrogen bond with Asp-170 (along with the transient proton transfer) during the acylation reaction might play a role in the relatively high specificity for the substrate with Glu at P(1) site. The role of some key residues in the catalysis is confirmed through free energy simulations. Glu-80 is found to act as a general base to accept a proton from Ser-287 during the nucleophilic attack and then as a general acid to protonate the leaving group (N-H of P(1')-Phe) during the cleavage of the scissile peptide bond. Another acidic residue, Asp-170, acts as a general acid catalyst to protonate the carbonyl of P(1)-Glu during the formation of the tetrahedral intermediate and as a general base for the formation of the acyl-enzyme. The energetic results from the free energy simulations support the importance of proton transfer from Asp-170 to the carbonyl of P(1)-Glu in the stabilization of the tetrahedral intermediate and the formation of a low-barrier hydrogen bond between the carboxyl group of P(1)-Glu and Asp-170 in the lowering of the free energy barrier for the cleavage of the peptide bond. Detailed analyses of the proton transfers during acylation are also given.

Understanding the mechanism of deacylation reaction catalyzed by the serine carboxyl peptidase kumamolisin-As: insights from QM/MM free energy simulations

Quantum mechanical/molecular mechanical (QM/MM) molecular dynamics and free energy simulations are performed to study the process of the deacylation reaction catalyzed by kumamolisin-As, a serine-carboxyl peptidase, and to elucidate the catalytic mechanism. The results given here suggest that Asp-164 acts as a general acid/base catalyst not only for the acylation reaction but also for the deacylation reaction. It is shown that the electrostatic oxyanion hole interactions may be less effective in transition state stabilization for the kumamolisin-As catalyzed reaction compared to the general acid/base mechanism involving the proton transfer from or to Asp-164. The dynamic substrate-assisted catalysis (DSAC) involving His at the P1 site of the substrate is found to be less important for the deacylation reaction than for the acylation reaction in the kumamolisin-As catalyzed reaction. The proton transfer processes during the enzyme-catalyzed process are examined and their role in the catalysis is discussed.

Combined QM/MM mechanistic study of the acylation process in furin complexed with the H5N1 avian influenza virus hemagglutinin's cleavage site

Combined quantum mechanical/molecular mechanical (QM/MM) techniques have been applied to investigate the detailed reaction mechanism of the first step of the acylation process by furin in which the cleavage site of the highly pathogenic avian influenza virus subtype H5N1 (HPH5) acts as its substrate. The energy profile shows a simultaneous mechanism, known as a concerted reaction, of the two subprocesses: the proton transfer from Ser368 to His194 and the nucleophilic attack on the carbonyl carbon of the scissile peptide of the HPH5 cleavage site with a formation of tetrahedral intermediate (INT). The calculated energy barrier for this reaction is 16.2 kcal.mol(-1) at QM/MM B3LYP/6-31+G*//PM3-CHARMM22 level of theory. Once the reaction proceeds, the ordering of the electrostatic stabilization by protein environment is of the enzyme-substrate < transition state < INT complexes. Asp153 was found to play the most important role in the enzymatic reaction by providing the highest degree of intermediate complex stabilization. In addition, the negatively charged carbonyl oxygen of INT is well stabilized by the oxyanion hole constructed by Asn295's carboxamide and Ser368's backbone.

Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme

Despite numerous structural and biochemical investigations, the catalytic mechanism of hairpin ribozyme self-cleavage remains elusive. To gain insight into the coupling of active site dynamics with activity of this small catalytic RNA, we analyzed a total of approximately 300 ns of molecular dynamics (MD) simulations. Our simulations predict improved global stability for an in vitro selected "gain of function" mutation, which is validated by native gel electrophoretic mobility shift assay. We observe that active site nucleobases and water molecules stabilize a geometry favorable to catalysis through a dynamic hydrogen bonding network. Simulations in which A38 is unprotonated show its N1 move into close proximity of the active site 2'-OH, indicating that A38 may act as a general base during cleavage, a role that has generally been discounted due to the longer distances observed in crystal structures involving inactivating substrate analogs. By contrast, simulations in which N1 of A38 is protonated place N1 in close proximity to the 5'-oxygen leaving group, which supports the proposal that A38 serves as a general acid. In analogy to protein enzymes, we discuss a plausible mechanism in which A38 acts bifunctionally and shuttles a proton directly from the 2'-OH to the 5'-oxygen. Furthermore, our simulations suggest an important role for protonation of N1 of A38 in promoting a favorable geometry similar to that observed in transition-state analog crystal structures, and support previously proposed roles of A38, G8, and long residency water molecules in transition-state stabilization.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Crystal structure and autoactivation pathway of the precursor form of human tripeptidyl-peptidase 1, the enzyme deficient in late infantile ceroid lipofuscinosis

Late infantile neuronal ceroid lipofuscinosis is a fatal childhood neurological disorder caused by a deficiency in the lysosomal protease tripeptidyl-peptidase 1 (TPP1). TPP1 represents the only known mammalian member of the S53 family of serine proteases, a group characterized by a subtilisin-like fold, a Ser-Glu-Asp catalytic triad, and an acidic pH optimum. TPP1 is synthesized as an inactive proenzyme (pro-TPP1) that is proteolytically processed into the active enzyme after exposure to low pH in vitro or targeting to the lysosome in vivo. In this study, we describe an endoglycosidase H-deglycosylated form of TPP1 containing four Asn-linked N-acetylglucosamines that is indistinguishable from fully glycosylated TPP1 in terms of autocatalytic processing of the proform and enzymatic properties of the mature protease. The crystal structure of deglycosylated pro-TPP1 was determined at 1.85 angstroms resolution. A large 151-residue C-shaped prodomain makes extensive contacts as it wraps around the surface of the catalytic domain with the two domains connected by a 24-residue flexible linker that passes through the substrate-binding groove. The proenzyme structure reveals suboptimal catalytic triad geometry with its propiece linker partially blocking the substrate-binding site, which together serve to prevent premature activation of the protease. Finally, we have identified numerous processing intermediates and propose a structural model that explains the pathway for TPP1 activation in vitro. These data provide new insights into TPP1 function and represent a valuable resource for constructing improved TPP1 variants for treatment of late infantile neuronal ceroid lipofuscinosis.

The expanding diversity of serine hydrolases

Serine hydrolases use a hydroxyl of a serine, assisted by one or more other residues, to cleave peptide bonds. They belong to several different families whose general mechanism is well known. However, the subtle structural differences that have recently been observed across a variety of families shed light on their functional diversity, including variations in mechanism of action, differences in the modes of substrate binding, and substrate-assisted orientation of catalytic residues. Of particular interest are the Rhomboid family serine proteinases that are active within the plasma membrane, for which several new structures have been reported. Because these enzymes are involved in biological and pathological processes, many are becoming important targets of drug design.