The drug binding sites and transport mechanism of the RND pumps from Mycobacterium tuberculosis: Insights from molecular dynamics simulations

Inhibition Mechanism of Anti-TB Drug SQ109: Allosteric Inhibition of TMM Translocation of Mycobacterium Tuberculosis MmpL3 Transporter

The mycolic acid transporter MmpL3 is driven by proton motive forces (PMF) and functions via an antiport mechanism. Although the crystal structures of the Mycobacterium smegmatis MmpL3 transporter alone and in complex with a trehalose monomycolate (TMM) substrate and an antituberculosis drug candidate SQ109 under Phase 2b-3 Clinical Trials are available, no water and no conformational change in MmpL3 were observed in these structures to explain SQ109's inhibition mechanism of proton and TMM transportation. In this study, molecular dynamics simulations of both apo form and inhibitor-bound MmpL3 in an explicit membrane were used to decipher the inhibition mechanism of SQ109. In the apo system, the close-open motion of the two TM domains, likely driven by the proton translocation, drives the close-open motion of the two PD domains, presumably allowing for TMM translocation. In contrast, in the holo system, the two PD domains are locked in a closed state, and the two TM domains are locked in an off pathway wider open state due to the binding of the inhibitor. Consistent with the close-open motion of the two PD domains, TMM entry size changes in the apo system, likely loading and moving the TMM, but does not vary much in the holo system and probably impair the movement of the TMM. Furthermore, we observed that water molecules passed through the central channel of the MmpL3 transporter to the cytoplasmic side in the apo system but not in the holo system, with a mean passing time of ∼135 ns. Because water wires play an essential role in transporting protons, our findings shed light on the importance of PMF in driving the close-open motion of the two TM domains. Interestingly, the key channel residues involved in water passage display considerable overlap with conserved residues within the MmpL protein family, supporting their critical function role.

Structure-based in silico approaches for drug discovery against Mycobacterium tuberculosis

Mycobacterium tuberculosis is the causative agent of TB and was estimated to cause 1.4 million death in 2019, alongside 10 million new infections. Drug resistance is a growing issue, with multi-drug resistant infections representing 3.3% of all new infections, hence novel antimycobacterial drugs are urgently required to combat this growing health emergency. Alongside this, increased knowledge of gene essentiality in the pathogenic organism and larger compound databases can aid in the discovery of new drug compounds. The number of protein structures, X-ray based and modelled, is increasing and now accounts for greater than > 80% of all predicted M. tuberculosis proteins; allowing novel targets to be investigated. This review will focus on structure-based in silico approaches for drug discovery, covering a range of complexities and computational demands, with associated antimycobacterial examples. This includes molecular docking, molecular dynamic simulations, ensemble docking and free energy calculations. Applications of machine learning onto each of these approaches will be discussed. The need for experimental validation of computational hits is an essential component, which is unfortunately missing from many current studies. The future outlooks of these approaches will also be discussed.

Auxin transport mechanism of membrane transporter encoded by AEC gene of Bacillus licheniformis isolated from metagenome of Tapta Kund Hotspring of Uttrakhand, India

Members of group Bacillus are most widely occurring microbes in agricultural soil and they affect crop health in various ways. They directly stimulate plant growth either by augmenting nutrients availability, invigorating plants' defence mechanisms; repressing soil-borne phytopathogens or by producing growth-regulating hormones like auxins and cytokinins. It is a well known fact that indole-3- acetic acid (a type of auxin) is a vital biologically active phytohormone excreted by certain Bacillus species, but its molecular mechanism has not yet been described. In this study, the auxin efflux carrier gene is isolated from the metagenome of the Tapta Kund hot spring, Uttrakhand, India. In addition, auxin efflux carrier (AEC) transporter protein of Bacillus licheniformis is modeled and the 318 amino acid residues long protein was found homologous to the apical sodium-dependent bile acid transporter (ASBT) of Yersinia frederiksnii, with 10 transmembrane segments (TM1-10) split into different domains: a panel domain defined by TM1, 2, 6 and 7; and a core domain defined by TM3-5 and 8-10. Finally, the predicted Bacillus licheniformis AEC protein has also been phylogenetically evaluated and its detailed molecular transport mechanism was worked out using molecular dynamics simulation analysis. Conclusively, this study demonstrates the efflux mechanism of the substrate, Indole 3- acetic acid by AEC transporter protein.

Efflux mechanism and pathway of verapamil pumping by human P-glycoprotein

Multidrug resistance (MDR) caused by overexpressed permeability-glycoprotein (P-gp) in cancer cells is the main barrier for the cure of cancers. P-gp can pump many chemotherapeutic drugs, which is a viable target to overcome P-gp-mediated MDR by efficient inhibitors of P-gp. However, limited understanding of the efflux mechanism by human P-gp hinders the development of efficient inhibitors. Herein, the transport of a P-gp inhibitor, verapamil, by human P-gp has been investigated using targeted molecular dynamics simulations and energetics analysis based on our previous research on the transport of a drug (doxorubicin). The energetics analysis identifies that the driving forces for the transport of verapamil are electrostatic repulsions contributed by the positively charged residues in the initial stage and then hydrophobic interactions contributed by the important residues in the later stage. This scenario is generally consistent with that in the transport of doxorubicin. However, the positively charged residues and the important residues for the transport of verapamil are incompletely consistent with the relative residues for the transport of doxorubicin. Moreover, the binding free energy contributions of the positively charged residues for the transport of verapamil are generally higher than them for the transport of doxorubicin, while the important residues constitute significantly different binding free energy compositions in the transports of the two substrates. Consequently, the pathway for the transport of verapamil is identified, which shares only two residues (F336 and M986) with the pathway of doxorubicin. This may imply the weak competitiveness of verapamil with doxorubicin in the substrate efflux. Taken together, this work provided new insights into the efflux mechanisms by human P-gp and would be beneficial in the design of potent P-gp inhibitors.

In Silico Strategies in Tuberculosis Drug Discovery

Tuberculosis (TB) remains a serious threat to global public health, responsible for an estimated 1.5 million mortalities in 2018. While there are available therapeutics for this infection, slow-acting drugs, poor patient compliance, drug toxicity, and drug resistance require the discovery of novel TB drugs. Discovering new and more potent antibiotics that target novel TB protein targets is an attractive strategy towards controlling the global TB epidemic. In silico strategies can be applied at multiple stages of the drug discovery paradigm to expedite the identification of novel anti-TB therapeutics. In this paper, we discuss the current TB treatment, emergence of drug resistance, and the effective application of computational tools to the different stages of TB drug discovery when combined with traditional biochemical methods. We will also highlight the strengths and points of improvement in in silico TB drug discovery research, as well as possible future perspectives in this field.

Structure-Guided Drug Design of 6-Substituted Adenosine Analogues as Potent Inhibitors of Mycobacterium tuberculosis Adenosine Kinase

Mycobacterium tuberculosis adenosine kinase (MtbAdoK) is an essential enzyme of Mtb and forms part of the purine salvage pathway within mycobacteria. Evidence suggests that the purine salvage pathway might play a crucial role in Mtb survival and persistence during its latent phase of infection. In these studies, we adopted a structural approach to the discovery, structure-guided design, and synthesis of a series of adenosine analogues that displayed inhibition constants ranging from 5 to 120 nM against the enzyme. Two of these compounds exhibited low micromolar activity against Mtb with half maximal effective inhibitory concentrations of 1.7 and 4.0 μM. Our selectivity and preliminary pharmacokinetic studies showed that the compounds possess a higher degree of specificity against MtbAdoK when compared with the human counterpart and are well tolerated in rodents, respectively. Finally, crystallographic studies showed the molecular basis of inhibition, potency, and selectivity and revealed the presence of a potentially therapeutically relevant cavity unique to the MtbAdoK homodimer.

Insights on Mycobacterium leprae Efflux Pumps and Their Implications in Drug Resistance and Virulence

Drug resistance in Mycobacterium leprae is assumed to be due to genetic alterations in the drug targets and reduced cell wall permeability. However, as observed in Mycobacterium tuberculosis, drug resistance may also result from the overactivity of efflux systems, which is mostly unexplored. In this perspective, we discuss known efflux pumps involved in M. tuberculosis drug resistance and virulence and investigate similar regions in the genome of M. leprae. In silico analysis reveals that the major M. tuberculosis efflux pumps known to be associated with drug resistance and virulence have been retained during the reductive evolutionary process that M. leprae underwent, e.g., RND superfamily, the ABC transporter BacA, and the MFS P55. However, some are absent (DinF, MATE) while others are derepressed (Mmr, SMR) in M. leprae reflecting the specific environment where M. leprae may live. The occurrence of several multidrug resistance efflux transporters shared between M. leprae and M. tuberculosis reveals potential implications in drug resistance and virulence. The conservation of the described efflux systems in M. leprae upon genome reduction indicates that these systems are potentially required for its intracellular survival and lifestyle. They potentially are involved in M. leprae drug resistance, which could hamper leprosy treatment success. Studying M. leprae efflux pumps as new drug targets is useful for future leprosy therapeutics, enhancing the global efforts to eradicate endemic leprosy, and prevent the emergence of drug resistance in afflicted countries.

Overcoming antibiotic resistance: Is siderophore Trojan horse conjugation an answer to evolving resistance in microbial pathogens?

Comparative study of siderophore biosynthesis pathway in pathogens provides potential targets for antibiotics and host drug delivery as a part of computationally feasible microbial therapy. Iron acquisition using siderophore models is an essential and well established model in all microorganisms and microbial infections a known to cause great havoc to both plant and animal. Rapid development of antibiotic resistance in bacterial as well as fungal pathogens has drawn us at a verge where one has to get rid of the traditional way of obstructing pathogen using single or multiple antibiotic/chemical inhibitors or drugs. 'Trojan horse' strategy is an answer to this imperative call where antibiotic are by far sneaked into the pathogenic cell via the siderophore receptors at cell and outer membrane. This antibiotic once gets inside, generates a 'black hole' scenario within the opportunistic pathogens via iron scarcity. For pathogens whose siderophore are not compatible to smuggle drug due to their complex conformation and stiff valence bonds, there is another approach. By means of the siderophore biosynthesis pathways, potential targets for inhibition of these siderophores in pathogenic bacteria could be achieved and thus control pathogenic virulence. Method to design artificial exogenous siderophores for pathogens that would compete and succeed the battle of intake is also covered with this review. These manipulated siderophore would enter pathogenic cell like any other siderophore but will not disperse iron due to which iron inadequacy and hence pathogens control be accomplished. The aim of this review is to offer strategies to overcome the microbial infections/pathogens using siderophore.

Evolution of structural fitness and multifunctional aspects of mycobacterial RND family transporters

Drug resistance is a major concern due to the evolution and emergence of pathogenic bacterial strains with novel strategies to resist the antibiotics in use. Mycobacterium tuberculosis (Mtb) is one of such pathogens with reported strains, which are not treatable with any of the available anti-TB drugs. This scenario has led to the need to look for some novel drug targets in Mtb, which may be exploited to design effective treatment strategies against the infection. The goal of this review is to discuss one such class of emerging drug targets in Mtb. MmpL (mycobacterial membrane protein large) proteins from Mtb are reported to be involved in multi-substrate transport including drug efflux and considered as one of the contributing factors for the emergence of multidrug-resistant strains. MmpL proteins belong to resistance nodulation division permeases superfamily of membrane transporters, which are viably and pathogenetically important and their inhibition could be lethal for the bacteria.

Excavating the surface-associated and secretory proteome of Mycobacterium leprae for identifying vaccines and diagnostic markers relevant immunodominant epitopes

For centuries, Mycobacterium leprae, etiological agent of leprosy, has been afflicting mankind regardless of extensive use of live-attenuated vaccines and antibiotics. Surface-associated and secretory proteins (SASPs) are attractive targets against bacteria. We have integrated biological knowledge with computational approaches and present a proteome-wide identification of SASPs. We also performed computational assignment of immunodominant epitopes as coordinates of prospective antigenic candidates in most important class of SASPs, the outer membrane proteins (OMPs). Exploiting the known protein sequence and structural characteristics shared by the SASPs from bacteria, 17 lipoproteins, 11 secretory and 19 novel OMPs (including 4 essential proteins) were identified in M. leprae As OMPs represent the most exposed antigens on the cell surface, their immunoinformatics analysis showed that the identified 19 OMPs harbor T-cell MHC class I epitopes and class II epitopes against HLA-DR alleles (54), while 15 OMPs present potential T-cell class II epitopes against HLA-DQ alleles (6) and 7 OMPs possess T-cell class II epitopes against HLA-DP alleles (5) of humans. Additionally, 11 M. leprae OMPs were found to have B-cell epitopes and these may be considered as prime candidates for the development of new immunotherapeutics against M. leprae.

Structural basis of transport function in major facilitator superfamily protein from Trichoderma harzianum

Trichothecenes are the sesquiterpenes secreted by Trichoderma spp. residing in the rhizosphere. These compounds have been reported to act as plant growth promoters and bio-control agents. The structural knowledge for the transporter proteins of their efflux remained limited. In this study, three-dimensional structure of Thmfs1 protein, a trichothecene transporter from Trichoderma harzianum, was homology modelled and further Molecular Dynamics (MD) simulations were used to decipher its mechanism. Fourteen transmembrane helices of Thmfs1 protein are observed contributing to an inward-open conformation. The transport channel and ligand binding sites in Thmfs1 are identified based on heuristic, iterative algorithm and structural alignment with homologous proteins. MD simulations were performed to reveal the differential structural behaviour occurring in the ligand free and ligand bound forms. We found that two discrete trichothecene binding sites are located on either side of the central transport tunnel running from the cytoplasmic side to the extracellular side across the Thmfs1 protein. Detailed analysis of the MD trajectories showed an alternative access mechanism between N and C-terminal domains contributing to its function. These results also demonstrate that the transport of trichodermin occurs via hopping mechanism in which the substrate molecule jumps from one binding site to another lining the transport tunnel.

Structure-Function Profile of MmpL3, the Essential Mycolic Acid Transporter from Mycobacterium tuberculosis

The MmpL family of proteins translocates complex (glyco)lipids and siderophores across the cell envelope of mycobacteria and closely related Corynebacteriaceae and plays important roles in the biogenesis of the outer membrane of these organisms. Despite their significance in the physiology and virulence of Mycobacterium tuberculosis, and from the perspective of developing novel antituberculosis agents, little is known about their structure and mechanism of translocation. In this study, the essential mycobacterial mycolic acid transporter, MmpL3, and its orthologue in Corynebacterium glutamicum, CmpL1, were investigated as prototypical MmpL proteins to gain insight into the transmembrane topology, tertiary and quaternary structures, and functional regions of this transporter family. The combined genetic, biochemical, and biophysical studies indicate that MmpL3 and CmpL1 are structurally similar to Gram-negative resistance-nodulation and division efflux pumps. They harbor 12 transmembrane segments interrupted by two large soluble periplasmic domains and function as homotrimers to export long-chain (C₂₂-C₉₀) mycolic acids, possibly in their acetylated form, esterified to trehalose. The mapping of a number of functional residues within the middle region of the transmembrane domain of MmpL3 shows a striking overlap with mutations associated with resistance to MmpL3 inhibitors. The results suggest that structurally diverse inhibitors of MmpL3 all target the proton translocation path of the transporter and that multiresistance to these inhibitors is enabled by conformational changes in MmpL3.