Araz Jakalian

Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation

We present the first global parameterization and validation of a novel charge model, called AM1‐BCC, which quickly and efficiently generates high‐quality atomic charges for computer simulations of organic molecules in polar media. The goal of the charge model is to produce atomic charges that emulate the HF/6‐31G* electrostatic potential (ESP) of a molecule. Underlying electronic structure features, including formal charge and electron delocalization, are first captured by AM1 population charges; simple additive bond charge corrections (BCCs) are then applied to these AM1 atomic charges to produce the AM1‐BCC charges. The parameterization of BCCs was carried out by fitting to the HF/6‐31G* ESP of a training set of >2700 molecules. Most organic functional groups and their combinations were sampled, as well as an extensive variety of cyclic and fused bicyclic heteroaryl systems. The resulting BCC parameters allow the AM1‐BCC charging scheme to handle virtually all types of organic compounds listed in The Merck Index and the NCI Database. Validation of the model was done through comparisons of hydrogen‐bonded dimer energies and relative free energies of solvation using AM1‐BCC charges in conjunction with the 1994 Cornell et al. forcefield for AMBER.13 Homo‐ and hetero‐dimer hydrogen‐bond energies of a diverse set of organic molecules were reproduced to within 0.95 kcal/mol RMS deviation from the ab initio values, and for DNA dimers the energies were within 0.9 kcal/mol RMS deviation from ab initio values. The calculated relative free energies of solvation for a diverse set of monofunctional isosteres were reproduced to within 0.69 kcal/mol of experiment. In all these validation tests, AMBER with the AM1‐BCC charge model maintained a correlation coefficient above 0.96. Thus, the parameters presented here for use with the AM1‐BCC method present a fast, accurate, and robust alternative to HF/6‐31G* ESP‐fit charges for general use with the AMBER force field in computer simulations involving organic small molecules.

Solvated interaction energy (SIE) for scoring protein-ligand binding affinities. 1. Exploring the parameter space.

We present a binding free energy function that consists of force field terms supplemented by solvation terms. We used this function to calibrate the solvation model along with the binding interaction terms in a self-consistent manner. The motivation for this approach was that the solute dielectric-constant dependence of calculated hydration gas-to-water transfer free energies is markedly different from that of binding free energies (J. Comput. Chem. 2003, 24, 954). Hence, we sought to calibrate directly the solvation terms in the context of a binding calculation. The five parameters of the model were systematically scanned to best reproduce the absolute binding free energies for a set of 99 protein-ligand complexes. We obtained a mean unsigned error of 1.29 kcal/mol for the predicted absolute binding affinity in a parameter space that was fairly shallow near the optimum. The lowest errors were obtained with solute dielectric values of Din = 20 or higher and scaling of the intermolecular van der Waals interaction energy by factors ranging from 0.03 to 0.15. The high apparent Din and strong van der Waals scaling may reflect the anticorrelation of the change in solvated potential energy and configurational entropy, that is, enthalpy-entropy compensation in ligand binding (Biophys. J. 2004, 87, 3035-3049). Five variations of preparing the protein-ligand data set were explored in order to examine the effect of energy refinement and the presence of bound water on the calculated results. We find that retaining water in the final protein structure used for calculating the binding free energy is not necessary to obtain good results; that is the continuum solvation model is sufficient. Virtual screening enrichment studies on estrogen receptor and thymidine kinase showed a good ability of the binding free energy function to recover true hits in a collection of decoys.

Revealing Atropisomer Axial Chirality in Drug Discovery

An often overlooked source of chirality is atropisomerism, which results from slow rotation along a bond axis due to steric hindrance and/or electronic factors. If undetected or not managed properly, this time‐dependent chirality has the potential to lead to serious consequences, because atropisomers can be present as distinct enantiomers or diastereoisomers with their attendant different properties. Herein we introduce a strategy to reveal and classify compounds that have atropisomeric chirality. Energy barriers to axial rotation were calculated using quantum mechanics, from which predicted high barriers could be experimentally validated. A calculated rotational energy barrier of 20 kcal mol−1 was established as a suitable threshold to distinguish between atropisomers and non‐atropisomers with a prediction accuracy of 86 %. This methodology was applied to subsets of drug databases in the course of which atropisomeric drugs were identified. In addition, some drugs were exposed that were not yet known to have this chiral attribute. The most valuable utility of this tool will be to predict atropisomerism along the drug discovery pathway. When used in concert with our compound classification scheme, decisions can be made during early discovery stages such as “hit‐to‐lead” and “lead optimization,” to foresee and validate the presence of atropisomers and to exercise options of removing, further stabilizing, or rendering the chiral axis of interest more freely rotatable via SAR design, thereby decreasing this potential liability within a compound series. The strategy can also improve drug development plans, such as determining whether a drug or series should be developed as a racemic mixture or as an isolated single compound. Moreover, the work described herein can be extended to other chemical fields that require the assessment of potential chiral axes.

Investigation on the role of the tetrazole in the binding of thiotetrazolylacetanilides with HIV-1 wild type and K103N/Y181C double mutant reverse transcriptases.

The role of the tetrazole moiety in the binding of aryl thiotetrazolylacetanilides with HIV-1 wild type and K103N/Y181C double mutant reverse transcriptases was explored. Different acyclic, cyclic and heterocyclic replacements were investigated in order to evaluate the conformational and electronic contribution of the tetrazole ring to the binding of the inhibitors in the NNRTI pocket. The replacement of the tetrazole by a pyrazolyl group led to reversal of selectivity, providing inhibitors with excellent potency against the double mutant reverse transcriptase.

Preclinical Profile of BI 224436, a Novel HIV-1 Non-Catalytic Site Integrase Inhibitor

ABSTRACT BI 224436 is an HIV-1 integrase inhibitor with effective antiviral activity that acts through a mechanism that is distinct from that of integrase strand transfer inhibitors (INSTIs). This 3-quinolineacetic acid derivative series was identified using an enzymatic integrase long terminal repeat (LTR) DNA 3′-processing assay. A combination of medicinal chemistry, parallel synthesis, and structure-guided drug design led to the identification of BI 224436 as a candidate for preclinical profiling. It has antiviral 50% effective concentrations (EC50s) of <15 nM against different HIV-1 laboratory strains and cellular cytotoxicity of >90 μM. BI 224436 also has a low, ∼2.1-fold decrease in antiviral potency in the presence of 50% human serum and, by virtue of a steep dose-response curve slope, exhibits serum-shifted EC95 values ranging between 22 and 75 nM. Passage of virus in the presence of inhibitor selected for either A128T, A128N, or L102F primary resistance substitutions, all mapping to a conserved allosteric pocket on the catalytic core of integrase. BI 224436 also retains full antiviral activity against recombinant viruses encoding INSTI resistance substitutions N155S, Q148H, and E92Q. In drug combination studies performed in cellular antiviral assays, BI 224436 displays an additive effect in combination with most approved antiretrovirals, including INSTIs. BI 224436 has drug-like in vitro absorption, distribution, metabolism, and excretion (ADME) properties, including Caco-2 cell permeability, solubility, and low cytochrome P450 inhibition. It exhibited excellent pharmacokinetic profiles in rat (clearance as a percentage of hepatic flow [CL], 0.7%; bioavailability [F], 54%), monkey (CL, 23%; F, 82%), and dog (CL, 8%; F, 81%). Based on the excellent biological and pharmacokinetic profile, BI 224436 was advanced into phase 1 clinical trials.

Development of Specific “Drug‐Like Property” Rules for Carboxylate‐Containing Oral Drug Candidates

The carboxylate moiety is an important pharmacophore in the medicinal chemists arsenal and is sometimes an irreplaceable functionality in drug–target interactions. Thus, practical guidance on its use in the most optimized manner would be a welcome addition to rational drug design. Key physicochemical and ADMET‐PK properties from a dataset of drugs containing a carboxylate (COOH) moiety were assembled and compared with those of a broader, general drug dataset. Our main objective was to identify features specific to COOH‐containing oral drugs that could be converted into simple rules delineating the boundaries within which prospective COOH‐containing chemical series and COOH‐containing drug candidates would be reasonably expected to possess properties suitable for oral administration. These specific “drug‐like” property rules include molecular weight, the number of rotatable bonds, the number of hydrogen bond donors and acceptors, predictions of lipophilic character (calculated log P and log D values), topological polar surface area (TPSA), and the pKa value of the carboxylate moiety. Similar to the various sets of criteria that have emerged over the past decade and which have significantly reshaped the way medicinal chemists think about preferred drug chemical space, we propose these specific COOH “drug‐like” property rules as a guide for the design of superior COOH‐containing drug candidates and as a tool to better manage the liabilities generally associated with the presence of a COOH moiety.

Enantiomeric Atropisomers Inhibit HCV Polymerase and/or HIV Matrix: Characterizing Hindered Bond Rotations and Target Selectivity.

An anthranilic acid series of allosteric thumb pocket 2 HCV NS5B polymerase inhibitors exhibited hindered rotation along a covalent bond axis, and the existence of atropisomer chirality was confirmed by NMR, HPLC analysis on chiral supports, and computational studies. A thorough understanding of the concerted rotational properties and the influence exerted by substituents involved in this steric phenomenon was attained through biophysical studies on a series of truncated analogues. The racemization half-life of a compound within this series was determined to be 69 min, which was consistent with a class 2 atropisomer (intermediate conformational exchange). It was further found by X-ray crystallography that one enantiomer of a compound bound to the intended HCV NS5B polymerase target whereas the mirror image atropisomer was able to bind to an unrelated HIV matrix target. Analogues were then identified that selectively inhibited the former. These studies highlight that atropisomer chirality can lead to distinct entities with specific properties, and the phenomenon of atropisomerism in drug discovery should be evaluated and appropriately managed.

Ligand bioactive conformation plays a critical role in the design of drugs that target the hepatitis C virus NS3 protease.

A ligand-focused strategy employed NMR, X-ray, modeling, and medicinal chemistry to expose the critical role that bioactive conformation played in the design of a variety of drugs that target the HCV protease. The bioactive conformation (bound states) were determined for key inhibitors identified along our drug discovery pathway from the hit to clinical compounds. All adopt similar bioactive conformations for the common core derived from the hit peptide DDIVPC. A carefully designed SAR analysis, based on the advanced inhibitor 1 in which the P1 to P3 side chains and the N-terminal Boc were sequentially truncated, revealed a correlation between affinity and the relative predominance of the bioactive conformation in the free state. Interestingly, synergistic conformation effects on potency were also noted. Comparisons with clinical and recently marketed drugs from the pharmaceutical industry showed that all have the same core and similar bioactive conformations. This suggested that the variety of appendages discovered for these compounds also properly satisfy the bioactive conformation requirements and allowed for a large variety of HCV protease drug candidates to be designed.

N-Acetamideindolecarboxylic acid allosteric ‘finger-loop’ inhibitors of the hepatitis C virus NS5B polymerase: discovery and initial optimization studies

SAR studies at the N(1)-position of allosteric indole-based HCV NS5B inhibitors has led to the discovery of acetamide derivatives with good cellular potency in subgenomic replicons (EC(50) <200 nM). This class of inhibitors displayed improved physicochemical properties and favorable ADME-PK profiles over previously described analogs in this class.

Conformation-Based Restrictions and Scaffold Replacements in the Design of Hepatitis C Virus Polymerase Inhibitors: Discovery of Deleobuvir (BI 207127)

Conformational restrictions of flexible torsion angles were used to guide the identification of new chemotypes of HCV NS5B inhibitors. Sites for rigidification were based on an acquired conformational understanding of compound binding requirements and the roles of substituents in the free and bound states. Chemical bioisosteres of amide bonds were explored to improve cell-based potency. Examples are shown, including the design concept that led to the discovery of the phase III clinical candidate deleobuvir (BI 207127). The structure-based strategies employed have general utility in drug design.