Ursolic acid as a potential inhibitor of Mycobacterium tuberculosis cytochrome bc1 oxidase-a molecular modelling perspective

The escalating burden of tuberculosis disease and drastic effects of current medicine has stimulated a search for alternative drugs. A medicinal plant Warburgia salutaris has been reported to possess inhibitory properties against M. tuberculosis. In this study, we apply computational methods to investigate the probability of W. salutaris compounds as potential inhibitors of M. tuberculosis QcrB protein. We performed molecular docking, molecular dynamics simulations, radius of gyration, principal component analysis (PCA), and molecular mechanics-generalized born surface area (MM-GBSA) binding-free energy calculations in explicit solvent to achieve our objective. The results suggested that ursolic acid (UA) and ursolic acid acetate (UAA) could serve as preferred potential inhibitors of mycobacterial QcrB compared to lansoprazole sulphide (LSPZ) and telacebec (Q203)-UA and UAA have a higher binding affinity to QcrB compared to LSPZ and Q203 drugs. UA binding affinity is attributed to hydrogen bond formation with Val120, Arg364 and Arg366, and largely resonated from van der Waals forces resulting from UA interactions with hydrophobic amino acids in its vicinity. UAA binds to the porphyrin ring binding site with higher binding affinity compared to LSPZ. The binding affinity results primarily from van der Waals forces between UAA and hydrophobic residues of QcrB in the porphyrin ring binding site where UAA binds competitively. UA and UAA formed stable complexes with the protein with reduced overall residue mobility, consequently supporting the magnitude of binding affinity of the respective ligands. UAA could potentially compete with the porphyrin ring for the binding site and deprive the mycobacterial cell from oxygen, consequently disturbing mycobacterial oxygen-dependent metabolic processes. Therefore, discovery of a compound that competes with porphyrin ring for the binding site may be useful in QcrB pharmocological studies. UA proved to be a superior compound, although its estimated toxicity profile revealed UA to be hepatotoxic within acceptable parameters. Although preliminary findings of this report still warrant experimental validation, they could serve as a baseline for the development of new anti-tubercular drugs from natural resources that target QcrB.

Impact of Drug Repurposing on SARS-Cov-2 Main Protease

The recent emergence of the severe acute respiratory disease caused by a novel coronavirus remains a concern posing many challenges to public health and the global economy. The resolved crystal structure of the main protease of SARS-CoV-2 or SCV2 (M^pro) has led to its identification as an attractive target for designing potent antiviral drugs. Herein, we provide a comparative molecular impact of hydroxychloroquine (HCQ), remdesivir, and β-D-N⁴-Hydroxycytidine (NHC) binding on SCV2 M^pro using various computational approaches like molecular docking and molecular dynamics (MD) simulation. Data analyses showed that HCQ, remdesivir, and NHC binding to SARS-CoV-2 M^pro decrease the protease loop capacity to fluctuate. These binding influences the drugs’ optimum orientation in the conformational space of SCV2 M^pro and produce noticeable steric effects on the interactive residues. An increased hydrogen bond formation was observed in SCV2 M^pro–NHC complex with a decreased receptor residence time during NHC binding. The binding mode of remdesivir to SCV2 M^pro differs from other drugs having van der Waals interaction as the force stabilizing protein–remdesivir complex. Electrostatic interaction dominates in the SCV2 M^pro−HCQ and SCV2 M^pro–NHC. Residue Glu166 was highly involved in the stability of remdesivir and NHC binding at the SCV2 M^pro active site, while Asp187 provides stability for HCQ binding.

Identifying the analogues of berberine as promising antitubercular drugs targeting Mtb-FtsZ polymerisation through ligand-based virtual screening and molecular dynamics simulations

Berberine, an active compound in the extract of golden seal (an age-long remedy for many infections) has been confirmed to be responsible for the extract's activity against multi-drug resistant strain of Mycobacterium tuberculosis. There is no available study that shows the exact target of berberine in M tuberculosis, although it is confirmed that berberine inhibits the polymerisation of filamentous temperature-sensitive mutant Z (FtsZ), an important bacteria cytokinesis protein, in Escherichia coli, suggesting that FtsZ could as well be the target of berberine in M tuberculosis. In this study, we carried out ligand-based virtual screening to identify analogues of berberine followed by molecular dynamics (MD) simulations of the complexes of Mtb-FtsZ with berberine (berb1) and the five selected analogues (berb9 [ZINC1709414], berb37 [ZINC238749993], berb38 [ZINC13509022], berb43 [ZINC14765594], and berb48 [ZINC238758595]). Post-MD analyses such as binding free energy, RMSD, RMSF, RoG and hydrogen bond lifetime analysis were used to understand the interactions between these ligands and the receptor. The results suggested that Mtb-FtsZ could likely be the target of berberine in M tuberculosis as it forms a stable complex coupled with a significantly high binding energy. The study also identified other potential inhibitors of MTB-FtsZ polymerisation. Berb38 specifically showed greater interaction with the residues at the binding site of the protein, forming a far more stable complex with the receptor than any of the other compounds under investigation, including berberine itself. ADME properties calculations also predicted all the ligands to be bioactive as orally administered drugs.

A consequence of drug targeting of aminoacyl-tRNA synthetases in Mycobacteriumtuberculosis

Drug-resistant Mycobacterium tuberculosis poses a great threat to public health and remains one of the red-flag tagged infectious diseases, with the tendency of comorbidity with other disease conditions such as HIV/AIDS. This perhaps is responsible for redoubling of effort in tuberculosis research and continuous change in patient management to optimize the drug therapy. Aminoacyl-tRNA synthetases are essential enzymes in M. tuberculosis that catalyse the transfer of a particular amino acid to its corresponding specific tRNA to form an aminoacyl-tRNA. These enzymes are believed to be novel antibacterial, antifungal and antiparasitic drug targets because of their role in the process of protein translation. Therefore, their existence as a compliment of M. tuberculosis has attracted a lot of research interest with the aim of curbing the scourge and provide the most effective drug in the treatment of tuberculosis. This leads to the discovery of a pool of aminoacyl-tRNA synthetases with their essential inhibitors. This review seeks to articulate the current advances in the development of new TB drugs exhibiting novelty in their mode of action with specific emphasis on aminoacyl-tRNA synthetases as drug targets.

Probing the dual inhibitory mechanisms of novel thiophenecarboxamide derivatives against Mycobacterium tuberculosis PyrG and PanK: an insight from biomolecular modeling study

The growing occurrence of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mycobacterium tuberculosis (Mtb) strains underscores an urgent need for new antibiotics. The development of more bioactive antibiotics against drug-resistant organisms with a different mode of action could be a game-changer for the cure and eradication of tuberculosis (TB). Pantothenate Kinase (PanK) and CTP synthetase (PyrG) are both essential for RNA, DNA, and Lipids biosynthesis pathways. Given the extensive knowledge on these biosynthesis pathways inhibition of Mtb growth and survival, these enzymes present a fascinating opportunity for anti-mycobacterial drug discovery. Recently, it was experimentally established that the active metabolite 11426026 of compound 7947882 (a prodrug activated by EthA monooxygenase, 5-methyl-N-(4-nitrophenyl) thiophene-2-carboxamide) inhibits the activities of PyrG and PanK to indicate novel multitarget therapy aimed at discontinuing Mtb growth. However, the molecular mechanisms of their selective inhibition remain subtle. In this work, molecular dynamics simulations were employed to investigate the inhibitory mechanism as well as the selectivity impact of the active metabolite inhibitor of these enzymes. Computational modeling of the studied protein-ligand systems reveals that the active metabolite can potentially inhibit both PanK and PyrG, thereby creating a pathway as a double target approach in tuberculosis treatment.Communicated by Ramaswamy H. Sarma.

Computational studies of the properties and activities of selected trisubstituted benzimidazoles as potential antitubercular drugs inhibiting MTB-FtsZ polymerization

Trisubstituted benzimidazoles (trisbenz) are significantly active against nonreplicating Mycobacterium tuberculosis (MTB) by inhibiting the polymerization of Filamentous Temperature Sensitive Mutant Z (FtsZ), an essential bacteria cell division protein. In-depth in-silico study of 5 of the most active trisubstituted benzimidazoles; trisbenz 1, 2, 3, 4 and 5, giving insight into their properties, such as stability, bioavailability, interactions with residues at the binding site of MTB-FtsZ and their influence on structural dynamics of the protein have been conducted. This was achieved through the application of in-silico methods including density functional theory (DFT) calculations, ADME properties calculations, molecular docking and molecular dynamics simulations. A DFT approach was applied to predict reactivity properties of potent FtsZ inhibitors, and the results reveal the relative reactivity of these inhibitors as bioactive moieties. The estimated ADME properties predicted all 5 compounds to be bioavailable and suitable for oral administration. Molecular docking, binding free energy, RMSD, RMSF, and hydrogen bond analysis confirmed these 5 compounds as potent MTB-FtsZ inhibitors. Although analyses proved these compounds to be bioactive and potent MTB-FtsZ inhibitors, however, trisbenz 1 appeared to be the most active against this protein while trisbenz 5 was the least active. This study further confirms the experimental study while also giving insight on the compounds mechanism of action and presents their bioavailability properties.

Molecular Insights Into Di(2-Picolyl) Amine-Induced Cytotoxicity and Apoptosis in Human Kidney (HEK293) Cells

Di(2-picolyl) amine (DPA) is a pyridine derivative known to chelate metal ions and thus has potential anticancer properties; however, its effect on normal cells remains unchartered necessitating further research. This study, therefore, investigated the mechanistic effects of DPA-induced cytotoxicity and apoptosis in the HEK293 cell line. Methods required that an half the maximum inhibition concentration (IC₅₀) was derived using the 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide (MTT) assay. Analyses aimed to assess oxidative stress, membrane damage, and DNA fragmentation by means of biochemical assays were performed. Luminometry analysis was carried out to understand the mechanism of apoptosis induction by determining the levels of adenosine triphosphate (ATP) and the activities of caspase-8, -9, and -3/7. Western blotting was used to ascertain the expression of apoptotic and stress-related proteins. An IC₅₀ of 1,079 µM DPA was obtained. Antioxidant effect correlated with a minimum increase in reactive oxygen species induced lipid peroxidation. The increase in initiator caspase-8 and -9 and executioner caspase-3/7 activities by DPA-induced apoptosis albeit prompting a decline in the levels of ATP. Furthermore, DPA brought about the following consequences on HEK293 cells: markedly elevated tail lengths of the comets, poly (ADP-ribose) polymerase 1 cleavage, and apoptotic body formation observed in the late stages. The cytotoxic effects of DPA in HEK293 cells may be mediated by induction of apoptosis via the caspase-dependent mechanism.

Understanding the Hsp90 N-terminal Dynamics: Structural and Molecular Insights into the Therapeutic Activities of Anticancer Inhibitors Radicicol (RD) and Radicicol Derivative (NVP-YUA922)

Heat shock protein 90 (Hsp90) is a crucial component in carcinogenesis and serves as a molecular chaperone that facilitates protein maturation whilst protecting cells against temperature-induced stress. The function of Hsp90 is highly dependent on adenosine triphosphate (ATP) binding to the N-terminal domain of the protein. Thus, inhibition through displacement of ATP by means of competitive binding with a suitable organic molecule is considered an attractive topic in cancer research. Radicicol (RD) and its derivative, resorcinylic isoxazole amine NVP-AUY922 (NVP), have shown promising pharmacodynamics against Hsp90 activity. To date, the underlying binding mechanism of RD and NVP has not yet been investigated. In this study, we provide a comprehensive understanding of the binding mechanism of RD and NVP, from an atomistic perspective. Density functional theory (DFT) calculations enabled the analyses of the compounds' electronic properties and results obtained proved to be significant in which NVP was predicted to be more favorable with solvation free energy value of -23.3 kcal/mol and highest stability energy of 75.5 kcal/mol for a major atomic delocalization. Molecular dynamic (MD) analysis revealed NVP bound to Hsp90 (NT-NVP) is more stable in comparison to RD (NT-RD). The Hsp90 protein exhibited a greater binding affinity for NT-NVP (-49.4 ± 3.9 kcal/mol) relative to NT-RD (-28.9 ± 4.5 kcal/mol). The key residues influential in this interaction are Gly 97, Asp 93 and Thr 184. These findings provide valuable insights into the Hsp90 dynamics and will serve as a guide for the design of potent novel inhibitors for cancer treatment.

Exploring the impact of H5N1 neuraminidase (H274Y) mutation on Peramivir: a bio-computational study from a molecular perspective

Communicated by Ramaswamy H. Sarma.

Induced Mutation Proves a Potential Target for TB Therapy: A Molecular Dynamics Study on LprG

Molecular dynamics (MD) simulations of wild-type and V91W mutant Mycobacterium tuberculosis-LprG (Mtb-LprG) were performed with the goal to provide a comprehensive understanding of the Mtb-LprG as a potential antimycobacterial target. A long-range MD simulations and post-MD analyzes led us to various results that plainly explained the impact of V91W mutation on Mtb-LprG. Herein, the results revealed that the wild-type is less stable compared to V91W mutant. This was further supported by root mean square fluctuation, where the V91W mutant showed a higher degree of flexibility compared to the wild-type. Dynamic cross-correlation analysis revealed that induced mutation leads to higher residual flexibility in the mutant structure as compared to the wild-type structure thus resulting in the existence of negatively correlated motions. The difference in principal component analysis scatter plot across the first two normal modes suggests a greater mobility of the V91W mutant conformation compared to the wild-type. Thermodynamic calculations revealed that the van der Waals (E_vdw) forces contribute the most towards binding free energy in a case of the V91W mutant as compared to the wild-type LprG complex. In addition, the residue interaction networks revealed more of E_vdw interaction existence among residues in case of the V91W mutant. This study supports the Mtb-LprG as a potential antimycobacterial target and also serves as a cornerstone to identifying new potential targets that have no inhibitors.

The impact of Thr91 mutation on c-Src resistance to UM-164: molecular dynamics study revealed a new opportunity for drug design

The emergence of a drug resistant non-receptor tyrosine kinase (c-Src) in triple-negative breast cancer (TNBC) remains a prime concern in relation to the burden of TNBC among people living with breast cancer and drug development. Thr91 mutation was found to induce a complete loss of protein conformation required for drug fitness. Herein, we provide the first account of the molecular impact of the Thr91 mutation on c-Src resistance to experimental drug UM-164 using various computational approaches, namely molecular dynamics simulation, principal component analysis (PCA), dynamic cross-correlation matrices (DCCM) analysis, hydrogen bond occupancy, thermodynamics calculation, ligand-residue interaction and residue interaction networks (RINs). Findings from this study revealed that Thr91 mutation leads to a steric conflict between UM-164 and the side chain of methionine (Met91); this mutation distorts the UM-164 optimum orientation on the conformational space of mutant c-Src compared to the wild-type; decreases hydrogen bond formation between the residues in the mutant protein structure; decreases the UM-164 binding energy in the mutant by -13.416 kcal mol^-1; reduces the residue correlation in the mutant protein structure; induces a change in the overall protein structure conformation from an inactive to active conformation; and distorts the ligand atomic interaction network and the residue interaction network. This report provides important insights that will assist in the further design of novel dual kinase inhibitors to minimise the chances of drug resistance in triple negative breast cancer.

Re-emergence of an orphan therapeutic target for the treatment of resistant prostate cancer - a thorough conformational and binding analysis for ROR-γ protein

Recent studies have linked a deadly form of prostate cancer known as metastatic castration-resistant prostate cancer to retinoic acid-related orphan-receptor gamma (ROR-γ). Most of these studies continued to place ROR-γ as orphan because of unidentifiable inhibitor. Recently identified inhibitors of ROR-γ and their therapeutic potential were evaluated, among which inhibitor XY018 was the potent. However, molecular understanding of the conformational features of XY018-ROR-γ complex is still elusive. Herein, molecular dynamics simulations were conducted on HC9-ROR-γ and XY018-ROR-γ complexes to understand their conformational features at molecular level and the influence of XY018 binding on the dynamics of ROR-γ with the aid of post-dynamic analytical tools. These include; principal component analysis, radius of gyration, binding free energy calculation (MM/GBSA), per-residue fluctuation and hydrogen bond occupancy. Findings from this study revealed that (1) hydrophobic packing contributes significantly to binding free energy, (2) Ile136 and Leu60 exhibited high hydrogen-bond occupancy in XY018-ROR-γ and HC9-ROR-γ, respectively, (3) XY018-ROR-γ displayed a relatively high loop region residue fluctuation compared to HC9-ROR-γ, (4) electrostatic interactions are a potential binding force in XY018-ROR-γ complex compared to HC9-ROR-γ, (5) XY018-ROR-γ assumes a rigid conformation which is highlighted by a decrease in residual fluctuation, (6) XY018 could potentially induce pseudoporphyria, nephritis and interstitial nephritis but potentially safe in renal failure. This study could serve as a base line for the design of new potential ROR-γ inhibitors.

Sliding Clamp of DNA Polymerase III as a Drug Target for TB Therapy: Comprehensive Conformational and Binding Analysis from Molecular Dynamic Simulations

Tuberculosis, caused by Mycobacterium tuberculosis, is one of the most common causes of death in the world. Mycobacterium tuberculosis -sliding clamp is a protein essential for many important DNA transactions including replication and DNA repair proteins, thus, a potential drug target for tuberculosis. Further investigation is needed in understanding DNA polymerase sliding clamp structure, especially from a computational perspective. In this study, we employ a wide-range of comparative molecular dynamic analyses on two systems: Mycobacterium tuberculosis - sliding clamp enzyme in its apo and bound form. The results reported in this study shows apo conformation to be less stable, as compared to bound conformation with an average radius of gyration of 25.812 and 25.459 Å, respectively. This was further supported by root mean square fluctuation, where an apo enzyme showed a higher degree of flexibility. However, the presence of the ligand lowers radius of gyration and root mean square fluctuation and also leads to an existence of negative correlated motions. Principal component analysis further justifies the same findings, whereby the apo enzyme exhibits a higher fluctuation compared to the bound complex. In addition, a stable 3₁₀ helix located at the binding site appears to be unstable in the presence of the ligand. Hence, it is possible that the binding of the ligand may have caused a rearrangement of the structure, leading to a change in the unwinding of 3₁₀ helix. Findings reported in this study further enhance the understanding of Mycobacterium tuberculosis -DnaN and also give a lead to the development of potent tuberculosis drugs.

Single H5N1 influenza A neuraminidase mutation develops resistance to oseltamivir due to distorted conformational and drug binding landscape: multiple molecular dynamics analyses

Clinical studies showed that a single mutation, I117V, develops severe resistance to oseltamivir, the first orally active influenza A neuraminidase inhibitor, in highly pathogenic H5N1 influenza A viruses.

Dynamics of the thumb-finger regions in a GH11 xylanase Bacillus circulans: comparison between the Michaelis and covalent intermediate

The structure and dynamics of B. circulans β-1,4-xylanase (BCX) were comparatively studied utilizing molecular dynamics.