Three Decades of Targeting Falcipains to Develop Antiplasmodial Agents: What have we Learned and What can be Done Next?

Malaria is a devastating infectious disease that affects large swathes of human populations across the planet's tropical regions. It is caused by parasites of the genus Plasmodium, with Plasmodium falciparum being responsible for the most lethal form of the disease. During the intraerythrocytic stage in the human hosts, malaria parasites multiply and degrade hemoglobin (Hb) using a battery of proteases, which include two cysteine proteases, falcipains 2 and 3 (FP-2 and FP-3). Due to their role as major hemoglobinases, FP-2 and FP-3 have been targeted in studies aiming to discover new antimalarials and numerous inhibitors with activity against these enzymes, and parasites in culture have been identified. Nonetheless, cross-inhibition of human cysteine cathepsins remains a serious hurdle to overcome for these compounds to be used clinically. In this article, we have reviewed key functional and structural properties of FP-2/3 and described different compound series reported as inhibitors of these proteases during decades of active research in the field. Special attention is also paid to the wide range of computer-aided drug design (CADD) techniques successfully applied to discover new active compounds. Finally, we provide guidelines that, in our understanding, will help advance the rational discovery of new FP-2/3 inhibitors.

Computational Study of Base-Catalyzed Thiohemiacetal Decomposition in Pseudomonas mevalonii HMG-CoA Reductase

Thiohemiacetals are key intermediates in the active sites of many enzymes catalyzing a variety of reactions. In the case of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase (PmHMGR), this intermediate connects the two hydride transfer steps where a thiohemiacetal is the product of the first hydride transfer and its breakdown forms the substrate of the second one, serving as the intermediate during cofactor exchange. Despite the many examples of thiohemiacetals in a variety of enzymatic reactions, there are few studies that detail their reactivity. Here, we present computational studies on the decomposition of the thiohemiacetal intermediate in PmHMGR using both QM-cluster and QM/MM models. This reaction mechanism involves a proton transfer from the substrate hydroxyl to an anionic Glu83 followed by a C-S bond elongation stabilized by a cationic His381. The reaction provides insight into the varying roles of the residues in the active site that favor this multistep mechanism.

An insight into the interaction between α-ketoamide- based inhibitor and coronavirus main protease: A detailed in silico study

The search for therapeutic drugs that can neutralize the effects of COVID-2019 (SARS-CoV-2) infection is the main focus of current research. The coronavirus main protease (M_pro) is an attractive target for anti-coronavirus drug design. Further, α-ketoamide is proved to be very effective as a reversible covalent-inhibitor against cysteine proteases. Herein, we report on the non-covalent to the covalent adduct formation mechanism of α-ketoamide-based inhibitor with the enzyme active site amino acids by QM/SQM model (QM = quantum mechanical, SQM = semi-empirical QM). To uncover the mechanism, we focused on two approaches: a concerted and a stepwise fashion. The concerted pathway proceeds via deprotonation of the thiol of cysteine (here, Cys₁₄₅ SγH) and simultaneous reversible nucleophilic attack of sulfur onto the α-ketoamide warhead. In this work, we propose three plausible concerted pathways. On the contrary, in a traditional two-stage pathway, the first step is proton transfer from Cys₁₄₅ SγH to His₄₁ Nδ forming an ion pair, and consecutively, in the second step, the thiolate ion attacks the α-keto group to form a thiohemiketal. In this reaction, we find that the stability of the tetrahedral intermediate oxyanion/hydroxyl group plays an important role. Moreover, as the α-keto group has two faces Si or Re for the nucleophilic attack, we considered both possibilities of attack leading to S- and R-thiohemiketal. We computed the structural, electronic, and energetic parameters of all stationary points including transition states via ONIOM and pure DFT method. Additionally, to characterize covalent, weak noncovalent interaction (NCI) and hydrogen-bonds, we applied NCI-reduced density gradient (NCI-RDG) methods along with Bader's Quantum Theory of Atoms-in-Molecules (QTAIM) and natural bonding orbital (NBO) analysis.

Azanitrile Inhibitors of the SmCB1 Protease Target Are Lethal to Schistosoma mansoni: Structural and Mechanistic Insights into Chemotype Reactivity

Azapeptide nitriles are postulated to reversibly covalently react with the active-site cysteine residue of cysteine proteases and form isothiosemicarbazide adducts. We investigated the interaction of azadipeptide nitriles with the cathepsin B1 drug target (SmCB1) from Schistosoma mansoni, a pathogen that causes the global neglected disease schistosomiasis. Azadipeptide nitriles were superior inhibitors of SmCB1 over their parent carba analogs. We determined the crystal structure of SmCB1 in complex with an azadipeptide nitrile and analyzed the reaction mechanism using quantum chemical calculations. The data demonstrate that azadipeptide nitriles, in contrast to their carba counterparts, undergo a change from E- to Z-configuration upon binding, which gives rise to a highly favorable energy profile of noncovalent and covalent complex formation. Finally, azadipeptide nitriles were considerably more lethal than their carba analogs against the schistosome pathogen in culture, supporting the further development of this chemotype as a treatment for schistosomiasis.

Activation and selectivity of OTUB-1 and OTUB-2 deubiquitinylases

The ovarian tumor domain (OTU) deubiquitinylating cysteine proteases OTUB1 and OTUB2 (OTU ubiquitin aldehyde binding 1 and 2) are representative members of the OTU subfamily of deubiquitinylases. Deubiquitinylation critically regulates a multitude of important cellular processes, such as apoptosis, cell signaling, and growth. Moreover, elevated OTUB expression has been observed in various cancers, including glioma, endometrial cancer, ovarian cancer, and breast cancer. Here, using molecular dynamics simulation approaches, we found that both OTUB1 and OTUB2 display a catalytic triad characteristic of proteases but differ in their configuration and protonation states. The OTUB1 protein had a prearranged catalytic site, with strong electrostatic interactions between the active-site residues His²⁶⁵ and Asp²⁶⁷ In OTUB2, however, the arrangement of the catalytic triad was different. In the absence of ubiquitin, the neutral states of the catalytic-site residues in OTUB2 were more stable, resulting in larger distances between these residues. Only upon ubiquitin binding did the catalytic triad in OTUB2 rearrange and bring the active site into a catalytically feasible state. An analysis of water access channels revealed only a few diffusion trajectories for the catalytically active form of OTUB1, whereas in OTUB2 the catalytic site was solvent-accessible, and a larger number of water molecules reached and left the binding pocket. Interestingly, in OTUB2, the catalytic residues His²²⁴ and Asn²²⁶ formed a stable hydrogen bond. We propose that the observed differences in activation kinetics, protonation states, water channels, and active-site accessibility between OTUB1 and OTUB2 may be relevant for the selective design of OTU inhibitors.

Distinct Roles of Catalytic Cysteine and Histidine in the Protease and Ligase Mechanisms of Human Legumain As Revealed by DFT-Based QM/MM Simulations

The cysteine protease enzyme legumain hydrolyzes peptide bonds with high specificity after asparagine and under more acidic conditions after aspartic acid [BakerE. N.J. Mol. Biol.1980, 141, 441−4847003158; BakerE. N.; J. Mol. Biol.1977, 111, 207–210859183; DrenthJ.; Biochemistry1976, 15, 3731–3738952885; MenardR.; J. Cell. Biochem.1994, 137; PolgarL.Eur. J. Biochem.1978, 88, 513–521689035; StorerA. C.; Methods Enzymol.1994, 244, 486–5007845227. Remarkably, legumain additionally exhibits ligase activity that prevails at pH > 5.5. The atomic reaction mechanisms including their pH dependence are only partly understood. Here we present a density functional theory (DFT)-based quantum mechanics/molecular mechanics (QM/MM) study of the detailed reaction mechanism of both activities for human legumain in solution. Contrasting the situation in other papain-like proteases, our calculations reveal that the active site Cys189 must be present in the protonated state for a productive nucleophilic attack and simultaneous rupture of the scissile peptide bond, consistent with the experimental pH profile of legumain-catalyzed cleavages. The resulting thioester intermediate (INT1) is converted by water attack on the thioester into a second intermediate, a diol (INT2), which is released by proton abstraction by Cys189. Surprisingly, we found that ligation is not the exact reverse of the proteolysis but can proceed via two distinct routes. Whereas the transpeptidation route involves aminolysis of the thioester (INT1), at pH 6 a cysteine-independent, histidine-assisted ligation route was found. Given legumain’s important roles in immunity, cancer, and neurodegenerative diseases, our findings open up possibilities for targeted drug design in these fields.

Modeling covalent-modifier drugs

In this review, we present a summary of how computer modeling has been used in the development of covalent-modifier drugs. Covalent-modifier drugs bind by forming a chemical bond with their target. This covalent binding can improve the selectivity of the drug for a target with complementary reactivity and result in increased binding affinities due to the strength of the covalent bond formed. In some cases, this results in irreversible inhibition of the target, but some targeted covalent inhibitor (TCI) drugs bind covalently but reversibly. Computer modeling is widely used in drug discovery, but different computational methods must be used to model covalent modifiers because of the chemical bonds formed. Structural and bioinformatic analysis has identified sites of modification that could yield selectivity for a chosen target. Docking methods, which are used to rank binding poses of large sets of inhibitors, have been augmented to support the formation of protein-ligand bonds and are now capable of predicting the binding pose of covalent modifiers accurately. The pK_a's of amino acids can be calculated in order to assess their reactivity towards electrophiles. QM/MM methods have been used to model the reaction mechanisms of covalent modification. The continued development of these tools will allow computation to aid in the development of new covalent-modifier drugs. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Crystallographic, DFT and docking (cathepsin B) studies on an organotellurium(IV) compound

Some biologically active organotellurium compounds exhibit inhibitory potency against cathepsin B. In this study, an alkyl derivative, viz. [CH3(CH2)2C(I)=C(H)](nBu)TeI2, 1, has been structurally characterised by X-ray crystallography and shown to be coordinated within a C2I2 donor set. When the stereochemically active lone pair of electrons is taken into account, a distorted trigonal bipyramidal geometry results with the iodide atoms in axial positions. Both intra- and inter-molecular Te···I interactions are also noted. If all interactions are considered, the coordination geometry is based on a Ψ-pentagonal bipyramidal geometry. An unusual feature of the structure is the curving of the functionalised C5 chain. This feature has been explored by DFT methods and shown to arise as a result of close C–H···I interactions. A docking study (cathepsin B) was performed to understand the inhibition mechanism and to compare the new results with previous observations. Notably, 1 has the same pose exhibited by analogous biologically active compounds with aryl groups. Thus, the present study suggests that (alkyl)2TeX2 compounds should also be evaluated for biological activity.

Cysteine cathepsins: their role in tumor progression and recent trends in the development of imaging probes

Papain-like cysteine proteases bear an enormous potential as drug discovery targets for both infectious and systemic human diseases. The considerable progress in this field over the last two decades has also raised interest in the visualization of these enzymes in their native context, especially with regard to tumor imaging. After a short introduction to structure and general functions of human cysteine cathepsins, we highlight their importance for drug discovery and development and provide a critical update on the current state of knowledge toward their involvement in tumor progression, with a special emphasis on their role in therapy response. In accordance with a radiopharmaceutical point of view, the main focus of this review article will be the discussion of recently developed fluorescence and radiotracer-based imaging agents together with related molecular probes.

Crystal structure of the lytic CHAP(K) domain of the endolysin LysK from Staphylococcus aureus bacteriophage K

BACKGROUND

Bacteriophages encode endolysins to lyse their host cell and allow escape of their progeny. Endolysins are also active against Gram-positive bacteria when applied from the outside and are thus attractive anti-bacterial agents. LysK, an endolysin from staphylococcal phage K, contains an N-terminal cysteine-histidine dependent amido-hydrolase/peptidase domain (CHAP(K)), a central amidase domain and a C-terminal SH3b cell wall-binding domain. CHAP(K) cleaves bacterial peptidoglycan between the tetra-peptide stem and the penta-glycine bridge.

METHODS

The CHAP(K) domain of LysK was crystallized and high-resolution diffraction data was collected both from a native protein crystal and a methylmercury chloride derivatized crystal. The anomalous signal contained in the derivative data allowed the location of heavy atom sites and phase determination. The resulting structures were completed, refined and analyzed. The presence of calcium and zinc ions in the structure was confirmed by X-ray fluorescence emission spectroscopy. Zymogram analysis was performed on the enzyme and selected site-directed mutants.

RESULTS

The structure of CHAP(K) revealed a papain-like topology with a hydrophobic cleft, where the catalytic triad is located. Ordered buffer molecules present in this groove may mimic the peptidoglycan substrate. When compared to previously solved CHAP domains, CHAP(K) contains an additional lobe in its N-terminal domain, with a structural calcium ion, coordinated by residues Asp45, Asp47, Tyr49, His51 and Asp56. The presence of a zinc ion in the active site was also apparent, coordinated by the catalytic residue Cys54 and a possible substrate analogue. Site-directed mutagenesis was used to demonstrate that residues involved in calcium binding and of the proposed active site were important for enzyme activity.

CONCLUSIONS

The high-resolution structure of the CHAP(K) domain of LysK was determined, suggesting the location of the active site, the substrate-binding groove and revealing the presence of a structurally important calcium ion. A zinc ion was found more loosely bound. Based on the structure, we propose a possible reaction mechanism. Future studies will be aimed at co-crystallizing CHAP(K) with substrate analogues and elucidating its role in the complete LysK protein. This, in turn, may lead to the design of site-directed mutants with altered activity or substrate specificity.

A new method for filtering of reactive "warheads" of transition-state analog protease inhibitors

In light of the major contribution of the reactive warhead to the binding energy trend in reversible covalent transition-state analog inhibitors of serine and cysteine hydrolases, would it be possible to rationally design and quickly filter such warheads, especially for large-scale screening? The previously defined W1 and W2 covalent descriptors quantitatively account for the energetic effect of the covalent bonds reorganization, accompanying enzyme-inhibitor covalent binding. The quantum mechanically calculated W1 and W2 reflect the warhead binding energy by modeling of the enzyme-inhibitor reaction core. Here, we demonstrate the use of these descriptors for warhead filtering, and examine its scope and limitations. The W1 and W2 descriptors provide a tool for rational design of various warheads as universal building blocks of real inhibitors without the requirement of 3D structural information about the target enzyme or QSAR studies. These warheads could then be used as hit structural templates in the subsequent optimization of inhibitors recognition sites.

Application of EMBM to Structure-Based Design of Warheads for Protease Inhibitors

Most CADD tools handle non-covalent enzyme inhibitors, despite the growing interest of the pharma industry in covalent inhibitors. We have recently introduced an enzyme mechanism-based method, EMBM, as a computational tool for binding trend analysis and prediction of chemical sites (CS) of reversible covalent enzyme inhibitors. In the current study we demonstrate the utility of EMBM to structure-based applications. In this mode, the energy of the enzyme-inhibitor covalent bond is accounted for by the W1 and W2 covalent descriptors we have developed, whereas the non-covalent interactions between the inhibitor CS and the enzyme active site can be estimated directly on the 3D structure of the enzyme-inhibitor complex.

Peptidyl cyclopropenones: reversible inhibitors, irreversible inhibitors, or substrates of cysteine proteases?

Peptidyl cyclopropenones were previously introduced as selective cysteine protease reversible inhibitors. In the present study we synthesized one such peptidyl cyclopropenone and investigated its interaction with papain, a prototype cysteine protease. A set of kinetics, biochemical, HPLC, MS, and (13)C-NMR experiments revealed that the peptidyl cyclopropenone was an irreversible inhibitor of the enzyme, alkylating the catalytic cysteine. In parallel, this cyclopropenone also behaved as an alternative substrate of the enzyme, providing a product that was tentatively suggested to be either a spiroepoxy cyclopropanone or a gamma-lactone. Thus, a single family of compounds exhibits an unusual variety of activities, being reversible inhibitors, irreversible inhibitors and alternative substrates towards enzymes of the same family.

Modeling methods for studying post-translational and transcriptional modifying enzymes

Biological catalysis is a complex chemical process that involves not only electronic reorganization in the substrate but also the reorganization of the catalyst. This complexity is even larger in the case of post-transcriptional and post-translational modifications which may involve the interaction between two biomacromolecules. However, the development over the past decades of new computational methods and strategies is offering a detailed molecular picture of the catalytic event and a deep understanding of the mechanisms of chemical reactions in biological environments. Here we review the efforts made in the last years to model catalysis in post-transcriptional and post-translational processes. We stress on the advantages and problems of the different computational strategies and their applicability in different cases.

Dynamics induced by β-lactam antibiotics in the active site of Bacillus subtilis L,D-transpeptidase

β-lactams inhibit peptidoglycan polymerization by acting as suicide substrates of essential d,d-transpeptidases. Bypass of these enzymes by unrelated l,d-transpeptidases results in β-lactam resistance, although carbapenems remain unexpectedly active. To gain insight into carbapenem specificity of l,d-transpeptidases (Ldts), we solved the nuclear magnetic resonance (NMR) structures of apo and imipenem-acylated Bacillus subtilis Ldt and show that the cysteine nucleophile is present as a neutral imidazole-sulfhydryl pair in the substrate-free enzyme. NMR relaxation dispersion does not reveal any preexisting conformational exchange in the apoenzyme, and change in flexibility is not observed upon noncovalent binding of β-lactams (K(D) > 37.5 mM). In contrast, covalent modification of active cysteine by both carbapenems and 2-nitro-5-thiobenzoate induces backbone flexibility that does not result from disruption of the imidazole-sulfhydryl proton interaction or steric hindrance. The chemical step of the reaction determines enzyme specificity since no differences in drug affinity were observed.

EMBM - a new enzyme mechanism-based method for rational design of chemical sites of covalent inhibitors

We introduce an enzyme mechanism-based method (EMBM) aimed at rational design of chemical sites (CS) of reaction coordinate analog inhibitors. The energy of valence reorganization of CS, caused by the formation of the enzyme-inhibitor covalent complex, is accounted for by new covalent descriptors W1 and W2. We considered CS fragments with a carbonyl reactivity center, like in native protease substrates. The W1 and W2 descriptors are calculated quantum mechanically on small molecular clusters simulating the reaction core of the formed covalent tetrahedral complex, anionic TC(O-) or neutral TC(OH). The modeling on a reaction core allows generation of various CS and corresponding TC(O-) and TC(OH) as universal building blocks of real inhibitors and their covalent complexes with serine or cysteine hydrolases. Moreover, the approach avoids the need for 3D structure of the target enzyme, so EMBM may be used for ligand-based design. We have built a chemical site of inhibitors (CSI) databank with pairs of W1 and W2 descriptors precalculated for both CH₃O(-) and CH₃S(-) nucleophiles for every collected CS fragment. We demonstrated that contribution of a CS fragment to the binding affinity of an inhibitor depends on both its covalent reorganization during the chemical transformation and its noncovalent interactions in the enzyme active site. Consequently, prediction of inhibitors binding trend can be done only by accounting for all of these factors, using W1 and W2 in combination with noncovalent QSAR descriptors.