Hezekiel M. Kumalo

Theory and Applications of Covalent Docking in Drug Discovery: Merits and Pitfalls

The present art of drug discovery and design of new drugs is based on suicidal irreversible inhibitors. Covalent inhibition is the strategy that is used to achieve irreversible inhibition. Irreversible inhibitors interact with their targets in a time-dependent fashion, and the reaction proceeds to completion rather than to equilibrium. Covalent inhibitors possessed some significant advantages over non-covalent inhibitors such as covalent warheads can target rare, non-conserved residue of a particular target protein and thus led to development of highly selective inhibitors, covalent inhibitors can be effective in targeting proteins with shallow binding cleavage which will led to development of novel inhibitors with increased potency than non-covalent inhibitors. Several computational approaches have been developed to simulate covalent interactions; however, this is still a challenging area to explore. Covalent molecular docking has been recently implemented in the computer-aided drug design workflows to describe covalent interactions between inhibitors and biological targets. In this review we highlight: (i) covalent interactions in biomolecular systems; (ii) the mathematical framework of covalent molecular docking; (iii) implementation of covalent docking protocol in drug design workflows; (iv) applications covalent docking: case studies and (v) shortcomings and future perspectives of covalent docking. To the best of our knowledge; this review is the first account that highlights different aspects of covalent docking with its merits and pitfalls. We believe that the method and applications highlighted in this study will help future efforts towards the design of irreversible inhibitors.

Investigation of flap flexibility of β-secretase using molecular dynamic simulations

Flap motif and its dynamics were extensively reported in aspartate proteases, e.g. HIV proteases and plasmepsins. Herein, we report the first account of flap dynamics amongst different conformations of β-secretase using molecular dynamics simulation. Various parameters were proposed and a selected few were picked which could appropriately describe the flap motion. Three systems were studied, namely Free (BACEFree) and two ligand-bound conformations, which belonged to space groups P6122 (BACEBound1) and C2221 (BACEBound2), respectively and four parameters (distance between the flaps tip residue, Thr72 and Ser325, d1; dihedral angle, ϕ (Thr72-Asp32-Asp228-Ser325); TriCα angles, θ1 (Thr72-Asp32-Ser325), and θ2 (Thr72-Asp228-Ser325)) were proposed to understand the change in dynamics of flap domain and the extent of flap opening and closing. Analysis of, θ2, d1, θ1 and ϕ confirmed that the BACEFree adopted semi-open, open and closed conformations with slight twisting during flap opening. However, BACEBound1 (P6122) showed an adaptation to open conformation due to lack of hydrogen bond interaction between the ligand and flap tip residue. A slight flap twisting, ϕ (lateral twisting) was observed for BACEBound1 during flap opening which correlates with the opening of BACEFree. Contradictory to the BACEBound1, the BACEBound2 locked the flap in a closed conformation throughout the simulation due to formation of a stable hydrogen bond interaction between the flap tip residue and ligand. Analyses of all three systems highlight that d1, θ2 and ϕ can be precisely used to describe the extent of flap opening and closing concurrently with snapshots along the molecular dynamics trajectory across several conformations of β-secretase.

Heat-Shock Protein 90 (Hsp90) as Anticancer Target for Drug Discovery: An Ample Computational Perspective

There are over 100 different types of cancer, and each is classified based on the type of cell that is initially affected. If left untreated, cancer can result in serious health problems and eventually death. Recently, the paradigm of cancer chemotherapy has evolved to use a combination approach, which involves the use of multiple drugs each of which targets an individual protein. Inhibition of heat‐shock protein 90 (Hsp90) is one of the novel key cancer targets. Because of its ability to target several signaling pathways, Hsp90 inhibition emerged as a useful strategy to treat a wide variety of cancers. Molecular modeling approaches and methodologies have become ‘close counterparts’ to experiments in drug design and discovery workflows. A wide range of molecular modeling approaches have been developed, each of which has different objectives and outcomes. In this review, we provide an up‐to‐date systematic overview on the different computational models implemented toward the design of Hsp90 inhibitors as anticancer agents. Although this is the main emphasis of this review, different topics such as background and current statistics of cancer, different anticancer targets including Hsp90, and the structure and function of Hsp90 from an experimental perspective, for example, X‐ray and NMR, are also addressed in this report. To the best of our knowledge, this review is the first account, which comprehensively outlines various molecular modeling efforts directed toward identification of anticancer drugs targeting Hsp90. We believe that the information, methods, and perspectives highlighted in this report would assist researchers in the discovery of potential anticancer agents.

A comparative molecular dynamics study on BACE1 and BACE2 flap flexibility

Abstract Beta-amyloid precursor protein cleavage enzyme1 (BACE1) and beta-amyloid precursor protein cleavage enzyme2 (BACE2), members of aspartyl protease family, are close homologs and have high similarity in their protein crystal structures. However, their enzymatic properties are different, which leads to different clinical outcomes. In this study, we performed sequence analysis and all-atom molecular dynamic (MD) simulations for both enzymes in their ligand-free states in order to compare their dynamical flap behaviors. This is to enhance our understanding of the relationship between sequence, structure and the dynamics of this protein family. Sequence analysis shows that in BACE1 and BACE2, most of the ligand-binding sites are conserved, indicative of their enzymatic property as aspartyl protease members. The other conserved residues are more or less unsystematically localized throughout the structure. Herein, we proposed and applied different combined parameters to define the asymmetric flap motion; the distance, d1, between the flap tip and the flexible region; the dihedral angle, φ, to account for the twisting motion and the TriCα angle, θ2 and θ1. All four combined parameters were found to appropriately define the observed “twisting” motion during the flaps different conformational states. Additional analysis of the parameters indicated that the flaps can exist in an ensemble of conformations, i.e. closed, semi-open and open conformations for both systems. However, the behavior of the flap tips during simulations is different between BACE1 and BACE2. The BACE1 active site cavity is more spacious as compared to that of BACE2. The analysis of 10S loop and 113S loop showed a similar trend to that of flaps, with the BACE1 loops being more flexible and less stable than those of BACE2. We believe that the results, methods and perspectives highlighted in this report would assist researchers in the discovery of BACE inhibitors as potential Alzheimer’s disease therapies.

Evaluation of 1, 4, 7-Triazacyclononane (TACN) as a potential Metallo-B-Lactamase inhibitor in Enterobacteriaceae: Restoring the Activity of B-lactams

Metallo-β-lactamase producing Enterobacteriaceae are of grave clinical concern particularly as there are no Metallo-β-lactamase (MBL) inhibitors approved for clinical use. The discovery and development of MBL inhibitors to restore the efficacy of available β-lactams are thus imperative. We investigated a zinc-chelating moiety, 1, 4, 7-triazacyclononane (TACN) for its inhibitory activity against clinical carbapenem-resistant Enterobacteriaceae. Minimum inhibitory concentrations (MICs), minimum bactericidal concentrations (MBCs), serum effect, fractional inhibitory concentrations index and time-kill kinetics were performed using broth microdilution techniques according to the Clinical Laboratory Standard Institute (CSLI) guidelines. Enzyme kinetic parameters and cytotoxicity effects of TACN were determined using spectrophotometric assays. The interactions of the enzyme-TACN complex were investigated by computational studies. Meropenem regained its activity against carbapenemase-producing Enterobacteriaceae, with the MIC decreasing to 0.03 mg/L in the presence of TACN. TACN-Meropenem combinations showed bactericidal effects with MIC/MBC ratio of ≤4, and synergistic activity was observed. Human serum effects on the MICs were insignificant, and TACN was found to be non-cytotoxic at concentrations above the MIC values. Computational studies predicted that TACN inhibits MBLs by targeting their catalytic active site pockets. This was supported by its inhibition constant Ki = 0.044 µM and inactivation constant kinact= 0.0406 (min-1) demonstrating that TACN inhibits MBLs efficiently and holds promise as a potential inhibitor. Importance Carbapenem-resistant Enterobacteriaceae (CRE)-mediated infections remain a significant public health concern and have been reported as critical in the World Health Organization’s Priority Pathogens List for the Research and Development of New Antibiotics. CRE produce enzymes such as Metallo-β-lactamases (MBLs), which inactivate β-lactam antibiotics. Combination therapies involving a β-lactam antibiotic and a β-lactamase inhibitor remain a major treatment option for infections caused by β-lactamase-producing organisms. Currently, no MBL inhibitor-β-lactam combination therapy is clinically available for MBL-positive bacterial infections. Hence, developing efficient molecules capable of inhibiting these enzymes could be a promising way to overcome this phenomenon. TACN played a significant role in the inhibitory activity of the tested molecules against CREs by potentiating the carbapenem. This study demonstrated that TACN inhibits MBLs efficiently and holds promises as a potential MBLs inhibitor to help curb the global health threat posed by MBL-producing CREs.Metallo-B-lactamase (MBL)-producing Enterobacteriaceae are of grave clinical concern particularly as there are no MBL inhibitors approved for clinical use. The discovery and development of MBL inhibitors to restore the efficacy of available B-lactams is thus imperative. We investigated a zinc-chelating moiety, 1, 4, 7-triazacyclononane (TACN) for its inhibitory activity against clinical carbapenem-resistant Enterobacteriaceae (CRE). Minimum inhibitory concentrations, minimum bactericidal concentrations, serum effect, fractional inhibitory concentrations index and time-kill kinetics were performed using broth micro-dilution techniques according to the CLSI guidelines. Enzyme kinetic parameters and cytotoxicity effects of TACN were determined using spectrophotometric assays. The interactions of the enzyme-TACN complex were investigated by computational studies. Meropenem (MEM) regained its activity against carbapenemase-producing Enterobacteriaceae, with the MIC decreasing to 0.03 mg/L in the presence of TACN. TACN-MEM combinations showed bactericidal effect with MIC/MBC ratio of less than or equal to 4 and synergistic activity was observed. Human serum effects on the MICs were insignificant and TACN was found to be non-cytotoxic at concentrations above the MIC values. Computational studies predicted that TACN inhibit MBLs by targeting their catalytic active site pockets. This was supported by its inhibition constant Ki = 0.044 uM and inactivation constant kinact = 0.0406 (min -1 ) demonstrating that TACN inhibits MBLs efficiently and holds promise as a potential inhibitor. Importance of study: Carbapenem-resistant Enterobacteriaceae (CRE)-mediated infections remain a huge public health concern and have been reported as critical in the World Health Organization Priority Pathogens List for the Research and Development of New Antibiotics. CRE produce enzymes such as Metallo-B-lactamases (MBLs), which inactivate B-lactam antibiotics. Combination therapies involving a B-lactam antibiotic and a B-lactamase inhibitor remain a major treatment option for infections caused by B-lactamase-producing organisms. Currently, no MBL inhibitor-B-lactam combination therapy is clinically available for MBL-positive bacterial infections. Hence, developing efficient molecules capable of inhibiting these enzymes remain a way forward to overcome this phenomenon. TACN played a significant role in the inhibitory activity of the tested molecules against CREs by potentiating the carbapenem. This study demonstrated that TACN inhibits MBLs efficiently and holds promise as a potential MBLs inhibitor to help curb the global health threat posed by MBL-producing CREs.

In Vitro Potentiation of Carbapenems with Tannic Acid Against Carbapenemase Producing Enterobacteriaceae: Exploring Natural Products as Potential Carbapenemase Inhibitors

We hypothesized and confirmed that tannic acid (TA) reverses carbapenem resistance by inhibiting carbapenemases in class A and B carbapenemase‐producing Enterobacteriaceae.

In Vitro Antibacterial Activity of Teixobactin Derivatives on Clinically Relevant Bacterial Isolates

Methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant enterococcus (VRE) are included on the WHO high priority list of pathogens that require urgent intervention. Hence emphasis needs to be placed on developing novel class of molecules to tackle these pathogens. Teixobactin is a new class of antibiotic that has demonstrated antimicrobial activity against common bacteria. Here we examined the antimicrobial properties of three Teixobactin derivatives against clinically relevant bacterial isolates taken from South African patients. The minimum inhibitory concentration (MIC), the minimal bactericidal concentration (MBC), the effect of serum on MICs and the time-kill kinetics studies of our synthesized Teixobactin derivatives (3, 4, and 5) were ascertained following the CLSI 2017 guidelines and using the broth microdilution method. Haemolysis on red blood cells (RBCs) and cytotoxicity on peripheral blood mononuclear cells (PBMCs) were performed to determine the safety of these compounds. The MICs of 3, 4, and 5 against reference strains were 4–64 μg/ml, 2–64 μg/ml, and 0.5–64 μg/ml, respectively. The MICs observed for MRSA were (3) 32 μg/ml, (4) 2–4 μg/ml and (5) 2–4 μg/ml whilst those for VRE were (3) 8–16 μg/ml, (4) 4 μg/ml and (5) 2–16 μg/ml, respectively. In the presence of 50% human serum, there was no significant effect on the MICs. The compounds did not exhibit any effect on cell viability at their effective concentrations. Teixobactin derivatives (3, 4, and 5) inhibited bacterial growth in drug-resistant bacteria and hence emerge as potential antimicrobial agents. Molecular dynamic simulations suggested that the most dominant binding mode of Lys10-teixobactin (4) to lipid II is through the amide protons of the cycle, which is identical to data described in the literature for the natural teixobactin hence predicting the possibility of a similar mechanism of action.

Mechanism of Inhibition of Hsp90 Dimerization by Gyrase B Inhibitor Coumermycin A1 (C–A1) Revealed by Molecular Dynamics Simulations and Thermodynamic Calculations

Heat shock protein (Hsp) 90 an emerging and attracting target in the anti-HIV drug discovery process due to the key role it plays in the pathogenicity of HIV-1 virus. In this research study, long-range all-atom molecular dynamics simulations were engaged for the bound and the unbound proteins to enhance the understanding of the molecular mechanisms of the Hsp90 dimerization and inhibition. Results evidently showed that coumermycin A1 (C–A1), a recently discovered Hsp90 inhibitor, binds at the dimer’s active site of the Hsp90 protein and leads to a substantial parting between dimeric opposed residues, which include Arg591.B, Lys594.A, Ser663.A, Thr653.B, Ala665.A, Thr649.B, Leu646.B and Asn669.A. Significant differences in magnitudes were observed in radius of gyration, root-mean-square deviation and root-mean-square fluctuation, which confirms a reasonably more flexible state in the apo conformation associated with it dimerization. In contrast, the bound conformer of Hsp90 showed less flexibility. This visibly highpoints the inhibition process resulting from the binding of the ligand. These findings were further validated by principal component analysis. We believe that the detailed dynamic analyses of Hsp90 presented in this study, would give an imperative insight and better understanding to the function and mechanisms of inhibition. Furthermore, information obtained from the binding mode of the inhibitor would be of great assistance in the design of more potent inhibitors against the HIV target Hsp90.