Evaluating the Contribution of Transition-State Destabilization to Changes in the Residence Time of Triazole-Based InhA Inhibitors

Accurate Characterization of Binding Kinetics and Allosteric Mechanisms for the HSP90 Chaperone Inhibitors Using AI-Augmented Integrative Biophysical Studies

The binding kinetics of drugs to their targets are gradually being recognized as a crucial indicator of the efficacy of drugs in vivo, leading to the development of various computational methods for predicting the binding kinetics in recent years. However, compared with the prediction of binding affinity, the underlying structure and dynamic determinants of binding kinetics are more complicated. Efficient and accurate methods for predicting binding kinetics are still lacking. In this study, quantitative structure-kinetics relationship (QSKR) models were developed using 132 inhibitors targeting the ATP binding domain of heat shock protein 90α (HSP90α) to predict the dissociation rate constant (koff), enabling a direct assessment of the drug-target residence time. These models demonstrated good predictive performance, where hydrophobic and hydrogen bond interactions significantly influence the koff prediction. In subsequent applications, our models were used to assist in the discovery of new inhibitors for the N-terminal domain of HSP90α (N-HSP90α), demonstrating predictive capabilities on an experimental validation set with a new scaffold. In X-ray crystallography experiments, the loop-middle conformation of apo N-HSP90α was observed for the first time (previously, the loop-middle conformation had only been observed in holo-N-HSP90α structures). Interestingly, we observed different conformations of apo N-HSP90α simultaneously in an asymmetric unit, which was also observed in a holo-N-HSP90α structure, suggesting an equilibrium of conformations between different states in solution, which could be one of the determinants affecting the binding kinetics of the ligand. Different ligands can undergo conformational selection or alter the equilibrium of conformations, inducing conformational rearrangements and resulting in different effects on binding kinetics. We then used molecular dynamics simulations to describe conformational changes of apo N-HSP90α in different conformational states. In summary, the study of the binding kinetics and molecular mechanisms of N-HSP90α provides valuable information for the development of more targeted therapeutic approaches.

Discovery of new diaryl ether inhibitors against Mycobacterium tuberculosis targeting the minor portal of InhA

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) affects 10 million people each year and the emergence of resistant TB augurs for a growing incidence. In the last 60 years, only three new drugs were approved for TB treatment, for which resistances are already emerging. Therefore, there is a crucial need for new chemotherapeutic agents capable of eradicating TB. Enzymes belonging to the type II fatty acid synthase system (FAS-II) are involved in the biosynthesis of mycolic acids, cell envelope components essential for mycobacterial survival. Among them, InhA is the primary target of isoniazid (INH), one of the most effective compounds to treat TB. INH acts as a prodrug requiring activation by the catalase-peroxidase KatG, whose mutations are the major cause for INH resistance. Herein, a new series of direct InhA inhibitors were designed based on a molecular hybridization approach. They exhibit potent inhibitory activities of InhA and, for some of them, good antitubercular activities. Moreover, they display a low toxicity on human cells. A study of the mechanism of action of the most effective molecules shows that they inhibit the biosynthesis of mycolic acids. The X-ray structures of two InhA/NAD+/inhibitor complexes have been obtained showing a binding mode of a part of the molecule in the minor portal, rarely seen in the InhA structures reported so far.

The translational value of ligand-receptor binding kinetics in drug discovery

The translation of in vitro potency of a candidate drug, as determined by traditional pharmacology metrics (such as EC50 /IC50 and KD /Ki values), to in vivo efficacy and safety is challenging. Residence time, which represents the duration of drug-target interaction, can be part of a more comprehensive understanding of the dynamic nature of drug-target interactions in vivo, thereby enabling better prediction of drug efficacy and safety. As a consequence, a prolonged residence time may help in achieving sustained pharmacological activity, while transient interactions with shorter residence times may be favourable for targets associated with side effects. Therefore, integration of residence time into the early stages of drug discovery and development has yielded a number of clinical candidates with promising in vivo efficacy and safety profiles. Insights from residence time research thus contribute to the translation of in vitro potency to in vivo efficacy and safety. Further research and advances in measuring and optimizing residence time will bring a much-needed addition to the drug discovery process and the development of safer and more effective drugs. In this review, we summarize recent research progress on residence time, highlighting its importance from a translational perspective.

Evaluating the Impact of the Tyr158 pKa on the Mechanism and Inhibition of InhA, the Enoyl-ACP Reductase from Mycobacterium tuberculosis

InhA, the Mycobacterium tuberculosis enoyl-ACP reductase, is a target for the tuberculosis (TB) drug isoniazid (INH). InhA inhibitors that do not require KatG activation avoid the most common mechanism of INH resistance, and there are continuing efforts to fully elucidate the enzyme mechanism to drive inhibitor discovery. InhA is a member of the short-chain dehydrogenase/reductase superfamily characterized by a conserved active site Tyr, Y158 in InhA. To explore the role of Y158 in the InhA mechanism, this residue has been replaced by fluoroTyr residues that increase the acidity of Y158 up to ∼3200-fold. Replacement of Y158 with 3-fluoroTyr (3-FY) and 3,5-difluoroTyr (3,5-F2Y) has no effect on kcatapp/KMapp nor on the binding of inhibitors to the open form of the enzyme (Kiapp), whereas both kcatapp/KMapp and Kiapp are altered by seven-fold for the 2,3,5-trifluoroTyr variant (2,3,5-F3Y158 InhA). 19F NMR spectroscopy suggests that 2,3,5-F3Y158 is ionized at neutral pH indicating that neither the acidity nor ionization state of residue 158 has a major impact on catalysis or on the binding of substrate-like inhibitors. In contrast, Ki*app is decreased 6- and 35-fold for the binding of the slow-onset inhibitor PT504 to 3,5-F2Y158 and 2,3,5-F3Y158 InhA, respectively, indicating that Y158 stabilizes the closed form of the enzyme adopted by EI*. The residence time of PT504 is reduced ∼four-fold for 2,3,5-F3Y158 InhA compared to wild-type, and thus, the hydrogen bonding interaction of the inhibitor with Y158 is an important factor in the design of InhA inhibitors with increased residence times on the enzyme.

The pathogenic mechanism of Mycobacterium tuberculosis: implication for new drug development

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), is a tenacious pathogen that has latently infected one third of the world's population. However, conventional TB treatment regimens are no longer sufficient to tackle the growing threat of drug resistance, stimulating the development of innovative anti-tuberculosis agents, with special emphasis on new protein targets. The Mtb genome encodes ~4000 predicted proteins, among which many enzymes participate in various cellular metabolisms. For example, more than 200 proteins are involved in fatty acid biosynthesis, which assists in the construction of the cell envelope, and is closely related to the pathogenesis and resistance of mycobacteria. Here we review several essential enzymes responsible for fatty acid and nucleotide biosynthesis, cellular metabolism of lipids or amino acids, energy utilization, and metal uptake. These include InhA, MmpL3, MmaA4, PcaA, CmaA1, CmaA2, isocitrate lyases (ICLs), pantothenate synthase (PS), Lysine-ε amino transferase (LAT), LeuD, IdeR, KatG, Rv1098c, and PyrG. In addition, we summarize the role of the transcriptional regulator PhoP which may regulate the expression of more than 110 genes, and the essential biosynthesis enzyme glutamine synthetase (GlnA1). All these enzymes are either validated drug targets or promising target candidates, with drugs targeting ICLs and LAT expected to solve the problem of persistent TB infection. To better understand how anti-tuberculosis drugs act on these proteins, their structures and the structure-based drug/inhibitor designs are discussed. Overall, this investigation should provide guidance and support for current and future pharmaceutical development efforts against mycobacterial pathogenesis.

Structure-Kinetic Relationship Studies for the Development of Long Residence Time LpxC Inhibitors

UDP-3-O-(R-3-hydroxymyristoyl)-N-acetylglucosamine deacetylase (LpxC) is a promising drug target in Gram-negative bacteria. Previously, we described a correlation between the residence time of inhibitors on Pseudomonas aeruginosa LpxC (paLpxC) and the post-antibiotic effect (PAE) caused by the inhibitors on the growth of P. aeruginosa. Given that drugs with prolonged activity following compound removal may have advantages in dosing regimens, we have explored the structure-kinetic relationship for paLpxC inhibition by analogues of the pyridone methylsulfone PF5081090 (1) originally developed by Pfizer. Several analogues have longer residence times on paLpxC than 1 (41 min) including PT913, which has a residence time of 124 min. PT913 also has a PAE of 4 h, extending the original correlation observed between residence time and PAE. Collectively, the studies provide a platform for the rational modulation of paLpxC inhibitor residence time and the potential development of antibacterial agents that cause prolonged suppression of bacterial growth.

Binding Thermodynamics and Dissociation Kinetics Analysis Uncover the Key Structural Motifs of Phenoxyphenol Derivatives as the Direct InhA Inhibitors and the Hotspot Residues of InhA

Given the current epidemic of multidrug-resistant tuberculosis, there is an urgent need to develop new drugs to combat drug-resistant tuberculosis. Direct inhibitors of the InhA target do not require activation and thus can overcome drug resistance caused by mutations in drug-activating enzymes. In this work, the binding thermodynamic and kinetic information of InhA to its direct inhibitors, phenoxyphenol derivatives, were explored through multiple computer-aided drug design (CADD) strategies. The results show that the van der Waals interactions were the main driving force for protein–ligand binding, among which hydrophobic residues such as Tyr158, Phe149, Met199 and Ile202 have high energy contribution. The AHRR pharmacophore model generated by multiple ligands demonstrated that phenoxyphenol derivatives inhibitors can form pi–pi stacking and hydrophobic interactions with InhA target. In addition, the order of residence time predicted by random acceleration molecular dynamics was consistent with the experimental values. The intermediate states of these inhibitors could form hydrogen bonds and van der Waals interactions with surrounding residues during dissociation. Overall, the binding and dissociation mechanisms at the atomic level obtained in this work can provide important theoretical guidance for the development of InhA direct inhibitors with higher activity and proper residence time.

Discovery of novel and potent InhA inhibitors by an in silico screening and pharmacokinetic prediction

Aim: In silico screening approaches were performed to discover novel InhA inhibitors. Methods: Candidate InhA inhibitors were obtained from the combination of virtual screening and pharmacokinetic prediction. In addition, molecular mechanics Poisson-Boltzmann surface area, molecular mechanics Generalized Born surface area and WaterSwap methods were performed to investigate the binding interactions and binding energy of candidate compounds. Results: Four candidate compounds with suitable physicochemical, pharmacokinetic and antibacterial properties are proposed. The crucial interactions of the candidate compounds were H-bond, pi-pi and sigma-pi interactions observed in the InhA binding site. The binding affinity of these compounds was improved by hydrophobic interactions with hydrophobic side chains in the InhA pocket. Conclusion: The four newly identified InhA inhibitors reported in this study could serve as promising hit compounds against Mycobacterium tuberculosis and may be considered for further experimental studies.

Biochemical and structural basis for Moraxella catarrhalis enoyl-acyl carrier protein reductase (FabI) inhibition by triclosan and estradiol

The enoyl-acyl carrier protein reductase (ENR) is an established drug target and catalyzes the last reduction step of the fatty acid elongation cycle. Here, we report the crystal structures of FabI from Moraxella catarrhalis (McFabI) in the apo form, binary complex with NAD+ and ternary complex with NAD ⁺ -triclosan (TCL) determined at 2.36, 2.12 and 2.22 Å resolutions, respectively. The comparative study of these three structures revealed three different conformational states for the substrate-binding loop (SBL), including an unstructured intermediate, a structured intermediate and a closed conformation in the apo, binary and ternary complex forms, respectively; indicating the flexibility of SBL during the ligand binding. Virtual screening has suggested that estradiol cypionate may be a potential inhibitor of McFabI. Subsequently, estradiol (EST), the natural form of estradiol cypionate, was assessed for its FabI-binding and -inhibition properties. In vitro studies demonstrated that TCL and EST bind to McFabI with high affinity (K_D = 0.038 ± 0.004 and 5 ± 0.06 μM respectively) and inhibit its activity (K_i = 62.93 ± 3.95 nM and 25.97 ± 1.93 μM respectively) and suppress the growth of M. catarrhalis. These findings reveal that TCL and EST inhibit the McFabI activity and thereby affect cell growth. This study suggests that estradiol may be exploited as a novel scaffold for the designing and development of more potential FabI inhibitors.

Mycobacterium enoyl acyl carrier protein reductase (InhA): A key target for antitubercular drug discovery

Enoyl acyl carrier protein reductase (InhA) is a key enzyme involved in fatty acid synthesis mainly mycolic acid biosynthesis that is a part of NADH dependent acyl carrier protein reductase family. The aim of the present literature is to underline the different scaffolds or enzyme inhibitors that inhibit mycolic acid biosynthesis mainly cell wall synthesis by inhibiting enzyme InhA. Various scaffolds were identified based on the screening technologies like high throughput screening, encoded library technology, fragment-based screening. The compounds studied include indirect inhibitors (Isoniazid, Ethionamide, Prothionamide) and direct inhibitors (Triclosan/Diphenyl ethers, Pyrrolidine Carboxamides, Pyrroles, Acetamides, Thiadiazoles, Triazoles) with better efficacy against drug resistance. Out of the several scaffolds studied, pyrrolidine carboxamides were found to be the best molecules targeting InhA having good bioavailability properties and better MIC. This review provides with a detailed information, analysis, structure activity relationship and useful insight on various scaffolds as InhA inhibitors.

Evolution of the drug-target residence time model

IntroductionThe pharmacological action of a drug is linked to its affinity for a specific molecular target as quantified by in vitro equilibrium measurements. However, it is clear that for many highly effective drugs, interactions with their molecular targets do not conform to simple, equilibrium conditions in vivo and this results in a temporal discordance between pharmacokinetics and pharmacodynamics. The drug-target residence time model was developed to provide a theoretical framework with which to understand cases in which very slow dissociation of the drug-target complex in vivo results in durable PD effects even after systemic concentrations of drug have waned.Area coveredIn this article, the author provides a brief description of the drug-target residence time model and focuses on the refinements that have been made to the original model to incorporate the influences of compound rebinding in cells and pharmacokinetic properties of drug molecules.Expert opinionThere is now overwhelming evidence for the utility of the drug-target residence time model as a framework for understanding in vivo drug action. The in vitro measured residence time (τ_R) must be used in concert with equilibrium measures of drug-target affinity (e.g. IC₅₀) and with in vivo measures of pharmacokinetic half-life, to afford the researcher a powerful approach to compound optimization for clinical effect. Despite the significant use and refinement of this model, continued studies are required to better understand the dynamic interplay between residence time, target pathobiology, drug distribution and drug pharmacokinetics.

Anti-tubercular activity and molecular docking studies of indolizine derivatives targeting mycobacterial InhA enzyme

A series of 1,2,3-trisubstituted indolizines (2a-2f, 3a-3d, and 4a-4c) were screened for in vitro whole-cell anti-tubercular activity against the susceptible H37Rv and multidrug-resistant (MDR) Mycobacterium tuberculosis (MTB) strains. Compounds 2b-2d, 3a-3d, and 4a-4c were active against the H37Rv-MTB strain with minimum inhibitory concentration (MIC) ranging from 4 to 32 µg/mL, whereas the indolizines 4a-4c, with ethyl ester group at the 4-position of the benzoyl ring also exhibited anti-MDR-MTB activity (MIC = 16-64 µg/mL). In silico docking study revealed the enoyl-acyl carrier protein reductase (InhA) and anthranilate phosphoribosyltransferase as potential molecular targets for the indolizines. The X-ray diffraction analysis of the compound 4b was also carried out. Further, a safety study (in silico and in vitro) demonstrated no toxicity for these compounds. Thus, the indolizines warrant further development and may represent a novel promising class of InhA inhibitors and multi-targeting agents to combat drug-sensitive and drug-resistant MTB strains.

A Long Residence Time Enoyl-Reductase Inhibitor Explores an Extended Binding Region with Isoenzyme-Dependent Tautomer Adaptation and Differential Substrate-Binding Loop Closure

The enoyl-acyl carrier protein (ACP) reductase (ENR) is a key enzyme within the bacterial fatty-acid synthesis pathway. It has been demonstrated that small-molecule inhibitors carrying the diphenylether (DPE) scaffold bear a great potential for the development of highly specific and effective drugs against this enzyme class. Interestingly, different substitution patterns of the DPE scaffold have been shown to lead to varying effects on the kinetic and thermodynamic behavior toward ENRs from different organisms. Here, we investigated the effect of a 4'-pyridone substituent in the context of the slow tight-binding inhibitor SKTS1 on the inhibition of the Staphylococcus aureus enoyl-ACP-reductase saFabI and the closely related isoenzyme from Mycobacterium tuberculosis, InhA, and explored a new interaction site of DPE inhibitors within the substrate-binding pocket. Using high-resolution crystal structures of both complexes in combination with molecular dynamics (MD) simulations, kinetic measurements, and quantum mechanical (QM) calculations, we provide evidence that the 4'-pyridone substituent adopts different tautomeric forms when bound to the two ENRs. We furthermore elucidate the structural determinants leading to significant differences in the residence time of SKTS1 on both enzymes.

Design, synthesis, in silico, and in vitro evaluation of 3-phenylpyrazole acetamide derivatives as antimycobacterial agents

Mycobacterium tuberculosis (Mtb) is one of the most dangerous pathogens affecting immunocompetent and immunocompromised patients worldwide. Novel molecules, which are efficient and can reduce the duration of therapy against drug-resistant strains, are an urgent unmet need of the hour. In our current study, a series of new 2-(3-phenyl-1H-pyrazol-1-yl)acetamide and N'-benzylidene-2-(3-phenyl-1H-pyrazol-1-yl)acetohydrazide derivatives were designed, synthesized, and evaluated for their antimycobacterial potential. The biological evaluation revealed that 6b, 6m, 6l, 7a, and 7k exhibited selective and potent inhibitory activity against Mtb. Furthermore, compounds 6m and 7h were found to be nontoxic to Vero cells with CC₅₀ of greater than 20 and 80 mg/ml, respectively, and exhibited promising selectivity indices (SI) of greater than 666 and 320, respectively. All derivatives exhibited excellent ADME (absorption, distribution, metabolism, and excretion) properties in silico. Also, all the derivatives were found compliant with Lipinski's rule of five, showing their druggability profile. Molecular docking insights of these derivatives have shown outstanding binding energies on the mycobacterial membrane protein large transporters. These results indicate that this scaffold may lead to a potential antimycobacterial drug candidate in the discovery of antitubercular agents.

Exploring the chemical space of 1,2,3-triazolyl triclosan analogs for discovery of new antileishmanial chemotherapeutic agents

Triclosan and isoniazid are known antitubercular compounds that have proven to be also active against Leishmania parasites. On these grounds, a collection of 37 diverse 1,2,3-triazoles based on the antitubercular molecules triclosan and 5-octyl-2-phenoxyphenol (8PP) were designed in search of novel structures with leishmanicidal activity and prepared using different alkynes and azides. The 37 compounds were assayed against Leishmania donovani, the etiological agent of leishmaniasis, yielding some analogs with activity at micromolar concentrations and against M. tuberculosis H37Rv resulting in scarce active compounds with an MIC of 20 μM. To study the mechanism of action of these catechols, we analyzed the inhibition activity of the library on the M. tuberculosis enoyl-ACP reductase (ENR) InhA, obtaining poor inhibition of the enzyme. The cytotoxicity against Vero cells was also tested, resulting in none of the compounds being cytotoxic at concentrations of up to 20 μM. Derivative 5f could be considered a valuable starting point for future antileishmanial drug development. The validation of a putative leishmanial InhA orthologue as a therapeutic target needs to be further investigated.

A structure-kinetic relationship study using matched molecular pair analysis

The lifetime of a binary drug–target complex is increasingly acknowledged as an important parameter for drug efficacy and safety. With a better understanding of binding kinetics and better knowledge about kinetic parameter optimization, intentionally induced prolongation of the drug–target residence time through structural changes of the ligand could become feasible. In this study we assembled datasets from 21 publications and the K4DD (Kinetic for Drug Discovery) database to conduct large scale data analysis. This resulted in 3812 small molecules annotated to 78 different targets from five protein classes (GPCRs: 273, kinases: 3238, other enzymes: 240, HSPs: 160, ion channels: 45). Performing matched molecular pair (MMP) analysis to further investigate the structure–kinetic relationship (SKR) in this data collection allowed us to identify a fundamental contribution of a ligand's polarity to its association rate, and in selected cases, also to its dissociation rate. However, we furthermore observed that the destabilization of the transition state introduced by increased polarity is often accompanied by simultaneous destabilization of the ground state resulting in an unaffected or even worsened residence time. Supported by a set of case studies, we provide concepts on how to alter ligands in ways to trigger on-rates, off-rates, or both.

A large-scale study employing matched molecular pair (MMP) analysis to uncover the contribution of a compound's polarity to its association and dissociation rates.

Efficient fungal and bacterial facilitated remediation of thiencarbazone methyl in the environment

Triazole herbicide, Thiencarbazone-methyl (TCM) applied on different crops for weedicidal activity is associated with an inherent toxicity towards bladder and urinary functionality. TCM has been first time explored for its biodegradative behavior utilizing microbes, previously isolated from soils. Simulated bio-transformation assemblies of five fungal strains; Aspergillus flavus (AF), Penicillium chrysogenum (PC), Aspergillus niger (AN), Aspergillus terrus (AT), Aspergillus fumigatus (AFu) and two bacterial strains: Xanthomonas citri (XC), Pseudomonassyringae (PS), were utilized. 10 mg/L TCM concentration was set up utilizing each microbe and analysed for 42 days. TCM bio-degradation was evaluated by UV-Visible spectrophotometery and gas chromatography mass spectroscopy. Aspergillus terrus (R² = 0.86) and Penicillium chrysogenum (R² = 0.88) exhibited highest capability to metabolize TCM while forming intermediate metabolites including; 2,4-dihydro-[1,2,4] triazol-3-one, semicarbazide and urea, methyl 4-isocyanatosulfonyl-5-methylthiophene-3-carboxylate. TCM degradation by all strains AF, AFu, AN, PC, AT, PS and XC was found to be 74, 74, 81, 95, 98, 90 and 95%, respectively after 42 days elucidating the effectiveness of all the utilized strains in degrading TCM. Current investigations can impact vital bioremediation approaches for pesticides mitigation from the ecological compartments. Furthermore, present research can be extended to the optimization of the bio-deteriorative assays to be employed on the practical scale for the successful management of environment through sustainable and cost effective ways.

FabI (enoyl acyl carrier protein reductase) - A potential broad spectrum therapeutic target and its inhibitors

Development of new anti-bacterial agents acting upon underexploited targets and thus evading known mechanisms of resistance is the need of the hour. The highly conserved and distinct bacterial fatty acid biosynthesis pathway (FAS-II), presents a validated and yet relatively underexploited target for drug discovery. FabI and its isoforms (FabL, FabK, FabV and InhA) are essential enoyl-ACP reductases present in several microorganisms. In addition, the components of the FAS-II pathway are distinct from the multi-enzyme FAS-I complex found in mammals. Thus, inhibition of FabI and its isoforms is anticipated to result in broad-spectrum antibacterial activity. Several research groups from industry and academic laboratories have devoted significant efforts to develop effective FabI-targeting antibiotics, which are currently in various stages of clinical development for the treatment of multi-drug resistant bacterial infections. This review summarizes all the natural as well as synthetic inhibitors of gram-positive and gram-negative enoyl ACP reductases (FabI). The knowledge of the reported inhibitors can aid in the development of broad-spectrum antibacterials specifically targeting FabI enzymes from S. aureus, S. epidermidis, B. anthracis, B. cereus, E. coli, P. aeruginosa, P. falciparum and M. tuberculosis.

On Calculating Free Energy Differences Using Ensembles of Transition Paths

The free energy of a process is the fundamental quantity that determines its spontaneity or propensity at a given temperature. In particular, the binding free energy of a drug candidate to its biomolecular target is used as an objective quantity in drug design. Recently, binding kinetics—rates of association (k_on) and dissociation (k_off)—have also demonstrated utility for their ability to predict efficacy and in some cases have been shown to be more predictive than the binding free energy alone. Some methods exist to calculate binding kinetics from molecular simulations, although these are typically more difficult to calculate than the binding affinity as they depend on details of the transition path ensemble. Assessing these rate constants can be difficult, due to uncertainty in the definition of the bound and unbound states, large error bars and the lack of experimental data. As an additional consistency check, rate constants from simulation can be used to calculate free energies (using the log of their ratio) which can then be compared to free energies obtained experimentally or using alchemical free energy perturbation. However, in this calculation it is not straightforward to account for common, practical details such as the finite simulation volume or the particular definition of the “bound” and “unbound” states. Here we derive a set of correction terms that can be applied to calculations of binding free energies using full reactive trajectories. We apply these correction terms to revisit the calculation of binding free energies from rate constants for a host-guest system that was part of a blind prediction challenge, where significant deviations were observed between free energies calculated with rate ratios and those calculated from alchemical perturbation. The correction terms combine to significantly decrease the error with respect to computational benchmarks, from 3.4 to 0.76 kcal/mol. Although these terms were derived with weighted ensemble simulations in mind, some of the correction terms are generally applicable to free energies calculated using physical pathways via methods such as Markov state modeling, metadynamics, milestoning, or umbrella sampling.

Hydrolysis of Amisulbrom in Buffer Solutions and Natural Water Samples: Kinetics and Products Identification

In this study, the hydrolysis of amisulbrom in buffer solutions and natural water samples were investigated. Effects of pH and temperature were tested in buffer solutions. Amisulbrom was stable in acidic and neutral aqueous solutions at 25°C, while quickly hydrolyzed with a half-life of 4.5 days (25°C) at pH 9.0. The kinetics rate equation was determined as k = 1.0234 × 10¹⁰ exp (-61.3760/R·T) (R² = 0.9642) for hydrolysis of amisulbrom at pH 9.0. The pH, ionic strength, and solubility were important factors influencing the hydrolysis of amisulbrom in natural water samples. Furthermore, three hydrolysis products were separated and identified in buffer solution (pH 9.0) and natural water samples. A tentative transformation mechanism of amisulbrom was proposed to rationalize the formation of HPs (hydrolysis products) based on their structural identification, DFT (density functional theory), and hydrolysis profiles. Toxicity prediction using the quantitative structure-activity relationship model revealed that the HP-I, and HP-II were more toxic than the parent amisulbrom. This investigation was the first to evaluate the behavior of amisulbrom hydrolysis in aquatic systems.

First triclosan-based macrocyclic inhibitors of InhA enzyme

Two macrocyclic derivatives based on the triclosan frame were designed and synthesized as inhibitors of Mycobacterium tuberculosis InhA enzyme. One of the two molecules M02 displayed promising inhibitory activity against InhA enzyme with an IC₅₀ of 4.7 μM. Molecular docking studies of these two compounds were performed and confirmed that M02 was more efficient as inhibitor of InhA activity. These molecules are the first macrocyclic direct inhibitors of InhA enzyme able to bind into the substrate pocket. Furthermore, these biaryl ether compounds exhibited antitubercular activities comparable to that of triclosan against M. tuberculosis H37Rv strain.

Overview of Recent Strategic Advances in Medicinal Chemistry

Introducing novel strategies, concepts, and technologies that speed up drug discovery and the drug development cycle is of great importance both in the highly competitive pharmaceutical industry as well as in academia. This Perspective aims to present a "big-picture" overview of recent strategic innovations in medicinal chemistry and drug discovery.

Ensemble-Based Steered Molecular Dynamics Predicts Relative Residence Time of A2A Receptor Binders

Drug-target residence time, the length of time for which a small molecule stays bound to its receptor target, has increasingly become a key property for optimization in drug discovery programs. However, its in silico prediction has proven difficult. Here we describe a method, using atomistic ensemble-based steered molecular dynamics (SMD), to observe the dissociation of ligands from their target G protein-coupled receptor in a time scale suitable for drug discovery. These dissociation simulations accurately, precisely, and reproducibly identify ligand-residue interactions and quantify the change in ligand energy values for both protein and water. The method has been applied to 17 ligands of the A_2A adenosine receptor, all with published experimental kinetic binding data. The residues that interact with the ligand as it dissociates are known experimentally to have an effect on binding affinities and residence times. There is a good correlation ( R² = 0.79) between the computationally calculated change in water-ligand interaction energy and experimentally determined residence time. Our results indicate that ensemble-based SMD is a rapid, novel, and accurate semi-empirical method for the determination of drug-target relative residence time.

Binding kinetics of ligands acting at GPCRs

The influence of drug-receptor binding kinetics has often been overlooked during the development of new therapeutics that target G protein-coupled receptors (GPCRs). Over the last decade there has been a growing understanding that an in-depth knowledge of binding kinetics at GPCRs is required to successfully target this class of proteins. Ligand binding to a GPCR is often not a simple single step process with ligand freely diffusing in solution. This review will discuss the experiments and equations that are commonly used to measure binding kinetics and how factors such as allosteric regulation, rebinding and ligand interaction with the plasma membrane may influence these measurements. We will then consider the molecular characteristics of a ligand and if these can be linked to association and dissociation rates.

Do new N-substituted 3-amino-4-phenyl-5-oxo-pyrazolinecarboxamide derivatives exhibit antitubercular potential?

As a continuation of previous tests concerning new N-substituted 3-amino-4-phenyl-5-oxo-pyrazolinecarboxamide derivatives (R3, R4 and R8) of notable antibacterial activity, their antitubercular potential against different mycobacterial strains was estimated. Tests performed on virulent (reference and clinical) strains of Mycobacterium bovis and Mycobacterium tuberculosis revealed the highest therapeutic potential of R8 derivative: MIC within the range 7.8-15.6 μg/ml and TI (therapeutic index) within the range 46.5-93. Moreover, the synergistic interaction was found between R3, R4 and R8 derivatives and rifampicin, one of the front-line antitubercular drugs. R8/rifampicin mixture in concentrations effective in inhibition of Mycobacterium tuberculosis strain was non-cytotoxic against GMK cells, displaying cell viability approximately 88-97% when compared to control. Molecular docking study enabled to conclude that enoyl acyl carrier protein reductase (InhA) can be considered as a potential molecular target of tested pyrazole derivatives. Although further modifications of chemical structure of the investigated pyrazole derivatives is required, in order to increase their antitubercular efficacy and therapeutic safety, these compounds, in particular R8 compound, can be promising for the treatment of human and bovine tuberculosis.

Design, synthesis, and biological evaluation of m-amidophenol derivatives as a new class of antitubercular agents

A series of m-amidophenol derivatives (6a-6l, 7a-7q, 9a, 9b, 12a-12c, 14 and 15) were designed and synthesized. Their antitubercular activities were evaluated in vitro against M. tuberculosis strains H37Ra and H37Rv and clinically isolated multidrug-resistant M. tuberculosis strains. Ten compounds displayed minimal inhibitory concentrations (MICs) against M. tuberculosis H37Ra below 2.5 μg mL-1 and 6g was the most active compound (MIC = 0.625 μg mL-1). Compounds 6g and 7a also showed potent inhibitory activity against M. tuberculosis H37Rv (MIC = 0.39 μg mL-1) and several clinically isolated multidrug-resistant M. tuberculosis strains (MIC = 0.39-3.125 μg mL-1). The compounds did not show inhibitory activity against normal Gram-positive and Gram-negative bacteria. They exhibited low cytotoxicity against HepG2 and RAW264.7 cell lines. The results demonstrated m-amidophenol as an attractive scaffold for the development of new antitubercular agents.

Structure-kinetic relationships that control the residence time of drug-target complexes: insights from molecular structure and dynamics

Time-dependent target occupancy is a function of both the thermodynamics and kinetics of drug-target interactions. However, while the optimization of thermodynamic affinity through approaches such as structure-based drug design is now relatively straight forward, less is understood about the molecular interactions that control the kinetics of drug complex formation and breakdown since this depends on both the ground and transition state energies on the binding reaction coordinate. In this opinion we highlight several recent examples that shed light on current approaches that are elucidating the factors that control the life-time of the drug-target complex.

Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery

Simple comparative correlation analyses and quantitative structure-kinetics relationship (QSKR) models highlight the interplay of kinetic rates and binding affinity as an essential feature in drug design and discovery. The choice of the molecular series, and their structural variations, used in QSKR modeling is fundamental to understanding the mechanistic implications of ligand and/or drug-target binding and/or unbinding processes. Here, we discuss the implications of linear correlations between kinetic rates and binding affinity constants and the relevance of the computational approaches to QSKR modeling.

Discovery of a cofactor-independent inhibitor of Mycobacterium tuberculosis InhA

AN12855 is a novel cofactor-independent inhibitor of Mycobacterium tuberculosis InhA. AN12855 has potent activity against M. tuberculosis, good oral bioavailability, and comparable efficacy to isoniazid in infection models.

New antitubercular agents are needed to combat the spread of multidrug- and extensively drug-resistant strains of Mycobacterium tuberculosis. The frontline antitubercular drug isoniazid (INH) targets the mycobacterial enoyl-ACP reductase, InhA. Resistance to INH is predominantly through mutations affecting the prodrug-activating enzyme KatG. Here, we report the identification of the diazaborines as a new class of direct InhA inhibitors. The lead compound, AN12855, exhibited in vitro bactericidal activity against replicating bacteria and was active against several drug-resistant clinical isolates. Biophysical and structural investigations revealed that AN12855 binds to and inhibits the substrate-binding site of InhA in a cofactor-independent manner. AN12855 showed good drug exposure after i.v. and oral delivery, with 53% oral bioavailability. Delivered orally, AN12855 exhibited dose-dependent efficacy in both an acute and chronic murine model of tuberculosis infection that was comparable with INH. Combined, AN12855 is a promising candidate for the development of new antitubercular agents.

Rational design of isonicotinic acid hydrazide derivatives with antitubercular activity: Machine learning, molecular docking, synthesis and biological testing

The problem of designing new antitubercular drugs against multiple drug-resistant tuberculosis (MDR-TB) was addressed using advanced machine learning methods. As there are only few published measurements against MDR-TB, we collected a large literature data set and developed models against the non-resistant H37Rv strain. The predictive accuracy of these models had a coefficient of determination q²  = .7-.8 (regression models) and balanced accuracies of about 80% (classification models) with cross-validation and independent test sets. The models were applied to screen a virtual chemical library, which was designed to have MDR-TB activity. The seven most promising compounds were identified, synthesized and tested. All of them showed activity against the H37Rv strain, and three molecules demonstrated activity against the MDR-TB strain. The docking analysis indicated that the discovered molecules could bind enoyl reductase, InhA, which is required in mycobacterial cell wall development. The models are freely available online (http://ochem.eu/article/103868) and can be used to predict potential anti-TB activity of new chemicals.

Genome-scale metabolic models as platforms for identification of novel genes as antimicrobial drug targets

The growing number of multidrug-resistant pathogenic bacteria is becoming a world leading challenge for the scientific community and for public health. However, advances in high-throughput technologies and whole-genome sequencing of bacterial pathogens make the construction of bacterial genome-scale metabolic models (GEMs) increasingly realistic. The use of GEMs as an alternative platforms will expedite identification of novel unconditionally essential genes and enzymes of target organisms with existing and forthcoming GEMs. This approach will follow the existing protocol for construction of high-quality GEMs, which could ultimately reduce the time, cost and labor-intensive processes involved in identification of novel antimicrobial drug targets in drug discovery pipelines. We discuss the current impact of existing GEMs of selected multidrug-resistant pathogenic bacteria for identification of novel antimicrobial drug targets and the challenges of closing the gap between genome-scale metabolic modeling and conventional experimental trial-and-error approaches in drug discovery pipelines.

An overview on crystal structures of InhA protein: Apo-form, in complex with its natural ligands and inhibitors

The enoyl-ACP reductase InhA from the mycobacterial fatty acid biosynthesis pathway has become a target of interest for the development of new anti-tubercular drugs. This protein has been identified as essential for the survival of Mycobacterium tuberculosis, the causative agent of tuberculosis, and as the main target of two pro-drugs: isoniazid, the frontline anti-tubercular drug, and ethionamide, a second-line medicine. Since most cases of resistance to isoniazid and ethionamide result from mutations in the mycobacterial activating enzyme (KatG for isoniazid and EthA for ethionamide), research of direct InhA inhibitors, avoiding the activation step, has emerged as a promising strategy for combating tuberculosis. Thereby, InhA is drawing much attention and its three-dimensional structure has been particularly studied. A better understanding of key sites of interactions responsible for InhA inhibition arises thus as an essential tool for the rational design of new potent inhibitors. In this paper, we propose an overview of the 80 available crystal structures of wild-type and mutant InhA, in its apo form, in complex with its cofactor, with an analogue of its natural ligands (C16 fatty acid and/or NADH) or with inhibitors. We will first discuss structural and mechanistic aspects in order to highlight key features of the protein before delivering thorough inventory of structures of InhA in the presence of synthetic ligands to underline the key interactions implicated in high affinity inhibition.

Drug-Target Kinetics in Drug Discovery

The development of therapies for the treatment of neurological cancer faces a number of major challenges including the synthesis of small molecule agents that can penetrate the blood-brain barrier (BBB). Given the likelihood that in many cases drug exposure will be lower in the CNS than in systemic circulation, it follows that strategies should be employed that can sustain target engagement at low drug concentration. Time dependent target occupancy is a function of both the drug and target concentration as well as the thermodynamic and kinetic parameters that describe the binding reaction coordinate, and sustained target occupancy can be achieved through structural modifications that increase target (re)binding and/or that decrease the rate of drug dissociation. The discovery and deployment of compounds with optimized kinetic effects requires information on the structure-kinetic relationships that modulate the kinetics of binding, and the molecular factors that control the translation of drug-target kinetics to time-dependent drug activity in the disease state. This Review first introduces the potential benefits of drug-target kinetics, such as the ability to delineate both thermodynamic and kinetic selectivity, and then describes factors, such as target vulnerability, that impact the utility of kinetic selectivity. The Review concludes with a description of a mechanistic PK/PD model that integrates drug-target kinetics into predictions of drug activity.

Protein conformational flexibility modulates kinetics and thermodynamics of drug binding

Structure-based drug design has often been restricted by the rather static picture of protein-ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery.

Intracellular rebinding of transition-state analogues provides extended in vivo inhibition lifetimes on human purine nucleoside phosphorylase

Purine nucleoside phosphorylase (PNP) is part of the human purine salvage pathway. Its deficiency triggers apoptosis of activated T-cells, making it a target for T-cell proliferative disorders. Transition-state analogues of PNP bind with picomolar (pm) dissociation constants. Tight-binding PNP inhibitors show exceptionally long lifetimes on the target enzyme. We solve the mechanism of the target residence time by comparing functional off-rates in vitro and in vivo We report in vitro PNP-inhibitor dissociation rates (t_½) from 3 to 31 min for seven Immucillins with dissociation constants of 115 to 6 pm Treatment of human erythrocytes with DADMe-Immucillin-H (DADMe-ImmH, 22 pm) causes complete inhibition of PNP. Loss of [¹⁴C]DADMe-ImmH from erythrocytes during multiple washes is slow and biphasic, resulting from inhibitor release and rebinding to PNP catalytic sites. The slow phase gave a t_½ of 84 h. Loss of [¹⁴C]DADMe-ImmH from erythrocytes in the presence of excess unlabeled DADMe-ImmH increased to a t_½ of 1.6 h by preventing rebinding. Thus, in human erythrocytes, rebinding of DADMe-ImmH is 50-fold more likely than diffusional loss of the inhibitor from the erythrocyte. Humans treated with a single oral dose of DADMe-ImmH in phase 1 clinical trials exhibit regain of PNP activity with a t_½ of 59 days, corresponding to the erythropoiesis rate in humans. Thus, the PNP catalytic site recapture of DADMe-ImmH is highly favored in vivo We conclude that transition-state analogues with picomolar dissociation constants exhibit long lifetimes on their targets in vivo because the probability of the target enzyme recapturing inhibitor molecules is greater than diffusional loss to the extracellular space.