Rapid development of two factor IXa inhibitors from hit to lead

Molecular Dynamics Simulation Study of the Selective Inhibition of Coagulation Factor IXa over Factor Xa

Thromboembolic disorders, arising from abnormal coagulation, pose a significant risk to human life in the modern world. The FDA has recently approved several anticoagulant drugs targeting factor Xa (FXa) to manage these disorders. However, these drugs have potential side effects, leading to bleeding complications in patients. To mitigate these risks, coagulation factor IXa (FIXa) has emerged as a promising target due to its selective regulation of the intrinsic pathway. Due to the high structural and functional similarities of these coagulation factors and their inhibitor binding modes, designing a selective inhibitor specifically targeting FIXa remains a challenging task. The dynamic behavior of protein-ligand interactions and their impact on selectivity were analyzed using molecular dynamics simulation, considering the availability of potent and selective compounds for both coagulation factors and the co-crystal structures of protein-ligand complexes. Throughout the simulations, we examined ligand movements in the binding site, as well as the contact frequencies and interaction fingerprints, to gain insights into selectivity. Interaction fingerprint (IFP) analysis clearly highlights the crucial role of strong H-bond formation between the ligand and D189 and A190 in the S1 subsite for FIXa selectivity, consistent with our previous study. This dynamic analysis also reveals additional FIXa-specific interactions. Additionally, the absence of polar interactions contributes to the selectivity for FXa, as observed from the dynamic profile of interactions. A contact frequency analysis of the protein-ligand complexes provides further confirmation of the selectivity criteria for FIXa and FXa, as well as criteria for binding and activity. Moreover, a ligand movement analysis reveals key interaction dynamics that highlight the tighter binding of selective ligands to the proteins compared to non-selective and inactive ligands.

HydraMap v.2: Prediction of Hydration Sites and Desolvation Energy with Refined Statistical Potentials

The complex network of water molecules within the binding pocket of a target protein undergoes alterations upon ligand binding, presenting a significant challenge for conventional molecular modeling methods to accurately characterize and compute the associated energy changes. We have previously developed an empirical method, HydraMap (J. Chem. Inf. Model. 2020, 60, 4359-4375), which employs statistical potentials to predict hydration sites and compute desolvation energy, achieving a reasonable balance between accuracy and speed. In this work, we present its improved version, namely, HydraMap v.2. We updated the statistical potentials for protein-water interactions through an analysis of 17 042 crystal protein structures. We also introduced a new feature to evaluate ligand-water interactions by incorporating statistical potentials derived from the solvated structures of 9878 small organic molecules produced by molecular dynamics simulations. By combining these potentials, HydraMap v.2 can predict and compare the hydration sites in a binding pocket before and after ligand binding, identifying key water molecules involved in the binding process, such as those forming bridging hydrogen bonds and unstable ones that can be replaced. We demonstrated the application of HydraMap v.2 in explaining the structure-activity relationship of a panel of MCL-1 inhibitors. The desolvation energies calculated by summing the energy change of each hydration site before and after ligand binding showed good correlation with known ligand binding affinities on six target proteins. In conclusion, HydraMap v.2 offers a cost-effective solution for estimating the desolvation energy during protein-ligand binding and also is practical in guiding lead optimization in structure-based drug discovery.

Cloning and characterization of Aedes aegypti blood downregulated chymotrypsin II

Aedes aegypti adult and larval blood downregulated chymotrypsin II was cloned, sequenced and its 3D conformation modeled. Cloning of the enzymes from adult and larval guts indicated that both genes sit at the same location on Chromosome 2. Genomic analyses showed that larval and adult genes are the same and both have four exons and three introns that are located on an 8.32 Kb DNA in direction with the Ae. aegypti genome. The adult and larval transcript synthesis is controlled by alternative splicing explaining small difference in the amino acids sequences. Chymotrypsin II that was extracted from guts of sugar-fed and at 48 after blood feeding showed a pH optimum of 4-5 with a broad shoulder of activity from pH 6 to 10. Dot blot analyses show that the enzyme's transcript is downregulated after females take a blood meal and upregulated at 48 h after the blood meal. A Chymotrypsin II transcript was also detected in the larval gut during different times of larval developmental stages, indication that Ae. aegypti chymotrypsin II is synthesized by adults and larval guts. The possibility that JH III and 20HE play an active role in the regulation is discussed.

A Structure Based Study of Selective Inhibition of Factor IXa over Factor Xa

Blood coagulation is an essential physiological process for hemostasis; however, abnormal coagulation can lead to various potentially fatal disorders, generally known as thromboembolic disorders, which are a major cause of mortality in the modern world. Recently, the FDA has approved several anticoagulant drugs for Factor Xa (FXa) which work via the common pathway of the coagulation cascade. A main side effect of these drugs is the potential risk for bleeding in patients. Coagulation Factor IXa (FIXa) has recently emerged as the strategic target to ease these risks as it selectively regulates the intrinsic pathway. These aforementioned coagulation factors are highly similar in structure, functional architecture, and inhibitor binding mode. Therefore, it remains a challenge to design a selective inhibitor which may affect only FIXa. With the availability of a number of X-ray co-crystal structures of these two coagulation factors as protein–ligand complexes, structural alignment, molecular docking, and pharmacophore modeling were employed to derive the relevant criteria for selective inhibition of FIXa over FXa. In this study, six ligands (three potent, two selective, and one inactive) were selected for FIXa inhibition and six potent ligands (four FDA approved drugs) were considered for FXa. The pharmacophore hypotheses provide the distribution patterns for the principal interactions that take place in the binding site. None of the pharmacophoric patterns of the FXa inhibitors matched with any of the patterns of FIXa inhibitors. Based on pharmacophore analysis, a selectivity of a ligand for FIXa over FXa may be defined quantitatively as a docking score of lower than −8.0 kcal/mol in the FIXa-grids and higher than −7.5 kcal/mol in the FXa-grids.

Probing activation-driven changes in coagulation factor IX by mass spectrometry

BACKGROUND

Activated factor IX (FIXa) is an inefficient enzyme that needs activated factor VIII (FVIII) for full activity. Recently, we identified a network of FVIII-driven changes in FIXa employing hydrogen-deuterium eXchange mass spectrometry (HDX-MS). Some changes also occurred in active-site inhibited FIXa, but others were not cofactor-driven, in particular those within the 220-loop (in chymotrypsin numbering).

OBJECTIVE

The aim of this work is to better understand the zymogen-to-enzyme transition in FIX, with specific focus on substrate-driven changes at the catalytic site.

METHODS

Footprinting mass spectrometry by HDX and Tandem-Mass Tags (TMT) labelling were used to explore changes occurring upon the conversion from FIX into FIXa. Mutagenesis and kinetic studies served to assess the role of the 220-loop.

RESULTS

HDX-MS displayed remarkably few differences between FIX and FIXa. In comparison with FIX, FIXa did exhibit decreased deuterium uptake at the N-terminus region. This was more prominent when the FIXa active site was occupied by an irreversible inhibitor. TMT-labelling showed that the N-terminus is largely protected from labelling, and that inhibitor binding increases protection to a minor extent. Occupation of the active site also reduced deuterium uptake within the 220-loop backbone. Mutagenesis within the 220-loop revealed that a putative H-bond network contributes to FIXa activity. TMT labeling of the N-terminus suggested that these 220-loop variants are more zymogen-like than wild-type FIXa.

CONCLUSION

In the absence of cofactor and substrate, FIXa is predominantly zymogen-like. Stabilization in its enzyme-like form involves, apart from FVIII-binding, also interplay between the 220-loop, N-terminus, and the substrate binding site.

Computer Design of Low-Molecular-Weight Inhibitors of Coagulation Factors

The review discusses main approaches to searching for new low-molecular-weight inhibitors of coagulation factors IIa, Xa, IXa, and XIa and the results of such studies conducted from 2015 to 2018. For each of these factors, several inhibitors with IC₅₀ < 10 nM have been found, some of which are now tested in clinical trials. However, none of the identified inhibitors meets the requirements for an "ideal" anticoagulant, so further studies are required.

Discovery of FIXa inhibitors by combination of pharmacophore modeling, molecular docking, and 3D-QSAR modeling

Human Coagulation Factor IXa (FIXa), specifically inhibited at the initiation stage of the blood coagulation cascade, is an excellent target for developing selective and safe anticoagulants. To explore this inhibitory mechanism, 86 FIXa inhibitors were selected to generate pharmacophore models and subsequently SAR models. Both best pharmacophore model and ROC curve were built through the Receptor-Ligand Pharmacophore Generation module. CoMFA model based on molecular docking and PLS factor analysis methods were developed. Model propagations values are q² = 0.709, r² = 0.949, and r²_pred = 0.905. The satisfactory q² value of 0.609, r² value of 0.962, and r²_pred value of 0.819 for CoMSIA indicated that the CoMFA and CoMSIA models are both available to predict the inhibitory activity on FIXa. On the basis of pharmacophore modeling, molecular docking, and 3D-QSAR modeling screening, six molecules are screened as potential FIXa inhibitors.

Releasing the brakes in coagulation Factor IXa by co-operative maturation of the substrate-binding site

Coagulation Factor IX is positioned at the merging point of the intrinsic and extrinsic blood coagulation cascades. Factor IXa (activated Factor IX) serves as the trigger for amplification of coagulation through formation of the so-called Xase complex, which is a ternary complex of Factor IXa, its substrate Factor X and the cofactor Factor VIIIa on the surface of activated platelets. Within the Xase complex the substrate turnover by Factor IXa is enhanced 200000-fold; however, the mechanistic and structural basis for this dramatic enhancement remains only partly understood. A multifaceted approach using enzymatic, biophysical and crystallographic methods to evaluate a key set of activity-enhanced Factor IXa variants has demonstrated a delicately balanced bidirectional network. Essential molecular interactions across multiple regions of the Factor IXa molecule co-operate in the maturation of the active site. This maturation is specifically facilitated by long-range communication through the Ile212–Ile213 motif unique to Factor IXa and a flexibility of the 170-loop that is further dependent on the conformation in the Cys168–Cys182 disulfide bond. Ultimately, the network consists of compensatory brakes (Val16 and Ile213) and accelerators (Tyr99 and Phe174) that together allow for a subtle fine-tuning of enzymatic activity.

Definitive Metabolite Identification Coupled with Automated Ligand Identification System (ALIS) Technology: A Novel Approach to Uncover Structure-Activity Relationships and Guide Drug Design in a Factor IXa Inhibitor Program

A potent and selective Factor IXa (FIXa) inhibitor was subjected to a series of liver microsomal incubations, which generated a number of metabolites. Using automated ligand identification system-affinity selection (ALIS-AS) methodology, metabolites in the incubation mixture were prioritized by their binding affinities to the FIXa protein. Microgram quantities of the metabolites of interest were then isolated through microisolation analytical capabilities, and structurally characterized using MicroCryoProbe heteronuclear 2D NMR techniques. The isolated metabolites recovered from the NMR experiments were then submitted directly to an in vitro FIXa enzymatic assay. The order of the metabolites' binding affinity to the Factor IXa protein from the ALIS assay was completely consistent with the enzymatic assay results. This work showcases an innovative and efficient approach to uncover structure-activity relationships (SARs) and guide drug design via microisolation-structural characterization and ALIS capabilities.

Development of a novel tricyclic class of potent and selective FIXa inhibitors

Using structure based drug design, a novel class of potent coagulation factor IXa (FIXa) inhibitors was designed and synthesized. High selectivity over FXa inhibition was achieved. Selected compounds were evaluated in rat IV/PO pharmacokinetic (PK) studies and demonstrated desirable oral PK profiles. Finally, the pharmacodynamics (PD) of this class of molecules were evaluated in thrombin generation assay (TGA) in Corn Trypsin Inhibitor (CTI) citrated human plasma and demonstrated characteristics of a FIXa inhibitor.