Exploring the Catalytic Reaction of Cysteine Proteases

Covalent small-molecule inhibitors of SARS-CoV-2 Mpro: Insights into their design, classification, biological activity, and binding interactions

Since 2020, many compounds have been investigated for their potential use in the treatment of SARS-CoV-2 infection. Among these agents, a huge number of natural products and FDA-approved drugs have been evaluated as potential therapeutics for SARS-CoV-2 using virtual screening and docking studies. However, the identification of the molecular targets involved in viral replication led to the development of rationally designed anti-SARS-CoV-2 agents. Among these targets, the main protease (Mpro) is one of the key enzymes needed in the replication of the virus. The data gleaned from the crystal structures of SARS-CoV-2 Mpro complexes with small-molecule covalent inhibitors has been used in the design and discovery of many highly potent and broad-spectrum Mpro inhibitors. The current review focuses mainly on the covalent type of SARS-CoV-2 Mpro inhibitors. The design, chemistry, and classification of these inhibitors were also in focus. The biological activity of these inhibitors, including their inhibitory activities against Mpro, their antiviral activities, and the SAR studies, were discussed. The review also describes the potential mechanism of the interaction between these inhibitors and the catalytic Cys145 residue in Mpro. Moreover, the binding modes and key binding interactions of these covalent inhibitors were also illustrated. The covalent inhibitors discussed in this review were of diverse chemical nature and origin. Their antiviral activity was mediated mainly by the inhibition of SARS-CoV-2 Mpro, with IC50 values in the micromolar to the nanomolar range. Many of these inhibitors exhibited broad-spectrum inhibitory activity against the Mpro enzymes of other coronaviruses (SARS-CoV-1 and MERS-CoV). The dual inhibition of the Mpro and PLpro enzymes of SARS-CoV-2 could also provide higher therapeutic benefits than Mpro inhibition. Despite the approval of nirmatrelvir by the FDA, many mutations in the Mpro enzyme of SARS-CoV-2 have been reported. Although some of these mutations did not affect the potency of nirmatrelvir, there is an urgent need to develop a second generation of Mpro inhibitors. We hope that the data summarized in this review could help researchers in the design of a new potent generation of SARS-CoV-2 Mpro inhibitors.

Exploring the extraction, functional properties, and industrial applications of papain from Carica papaya

Papain a protease enzyme naturally present in the Carica papaya has gained significant interest across several industries due to its unique properties and versatility. The unique structure of papain imparts the functionality that assists in elucidating how papain enzyme works and making it beneficial for a variety of purposes. This review highlights recent advancements in papain extraction techniques to enhance production efficiency to meet market demand. The extraction of papain from the Carica papaya plant offers various advantages such as cost-effectiveness, biodegradability, safety, and the ability to withstand a wide range of pH and temperature conditions. Key findings reveal that non-conventional papain extraction techniques offer significant advantages in terms of efficiency, product quality, and environmental sustainability. Furthermore, papain treatment enhances the value of final products due to its anti-bacterial, anti-oxidant, and anti-obesity properties. The ability of papain to hydrolyze a wide range of proteins across various conditions makes it a suitable protease enzyme. While the study emphasizes the advantages of papain, the study also acknowledges limitations such as the continuous research and development to optimize extraction processes which will help unlock papain's potential and meet the growing demand. © 2024 Society of Chemical Industry.

Computational Insights into SARS-CoV-2 Main Protease Mutations and Nirmatrelvir Efficacy: The Effects of P132H and P132H-A173V

Nirmatrelvir, a pivotal component of the oral antiviral Paxlovid for COVID-19, targets the SARS-CoV-2 main protease (Mpro) as a covalent inhibitor. Here, we employed combined computational methods to explore how the prevalent Omicron variant mutation P132H, alone and in combination with A173V (P132H-A173V), affects nirmatrelvir's efficacy. Our findings suggest that P132H enhances the noncovalent binding affinity of Mpro for nirmatrelvir, whereas P132H-A173V diminishes it. Although both mutants catalyze the rate-limiting step more efficiently than the wild-type (WT) Mpro, P132H slows the overall rate of covalent bond formation, whereas P132H-A173V accelerates it. Comprehensive analysis of noncovalent and covalent contributions to the overall binding free energy of the covalent complex suggests that P132H likely enhances Mpro sensitivity to nirmatrelvir, while P132H-A173V may confer resistance. Per-residue decompositions of the binding and activation free energies pinpoint key residues that significantly affect the binding affinity and reaction rates, revealing how the mutations modulate these effects. The mutation-induced conformational perturbations alter drug-protein local contact intensities and the electrostatic preorganization of the protein, affecting noncovalent binding affinity and the stability of key reaction states, respectively. Our findings inform the mechanisms of nirmatrelvir resistance and sensitivity, facilitating improved drug design and the detection of resistant strains.

Bacterial extracellular vesicles: Vital contributors to physiology from bacteria to host

Bacterial extracellular vesicles (bEVs) represent spherical particles with diameters ranging from 20 to 400 nm filled with multiple parental bacteria-derived components, including proteins, nucleic acids, lipids, and other biomolecules. The production of bEVs facilitates bacteria interacting with their environment and exerting biological functions. It is increasingly evident that the bEVs play integral roles in both bacterial and host physiology, contributing to environmental adaptations to functioning as health promoters for their hosts. This review highlights the current state of knowledge on the composition, biogenesis, and diversity of bEVs and the mechanisms by which different bEVs elicit effects on bacterial physiology and host health. We posit that an in-depth exploration of the mechanistic aspects of bEVs activity is essential to elucidate their health-promoting effects on the host and may facilitate the translation of bEVs into applications as novel natural biological nanomaterials.

Charge Relay Without Proton Transfer: Coupling of Two Short Hydrogen Bonds via Imidazole in Models of Catalytic Triad of Serine Protease Active Site

A homologous series of 20 substituted alcohol-imidazole-acetate model complexes imitating the charge relay system in Ser-His-Asp catalytic triad of serine proteases is considered quantum-chemically. We show qualitatively that the geometries of alcohol-imidazole and imidazole-acetate short hydrogen bonds are strongly coupled via the central imidazole and such complexes are capable of effectively relaying the charge from acetate to alcohol moiety upon relatively small concerted proton displacements. We hypothesize an alternative catalytic mechanism of serine proteases that does not require two complete proton transfers or hydrogen bond breakage between Ser and His residues.

Extracellular proteolysis in cancer: Proteases, substrates, and mechanisms in tumor progression and metastasis

A vast ensemble of extracellular proteins influences the development and progression of cancer, shaped and reshaped by a complex network of extracellular proteases. These proteases, belonging to the distinct classes of metalloproteases, serine proteases, cysteine proteases, and aspartic proteases, play a critical role in cancer. They often become dysregulated in cancer, with increases in pathological protease activity frequently driven by the loss of normal latency controls, diminished regulation by endogenous protease inhibitors, and changes in localization. Dysregulated proteases accelerate tumor progression and metastasis by degrading protein barriers within the extracellular matrix (ECM), stimulating tumor growth, reactivating dormant tumor cells, facilitating tumor cell escape from immune surveillance, and shifting stromal cells toward cancer-promoting behaviors through the precise proteolysis of specific substrates to alter their functions. These crucial substrates include ECM proteins and proteoglycans, soluble proteins secreted by tumor and stromal cells, and extracellular domains of cell surface proteins, including membrane receptors and adhesion proteins. The complexity of the extracellular protease web presents a significant challenge to untangle. Nevertheless, technological strides in proteomics, chemical biology, and the development of new probes and reagents are enabling progress and advancing our understanding of the pivotal importance of extracellular proteolysis in cancer.

Landscape of In Silico Tools for Modeling Covalent Modification of Proteins: A Review on Computational Covalent Drug Discovery

Covalent drug discovery has been a challenging research area given the struggle of finding a sweet balance between selectivity and reactivity for these drugs, the lack of which often leads to off-target activities and hence undesirable side effects. However, there has been a resurgence in covalent drug design following the success of several covalent drugs such as boceprevir (2011), ibrutinib (2013), neratinib (2017), dacomitinib (2018), zanubrutinib (2019), and many others. Design of covalent drugs includes many crucial factors, where "evaluation of the binding affinity" and "a detailed mechanistic understanding on covalent inhibition" are at the top of the list. Well-defined experimental techniques are available to elucidate these factors; however, often they are expensive and/or time-consuming and hence not suitable for high throughput screens. Recent developments in in silico methods provide promise in this direction. In this report, we review a set of recent publications that focused on developing and/or implementing novel in silico techniques in "Computational Covalent Drug Discovery (CCDD)". We also discuss the advantages and disadvantages of these approaches along with what improvements are required to make it a great tool in medicinal chemistry in the near future.

Unraveling the Reaction Mechanism of Russell's Viper Venom Factor X Activator: A Paradigm for the Reactivity of Zinc Metalloproteinases?

Snake venom metalloproteinases (SVMPs) are important drug targets against snakebite envenoming, the neglected tropical disease with the highest mortality worldwide. Here, we focus on Russell's viper (Daboia russelii), one of the "big four" snakes of the Indian subcontinent that, together, are responsible for ca. 50,000 fatalities annually. The "Russell's viper venom factor X activator" (RVV-X), a highly toxic metalloproteinase, activates the blood coagulation factor X (FX), leading to the prey's abnormal blood clotting and death. Given its tremendous public health impact, the WHO recognized an urgent need to develop efficient, heat-stable, and affordable-for-all small-molecule inhibitors, for which a deep understanding of the mechanisms of action of snake's principal toxins is fundamental. In this study, we determine the catalytic mechanism of RVV-X by using a density functional theory/molecular mechanics (DFT:MM) methodology to calculate its free energy profile. The results showed that the catalytic process takes place via two steps. The first step involves a nucleophilic attack by an in situ generated hydroxide ion on the substrate carbonyl, yielding an activation barrier of 17.7 kcal·mol-1, while the second step corresponds to protonation of the peptide nitrogen and peptide bond cleavage with an energy barrier of 23.1 kcal·mol-1. Our study shows a unique role played by Zn2+ in catalysis by lowering the pKa of the Zn2+-bound water molecule, enough to permit the swift formation of the hydroxide nucleophile through barrierless deprotonation by the formally much less basic Glu140. Without the Zn2+ cofactor, this step would be rate-limiting.

High-throughput computational investigation of protein electrostatics and cavity for SAM-dependent methyltransferases

S-adenosyl methionine (SAM)-dependent methyl transferases (MTases) are a ubiquitous class of enzymes catalyzing dozens of essential life processes. Despite targeting a large space of substrates with diverse intrinsic reactivity, SAM MTases have similar catalytic efficiency. While understanding of MTase mechanism has grown tremendously through the integration of structural characterization, kinetic assays, and multiscale simulations, it remains elusive how these enzymes have evolved to fit the diverse chemical needs of their respective substrates. In this work, we performed a high-throughput molecular modeling analysis of 91 SAM MTases to better understand how their properties (i.e., electric field [EF] strength and active site volumes) help achieve similar catalytic efficiency toward substrates of different reactivity. We found that EF strengths have largely adjusted to make the target atom a better methyl acceptor. For MTases that target RNA/DNA and histone proteins, our results suggest that EF strength accommodates formal hybridization state and variation in cavity volume trends with diversity of substrate classes. Metal ions in SAM MTases contribute negatively to EF strength for methyl donation and enzyme scaffolds tend to offset these contributions.

Multiscale Simulations of the Covalent Inhibition of the SARS-CoV-2 Main Protease: Four Compounds and Three Reaction Mechanisms

We report the results of computational modeling of the reactions of the SARS-CoV-2 main protease (MPro) with four potential covalent inhibitors. Two of them, carmofur and nirmatrelvir, have shown experimentally the ability to inhibit MPro. Two other compounds, X77A and X77C, were designed computationally in this work. They were derived from the structure of X77, a non-covalent inhibitor forming a tight surface complex with MPro. We modified the X77 structure by introducing warheads capable of reacting with the catalytic cysteine residue in the MPro active site. The reaction mechanisms of the four molecules with MPro were investigated by quantum mechanics/molecular mechanics (QM/MM) simulations. The results show that all four compounds form covalent adducts with the catalytic cysteine Cys 145 of MPro. From the chemical perspective, the reactions of these four molecules with MPro follow three distinct mechanisms. The reactions are initiated by a nucleophilic attack of the thiolate group of the deprotonated cysteine residue from the catalytic dyad Cys145-His41 of MPro. In the case of carmofur and X77A, the covalent binding of the thiolate to the ligand is accompanied by the formation of the fluoro-uracil leaving group. The reaction with X77C follows the nucleophilic aromatic substitution SNAr mechanism. The reaction of MPro with nirmatrelvir (which has a reactive nitrile group) leads to the formation of a covalent thioimidate adduct with the thiolate of the Cys145 residue in the enzyme active site. Our results contribute to the ongoing search for efficient inhibitors of the SARS-CoV-2 enzymes.

Elucidation of novel compounds and epitope-based peptide vaccine design against C30 endopeptidase regions of SARS-CoV-2 using immunoinformatics approaches

There has been progressive improvement in immunoinformatics approaches for epitope-based peptide design. Computational-based immune-informatics approaches were applied to identify the epitopes of SARS-CoV-2 to develop vaccines. The accessibility of the SARS-CoV-2 protein surface was analyzed, and hexa-peptide sequences (KTPKYK) were observed having a maximum score of 8.254, located between amino acids 97 and 102, whereas the FSVLAC at amino acids 112 to 117 showed the lowest score of 0.114. The surface flexibility of the target protein ranged from 0.864 to 1.099 having amino acid ranges of 159 to 165 and 118 to 124, respectively, harboring the FCYMHHM and YNGSPSG hepta-peptide sequences. The surface flexibility was predicted, and a 0.864 score was observed from amino acids 159 to 165 with the hepta-peptide (FCYMHHM) sequence. Moreover, the highest score of 1.099 was observed between amino acids 118 and 124 against YNGSPSG. B-cell epitopes and cytotoxic T-lymphocyte (CTL) epitopes were also identified against SARS-CoV-2. In molecular docking analyses, -0.54 to -26.21 kcal/mol global energy was observed against the selected CTL epitopes, exhibiting binding solid energies of -3.33 to -26.36 kcal/mol. Based on optimization, eight epitopes (SEDMLNPNY, GSVGFNIDY, LLEDEFTPF, DYDCVSFCY, GTDLEGNFY, QTFSVLACY, TVNVLAWLY, and TANPKTPKY) showed reliable findings. The study calculated the associated HLA alleles with MHC-I and MHC-II and found that MHC-I epitopes had higher population coverage (0.9019% and 0.5639%) than MHC-II epitopes, which ranged from 58.49% to 34.71% in Italy and China, respectively. The CTL epitopes were docked with antigenic sites and analyzed with MHC-I HLA protein. In addition, virtual screening was conducted using the ZINC database library, which contained 3,447 compounds. The 10 top-ranked scrutinized molecules (ZINC222731806, ZINC077293241, ZINC014880001, ZINC003830427, ZINC030731133, ZINC003932831, ZINC003816514, ZINC004245650, ZINC000057255, and ZINC011592639) exhibited the least binding energy (-8.8 to -7.5 kcal/mol). The molecular dynamics (MD) and immune simulation data suggest that these epitopes could be used to design an effective SARS-CoV-2 vaccine in the form of a peptide-based vaccine. Our identified CTL epitopes have the potential to inhibit SARS-CoV-2 replication.

Nitriles: an attractive approach to the development of covalent inhibitors

Nitriles have broad applications in medicinal chemistry, with more than 60 small molecule drugs on the market containing the cyano functional group. In addition to the well-known noncovalent interactions that nitriles can perform with macromolecular targets, they are also known to improve drug candidates' pharmacokinetic profiles. Moreover, the cyano group can be used as an electrophilic warhead to covalently bind an inhibitor to a target of interest, forming a covalent adduct, a strategy that can present benefits over noncovalent inhibitors. This approach has gained much notoriety in recent years, mainly with diabetes and COVID-19-approved drugs. Nevertheless, the application of nitriles in covalent ligands is not restricted to it being the reactive center, as it can also be employed to convert irreversible inhibitors into reversible ones, a promising strategy for kinase inhibition and protein degradation. In this review, we introduce and discuss the roles of the cyano group in covalent inhibitors, how to tune its reactivity and the possibility of achieving selectivity only by replacing the warhead. Finally, we provide an overview of nitrile-based covalent compounds in approved drugs and inhibitors recently described in the literature.

Exploring the Phospholipid Transport Mechanism of ATP8A1-CDC50

P4-ATPase translocates lipids from the exoplasmic to the cytosolic plasma membrane leaflet to maintain lipid asymmetry distribution in eukaryotic cells. P4-ATPase is associated with severe neurodegenerative and metabolic diseases such as neurological and motor disorders. Thus, it is important to understand its transport mechanism. However, even with progress in X-ray diffraction and cryo-electron microscopy techniques, it is difficult to obtain the dynamic information of the phospholipid transport process in detail. There are still some problems required to be resolved: (1) when does the lipid transport happen? (2) How do the key residues on the transmembrane helices contribute to the free energy of important states? In this work, we explore the phospholipid transport mechanism using a coarse-grained model and binding free energy calculations. We obtained the free energy landscape by coupling the protein conformational changes and the phospholipid transport event, taking ATP8A1-CDC50 (the typical subtype of P4-ATPase) as the research object. According to the results, we found that the phospholipid would bind to the ATP8A1-CDC50 at the early stage when ATP8A1-CDC50 changes from E2P to E2Pi-PL state. We also found that the electrostatic effects play crucial roles in the phospholipid transport process. The information obtained from this work could help us in designing novel drugs for P-type flippase disorders.

Antibacterial Activity of Ebselen

Ebselen is a low-molecular-weight organoselenium compound that has been broadly studied for its antioxidant, anti-inflammatory, and cytoprotective properties. These advantageous properties were initially associated with mimicking the activity of selenoprotein glutathione peroxidase, but the biomedical impact of this compound appear to be far more complex. Ebselen serves as a substrate or inhibitor with multiple protein/enzyme targets, whereas inhibition typically originates from the covalent modification of cysteine residues by opening the benzisoselenazolone ring and S-Se bond formation. The inhibition of enzymes of various classes and origins has been associated with substantial antimicrobial potential among other activities. In this contribution, we summarize the current state of the art regarding the antibacterial activity of ebselen. This activity, alone and in combination with commercial pharmaceuticals, against pathogens, including those resistant to drugs, is presented, together with the molecular mechanism behind the reactivity. The specific inactivation of thioredoxin reductase, bacterial toxins, and other resistance factors is considered to have certain therapeutic implications. Synergistic action and sensitization to common antibiotics assisted with the use of ebselen appear to be promising directions in the treatment of persistent infections.

The catalytic mechanism of the mitochondrial methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2)

Methylenetetrahydrofolate dehydrogenase/cyclohydrolase (MTHFD2) is a new drug target that is expressed in cancer cells but not in normal adult cells, which provides an Achilles heel to selectively kill cancer cells. Despite the availability of crystal structures of MTHFD2 in the inhibitor- and cofactor-bound forms, key information is missing due to technical limitations, including (a) the location of absolutely required Mg2+ ion, and (b) the substrate-bound form of MTHFD2. Using computational modeling and simulations, we propose that two magnesium ions are present at the active site whereby (i) Arg233, Asp225, and two water molecules coordinate MgA2+, while MgA2+ together with Arg233 stabilize the inorganic phosphate (Pi); (ii) Asp168 and three water molecules coordinate MgB2+, and MgB2+ further stabilizes Pi by forming a hydrogen bond with two oxygens of Pi; (iii) Arg201 directly coordinates the Pi; and (iv) through three water-mediated interactions, Asp168 contributes to the positioning and stabilization of MgA2+, MgB2+ and Pi. Our computational study at the empirical valence bond level allowed us also to elucidate the detailed reaction mechanisms. We found that the dehydrogenase activity features a proton-coupled electron transfer with charge redistribution connected to the reorganization of the surrounding water molecules which further facilitates the subsequent cyclohydrolase activity. The cyclohydrolase activity then drives the hydration of the imidazoline ring and the ring opening in a concerted way. Furthermore, we have uncovered that two key residues, Ser197/Arg233, are important factors in determining the cofactor (NADP+/NAD+) preference of the dehydrogenase activity. Our work sheds new light on the structural and kinetic framework of MTHFD2, which will be helpful to design small molecule inhibitors that can be used for cancer treatment.

Thiosemicarbazone derivatives: Evaluation as cruzipain inhibitors and molecular modeling study of complexes with cruzain

The development of cruzipain inhibitors represents one of the most attractive challenges in the search for drugs for the treatment of Chagas disease. A recombinant form of this enzyme, cruzain, has been crystallized with numerous inhibitors, excluding thiosemicarbazones. These compounds have been established as potent inhibitors of cruzain, although there is very little data in the literature of thiosemicarbazones tested on cruzipain. In this work, we present the results of the evaluation of eleven thiosemicarbazones on cruzipain, isolated from T. cruzi epimastigotes, six of them previously evaluated on cruzain. For these latter, we studied through computational methods, the mode of interaction with the active site of cruzain and the contribution of geometric parameters to the possible mechanism of action involved in the observed inhibition. Finally, from some geometric parameters analyzed on modeled TSC-cruzain complexes, a semi-quantitative relationship was established that could explain the inhibitory activity of thiosemicarbazones on cruzipain, the enzyme actually present in the parasite.

Pairing structural reconstruction with catalytic competence to evaluate the mechanisms of key enzymes in the folate-mediated one-carbon pathway

Mammalian metabolism comprises a series of interlinking pathways that include two major cycles: the folate and methionine cycles. The folate-mediated metabolic cycle uses several oxidation states of tetrahydrofolate to carry activated one-carbon units to be readily used and interconverted within the cell. They are required for nucleotide synthesis, methylation and metabolism, and particularly for proliferation of cancer cells. Based on the latest progress in genome-wide CRISPR loss-of-function viability screening of 789 cell lines, we focus on the most cancer-dependent enzymes in this pathway, especially those that are hyperactivated in cancer, to provide new insight into the chemical basis for cancer drug development. Since the complete 3D structure of several of these enzymes of the one-carbon pathway in their active form are not available, we used homology modelling integrated with the interpretation of the reaction mechanism. In addition, have reconstructed the most likely scenario for the reactions taking place paired with their catalytic competence that provides a testable framework for this pathway.

Structural implication of substrate binding by peptidoglycan remodeling enzyme MepS

Constant remodeling is necessary for bacterial cell growth and bacterial morphogenesis; peptidoglycan (PG) is a crucial component in this process. Murein DD-endopeptidase (MepS), initially annotated as Spr from E. coli K12, is a NlpC/P60 family endopeptidase, which cleaves the meso-diaminopimelate (DAP)-D-Ala peptide bond of PG. The Cys68, His119, His131 triad form the active site residues. MepS has autolytic activity, which is strictly regulated by a periplasmic degradation system comprising the NlpI/Prc protease complex. MepS is essential for maintaining the cell viability, and therefore, it is a potential target for developing antibiotics. This study aimed to understand the structural basis of substrate recognition and degradation. We determined the high-resolution structures of MepS, after mutating Cys68 to serine (MepS-C68S) to improve stability. We further found that citrate and L-malate molecules bind to the active site of MepS-C68S; this is in line with the recurrent observation of organic acids binding to PG endopeptidases. The presence of conserved residues on the surface revealed the potential peptide binding sites of MepS. We modelled a cross-linked peptide model of meso-DAP-D-Ala-meso-DAP, bound to the active site groove of MepS-C68S. Two conserved tyrosine residues, Tyr56 and Tyr147 appeared to be essential for the recognition of peptides. Our structural discoveries could provide insights for the design of novel antibiotics targeting MepS.

Predicting the Relative Binding Affinity for Reversible Covalent Inhibitors by Free Energy Perturbation Calculations

Covalent inhibitors are assuming central importance in drug discovery projects, especially in this pandemic scenario. Many research groups have focused their attention on inhibiting viral proteases or human proteases such as cathepsin L (hCatL). The inhibition of these critical enzymes may impair viral replication. However, molecular modeling of covalent ligands is challenging since covalent and noncovalent ligand-bound states must be considered in the binding process. In this work, we evaluated the suitability of free energy perturbation (FEP) calculations as a tool for predicting the binding affinity of reversible covalent inhibitors of hCatL. Our strategy relies on the relative free energy calculated for both covalent and noncovalent complexes and the free energy changes have been compared with experimental data for eight nitrile-based inhibitors, including three new inhibitors of hCatL. Our results demonstrate that the covalent complex can be employed to properly rank the inhibitors. Nevertheless, a comparison of the free energy changes in both noncovalent and covalent states is valuable to interpret the effect triggered by the formation of the covalent bond on the interactions played by functional groups distant from the warhead. Overall, FEP can be employed as a powerful predictor tool in developing and understanding the activity of reversible covalent inhibitors.

Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design

The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.

Synthetic Biology and Biocomputational Approaches for Improving Microbial Endoglucanases toward Their Innovative Applications

Microbial endoglucanases belonging to the β-1-4 glycosyl hydrolase family are useful enzymes due to their vast industrial applications in pulp and paper industries and biorefinery. They convert lignocellulosic substrates to soluble sugars and help in the biodegradation process. Various biocomputational tools can be utilized to understand the catalytic activity, reaction kinetics, complexity of active sites, and chemical behavior of enzyme complexes in reactions. This might be helpful in increasing productivity and cost reduction in industries. The present review gives an overview of some interesting aspects of enzyme design, including computational techniques such as molecular dynamics simulation, homology modeling, mutational analysis, etc., toward enhancing the quality of these enzymes. Moreover, the review also covers the aspects of synthetic biology, which could be helpful in faster and reliable development of useful enzymes with desired characteristics and applications. Finally, the review also deciphers the utilization of endoglucanases in biodegradation and emphasizes the use of diversified protein engineering tools and the modification of metabolic pathways for enzyme engineering.