Structure-Based Virtual Screening of Benzaldehyde Thiosemicarbazone Derivatives against DNA Gyrase B of Mycobacterium tuberculosis

Emergence of antibiotic-resistant Mycobacterium tuberculosis (M. tuberculosis) restricts the availability of drugs for the treatment of tuberculosis, which leads to the increased morbidity and mortality of the disease worldwide. There are many intrinsic and extrinsic factors that have been reported for the resistance mechanism. To overcome such mechanisms, chemically synthesized benzaldehyde thiosemicarbazone derivatives were screened against M. tuberculosis to find potential inhibitor for tuberculosis. Such filtering process resulted in compound 13, compound 21, and compound 20 as the best binding energy compounds against DNA gyrase B, an important protein in the replication process. The ADMET prediction has shown the oral bioavailability of the novel compounds.


Introduction
Tuberculosis (TB) is a potentially serious communicable disease caused by the bacillus Mycobacterium tuberculosis (M. tuberculosis) [1]. TB was found with increased mortality, and it spreads from person to person through tiny droplets released into the air via coughs and sneezes [2]. According to the World Health Organization (WHO) report, worldwide, around 10 million people were affected with TB [3]. Here, 56% were men, 32% were women, and children accounted for 12%. Also, 1.2 million deaths were reported in HIVnegative people and 208000 deaths were reported among HIV-positive people. e largest number of new TB cases occurred in the Southeast Asian region with 44% of new cases, followed by the African region with 25% of new cases and Western Pacific with 18%. Eight countries accounted for two-thirds of the new TB cases including India (26%), Indonesia (8.5%), China (8.4%), Philippines (6.0%), Pakistan (5.7%), Nigeria (4.4%), Bangladesh (3.6%), and South Africa (3. 6) [3]. M. tuberculosis does not have any particular virulence factor, but it perseveres long-term in the human body without causing any significant damage and transmission, if not, the immune system of the host is compromised. It secrets effector proteins to complicate the immune system, thereby, stimulating its intracellular survival in granulomas throughout the latency period of infection [4]. M. tuberculosis can develop in different conditions such as pulmonary and extrapulmonary TB (pleural, lymphadenitis, skeletal, gastrointestinal, and ocular) [5]. Once M. tuberculosis entered into the lungs, alveolar macrophages engulf the organism, and it was captured by phagosomes. Finally, it was delivered to lysosomes and get degraded. However, in most of the instances, to persist in human alveolar macrophages, the organism inhibits the acidification and phagosomes maturation. Also, to escape from the immune system recognition and hypoxic condition, the organism remains in a "quiescent status" (in nongrowing state but metabolically active) in 90% of the diseased individuals [6]. e innate (antimycobacterial elements, IFNc and TNFα) as well as adaptive immunity (T cells, 17, CD4 + , and CD8 + ) controls M. tuberculosis activity in the latent phase of the infection. Hence, it is known as the conditional pathogenic bacterium [7]. If the host immune system is compromised, the organism gets activated and initiate replication. is promotes the diseased macrophage necrosis and, thereby, discharges the intracellular bacteria, which infect novel cells and invade other tissues [8]. Recent reports suggest that TB is also related with many other human complications, namely, autoimmune diseases (sarcoidosis: M. tuberculosis activates toll-like receptors, thereby, promoting pulmonary sarcoidosis [9]; systemic lupus erythematosus: immunosuppressive therapy and several immune abnormalities cause reactivation and diffusion of TB [10]), metabolic syndromes (diabetes mellitus: promotes the proliferation of mycobacterium [11]; atherosclerosis: mycobacterium rushes the development of atherosclerosis [12]), and pulmonary complications (pneumonia: TB infection rises the risk of secondary bacterial infection in children [13]; chronic obstructive pulmonary disease: pulmonary TB may alter the lung architecture [14]; lung cancer: TB is one of the risk factors for the development of lung cancer [15,16]).
Treatment options available for drug-susceptible TB involve first-line drugs (isoniazid, rifampicin, ethambutol, and pyrazinamide) for six months. e indiscriminate use of the antibiotics leads to the development of resistance. e WHO classifies the resistance as multidrug resistance (MDR), extensive drug resistance (XDR), and total drug resistance (TDR). e rise of resistance to the first-line drugs termed as multidrug resistance (MDR) and resistance to the fluoroquinolones, one of the second-line drugs (capreomycin, kanamycin, and amikacin) tends to the occurrence of XDR. Both resistance to the first and second-line drugs leads to the development of TDR [3]. e U.S Food and Drug Administration (FDA) has developed bedaquiline, a novel drug against MDR M. tuberculosis. Over the 40 years, bedaquiline have been in use to treat TB. Unfortunately, resistance to bedaquiline antibiotic has been reported in recent years [17,18]. e extrinsic factors associated with antibiotic resistance are social elements of TB in inhabitants and the eminence of prevention and control services of TB. Gygli et al., in 2017, explained the mechanism of intrinsic drug resistance in M. tuberculosis, such as decreased cell wall permeation, secretion of drug inactivating enzymes, mutation in the cell wall efflux system, bacterial drug target modification, and overexpression of drug targets. Here, isoniazid prodrug is activated by catalase or peroxidase (gene: katG), and it functions through enoyl acyl carrier protein reductase (gene: inhA). Rifampicin acts against M. tuberculosis through binding to the RNA polymerase β subunit (gene: rpoB). Ethambutol inhibits arabinosyl transferase (gene: embB) involved in the biosynthesis of cell wall arabinogalactan. Pyrazinamide, a nicotinamide analog, requires pyrazinamidase or nicotinamidase (gene: pncA) to get converted into pyrazinoic acid, an active form. For second-line drugs, fluoroquinolones inhibit the topoisomerase II (DNA gyrase), an essential enzyme involved in the replication process (genes: gyrA and gyrB). e role of kanamycin and amikacin modifies the range of 16S rRNA (gene: rrs) and, thereby, inhibits the protein synthesis. e capreomycin is found to inhibit the translation process in mycobacterium. e target gene tlyA participates in the rRNA ribose specific 2′-O-methylation. Here, the development of resistance to these first-line and second-line drugs is due to the mutations in the drug-targeting genes [19][20][21][22][23]. Such situation necessitates the design and development of novel drug with high antitubercular activity. In the recent study, Volynets et al., in 2019, have experimentally proved the activity of benzaldehyde thiosemicarbazone against M. tuberculosis. Hence, in the present study, we have retrieved the synthesized benzaldehyde thiosemicarbazone derivatives from the literature and screened against one of the main target DNA gyrases to find a potential inhibitor for M. tuberculosis [24] through molecular modelling methods. Our study may give an idea to the researchers who are designing drug against M. tuberculosis in the molecular level.

Protein Structure Preparation.
e 3-dimensional (3D) structure of the target protein DNA gyrase was obtained from the protein data bank (PDB: 6GAU) [25,26]. e PDB structure is a homodimer, which has two chains, chain A (DNA gyrase subunit B) and B (DNA gyrase subunit A) ( Figure 1). For our virtual screening and molecular docking studies, we have utilized chain A. e sequence length of the subunit is 1179 amino acids. e 3D structure was crystalized by the X-ray diffraction method. e nonamino acid structures cocrystallized with DNA gyrase, phosphoaminophosphonic acid-adenylate ester, magnesium ions, and water molecules were detached from the structure.

Active Site Prediction.
e active site of the DNA gyrase subunit B was predicted through the DoGSiteScorer [27]. Finding the low molecular weight ligand molecule by the target protein is the source for maintenance of the biological system. Here, active site is crucial for the function of an enzyme [28]. Usually, the active site prediction helps to ensure the protein function, druggability, and family classification. Such prediction is carried out by the DoGSite server, a structure-based prediction method, which functions on the basis of the difference of Gaussian (DoG) method. Most of the computational techniques available are geometry-based, energy-based, evolutionary-based, and combine method predictions. All of these approaches have disadvantages [29]. In contrast, DoGSite uses the pattern recognition method to predict the active sites. is locates the active site regions by sieving the grid of the protein using the DoG filter, and it finds the spherically designed structures called DoG cores in the grid. Finally, these detected DoG cores are gathered to form pockets, where the ligand molecules can fit firmly.

Ligand Preparation.
e synthesized 30 benzaldehyde thiosemicarbazone ( Figure 2) derivatives were obtained from the literature (9). e structures were drawn with the help of ChemSketch software and optimized (10) ( Table 1).

Virtual Screening.
A total of 30 benzaldehyde thiosemicarbazone derivatives were screened against gyrB using the python prescription virtual screening tool (PyRx) [30], a structure-based virtual screening process, which usually screens compound libraries against protein targets. Using the OpenBabel tool combined with the PyRx server, the ligand molecules were added and subjected to energy minimization with the help of the universal force field (UFF) by the conjugate gradient algorithm. Both protein and ligand structures were converted into PDBQT format. e docking procedure was carried out through the AutoDock Vina of PyRx Tool. From the 30 compounds, best binding energy was selected and passed to the next procedure, docking studies to get final best compounds.

Molecular Docking Studies.
Prediction of binding orientation of ligand with the protein is carried out through the molecular docking studies; in the current study, AutoDock version 4.2.6 was utilized [31]. It works on the basis of the Lamarckian genetic algorithm [32]. e target protein DNA gyrase B 3D structure was added with polar hydrogen atoms using AutoDock Tool (ADT). Both ligand and the protein structures were converted into PDBQT format. e grid was fixed with 80 × 80 x 80 size, which covers the active site of the protein. e desolvation map and grid map was created by the AutoGrid program. e algorithm was set with energy evaluation of 2500000, population of 150, 27000 generations, 0.02 mutation rate, and 0.8 crossover rate. Next, the AutoDock run was simulated, and the binding energy of protein-ligand affinity was checked to identify the best binding affinity. AutoDock produces empirical scoring functions; it is the combination of hydrogen bonding, van der Waals, electrostatic, hydrophobicity, entropy, and desolvation energies. e unit of binding energy is kcal/mol. e higher binding affinity between the complexes reflects the increased intermolecular forces. e lowest binding energy stabilizes the complex. e hydrophobic interactions were also explored using LigPlot+ v.2.2 analysis.

Absorption, Distribution, Metabolism, Excretion, and Toxicity Prediction.
e pharmacokinetic properties such as absorption (A), distribution (D), metabolism (M), excretion (E), and toxicity (T) were predicted through the online server ADMET structure activity relationship (admetSAR). e database has structure and text search options, and also, the database collects, curates, and holds ADMET-related properties data from the available literature. e database Evidence-Based Complementary and Alternative Medicine

Results and Discussion
TB, malaria, and acquired immune deficiency syndrome (AIDS) are the top most fatal infectious disease becoming a worldwide public health threat. e development of TB causes increased morbidity and mortality [34]. Different combinations of drugs have been given in treatment. ere are several factors that have been reported for the failure of TB treatment, namely, delay in diagnosis, lack of effective drug administration, decreased accessibility of inexpensive, low toxic and effective drugs, prolong intake of drugs, lack of adherence to drug regimen, and rise of drug-resistant TB strains. Hence, to overcome the resistance mechanism, novel compounds were synthesized against the M. tuberculosis, and molecular modelling studies were carried out to find the potential drug candidate. e DoGSite server predicted the binding site pocket on the DNA gyrase B, which covered the area of 2834.65Å 3 , surface 3418.38Å 2 , and the drug score is 0.81. e drug score falls between 0 and 1, and the higher drug score reveals the high potential of the predicted pocket on the protein (Figure 3).

Evidence-Based Complementary and Alternative Medicine
Similar computational studies have been reported in the screening of antimicrobial compounds. Pertersen et al. [35] have used the virtual screening techniques to identify novel inhibitor for M. tuberculosis 3-dehydroquinate. Similarly, another study has employed pharmacophore-based virtual screening to identify the pyrazolo [1,5-a] pyrimidine derivatives against InhA of M. tuberculosis [36,37].
Here, the first-line antimicrobial drugs rifampicin and isoniazid have lower energies −6.7 kcal/mol and −5.0 kcal/ mol than the benzaldehyde thiosemicarbazone derivatives     e second-line drug levofloxacin is −7.0 kcal/mol.
From the virtual screening and molecular docking studies, it is confirmed that compound 13 could serve as a potential inhibitor for M. tuberculosis. Further preclinical studies have to be conducted to confirm the antimicrobial activity of the compound.

ADMET Prediction.
ADMET prediction of benzaldehyde thiosemicarbazone derivatives using admetSAR server results in intestinal absorption capacity of the compounds 13, 21, and 20. ese compounds also act as a nonsubstrate and noninhibitor for P-glycoprotein. Compound 13 serves as an inhibitor for renal organic cation transporter and other two compounds 20 and 21 as a noninhibitor. In the metabolism, for cytochrome P450 2C9, 2D6, and 3A4, the compounds 13, 21, and 20 possess nonsubstrate property and for cytochrome P450 1A2, 2C9, 2D6, 2C19, and 3A4, it results in noninhibitory activity. Toxicity prediction reveals the non-AMES toxicity and noncarcinogenic activity. Overall, the ADMET prediction further confirms the bioavailability of the compounds (Table 3).

Conclusion
e molecular modelling studies on 30 benzaldehyde thiosemicarbazone derivatives reveals the best binding energy of compound 13, compound 21, and compound 20 against DNA gyrase B of M. tuberculosis. To further confirm the activity of these compounds, preclinical studies have to be conducted.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest.