Design and synthesis of a series of potent and orally bioavailable noncovalent thrombin inhibitors that utilize nonbasic groups in the P1 position

Design and Preclinical Characterization Program toward Asundexian (BAY 2433334), an Oral Factor XIa Inhibitor for the Prevention and Treatment of Thromboembolic Disorders

Activated coagulation factor XI (FXIa) is a highly attractive antithrombotic target as it contributes to the development and progression of thrombosis but is thought to play only a minor role in hemostasis so that its inhibition may allow for decoupling of antithrombotic efficacy and bleeding time prolongation. Herein, we report our major efforts to identify an orally bioavailable, reversible FXIa inhibitor. Using a protein structure-based de novo design approach, we identified a novel micromolar hit with attractive physicochemical properties. During lead modification, a critical problem was balancing potency and absorption by focusing on the most important interactions of the lead series with FXIa while simultaneously seeking to improve metabolic stability and the cytochrome P450 interaction profile. In clinical trials, the resulting compound from our extensive research program, asundexian (BAY 2433334), proved to possess the desired DMPK properties for once-daily oral dosing, and even more importantly, the initial pharmacological hypothesis was confirmed.

Experimental and computational evidence for a stabilising C–Cl(lone-pair)⋯π(chelate-ring) interaction

Stabilising C–Cl(lone-pair)⋯π(chelate ring) interactions are described.

Free Fatty Acids Modulate Thrombin Mediated Fibrin Generation Resulting in Less Stable Clots

Upon platelet activation, free fatty acids are released at the stage of thrombus formation, but their effects on fibrin formation are largely unexplored. Our objective was to characterize the kinetic effects of fatty acids on thrombin activity, as well as the structural and mechanical properties of the resultant fibrin clots. Thrombin activity on fibrinogen was followed by turbidimetry and detailed kinetic characterization was performed using a fluorogenic short peptide substrate. The viscoelastic properties of fibrin were measured with rotatory oscillation rheometer, whereas its structure was analyzed with scanning electron microscopy (SEM). In turbidimetric assays of fibrin generation, oleate and stearate at physiologically relevant concentrations (60–600 μM) produced a bell-shaped inhibitory dose response, increasing 10- to 30-fold the time to half-maximal clotting. Oleate inhibited thrombin activity on a short peptide substrate according to a mixed-type inhibitor pattern (a 9-fold increase of the Michaelis constant, K_m and a 20% decrease of the catalytic constant), whereas stearate resulted in only a minor (15%) drop in the catalytic constant without any change in the K_m. Morphometric analysis of SEM images showed a 73% increase in the median fiber diameter in the presence of stearate and a 20% decrease in the presence of oleate. Concerning the viscoelastic parameters of the clots, storage and loss moduli, maximal viscosity and critical shear stress decreased by 32–65% in the presence of oleate or stearate, but loss tangent did not change indicating decreased rigidity, higher deformability and decreased internal resistance to shear stress. Our study provides evidence that free fatty acids (at concentrations comparable to those reported in thrombi) reduce the mechanical stability of fibrin through modulation of thrombin activity and the pattern of fibrin assembly.

4-Fluoroprolines: Conformational Analysis and Effects on the Stability and Folding of Peptides and Proteins

Proline is unique among proteinogenic amino acids because a pyrrolidine ring links its amino group to its side chain. This heterocycle constrains the conformations of the main chain and thus templates particular secondary structures. Proline residues undergo post-translational modification at the 4-position to yield 4-hydroxyproline, which is especially prevalent in collagen. Interest in characterizing the effects of this modification led to the use of 4-fluoroprolines to enhance inductive properties relative to the hydroxyl group of 4-hydroxyproline and to eliminate contributions from hydrogen bonding. The strong inductive effect of the fluoro group has three main consequences: enforcing a particular pucker upon the pyrrolidine ring, biasing the conformation of the preceding peptide bond, and accelerating cis/trans prolyl peptide bond isomerization. These subtle, yet reliable modulations make 4-fluoroproline-incorporation a complement to traditional genetic approaches for exploring structure-function relationships in peptides and proteins, as well as for endowing peptides and proteins with conformational stability.

Discovery of AZD8165 – a clinical candidate from a novel series of neutral thrombin inhibitors

A novel series of neutral thrombin inhibitors has been developed using a selection process based on docking experiments and property calculations and predictions.

Discovery of novel P1 groups for coagulation factor VIIa inhibition using fragment-based screening

A multidisciplinary, fragment-based screening approach involving protein ensemble docking and biochemical and NMR assays is described. This approach led to the discovery of several structurally diverse, neutral surrogates for cationic factor VIIa P1 groups, which are generally associated with poor pharmacokinetic (PK) properties. Among the novel factor VIIa inhibitory fragments identified were aryl halides, lactams, and heterocycles. Crystallographic structures for several bound fragments were obtained, leading to the successful design of a potent factor VIIa inhibitor with a neutral lactam P1 and improved permeability.

3D-QSAR based on quantum-chemical molecular fields: toward an improved description of halogen interactions

Current 3D-QSAR methods such as CoMFA or CoMSIA make use of classical force-field approaches for calculating molecular fields. Thus, they can not adequately account for noncovalent interactions involving halogen atoms like halogen bonds or halogen-π interactions. These deficiencies in the underlying force fields result from the lack of treatment of the anisotropy of the electron density distribution of those atoms, known as the "σ-hole", although recent developments have begun to take specific interactions such as halogen bonding into account. We have now replaced classical force field derived molecular fields by local properties such as the local ionization energy, local electron affinity, or local polarizability, calculated using quantum-mechanical (QM) techniques that do not suffer from the above limitation for 3D-QSAR. We first investigate the characteristics of QM-based local property fields to show that they are suitable for statistical analyses after suitable pretreatment. We then analyze these property fields with partial least-squares (PLS) regression to predict biological affinities of two data sets comprising factor Xa and GABA-A/benzodiazepine receptor ligands. While the resulting models perform equally well or even slightly better in terms of consistency and predictivity than the classical CoMFA fields, the most important aspect of these augmented field-types is that the chemical interpretation of resulting QM-based property field models reveals unique SAR trends driven by electrostatic and polarizability effects, which cannot be extracted directly from CoMFA electrostatic maps. Within the factor Xa set, the interaction of chlorine and bromine atoms with a tyrosine side chain in the protease S1 pocket are correctly predicted. Within the GABA-A/benzodiazepine ligand data set, PLS models of high predictivity resulted for our QM-based property fields, providing novel insights into key features of the SAR for two receptor subtypes and cross-receptor selectivity of the ligands. The detailed interpretation of regression models derived using improved QM-derived property fields thus provides a significant advantage by revealing chemically meaningful correlations with biological activity and helps in understanding novel structure-activity relationship features. This will allow such knowledge to be used to design novel molecules on the basis of interactions additional to steric and hydrogen-bonding features.

The application of bioisosteres in drug design for novel drug discovery: focusing on acid protease inhibitors

INTRODUCTION

A bioisostere is a powerful concept for medicinal chemistry. It allows the improvement of the stability; oral absorption; membrane permeability; and absorption, distribution, metabolism and excretion (ADME) of drug candidate, while retaining their biological properties. The term 'bioisostere' is derived from 'isostere', whose physical and chemical properties, such as steric size, hydrophobicity, and electronegativity, are similar to those of a functional or atomic group, and is considered to possess biological properties. Here, the authors highlight the recent applications of bioisosteres in drug design, mainly based on our drug discovery studies.

AREAS COVERED

This review discusses the application of bioisosteres for novel drug discovery with focus on the authors' drug discovery studies such as renin, HIV-protease, and β-secretase inhibitors. The authors highlight that some bioisosteres can form the scaffolding for drug candidates, namely substrate transition state, amide/ester, and carboxylic acid bioisosteres. Moreover, the authors propose the new terms 'electron-donor bioisostere' and 'conformational bioisostere' for drug discovery.

EXPERT OPINION

The authors discuss the importance of bioisostere's design concept based on specific interaction with the corresponding biomolecule. In addition, some strategies for drug discovery based on the bioisostere concept are introduced. Many bioisosteres, which are recognized by corresponding target biomolecules as exhibiting similar biological properties, have been reported to date; most of the recently developed bioisosteres were designed by cheminformatics approaches. Some molecular design softwares and databases are introduced.

Designing allosteric regulators of thrombin. Monosulfated benzofuran dimers selectively interact with Arg173 of exosite 2 to induce inhibition

Earlier, we reported on the design of sulfated benzofuran dimers (SBDs) as allosteric inhibitors of thrombin (Sidhu et al. J. Med. Chem.201154 5522-5531). To identify the site of binding of SBDs, we studied thrombin inhibition in the presence of exosite 1 and 2 ligands. Whereas hirudin peptide and heparin octasaccharide did not affect the IC(50) of thrombin inhibition by a high affinity SBD, the presence of full-length heparin reduced inhibition potency by 4-fold. The presence of γ' fibrinogen peptide, which recognizes Arg93, Arg97, Arg173, Arg175, and other residues, resulted in a loss of affinity that correlated with the ideal Dixon-Webb competitive profile. Replacement of several arginines and lysines of exosite 2 with alanine did not affect thrombin inhibition potency, except for Arg173, which displayed a 22-fold reduction in IC(50). Docking studies suggested a hydrophobic patch around Arg173 as a plausible site of SBD binding to thrombin. The absence of the Arg173-like residue in factor Xa supported the observed selectivity of inhibition by SBDs. Cellular toxicity studies indicated that SBDs are essentially nontoxic to cells at concentrations as high as 250 mg/kg. Overall, the work presents the localization of the SBD binding site, which could lead to allosteric modulators of thrombin that are completely different from all clinically used anticoagulants.

The role of structural information in the discovery of direct thrombin and factor Xa inhibitors

The quest for novel medications to treat thromboembolic disorders such as venous thrombosis, pulmonary embolism and stroke received a boost when the 3D structures of two major players in the blood coagulation cascade were determined in 1989 and 1993. Structure-guided design of inhibitors of thrombin (factor IIa, fIIa) and factor Xa (fXa) eventually led to the discovery of potent, selective, efficacious, orally active and safe compounds that proved successful in clinical studies. In 2008, the direct thrombin inhibitor dabigatran etexilate developed by Boehringer Ingelheim became the first novel antithrombotic molecular entity to enter the market in 50 years. Additional compounds targeting factor Xa were subsequently granted marketing authorization or are in late-stage clinical studies. In this review, I use selected case studies to describe the discovery of novel fIIa and fXa inhibitors, with a particular emphasis on the pre-eminent role that structural information played in this process.

Site-Specific Fragment Identification Guided by Single-Step Free Energy Perturbation Calculations

The in-silico Site Identification by Ligand Competitive Saturation (SILCS) approach identifies the binding sites of representative chemical entities on the entire protein surface, information that can be applied for computational fragment-based drug design. In this study, we report an efficient computational protocol that uses sampling of the protein-fragment conformational space obtained from the SILCS simulations and performs single step free energy perturbation (SSFEP) calculations to identify site-specific favorable chemical modifications of benzene involving substitutions of ring hydrogens with individual non-hydrogen atoms. The SSFEP method is able to capture the experimental trends in relative hydration free energies of benzene analogues and for two datasets of experimental relative binding free energies of congeneric series of ligands of the proteins α-thrombin and P38 MAP kinase. The approach includes a protocol in which data obtained from SILCS simulations of the proteins is first analyzed to identify favorable benzene binding sites following which an ensemble of benzene-protein conformations for that site is obtained. The SSFEP protocol applied to that ensemble results in good reproduction of experimental free energies of the α-thrombin ligands, but not for P38 MAP kinase ligands. Comparison with results from a P38 full-ligand simulation and analysis of conformations reveals the reason for the poor agreement being the connectivity with the remainder of the ligand, a limitation inherent in fragment-based methods. Since the SSFEP approach can identify favorable benzene modifications as well as identify the most favorable fragment conformations, the obtained information can be of value for fragment linking or structure-based optimization.

Synthesis and biochemical evaluation of triazole/tetrazole-containing sulfonamides against thrombin and related serine proteases

A small library of 25 triazole/tetrazole-based sulfonamides have been synthesized and further evaluated for their inhibitory activity against thrombin, trypsin, tryptase and chymase. In general, the triazole-based sulfonamides inhibited thrombin more efficiently than the tetrazole counterparts. Particularly, compound 26 showed strong thrombin inhibition (K(i)=880 nM) and significant selectivity against other human related serine proteases like trypsin (K(i)=729 μM). Thrombin binding affinity of the same compound was determined by ITC and demonstrated that the binding of this new triazole-based scaffold is enthalpically driven, making it a good candidate for further development.

Rational design of potent, small, synthetic allosteric inhibitors of thrombin

Thrombin is a key enzyme targeted by the majority of current anticoagulants that are direct inhibitors. Allosteric inhibition of thrombin may offer a major advantage of finely tuned regulation. We present here sulfated benzofurans as the first examples of potent, small allosteric inhibitors of thrombin. A sulfated benzofuran library of 15 sulfated monomers and 13 sulfated dimers with different charged, polar, and hydrophobic substituents was studied in this work. Synthesis of the sulfated benzofurans was achieved through a multiple step, highly branched strategy, which culminated with microwave-assisted chemical sulfation. Of the 28 potential inhibitors, 11 exhibited reasonable inhibition of human α-thrombin at pH 7.4. Structure-activity relationship analysis indicated that sulfation at the 5-position of the benzofuran scaffold was essential for targeting thrombin. A tert-butyl 5-sulfated benzofuran derivative was found to be the most potent thrombin inhibitor with an IC(50) of 7.3 μM under physiologically relevant conditions. Michaelis-Menten studies showed an allosteric inhibition phenomenon. Plasma clotting assays indicate that the sulfated benzofurans prolong both the activated partial thromboplastin time and prothrombin time. Overall, this work puts forward sulfated benzofurans as the first small, synthetic molecules as powerful lead compounds for the design of a new class of allosteric inhibitors of thrombin.

Contribution of explicit solvent effects to the binding affinity of small-molecule inhibitors in blood coagulation factor serine proteases

The prevention of blood coagulation is important in treating thromboembolic disorders, and several serine proteases involved in the coagulation cascade have been classified as pharmaceutically relevant. Whereas structure-based drug design has contributed to the development of some serine protease inhibitors, traditional computational methods have not been able to fully describe structure-activity relationships (SAR). Here, we study the SAR for a number of serine proteases by using a method that calculates the thermodynamic properties (enthalpy and entropy) of the water that solvates the active site. We show that the displacement of water from specific subpockets (such as S1-4 and the ester binding pocket) of the active site by the ligand can govern potency, especially for cases in which small chemical changes (i.e., a methyl group or halogen) result in a substantial increase in potency. Furthermore, we describe how relative binding free energies can be estimated by combining the water displacement energy with complementary terms from an implicit solvent molecular mechanics description binding.

A network-based multi-target computational estimation scheme for anticoagulant activities of compounds

Background

Traditional virtual screening method pays more attention on predicted binding affinity between drug molecule and target related to a certain disease instead of phenotypic data of drug molecule against disease system, as is often less effective on discovery of the drug which is used to treat many types of complex diseases. Virtual screening against a complex disease by general network estimation has become feasible with the development of network biology and system biology. More effective methods of computational estimation for the whole efficacy of a compound in a complex disease system are needed, given the distinct weightiness of the different target in a biological process and the standpoint that partial inhibition of several targets can be more efficient than the complete inhibition of a single target.

Methodology

We developed a novel approach by integrating the affinity predictions from multi-target docking studies with biological network efficiency analysis to estimate the anticoagulant activities of compounds. From results of network efficiency calculation for human clotting cascade, factor Xa and thrombin were identified as the two most fragile enzymes, while the catalytic reaction mediated by complex IXa:VIIIa and the formation of the complex VIIIa:IXa were recognized as the two most fragile biological matter in the human clotting cascade system. Furthermore, the method which combined network efficiency with molecular docking scores was applied to estimate the anticoagulant activities of a serial of argatroban intermediates and eight natural products respectively. The better correlation (r = 0.671) between the experimental data and the decrease of the network deficiency suggests that the approach could be a promising computational systems biology tool to aid identification of anticoagulant activities of compounds in drug discovery.

Conclusions

This article proposes a network-based multi-target computational estimation method for anticoagulant activities of compounds by combining network efficiency analysis with scoring function from molecular docking.

P3 optimization of functional potency, in vivo efficacy and oral bioavailability in 3-aminopyrazinone thrombin inhibitors bearing non-charged groups at the P1 position

Although the S3 pocket of the thrombin active site is lined with lipophilic amino acid residues, the accommodation of polarity within the lipophilic P3 moiety of small molecule inhibitors is possible provided that the polar functionality is capable of pointing away from the binding pocket outwards toward solvent while simultaneously allowing the lipophilic portion of the P3 ligand to interact with the S3 amino acid residues. Manipulation of this motif provided the means to effect optimization of functional potency, in vivo antithrombotic efficacy and oral bioavailability in a series of 3-aminopyrazinone thrombin inhibitors which contained non-charged groups at the P1 position.

Design, synthesis and SAR of a series of 1,3,5-trisubstituted benzenes as thrombin inhibitors

A novel 1,3,5-trisubstituted benzamide thrombin inhibitor template was designed via hybridization of a known aminopyridinoneacetamide and a known 1,3,5-trisubstituted phenyl ether. Optimization of this lead afforded a novel potent series of biaryl 1,3,5-trisubstituted benzenes with excellent functional anticoagulant potency.

Discovery of a tetrahydropyrimidin-2(1H)-one derivative (TAK-442) as a potent, selective, and orally active factor Xa inhibitor

Coagulation enzyme factor Xa (FXa) is a particularly promising target for the development of new anticoagulant agents. We previously reported the imidazo[1,5-c]imidazol-3-one derivative 1 as a potent and orally active FXa inhibitor. However, it was found that 1 predominantly undergoes hydrolysis upon incubation with human liver microsomes, and the human specific metabolic pathway made it difficult to predict the human pharmacokinetics. To address this issue, our synthetic efforts were focused on modification of the imidazo[1,5-c]imidazol-3-one moiety of the active metabolite 3a, derived from 1, which resulted in the discovery of the tetrahydropyrimidin-2(1H)-one derivative 5k as a highly potent and selective FXa inhibitor. Compound 5k showed no detectable amide bond cleavage in human liver microsomes, exhibited a good pharmacokinetic profile in monkeys, and had a potent antithrombotic efficacy in a rabbit model without prolongation of bleeding time. Compound 5k is currently under clinical development with the code name TAK-442.

Enhancement of hydrophobic interactions and hydrogen bond strength by cooperativity: synthesis, modeling, and molecular dynamics simulations of a congeneric series of thrombin inhibitors

Accurately predicting the binding affinity of ligands to their receptors by computational methods is one of the major challenges in structure-based drug design. One of the potentially significant errors in these predictions is the common assumption that the ligand binding affinity contributions of noncovalent interactions are additive. Herein we present data obtained from two separate series of thrombin inhibitors containing hydrophobic side chains of increasing size that bind in the S3 pocket and with, or without, an adjacent amine that engages in a hydrogen bond with Gly 216. The first series of inhibitors has a m-chlorobenzyl moiety binding in the S1 pocket, and the second has a benzamidine moiety. When the adjacent hydrogen bond is present, the enhanced binding affinity per A(2) of hydrophobic contact surface in the S3 pocket improves by 75% and 59%, respectively, over the inhibitors lacking this hydrogen bond. This improvement of the binding affinity per A(2) demonstrates cooperativity between the hydrophobic interaction and the hydrogen bond.

Think twice: understanding the high potency of bis(phenyl)methane inhibitors of thrombin

Successful design of potent and selective protein inhibitors, in terms of structure-based drug design, strongly relies on the correct understanding of the molecular features determining the ligand binding to the target protein. We present a case study of serine protease inhibitors with a bis(phenyl)methane moiety binding into the S3 pocket. These inhibitors bind with remarkable potency to the active site of thrombin, the blood coagulation factor IIa. A combination of X-ray crystallography and isothermal titration calorimetry provides conclusive insights into the driving forces responsible for the surprisingly high potency of these inhibitors. Analysis of six well-resolved crystal structures (resolution 1.58-2.25 A) along with the thermodynamic data allows an explanation of the tight binding of the bis(phenyl)methane inhibitors. Interestingly, the two phenyl rings contribute to binding affinity for very different reasons - a fact that can only be elucidated by a structure-based approach. The first phenyl moiety occupies the hydrophobic S3 pocket, resulting in a mainly entropic advantage of binding. This observation is based on the displacement of structural water molecules from the S3 pocket that are observed in complexes with inhibitors that do not bind in the S3 pocket. The same classic hydrophobic effect cannot explain the enhanced binding affinity resulting from the attachment of the second, more solvent-exposed phenyl ring. For the bis(phenyl)methane inhibitors, an observed adaptive rotation of a glutamate residue adjacent to the S3 binding pocket attracted our attention. The rotation of this glutamate into salt-bridging distance with a lysine moiety correlates with an enhanced enthalpic contribution to binding for these highly potent thrombin binders. This explanation for the magnitude of the attractive force is confirmed by data retrieved by a Relibase search of several thrombin-inhibitor complexes deposited in the Protein Data Bank exhibiting similar molecular features. Special attention was attributed to putative changes in the protonation states of the interaction partners. For this purpose, two analogous inhibitors differing mainly in their potential to change the protonation state of a hydrogen-bond donor functionality were compared. Buffer dependencies of the binding enthalpy associated with complex formation could be traced by isothermal titration calorimetry, which revealed, along with analysis of the crystal structures (resolution 1.60 and 1.75 A), that a virtually compensating proton interchange between enzyme, inhibitor and buffer is responsible for the observed buffer-independent thermodynamic signatures.

More than a simple lipophilic contact: a detailed thermodynamic analysis of nonbasic residues in the s1 pocket of thrombin

The field of medicinal chemistry aims to design and optimize small molecule leads into drug candidates that may positively interfere with pathological disease situations in humans or combat the growth of infective pathogens. From the plethora of crystal structures of protein-inhibitor complexes we have learned how molecules recognize each other geometrically, but we still have rather superficial understanding of why they bind to each other. This contribution surveys a series of 26 thrombin inhibitors with small systematic structural differences to elucidate the rationale for their widely deviating binding affinity from 185 microM to 4 nM as recorded by enzyme kinetic measurements. Five well-resolved (resolution 2.30 - 1.47 A) crystal structures of thrombin-inhibitor complexes and an apo-structure of the uncomplexed enzyme (1.50 A) are correlated with thermodynamic data recorded by isothermal titration calorimetry with 12 selected inhibitors from the series. Taking solubility data into account, the variation in physicochemical properties allows conclusions to be reached about the relative importance of the enthalpic binding features as well as to estimate the importance of the parameters more difficult to capture, such as residual ligand entropy and desolvation properties. The collected data reveal a comprehensive picture of the thermodynamic signature that explains the so far poorly understood attractive force experienced by m-chloro-benzylamides to thrombin.

Design, structure-activity relationships, X-ray crystal structure, and energetic contributions of a critical P1 pharmacophore: 3-chloroindole-7-yl-based factor Xa inhibitors

An indole-based P1 moiety was incorporated into a previously established factor Xa inhibitor series. The indole group was designed to hydrogen-bond with the carbonyl of Gly218, while its 3-methyl or 3-chloro substituent was intended to interact with Tyr228. These interactions were subsequently observed in the X-ray crystal structure of compound 18. SAR studies led to the identification of compound 20 as the most potent FXa inhibitor in this series (IC(50) = 2.4 nM, EC(2xPT) = 1.2 microM). An in-depth energetic analysis suggests that the increased binding energy of 3-chloroindole-versus 3-methylindole-containing compounds in this series is due primarily to (a) the more hydrophobic nature of chloro- versus methyl-containing compounds and (b) an increased interaction of 3-chloroindole versus 3-methylindole with Gly218 backbone. The stronger hydrophobicity of chloro- versus methyl-substituted aromatics may partly explain the general preference for chloro- versus methyl-substituted P1 groups in FXa, which extends beyond the current series.

Discovery of imidazo[1,5-c]imidazol-3-ones: weakly basic, orally active factor Xa inhibitors

The coagulation enzyme factor Xa (FXa) has been recognized as a promising target for the development of new antithrombotic agents. We previously found compound 1 to be an orally bioavailable FXa inhibitor in fasted monkeys; however, 1 showed poor bioavailability in rats and fed monkeys. To work out the pharmacokinetic problems, we focused our synthetic efforts on the chemical conversion of the 4-(imidazo[1,2- a]pyridin-5-yl)piperazine moiety of 1 to imidazolylpiperidine derivatives (fused and nonfused), which resulted in the discovery of the weakly basic imidazo[1,5- c]imidazol-3-one 3q as a potent and selective FXa inhibitor. Compound 3q showed favorable oral bioavailability in rats and monkeys under both fasted and fed conditions and antithrombotic efficacy in a rat model of venous thrombosis after oral administration, without a significant increase in bleeding time (unlike warfarin). On the basis of these promising properties, compound 3q was selected for further evaluation.

Cl-pi interactions in protein-ligand complexes

During systematic analysis of nonbonded contacts in protein-ligand complexes derived from crystal structures in the Protein Data Bank, Cl-pi interactions have been found, not only in the well-documented serine proteases but also, to a lesser extent, in other proteins. From geometric analysis of such Cl-pi interactions in the crystal structures, two distinct geometries were found: the "edge-on" approach of a Cl atom to a ring atom or C-C bond and the "face-on" approach toward the ring centroid with an average interatomic distance of 3.6 A. High-level ab initio calculations using benzene-chlorohydrocarbon model systems elucidated that the calculated Cl-pi interaction energy is -2.01 kcal/mol, and the dispersion force is the major source of attraction. We also discussed the geometric flexibility in Cl-pi interactions and a relationship between the intensity of the pi density in an aromatic ring and the interaction position of the Cl atom.

Evaluating docking programs: keeping the playing field level

Over recent years many enrichment studies have been published which purport to rigorously compare the performance of two or more docking protocols. It has become clear however that such studies often have flaws within their methodologies, which cast doubt on the rigour of the conclusions. Setting up such comparisons is fraught with difficulties and no best mode of practice is available to guide the experimenter. Careful choice of structural models and ligands appropriate to those models is important. The protein structure should be representative for the target. In addition the set of active ligands selected should be appropriate to the structure in cases where different forms of the protein bind different classes of ligand. Binding site definition is also an area in which errors arise. Particular care is needed in deciding which crystallographic waters to retain and again this may be predicated by knowledge of the likely binding modes of the ligands making up the active ligand list. Geometric integrity of the ligand structures used is clearly important yet it is apparent that published sets of actives + decoys may contain sometimes high proportions of incorrect structures. Choice of protocol for docking and analysis needs careful consideration as many programs can be tweaked for optimum performance. Should studies be run using 'black box' protocols supplied by the software provider? Lastly, the correct method of analysis of enrichment studies is a much discussed topic at the moment. However currently promoted approaches do not consider a crucial aspect of a successful virtual screen, namely that a good structural diversity of hits be returned. Overall there is much to consider in the experimental design of enrichment studies. Hopefully this study will be of benefit in helping others plan such experiments.

Mechanism-Based Thrombin Inhibitors: Design, Synthesis, and Molecular Docking of a New Selective 2-Oxo-2H-1-benzopyran Derivative

New 2-oxo-2H-1-benzopyran derivatives were prepared to optimize 2a,b, initially developed as mechanism-based alpha-chymotrypsin (alpha-CT) inhibitors, into potent and selective thrombin (THR) inhibitors. From this study, 22, characterized by a 2-(N-ethyl-2'-oxoacetamide)-5'-chlorophenyl ester side chain, was shown to be a good THR inhibitor (ki/KI = 3455 M(-1) x s(-1)), displaying an excellent selectivity profile against other serine proteases such as factor Xa, trypsin, and alpha-CT. Docking analysis of this compound into the different protein structures revealed the molecular basis responsible for its potency and selectivity.

Quantum and molecular dynamics study for binding of macrocyclic inhibitors to human alpha-thrombin

Molecular dynamics simulations followed by quantum mechanical calculation and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analysis have been carried out to study binding of proline- and pyrazinone-based macrocyclic inhibitors (L86 and T76) to human alpha-thrombin. Detailed binding interaction energies between these inhibitors and individual protein fragments are calculated using DFT method based on a new quantum mechanical approach for computing protein-ligand interaction energy. The analysis of detailed interaction energies provides insight on the protein-ligand binding mechanism. Study shows that T76 and L86 bind to thrombin in a very similar "inhibition mode" except that T76 has relatively weaker binding interaction with Glu(217). The analysis from quantum calculation of binding interaction is consistent with the MM-PBSA calculation of binding free energy, and the calculated free energies for L86/T76-thrombin binding agree well with the experimental data.