Ligand Transport Research Articles

Large-scale annotation of biochemically relevant pockets and tunnels in cognate enzyme–ligand complexes

Tunnels in enzymes with buried active sites are key structural features allowing the entry of substrates and the release of products, thus contributing to the catalytic efficiency. Targeting the bottlenecks of protein tunnels is also a powerful protein engineering strategy. However, the identification of functional tunnels in multiple protein structures is a non-trivial task that can only be addressed computationally. We present a pipeline integrating automated structural analysis with an in-house machine-learning predictor for the annotation of protein pockets, followed by the calculation of the energetics of ligand transport via biochemically relevant tunnels. A thorough validation using eight distinct molecular systems revealed that CaverDock analysis of ligand un/binding is on par with time-consuming molecular dynamics simulations, but much faster. The optimized and validated pipeline was applied to annotate more than 17,000 cognate enzyme–ligand complexes. Analysis of ligand un/binding energetics indicates that the top priority tunnel has the most favourable energies in 75% of cases. Moreover, energy profiles of cognate ligands revealed that a simple geometry analysis can correctly identify tunnel bottlenecks only in 50% of cases. Our study provides essential information for the interpretation of results from tunnel calculation and energy profiling in mechanistic enzymology and protein engineering. We formulated several simple rules allowing identification of biochemically relevant tunnels based on the binding pockets, tunnel geometry, and ligand transport energy profiles.Scientific contributionsThe pipeline introduced in this work allows for the detailed analysis of a large set of protein–ligand complexes, focusing on transport pathways. We are introducing a novel predictor for determining the relevance of binding pockets for tunnel calculation. For the first time in the field, we present a high-throughput energetic analysis of ligand binding and unbinding, showing that approximate methods for these simulations can identify additional mutagenesis hotspots in enzymes compared to purely geometrical methods. The predictor is included in the supplementary material and can also be accessed at https://github.com/Faranehhad/Large-Scale-Pocket-Tunnel-Annotation.git. The tunnel data calculated in this study has been made publicly available as part of the ChannelsDB 2.0 database, accessible at https://channelsdb2.biodata.ceitec.cz/.

Redox Signaling in Endosomes Using the Example of EGF Receptors: A Graphical Review.

Early endosomes represent first-line sorting compartments or even organelles for internalized molecules. They enable the transport of molecules or ligands to other compartments of the cell, such as lysosomes, for degradation or recycle them back to the membrane by various mechanisms. Moreover, early endosomes function as signaling and scaffolding platforms to initiate or prolong distinct signaling pathways. Accordingly, early endosomes have to be recognized as either part of a degradation or recycling pathway. The physical proximity of many ligand-binding receptors with other membrane-bound proteins or complexes such as NADPH oxidases may result in an interaction of second messengers, like reactive oxygen species (ROS) and early endosomes, that promote the correct recognition of individual early endosomes. In fact, redoxosomes comprise an endosomal subsection of signaling endosomes. One example of such potential interaction is epidermal growth factor receptor (EGFR) signaling. Here we summarize recent findings on EGFR signaling as a well-studied example for receptor trafficking and trans-activation and illustrate the interplay between cellular and endosomal ROS.

Structural basis of thiamine transport and drug recognition by SLC19A3

Thiamine (vitamin B1) functions as an essential coenzyme in cells. Humans and other mammals cannot synthesise this vitamin de novo and thus have to take it up from their diet. Eventually, every cell needs to import thiamine across its plasma membrane, which is mainly mediated by the two specific thiamine transporters SLC19A2 and SLC19A3. Loss of function mutations in either of these transporters lead to detrimental, life-threatening metabolic disorders. SLC19A3 is furthermore a major site of drug interactions. Many medications, including antidepressants, antibiotics and chemotherapeutics are known to inhibit this transporter, with potentially fatal consequences for patients. Despite a thorough functional characterisation over the past two decades, the structural basis of its transport mechanism and drug interactions has remained elusive. Here, we report seven cryo-electron microscopy (cryo-EM) structures of the human thiamine transporter SLC19A3 in complex with various ligands. Conformation-specific nanobodies enable us to capture different states of SLC19A3’s transport cycle, revealing the molecular details of thiamine recognition and transport. We identify seven previously unknown drug interactions of SLC19A3 and present structures of the transporter in complex with the inhibitors fedratinib, amprolium and hydroxychloroquine. These data allow us to develop an understanding of the transport mechanism and ligand recognition of SLC19A3.

Structural visualization of inhibitor binding in prolyl oligopeptidase.

The association and dissociation of proteins and ligands are crucial in biophysics for potential drug development [Baron and McCammon, Annu. Rev. Phys. Chem. 64, 151-175 (2013)]. However, identifying and characterizing the reaction pathways for these rare events has been a long-standing challenge. Molecular dynamics (MD) simulations are limited in exploring biophysical processes on experimental timescales, so ligand transport processes through complex transient tunnels formed by proteins during dynamics are often simulated using enhanced sampling MD [Rydzewski and Nowak, Phys. Life Rev. 22-23, 58-74 (2017)]. Erroneously identified ligand binding pathways can affect thermodynamic and kinetic characteristics calculated from MD trajectories. A system that has the potential to be a therapeutic target for neurodegenerative diseases is prolyl oligopeptidase (PREP). This is due to its involvement in promoting protein aggregation and disrupting cellular function through affecting protein-protein interactions (PPI). The recent discovery of a new type of PREP inhibitor that targets PPI raises important questions about the diversity of ligand binding pathways in PREP and their impact on protein dynamics [Pätsi et al., J. Med. Chem. 67, 5421-5436 (2024); Kilpeläinen et al., J. Med. Chem. 66, 7475-7496 (2023); and Walczewska-Szewc et al., Phys. Chem. Chem. Phys. 24, 4366-4373 (2022)]. In this article, using results from enhanced sampling MD, we visually present how the binding process in PREP depends on subtle changes in inhibitors, which could be crucial in treating neurodegenerative disorders.

Recent Updates on Interaction Studies and Drug Delivery of Antimalarials with Serum Albumin Proteins.

This review focuses on recent trends in the binding study of various antimalarial agents with serum albumins in detail. Serum albumin has a significant role in the transport of drugs and endogenous ligands. The nature and magnitude of serum albumin and drug interactions have a tremendous impact on the pharmacological behavior and toxicity of that drug. Binding of drug to serum albumin not only controls its free and active concentration, but also provides a reservoir for a long duration of action. This ultimately affects drug absorption, distribution, metabolism, and excretion. Such interaction determines the actual drug efficacy as the drug action can be correlated with the amount of unbound drug. With the advancement in spectroscopic techniques and simulation studies, binding studies play an increasingly important role in biophysical and biomedical science, especially in the field of drug delivery and development. This review assesses the insight we have gained so far to improve drug delivery and discovery of antimalarials on the basis of a plethora of drug-serum protein interaction studies done so far.

Characterization of Intrinsically Disordered Proteins in Healthy and Diseased States by Nuclear Magnetic Resonance.

Intrinsically Disordered Proteins (IDPs) are active in different cellular procedures like ordered assembly of chromatin and ribosomes, interaction with membrane, protein, and ligand binding, molecular recognition, binding, and transportation via nuclear pores, microfilaments and microtubules process and disassembly, protein functions, RNA chaperone, and nucleic acid binding, modulation of the central dogma, cell cycle, and other cellular activities, post-translational qualification and substitute splicing, and flexible entropic linker and management of signaling pathways. The intrinsic disorder is a precise structural characteristic that permits IDPs/IDPRs to be involved in both one-to-many and many-to-one signaling. IDPs/IDPRs also exert some dynamical and structural ordering, being much less constrained in their activities than folded proteins. Nuclear magnetic resonance (NMR) spectroscopy is a major technique for the characterization of IDPs, and it can be used for dynamic and structural studies of IDPs. This review was carried out to discuss intrinsically disordered proteins and their different goals, as well as the importance and effectiveness of NMR in characterizing intrinsically disordered proteins in healthy and diseased states.

Incorporating Prior Knowledge in the Seeds of Adaptive Sampling Molecular Dynamics Simulations of Ligand Transport in Enzymes with Buried Active Sites.

Because most proteins have buried active sites, protein tunnels or channels play a crucial role in the transport of small molecules into buried cavities for enzymatic catalysis. Tunnels can critically modulate the biological process of protein-ligand recognition. Various molecular dynamics methods have been developed for exploring and exploiting the protein-ligand conformational space to extract high-resolution details of the binding processes, a recent example being energetically unbiased high-throughput adaptive sampling simulations. The current study systematically contrasted the role of integrating prior knowledge while generating useful initial protein-ligand configurations, called seeds, for these simulations. Using a nontrivial system of a haloalkane dehalogenase mutant with multiple transport tunnels leading to a deeply buried active site, simulations were employed to derive kinetic models describing the process of association and dissociation of the substrate molecule. The most knowledge-based seed generation enabled high-throughput simulations that could more consistently capture the entire transport process, explore the complex network of transport tunnels, and predict equilibrium dissociation constants, koff/kon, on the same order of magnitude as experimental measurements. Overall, the infusion of more knowledge into the initial seeds of adaptive sampling simulations could render analyses of transport mechanisms in enzymes more consistent even for very complex biomolecular systems, thereby promoting drug development efforts and the rational design of enzymes with buried active sites.

Mechanistic insights into P-glycoprotein ligand transport and inhibition revealed by enhanced molecular dynamics simulations

P-glycoprotein (P-gp) plays a crucial role in cellular detoxification and drug efflux processes, transitioning between inward-facing (IF) open, occluded, and outward-facing (OF) states to facilitate substrate transport. Its role is critical in cancer therapy, where P-gp contributes to the multidrug resistance phenotype. In our study, classical and enhanced molecular dynamics (MD) simulations were conducted to dissect the structural and functional features of the P-gp conformational states. Our advanced MD simulations, including kinetically excited targeted MD (ketMD) and adiabatic biasing MD (ABMD), provided deeper insights into state transition and translocation mechanisms. Our findings suggest that the unkinking of TM4 and TM10 helices is a prerequisite for correctly achieving the outward conformation. Simulations of the IF-occluded conformations, characterized by kinked TM4 and TM10 helices, consistently demonstrated altered communication between the transmembrane domains (TMDs) and nucleotide binding domain 2 (NBD2), suggesting the implication of this interface in inhibiting P-gp's efflux function. A particular emphasis was placed on the unstructured linker segment connecting the NBD1 to TMD2 and its role in the transporter's dynamics. With the linker present, we specifically noticed a potential entrance of cholesterol (CHOL) through the TM4-TM6 portal, shedding light on crucial residues involved in accommodating CHOL. We therefore suggest that this entry mechanism could be employed for some P-gp substrates or inhibitors. Our results provide critical data for understanding P-gp functioning and developing new P-gp inhibitors for establishing more effective strategies against multidrug resistance.

Revealing the interaction mechanism between bovine serum albumin (BSA) and a fluorescent coumarin derivative: A multispectroscopic and in silico approach

The interaction of biomacromolecules with each other or with the ligands is essential for biological activity. In this context, the molecular recognition of bovine serum albumin (BSA) with 4-(Benzo[1,3]dioxol-5-yloxymethyl)-7-hydroxy-chromen-2-one (4BHC) is explored using multispectroscopic and computational techniques. UV–Vis spectroscopy helped in predicting the conformational variations in BSA. Using fluorescence spectroscopy, the quenching behaviour of the fluorophore upon interaction with the ligand is examined, which is found to be a static type of quenching; fluorescence lifetime studies further verify this. The binding constant is discovered to be in the range of 104 M−1, which indicates the moderate type of association that results in reversible binding, where the transport and release of ligands in the target tissue takes place. Fourier-transform infrared spectroscopy (FT-IR) measurements validate the secondary structure conformational changes of BSA after complexing with 4BHC. The thermodynamic factors obtained through temperature-dependent fluorescence studies suggest that the dominant kind of interaction force is hydrophobic in nature, and the interaction process is spontaneous. The alterations in the surrounding microenvironment of the binding site and conformational shifts in the structure of the protein are studied through 3D fluorescence and synchronous fluorescence studies. Molecular docking and molecular dynamics (MD) simulations agree with experimental results and explain the structural stability throughout the discussion. The outcomes might have possible applications in the field of pharmacodynamics and pharmacokinetics.

Garnet geochronology, major and trace element geochemistry of the Huanggangliang Fe-Sn polymetallic deposit, NE China

Cu-Fe (Au), W-Sn (Mo), Pb-Zn (Ag) and other element assemblages are common in many skarn deposits, but examples of Fe and Sn coexistence in nature are relatively rare. The Huanggangliang deposit in the southern Great Xing'an Range is the largest Fe-Sn polymetallic deposit in northern China, and there are five main mining areas (SK Ⅰ–SK Ⅴ) distributed along a SW–NE trend. SK Ⅰ and SK Ⅲ are the two larger mining areas in the Huanggangliang deposit, and their mineral assemblages and wall-rocks are significantly different: cassiterite and metal sulfides are duncommon in SK Ⅰ (andesite), while they are widely distributed in SK Ⅲ (marble). Garnet is widely distributed in SK Ⅰ and SK Ⅲ mining areas, and provides a good opportunity to explore the mechanism of Fe-Sn assemblage differentiation. Three types of garnet are distinguishable (Grt 1, Grt 2 and Grt 3) based on their mineral associations, textures and locations, and they mainly belong to the andradite-grossular solid solution series. The LA–ICP–MS U-Pb ages of the garnets from SK Ⅰ and SK Ⅲ are 136.0 ± 1.2 Ma (n = 25, MSWD = 1.9) and 136.2 ± 2.6 Ma (n = 37, MSWD = 1.8), respectively, suggesting that magmatic-hydrothermal activity in these two mining areas began at the same time. Grt 1 (located deep in SK Ⅰ and SK Ⅲ) and Grt 2 (located shallowly in SK Ⅰ) are mainly andradite, while Grt 3 (located shallowly in SK Ⅲ) is mainly grossular. They are all enriched in high field strength elements (HFSEs) and depleted in large ion lithophile elements (LILEs), and their REE patterns are enriched in light rare earth elements (LREEs) and depleted in heavy rare earth elements (HREEs). However, Grt Ⅰ-1a and 1b (deep garnets in SK Ⅰ) exhibit obvious positive Eu anomalies, while the other garnet types exhibit negative or no obvious Eu anomalies. Additionally, Grt 3 has high contents of V, Sn, U, Th, Zn and Y, and low contents of W and Mo. This pattern implies that the garnets in deep SK Ⅰ may have formed under oxidizing conditions with high temperatures, high oxygen fugacity (fO2) and weak acidity, while the garnets in shallow SK Ⅲ formed under reducing conditions with moderate temperatures, moderate fO2 and near-neutral pH. The garnets in shallow SK Ⅰ and deep SK Ⅲ may have formed under weak oxidation–reduction conditions at moderate to high temperatures, moderate fO2 and weakly acidic to neutral pH conditions. Our newly reported field geological, mineralogical, and geochemical data from garnet suggest that the ore-forming fluids of the Huanggangliang were magmatic-hydrothermal fluids mixed with slightly acidic external fluids, chloride complexes may have been the main transport ligands of REEs. The metallogenic environment may be an important factor leading to the difference of mineral assemblages in the Huanggangliang deposit.

Structure and mechanisms of transport of human Asc1/CD98hc amino acid transporter

Recent cryoEM studies elucidated details of the structural basis for the substrate selectivity and translocation of heteromeric amino acid transporters. However, Asc1/CD98hc is the only neutral heteromeric amino acid transporter that can function through facilitated diffusion, and the only one that efficiently transports glycine and D-serine, and thus has a regulatory role in the central nervous system. Here we use cryoEM, ligand-binding simulations, mutagenesis, transport assays, and molecular dynamics to define human Asc1/CD98hc determinants for substrate specificity and gain insights into the mechanisms that govern substrate translocation by exchange and facilitated diffusion. The cryoEM structure of Asc1/CD98hc is determined at 3.4–3.8 Å resolution, revealing an inward-facing semi-occluded conformation. We find that Ser 246 and Tyr 333 are essential for Asc1/CD98hc substrate selectivity and for the exchange and facilitated diffusion modes of transport. Taken together, these results reveal the structural bases for ligand binding and transport features specific to human Asc1.

Hedgehog on the Move: Glypican-Regulated Transport and Gradient Formation in Drosophila.

Glypicans (Glps) are a family of heparan sulphate proteoglycans that are attached to the outer plasma membrane leaflet of the producing cell by a glycosylphosphatidylinositol anchor. Glps are involved in the regulation of many signalling pathways, including those that regulate the activities of Wnts, Hedgehog (Hh), Fibroblast Growth Factors (FGFs), and Bone Morphogenetic Proteins (BMPs), among others. In the Hh-signalling pathway, Glps have been shown to be essential for ligand transport and the formation of Hh gradients over long distances, for the maintenance of Hh levels in the extracellular matrix, and for unimpaired ligand reception in distant recipient cells. Recently, two mechanistic models have been proposed to explain how Hh can form the signalling gradient and how Glps may contribute to it. In this review, we describe the structure, biochemistry, and metabolism of Glps and their interactions with different components of the Hh-signalling pathway that are important for the release, transport, and reception of Hh.

Molecular insights into ferric-siderophore transport by the putative TonB-dependent transporter in Mycobacterium tuberculosis

Iron acquisition is critical to the virulence of Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis. To acquire iron within the host, Mtb secretes siderophores that chelate iron with high affinity. Siderophores scavenge iron from host cells using TonB-dependent transporters like FecA. We investigated molecular mechanisms of FecA-mediated ferric-siderophore transport in Mtb. Molecular docking and molecular dynamics simulations revealed a series of interactions between ferric siderophores and FecA. The initial binding occurs at a pocket located on the extracellular surface of FecA. The ligand then migrates deeper through the transport tunnel to a subsequent binding site, aided by conformational changes in FecA that expand the tunnel diameter. We observed the key roles of precise positioning of extracellular loops in the outer membrane barrel and plug domains in the optimal ligand binding and transport. Transport of ferric–siderophore complex into Mtb follows an induced fit model, with ligand interaction eliciting 2–10 Å shifts in the barrel and plug regions. By revealing the conformational dynamics enabling iron import, these findings provide molecular-level insights into a metal ion uptake mechanism in Mtb. Iron acquisition is essential for Mtb pathogenesis, so this work may inform novel therapeutic strategies that disrupt siderophore uptake pathways. Communicated by Ramaswamy H. Sarma

Fatty-Acid-Binding Proteins: From Lipid Transporters to Disease Biomarkers.

Fatty-acid-binding proteins (FABPs) serve a crucial role in the metabolism and transport of fatty acids and other hydrophobic ligands as an intracellular protein family. They are also recognized as a critical mediator in the inflammatory and ischemic pathways. FABPs are found in a wide range of tissues and organs, allowing them to contribute to various disease/injury developments that have not been widely discussed. We have collected and analyzed research journals that have investigated the role of FABPs in various diseases. Through this review, we discuss the findings on the potential of FABPs as biomarkers for various diseases in different tissues and organs, looking at their expression levels and their roles in related diseases according to available literature data. FABPs have been reported to show significantly increased expression levels in various tissues and organs associated with metabolic and inflammatory diseases. Therefore, FABPs are a promising novel biomarker that needs further development to optimize disease diagnosis and prognosis methods along with previously discovered markers.

Molecular Insight into Ligand Binding and Transport by the Lentil Lipid Transfer Protein Lc-LTP2: The Role of Basic Amino Acid Residues at Opposite Entrances to the Hydrophobic Cavity.

Lipid transfer proteins (LTPs) realize their functions in plants due to their ability to bind and transport various ligands. Structures of many LTPs have been studied; however, the mechanism of ligand binding and transport is still not fully understood. In this work, we studied the role of Lys61 and Lys81 located near the "top" and "bottom" entrances to the hydrophobic cavity of the lentil lipid transfer protein Lc-LTP2, respectively, in these processes. Using site-directed mutagenesis, we showed that both amino acid residues played a key role in lipid binding to the protein. In experiments with calcein-loaded liposomes, we demonstrated that both the above-mentioned lysine residues participated in the protein interaction with model membranes. According to data obtained from fluorescent spectroscopy and TNS probe displacement, both amino acid residues are necessary for the ability of the protein to transfer lipids between membranes. Thus, we hypothesized that basic amino acid residues located at opposite entrances to the hydrophobic cavity of the lentil Lc-LTP2 played an important role in initial protein-ligand interaction in solution as well as in protein-membrane docking.

Nanobody-based VSR7 tracing shows clathrin-dependent TGN to Golgi recycling

Receptor-mediated transport of soluble proteins is nature’s key to empowering eukaryotic cells to access a plethora of macromolecules, either by direct accumulation or as products from resulting biochemical pathways. The transport efficiency of these mechanisms results from the receptor’s capability to capture, transport, and release ligands on the one hand and the cycling ability that allows for performing multiple rounds of ligand transport on the other. However, the plant VACUOLAR SORTING RECEPTOR (VSR) protein family is diverse, and their ligand-specificity and bidirectional trafficking routes and transport mechanisms remain highly controversial. Here we employ nanobody-epitope interaction-based molecular tools to assess the function of the VSR 7 in vivo. We demonstrate the specificity of the VSR7 for sequence-specific vacuolar sorting signals, and we trace its anterograde transport and retrograde recycling route. VSR7 localizes at the cis-Golgi apparatus at steady state conditions and transports ligands downstream to release them in the trans-Golgi network/early endosome (TGN/EE) before undergoing clathrin-dependent recycling from the TGN/EE back to the cis-Golgi.

Phylogenetic Analysis Guides Transporter Protein Deorphanization: A Case Study of the SLC25 Family of Mitochondrial Metabolite Transporters.

Homology search and phylogenetic analysis have commonly been used to annotate gene function, although they are prone to error. We hypothesize that the power of homology search in functional annotation depends on the coupling of sequence variation to functional diversification, and we herein focus on the SoLute Carrier (SLC25) family of mitochondrial metabolite transporters to survey this coupling in a family-wide manner. The SLC25 family is the largest family of mitochondrial metabolite transporters in eukaryotes that translocate ligands of different chemical properties, ranging from nucleotides, amino acids, carboxylic acids and cofactors, presenting adequate experimentally validated functional diversification in ligand transport. Here, we combine phylogenetic analysis to profile SLC25 transporters across common eukaryotic model organisms, from Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Danio rerio, to Homo sapiens, and assess their sequence adaptations to the transported ligands within individual subfamilies. Using several recently studied and poorly characterized SLC25 transporters, we discuss the potentials and limitations of phylogenetic analysis in guiding functional characterization.

PyCaverDock: Python implementation of the popular tool for analysis of ligand transport with advanced caching and batch calculation support.

Access pathways in enzymes are crucial for the passage of substrates and products of catalysed reactions. The process can be studied by computational means with variable degrees of precision. Our in-house approximative method CaverDock provides a fast and easy way to set up and run ligand binding and unbinding calculations through protein tunnels and channels. Here we introduce pyCaverDock, a Python3 API designed to improve user experience with the tool and further facilitate the ligand transport analyses. The API enables users to simplify the steps needed to use CaverDock, from automatizing setup processes to designing screening pipelines. pyCaverDock API is implemented in Python 3 and is freely available with detailed documentation and practical examples at https://loschmidt.chemi.muni.cz/caverdock/.

Interaction at a distance: Xenon migration in Mb

The transport of ligands, such as NO or O2, through internal cavities is essential for the function of globular proteins, including hemoglobin, myoglobin (Mb), neuroglobin, truncated hemoglobins, or cytoglobin. For Mb, several internal cavities (Xe1 through Xe4) were observed experimentally and they were linked to ligand storage. The present work determines barriers for xenon diffusion and relative stabilization energies for the ligand in the initial and final pocket, linking a transition depending on the occupancy state of the remaining pockets from both biased and unbiased molecular dynamics simulations. It is found that the energetics of a particular ligand migration pathway may depend on the direction in which the transition is followed and the occupancy state of the other cavities. Furthermore, the barrier height for a particular transition can depend in a non-additive fashion on the occupancy of either cavity A or B or simultaneous population of both cavities, A and B. Multiple repeats for the Xe1 → Xe2 transition reveal that the activation barrier is a distribution of barrier heights rather than one single value, which is confirmed by a distribution of transition times for the same transition from unbiased simulations. Dynamic cross correlation maps demonstrate that correlated motions occur between adjacent residues or through space, residue Phe138 is found to be a gate for the Xe1 → Xe2 transition, and the volumes of the internal cavities vary along the diffusion pathway, indicating that there is dynamic communication between the ligand and the protein. These findings suggest that Mb is an allosteric protein.

In SilicoIdentification of Triclosan Derivatives as Potential Inhibitors of MutantMycobacterium tuberculosisInhA

Enoyl acyl carrier protein reductase (InhA) is a crucial enzyme for the biosynthesis of mycolic acids which are major compartments of the Mycobacterium tuberculosis (Mtb) cell wall. Direct inhibition of InhA without activation by drug-NADH adduct has clinical utility to overcome drug resistance. We aimed at the in silico identification of triclosan derivatives with the potential inhibitory effect of S94A-InhA as a clinically important mutant form. Caver Web 1.0 server was used to analyze the ligand transport through access tunnels. Two macrocyclic triclosan derivatives ( 4 and 6) could be identified with more energy-favorable transfer routes toward the enzyme active site. Molecular dynamics (MD) simulations (50 ns) of the best-scored compounds revealed the stability of enzyme structure upon binding to 4 and 6. Compound 4 could better retain enzyme stability upon target binding. Results of intermolecular H-bond analysis indicated that both complexes were mediated through hydrophobic contacts. Declined solvent accessible surface area (SASA) for the apo and bound enzyme states indicated non-destabilizing behavior and no structural relaxation. Electrostatic and van der Waals interactions between triclosan derivatives and their surroundings were used to acquire binding free energies through the linear interaction energy (LIE) method based on MD simulations (Average [Formula: see text], [Formula: see text] kcal/mol and [Formula: see text] kcal/mol). Both of the triclosan derivatives showed relatively stable energy variations and their steady accommodation inside enzyme active site could be confirmed during 50 ns. These results may be implicated in further structure-guided approaches against drug-resistant Mtb.