Alanine mutation of the catalytic sites of Pantothenate Synthetase causes distinct conformational changes in the ATP binding region

Efficient carbon flux allocation towards D-pantothenic acid production via growth-decoupled strategy in Escherichia coli

For industrial strain construction, rational allocation of carbon flux is of paramount importance especially for decoupling cell growth and chemical productions to get maximum titer, rate, yield (TRY), which become Gordian Knot. Here, a temperature-sensitive switch and genetic circuits was used for effectively decoupling cell growth from D-pantothenic acid (DPA) production, along with systematically metabolic engineering including blocking redundant pathways of pyruvate and enhancing DPA driving force. Afterwards, rapid biomass accumulation only happened during growth stage, and subsequent high-efficient DPA production was initiated with reducing fermentation temperature. Finally, 97.20 g/L DPA and 0.64 g/g glucose conversion rate were achieved in 5-liter fed-batch fermentation. These undisputedly represent a milestone for the biosynthesis of DPA. With using strategies for decoupling cell growth from chemical productions, it would serve as "Alexander's sword" to cut Gordian Knot to get industrial chassis cells with excellent TRY for de novo biosynthesis of valuable chemicals.

Exploring optimal drug targets through subtractive proteomics analysis and pangenomic insights for tailored drug design in tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis, ranks among the top causes of global human mortality, as reported by the World Health Organization's 2022 TB report. The prevalence of M. tuberculosis strains that are multiple and extensive-drug resistant represents a significant barrier to TB eradication. Fortunately, having many completely sequenced M. tuberculosis genomes available has made it possible to investigate the species pangenome, conduct a pan-phylogenetic investigation, and find potential new drug targets. The 442 complete genome dataset was used to estimate the pangenome of M. tuberculosis. This study involved phylogenomic classification and in-depth analyses. Sequential filters were applied to the conserved core genome containing 2754 proteins. These filters assessed non-human homology, virulence, essentiality, physiochemical properties, and pathway analysis. Through these intensive filtering approaches, promising broad-spectrum therapeutic targets were identified. These targets were docked with FDA-approved compounds readily available on the ZINC database. Selected highly ranked ligands with inhibitory potential include dihydroergotamine and abiraterone acetate. The effectiveness of the ligands has been supported by molecular dynamics simulation of the ligand-protein complexes, instilling optimism that the identified lead compounds may serve as a robust basis for the development of safe and efficient drugs for TB treatment, subject to further lead optimization and subsequent experimental validation.

Lead identification against Mycobacterium tuberculosis using highly enriched active molecules against pantothenate synthetase

The Pantothenate synthetase (PS) from the Mycobacterium tuberculosis (Mtb) holds a crucial role in the survival and robust proliferation of bacteria through its catalysis of coenzyme A and acyl carrier protein synthesis. The present study undertook the PS drug target in complex with a co-crystallized ligand and subjected it to docking and virtual screening approaches. The experimental design encompassed three discrete datasets: an active dataset featuring 136 compounds, an inactive dataset comprising 56 compounds, and a decoys dataset curated from the zinc library, comprising an extensive compilation of approximately 53,000 compounds. The compounds' binding energies were observed to be in the range of -5 to ∼-14 kcal/mol. Additionally, binding energy results were further refined through Enrichment Factor analysis (EF). EF is a new statistical approach which uses the scores obtained from docking-based virtual screening and predicts the precision of the scoring function. Remarkably, the Enrichment Factor (EF) analysis produced exceptionally favorable outcomes, attaining an EF of approximately 49% within the uppermost 1% fraction of the compound distribution. Finally, a total of eight compounds, evenly distributed between the active dataset and the decoys dataset, emerged as potent inhibitors of the Pantothenate synthetase (PS) enzyme. The analysis of inhibition constants and binding energy revealed a notable correlation, with an r-squared value (r2) of 0.912 between the two parameters. Furthermore, the shortlisted compounds were subjected to 100 ns MD simulation to determine their stability and dynamics behavior. The decoy compounds that have been identified, exhibiting properties comparable to the active compounds, are postulated as potential candidates for targeting the Pantothenate synthetase (PS) enzyme to treat Mtb infection. Nevertheless, in the pursuit of a comprehensive investigation, it is advisable to undertake additional experimental validation as a component of the subsequent study.Communicated by Ramaswamy H. Sarma.

Dynamozones are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease

Proteins are not static but are flexible molecules that can adopt many different conformations. The HIV-1 protease is an important target for the development of therapies to treat AIDS, due to its critical role in the viral life cycle. We investigated several dynamics studies on the HIV-1 protease families to illustrate the significance of examining the dynamic behaviors and molecular motions for an entire understanding of their dynamics-structure-function relationships. Using computer simulations and principal component analysis approaches, the dynamics data obtained revealed that: (i) The flap regions are the most obvious sign of the evolution of conformational dynamics in HIV-1 protease; (ii) There are dynamic structural regions in some proteins that contribute to the biological function and allostery of proteins via appropriate flexibility. These regions are a clear sign of the evolution of conformational dynamics of proteins, which we call dynamozones. The flap regions are one of the most important dynamozones members that are critical for HIV-1 protease function. Due to the existence of other members of dynamozones in different proteins, we propose to consider dynamozones as a footprint of the evolution of the conformational dynamics of proteins.

Calendulaglycoside A showing potential activity against SARS-CoV-2 main protease: Molecular docking, molecular dynamics, and SAR studies

BACKGROUND AND AIM

The discovery of drugs capable of inhibiting SARS-CoV-2 is a priority for human beings due to the severity of the global health pandemic caused by COVID-19. To this end, natural products can provide therapeutic alternatives that could be employed as an effective safe treatment for COVID-19.

EXPERIMENTAL PROCEDURE

Twelve compounds were isolated from the aerial parts of C. officinalis L. and investigated for their inhibitory activities against SARS-CoV-2 Mpro compared to its co-crystallized N3 inhibitor using molecular docking studies. Furthermore, a 100 ns MD simulation was performed for the most active two promising compounds, Calendulaglycoside A (SAP5) and Osteosaponin-I (SAP8).

RESULTS AND CONCLUSION

At first, molecular docking studies showed interesting binding scores as compared to the N3 inhibitor. Calendulaglycoside A (SAP5) achieved a superior binding than the co-crystallized inhibitor indicating promising affinity and intrinsic activity towards the Mpro of SARS-CoV-2 as well. Moreover, findings illustrated preferential stability for SAP5 within the Mpro pocket over that of N3 beyond the 40 ns MD simulation course. Structural preferentiality for triterpene-Mpro binding highlights the significant role of 17β-glucosyl and carboxylic 3α-galactosyl I moieties through high electrostatic interactions across the MD simulation trajectories. Furthermore, this study clarified a promising SAR responsible for the antiviral activity against the SARS-CoV-2 Mpro and the design of new drug candidates targeting it as well. The above findings could be promising for fast examining the previously isolated triterpenes both pre-clinically and clinically for the treatment of COVID-19.

Crocin Inhibits the Fibrillation of Human α-synuclein and Disassembles Mature Fibrils: Experimental Findings and Mechanistic Insights from Molecular Dynamics Simulation

The aggregation of human alpha-synuclein (hαS) is pivotally implicated in the development of most types of synucleinopathies. Molecules that can inhibit or reverse the aggregation process of amyloidogenic proteins have potential therapeutic value. The anti-aggregating activity of multiple carotenoid compounds has been reported over the past decades against a growing list of amyloidogenic polypeptides. Here, we aimed to determine whether crocin, the main carotenoid glycoside component of saffron, would inhibit hαS aggregation or could disassemble its preformed fibrils. By employing a series of biochemical and biophysical techniques, crocin was exhibited to inhibit hαS fibrillation in a dose-dependent fashion by stabilizing very early aggregation intermediates in off-pathway non-toxic conformations with little β-sheet content. We also observed that crocin at high concentrations could efficiently destabilize mature fibrils and disassemble them into seeding-incompetent intermediates by altering their β-sheet conformation and reshaping their structure. Our atomistic molecular dynamics (MD) simulations demonstrated that crocin molecules bind to both the non amyloid-β component (NAC) region and C-terminal domain of hαS. These interactions could thereby stabilize the autoinhibitory conformation of the protein and prevent it from adopting aggregation-prone structures. MD simulations further suggested that ligand molecules prefer to reside longitudinally along the fibril axis onto the edges of the inter-protofilament interface where they establish hydrogen and hydrophobic bonds with steric zipper stabilizing residues. These interactions turned out to destabilize hαS fibrils by altering the interstrand twist angles, increasing the rigidity of the fibril core, and elevating its radius of gyration. Our findings suggest the potential pharmaceutical implication of crocin in synucleinopathies.

Investigating the structure-activity relationship of marine natural polyketides as promising SARS-CoV-2 main protease inhibitors

Since its first report in December 2019, the novel coronavirus virus, SARS-CoV-2, has caused an unprecedented global health crisis and economic loss imposing a tremendous burden on the worldwide finance, healthcare system, and even daily life. Even with the introduction of different preventive vaccines, there is still a dire need for effective antiviral therapeutics. Nature has been considered as the historical trove of drug discovery and development, particularly in cases of worldwide crises. Herein, a comprehensive in silico investigation of a highly focused chemical library of 34 pederin-structurally related marine compounds, belonging to four polyketides families, was initiated against the SARS-CoV-2 main protease, Mpro, being the key replicating element of the virus and main target in many drugs development programs. Two of the most potent SARS-CoV-2 Mpro co-crystallized inhibitors, O6K and N3, were added to the tested database as reference standards. Through molecular docking simulation, promising compounds including Pederin (1), Dihydro-onnamide A (11), Onnamide C (14), Pseudo-onnamide A (17), and Theopederin G (29) have been identified from different families based on their superior ligand–protein energies and relevant binding profiles with the key Mpro pocket residues. Thermodynamic behaviors of the identified compounds were investigated through 200 ns all-atom molecular dynamics simulation illustrating their significant stability and pocket accommodation. Furthermore, structural activity preferentiality was identified for the pederin-based marine compounds highlighting the importance of the terminal guanidine and cyclic hemiacetal linker, and the length of the sidechain. Our findings highlight the challenges of targeting SARS-CoV-2 Mpro as well as recommending further in vitro and in vivo studies regarding the examined marine products either alone or in combination paving the way for promising lead molecules.

Marine natural polyketides showed promising SARS-CoV-2 main protease inhibitory activities.

Molecular Docking and Dynamics Simulation Study of Hyrtios erectus Isolated Scalarane Sesterterpenes as Potential SARS-CoV-2 Dual Target Inhibitors

Presently, the world is under the toll of pandemic coronavirus disease-2019 (COVID-19) outbreak caused by SARS-CoV-2. Lack of effective and safe therapeutics has stressed the scientific community for developing novel therapeutics capable of alleviating and stopping this pandemic. Within the presented study, molecular docking, ADME properties and all-atom molecular dynamic (MD) simulation, along with two standard antiviral agents (lopinavir and benzopurpurin-4B), were applied to investigate 15 scalaranes sesterterpenes natural compounds, purified from the Red Sea marine sponge Hyrtios erectus, as potential COVID-19 dual-target inhibitors. Following multi-step docking within COVID-19 main protease and Nsp15 endoribonuclease cavities, nine promising drug-like compounds exhibited higher docking scores as well as better interactions with the target's crucial residues than those of reference ligands. Compounds 2, 6, 11, and 15, were predicted to simultaneously subdue the activity of the two COVID-19 targets. Dynamics behavior of the best-docked molecules, compounds 15 and 6, within COVID-19 target pockets showed substantial stability of ligand-protein complexes as presented via several MD simulation parameters. Furthermore, calculated free-binding energies from MD simulation illustrated significant ligand's binding affinity towards respective target pockets. All provided findings supported the utility of scalarane-based sesterterpenes, particularly compounds 15 and 6, as promising lead candidates guiding the development of effective therapeutics against SARS-CoV-2.

In a search for potential drug candidates for combating COVID-19: computational study revealed salvianolic acid B as a potential therapeutic targeting 3CLpro and spike proteins

The global prevalence of COVID-19 disease and the overwhelming increase in death toll urge scientists to discover new effective drugs. Although the drug discovery process is a challenging and time-consuming, fortunately, the plant kingdom was found to have many active therapeutics possessing broad-spectrum antiviral activity including those candidates active against severe acute respiratory syndrome coronaviruses (SARS-CoV). Herein, nine traditional Chinese medicinal plant constituents from different origins (Glycyrrhizin 1, Lycorine 2, Puerarin 3, Daidzein 4, Daidzin 5, Salvianolic acid B 6, Dihydrotanshinone I 7, Tanshinone I 8, Tanshinone IIa 9) previously reported to exhibit antiviral activity against SARS-CoV were virtually screened in silico (molecular docking) as potential inhibitors of SARS-CoV-2 target proteins. The tested medicinal plant compounds were in silico screened for their activity against two key SARS-CoV-2 target proteins; 3CLpro, and Spike binding-domain proteins. Among the tested medicinal plant compounds, Salvianolic acid B 6 (Sal-B) showed promising binding affinities against the two specified SARS-CoV-2 target proteins compared to the reference standards used. Hence molecular dynamics simulations followed by calculating the free-binding energy were carried out for Sal-B providing information on its affinity, stability, and thermodynamic behavior within the two SARS-CoV-2 target proteins as well as key ligand-protein binding aspects. Besides, the quantum mechanical calculations showed that Sal-B can adopt different conformations due to the existence of various rotatable bonds. Therefore, the enhanced antiviral activity of Sal-B among other studied compounds can be also attributed to the structural flexibility of Sal-B. Our study gives an explanation of the structure activity relationship required for targeting SARS-CoV-2 3CLpro and Spike proteins and also facilitates the future design and synthesis of new potential drugs exhibiting better affinity and specificity. Besides, an ADME study was carried out on screened compounds and reference controls revealing their pharmacokinetics properties.Communicated by Ramaswamy H. Sarma.

Vitamin in the Crosshairs: Targeting Pantothenate and Coenzyme A Biosynthesis for New Antituberculosis Agents

Despite decades of dedicated research, there remains a dire need for new drugs against tuberculosis (TB). Current therapies are generations old and problematic. Resistance to these existing therapies results in an ever-increasing burden of patients with disease that is difficult or impossible to treat. Novel chemical entities with new mechanisms of action are therefore earnestly required. The biosynthesis of coenzyme A (CoA) has long been known to be essential in Mycobacterium tuberculosis (Mtb), the causative agent of TB. The pathway has been genetically validated by seminal studies in vitro and in vivo. In Mtb, the CoA biosynthetic pathway is comprised of nine enzymes: four to synthesize pantothenate (Pan) from l-aspartate and α-ketoisovalerate; five to synthesize CoA from Pan and pantetheine (PantSH). This review gathers literature reports on the structure/mechanism, inhibitors, and vulnerability of each enzyme in the CoA pathway. In addition to traditional inhibition of a single enzyme, the CoA pathway offers an antimetabolite strategy as a promising alternative. In this review, we provide our assessment of what appear to be the best targets, and, thus, which CoA pathway enzymes present the best opportunities for antitubercular drug discovery moving forward.

Computational SNP Analysis and Molecular Simulation Revealed the Most Deleterious Missense Variants in the NBD1 Domain of Human ABCA1 Transporter

The ATP-binding cassette transporter A1 (ABCA1) is a membrane-bound exporter protein involved in regulating serum HDL level by exporting cholesterol and phospholipids to load up in lipid-poor ApoA-I and ApoE, which allows the formation of nascent HDL. Mutations in the ABCA1 gene, when presents in both alleles, disrupt the canonical function of ABCA1, which associates with many disorders related to lipid transport. Although many studies have reported the phenotypic effects of a large number of ABCA1 variants, the pathological effect of non-synonymous polymorphisms (nsSNPs) in ABCA1 remains elusive. Therefore, aiming at exploring the structural and functional consequences of nsSNPs in ABCA1, in this study, we employed an integrated computational approach consisting of nine well-known in silico tools to identify damaging SNPs and molecular dynamics (MD) simulation to get insights into the magnitudes of the damaging effects. In silico tools revealed four nsSNPs as being most deleterious, where the two SNPs (G1050V and S1067C) are identified as the highly conserved and functional disrupting mutations located in the NBD1 domain. MD simulation suggested that both SNPs, G1050V and S1067C, changed the overall structural flexibility and dynamics of NBD1, and induced substantial alteration in the structural organization of ATP binding site. Taken together, these findings direct future studies to get more insights into the role of these variants in the loss of the ABCA1 function.

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.

Pyrazinamide drug resistance in RpsA mutant (∆438A) of Mycobacterium tuberculosis: Dynamics of essential motions and free-energy landscape analysis

Pyrazinamide is an essential first-line antitubercular drug which plays pivotal role in tuberculosis treatment. It is a prodrug that requires amide hydrolysis by mycobacterial pyrazinamidase enzyme for conversion into pyrazinoic acid (POA). POA is known to target ribosomal protein S1 (RpsA), aspartate decarboxylase (PanD), and some other mycobacterial proteins. Spontaneous chromosomal mutations in RpsA have been reported for phenotypic resistance against pyrazinamide. We have constructed and validated 3D models of the native and Δ438A mutant form of RpsA protein. RpsA protein variants were then docked to POA and long range molecular dynamics simulations were carried out. Per residue binding free-energy calculations, free-energy landscape analysis, and essential dynamics analysis were performed to outline the mechanism underlying the high-level PZA resistance conferred by the most frequently occurring deletion mutant of RpsA. Our study revealed the conformational modulation of POA binding site due to the disruptive collective modes of motions and increased conformational flexibility in the mutant than the native form. Residue wise MMPBSA decomposition and protein-drug interaction pattern revealed the difference of energetically favorable binding site in the wild-type (WT) protein in comparison with the mutant. Analysis of size and shape of minimal energy landscape area delineated higher stability of the WT complex than the mutant form. Our study provides mechanistic insights into pyrazinamide resistance in Δ438A RpsA mutant, and the results arising out of this study will pave way for design of novel and effective inhibitors targeting the resistant strains of Mycobacterium tuberculosis.

Molecular mechanism of acetoacetyl-CoA enhanced kinetics for increased bioplastic production from Cupriavidus necator 428

Polyhydroxyalkanoates are gaining importance due to their biodegradable nature and close analogy to plastics. Polyhydroxybutyrate (PHB) is the most widely used bioplastic from polyalkanoate family, which is produced by a legion of bacterial species via phbCAB operon encoding β-ketothiolase (PhaA), NADPH-dependent acetoacetyl-coenzyme A (acetoacetyl-CoA) reductase (PhaB) and polyhydroxyalkanoate synthase (PhaC). Augmentation in the activity of these enzymes is promising for increased PHB production which is achieved by enzyme engineering strategies including non-structural and structural approaches. Our study is deployed on directed evolution-based experimentally reported mutants of PhaB enzyme with increased efficiency due to impact on critical structural factors. We have analyzed and compared the native PhaB with two of its variants Q47L and T173S in complex with their cofactor i.e. NADPH as well as the substrate i.e. acetoacetyl-CoA, via long range molecular dynamics simulations. Interaction profile, MMPBSA, essential dynamics, and free energy landscape analysis revealed that the enzyme efficiency is critically affected by cofactor interactions. It was also observed that mutants have higher equilibrium constant with lesser but optimal affinity for substrate and cofactor than the wild type, which might be the reason for increased efficiency of the mutants via enhanced substrate and cofactor exchange rate. Our study provides insights into the cofactor and substrate binding affinities to PhaB enzyme at atomistic level, which will facilitate designing of highly efficient PhaB enzymes for increased PHB production. Communicated by Ramaswamy H. Sarma.

Drug repurposing against arabinosyl transferase (EmbC) of Mycobacterium tuberculosis: Essential dynamics and free energy minima based binding mechanics analysis

Arabinosyl tranferases (embA, embB, embC) are the key enzymes responsible for biogenesis of arabinan domain of arabinogalactan (AG) and lipoarabinomannan (LAM), two major heteropolysaccharide constituents of the peculiar mycobacterial cell envelope. EmbC is predominantly responsible for LAM synthesis and has been commonly associated with Ethambutol resistance. We have screened the FDA library against EmbC to reposition a drug better than Ethambutol with higher binding affinity to Embc. High throughput virtual screening followed by extra precision docking using Glide gave two best leads i.e. Terlipressin and Amikacin with docking score of -11.39 kcal/mol and -10.71 kcal/mol, respectively. Binding mechanics of the selected drugs was elucidated through long range molecular dynamics simulations (100 ns) using binding free energy rescoring, essential dynamics and free energy minima based approaches, thus revealing the most stable binding modes of Terlipressin with EmbC. Our study establishes the EmbC binding potential of the repurposed drugs Terlipressin and Amikacin.

The PGRS Domain of Mycobacterium tuberculosis PE_PGRS Protein Rv0297 Is Involved in Endoplasmic Reticulum Stress-Mediated Apoptosis through Toll-Like Receptor 4

The genome of Mycobacterium tuberculosis, the causal organism of tuberculosis (TB), encodes a unique protein family known as the PE/PPE/PGRS family, present exclusively in the genus Mycobacterium and nowhere else in the living kingdom, with largely unexplored functions. We describe the functional significance of the PGRS domain of Rv0297, a member of this family. In silico analyses revealed the presence of intrinsically disordered stretches and putative endoplasmic reticulum (ER) localization signals in the PGRS domain of Rv0297 (Rv0297PGRS). The PGRS domain aids in ER localization, which was shown by infecting macrophage cells with M. tuberculosis and by overexpressing the protein by transfection in macrophage cells followed by activation of the unfolded protein response, as evident from increased expression of GRP78/GRP94 and CHOP/ATF4, leading to disruption of intracellular Ca²⁺ homeostasis and increased nitric oxide (NO) and reactive oxygen species (ROS) production. The consequent activation of the effector caspase-8 resulted in apoptosis of macrophages, which was Toll-like receptor 4 (TLR4) dependent. Administration of recombinant Rv0297PGRS (rRv0297PGRS) also exhibited similar effects. These results implicate a hitherto-unknown role of the PGRS domain of the PE_PGRS protein family in ER stress-mediated cell death through TLR4. Since this protein is already known to be present at later stages of infection in human granulomas it points to the possibility of it being employed by M. tuberculosis for its dissemination via an apoptotic mechanism.

Apoptosis is generally thought to be a defense mechanism in protecting the host against Mycobacterium tuberculosis in early stages of infection. However, apoptosis during later stages in lung granulomas may favor the bacterium in disseminating the disease. ER stress has been found to induce apoptosis in TB granulomas, in zones where apoptotic macrophages accumulate in mice and humans. In this study, we report ER stress-mediated apoptosis of host cells by the Rv0297-encoded PE_PGRS5 protein of M. tuberculosis exceptionally present in the pathogenic Mycobacterium genus. The PGRS domain of Rv0297 aids the protein in localizing to the ER and induces the unfolded protein response followed by apoptosis of macrophages. The effect of the Rv0297PGRS domain was found to be TLR4 dependent. This study presents novel insights on the strategies employed by M. tuberculosis to disseminate the disease.

Drug targeted virtual screening and molecular dynamics of LipU protein of Mycobacterium tuberculosis and Mycobacterium leprae

The lipolytic protein LipU was conserved in mycobacterium sp. including M. tuberculosis (MTB LipU) and M. leprae (MLP LipU). The MTB LipU was identified in extracellular fraction and was reported to be essential for the survival of mycobacterium. Therefore to address the problem of drug resistance in pathogen, LipU was selected as a drug target and the viability of finding out some FDA approved drugs as LipU inhibitors in both the cases was explored. Three-dimensional (3D) model structures of MTB LipU and MLP LipU were generated and stabilized through molecular dynamics (MD). FDA approved drugs were screened against these proteins. The result showed that the top-scoring compounds for MTB LipU were Diosmin, Acarbose and Ouabain with the Glide XP score of -12.8, -11.9 and -11.7 kcal/mol, respectively, whereas for MLP LipU protein, Digoxin (-9.2 kcal/mol), Indinavir (-8.2 kcal/mol) and Travoprost (-8.2 kcal/mol) showed highest affinity. These drugs remained bound in the active site pocket of MTB LipU and MLP LipU structure and interaction grew stronger after dynamics. RMSD, RMSF and Rg were found to be persistent throughout the simulation period. Hydrogen bonds along with large number of hydrophobic interactions stabilized the complex structures. Binding free energies obtained through Prime/MM-GBSA were found in the significant range from -63.85 kcal/mol to -34.57 kcal/mol for MTB LipU and -71.33 kcal/mol to -23.91 kcal/mol for MLP LipU. The report suggested high probability of these drugs to demolish the LipU activity and could be probable drug candidates to combat TB and leprosy disease.