Factor Xa Inhibitors: S1 Binding Interactions of a Series of N-{(3S)-1-[(1S)-1-Methyl-2-morpholin-4-yl-2-oxoethyl]-2-oxopyrrolidin-3-yl}sulfonamides

Toward non-factor therapy in hemophilia: an antithrombin insensitive Gla-domainless factor Xa as tissue factor pathway inhibitor bait

Background

Gla-domainless factor (F) Xa (GD-FXa) was proposed as a trap to endogenous anticoagulant tissue factor pathway inhibitor (TFPI) to restore thrombin generation in hemophilia. Using computational chemistry and experimental approaches, we previously showed that S195A GD-FXa also binds TFPI and restores ex vivo coagulation in plasma obtained from person(s) with hemophilia.

Methods

To design a GD-FXa variant with improved anti-TFPI affinity, we performed molecular dynamics simulations and identified suitable sites for mutagenesis. The calculations identified residues R150FXa and K96Fxa as cold-spots of interaction between GD-FXa and the K2 domain of TFPI. In the three-dimensional model, both residues face toward TFPI hydrophobic residues and are thus potential candidates for mutagenesis into hydrophobic residues to favor an improved protein-protein interaction.

Results

Catalytically inactive GD-FXa variants containing the S195A mutation and the additional mutations K96Y, R150I, R150G, R150F, and K96YR150F, were produced to experimentally confirm these computational hypotheses. Among these mutants, the R150FFXa and the K96YR150FFXa were slightly more effective than S195A GD-FXa in restoring coagulation in FVIII deficient plasmas. However, in surface plasmon resonance experiments, they showed TFPI binding affinities in the same range and acted similarly as S195A GD-FXa in FXa/TFPI competition assays. In contrast, the R150 mutants completely lost their interactions with antithrombin as observed in the surface plasmon resonance experiments.

Conclusions

We therefore conclude that their antithrombin resistance is responsible for their improved thrombin generation, through an extension of their half-lives.

Unified Approach to Benzo[d]thiazol-2-yl-Sulfonamides

In this paper, we report a unified approach to N-substituted and N,N-disubstituted benzothiazole (BT) sulfonamides. Our approach to BT-sulfonamides starts from simple commercially available building blocks (benzo[d]thiazole-2-thiol and primary and secondary amines) that are connected via (a) a S oxidation/S-N coupling approach, (b) a S-N coupling/S-oxidation sequence, or via (c) a S-oxidation/S-F bond formation/SuFEx approach. The labile N-H bond in N-monoalkylated BT-sulfonamides (pK_a (BTSO₂N(H)Bn) = 3.34 ± 0.05) further allowed us to develop a simple weak base-promoted N-alkylation method and a stereoselective microwave-promoted Fukuyama-Mitsunobu reaction. N-Alkyl-N-aryl BT-sulfonamides were accessed with the help of the Chan-Lam coupling reaction. Developed methods were further used in stereo and chemoselective transformations of podophyllotoxin and several amino alcohols.

Sensitivity of Binding Free Energy Calculations to Initial Protein Crystal Structure

Binding free energy calculations using alchemical free energy (AFE) methods are widely considered to be the most rigorous tool in the computational drug discovery arsenal. Despite this, the calculations suffer from accuracy, precision, and reproducibility issues. In this publication, we perform a high-throughput study of more than a thousand AFE calculations, utilizing over 220 μs of total sampling time, on three different protein systems to investigate the impact of the initial crystal structure on the resulting binding free energy values. We also consider the influence of equilibration time and discover that the initial crystal structure can have a significant effect on free energy values obtained at short timescales that can manifest itself as a free energy difference of more than 1 kcal/mol. At longer timescales, these differences are largely overtaken by important rare events, such as torsional ligand motions, typically resulting in a much higher uncertainty in the obtained values. This work emphasizes the importance of rare event sampling and long-timescale dynamics in free energy calculations even for routinely performed alchemical perturbations. We conclude that an optimal protocol should not only concentrate computational resources on achieving convergence in the alchemical coupling parameter (λ) space but also on longer simulations and multiple repeats.

Tailor-Made Amino Acids and Fluorinated Motifs as Prominent Traits in Modern Pharmaceuticals

Structural analysis of modern pharmaceutical practices allows for the identification of two rapidly growing trends: the introduction of tailor-made amino acids and the exploitation of fluorinated motifs. Curiously, the former represents one of the most ubiquitous classes of naturally occurring compounds, whereas the latter is the most xenobiotic and comprised virtually entirely of man-made derivatives. Herein, 39 selected compounds, featuring both of these traits in the same molecule, are profiled. The total synthesis, source of the corresponding amino acids and fluorinated residues, and medicinal chemistry aspects and biological properties of the molecules are discussed.

The selectivity and bioavailability improvement of novel oral anticoagulants: An overview

Anticoagulants have exhibited a critical role in the prevention and/or treatment of thrombotic diseases. Up to now, kinds of novel oral anticoagulants, inhibiting plasma serine proteases in the coagulation cascade, have been developed to overcome the clinical limitations of classical anticoagulants (like warfarin and heparins). Some of them, such as Apixaban, Rivaroxaban, Edoxaban, and Dabigatran, have been approved by FDA in recent years. This review summarizes the discovery and optimization of representative novel oral anticoagulants with the aim to improve selectivity and bioavailability of compounds. The impact of different targets in the cascade on bleeding risk also is discussed. We hope some more effective, selective, and safer anticoagulants can be developed in the future on the basis of these design experiences.

Theoretical and experimental characterization of 1,4-N⋯S σ-hole intramolecular interactions in bioactive N-acylhydrazone derivatives

Sigma-hole (σ-hole) bonds are interactions that are gaining special attention in medicinal chemistry.

CSAR 2014: A Benchmark Exercise Using Unpublished Data from Pharma

The 2014 CSAR Benchmark Exercise was the last community-wide exercise that was conducted by the group at the University of Michigan, Ann Arbor. For this event, GlaxoSmithKline (GSK) donated unpublished crystal structures and affinity data from in-house projects. Three targets were used: tRNA (m1G37) methyltransferase (TrmD), Spleen Tyrosine Kinase (SYK), and Factor Xa (FXa). A particularly strong feature of the GSK data is its large size, which lends greater statistical significance to comparisons between different methods. In Phase 1 of the CSAR 2014 Exercise, participants were given several protein-ligand complexes and asked to identify the one near-native pose from among 200 decoys provided by CSAR. Though decoys were requested by the community, we found that they complicated our analysis. We could not discern whether poor predictions were failures of the chosen method or an incompatibility between the participant's method and the setup protocol we used. This problem is inherent to decoys, and we strongly advise against their use. In Phase 2, participants had to dock and rank/score a set of small molecules given only the SMILES strings of the ligands and a protein structure with a different ligand bound. Overall, docking was a success for most participants, much better in Phase 2 than in Phase 1. However, scoring was a greater challenge. No particular approach to docking and scoring had an edge, and successful methods included empirical, knowledge-based, machine-learning, shape-fitting, and even those with solvation and entropy terms. Several groups were successful in ranking TrmD and/or SYK, but ranking FXa ligands was intractable for all participants. Methods that were able to dock well across all submitted systems include MDock,1 Glide-XP,2 PLANTS,3 Wilma,4 Gold,5 SMINA,6 Glide-XP2/PELE,7 FlexX,8 and MedusaDock.9 In fact, the submission based on Glide-XP2/PELE7 cross-docked all ligands to many crystal structures, and it was particularly impressive to see success across an ensemble of protein structures for multiple targets. For scoring/ranking, submissions that showed statistically significant achievement include MDock1 using ITScore1,10 with a flexible-ligand term,11 SMINA6 using Autodock-Vina,12,13 FlexX8 using HYDE,14 and Glide-XP2 using XP DockScore2 with and without ROCS15 shape similarity.16 Of course, these results are for only three protein targets, and many more systems need to be investigated to truly identify which approaches are more successful than others. Furthermore, our exercise is not a competition.

OPLS3: A Force Field Providing Broad Coverage of Drug-like Small Molecules and Proteins

The parametrization and validation of the OPLS3 force field for small molecules and proteins are reported. Enhancements with respect to the previous version (OPLS2.1) include the addition of off-atom charge sites to represent halogen bonding and aryl nitrogen lone pairs as well as a complete refit of peptide dihedral parameters to better model the native structure of proteins. To adequately cover medicinal chemical space, OPLS3 employs over an order of magnitude more reference data and associated parameter types relative to other commonly used small molecule force fields (e.g., MMFF and OPLS_2005). As a consequence, OPLS3 achieves a high level of accuracy across performance benchmarks that assess small molecule conformational propensities and solvation. The newly fitted peptide dihedrals lead to significant improvements in the representation of secondary structure elements in simulated peptides and native structure stability over a number of proteins. Together, the improvements made to both the small molecule and protein force field lead to a high level of accuracy in predicting protein-ligand binding measured over a wide range of targets and ligands (less than 1 kcal/mol RMS error) representing a 30% improvement over earlier variants of the OPLS force field.

Integration of Ligand and Structure Based Approaches for CSAR-2014

The prediction of binding poses and affinities is an area of active interest in computer-aided drug design (CADD). Given the documented limitations with either ligand or structure based approaches, we employed an integrated approach and developed a rapid protocol for binding mode and affinity predictions. This workflow was applied to the three protein targets of Community Structure-Activity Resource-2014 (CSAR-2014) exercise: Factor Xa (FXa), Spleen Tyrosine Kinase (SYK), and tRNA (guanine-N(1))-methyltransferase (TrmD). Our docking and scoring workflow incorporates compound clustering and ligand and protein structure based pharmacophore modeling, followed by local docking, minimization, and scoring. While the former part of the protocol ensures high-quality ligand alignments and mapping, the subsequent minimization and scoring provides the predicted binding modes and affinities. We made blind predictions of docking pose for 1, 5, and 14 ligands docked into 1, 2, and 12 crystal structures of FXa, SYK, and TrmD, respectively. The resulting 174 poses were compared with cocrystallized structures (1, 5, and 14 complexes) made available at the end of CSAR. Our predicted poses were related to the experimentally determined structures with a mean root-mean-square deviation value of 3.4 Å. Further, we were able to classify high and low affinity ligands with the area under the curve values of 0.47, 0.60, and 0.69 for FXa, SYK, and TrmD, respectively, indicating the validity of our approach in at least two of the three systems. Detailed critical analysis of the results and CSAR methodology ranking procedures suggested that a straightforward application of our workflow has limitations, as some of the performance measures do not reflect the actual utility of pose and affinity predictions in the biological context of individual systems.

Arylsulfonamidopiperidone derivatives as a novel class of factor Xa inhibitors

The design, synthesis and SAR of a novel class of valerolactam-based arylsulfonamides as potent and selective FXa inhibitors is reported. The arylsulfonamide-valerolactam scaffold was derived based on the proposed bioisosterism to the arylcyanoguanidine-caprolactam core in known FXa inhibitors. The SAR study led to compound 46 as the most potent FXa inhibitor in this series, with an IC(50) of 7 nM and EC(2×PT) of 1.7 μM. The X-ray structure of compound 40 bound to FXa shows that the sulfonamide-valerolactam scaffold anchors the aryl group in the S1 and the novel acylcytisine pharmacophore in the S4 pockets.

Structure and property based design of factor Xa inhibitors: pyrrolidin-2-ones with aminoindane and phenylpyrrolidine P4 motifs

The rational design, syntheses and evaluation of potent sulfonamidopyrrolidin-2-one-based factor Xa inhibitors incorporating aminoindane and phenylpyrrolidine P4 motifs are described. These series delivered highly potent anticoagulant compounds with excellent oral pharmacokinetic profiles; however, significant time dependant P450 inhibition was an issue for the aminoindane series, but this was not observed with the phenylpyrrolidine motif, which produced candidate quality molecules with potential for once-daily oral dosing in humans.

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.

Novel factor Xa inhibitors: a patent review

IMPORTANCE OF THE FIELD

New oral anticoagulants with favorable safety profiles and fixed doses are required for the management of thromboembolism and stroke prevention in patients with atrial fibrillation. Among them, fXa inhibitors (the so-called xabans) are attractive options that can overcome limitations (e.g., bleeding) of the current oral antithrombotic therapy. The rational design of small-molecule direct fXa inhibitors, whose importance is testified by the growing number of publications and patents recently registered, has been fully supported by the X-ray crystallography of enzyme-ligand complexes.

AREAS COVERED IN THIS REVIEW

Pubmed, SciFinder Scholar, ISI web of knowledge(SM), http://ep.espacenet.com/ and Google websites were used as the main sources for literature retrieving, and > 100 patents filed between 2006 and April 2009, reviewed and discussed herein, highlight the variety among the P1 and P4 moieties on suitable scaffolds.

WHAT THE READER WILL GAIN

The replacement of the benzamidine P1 moiety, which characterizes the first generation, with less basic bioisosteric or nonpolar neutral P1 groups led to the disclosure of numerous fXa inhibitors with high potency, selectivity and oral bioavailability. Novel selective fXa inhibitors with stable pharmacokinetics, better therapeutic windows and ease-of-use than the existing anticoagulants are currently under advanced stage clinical trials.

TAKE-HOME MESSAGE

Available data from Phase II and Phase III studies reflect the drive towards fXa inhibitors as potentially more effective and safer antithrombotic drugs. Their development is expected to address two major needs for anticoagulation, namely safety and ease-of-use, and to significantly affect the anticoagulant market.

Structure and property based design of factor Xa inhibitors: pyrrolidin-2-ones with monoaryl P4 motifs

Structure and property based drug design was exploited in the synthesis of sulfonamidopyrrolidin-2-one-based factor Xa inhibitors, incorporating neutral and basic monoaryl P4 groups, ultimately producing potent inhibitors with effective levels of anticoagulant activity and extended oral pharmacokinetic profiles. However, time dependant inhibition of Cytochrome P450 3A4 was a particular issue with this series.

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.

Antithrombotic potential of GW813893: a novel, orally active, active-site directed factor Xa inhibitor

BACKGROUND

Factor Xa (FXa) has been a target of considerable interest for drug development efforts aimed at suppressing thrombosis. In this report, a new orally active, small molecule, active-site directed FXa inhibitor, GW813893, has been profiled in a succession of in vitro and in vivo assays involved in its preclinical characterization as a potential antithrombotic therapeutic.

METHODS

In vitro profiling of GW813893 consisted of assessing its inhibitory potential against FXa and a broad panel of related and unrelated enzymes and receptors. Additionally, the FXa inhibition potential of GW813893 was assessed in prothrombinase and plasma-based clotting assays. In vivo characterization of GW813893 consisted of thrombosis studies in a rat inferior vena cava model, a rat carotid artery thrombosis model, and a rabbit jugular thrombosis model. Bleeding studies were conducted in a rat tail transection model. Ex vivo determinations of compound effects on FX and clotting activity were also undertaken.

RESULTS

GW813893 was more than 90-fold selective over all enzymes tested, and it inhibited FXa and prothrombinase activity with a Ki of 4.0 nM and 9.7 nM, respectively. In vivo, GW813893 concentration-dependently suppressed thrombotic activity in all models tested. The antithrombotic activity correlated with the suppression of plasma-based clotting activity and the inhibition of plasma FX activity (P < 0.02). Over the antithrombotic dose-range, an increased bleeding diathesis was not observed.

CONCLUSION

These experiments demonstrate that GW813893 is a potent, selective, orally active inhibitor of FXa. The data suggest that GW813893 has robust antithrombotic potential at doses that have no detectable hemostasis liability. Collectively, the profile suggests that GW813893 has the preclinical pharmacology underpinnings of an oral antithrombotic therapeutic.

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.

Selective and dual action orally active inhibitors of thrombin and factor Xa

The synthetic entry to new classes of dual fXa/thrombin and selective thrombin inhibitors with significant oral bioavailability is described. This was achieved through minor modifications to the sulfonamide group in our potent and selective fXa inhibitor (E)-2-(5-chlorothien-2-yl)-N-{(3S)-1-[(1S)-1-methyl-2-(morpholin-4-yl)-2-oxoethyl]-2-oxopyrrolidin-3-yl}ethenesulfonamide and these observed activity changes have been rationalised using structural studies.