In silico study directed towards identification of the key structural features of GyrB inhibitors targeting MTB DNA gyrase: HQSAR, CoMSIA and molecular dynamics simulations

Revealing the Interaction Mechanism between Mycobacterium tuberculosis GyrB and Novobiocin, SPR719 through Binding Thermodynamics and Dissociation Kinetics Analysis

With the rapid emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), various levels of resistance against existing anti-tuberculosis (TB) drugs have developed. Consequently, the identification of new anti-TB targets and drugs is critically urgent. DNA gyrase subunit B (GyrB) has been identified as a potential anti-TB target, with novobiocin and SPR719 proposed as inhibitors targeting GyrB. Therefore, elucidating the molecular interactions between GyrB and its inhibitors is crucial for the discovery and design of efficient GyrB inhibitors for combating multidrug-resistant TB. In this study, we revealed the detailed binding mechanisms and dissociation processes of the representative inhibitors, novobiocin and SPR719, with GyrB using classical molecular dynamics (MD) simulations, tau-random acceleration molecular dynamics (τ-RAMD) simulations, and steered molecular dynamics (SMD) simulations. Our simulation results demonstrate that both electrostatic and van der Waals interactions contribute favorably to the inhibitors' binding to GyrB, with Asn52, Asp79, Arg82, Lys108, Tyr114, and Arg141 being key residues for the inhibitors' attachment to GyrB. The τ-RAMD simulations indicate that the inhibitors primarily dissociate from the ATP channel. The SMD simulation results reveal that both inhibitors follow a similar dissociation mechanism, requiring the overcoming of hydrophobic interactions and hydrogen bonding interactions formed with the ATP active site. The binding and dissociation mechanisms of GyrB with inhibitors novobiocin and SPR719 obtained in our work will provide new insights for the development of promising GyrB inhibitors.

Molecular docking-based interaction studies on imidazo[1,2-a] pyridine ethers and squaramides as anti-tubercular agents

Development of new anti-tubercular agents is required in the wake of resistance to the existing and newly approved drugs through novel-validated targets like ATP synthase, etc. The major limitation of poor correlation between docking scores and biological activity by SBDD was overcome by a novel approach of quantitatively correlating the interactions of different amino acid residues present in the target protein structure with the activity. This approach well predicted the ATP synthase inhibitory activity of imidazo[1,2-a] pyridine ethers and squaramides (r = 0.84) in terms of Glu65b interactions. Hence, the models were developed on combined (r = 0.78), and training (r = 0.82) sets of 52, and 27 molecules, respectively. The training set model well predicted the diverse dataset (r = 0.84), test set (r = 0.755), and, external dataset (r_ext = 0.76). This model predicted three compounds from a focused library generated by incorporating the essential features of the ATP synthase inhibition with the pIC₅₀ values in the range of 0.0508-0.1494 µM. Molecular dynamics simulation studies ascertain the stability of the protein structure and the docked poses of the ligands. The developed model(s) may be useful in the identification and optimization of novel compounds against TB.

Molecular docking-based interactions in QSAR studies on Mycobacterium tuberculosis ATP synthase inhibitors

Tuberculosis (TB) is a global threat with a large burden across the continents in terms of mortality, morbidity, and financial losses. The disease has evolved into multi-drug-resistant (MDR-TB) and extensively drug-resistant (XDR-TB) tuberculosis owing to numerous factors ranging from patients' non-compliance to demographical implications. There have been very few new drugs for resistant TB. Resistance has already been reported even for the newly introduced drug bedaquiline. An attempt has been made to integrate both structure-based and QSAR drug design techniques (QSAR-SBDD) for the identification of novel leads. The docking scores normally do not correlate with the activity. Hence, the docking results have been analysed in terms of the number of interactions rather than docking scores. The parameters derived from interactions have been used in developing the QSAR models. The best model shows a good correlation (r = 0.908) between the activity and interaction parameter 'C' describing the sum of all the interactions with each amino acid residue. This model also predicts external dataset with a good correlation (r_ext = 0.851) and can be used for the identification of novel chemical entities (NCEs) and repurposed drugs for TB therapeutics.

Ligand-based design of anticancer MMP2 inhibitors: a review

MMP2, a Zn²⁺-dependent metalloproteinase, is related to cancer and angiogenesis. Inhibition of this enzyme might result in a potential antimetastatic drug to leverage the anticancer drug armory. In silico or computer-aided ligand-based drug design is a method of rational drug design that takes multiple chemometrics (i.e., multi-quantitative structure-activity relationship methods) into account for virtually selecting or developing a series of probable selective MMP2 inhibitors. Though existing matrix metalloproteinase inhibitors have shown plausible pan-matrix metalloproteinase (MMP) activity, they have resulted in various adverse effects leading to their being rescinded in later phases of clinical trials. Therefore a review of the ligand-based designing methods of MMP2 inhibitors would result in an explicit route map toward successfully designing and synthesizing novel and selective MMP2 inhibitors.

Insights into interaction mechanism of inhibitors E3T, E3H and E3B with CREB binding protein by using molecular dynamics simulations and MM-GBSA calculations

CREB binding protein (CBP) and its paralog E1A binding protein (p300) are related to the development of inflammatory diseases, cancers and other diseases, and have been potential targets for the treatment of human diseases. In this work, interaction mechanism of three small molecules E3T, E3H, and E3B with CBP was investigated by employing molecular dynamics (MD) simulations, principal component analysis (PCA), and molecular mechanics/generalized born surface area (MM-GBSA) method. The results indicate that inhibitor bindings cause the changes of movement modes and structural flexibility of CBP, and van der Waals interactions mostly drive associations of inhibitors with CBP. In the meantime, the results based on inhibitor-residue interactions not only show that eight residues of CBP can strongly interact with E3T, E3H and E3B but also verify that the CH-CH, CH-π, and π-π interactions are responsible for vital contributions in associations of E3T, E3H and E3B with CBP. In addition, the H-O radial distribution functions (RDFs) were computed to assess the stability of hydrogen bonding interactions between inhibitors and CBP, and the obtained information identifies several key hydrogen bonds playing key roles in bindings of E3T, E3H and E3B to CBP. Potential new inhibitors have been proposed.

Structure-Based Screening of DNA GyraseB Inhibitors for Therapeutic Applications in Tuberculosis: a Pharmacoinformatics Study

Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis (MTB) and considered as serious public health concern worldwide which kills approximately five thousand people every day. Therefore, TB drug development efforts are in gigantic need for identification of new potential chemical agents to eradicate TB from the society. The bacterial DNA gyrase B (GyrB) protein as an experimentally widely accepted effective drug target for the development of TB chemotherapeutics. In the present study, advanced pharmacoinformatics approaches were used to screen the Mcule database against the GyrB protein. Based on a number of chemometric parameters, five molecules were found to be crucial to inhibit the GyrB. A number of molecular binding interactions between the proposed inhibitors and important active site residues of GyrB were observed. The predicted drug-likeness properties of all molecules were indicated that compounds possess characteristics to be the drug-like candidates. The dynamic nature of each molecule was explored through the molecular dynamics (MD) simulation study. Various analyzing parameters from MD simulation trajectory have suggested rationality of the molecules to be potential GyrB inhibitor. Moreover, the binding free energy was calculated from the entire MD simulation trajectories highlighted greater binding free energy values for all newly identified compounds also substantiated the strong binding affection towards the GyrB in comparison to the novobiocin. Therefore, the proposed molecules might be considered as potential anti-TB chemical agents for future drug discovery purposes subjected to experimental validation. Graphical Abstract.

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

Advancement of GyrB Inhibitors for Treatment of Infections Caused by Mycobacterium tuberculosis and Non-tuberculous Mycobacteria

The prospect of ever increasing antibiotic resistance eroding currently available treatment options for bacterial infections underscores the need to continue to identify new antibiotics, preferably those that act on novel targets or with novel mechanisms of action. Bacterial gyrase B subunit (GyrB), an essential component of bacterial gyrase required for successful DNA replication, represents such a target. We describe recent examples of GyrB inhibitors and point out their potential utility for treatment of mycobacterial diseases caused by Mycobacterium tuberculosis (TB) and non-tuberculous mycobacteria (NTM). Current therapeutic options for these diseases are often suboptimal due to resistance to current standard of care antibiotics. A future GyrB inhibitor-based antibiotic could offer a new and effective addition to the armamentarium for treatment of mycobacterial diseases and possibly for infections caused by other bacterial pathogens. One GyrB inhibitor, SPR720, has recently completed a first-in-human clinical trial and is in clinical development for the treatment of NTM and TB infections.