Design and structure-activity relationships of potent and selective inhibitors of blood coagulation factor Xa

Rapid and Mild Lactamization Using Highly Electrophilic Triphosgene in a Microflow Reactor

Lactams are cyclic amides that are indispensable as drugs and as drug candidates. Conventional lactamization includes acid-mediated and coupling-agent-mediated approaches that suffer from narrow substrate scope, much waste, and/or high cost. Inexpensive, less-wasteful approaches mediated by highly electrophilic reagents are attractive, but there is an imminent risk of side reactions. Herein, a methods using highly electrophilic triphosgene in a microflow reactor that accomplishes rapid (0.5-10 s), mild, inexpensive, and less-wasteful lactamization are described. Methods A and B, which use N-methylmorpholine and N-methylimidazole, respectively, were developed. Various lactams and a cyclic peptide containing acid- and/or heat-labile functional groups were synthesized in good to high yields without the need for tedious purification. Undesired reactions were successfully suppressed, and the risk of handling triphosgene was minimized by the use of microflow technology.

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.

Realizing private and practical pharmacological collaboration

Although combining data from multiple entities could power life-saving breakthroughs, open sharing of pharmacological data is generally not viable because of data privacy and intellectual property concerns. To this end, we leverage modern cryptographic tools to introduce a computational protocol for securely training a predictive model of drug-target interactions (DTIs) on a pooled dataset that overcomes barriers to data sharing by provably ensuring the confidentiality of all underlying drugs, targets, and observed interactions. Our protocol runs within days on a real dataset of more than 1 million interactions and is more accurate than state-of-the-art DTI prediction methods. Using our protocol, we discover previously unidentified DTIs that we experimentally validated via targeted assays. Our work lays a foundation for more effective and cooperative biomedical research.

Quantitative Correlation of Conformational Binding Enthalpy with Substrate Specificity of Serine Proteases

Members of the same protease family show different substrate specificity, even if they share identical folds, depending on the physiological processes they are part of. Here, we investigate the key factors for subpocket and global specificity of factor Xa, elastase, and granzyme B which despite all being serine proteases and sharing the chymotrypsin-fold show distinct substrate specificity profiles. We determined subpocket interaction potentials with GRID for static X-ray structures and an in silico generated ensemble of conformations. Subpocket interaction potentials determined for static X-ray structures turned out to be insufficient to explain serine protease specificity for all subpockets. Therefore, we generated conformational ensembles using molecular dynamics simulations. We identified representative binding site conformations using distance-based hierarchical agglomerative clustering and determined subpocket interaction potentials for each representative conformation of the binding site. Considering the differences in subpocket interaction potentials for these representative conformations as well as their abundance allowed us to quantitatively explain subpocket specificity for the nonprime side for all three example proteases on a molecular level. The methods to identify key regions determining subpocket specificity introduced in this study are directly applicable to other serine proteases, and the results provide starting points for new strategies in rational drug design.

Prediction of factor Xa inhibitors by machine learning methods

Factor Xa (FXa) inhibitors have been explored as anticoagulants for treatment and prevention of thrombotic diseases. Molecular docking, pharmacophore, quantitative structure-activity relationships, and support vector machines (SVM) have been used for computer prediction of FXa inhibitors. These methods achieve promising prediction accuracies of 69-80% for FXa inhibitors and 85-99% for non-inhibitors. Prediction performance, particularly for inhibitors, may be further improved by exploring methods applicable to more diverse range of compounds and by using more appropriate set of molecular descriptors. We tested the capability of several machine learning methods (C4.5 decision tree, k-nearest neighbor, probabilistic neural network, and support vector machine) by using a much more diverse set of 1098 compounds (360 inhibitors and 738 non-inhibitors) than those in other studies. A feature selection method was used for selecting molecular descriptors appropriate for distinguishing FXa inhibitors and non-inhibitors. The prediction accuracies of these methods are 89.1-97.5% for FXa inhibitors and 92.3-98.1% for non-inhibitors. In particular, compared to other studies, support vector machine gives a substantially improved accuracy of 94.6% for FXa non-inhibitors and maintains a comparable accuracy of 98.1% for inhibitors, based-on a more rigorous test with more diverse range of compounds. Our study suggests that machine learning methods such as SVM are useful for facilitating the prediction of FXa inhibitors.

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

Factor Xa inhibitory activities for a series of N-{(3S)-1-[(1S)-1-methyl-2-morpholin-4-yl-2-oxoethyl]-2-oxopyrrolidin-3-yl}sulfonamides with different P1 groups are described. These data provide insight into binding interactions within the S1 primary specificity pocket; rationales are presented for the derived SAR on the basis of electronic interactions through crystal structures of fXa-ligand complexes and molecular modeling studies. A good correlation between in vitro anticoagulant activities with lipophilicity and the extent of human serum albumin binding is observed within this series of potent fXa inhibitors. Pharmacokinetic profiles in rat and dog, together with selectivity over other trypsin-like serine proteases, identified 1f as a candidate for further evaluation.

Design and synthesis of orally active pyrrolidin-2-one-based factor Xa inhibitors

A series of novel, non-basic 3-(6-chloronaphth-2-ylsulfonyl)aminopyrrolidin-2-one-based factor Xa (fXa) inhibitors, incorporating an alanylamide P4 group, was designed and synthesised. Within this series, the N-2-(morpholin-4-yl)-2-oxoethyl derivative 24 was shown to be a potent, selective fXa inhibitor with good anticoagulant activity. Moreover, 24 possessed highly encouraging rat and dog pharmacokinetic profiles with excellent oral bioavailabilities in both species.

Unusually Weak Binding Interactions in Tetrazole−Amidine Complexes

[reaction: see text] Tetrazoles frequently replace carboxylic acids in pharmaceutical drugs. However, while the binding modes of tetrazolate and carboxylate anions in amidinium complexes turns out to be similar, the association constant of the former is 2-3 orders of magnitude smaller in DMSO. Crystal structures revealed that the N...H-N hydrogen bonds in amidinium tetrazolates are bent (162 degrees and 169 degrees ) and noticeably longer (N...N 2.96 A) than corresponding hydrogen bonds in both amidinium carboxylates and ammonium tetrazolates.

Synthesis and aminoacyl-tRNA synthetase inhibitory activity of aspartyl adenylate analogs

Three nonhydrolyzable aspartyl adenylate analogs have been prepared and tested as inhibitors of E. coli aspartyl-tRNA synthetase. 5'-O-[N-(L-Aspartyl)sulfamoyl]adenosine is a potent competitive inhibitor (K(i) = 15 nM) whereas L-aspartol adenylate is a weaker inhibitor (K(i) = 45 microM) with respect to aspartic acid. The corresponding ketomethylphosphonate (a novel isosteric replacement) is also a strong inhibitor (K(i) = 123 nM).

Design and Synthesis of a Potent and Selective Peptidomimetic Inhibitor of Caspase-3

In this paper we report the synthesis and characterization of a novel potent and selective inhibitor of caspase-3, a member of the caspase family of cysteine proteases which plays an important role in many human disorders. This molecule represents 3(S)-acetylamino-N-[1-[(((3S)-2-hydroxy-5-oxo-tetrahydrofuran-3-yl)carbamoyl)methyl]-2-oxo-5-phenyl-2,3-dihydro-1H-benzo[e][1,4]diazepin-3-yl]succinamic acid, a monocyclic conformationally constrained form of the tetrapeptide Ac-DEVD-H, in which a 1,4-benzodiazepine nucleus is introduced internally to the peptidic sequence.

Virtual screening with flexible docking and COMBINE-based models. Application to a series of factor Xa inhibitors

A two-step, fully automatic virtual screening procedure consisting of flexible docking followed by activity prediction by COMparative BINding Energy (COMBINE) analysis is presented. This novel approach has been successfully applied, as an example with medicinal chemistry interest, to a recently reported series of 133 factor Xa (fXa)(1) inhibitors whose activities encompass 4 orders of magnitude. The docking algorithm is linked to the COMBINE analysis program and used to derive independent regression models of the 133 inhibitors docked within three different fXa structures (PDB entries 1fjs, 1f0r, and 1xka), so as to explore the effect of receptor conformation on the overall results. Reliable docking conformations and predictive regression models requiring eight latent variables could be derived for two of the fXa structures, with the best model achieving a Q(2) of 0.63 and a standard deviation of errors of prediction (SDEP) of 0.51 (leave-one-out). The two-step procedure was then employed to screen a designed virtual library of 112 ligands, containing both active and inactive compounds. While docking energies alone could show a good performance for selecting hits, including structurally diverse ones, inclusion of COMBINE analysis regression models provided improved rankings for the identification of structurally related molecules in external sets. In our best case, a recognition rate of approximately 80% of known binders at approximately 15% false positives rate was achieved, corresponding to an enrichment factor of approximately 450% over random.

Synthesis and antitumor evaluation of some new substituted amidino-benzimidazolyl-furyl-phenyl-acrylates and naphtho[2,1-b]furan-carboxylates

The multistep synthesis of a series of substituted amidino-benzimidazolyl-furyl-phenyl-acrylic acid's esters and substituted amidino-benzimidazolyl-naphtho[2,1-b]furan-carboxylic acid's esters is described starting from corresponding 3-(2-furyl)-2-phenyl-acrylic acids. The new compounds were tested on the cytostatic activities against malignant cell lines: pancreatic carcinoma (MiaPaCa2), breast carcinoma (MCF7), cervical carcinoma (HeLa), laryngeal carcinoma (Hep2), colon carcinoma (HT 29), melanoma (HBL), and human fibroblasts cell line (WI38). All compounds inhibited the proliferation of tumor cell lines. Inhibitory effect of examined compounds depended on concentration, but without significant difference among the type of tumor cells. The compounds 2 and 5 exerted very low inhibitory effect on the growth of human fibroblasts. Unsubstituted derivative 8 has not inhibited any tested cell lines.

Single-step perturbations to calculate free energy differences from unphysical reference states: limits on size, flexibility, and character

Relative free energies for a series of not too different compounds can be estimated accurately from a single simulation of an unphysical reference state that encompasses the characteristic molecular features of the compounds. Previously, this method has been applied to the calculation of free energies of solvation and of ligand binding for small molecules. In the present study we investigate the limits to the accuracy of the method by applying it to a realistic model of the binding of a set of rather large ligands to the protein factor Xa, a key protein in current efforts to design anticoagulation drugs. The evaluation of the binding free energies and conformations of nine derivatives of a biphenylamidino inhibitor leads to insights regarding the effect of the size, flexibility, and character of the unphysical part of the ligand in the reference state on the accuracy of the predicted binding free energies.

Efficacious and orally bioavailable thrombin inhibitors based on a 2,5-thienylamidine at the P1 position: discovery of N-carboxymethyl-d-diphenylalanyl-l-prolyl[(5-amidino-2-thienyl)methyl]amide

Thrombin, a crucial enzyme in the blood coagulation, has been a target for antithrombotic therapy. Orally active thrombin inhibitors would provide effective and safe prophylaxis for venous and arterial thrombosis. We conducted optimization of a highly efficacious benzamidine-based thrombin inhibitor LB30812 (3, K(i) = 3 pM) to improve oral bioavailability. Of a variety of arylamidines investigated at the P1 position, 2,5-thienylamidine effectively replaced the benzamidine without compromising the thrombin inhibitory potency and oral absorption. The sulfamide and sulfonamide derivatization at the N-terminal position in general afforded highly potent thrombin inhibitors but with moderate oral absorption, while the well-absorbable N-carbamate derivatives exhibited limited metabolic stability in S9 fractions. The present work culminated in the discovery of the N-carboxymethyl- and 2,5-thienylamidine-containing compound 22 that exhibits the most favorable profiles of anticoagulant and antithrombotic activities as well as oral bioavilability (K(i) = 15 pM; F = 43%, 42%, and 15% in rats, dogs, and monkeys, respectively). This compound on a gravimetric basis was shown to be more effective than a low molecular weight heparin, enoxaparin, in the venous thrombosis models of rat and rabbit. Compound 22 (LB30870) was therefore selected for further preclinical and clinical development.

Design, synthesis, and SAR of substituted acrylamides as factor Xa inhibitors

Substituted acrylamides were used as templates that bridge P1 and P4 binding elements, resulting in a series of potent (sub-nanomolar) and selective factor Xa inhibitors. In this template, cis-geometry of P1 and P4 ligands is highly preferred. SAR on the substituting groups, as well as on modification of P1 and P4 moieties is described. Compounds in this series show good in vivo efficacy in animal models.

Noncovalent tripeptidic thrombin inhibitors incorporating amidrazone, amine and amidine functions at P1

A series of noncovalent tripeptidic thrombin inhibitors incorporating amidrazone, amine and amidine functions at P1 was investigated. While the amidrazone and the amine series displayed limited oral absorption, the amidine series demonstrated generally good oral absorption and strong antithrombotic activity; the single-digit picomolar K(i) achieved from this series is among the best yet reported. The present work highlights the benzamidine compound 11f (LB30812) that exhibits excellent overall profiles of potency, oral absorption and antithrombotic efficacy.

Exploiting subsite S1 of trypsin-like serine proteases for selectivity: potent and selective inhibitors of urokinase-type plasminogen activator

A nonselective inhibitor of trypsin-like serine proteases, 2-(2-hydroxybiphenyl-3-yl)-1H-indole-5-carboxamidine (1) (Verner, E.; Katz, B. A.; Spencer, J.; Allen, D.; Hataye, J.; Hruzewicz, W.; Hui, H. C.; Kolesnikov, A.; Li, Y.; Luong, C.; Martelli, A.; Radika. K.; Rai, R.; She, M.; Shrader, W.; Sprengeler, P. A.; Trapp, S.; Wang, J.; Young, W. B.; Mackman, R. L. J. Med. Chem. 2001, 44, 2753-2771) has been optimized through minor structural changes on the S1 binding group to afford remarkably selective and potent inhibitors of urokinase-type plasminogen activator (uPA). The trypsin-like serine proteases(1) that comprise drug targets can be broadly categorized into two subfamilies, those with Ser190 and those with Ala190. A single-atom modification, for example, replacement of hydrogen for chlorine at the 6-position of the 5-amidinoindole P1 group on 1, generated up to 6700-fold selectivity toward the Ser190 enzymes and against the Ala190 enzymes. The larger chlorine atom displaces a water molecule (H(2)O1(S1)) that binds near residue 190 in all the complexes of 1, and related inhibitors, in uPA, thrombin, and trypsin. The water molecule, H(2)O1(S1), in both the Ser190 or Ala190 enzymes, hydrogen bonds with the amidine N1 nitrogen of the inhibitor. When it is displaced, a reduction in affinity toward the Ala190 enzymes is observed due to the amidine N1 nitrogen of the bound inhibitor being deprived of a key hydrogen-bonding partner. In the Ser190 enzymes the affinity is maintained since the serine hydroxyl oxygen O gamma(Ser190) compensates for the displaced water molecule. High-resolution crystallography provided evidence for the displacement of the water molecule and validated the design rationale. In summation, a novel and powerful method for engineering selectivity toward Ser190 proteases and against Ala190 proteases without substantially increasing molecular weight is described.

Solid-phase optimisation of achiral amidinobenzyl indoles as potent and selective factor Xa inhibitors

Starting from the achiral and potent factor Xa inhibitor 1, a new and flexible solid-phase optimisation strategy is described to reduce its cationic character. By replacing one positively charged side chain by a lipophilic substituent, a novel series of highly potent and selective achiral factor Xa inhibitors was discovered. The identified lipophilic replacements in the S4 pocket might be valuable for other approaches towards fXa inhibitors.

GRID/CPCA: a new computational tool to design selective ligands

We present a computational procedure aimed at understanding enzyme selectivity and guiding the design of drugs with respect to selectivity. It starts from a set of 3D structures of the target proteins characterized by the program GRID. In the multivariate description proposed, the variables are organized and scaled in a different way than previously published methodologies. Then, consensus principal component analysis (CPCA) is used to analyze the GRID descriptors, allowing the straightforward identification of possible modifications in the ligand to improve its selectivity toward a chosen target. As an important new feature the computational method is able to work with more than two target proteins and with several 3D structures for each protein. Additionally, the use of a 'cutout tool' allows to focus on the important regions around the active site. The method is validated for a total number of nine structures of the three homologous serine proteases thrombin, trypsin, and factor Xa. The regions identified by the method as being important for selectivity are in excellent agreement with available experimental data and inhibitor structure-activity relationships.

Factor Xa: simulation studies with an eye to inhibitor design

Factor Xa is a serine protease which activates thrombin and plays a key regulatory role in the blood-coagulation cascade. Factor Xa is at the crossroads of the extrinsic and intrinsic pathways of coagulation and, hence, has become an important target for the design of anti-thrombotics (inhibitors). It is not known to be involved in other processes than hemostasis and its binding site is different to that of other serine proteases, thus facilitating selective inhibition. The design of high-affinity selective inhibitors of factor Xa requires knowledge of the structural and dynamical characteristics of its active site. The three-dimensional structure of factor Xa was resolved by X-ray crystallography and refined at 2.2 A resolution by Padmanabhan and collaborators. In this article we present results from molecular dynamics simulations of the catalytic domain of factor Xa in aqueous solution. The simulations were performed to characterise the mobility and flexibility of the residues delimiting the unoccupied binding site of the enzyme, and to determine hydrogen bonding propensities (with protein and with solvent atoms) of those residues in the active site that could interact with a substrate or a potential inhibitor. The simulation data is aimed at facilitating the design of high-affinity selective inhibitors of factor Xa.

Solid-phase parallel synthesis of azarene pyrrolidinones as factor Xa inhibitors

A focused library (4 x 14) prepared from 4-aminopyridine and 4-, 5-, and 6-azoindole templates was synthesized using 14 polymer supported 4-amido-2,3,5,6-tetrafluorophenyl (TFP) sulfonate esters inputs. Several compounds were identified as factor Xa inhibitors (IC50< or =0.1 microM) helping to establish the SAR among these four series of azarene pyrrolidinones.

Sulfonamidopyrrolidinone factor Xa inhibitors: potency and selectivity enhancements via P-1 and P-4 optimization

Sulfonamidopyrrolidinones were previously disclosed as a selective class of factor Xa (fXa) inhibitors, culminating in the identification of RPR120844 as a potent member with efficacy in vivo. Recognizing the usefulness of the central pyrrolidinone template for the presentation of ligands to the S-1 and S-4 subsites of fXa, studies to optimize the P-1 and P-4 groups were initiated. Sulfonamidopyrrolidinones containing 4-hydroxy- and 4-aminobenzamidines were discovered to be effective inhibitors of fXa. X-ray crystallographic experiments in trypsin and molecular modeling studies suggest that our inhibitors bind by insertion of the 4-hydroxybenzamidine moiety into the S-1 subsite of the fXa active site. Of the P-4 groups examined, the pyridylthienyl sulfonamides were found to confer excellent potency and selectivity especially in combination with 4-hydroxybenzamidine. Compound 20b (RPR130737) was shown to be a potent fXa inhibitor (K(i) = 2 nM) with selectivity against structurally related serine proteinases (>1000 times). Preliminary biological evaluation demonstrates the effectiveness of this inhibitor in common assays of thrombosis in vitro (e.g. activated partial thromboplastin time) and in vivo (e.g. rat FeCl(2)-induced carotid artery thrombosis model).