Anders Karlén

Hydrogen bonding descriptors in the prediction of human in vivo intestinal permeability

Hydrogen bonding has been identified as an important parameter for describing drug permeability. Recently, we derived models for predicting intestinal permeability using the hydrogen bonding descriptors polar surface area (PSA) and number of hydrogen bond donors (HBD), and a lipophilicity descriptor [J. Med. Chem. 41 (1998) 4939]. We have now explored other types of hydrogen bonding descriptors to see if these improve the models. Both an experimental hydrogen bonding descriptor, deltalogP, and calculated descriptors, based either on semiempirical calculations or on experimentally derived hydrogen bond strength values of small molecules, were used. Principal component analyses (PCA) were performed in order to characterize the different parameters, using both a drug data set and a data set of small non-drug-like molecules for which deltalogP-values had been published. For a set of diverse drug molecules, for which human intestinal permeability data was available, a PLS-analysis was performed to study the correlation of permeability to the different hydrogen bonding parameters. No correlation could be identified between deltalogP and human intestinal permeability in this data set. However, the combination of a hydrogen bond donor descriptor, a general hydrogen bonding descriptor and a lipophilicity descriptor enabled the prediction of human intestinal permeability, whereas hydrogen bond acceptor descriptors were found to be less important. The obtained models successfully predicted the intestinal permeability values of two external data sets.

Functionalized 3-amino-imidazo[1,2-a]pyridines: A novel class of drug-like Mycobacterium tuberculosis glutamine synthetase inhibitors

3-Amino-imidazo[1,2-a]pyridines have been identified as a novel class of Mycobacterium tuberculosis glutamine synthetase inhibitors. Moreover, these compounds represent the first drug-like inhibitors of this enzyme. A series of compounds exploring structural diversity in the pyridine and phenyl rings have been synthesized and biologically evaluated. Compound 4n was found to be the most potent inhibitor (IC(50)=0.38+/-0.02 microM). This compound was significantly more potent than the known inhibitors, l-methionine-SR-sulfoximine and phosphinothricin.

Virtual Screening Data Fusion Using Both Structure- and Ligand-Based Methods

Virtual screening is widely applied in drug discovery, and significant effort has been put into improving current methods. In this study, we have evaluated the performance of compound ranking in virtual screening using five different data fusion algorithms on a total of 16 data sets. The data were generated by docking, pharmacophore search, shape similarity, and electrostatic similarity, spanning both structure- and ligand-based methods. The algorithms used for data fusion were sum rank, rank vote, sum score, Pareto ranking, and parallel selection. None of the fusion methods require any prior knowledge or input other than the results from the single methods and, thus, are readily applicable. The results show that compound ranking using data fusion improves the performance and consistency of virtual screening compared to the single methods alone. The best performing data fusion algorithm was parallel selection, but both rank voting and Pareto ranking also have good performance.

Design, Synthesis and X-Ray Crystallographic Studies of Alpha-Aryl Substituted Fosmidomycin Analogues as Inhibitors of Mycobacterium Tuberculosis 1-Deoxy-D-Xylulose-5-Phosphate Reductoisomerase

The natural antibiotic fosmidomycin acts via inhibition of 1-deoxy-d-xylulose 5-phosphate reductoisomerase (DXR), an essential enzyme in the non-mevalonate pathway of isoprenoid biosynthesis. Fosmidomycin is active on Mycobacterium tuberculosis DXR (MtDXR), but it lacks antibacterial activity probably because of poor uptake. α-Aryl substituted fosmidomycin analogues have more favorable physicochemical properties and are also more active in inhibiting malaria parasite growth. We have solved crystal structures of MtDXR in complex with 3,4-dichlorophenyl substituted fosmidomycin analogues; these show important differences compared to our previously described forsmidomycin-DXR complex. Our best inhibitor has an IC(50) = 0.15 μM on MtDXR but still lacked activity in a mycobacterial growth assay (MIC > 32 μg/mL). The combined results, however, provide insights into how DXR accommodates the new inhibitors and serve as an excellent starting point for the design of other novel and more potent inhibitors, particularly against pathogens where uptake is less of a problem, such as the malaria parasite.

Is it possible to increase hit rates in structure-based virtual screening by pharmacophore filtering? An investigation of the advantages and pitfalls of post-filtering

We have investigated the influence of post-filtering virtual screening results, with pharmacophoric features generated from an X-ray structure, on enrichment rates. This was performed using three docking softwares, zdock+, Surflex and FRED, as virtual screening tools and pharmacophores generated in UNITY from co-crystallized complexes. Sets of known actives along with 9997 pharmaceutically relevant decoy compounds were docked against six chemically diverse protein targets namely CDK2, COX2, ERalpha, fXa, MMP3, and NA. To try to overcome the inherent limitations of the well-known docking problem, we generated multiple poses for each compound. The compounds were first ranked according to their scores alone and enrichment rates were calculated using only the top scoring pose of each compound. Subsequently, all poses for each compound were passed through the different pharmacophores generated from co-crystallized complexes and the enrichment factors were re-calculated based on the top-scoring passing pose of each compound. Post-filtering with a pharmacophore generated from only one X-ray complex was shown to increase enrichment rates in all investigated targets compared to docking alone. This indicates that this is a general method, which works for diverse targets and different docking softwares.

Trisubstituted Imidazoles as Mycobacterium tuberculosis Glutamine Synthetase Inhibitors

Mycobacterium tuberculosis glutamine synthetase (MtGS) is a promising target for antituberculosis drug discovery. In a recent high-throughput screening study we identified several classes of MtGS inhibitors targeting the ATP-binding site. We now explore one of these classes, the 2-tert-butyl-4,5-diarylimidazoles, and present the design, synthesis, and X-ray crystallographic studies leading to the identification of MtGS inhibitors with submicromolar IC50 values and promising antituberculosis MIC values.

Discriminant and quantitative PLS analysis of competitive CYP2C9 inhibitors versus non-inhibitors using alignment independent GRIND descriptors

This study describes the use of alignment-independent descriptors for obtaining qualitative and quantitative predictions of the competitive inhibition of CYP2C9 on a serie of highly structurally diverse compounds. This was accomplished by calculating alignment independent descriptors in ALMOND. These GRid INdependent Descriptors (GRIND) represent the most important GRID-interactions as a function of the distance instead of the actual position of each grid-point. The experimental data was determined under uniform conditions. The inhibitor data set consists of 35 structurally diverse competitive stereospecific inhibitors of the cytochrome P450 2C9 and the non -inhibitor data set of 46 compounds. In a PLS discriminant analysis 21 inhibitors and 21 non-inhibitors (1 and 0 as activities) were analyzed using the ALMOND program obtaining a model with an r2 of 0.74 and a cross-validation value (q2) of 0.64. The model was externally validated with 39 compounds (14 inhibitors/25 non-inhibitors). 74% of the compounds were correctly predicted and an additional 13% was assigned to a borderline cluster. Thereafter, a model for quantitative predictions was generated by a PLS analysis of the GRIND descriptors using the experimental Ki-value for 21 of the competitive inhibitors (r2=0.77, q2=0.60). The model was externally validated using 12 compounds and predicted 11 out of 12 of the Ki-values within 0.5 log units. The discriminant model will be useful in screening for CYP2C9 inhibitors from large compound collections. The 3D-QSAR model will be used during lead optimization to avoid chemistry that result in inhibition of CYP2C9.

Characterization of a set of HIV-1 protease inhibitors using binding kinetics data from a biosensor-based screen

The interaction between 290 structurally diverse human immunodeficiency virus type 1 (HIV-1) protease inhibitors and the immobilized enzyme was analyzed with an optical biosensor. Although only a single concentration of inhibitor was used, information about the kinetics of the interaction could be obtained by extracting binding signals at discrete time points. The statistical correlation between the biosensor binding data, inhibition of enzyme activity (K;), and viral replication (EC50) revealed that the association and dissociation rates for the interaction could be resolved and that they were characteristic for the compounds. The most potent inhibitors, with respect to K; and EC50 values, including the clinically used drugs, all exhibited fast association and slow dissociation rates. Selective or partially selective binders for HIV-1 protease could be distinguished from compounds that showed a general protein-binding tendency by using three reference target proteins. This biosensor-based direct binding assay revealed a capacity to efficiently provide high-resolution information on the interaction kinetics and specificity of the interaction of a set of compounds with several targets simultaneously.

Synthesis of Functionalized Cinnamaldehyde Derivatives by an Oxidative Heck Reaction and Their Use as Starting Materials for Preparation of Mycobacterium tuberculosis 1-Deoxy-d-xylulose-5-phosphate Reductoisomerase Inhibitors

Cinnamaldehyde derivatives were synthesized in good to excellent yields in one step by a mild and selective, base-free palladium(II)-catalyzed oxidative Heck reaction starting from acrolein and various arylboronic acids. Prepared α,β-unsaturated aldehydes were used for synthesis of novel α-aryl substituted fosmidomycin analogues, which were evaluated for their inhibition of Mycobacterium tuberculosis 1-deoxy-d-xylulose 5-phosphate reductoisomerase. IC50 values between 0.8 and 27.3 μM were measured. The best compound showed activity comparable to that of the most potent previously reported α-aryl substituted fosmidomycin-class inhibitor.

Disulfide cyclized tripeptide analogues of angiotensin IV as potent and selective inhibitors of insulin-regulated aminopeptidase (IRAP)

The insulin-regulated aminopeptidase (IRAP) localized in areas of the brain associated with memory and learning is emerging as a new promising therapeutic target for the treatment of memory dysfunctions. The angiotensin II metabolite angiotensin IV (Ang IV, Val(1)-Tyr(2)-Ile(3)-His(4)-Pro(5)-Phe(6)) binds with high affinity to IRAP and inhibits this aminopeptidase (K(i) = 62.4 nM). Furthermore, Ang IV has been demonstrated to enhance cognition in animal models and is believed to play an important role in cognitive processes. It is herein reported that displacement of the C-terminal tripeptide His(4)-Pro(5)-Phe(6) with a phenylacetic acid functionality combined with a constrained macrocyclic system in the N-terminal affords potent IRAP inhibitors that are less peptidic in character than the hexapeptide Ang IV. Configurational analysis of three pairs of diastereomeric Ang IV analogues was performed using a combination of solution NMR spectroscopic methods, Monte Carlo conformational searches, and NAMFIS calculations. The compounds encompassing l-amino acids only (4, 8, and 12) showed significantly higher bioactivity compared to their lld-epimers (5, 9, and 13). The best inhibitors in the series, compounds 8 and 12, incorporating a 13- and 14-membered disulfide ring system, respectively, and both with a β(3)-homotyrosine residue (β(3)hTyr) replacing Tyr(2), exhibit K(i) values of 3.3 and 5.2 nM, respectively.