Shedding Light on Important Waters for Drug Design: Simulations versus Grid-Based Methods.

Water molecules play an important role in the association of drugs with their pharmaceutical targets. For this reason, calculating the energetic contribution of water is essential to make accurate predictions of compounds' affinity and selectivity. Water molecules can also modify the binding mode of compounds by forming water bridges, or clusters, that stabilize a particular orientation of the ligand. Several computational methods have been developed for solvent mapping, but few studies have attempted to compare them in a drug design context. In this paper, four commercially available solvent mapping tools (SZMAP, WaterFLAP, 3D-RISM, and WaterMap) are evaluated on three different protein targets. The methods were compared by looking at their ability to predict the structure-activity relations of lead compounds. All methods were found to be useful to some degree and to improve the predictions from docking alone. However, the only simulation-based approach tested, WaterMap, was found in some cases to be more accurate than grid-based methods.

Discovery of GLPG1972/S201086, a Potent, Selective, and Orally Bioavailable ADAMTS-5 Inhibitor for the Treatment of Osteoarthritis.

There are currently no approved disease-modifying osteoarthritis (OA) drugs (DMOADs). The aggrecanase ADAMTS-5 is key in the degradation of human aggrecan (AGC), a component of cartilage. Therefore, ADAMTS-5 is a promising target for the identification of DMOADs. We describe the discovery of GLPG1972/S201086, a potent and selective ADAMTS-5 inhibitor obtained by optimization of a promising hydantoin series following an HTS. Biochemical activity against rat and human ADAMTS-5 was assessed via a fluorescence-based assay. ADAMTS-5 inhibitory activity was confirmed with human aggrecan using an AGC ELISA. The most promising compounds were selected based on reduction of glycosaminoglycan release after interleukin-1 stimulation in mouse cartilage explants and led to the discovery of GLPG1972/S201086. The anticatabolic activity was confirmed in mouse cartilage explants (IC50 < 1.5 µM). The cocrystal structure of GLPG1972/S201086 with human recombinant ADAMTS-5 is discussed. GLPG1972/S201086 has been investigated in a phase 2 clinical study in patients with knee OA (NCT03595618).

Discovery and Optimization of Orally Bioavailable Phthalazone and Cinnolone Carboxylic Acid Derivatives as S1P2 Antagonists against Fibrotic Diseases.

Idiopathic pulmonary fibrosis (IPF) is a chronic and progressive lung disease. Current treatments only slow down disease progression, making new therapeutic strategies compelling. Increasing evidence suggests that S1P2 antagonists could be effective agents against fibrotic diseases. Our compound collection was mined for molecules possessing substructure features associated with S1P2 activity. The weakly potent indole hit 6 evolved into a potent phthalazone series, bearing a carboxylic acid, with the aid of a homology model. Suboptimal pharmacokinetics of a benzimidazole subseries were improved by modifications targeting potential interactions with transporters, based on concepts deriving from the extended clearance classification system (ECCS). Scaffold hopping, as a part of a chemical enablement strategy, permitted the rapid exploration of the position adjacent to the carboxylic acid. Compound 38, with good pharmacokinetics and in vitro potency, was efficacious at 10 mg/kg BID in three different in vivo mouse models of fibrotic diseases in a therapeutic setting.

Identification of a novel NR2B-selective NMDA receptor antagonist using a virtual screening approach.

We report the identification of a novel NR2B-selective NMDAR antagonist with an original scaffold, LSP10-0500. This compound was identified by a virtual high-throughput screening approach on the basis of a quantitative pharmacophore model of NR2B-specific NMDAR antagonists. A SAR study around LSP10-0500 is also described.

Rose Bengal analogs and vesicular glutamate transporters (VGLUTs).

Vesicular glutamate transporters (VGLUTs) allow the loading of presynaptic glutamate vesicles and thus play a critical role in glutamatergic synaptic transmission. Rose Bengal (RB) is the most potent known VGLUT inhibitor (Ki 25 nM); therefore we designed, synthesized and tested in brain preparations, a series of analogs based on this scaffold. We showed that among the two tautomers of RB, the carboxylic and not the lactonic form is active against VGLUTs and generated a pharmacophore model to determine the minimal structure requirements. We also tested RB specificity in other neurotransmitter uptake systems. RB proved to potently inhibit VMAT (Ki 64 nM) but weakly VACHT (Ki>9.7 microM) and may be a useful tool in glutamate/acetylcholine co-transmission studies.

Discovery, Structure-Activity Relationship, and Binding Mode of an Imidazo[1,2-a]pyridine Series of Autotaxin Inhibitors.

Autotaxin (ATX) is a secreted enzyme playing a major role in the production of lysophosphatidic acid (LPA) in blood through hydrolysis of lysophosphatidyl choline (LPC). The ATX-LPA signaling axis arouses a high interest in the drug discovery industry as it has been implicated in several diseases including cancer, fibrotic diseases, and inflammation, among others. An imidazo[1,2-a]pyridine series of ATX inhibitors was identified out of a high-throughput screening (HTS). A cocrystal structure with one of these compounds and ATX revealed a novel binding mode with occupancy of the hydrophobic pocket and channel of ATX but no interaction with zinc ions of the catalytic site. Exploration of the structure-activity relationship led to compounds displaying high activity in biochemical and plasma assays, exemplified by compound 40. Compound 40 was also able to decrease the plasma LPA levels upon oral administration to rats.

Discovery of 2-[[2-Ethyl-6-[4-[2-(3-hydroxyazetidin-1-yl)-2-oxoethyl]piperazin-1-yl]-8-methylimidazo[1,2-a]pyridin-3-yl]methylamino]-4-(4-fluorophenyl)thiazole-5-carbonitrile (GLPG1690), a First-in-Class Autotaxin Inhibitor Undergoing Clinical Evaluation for the Treatment of Idiopathic Pulmonary Fibrosis.

Autotaxin is a circulating enzyme with a major role in the production of lysophosphatic acid (LPA) species in blood. A role for the autotaxin/LPA axis has been suggested in many disease areas including pulmonary fibrosis. Structural modifications of the known autotaxin inhibitor lead compound 1, to attenuate hERG inhibition, remove CYP3A4 time-dependent inhibition, and improve pharmacokinetic properties, led to the identification of clinical candidate GLPG1690 (11). Compound 11 was able to cause a sustained reduction of LPA levels in plasma in vivo and was shown to be efficacious in a bleomycin-induced pulmonary fibrosis model in mice and in reducing extracellular matrix deposition in the lung while also reducing LPA 18:2 content in bronchoalveolar lavage fluid. Compound 11 is currently being evaluated in an exploratory phase 2a study in idiopathic pulmonary fibrosis patients.

Reinvestigation of the Branimycin Stereochemistry at Position 17-C.

A conformational study of branimycin was performed using single-crystal X-ray crystallography to characterize the solid-state form, while a combination of NMR spectroscopy and molecular modeling was employed to gain information about the solution structure. Comparison of the crystal structure with its solution counterpart showed no significant differences in conformation, confirming the relative rigidity of the tricyclic system. However, these experiments revealed that the formerly proposed stereochemistry of branimycin at 17-C should be revised.

Virtual screening workflow development guided by the "receiver operating characteristic" curve approach. Application to high-throughput docking on metabotropic glutamate receptor subtype 4.

The "receiver operating characteristic" (ROC) curve method is a well-recognized metric used as an objective way to evaluate the ability of a given test to discriminate between two populations. This facilitates decision-making in a plethora of fields in which a wrong judgment may have serious consequences including clinical diagnosis, public safety, travel security, and economic strategies. When virtual screening is used to speed-up the drug discovery process in pharmaceutical research, taking the right decision upon selecting or discarding a molecule prior to in vitro evaluation is of paramount importance. Characterizing both the ability of a virtual screening workflow to select active molecules and the ability to discard inactive ones, the ROC curve approach is well suited for this critical decision gate. As a case study, the first virtual screening workflow focused on metabotropic glutamate receptor subtype 4 (mGlu4R) agonists is reported here. Six compounds out of 38 selected and tested in vitro were shown to have agonist activity on this target of therapeutic interest.

Discovery and optimization of an azetidine chemical series as a free fatty acid receptor 2 (FFA2) antagonist: from hit to clinic.

FFA2, also called GPR43, is a G-protein coupled receptor for short chain fatty acids which is involved in the mediation of inflammatory responses. A class of azetidines was developed as potent FFA2 antagonists. Multiparametric optimization of early hits with moderate potency and suboptimal ADME properties led to the identification of several compounds with nanomolar potency on the receptor combined with excellent pharmacokinetic (PK) parameters. The most advanced compound, 4-[[(R)-1-(benzo[b]thiophene-3-carbonyl)-2-methyl-azetidine-2-carbonyl]-(3-chloro-benzyl)-amino]-butyric acid 99 (GLPG0974), is able to inhibit acetate-induced neutrophil migration strongly in vitro and demonstrated ability to inhibit a neutrophil-based pharmacodynamic (PD) marker, CD11b activation-specific epitope [AE], in a human whole blood assay. All together, these data supported the progression of 99 toward next phases, becoming the first FFA2 antagonist to reach the clinic.

Discovery of diarylhydantoins as new selective androgen receptor modulators.

A novel selective androgen receptor modulator scaffold has been discovered through structural modifications of hydantoin antiandrogens. Several 4-(4-hydroxyphenyl)-N-arylhydantoins displayed partial agonism with nanomolar in vitro potency in transactivation experiments using androgen receptor (AR) transfected cells. In a standard castrated male rat model, several compounds showed good anabolic activity on levator ani muscle, dissociated from the androgenic activity on ventral prostate, after oral dosing at 30 mg/kg. (+)-4-[3,4-Dimethyl-2,5-dioxo-4-(4-hydroxyphenyl)imidazolidin-1-yl]-2-(trifluoromethyl)benzonitrile ((+)-11b) displayed anabolic potency with a strong dissociation between levator ani muscle and ventral prostate (A(50) = 0.5 mg/kg vs 70 mg/kg). The binding modes of two compounds, including (+)-11b, within the AR ligand-binding domain have been studied by cocrystallization experiments using a coactivator-like peptide. Both compounds bound to the same site, and the overall structures of the AR were very similar.

Identification of a 4-(hydroxymethyl)diarylhydantoin as a selective androgen receptor modulator.

Structural modification performed on a 4-methyl-4-(4-hydroxyphenyl)hydantoin series is described which resulted in the development of a new series of 4-(hydroxymethyl)diarylhydantoin analogues as potent, partial agonists of the human androgen receptor. This led to the identification of (S)-(-)-4-(4-(hydroxymethyl)-3-methyl-2,5-dioxo-4-phenylimidazolidin-1-yl)-2-(trifluoromethyl)benzonitrile ((S)-(-)-18a, GLPG0492) evaluated in vivo in a classical model of orchidectomized rat. In this model, (-)-18a exhibited anabolic activity on muscle, strongly dissociated from the androgenic activity on prostate after oral dosing. (-)-18a has very good pharmacokinetic properties, including bioavailability in rat (F > 50%), and is currently under evaluation in phase I clinical trials.

A virtual screening hit reveals new possibilities for developing group III metabotropic glutamate receptor agonists.

(R)-PCEP (3-amino-3-carboxypropyl-2'-carboxyethyl phosphinic acid, 1), a new metabotropic glutamate receptor 4 (mGlu4R) agonist, was discovered in a previously reported virtual screening. The (S)-enantiomer and a series of derivatives were synthesized and tested on recombinant mGlu4 receptors. A large number of derivatives activated this receptor but was not able to discriminate between mGlu4 and mGlu8 receptors. The most potent ones 6 and 12 displayed an EC(50) of 1.0 +/- 0.2 microM at mGlu4R. Interestingly these agonists with longer alkyl chains revealed a new binding pocket adjacent to the glutamate binding site, which is lined with residues that differ among the mGluR subtypes and that will allow the design of more selective compounds. Additionally 6 was able to activate mGlu7 receptor with an EC(50) of 43 +/- 16 microM and is thus significantly more potent than L-AP4 (EC(50) of 249 +/- 106 microM).

High-potency olfactory receptor agonists discovered by virtual high-throughput screening: molecular probes for receptor structure and olfactory function.

The detection of diverse chemical structures by the vertebrate olfactory system is accomplished by the recognition of odorous ligands by their cognate receptors. In the present study, we used computational screening to discover novel high-affinity agonists of an olfactory G protein-coupled receptor that recognizes amino acid ligands. Functional testing of the top candidates validated several agonists with potencies higher than any of the receptor's known natural ligands. Computational modeling revealed molecular interactions involved in ligand binding and further highlighted interactions that have been conserved in evolutionarily divergent amino acid receptors. Significantly, the top compounds display robust activities as odorants in vivo and include a natural product that may be used to signal the presence of bacteria in the environment. Our virtual screening approach should be applicable to the identification of new bioactive molecules for probing the structure of chemosensory receptors and the function of chemosensory systems in vivo.

Amino-pyrrolidine tricarboxylic acids give new insight into group III metabotropic glutamate receptor activation mechanism.

Like most class C G-protein-coupled receptors, metabotropic glutamate (mGlu) receptors possess a large extracellular domain where orthosteric ligands bind. Crystal structures revealed that this domain, called Venus FlyTrap (VFT), adopts a closed or open conformation upon agonist or antagonist binding, respectively. We have described amino-pyrrolidine tricarboxylic acids (APTCs), including (2S,4S)-4-amino-1-[(E)-3-carboxyacryloyl]pyrrolidine-2,4-dicarboxylic acid (FP0429), as new selective group III mGlu agonists. Whereas FP0429 is an almost full mGlu4 agonist, it is a weak and partial agonist of the closely related mGlu8 subtype. To get more insight into the activation mechanism of mGlu receptors, we aimed to elucidate why FP0429 behaves differently at these two highly homologous receptors by focusing on two residues within the binding site that differ between mGlu4 and mGlu8. Site-directed mutagenesis of Ser157 and Gly158 of mGlu4 into their mGlu8 homologs (Ala) turned FP0429 into a weak partial agonist. Conversely, introduction of Ser and Gly residues into mGlu8 increased FP0429 efficacy. Docking of FP0429 in mGlu4 VFT 3D model helped us characterize the role of each residue. Indeed, mGlu4 Ser157 seems to have an important role in FP0429 binding, whereas Gly158 may allow a deeper positioning of this agonist in the cavity of lobe I, thereby ensuring optimal interactions with lobe II residues in the fully closed state of the VFT. In contrast, the presence of a methyl group in mGlu8 (Ala instead of Gly) weakens the interactions with the lobe II residues. This probably results in a less stable or a partially closed form of the mGlu8 VFT, leading to partial receptor activation.

Design and synthesis of APTCs (aminopyrrolidinetricarboxylic acids): identification of a new group III metabotropic glutamate receptor selective agonist.

A new family of mGlu receptor orthosteric ligands called APTCs was designed and synthesized using a parallel chemistry approach. Amongst 65 molecules tested on mGlu4, mGlu6 and mGlu8 subtypes, (2S,4S)-4-amino-1-[(E)-3-carboxyacryloyl]pyrrolidine-2,4-dicarboxylic acid (8a06-FP0429) has been shown to be a full mGlu4 agonist and a partial mGlu8 agonist. In addition, 8a06 was shown to be selective versus group I and II mGlu subtypes. A possible explanation for this efficacy difference is proposed by docking experiment performed with molecular model of the receptor.

Closure of the Venus flytrap module of mGlu8 receptor and the activation process: Insights from mutations converting antagonists into agonists.

Ca2+, pheromones, sweet taste compounds, and the main neurotransmitters glutamate and gamma-aminobutyric acid activate G protein-coupled receptors (GPCRs) that constitute the GPCR family 3. These receptors are dimers, and each subunit has a large extracellular domain called a Venus flytrap module (VFTM), where agonists bind. This module is connected to a heptahelical domain that activates G proteins. Recently, the structure of the dimer of mGlu1 VFTMs revealed two important conformational changes resulting from glutamate binding. First, agonists can stabilize a closed state of at least one VFTM in the dimer. Second, the relative orientation of the two VFTMs in the dimer is different in the presence of glutamate, such that their C-terminal ends (which are connected to the G protein-activating heptahelical domain) become closer by more than 20 A. This latter change in orientation has been proposed to play a key role in receptor activation. To elucidate the respective role of VFTM closure and the change in orientation of the VFTMs in family 3 GPCR activation, we analyzed the mechanism of action of the mGlu8 receptor antagonists ACPT-II and MAP4. Molecular modeling studies suggest that these two compounds prevent the closure of the mGlu8 VFTM because of ionic and steric hindrance, respectively. We show here that the replacement of the residues responsible for these hindrances (Asp-309 and Tyr-227, respectively) by Ala allows ACPT-II or MAP4 to fully activate the receptors. These data are consistent with the requirement of the VFTM closure for family 3 GPCR activation.