Smallest Maximum Intramolecular Distance: A Novel Method to Mitigate Pregnane Xenobiotic Receptor Activation.

Induction of cytochrome P450 isoform 3A4 via activation of the pregnane xenobiotic receptor (PXR) is a concern for pharmaceutical discovery and development, as it can lead to drug-drug interactions. We present a novel molecular descriptor, the smallest maximum intramolecular distance (SMID), which is correlated with PXR activation, and a method for using the SMID descriptor to guide discovery chemists in modifying lead compounds to decrease PXR activation.

D3R grand challenge 4: blind prediction of protein-ligand poses, affinity rankings, and relative binding free energies.

The Drug Design Data Resource (D3R) aims to identify best practice methods for computer aided drug design through blinded ligand pose prediction and affinity challenges. Herein, we report on the results of Grand Challenge 4 (GC4). GC4 focused on proteins beta secretase 1 and Cathepsin S, and was run in an analogous manner to prior challenges. In Stage 1, participant ability to predict the pose and affinity of BACE1 ligands were assessed. Following the completion of Stage 1, all BACE1 co-crystal structures were released, and Stage 2 tested affinity rankings with co-crystal structures. We provide an analysis of the results and discuss insights into determined best practice methods.

A Brief Recent History of Women in Computational Chemistry at Vertex Pharmaceuticals.

Countless reports cite the importance of diversity in the academic, industrial, and government workplace. This article shares the different perspective on gender diversity from five women who have recently joined Vertex's computational chemistry group. It is written with the hope that other scientists will take the themes which resonant and adopt them to their own institutions to inspire the fostering of an inclusive environment while in pursuit of scientific discoveries.

Indole acids as a novel PDE2 inhibitor chemotype that demonstrate pro-cognitive activity in multiple species.

An internal HTS effort identified a novel PDE2 inhibitor series that was subsequently optimized for improved PDE2 activity and off-target selectivity. The optimized lead, compound 4, improved cognitive performance in a rodent novel object recognition task as well as a non-human primate object retrieval task. In addition, co-crystallization studies of close analog of 4 in the PDE2 active site revealed unique binding interactions influencing the high PDE isoform selectivity.

Computational Chemistry: A Rising Tide of Women.

The authors were inspired to explore the topic of gender diversity in computational chemistry on the basis of similar recent publications in the related fields of medicinal chemistry ( Huryn , D. M. ; et al. ACS Med. Chem. Lett. 2017 , 8 , 900 ) and computational biology ( Bonham , K. S. ; Stefan , M. I. PLoS Comput. Biol. 2017 , 13 , e1005134 ). To do so, we examined historical demographics in two different professional settings, i.e., attendance/participation at the Gordon Research Conferences on Computer-Aided Drug Design and Computational Chemistry and membership in the Computers in Chemistry Division of the American Chemical Society. We conclude that female representation in computational chemistry has risen steadily over the last 40 years and likely stands at around 25%, which appears to slightly exceed that of the neighboring fields of computer science and medicinal chemistry. In accordance with the old slogan that "a rising tide lifts all boats", here a rising tide of women scientists is having an impact on the field of computational chemistry. Tactics to ensure that this number continues to improve are highlighted.

Modeling & Informatics at Vertex Pharmaceuticals Incorporated: our philosophy for sustained impact.

Molecular modelers and informaticians have the unique opportunity to integrate cross-functional data using a myriad of tools, methods and visuals to generate information. Using their drug discovery expertise, information is transformed to knowledge that impacts drug discovery. These insights are often times formulated locally and then applied more broadly, which influence the discovery of new medicines. This is particularly true in an organization where the members are exposed to projects throughout an organization, such as in the case of the global Modeling & Informatics group at Vertex Pharmaceuticals. From its inception, Vertex has been a leader in the development and use of computational methods for drug discovery. In this paper, we describe the Modeling & Informatics group at Vertex and the underlying philosophy, which has driven this team to sustain impact on the discovery of first-in-class transformative medicines.

Discovery of pyrazolopyrimidine phosphodiesterase 10A inhibitors for the treatment of schizophrenia.

Herein, we present the identification of a novel class of pyrazolopyrimidine phosphodiesterase 10A (PDE10A) inhibitors. Beginning with a lead molecule (1) identified through a fragment-based drug discovery (FBDD) effort, lead optimization was enabled by rational design, X-ray crystallography, metabolic and off-target profiling, and fragment scaffold-hopping. We highlight the discovery of PyP-1, a potent, highly selective, and orally bioavailable pyrazolopyrimidine inhibitor of PDE10A. PyP-1 exhibits sub-nanomolar potency (PDE10A Ki=0.23nM), excellent pharmacokinetic (PK) and physicochemical properties, and a clean off-target profile. It displays dose-dependent efficacy in numerous pharmacodynamic (PD) assays that measure potential for anti-psychotic activity and cognitive improvement. PyP-1 also has a clean preclinical profile with respect to cataleptic potential in rats, prolactin secretion, and weight gain, common adverse events associated with currently marketed therapeutics. Further, PyP-1 displays in vivo preclinical target engagement as measured by PET enzyme occupancy in concert with [(11)C]MK-8193, a novel PDE10A PET tracer.

Prediction of Protein Pairs Sharing Common Active Ligands Using Protein Sequence, Structure, and Ligand Similarity.

We benchmarked the ability of comparative computational approaches to correctly discriminate protein pairs sharing a common active ligand (positive protein pairs) from protein pairs with no common active ligands (negative protein pairs). Since the target and the off-targets of a drug share at least a common ligand, i.e., the drug itself, the prediction of positive protein pairs may help identify off-targets. We evaluated representative protein-centric and ligand-centric approaches, including (1) 2D and 3D ligand similarity, (2) several measures of protein sequence similarity in conjunction with different sequence sources (e.g., full protein sequence versus binding site residues), and (3) a newly described pocket shape similarity and alignment program called SiteHopper. While the sequence-based alignment of pocket residues achieved the best overall performance, SiteHopper outperformed sequence-based approaches for unrelated proteins with only 20-30% pocket residue identity. Analogously, among ligand-centric approaches, path-based fingerprints achieved the best overall performance, but ROCS-based ligand shape similarity outperformed path-based fingerprints for structurally dissimilar ligands (Tanimoto 25%-40%). A significant drop in recognition performance was observed for ligand-centric approaches when PDB ligands were used instead of ChEMBL ligands. Finally, we analyzed the relationship between pocket shape and ligand shape in our data set and found that similar ligands tend to bind to similar pockets while similar pockets may accept a range of different-shaped ligands.

Collaborating to improve the use of free-energy and other quantitative methods in drug discovery.

In May and August, 2016, several pharmaceutical companies convened to discuss and compare experiences with Free Energy Perturbation (FEP). This unusual synchronization of interest was prompted by Schrödinger's FEP+ implementation and offered the opportunity to share fresh studies with FEP and enable broader discussions on the topic. This article summarizes key conclusions of the meetings, including a path forward of actions for this group to aid the accelerated evaluation, application and development of free energy and related quantitative, structure-based design methods.

Discovery of [¹¹C]MK-8193 as a PET tracer to measure target engagement of phosphodiesterase 10A (PDE10A) inhibitors.

Phosphodiesterase 10A (PDE10A) inhibition has recently been identified as a potential mechanism to treat multiple symptoms that manifest in schizophrenia. In order to facilitate preclinical development and support key proof-of-concept clinical trials of novel PDE10A inhibitors, it is critical to discover positron emission tomography (PET) tracers that enable plasma concentration/PDE10A occupancy relationships to be established across species with structurally diverse PDE10A inhibitors. In this Letter, we describe how a high-throughput screening hit was optimized to provide [(11)C]MK-8193 (8j), a PET tracer that supports the determination of plasma concentration/PDE10A occupancy relationships for structurally diverse series of PDE10A inhibitors in both rat and rhesus monkey.

Methyl-substitution of an iminohydantoin spiropiperidine ß-secretase (BACE-1) inhibitor has a profound effect on its potency.

The IC50 of a beta-secretase (BACE-1) lead compound was improved ∼200-fold from 11 µM to 55 nM through the addition of a single methyl group. Computational chemistry, small molecule NMR, and protein crystallography capabilities were used to compare the solution conformation of the ligand under varying pH conditions to its conformation when bound in the active site. Chemical modification then explored available binding pockets adjacent to the ligand. A strategically placed methyl group not only maintained the required pKa of the piperidine nitrogen and filled a small hydrophobic pocket, but more importantly, stabilized the conformation best suited for optimized binding to the receptor.

Shaping suvorexant: application of experimental and theoretical methods for driving synthetic designs.

Dual Orexin Receptor Antagonists (DORA) bind to both the Orexin 1 and 2 receptors. High resolution crystal structures of the Orexin 1 and 2 receptors, both class A GPCRs, were not available at the time of this study, and thus, ligand-based analyses were invoked and successfully applied to the design of DORAs. Computational analysis, ligand based superposition, unbound small-molecule X-ray crystal structures and NMR analysis were utilized to understand the conformational preferences of key DORAs and excellent agreement between these orthogonal approaches was seen in the majority of compounds examined. The predominantly face-to-face (F2F) interaction observed between the distal aromatic rings was the core 3D shape motif in our design principle and was used in the development of compounds. A notable exception, however, was seen between computation and experiment for suvorexant where the molecule exhibits an extended conformation in the unbound small-molecule X-ray structure. Even taking into account solvation effects explicitly in our calculations, we nevertheless find support that the F2F conformation is the bioactive conformation. Using a dominant states approximation for the partition function, we made a comprehensive assessment of the free energies required to adopt both an extended and a F2F conformation of a number of DORAs. Interestingly, we find that only a F2F conformation is consistent with the activities reported.

Discovery of MK-8189, a Highly Potent and Selective PDE10A Inhibitor for the Treatment of Schizophrenia.

PDE10A is an important regulator of striatal signaling that, when inhibited, can normalize dysfunctional activity. Given the involvement of dysfunctional striatal activity with schizophrenia, PDE10A inhibition represents a potentially novel means for its treatment. With the goal of developing PDE10A inhibitors, early optimization of a fragment hit through rational design led to a series of potent pyrimidine PDE10A inhibitors that required further improvements in physicochemical properties, off-target activities, and pharmacokinetics. Herein we describe the discovery of an isomeric pyrimidine series that addresses the liabilities seen with earlier compounds and resulted in the invention of compound 18 (MK-8189), which is currently in Phase 2b clinical development for the treatment of schizophrenia.

Antiviral activity and in vitro mutation development pathways of MK-6186, a novel nonnucleoside reverse transcriptase inhibitor.

MK-6186 is a novel nonnucleoside reverse transcriptase inhibitor (NNRTI) which displays subnanomolar potency against wild-type (WT) virus and the two most prevalent NNRTI-resistant RT mutants (K103N and Y181C) in biochemical assays. In addition, it showed excellent antiviral potency against K103N and Y181C mutant viruses, with fold changes (FCs) of less than 2 and 5, respectively. When a panel of 12 common NNRTI-associated mutant viruses was tested with MK-6186, only 2 relatively rare mutants (Y188L and V106I/Y188L) were highly resistant, with FCs of >100, and the remaining viruses showed FCs of <10. Furthermore, a panel of 96 clinical virus isolates with NNRTI resistance mutations was evaluated for susceptibility to NNRTIs. The majority (70%) of viruses tested displayed resistance to efavirenz (EFV), with FCs of >10, whereas only 29% of the mutant viruses displayed greater than 10-fold resistance to MK-6186. To determine whether MK-6186 selects for novel resistance mutations, in vitro resistance selections were conducted with one isolate each from subtypes A, B, and C under low-multiplicity-of-infection (MOI) conditions. The results showed a unique mutation development pattern in which L234I was the first mutation to emerge in the majority of the experiments. In resistance selection under high-MOI conditions with subtype B virus, V106A was the dominant mutation detected in the breakthrough viruses. More importantly, mutant viruses selected by MK-6186 showed FCs of <10 against EFV or etravirine (ETR), and the mutant viruses containing mutations selected by EFV or ETR were sensitive to MK-6186 (FCs of <10).

Discovery of pyrrolidine-based ß-secretase inhibitors: lead advancement through conformational design for maintenance of ligand binding efficiency.

We have developed a novel series of pyrrolidine derived BACE-1 inhibitors. The potency of the weak initial lead structure was enhanced using library-based SAR methods. The series was then further advanced by rational design while maintaining a minimal ligand binding efficiency threshold. Ultimately, the co-crystal structure was obtained revealing that these inhibitors interacted with the enzyme in a unique fashion. In all, the potency of the series was enhanced by 4 orders of magnitude from the HTS lead with concomitant increases in physical properties needed for series advancement. The progression of these developments in a systematic fashion is described.

Discovery of tetrahydropyridopyrimidine phosphodiesterase 10A inhibitors for the treatment of schizophrenia.

We describe the discovery of potent and orally bioavailable tetrahydropyridopyrimidine inhibitors of phosphodiesterase 10A by systematic optimization of a novel HTS lead. Lead compound THPP-1 exhibits nanomolar potencies, excellent pharmacokinetic properties, and a clean off-target profile. It displays in vivo target engagement as measured by increased rat striatal cGMP levels upon oral dosing. It shows dose-dependent efficacy in a key pharmacodynamic assay predictive of antipsychotic activity, the psychostimulant-induced rat hyperlocomotion assay. Further, THPP-1 displays significantly fewer preclinical adverse events in assays measuring prolactin secretion, catalepsy, and weight gain, in contrast to the typical and atypical antipsychotics haloperidol and olanzapine.

Small molecule mimetics of an HIV-1 gp41 fusion intermediate as vaccine leads.

We describe here a novel platform technology for the discovery of small molecule mimetics of conformational epitopes on protein antigens. As a model system, we selected mimetics of a conserved hydrophobic pocket within the N-heptad repeat region of the HIV-1 envelope protein, gp41. The human monoclonal antibody, D5, binds to this target and exhibits broadly neutralizing activity against HIV-1. We exploited the antigen-binding property of D5 to select complementary small molecules using a high throughput screen of a diverse chemical collection. The resulting small molecule leads were rendered immunogenic by linking them to a carrier protein and were shown to elicit N-heptad repeat-binding antibodies in a fraction of immunized mice. Plasma from HIV-1-infected subjects shown previously to contain broadly neutralizing antibodies was found to contain antibodies capable of binding to haptens represented in the benzylpiperidine leads identified as a result of the high throughput screen, further validating these molecules as vaccine leads. Our results suggest a new paradigm for vaccine discovery using a medicinal chemistry approach to identify lead molecules that, when optimized, could become vaccine candidates for infectious diseases that have been refractory to conventional vaccine development.

Design and synthesis of pyridone inhibitors of non-nucleoside reverse transcriptase.

Next generation NNRTIs are sought which possess both broad spectrum antiviral activity against key mutant strains and a high genetic barrier to the selection of new mutant viral strains. Pyridones were evaluated as an acyclic conformational constraint to replace the aryl ether core of MK-4965 (1) and the more rigid indazole constraint of MK-6186 (2). The resulting pyridone compounds are potent inhibitors of HIV RT and have antiviral activity in cell culture that is superior to other next generation NNRTI's.

Pyridyl amides as potent inhibitors of T-type calcium channels.

A novel series of amide T-type calcium channel antagonists were prepared and evaluated using in vitro and in vivo assays. Optimization of the screening hit 3 led to identification of the potent and selective T-type antagonist 37 that displayed in vivo efficacy in rodent models of epilepsy and sleep.

Distinct mutation pathways of non-subtype B HIV-1 during in vitro resistance selection with nonnucleoside reverse transcriptase inhibitors.

Studies were conducted to investigate mutation pathways among subtypes A, B, and C of human immunodeficiency virus type 1 (HIV-1) during resistance selection with nonnucleoside reverse transcriptase inhibitors (NNRTIs) in cell culture under low-multiplicity of infection (MOI) conditions. The results showed that distinct pathways were selected by different virus subtypes under increasing selective pressure of NNRTIs. F227C and Y181C were the major mutations selected by MK-4965 in subtype A and C viruses during resistance selection. With efavirenz (EFV), F227C and V106M were the major mutations responsible for viral breakthrough in subtype A viruses, whereas a single pathway (G190A/V106M) accounted for mutation development in subtype C viruses. Y181C was the dominant mutation in the resistance selection with etravirine (ETV) in subtype A, and E138K/H221Y were the mutations detected in the breakthrough viruses from subtype C viruses with ETV. In subtype B viruses, on the other hand, known NNRTI-associated mutations (e.g., Y181C, P236L, L100I, V179D, and K103N) were selected by the NNRTIs. The susceptibility of the subtype A and B mutant viruses to NNRTIs was determined in order to gain insight into the potential mechanisms of mutation development. Collectively, these results suggest that minor differences may exist in conformation of the residues within the NNRTI binding pocket (NNRTIBP) of reverse transcriptase (RT) among the three subtypes of viruses. Thus, the interactions between NNRTIs and the residues in the NNRTIBPs of different subtypes may not be identical, leading to distinct mutation pathways during resistance selection in cell culture.