Apigenin analogues as SARS-CoV-2 main protease inhibitors: In-silico screening approach.

The COVID-19 new variants spread rapidly all over the world, and until now scientists strive to find virus-specific antivirals for its treatment. The main protease of SARS-CoV-2 (Mpro) exhibits high structural and sequence homology to main protease of SARS-CoV (93.23% sequence identity), and their sequence alignment indicated 12 mutated/variant residues. The sequence alignment of SARS-CoV-2 main protease led to identification of only one mutated/variant residue with no significant role in its enzymatic process. Therefore, Mpro was considered as a high-profile drug target in anti-SARS-CoV-2 drug discovery. Apigenin analogues to COVID-19 main protease binding were evaluated. The detailed interactions between the analogues of Apigenin and SARS-CoV-2 Mpro inhibitors were determined as hydrogen bonds, electronic bonds and hydrophobic interactions. The binding energies obtained from the molecular docking of Mpro with Boceprevir, Apigenin, Apigenin 7-glucoside-4'-p-coumarate, Apigenin 7-glucoside-4'-trans-caffeate and Apigenin 7-O-beta-d-glucoside (Cosmosiin) were found to be -6.6, -7.2, -8.8, -8.7 and -8.0 kcal/mol, respectively. Pharmacokinetic parameters and toxicological characteristics obtained by computational techniques and Virtual ADME studies of the Apigenin analogues confirmed that the Apigenin 7-glucoside-4'-p-coumarate is the best candidate for SARS-CoV-2 Mpro inhibition.

Design of SARS-CoV-2 PLpro Inhibitors for COVID-19 Antiviral Therapy Leveraging Binding Cooperativity.

Antiviral agents that complement vaccination are urgently needed to end the COVID-19 pandemic. The SARS-CoV-2 papain-like protease (PLpro), one of only two essential cysteine proteases that regulate viral replication, also dysregulates host immune sensing by binding and deubiquitination of host protein substrates. PLpro is a promising therapeutic target, albeit challenging owing to featureless P1 and P2 sites recognizing glycine. To overcome this challenge, we leveraged the cooperativity of multiple shallow binding sites on the PLpro surface, yielding novel 2-phenylthiophenes with nanomolar inhibitory potency. New cocrystal structures confirmed that ligand binding induces new interactions with PLpro: by closing of the BL2 loop of PLpro forming a novel "BL2 groove" and by mimicking the binding interaction of ubiquitin with Glu167 of PLpro. Together, this binding cooperativity translates to the most potent PLpro inhibitors reported to date, with slow off-rates, improved binding affinities, and low micromolar antiviral potency in SARS-CoV-2-infected human cells.

NMR Spectroscopy of the Main Protease of SARS-CoV-2 and Fragment-Based Screening Identify Three Protein Hotspots and an Antiviral Fragment.

The main protease (3CLp) of the SARS-CoV-2, the causative agent for the COVID-19 pandemic, is one of the main targets for drug development. To be active, 3CLp relies on a complex interplay between dimerization, active site flexibility, and allosteric regulation. The deciphering of these mechanisms is a crucial step to enable the search for inhibitors. In this context, using NMR spectroscopy, we studied the conformation of dimeric 3CLp from the SARS-CoV-2 and monitored ligand binding, based on NMR signal assignments. We performed a fragment-based screening that led to the identification of 38 fragment hits. Their binding sites showed three hotspots on 3CLp, two in the substrate binding pocket and one at the dimer interface. F01 is a non-covalent inhibitor of the 3CLp and has antiviral activity in SARS-CoV-2 infected cells. This study sheds light on the complex structure-function relationships of 3CLp and constitutes a strong basis to assist in developing potent 3CLp inhibitors.

Inhibition of Cysteine Proteases by 6,6'-Dihydroxythiobinupharidine (DTBN) from Nuphar lutea.

The specificity of inhibition by 6,6'-dihydroxythiobinupharidine (DTBN) on cysteine proteases was demonstrated in this work. There were differences in the extent of inhibition, reflecting active site structural-steric and biochemical differences. Cathepsin S (IC50 = 3.2 µM) was most sensitive to inhibition by DTBN compared to Cathepsin B, L and papain (IC50 = 1359.4, 13.2 and 70.4 µM respectively). DTBN is inactive for the inhibition of Mpro of SARS-CoV-2. Docking simulations suggested a mechanism of interaction that was further supported by the biochemical results. In the docking results, it was shown that the cysteine sulphur of Cathepsin S, L and B was in close proximity to the DTBN thiaspirane ring, potentially forming the necessary conditions for a nucleophilic attack to form a disulfide bond. Covalent docking and molecular dynamic simulations were performed to validate disulfide bond formation and to determine the stability of Cathepsins-DTBN complexes, respectively. The lack of reactivity of DTBN against SARS-CoV-2 Mpro was attributed to a mismatch of the binding conformation of DTBN to the catalytic binding site of Mpro. Thus, gradations in reactivity among the tested Cathepsins may be conducive for a mechanism-based search for derivatives of nupharidine against COVID-19. This could be an alternative strategy to the large-scale screening of electrophilic inhibitors.

Establishing an Analogue Based In Silico Pipeline in the Pursuit of Novel Inhibitory Scaffolds against the SARS Coronavirus 2 Papain-Like Protease.

The ongoing coronavirus pandemic has been a burden on the worldwide population, with mass fatalities and devastating socioeconomic consequences. It has particularly drawn attention to the lack of approved small-molecule drugs to inhibit SARS coronaviruses. Importantly, lessons learned from the SARS outbreak of 2002-2004, caused by severe acute respiratory syndrome coronavirus 1 (SARS-CoV-1), can be applied to current drug discovery ventures. SARS-CoV-1 and SARS-CoV-2 both possess two cysteine proteases, the main protease (Mpro) and the papain-like protease (PLpro), which play a significant role in facilitating viral replication, and are important drug targets. The non-covalent inhibitor, GRL-0617, which was found to inhibit replication of SARS-CoV-1, and more recently SARS-CoV-2, is the only PLpro inhibitor co-crystallised with the recently solved SARS-CoV-2 PLpro crystal structure. Therefore, the GRL-0617 structural template and pharmacophore features are instrumental in the design and development of more potent PLpro inhibitors. In this work, we conducted scaffold hopping using GRL-0617 as a reference to screen over 339,000 ligands in the chemical space using the ChemDiv, MayBridge, and Enamine screening libraries. Twenty-four distinct scaffolds with structural and electrostatic similarity to GRL-0617 were obtained. These proceeded to molecular docking against PLpro using the AutoDock tools. Of two compounds that showed the most favourable predicted binding affinities to the target site, as well as comparable protein-ligand interactions to GRL-0617, one was chosen for further analogue-based work. Twenty-seven analogues of this compound were further docked against the PLpro, which resulted in two additional hits with promising docking profiles. Our in silico pipeline consisted of an integrative four-step approach: (1) ligand-based virtual screening (scaffold-hopping), (2) molecular docking, (3) an analogue search, and, (4) evaluation of scaffold drug-likeness, to identify promising scaffolds and eliminate those with undesirable properties. Overall, we present four novel, and lipophilic, scaffolds obtained from an exhaustive search of diverse and uncharted regions of chemical space, which may be further explored in vitro through structure-activity relationship (SAR) studies in the search for more potent inhibitors. Furthermore, these scaffolds were predicted to have fewer off-target interactions than GRL-0617. Lastly, to our knowledge, this work contains the largest ligand-based virtual screen performed against GRL-0617.

EGCG, a Green Tea Catechin, as a Potential Therapeutic Agent for Symptomatic and Asymptomatic SARS-CoV-2 Infection.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has emerged to be the greatest threat to humanity in the modern world and has claimed nearly 2.2 million lives worldwide. The United States alone accounts for more than one fourth of 100 million COVID-19 cases across the globe. Although vaccination against SARS-CoV-2 has begun, its efficacy in preventing a new or repeat COVID-19 infection in immunized individuals is yet to be determined. Calls for repurposing of existing, approved, drugs that target the inflammatory condition in COVID-19 are growing. Our initial gene ontology analysis predicts a similarity between SARS-CoV-2 induced inflammatory and immune dysregulation and the pathophysiology of rheumatoid arthritis. Interestingly, many of the drugs related to rheumatoid arthritis have been found to be lifesaving and contribute to lower COVID-19 morbidity. We also performed in silico investigation of binding of epigallocatechin gallate (EGCG), a well-known catechin, and other catechins on viral proteins and identified papain-like protease protein (PLPro) as a binding partner. Catechins bind to the S1 ubiquitin-binding site of PLPro, which might inhibit its protease function and abrogate SARS-CoV-2 inhibitory function on ubiquitin proteasome system and interferon stimulated gene system. In the realms of addressing inflammation and how to effectively target SARS-CoV-2 mediated respiratory distress syndrome, we review in this article the available knowledge on the strategic placement of EGCG in curbing inflammatory signals and how it may serve as a broad spectrum therapeutic in asymptomatic and symptomatic COVID-19 patients.

In Silico Screening of the DrugBank Database to Search for Possible Drugs against SARS-CoV-2.

Coronavirus desease 2019 (COVID-19) is responsible for more than 1.80 M deaths worldwide. A Quantitative Structure-Activity Relationships (QSAR) model is developed based on experimental pIC50 values reported for a structurally diverse dataset. A robust model with only five descriptors is found, with values of R2 = 0.897, Q2LOO = 0.854, and Q2ext = 0.876 and complying with all the parameters established in the validation Tropsha's test. The analysis of the applicability domain (AD) reveals coverage of about 90% for the external test set. Docking and molecular dynamic analysis are performed on the three most relevant biological targets for SARS-CoV-2: main protease, papain-like protease, and RNA-dependent RNA polymerase. A screening of the DrugBank database is executed, predicting the pIC50 value of 6664 drugs, which are IN the AD of the model (coverage = 79%). Fifty-seven possible potent anti-COVID-19 candidates with pIC50 values > 6.6 are identified, and based on a pharmacophore modelling analysis, four compounds of this set can be suggested as potent candidates to be potential inhibitors of SARS-CoV-2. Finally, the biological activity of the compounds was related to the frontier molecular orbitals shapes.

Protease targeted COVID-19 drug discovery and its challenges: Insight into viral main protease (Mpro) and papain-like protease (PLpro) inhibitors.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) brutally perils physical and mental health worldwide. Unavailability of effective anti-viral drug rendering global threat of COVID-19 caused by SARS-CoV-2. In this scenario, viral protease enzymes are crucial targets for drug discovery. This extensive study meticulously focused on two viral proteases such as main protease (Mpro) and papain-like protease (PLpro), those are essential for viral replication. This review provides a detail overview of the targets (Mpro and PLpro) from a structural and medicinal chemistry point of view, together with recently reported protease inhibitors. An insight into the challenges in the development of effective as well as drug like protease inhibitors is discussed. Peptidomimetic and/or covalent coronavirus protease inhibitors possessed potent and selective active site inhibition but compromised in pharmacokinetic parameters to be a drug/drug like molecule. Lead optimization of non-peptidomimetic and/or low molecular weight compounds may be a better option for oral delivery. A masterly combination of adequate pharmacokinetic properties with coronavirus protease activity as well as selectivity will provide potential drug candidates in future. This study is a part of our endeavors which surely dictates medicinal chemistry efforts to discover effective anti-viral agent for this devastating disease.

Preclinical evaluation of [11 C]L-235 as a radioligand for Positron Emission Tomography cathepsin K imaging in bone.

The cathepsin K (CatK) enzyme is abundantly expressed in osteoclasts, and CatK inhibitors have been developed for the treatment of osteoporosis. In our effort to support discovery and clinical evaluations of a CatK inhibitor, we sought to discover a radioligand to determine target engagement of the enzyme by therapeutic candidates using positron emission tomography (PET). L-235, a potent and selective CatK inhibitor, was labeled with carbon-11. PET imaging studies recording baseline distribution of [11 C]L-235, and chase and blocking studies using the selective CatK inhibitor MK-0674 were performed in juvenile and adult nonhuman primates (NHP) and ovariectomized rabbits. Retention of the PET tracer in regions expected to be osteoclast-rich compared with osteoclast-poor regions was examined. Increased retention of the radioligand was observed in osteoclast-rich regions of juvenile rabbits and NHP but not in the adult monkey or adult ovariectomized rabbit. Target engagement of CatK was observed in blocking studies with MK-0674, and the radioligand retention was shown to be sensitive to the level of MK-0674 exposure. [11 C]L-235 can assess target engagement of CatK in bone only in juvenile animals. [11 C]L-235 may be a useful tool for guiding the discovery of CatK inhibitors.

Structures of the free and inhibitors-bound forms of bromelain and ananain from Ananas comosus stem and in vitro study of their cytotoxicity.

The Ananas comosus stem extract is a complex mixture containing various cysteine ​​proteases of the C1A subfamily, such as bromelain and ananain. This mixture used for centuries in Chinese medicine, has several potential therapeutic applications as anti-cancer, anti-inflammatory and ecchymosis degradation agent. In the present work we determined the structures of bromelain and ananain, both in their free forms and in complex with the inhibitors E64 and TLCK. These structures combined with protease-substrate complexes modeling clearly identified the Glu68 as responsible for the high discrimination of bromelain in favor of substrates with positively charged residues at P2, and unveil the reasons for its weak inhibition by cystatins and E64. Our results with purified and fully active bromelain, ananain and papain show a strong reduction of cell proliferation with MDA-MB231 and A2058 cancer cell lines at a concentration of about 1 µM, control experiments clearly emphasizing the need for proteolytic activity. In contrast, while bromelain and ananain had a strong effect on the proliferation of the OCI-LY19 and HL-60 non-adherent cell lines, papain, the archetypal member of the C1A subfamily, had none. This indicates that, in this case, sequence/structure identity beyond the active site of bromelain and ananain is more important than substrate specificity.

Benzothiophene-2-carboxamide derivatives as SENPs inhibitors with selectivity within SENPs family.

The SUMO (small ubiquitin-related modifier)-specific proteases (SENPs) are responsible for the cleavage of SUMO from its target proteins, thus play important roles in the dynamic SUMOylation and deSUMOylation processes. SENPs are related to a variety of human diseases including cancer and represent a new class of potential therapeutic targets with mechanism of action that is likely to be different from that of current clinically used drugs. However, potent inhibitors that are selective within the SENPs family members still remain a challenge due to their high homology. In order to demonstrate the feasibility of developing selective inhibitors within the SENPs family, we chose SENP1/2/5 as representatives, aiming to identify inhibitors with selectivity among the members. Starting from a hit compound ZCL951 from virtual screening, a series of benzothiophene-2-carboxamide inhibitors were designed based on the protein structures of SENP1, 2, and 5. First, an unoccupied hydrophobic pocket was first identified which led to IC50 as low as 0.56 µM. Furthermore, the ethylacetate 77 gave both submicromolar inhibitory activity and 33-fold selectivity for SENP2 versus SENP5. They are the most potent and selective nonpeptidic inhibitor reported so far for the SENPs family, as far as we are aware. Their structure-activity relationship was also discussed.

Virtual screening based on molecular docking of possible inhibitors of Covid-19 main protease.

Coronavirus (COVID-19) is an enveloped RNA virus that is diversely found in humans and that has now been declared a global pandemic by the World Health Organization. Thus, there is an urgent need to develop effective therapies and vaccines against this disease. In this context, this study aimed to evaluate in silico the molecular interactions of drugs with therapeutic indications for treatment of COVID-19 (Azithromycin, Baricitinib and Hydroxychloroquine) and drugs with similar structures (Chloroquine, Quinacrine and Ruxolitinib) in docking models from the SARS-CoV-2 main protease (M-pro) protein. The results showed that all inhibitors bound to the same enzyme site, more specifically in domain III of the SARS-CoV-2 main protease. Therefore, this study allows proposing the use of baricitinib and quinacrine, in combination with azithromycin; however, these computer simulations are just an initial step for conceiving new projects for the development of antiviral molecules.

Identification of Potential Inhibitors of Cathepsin-B using Shape & Pharmacophore-based Virtual Screening, Molecular Docking and Explicit Water Thermodynamics.

Lysosome has been long understood as a vital digestive organelle. Increasing reports indicate that the lysosome also plays a crucial role in the pathogenesis of a variety of neurodegenerative diseases, including Huntington's disease, Alzheimer's disease, Parkinson's disease and amyotrophic lateral sclerosis. Abnormal protein degradation and deposition stimulated by lysosomal dysfunction may cause age-related neurodegeneration. Enormous efforts have been devoted to the development of effective therapeutics against Alzheimer's disease, the most debilitating neurodegenerative disease. Endopeptidase activity of the Cathepsin-B is associated with the pathological processes. Work presented here focuses on identification of new inhibitors against Cathepsin-B protein using diverse computational approaches together. The inhibitors identified were further tested for in-vitro activity using enzyme based assay method. The identified inhibitors provided interesting understanding on how the water thermodynamic properties along with hydrophobic, steric, electronic, and structural requirements contribute to cathepsin-B inhibitory activity. These water thermodynamic studies, may further be used in computer aided drug discovery pipeline to design and predict more potent derivatives of various scaffolds as cathepsin-B inhibitors.

Discovery of covalent enzyme inhibitors using virtual docking of covalent fragments.

Here we present a virtual docking screen of 1648 commercially available covalent fragments, and identified covalent inhibitors of cysteine protease cathepsin L. These inhibitors did not inhibit closely related protease cathepsin B. Thus, we have established virtual docking of covalent fragments as an approach to discover covalent enzyme inhibitors.

Discovery of Non-Peptidic Compounds against Chagas Disease Applying Pharmacophore Guided Molecular Modelling Approaches.

Chagas disease is one of the primary causes of heart diseases accounting to 50,000 lives annually and is listed as the neglected tropical disease. Because the currently available therapies have greater toxic effects with higher resistance, there is a dire need to develop new drugs to combat the disease. In this pursuit, the 3D QSAR ligand-pharmacophore (pharm 1) and receptor-based pharmacophore (pharm 2) search was initiated to retrieve the candidate compounds from universal natural compounds database. The validated models were allowed to map the universal natural compounds database. The obtained lead candidates were subjected to molecular docking against cysteine protease (PDB code: 1ME3) employing -Cdocker available on the discovery studio. Subsequently, two Hits have satisfied the selection criteria and were escalated to molecular dynamics simulation and binding free energy calculations. These Hits have demonstrated higher dock scores, displayed interactions with the key residues portraying an ideal binding mode complemented by mapping to all the features of pharm 1 and pharm 2. Additionally, they have rendered stable root mean square deviation (RMSD) and potential energy profiles illuminating their potentiality as the prospective antichagastic agents. The study further demonstrates the mechanism of inhibition by tetrad residues compromising of Gly23 and Asn70 holding the ligand at each ends and the residues Gly65 and Gly160 clamping the Hits at the center. The notable feature is that the Hits lie in close proximity with the residues Glu66 and Leu67, accommodating within the S1, S2 and S3 subsites. Considering these findings, the study suggests that the Hits may be regarded as effective therapeutics against Chagas disease.

Bromein, a Bromelain Inhibitor from Pineapple Stem: Structural and Functional Characteristics.

Bromelain inhibitor, "bromein", is a proteinase-inhibitor specific to the cysteine proteinase bromelain from pineapple stem. In the stem, eight bromein isoforms are known to exist, and each isoform has a short peptide (light chain) and a long one (heavy chain) with five disulfide bonds. The three-dimensional structure of the sixth isoform (bromein-6) is composed of inhibitory and stabilizing domains, and each domain contains a three-stranded antiparallel ß-sheet. The genomic sequence of a bromein precursor encodes three homologous bromein isoform domains, and each isoform domain has a signal peptide, three interchain peptides between the light chain and heavy chain, two interdomain peptides and a propeptide. Interestingly, at the protein level, bromein- 6 appears to share a similar folding and disulfide-bonding connectivity with Bowman-Birk serine proteinase inhibitors and shows weak inhibition toward chymotrypsin and trypsin. However, no significant similarity was found between them at the genomic level. This indicates that they have evolved convergently to possess such a structural similarity. To identify the essential reactive site(s) with bromelain, we investigated the inhibitory activity of 44 kinds of the single/double and insertion/ deletion mutants of bromein-6 towards stem bromelain. As a result, it was shown that both the appropriate positioning and the complete side-chain structure of Leu10 in the light chain are absolutely crucial for the inhibition, with an additional measure of importance for the preceding Pro9. Bromein and stem bromelain coexist in the acidic vacuoles of the stem tissue, and one of the key role of bromein appears to be the regulation of the bromelain activity.

Discovery of benzimidazole-based Leishmania mexicana cysteine protease CPB2.8ΔCTE inhibitors as potential therapeutics for leishmaniasis.

Chemotherapy is currently the only effective approach to treat all forms of leishmaniasis. However, its effectiveness is severely limited due to high toxicity, long treatment length, drug resistance, or inadequate mode of administration. As a consequence, there is a need to identify new molecular scaffolds and targets as potential therapeutics for the treatment of this disease. We report a small series of 1,2-substituted-1H-benzo[d]imidazole derivatives (9a-d) showing affinity in the submicromolar range (Ki  = 0.15-0.69 µM) toward Leishmania mexicanaCPB2.8ΔCTE, one of the more promising targets for antileishmanial drug design. The compounds confirmed activity in vitro against intracellular amastigotes of Leishmania infantum with the best result being obtained with derivative 9d (IC50  = 6.8 µM), although with some degree of cytotoxicity (CC50  = 8.0 µM on PMM and CC50  = 32.0 µM on MCR-5). In silico molecular docking studies and ADME-Tox properties prediction were performed to validate the hypothesis of the interaction with the intended target and to assess the drug-likeness of these derivatives.

Synthesis and biological evaluation of potential inhibitors of the cysteine proteases cruzain and rhodesain designed by molecular simplification.

Analogues of 8-chloro-N-(3-morpholinopropyl)-5H-pyrimido[5,4-b]indol-4-amine 1, a known cruzain inhibitor, were synthesized using a molecular simplification strategy. Five series of analogues were obtained: indole, pyrimidine, quinoline, aniline and pyrrole derivatives. The activity of the compounds was evaluated against the enzymes cruzain and rhodesain as well as against Trypanosoma cruzi amastigote and trypomastigote forms. The 4-aminoquinoline derivatives showed promising activity against both enzymes, with IC50 values ranging from 15 to 125µM. These derivatives were selective inhibitors for the parasitic proteases, being unable to inhibit mammalian cathepsins B and S. The most active compound against cruzain (compound 5a; IC50=15µM) is considerably more synthetically accessible than 1, while retaining its ligand efficiency. As observed for the original lead, compound 5a was shown to be a competitive enzyme inhibitor. In addition, it was also active against T. cruzi (IC50=67.7µM). Interestingly, the pyrimidine derivative 4b, although inactive in enzymatic assays, was highly active against T. cruzi (IC50=3.1µM) with remarkable selectivity index (SI=128) compared to uninfected fibroblasts. Both 5a and 4b exhibit drug-like physicochemical properties and are predicted to have a favorable ADME profile, therefore having great potential as candidates for lead optimization in the search for new drugs to treat Chagas disease.

Structural and functional studies on a variant of cystatin purified from brain of Capra hircus.

Cystatins, known for their ubiquitous presence in mammalian system are thiol protease inhibitors serving important physiological functions. Here, we present a variant of cystatin isolated from brain of Capra hircus (goat) which is glycosylated but lacks disulphide bonds. Caprine brain cystatin (CBC) was isolated using alkaline treatment, ammonium sulphate fractionation (40-60%) and gel filtration chromatography on Sephacryl S-100HR column with an overall yield of 26.29% and 322-fold purification. The inhibitor gave a molecular mass of ~44 kDa as determined by SDS-PAGE and gel filtration behaviour. The Stokes radius and diffusion coefficient of CBC were 27.14 Å and 8.18 × 10-7 cm2 s-1, respectively. Kinetic data revealed that CBC inhibited thiol proteases reversibly and competitively, with the highest inhibition towards papain (Ki = 4.10 nM) followed by ficin and bromelain. CBC possessed 34.7% α-helical content as observed by CD spectroscopy. UV, fluorescence, CD and FTIR spectroscopy revealed significant conformational change upon CBC-papain complex formation. Isothermal titration calorimetry (ITC) was used to measure the thermodynamic parameters - ΔH, ΔS, ΔG along with N (binding stoichiometry) for CBC-papain complex formation. Binding stoichiometry (N = .97 ± .07 sites) for the CBC-papain complex indicates that cystatin is surrounded by nearly one papain molecule. Negative ΔH (-5.78 kcal mol-1) and positive ΔS (11.01 cal mol-1 deg-1) values suggest that the interaction between CBC and papain is enthalpically as well as entropically favoured process. The overall negative ΔG (-9.19 kcal mol-1) value implies a spontaneous CBC-papain interaction.