A General Acid−Base Mechanism for the Stabilization of a Tetrahedral Adduct in a Serine−Carboxyl Peptidase: A Computational Study

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).

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

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.

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 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.

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.

Degenerate intermolecular and intramolecular proton-transfer reactions: electronic structure of the transition states

Density functional theory (DFT) calculations are performed on a series of double and single proton-transfer reactions to study the variation in polarizations in complexes during the dynamics of proton transfer from one isoenergetic, hydrogen-bonded ground-state structure to the other. The isotropic average polarizability (alpha(av)) shows an interesting single-humped profile with a maxima coinciding with the transition state of the reaction. Similar profiles are also computed at Nd:YaG frequencies. The origin of the maximal polarizability at the transition state is traced to maximal charge separation and large D (donor)-A (acceptor) distances. Maximal polarizability for the transition state suggests an interesting, novel, and less memory extensive computational tool to locate the transition state for hydrogen-transfer reactions in hydrogen-bonded complexes.

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.

The importance of dynamics in substrate-assisted catalysis and specificity

The QM/MM MD and free energy simulations show that the dynamics involving a His residue at the P1 site of the substrate may play an important role in substrate-assisted catalysis and specificity for a serine-carboxyl peptidase.

Evolution of prokaryotic subtilases: genome-wide analysis reveals novel subfamilies with different catalytic residues

Subtilisin-like serine proteases (subtilases) are a very diverse family of serine proteases with low sequence homology, often limited to regions surrounding the three catalytic residues. Starting with different Hidden Markov Models (HMM), based on sequence alignments around the catalytic residues of the S8 family (subtilisins) and S53 family (sedolisins), we iteratively searched all ORFs in the complete genomes of 313 eubacteria and archaea. In 164 genomes we identified a total of 567 ORFs with one or more of the conserved regions with a catalytic residue. The large majority of these contained all three regions around the "classical" catalytic residues of the S8 family (Asp-His-Ser), while 63 proteins were identified as S53 (sedolisin) family members (Glu-Asp-Ser). More than 30 proteins were found to belong to two novel subsets with other evolutionary variations in catalytic residues, and new HMMs were generated to search for them. In one subset the catalytic Asp is replaced by an equivalent Glu (i.e. Glu-His-Ser family). The other subset resembles sedolisins, but the conserved catalytic Asp is not located on the same helix as the nucleophile Glu, but rather on a beta-sheet strand in a topologically similar position, as suggested by homology modeling. The Prokaryotic Subtilase Database (www.cmbi.ru.nl/subtilases) provides access to all information on the identified subtilases, the conserved sequence regions, the proposed family subdivision, and the appropriate HMMs to search for them. Over 100 proteins were predicted to be subtilases for the first time by our improved searching methods, thereby improving genome annotation.

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

Quantum mechanical/molecular mechanical molecular dynamics and free energy simulations are performed to study the acylation reaction catalyzed by kumamolisin-As, a serine-carboxyl peptidase, and to elucidate the catalytic mechanism and the origin of substrate specificity. It is demonstrated that the nucleophilic attack by the serine residue on the substrate may not be the rate-limiting step for the acylation of the GPH*FF substrate. The present study also confirms the earlier suggestions that Asp164 acts as a general acid during the catalysis and that the electrostatic oxyanion hole interactions may not be sufficient to lead a stable tetrahedral intermediate along the reaction pathway. Moreover, Asp164 is found to act as a general base during the formation of the acyl-enzyme from the tetrahedral intermediate. The role of dynamic substrate assisted catalysis (DSAC) involving His at the P1 site of the substrate is examined for the acylation reaction. It is demonstrated that the bond-breaking and -making events at each stage of the reaction trigger a change of the position for the His side chain and lead to the formation of the alternative hydrogen bonds. The back and forth movements of the His side chain between the C=O group of Pro at P2 and Odelta2 of Asp164 in a ping-pong-like mechanism and the formation of the alternative hydrogen bonds effectively lower the free energy barriers for both the nucleophilic attack and the acyl-enzyme formation and may therefore contribute to the relatively high activity of kumamolisin-As toward the substrates with His at the P1 site.

Processing, catalytic activity and crystal structures of kumamolisin-As with an engineered active site

Kumamolisin-As is an acid collagenase with a subtilisin-like fold. Its active site contains a unique catalytic triad, Ser278-Glu78-Asp82, and a putative transition-state stabilizing residue, Asp164. In this study, the mutants D164N and E78H/D164N were engineered in order to replace parts of the catalytic machinery of kumamolisin-As with the residues found in the equivalent positions in subtilisin. Unlike the wild-type and D164N proenzymes, which undergo instantaneous processing to produce their 37-kDa mature forms, the expressed E78H/D164N proenzyme exists as an equilibrated mixture of the nicked and intact forms of the precursor. X-ray crystallographic structures of the mature forms of the two mutants showed that, in each of them, the catalytic Ser278 makes direct hydrogen bonds with the side chain of Asn164. In addition, His78 of the double mutant is distant from Ser278 and Asp82, and the catalytic triad no longer exists. Consistent with these structural alterations around the active site, these mutants showed only low catalytic activity (relative k(cat) at pH 4.0 1.3% for D164N and 0.0001% for E78H/D164N). pH-dependent kinetic studies showed that the single D164N substitution did not significantly alter the logk(cat) vs. pH and log(k(cat)/Km) vs. pH profiles of the enzyme. In contrast, the double mutation resulted in a dramatic switch of the logk(cat) vs. pH profile to one that was consistent with catalysis by means of the Ser278-His78 dyad and Asn164, which may also account for the observed ligation/cleavage equilibrium of the precursor of E78H/D164N. These results corroborate the mechanistic importance of the glutamate-mediated catalytic triad and oxyanion-stabilizing aspartic acid residue for low-pH peptidase activity of the enzyme.

Tripeptidyl-peptidase I in health and disease

The lysosomal lumen contains numerous acidic hydrolases involved in the degradation of carbohydrates, lipids, proteins, and nucleic acids, which are basic cell components that turn over continuously within the cell and/or are ingested from outside of the cell. Deficiency in almost any of these hydrolases causes accumulation of the undigested material in secondary lysosomes, which manifests itself as a form of lysosomal storage disorder (LSD). Mutations in tripeptidyl-peptidase I (TPP I) underlie the classic late-infantile form of neuronal ceroid lipofuscinoses (CLN2), the most common neurodegenerative disorders of childhood. TPP I is an aminopeptidase with minor endopeptidase activity and Ser475 serving as an active-site nucleophile. The enzyme is synthesized as a highly glycosylated precursor transported by mannose-6-phosphate receptors to lysosomes, where it undergoes proteolytic maturation. This review summarizes recent progress in understanding of TPP I biology and molecular pathology of the CLN2 disease process, including distribution of the enzyme, its biosynthesis, glycosylation, transport and activation, as well as catalytic mechanisms and their potential implications for pathogenesis and treatment of the underlying disease. Promising data from gene and stem cell therapy in laboratory animals raise hope that CLN2 will be the first neurodegenerative LSD for which causative treatment will become available for humans.