Molecular Dynamics Assisted Mechanistic Study of Isoniazid-Resistance against Mycobacterium tuberculosis InhA

Exploring natural products library to identify potential inhibitors targeting isoniazid-resistant tuberculosis

Mycobacterium tuberculosis (MTB) causing tuberculosis (TB) infection is a leading source of illness and death in developing nations, and the emergence of drug-resistant TB remains a significant global threat and a challenge in treating the disease. Mutations in the inhA and katG genes are connected to the principal molecular mechanism of isoniazid (INH) resistance, and continuous treatment of INH for more than a decade led to the evolution of INH resistant-TB (inhR-TB). Structure-based drug discovery approaches on traditional natural compounds are the contemporary source to identify significant lead molecules. This work focuses on discovering effective small compounds from natural compound libraries and applying pharmacophore-based virtual screening to filter out the molecules. The best-identified hit complexes were used for molecular dynamics simulations (MDS) to observe their stability and compactness. A three-dimensional e-pharmacophore hypothesis and screening generated 62 hits based on phase fitness scores from the pharmacophore-based virtual screening. Molecular docking experiments in Maestro's GLIDE module indicated that ZINC000002383126 and ASN22022 may be potential inhibitors of inhA and katG (native, inhA mutants S94A, Y158A, Y158F and Y158S and D137S, Y229F, S315T, W321F, and R418L mutants of katG). In addition, MDS analysis indicated that the native and mutant docked complexes of inhA and katG had good stability and remained compact in the binding pocket of the targets. In vitro studies can further validate the compounds that can act as INH competitive inhibitors.Communicated by Ramaswamy H. Sarma.

Structural dynamics insights into the M306L, M306V, and D1024N mutations in Mycobacterium tuberculosis inducing resistance to ethambutol

Resistance to anti-tuberculosis drugs, especially ethambutol (EMB), has been widely reported worldwide. EMB resistance is caused by mutations in the embB gene, which encodes the arabinosyl transferase enzyme. This study aimed to detect mutations in the embB gene of Mycobacterium tuberculosis from Papua and to evaluate their impact on the effectiveness of EMB. We analyzed 20 samples of M. tuberculosis culture that had undergone whole-genome sequencing, of which 19 samples were of sufficient quality for further bioinformatics analysis. Mutation analysis was performed using TBProfiler, which identified M306L, M306V, D1024N, and E378A mutations. In sample TB035, the M306L mutation was present along with E378A. The binding affinity of EMB to arabinosyl transferase was calculated using AutoDock Vina. The molecular docking results revealed that all mutants demonstrated an increased binding affinity to EMB compared to the native protein (-0.948 kcal/mol). The presence of the M306L mutation, when coexisting with E378A, resulted in a slight increase in binding affinity compared to the M306L mutation alone. The molecular dynamics simulation results indicated that the M306L, M306L + E378A, M306V, and E378A mutants decreased protein stability. Conversely, the D1024N mutant exhibited stability comparable to the native protein. In conclusion, this study suggests that the M306L, M306L + E378A, M306V, and E378A mutations may contribute to EMB resistance, while the D1024N mutation may be consistent with continued susceptibility to EMB.

Facing Antitubercular Resistance: Identification of Potential Direct Inhibitors Targeting InhA Enzyme and Generation of 3D-pharmacophore Model by in silico Approach

Purpose

The enoyl-acyl carrier protein reductase (InhA) is one of the important key enzymes employed in mycolic acids biosynthesis pathway and an important component of mycobacterial cell walls. This enzyme has also been identified as major target of isoniazid drug, except that isoniazid needs to be activated first by the catalase peroxidase (KatG) protein to form the isonicotinoyl-NAD (INH-NAD) adduct that inhibits the action of InhA enzyme. However, this activation becomes more difficult and unreachable with the problem of mutation-related resistance caused mainly by acquired mutations in KatG and InhA protein. Our main interest in this study is to identify direct InhA inhibitors using computer-aided drug design.

Methods

Computer-aided drug design was used to solve this problem by applying three different approaches including mutation impact modelling, virtual screening and 3D-pharmacophore search.

Results

A total of 15 mutations were collected from the literature, then a 3D model was generated for each of them and their impact was predicted. Of the 15 mutations, 10 were found to be deleterious and have a direct effect on flexibility, stability and SASA of the protein. In virtual screening, from 1,000 similar INH-NAD analogues obtained by the similarity search method, 823 compounds passed toxicity filter and drug likeness rules, which were then docked to the wild-type of InhA protein. Subsequently, 34 compounds with binding energy score better than that of INH-NAD were selected and docked against the 10 generated mutated models of InhA. Only three leads showed a lower binding affinity better than the reference. The 3D-pharmacophore model approach was used to identify the common features between those three compounds by generating a pharmacophoric map.

Conclusion

The result of this study may pave the way to develop more potent mutant-specific inhibitors to overcome this resistance.

Applications of Molecular Simulation in the Discovery of Antituberculosis Drugs: A Review

After decades of efforts, tuberculosis has been well controlled in most places. The existing drugs are no longer sufficient for the treatment of drug-resistant Mycobacterium tuberculosis due to significant toxicity and selective pressure, especially for XDR-TB. In order to accelerate the development of high-efficiency, low-toxic antituberculosis drugs, it is particularly important to use Computer Aided Drug Design (CADD) for rational drug design. Here, we systematically reviewed the specific role of molecular simulation in the discovery of new antituberculosis drugs. The purpose of this review is to overview current applications of molecular simulation methods in the discovery of antituberculosis drugs. Furthermore, the unique advantages of molecular simulation was discussed in revealing the mechanism of drug resistance. The comprehensive use of different molecular simulation methods will help reveal the mechanism of drug resistance and improve the efficiency of rational drug design. With the help of molecular simulation methods such as QM/MM method, the mechanisms of biochemical reactions catalyzed by enzymes at atomic level in Mycobacterium tuberculosis has been deeply analyzed. QSAR and virtual screening both accelerate the development of highefficiency, low-toxic potential antituberculosis drugs. Improving the accuracy of existing algorithms and developing more efficient new methods for CADD will always be a hot topic in the future. It is of great value to utilize molecular dynamics simulation to investigate complex systems that cannot be studied in experiments, especially for drug resistance of Mycobacterium tuberculosis.

Analysis of mutations leading to para-aminosalicylic acid resistance in Mycobacterium tuberculosis

Thymidylate synthase A (ThyA) is the key enzyme involved in the folate pathway in Mycobacterium tuberculosis. Mutation of key residues of ThyA enzyme which are involved in interaction with substrate 2′-deoxyuridine-5′-monophosphate (dUMP), cofactor 5,10-methylenetetrahydrofolate (MTHF), and catalytic site have caused para-aminosalicylic acid (PAS) resistance in TB patients. Focusing on R127L, L143P, C146R, L172P, A182P, and V261G mutations, including wild-type, we performed long molecular dynamics (MD) simulations in explicit solvent to investigate the molecular principles underlying PAS resistance due to missense mutations. We found that these mutations lead to (i) extensive changes in the dUMP and MTHF binding sites, (ii) weak interaction of ThyA enzyme with dUMP and MTHF by inducing conformational changes in the structure, (iii) loss of the hydrogen bond and other atomic interactions and (iv) enhanced movement of protein atoms indicated by principal component analysis (PCA). In this study, MD simulations framework has provided considerable insight into mutation induced conformational changes in the ThyA enzyme of Mycobacterium.

An integrated computational approach of molecular dynamics simulations, receptor binding studies and pharmacophore mapping analysis in search of potent inhibitors against tuberculosis

Tuberculosis is an infectious chronic disease caused by obligate pathogen Mycobacterium tuberculosis that affects millions of people worldwide. Although many first and second line drugs are available for its treatment, but their irrational use has adversely lead to the emerging cases of multiple drug resistant and extensively drug-resistant tuberculosis. Therefore, there is an intense need to develop novel potent analogues for its treatment. This has prompted us to develop potent analogues against TB. The Mycobacterium tuberculosis genome provides us with number of validated targets to combat against TB. Study of Mtb genome disclosed six epoxide hydrolases (A to F) which convert harmful epoxide into diols and act as a potential drug target for rational drug design. Our current strategy is to develop such analogues which inhibits epoxide hydrolase enzyme present in Mtb genome. To achieve this, we adopted an integrated computational approach involving QSAR, pharmacophore mapping, molecular docking and molecular dynamics simulation studies. The approach envisaged vital information about the role of molecular descriptors, essential pharmacophoric features and binding energy for compounds to bind into the active site of epoxide hydrolase. Molecular docking analysis revealed that analogues exhibited significant binding to Mtb epoxide hydrolase. Further, three docked complexes 2s, 37s and 15s with high, moderate and low docking scores respectively were selected for molecular dynamics simulation studies. RMSD analysis revealed that all complexes are stable with average RMSD below 2 Å throughout the 10 ns simulations. The B-factor analysis showed that the active site residues of epoxide hydrolase are flexible enough to interact with inhibitor. Moreover, to confirm the binding of these urea derivatives, MM-GBSA binding energy analysis were performed. The calculations showed that 37s has more binding affinity (ΔGtotal = -52.24 kcal/mol) towards epoxide hydrolase compared to 2s (ΔGtotal = -51.70 kcal/mol) and 15s (ΔGtotal = -49.97 kcal/mol). The structural features inferred in our study may provide the future directions to the scientists towards the discovery of new chemical entity exhibiting anti-TB property.

The two mutations of actin-myosin interface and their effect on the dynamics, structures, and functions of skeletal muscle actin

Congenital myopathy is a broad category of muscular diseases with symptoms appearing at the time of birth. One type of congenital myopathy is Congenital Fiber Type Disproportion (CFTD), a severely debilitating disease. The G48D and G48C mutations in the D-loop and the actin-myosin interface are the two causes of CFTD. These mutations have been shown to significantly affect the structure and function of muscle fibers. To the author's knowledge, the effects of these mutations have not yet been studied. In this work, the power stroke structure of the head domain of myosin and the wild and mutated types of actin were modeled. Then, a MD simulation was run for the modeled structures to study the effects of these mutations on the structure, function, and molecular dynamics of actin. The wild and mutated actins docked with myosin showed differences in hydrogen bonding patterns, free binding energies, and hydrogen bond occupation frequencies. The G48D and G48C mutations significantly impacted the conformation of D-loops because of their larger size compared to Glycine and their ability to interfere with the polarity or hydrophobicity of this neutralized and hydrophobic loop. Therefore, the mutated loops were unable to fit properly into the hydrophobic groove of the adjacent G-actin. The abnormal structure of D-loops seems to result in the abnormal assembly of F-actins, giving rise to the symptoms of CFTD. It was also noted that G48C and G48D did not form hydrogen bonds with myosin in the residue 48 location. Nevertheless, in this case, muscles are unable to contract properly due to muscle atrophy.

Cell wall permeability assisted virtual screening to identify potential direct InhA inhibitors of Mycobacterium tuberculosis and their biological evaluation

The arising cases of isoniazid-resistance have motivated research interests toward new class of molecules known as direct InhA inhibitors. Here, a combine approach of shape-based pharmacophore and descriptor-based 2D QSAR was used to identify the potential direct InhA inhibitors. The approach is duly assisted with in vitro testing and molecular dynamics simulations. A combination of empirical parameters was derived to use as a filter for cell wall permeability while 2D QSAR was used as another filter to predict the biological activity. Both filters were applied to prioritize the molecules for biological evaluation against anti-TB activity. It led to 6 potential molecules which showed > 90% inhibition of H37Rv strain of Mycobacterium tuberculosis in BACTEC assay. Further, MMGBSA binding free energy of identified molecules was compared with available highly potent molecule, 5-hexyl-2-(2-methylphenoxy) phenol (IC₅₀ = 5nM) using molecular dynamics simulations. It showed two molecules with comparatively higher affinity toward InhA as compared to potent molecule. It indicated the candidature of identified molecules to be further considered in anti-TB drug development pipeline.

S. S, Kulkarni M. Design, synthesis and evaluation of diphenyl ether analogues as antitubercular agents

We herein report the investigation of new diphenyl ethers asMycobacterium tuberculosisenoyl-acyl carrier protein reductase (InhA) inhibitors by structure-based drug design approach.