Fluoroquinolone Derivatives in the Treatment of Mycobacterium tuberculosis Infection.

Tuberculosis (TB) is currently one of the leading causes of death due to infective agents, and the growing rate of multidrug-resistant tuberculosis (MDR TB) cases poses an emergent public health threat. Fluoroquinolones are commonly used in the treatment of both MDR TB and drug-sensitive tuberculosis patients who are intolerant to first-line antitubercular agents. Unfortunately, these drugs have mild side effects, relevant to the prolonged treatment regimens and diminished bioavailability due to binding of metal ions. Moreover, the resistance to fluoroquinolones is also on the rise, a characteristic of extensively drug-resistant TB (XDR TB). Here, we developed esters as prodrugs of the fluoroquinolones levofloxacin and ciprofloxacin, with long-chain fatty alcohols. Both the alcohols and the quinolone have previously shown antimycobacterial activity and the aim was to develop esters with improved lipophilicity and capable of delivering the free acid inside mycobacterial cells. The carboxylic acid group of fluoroquinolones is essential to the mode of action but is also responsible for many of its side effects and metal-chelating properties. The synthesis, stability in biological media, and antibacterial activity were evaluated, the latter not only against Mycobacterium tuberculosis but also against other clinically relevant bacterial species, since the parent compounds display a broad spectrum of activity. The biological results show a reduction in the antitubercular activity of the synthesized derivatives, probably due to deficient activation of the ester prodrug. Despite this, it was found that the derivatives exhibit bioactivity against other fluoroquinolone-resistant bacteria, indicating a different mode of action and suggesting that it may be worthwhile to research further modifications to the carboxylic acid group. This might lead to new compounds that are efficient against resistant strains. This idea that the compounds may act by a different mechanism of action was further supported by a brief computer investigation that demonstrated the potential lack of selectivity of the esters to the fluoroquinolone target.

Benzoic Acid Derivatives as Prodrugs for the Treatment of Tuberculosis.

One interesting approach to fight tuberculosis is the use of prodrugs that often have shown improved biological activities over drugs with poor absorption or difficulty to cross membranes. Previous studies demonstrate that weak acids such as benzoic acid, present antimycobacterial activity. Moreover, esters of those acids revealed to be a viable alternative since they may diffuse more easily through the cell membranes. Previously we showed that mycobacteria can easily activate benzoic acid esters by conversion to the corresponding acid. Since Zhang postulated that the activity of the acids can be dependent on their pKa, we set up to synthesize a library of benzoates with different electron withdrawing groups (4-chloro, 2,6-dichloro, 3,5-dichloro, 4-nitro, and 3,5 dinitro), to modulate pKa of the liberated acid and different alkoxy substituents (propyl, hexyl, and phenyl) to modulate their lipophilicity, and tested the activity of the esters and the corresponding free acids against mycobacteria. We also studied the activation of the esters by mycobacterial enzymes and the stability of the compounds in buffer and plasma. We concluded that all the benzoates in our study can be activated by mycobacterial enzymes and that the phenyl and hexyl esters presented higher activity than the corresponding free acids, with the nitrobenzoates, and especially the dinitrobenzoates, showing very interesting antitubercular activity that deserve further exploration. Our results did not show a correlation between the activity and the pKa of the acids.

Activation of benzoate model prodrugs by mycobacteria. Comparison with mammalian plasma and liver hydrolysis.

Due to difficulties in drug penetration in M. tuberculosis, a prodrug approach based on mycobacterial activation appears as a promising strategy to increase the delivery of antitubercular drugs to the target microorganisms. Esters have been successful used by us and others to deliver drugs to mycobacteria, however because very little is known about the metabolic hydrolysis of esters by mycobacteria in connection with prodrug activation, we decided to study the process further. For that we selected a series of 13 benzoates with different chain lengths and ramifications in the alkoxy side as model prodrugs and examined their hydrolysis by a mycobacterial homogenate, comparing the results with those obtained parallelly in human plasma and in total rat liver homogenate. In all biological media, the benzoates with a linear alkyl group showed a parabolic dependence between log(k) and logP (or the number of carbons of the linear alkyl chain) that reached a maximal value for the n-butyl chain. Considering linear correlations for the total number of compounds between log(k) and chosen descriptors, for mycobacterial esterases, pKa of the leaving alcohol (pKaLG) seem to be the most important descriptor. Plasma esterases seem to be quite sensitive to the Taft polarity parameter σ* and also to pKaLG and less sensitive to steric effects. Liver esterases seem to be more sensitive to the Taft steric descriptor ESc. Lipophilicity correlates weakly with log(k) in all the 3 media, however, is more important when one looks for mycobacterial activation selectivity in relation to plasma metabolism or in relation to liver homogenate metabolism. The importance of lipophilicity increases further when biparametric expressions are considered. We showed that it is easy to activate a wide variety of benzoate esters using a mycobacterial homogenate. The data also suggest that with careful design is possible to obtain tuberculostatic prodrug esters sensitive to mycobacterial hydrolases while reasonably resistant to plasma and liver hydrolysis. One important observation is that mycobacterial hydrolysis is less affected by bulky substituents than liver homogenate or plasma hydrolysis. tert-Butyl is probably the substituent in the alkoxy side that seems more adequate to resist simultaneously plasma and liver metabolism, while allowing activation by mycobacterial esterases. Hexyl is also a good option for the medicinal chemist if a linear alkoxy chain is needed.

FAME 3: Predicting the Sites of Metabolism in Synthetic Compounds and Natural Products for Phase 1 and Phase 2 Metabolic Enzymes.

In this work we present the third generation of FAst MEtabolizer (FAME 3), a collection of extra trees classifiers for the prediction of sites of metabolism (SoMs) in small molecules such as drugs, druglike compounds, natural products, agrochemicals, and cosmetics. FAME 3 was derived from the MetaQSAR database ( Pedretti et al. J. Med. Chem. 2018 , 61 , 1019 ), a recently published data resource on xenobiotic metabolism that contains more than 2100 substrates annotated with more than 6300 experimentally confirmed SoMs related to redox reactions, hydrolysis and other nonredox reactions, and conjugation reactions. In tests with holdout data, FAME 3 models reached competitive performance, with Matthews correlation coefficients (MCCs) ranging from 0.50 for a global model covering phase 1 and phase 2 metabolism, to 0.75 for a focused model for phase 2 metabolism. A model focused on cytochrome P450 metabolism yielded an MCC of 0.57. Results from case studies with several synthetic compounds, natural products, and natural product derivatives demonstrate the agreement between model predictions and literature data even for molecules with structural patterns clearly distinct from those present in the training data. The applicability domains of the individual models were estimated by a new, atom-based distance measure (FAMEscore) that is based on a nearest-neighbor search in the space of atom environments. FAME 3 is available via a public web service at https://nerdd.zbh.uni-hamburg.de/ and as a self-contained Java software package, free for academic and noncommercial research.

Prediction of UGT-mediated Metabolism Using the Manually Curated MetaQSAR Database.

Even though glucuronidations are the most frequent metabolic reactions of conjugation, both in quantitative and qualitative terms, they have rather seldom been investigated using computational approaches. To fill this gap, we have used the manually collected MetaQSAR metabolic reaction database to generate two models for the prediction of UGT-mediated metabolism, both based on molecular descriptors and implementing the Random Forest algorithm. The first model predicts the occurrence of the reaction and was internally validated with a Matthew correlation coefficient (MCC) of 0.76 and an area under the ROC curve (AUC) of 0.94, and further externally validated using a test set composed of 120 additional xenobiotics (MCC of 0.70 and AUC of 0.90). The second model distinguishes between O- and N-glucuronidations and was optimized by the random undersampling procedure to improve the predictive accuracy during the internal validation, with the recall measure of the minority class increasing from 0.55 to 0.78.

Prediction of the Formation of Reactive Metabolites by A Novel Classifier Approach Based on Enrichment Factor Optimization (EFO) as Implemented in the VEGA Program.

The study is aimed at developing linear classifiers to predict the capacity of a given substrate to yield reactive metabolites. While most of the hitherto reported predictive models are based on the occurrence of known structural alerts (e.g., the presence of toxophoric groups), the present study is focused on the generation of predictive models involving linear combinations of physicochemical and stereo-electronic descriptors. The development of these models is carried out by using a novel classification approach based on enrichment factor optimization (EFO) as implemented in the VEGA suite of programs. The study took advantage of metabolic data as collected by manually curated analysis of the primary literature and published in the years 2004⁻2009. The learning set included 977 substrates among which 138 compounds yielded reactive first-generation metabolites, plus 212 substrates generating reactive metabolites in all generations (i.e., metabolic steps). The results emphasized the possibility of developing satisfactory predictive models especially when focusing on the first-generation reactive metabolites. The extensive comparison of the classifier approach presented here using a set of well-known algorithms implemented in Weka 3.8 revealed that the proposed EFO method compares with the best available approaches and offers two relevant benefits since it involves a limited number of descriptors and provides a score-based probability thus allowing a critical evaluation of the obtained results. The last analyses on non-cheminformatics UCI datasets emphasize the general applicability of the EFO approach, which conveniently performs using both balanced and unbalanced datasets.

Emerging Chemodiversity and Stereoisomerism in Acyclic Mono- and Dichloroalkanes.

The objective of this pen-and-paper study is to monitor the growth of structural complexity as a function of the number of heavy atoms. This is done here using simple chloroalkanes, the indices of structural complexity being the number of regioisomers and mainly stereoisomeric features. Following a previous work dealing with acyclic alkanes, the present study focuses on all acyclic monochloroalkanes from four heavy atoms (four regioisomers) to 10 heavy atoms (2 to 211 regioisomers, respectively). Similarly, all acyclic dichloroalkanes were drawn, ranging from five heavy atoms to nine heavy atoms (4 to 130 regioisomers, respectively). The first part of the work describes the nature and number of stereoisomeric features seen to emerge in molecular structures featuring one or more stereogenic centers. These molecular features include enantiomerism, diastereoisomerism, pseudoasymmetry, and meso-isomerism. The first appearance of chirality among the monochloro and dichloro homologs occurred in the families with five heavy atoms. Diastereoisomerism, meso-isomerism, and pseudoasymmetry appear in higher families of homologs. The second part of the article is an attempt to rationalize the appearance of regioisomers, stereogenic centers, and stereoisomers as their number grows along homologous series, from one family of regioisomers to the higher ones. Far from being random as anticipated, such increases appear regular and ordered. Indeed, the data available demonstrate exponential increases and suggest the involvement of geometric order in the growth of chemical complexity.

Approaching Pharmacological Space: Events and Components.

With a view to introducing the concept of pharmacological space and its potential applications in investigating and predicting the toxic mechanisms of xenobiotics, this opening chapter describes the logical relations between conformational behavior, physicochemical properties and binding spaces, which are seen as the three key elements composing the pharmacological space. While the concept of conformational space is routinely used to encode molecular flexibility, the concepts of property spaces and, particularly, of binding spaces are more innovative. Indeed, their descriptors can find fruitful applications (a) in describing the dynamic adaptability a given ligand experiences when inserted into a specific environment, and (b) in parameterizing the flexibility a ligand retains when bound to a biological target. Overall, these descriptors can conveniently account for the often disregarded entropic factors and as such they prove successful when inserted in ligand- or structure-based predictive models. Notably, and although binding space parameters can clearly be derived from MD simulations, the chapter will illustrate how docking calculations, despite their static nature, are able to evaluate ligand's flexibility by analyzing several poses for each ligand. Such an approach, which represents the founding core of the binding space concept, can find various applications in which the related descriptors show an impressive enhancing effect on the statistical performances of the resulting predictive models.

MetaQSAR: An Integrated Database Engine to Manage and Analyze Metabolic Data.

The study describes the MetaQSAR tool, a new database engine specifically tailored to collect and analyze metabolic data. This is a plug-in embedded in the VEGA suite of programs (freely downloadable at www.vegazz.net ) and takes advantage from all cheminformatics features implemented in the software with additional tools aimed to perform statistical analyses, similarity searches, and physicochemical profiling of the stored molecules. MetaQSAR also implements a novel metabolism classification, which groups the metabolic reactions in 101 classes and can find numerous applications in metabolic analyses. The potentials of MetaQSAR are here assessed by using it to store and analyze an extended database focused on metabolism of xenobiotics, which was collected by manually curated meta-analysis of the recent literature. The database includes 1890 substrates taken from about 1500 original papers in the years 2004-2015. The database was utilized in both physicochemical analyses and similarity searches, so evidencing the notable potentialities of MetaQSAR, which can find particularly fruitful applications in developing targeted predictive approaches.

Emergent chemodiversity: The case of stereoisomerism in acyclic alkanes.

The objective of this pen-and-paper study is to witness the emergence of stereoisomeric properties when comparing lower to higher families of homologs. Specifically, the study compares all acyclic hexanes (five constitutional isomers, a.k.a. regioisomers), all nine heptanes, all 18 octanes, all 35 nonanes, and all 75 decanes. The first part of the work examines the nature and number of stereoisomeric properties seen to emerge in chemical structures featuring one chiral center (i.e., enantiomerism) or two such centers, in which case more complex stereoisomeric features emerge (enantiomerism, diastereoisomerism, pseudoasymmetry, and meso-isomers). The first emergence of chirality (i.e., enantiomerism) occurred in some heptanes. Diastereoisomerism and meso-isomers appear with some octanes, while a pseudoasymmetric center exists in a decane regioisomer. The second part of the work is an attempt to rationalize the numbers of regioisomers, chiral centers, and stereoisomers as these numbers grow from one family of regioisomers to the higher ones. Far from being random, such increases prove regular and ordered.

Binding Space Concept: A New Approach To Enhance the Reliability of Docking Scores and Its Application to Predicting Butyrylcholinesterase Hydrolytic Activity.

Docking simulations are very popular approaches able to assess the capacity of a given ligand to interact with a target. Docking simulations are usually focused on a single best complex even though many studies showed that ligands retain a significant mobility within a binding pocket by assuming different binding modes all of which may contribute to the monitored ligand affinity. The present study describes an innovative concept, the binding space, which allows an exploration of the ligand mobility within the binding pocket by simultaneously considering several ligand poses as generated by docking simulations. The multiple poses and the relative docking scores can then be analyzed by taking advantage of the same concepts already used in the property space analysis; hence the binding space can be parametrized by (a) mean scores, (b) score ranges, and (c) score sensitivity values. The first parameter represents a very simple procedure to account for the contribution of the often neglected alternative binding modes, while the last two descriptors encode the degree of mobility which a given ligand retains within the binding cavity (score range) as well as the ease with which a ligand explores such a mobility (score sensitivity). Here, the binding space concept is applied to the prediction of the hydrolytic activity of BChE by synergistically considering multiple poses and multiple protein structures. The obtained results shed light on the remarkable potential of the binding space concept, whose parameters allow a significant increase of the predictive power of the docking results as revealed by extended correlative analyses. Mean scores are the parameters affording the largest statistical improvement, and all the here proposed docking-based descriptors show enhancing effects in developing predictive models. Finally, the study describes a new score function (Contacts score) simply based on the number of surrounding residues which appears to be particularly productive in the framework of the binding space.

Mechanisms and pharmaceutical consequences of processes of stereoisomerisation - A didactic excursion.

This writing proposes an overview of chemical processes of isomerisation involving conformations and configurational isomers. Molecular structure is depicted as consisting in A) core features; and B) fluctuating properties. To avoid confusion and misuses in terminology, attention is paid to IUPAC recommendations and definitions, beginning with a focus on conformers versus conformations. Two sections consider rotations, the behaviour of exocyclic single bonds being kept separate from that of endocyclic single bonds. Processes of pyramidal inversions are discussed next, alone and combined with rotations. Given its sometimes unsuspected occurrence in some drug molecules, stereolability at tetracoordinated carbon centres receives particular attention. Stereolabile and stereostable nitrogen-containing groups come next, with a focus on the simultaneous contributions of inversion and rotation as occurring for example in the conformational reversal of saturated aza-heterocycles. The last type of isomerisation to be examined involves complex heterocyclic systems of medicinal relevance. The conclusion points to the mechanistic variety of stereoisomerisation processes, to a continuum in their energy barriers, to their relevance in all fields of stereochemistry, and to the need for junior medicinal chemists to be adequately trained in stereochemistry.

Esters of Pyrazinoic Acid Are Active against Pyrazinamide-Resistant Strains of Mycobacterium tuberculosis and Other Naturally Resistant Mycobacteria In Vitro and Ex Vivo within Macrophages.

Pyrazinamide (PZA) is active against major Mycobacterium tuberculosis species (M. tuberculosis, M. africanum, and M. microti) but not against M. bovis and M. avium. The latter two are mycobacterial species involved in human and cattle tuberculosis and in HIV coinfections, respectively. PZA is a first-line agent for the treatment of human tuberculosis and requires activation by a mycobacterial pyrazinamidase to form the active metabolite pyrazinoic acid (POA). As a result of this mechanism, resistance to PZA, as is often found in tuberculosis patients, is caused by point mutations in pyrazinamidase. In previous work, we have shown that POA esters and amides synthesized in our laboratory were stable in plasma (M. F. Simões, E. Valente, M. J. Gómez, E. Anes, and L. Constantino, Eur J Pharm Sci 37:257-263, 2009, http://dx.doi.org/10.1016/j.ejps.2009.02.012). Although the amides did not present significant activity, the esters were active against sensitive mycobacteria at concentrations 5- to 10-fold lower than those of PZA. Here, we report that these POA derivatives possess antibacterial efficacy in vitro and ex vivo against several species and strains of Mycobacterium with natural or acquired resistance to PZA, including M. bovis and M. avium. Our results indicate that the resistance probably was overcome by cleavage of the prodrugs into POA and a long-chain alcohol. Although it is not possible to rule out that the esters have intrinsic activity per se, we bring evidence here that long-chain fatty alcohols possess a significant antimycobacterial effect against PZA-resistant species and strains and are not mere inactive promoieties. These findings may lead to candidate dual drugs having enhanced activity against both PZA-susceptible and PZA-resistant isolates and being suitable for clinical development.

Predicting drug metabolism: experiment and/or computation?

Drug metabolism can produce metabolites with physicochemical and pharmacological properties that differ substantially from those of the parent drug, and consequently has important implications for both drug safety and efficacy. To reduce the risk of costly clinical-stage attrition due to the metabolic characteristics of drug candidates, there is a need for efficient and reliable ways to predict drug metabolism in vitro, in silico and in vivo. In this Perspective, we provide an overview of the state of the art of experimental and computational approaches for investigating drug metabolism. We highlight the scope and limitations of these methods, and indicate strategies to harvest the synergies that result from combining measurement and prediction of drug metabolism.

Types of stereoselectivity in drug metabolism: a heuristic approach.

Stereochemical factors are known to play a significant role in the metabolism of drugs and other xenobiotics. Following Prelog's lead, types of metabolic stereoselectivity can be categorized as (i) substrate stereoselectivity (the differential metabolism of two or more stereoisomeric substrates) and (ii) product stereoselectivity (the differential formation of two or more stereoisomeric metabolites from a single substrate). Combinations of the two categories exist as (iii) substrate-product stereoselectivities, meaning that product stereoselectivity itself is substrate stereoselective. Here, published examples of metabolic stereoselectivities are examined in the light of these concepts. In parallel, a graphical scheme is presented with a view to facilitate learning and help researchers to solve classification problems.

Passive lipoidal diffusion and carrier-mediated cell uptake are both important mechanisms of membrane permeation in drug disposition.

Recently, it has been proposed that drug permeation is essentially carrier-mediated only and that passive lipoidal diffusion is negligible. This opposes the prevailing hypothesis of drug permeation through biological membranes, which integrates the contribution of multiple permeation mechanisms, including both carrier-mediated and passive lipoidal diffusion, depending on the compound's properties, membrane properties, and solution properties. The prevailing hypothesis of drug permeation continues to be successful for application and prediction in drug development. Proponents of the carrier-mediated only concept argue against passive lipoidal diffusion. However, the arguments are not supported by broad pharmaceutics literature. The carrier-mediated only concept lacks substantial supporting evidence and successful applications in drug development.

Quantitative monitoring of tamoxifen in human plasma extended to 40 metabolites using liquid-chromatography high-resolution mass spectrometry: new investigation capabilities for clinical pharmacology.

Liquid-chromatography (LC) high-resolution (HR) mass spectrometry (MS) analysis can record HR full scans, a technique of detection that shows comparable selectivity and sensitivity to ion transitions (SRM) performed with triple-quadrupole (TQ)-MS but that allows de facto determination of "all" ions including drug metabolites. This could be of potential utility in in vivo drug metabolism and pharmacovigilance studies in order to have a more comprehensive insight in drug biotransformation profile differences in patients. This simultaneous quantitative and qualitative (Quan/Qual) approach has been tested with 20 patients chronically treated with tamoxifen (TAM). The absolute quantification of TAM and three metabolites in plasma was realized using HR- and TQ-MS and compared. The same LC-HR-MS analysis allowed the identification and relative quantification of 37 additional TAM metabolites. A number of new metabolites were detected in patients' plasma including metabolites identified as didemethyl-trihydroxy-TAM-glucoside and didemethyl-tetrahydroxy-TAM-glucoside conjugates corresponding to TAM with six and seven biotransformation steps, respectively. Multivariate analysis allowed relevant patterns of metabolites and ratios to be associated with TAM administration and CYP2D6 genotype. Two hydroxylated metabolites, α-OH-TAM and 4'-OH-TAM, were newly identified as putative CYP2D6 substrates. The relative quantification was precise (<20 %), and the semiquantitative estimation suggests that metabolite levels are non-negligible. Metabolites could play an important role in drug toxicity, but their impact on drug-related side effects has been partially neglected due to the tremendous effort needed with previous MS technologies. Using present HR-MS, this situation should evolve with the straightforward determination of drug metabolites, enlarging the possibilities in studying inter- and intra-patients drug metabolism variability and related effects.

Exploring the space of histidine containing dipeptides in search of novel efficient RCS sequestering agents.

The study reports a set of forty proteinogenic histidine-containing dipeptides as potential carbonyl quenchers. The peptides were chosen to cover as exhaustively as possible the accessible chemical space, and their quenching activities toward 4-hydroxy-2-nonenal (HNE) and pyridoxal were evaluated by HPLC analyses. The peptides were capped at the C-terminus as methyl esters or amides to favor their resistance to proteolysis and diastereoisomeric pairs were considered to reveal the influence of configuration on quenching. On average, the examined dipeptides are less active than the parent compound carnosine (ßAla + His) thus emphasizing the unfavorable effect of the shortening of the ßAla residue as confirmed by the control dipeptide Gly-His. Nevertheless, some peptides show promising activities toward HNE combined with a remarkable selectivity. The results emphasize the beneficial role of aromatic and positively charged residues, while negatively charged and H-bonding side chains show a detrimental effect on quenching. As a trend, ester derivatives are slightly more active than amides while heterochiral peptides are more active than their homochiral diastereoisomer. Overall, the results reveal that quenching activity strongly depends on conformational effects and vicinal residues (as evidenced by the reported QSAR analysis), offering insightful clues for the design of improved carbonyl quenchers and to rationalize the specific reactivity of histidine residues within proteins.

Simulations in evolution. III. Randomness as a generator of opportunities.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous reports presented a histogram model to simulate the evolution of populations of individuals classified into bins according to an unspecified, quantifiable phenotypic character, and whose number in each bin changed generation after generation under the influence of fitness, while the total population was maintained constant. The histogram model also allowed Shannon entropy (SE) to be monitored continuously as the information content of the total population decreased or increased. Here, a simple Perl (Practical Extraction and Reporting Language) application was developed to carry out these computations, with the critical feature of an added random factor in the percent of individuals whose offspring moved to a vicinal bin. The results of the simulations demonstrate that the random factor mimicking variation increased considerably the range of values covered by Shannon entropy, especially when the percentage of changed offspring was high. This increase in information content is interpreted as facilitated adaptability of the population.