Selective inhibitors of the PSEN1-gamma-secretase complex.

Clinical development of γ-secretases, a family of intramembrane cleaving proteases, as therapeutic targets for a variety of disorders including cancer and Alzheimer's disease was aborted because of serious mechanism-based side effects in the phase III trials of unselective inhibitors. Selective inhibition of specific γ-secretase complexes, containing either PSEN1 or PSEN2 as the catalytic subunit and APH1A or APH1B as supporting subunits, does provide a feasible therapeutic window in preclinical models of these disorders. We explore here the pharmacophoric features required for PSEN1 versus PSEN2 selective inhibition. We synthesized a series of brain penetrant 2-azabicyclo[2,2,2]octane sulfonamides and identified a compound with low nanomolar potency and high selectivity (>250-fold) toward the PSEN1-APH1B subcomplex versus PSEN2 subcomplexes. We used modeling and site-directed mutagenesis to identify critical amino acids along the entry part of this inhibitor into the catalytic site of PSEN1. Specific targeting one of the different γ-secretase complexes might provide safer drugs in the future.

Current State of Open Source Force Fields in Protein-Ligand Binding Affinity Predictions.

In drug discovery, the in silico prediction of binding affinity is one of the major means to prioritize compounds for synthesis. Alchemical relative binding free energy (RBFE) calculations based on molecular dynamics (MD) simulations are nowadays a popular approach for the accurate affinity ranking of compounds. MD simulations rely on empirical force field parameters, which strongly influence the accuracy of the predicted affinities. Here, we evaluate the ability of six different small-molecule force fields to predict experimental protein-ligand binding affinities in RBFE calculations on a set of 598 ligands and 22 protein targets. The public force fields OpenFF Parsley and Sage, GAFF, and CGenFF show comparable accuracy, while OPLS3e is significantly more accurate. However, a consensus approach using Sage, GAFF, and CGenFF leads to accuracy comparable to OPLS3e. While Parsley and Sage are performing comparably based on aggregated statistics across the whole dataset, there are differences in terms of outliers. Analysis of the force field reveals that improved parameters lead to significant improvement in the accuracy of affinity predictions on subsets of the dataset involving those parameters. Lower accuracy can not only be attributed to the force field parameters but is also dependent on input preparation and sampling convergence of the calculations. Especially large perturbations and nonconverged simulations lead to less accurate predictions. The input structures, Gromacs force field files, as well as the analysis Python notebooks are available on GitHub.

To Design Scalable Free Energy Perturbation Networks, Optimal Is Not Enough.

Drug discovery is accelerated with computational methods such as alchemical simulations to estimate ligand affinities. In particular, relative binding free energy (RBFE) simulations are beneficial for lead optimization. To use RBFE simulations to compare prospective ligands in silico, researchers first plan the simulation experiment, using graphs where nodes represent ligands and graph edges represent alchemical transformations between ligands. Recent work demonstrated that optimizing the statistical architecture of these perturbation graphs improves the accuracy of the predicted changes in the free energy of ligand binding. Therefore, to improve the success rate of computational drug discovery, we present the open-source software package High Information Mapper (HiMap)─a new take on its predecessor, Lead Optimization Mapper (LOMAP). HiMap removes heuristics decisions from design selection and instead finds statistically optimal graphs over ligands clustered with machine learning. Beyond optimal design generation, we present theoretical insights for designing alchemical perturbation maps. Some of these results include that for n number of nodes, the precision of perturbation maps is stable at n·ln(n) edges. This result indicates that even an "optimal" graph can result in unexpectedly high errors if a plan includes too few alchemical transformations for the given number of ligands and edges. And, as a study compares more ligands, the performance of even optimal graphs will deteriorate with linear scaling of the edge count. In this sense, ensuring an A- or D-optimal topology is not enough to produce robust errors. We additionally find that optimal designs will converge more rapidly than radial and LOMAP designs. Moreover, we derive bounds for how clustering reduces cost for designs with a constant expected relative error per cluster, invariant of the size of the design. These results inform how to best design perturbation maps for computational drug discovery and have broader implications for experimental design.

MELD-Bracket Ranks Binding Affinities of Diverse Sets of Ligands.

Affinity ranking of structurally diverse small-molecule ligands is a challenging problem with important applications in structure-based drug discovery. Absolute binding free energy methods can model diverse ligands, but the high computational cost of the current methods limits application to data sets with few ligands. We recently developed MELD-Bracket, a Molecular Dynamics method for efficient affinity ranking of ligands [ JCTC 2022, 18 (1), 374-379]. It utilizes a Bayesian framework to guide sampling to relevant regions of phase space, and it couples this with a bracket-like competition on a pool of ligands. Here we find that 6-competitor MELD-Bracket can rank dozens of diverse ligands that have low structural similarity and different net charges. We benchmark it on four protein systems─PTB1B, Tyk2, BACE, and JAK3─having varied modes of interactions. We also validated 8-competitor and 12-competitor protocols. The MELD-Bracket protocols presented here may have the appropriate balance of accuracy and computational efficiency to be suitable for ranking diverse ligands from typical drug discovery campaigns.

The Impact of Data on Structure-Based Binding Affinity Predictions Using Deep Neural Networks.

Artificial intelligence (AI) has gained significant traction in the field of drug discovery, with deep learning (DL) algorithms playing a crucial role in predicting protein-ligand binding affinities. Despite advancements in neural network architectures, system representation, and training techniques, the performance of DL affinity prediction has reached a plateau, prompting the question of whether it is truly solved or if the current performance is overly optimistic and reliant on biased, easily predictable data. Like other DL-related problems, this issue seems to stem from the training and test sets used when building the models. In this work, we investigate the impact of several parameters related to the input data on the performance of neural network affinity prediction models. Notably, we identify the size of the binding pocket as a critical factor influencing the performance of our statistical models; furthermore, it is more important to train a model with as much data as possible than to restrict the training to only high-quality datasets. Finally, we also confirm the bias in the typically used current test sets. Therefore, several types of evaluation and benchmarking are required to understand models' decision-making processes and accurately compare the performance of models.

The Impact of Experimental and Calculated Error on the Performance of Affinity Predictions.

The accurate prediction of binding affinity between protein and small molecules with free energy methods, particularly the difference in binding affinities via relative binding free energy calculations, has undergone a dramatic increase in use and impact over recent years. The improvements in methodology, hardware, and implementation can deliver results with less than 1 kcal/mol mean unsigned error between calculation and experiment. This is a remarkable achievement and beckons some reflection on the significance of calculation approaching the accuracy of experiment. In this article, we describe a statistical analysis of the implications of variance (standard deviation) of both experimental and calculated binding affinities with respect to the unknown true binding affinity. We reveal that plausible ratios of standard deviation in experiment and calculation can lead to unexpected outcomes for assessing the performance of predictions. The work extends beyond the case of binding free energies to other affinity or property prediction methods.

Divide and Conquer. Pocket-Opening Mixed-Solvent Simulations in the Perspective of Docking Virtual Screening Applications for Drug Discovery.

The existence of a druggable binding pocket is a prerequisite for computational drug-target interaction studies including virtual screening. Retrospective studies have shown that extended sampling methods like Markov State Modeling and mixed-solvent simulations can identify cryptic pockets relevant for drug discovery. Here, we apply a combination of mixed-solvent molecular dynamics (MD) and time-structure independent component analysis (TICA) to four retrospective case studies: NPC2, the CECR2 bromodomain, TEM-1, and MCL-1. We compare previous experimental and computational findings to our results. It is shown that the successful identification of cryptic pockets depends on the system and the cosolvent probes. We used alternative TICA internal features such as the unbiased backbone coordinates or backbone dihedrals versus biased interatomic distances. We found that in the case of NPC2, TEM-1, and MCL-1, the use of unbiased features is able to identify cryptic pockets, although in the case of the CECR2 bromodomain, more specific features are required to properly capture a pocket opening. In the perspective of virtual screening applications, it is shown how docking studies with the parent ligands depend critically on the conformational state of the targets.

Pre-Exascale Computing of Protein-Ligand Binding Free Energies with Open Source Software for Drug Design.

Nowadays, drug design projects benefit from highly accurate protein-ligand binding free energy predictions based on molecular dynamics simulations. While such calculations have been computationally expensive in the past, we now demonstrate that workflows built on open source software packages can efficiently leverage pre-exascale computing resources to screen hundreds of compounds in a matter of days. We report our results of free energy calculations on a large set of pharmaceutically relevant targets assembled to reflect industrial drug discovery projects.

Collaborative Assessment of Molecular Geometries and Energies from the Open Force Field.

Force fields form the basis for classical molecular simulations, and their accuracy is crucial for the quality of, for instance, protein-ligand binding simulations in drug discovery. The huge diversity of small-molecule chemistry makes it a challenge to build and parameterize a suitable force field. The Open Force Field Initiative is a combined industry and academic consortium developing a state-of-the-art small-molecule force field. In this report, industry members of the consortium worked together to objectively evaluate the performance of the force fields (referred to here as OpenFF) produced by the initiative on a combined public and proprietary dataset of 19,653 relevant molecules selected from their internal research and compound collections. This evaluation was important because it was completely blind; at most partners, none of the molecules or data were used in force field development or testing prior to this work. We compare the Open Force Field "Sage" version 2.0.0 and "Parsley" version 1.3.0 with GAFF-2.11-AM1BCC, OPLS4, and SMIRNOFF99Frosst. We analyzed force-field-optimized geometries and conformer energies compared to reference quantum mechanical data. We show that OPLS4 performs best, and the latest Open Force Field release shows a clear improvement compared to its predecessors. The performance of established force fields such as GAFF-2.11 was generally worse. While OpenFF researchers were involved in building the benchmarking infrastructure used in this work, benchmarking was done entirely in-house within industrial organizations and the resulting assessment is reported here. This work assesses the force field performance using separate benchmarking steps, external datasets, and involving external research groups. This effort may also be unique in terms of the number of different industrial partners involved, with 10 different companies participating in the benchmark efforts.

Non-equilibrium approach for binding free energies in cyclodextrins in SAMPL7: force fields and software.

In the current work we report on our participation in the SAMPL7 challenge calculating absolute free energies of the host-guest systems, where 2 guest molecules were probed against 9 hosts-cyclodextrin and its derivatives. Our submission was based on the non-equilibrium free energy calculation protocol utilizing an averaged consensus result from two force fields (GAFF and CGenFF). The submitted prediction achieved accuracy of [Formula: see text] in terms of the unsigned error averaged over the whole dataset. Subsequently, we further report on the underlying reasons for discrepancies between our calculations and another submission to the SAMPL7 challenge which employed a similar methodology, but disparate ligand and water force fields. As a result we have uncovered a number of issues in the dihedral parameter definition of the GAFF 2 force field. In addition, we identified particular cases in the molecular topologies where different software packages had a different interpretation of the same force field. This latter observation might be of particular relevance for systematic comparisons of molecular simulation software packages. The aforementioned factors have an influence on the final free energy estimates and need to be considered when performing alchemical calculations.

Accurate Prediction of GPCR Ligand Binding Affinity with Free Energy Perturbation.

The computational prediction of relative binding free energies is a crucial goal for drug discovery, and G protein-coupled receptors (GPCRs) are arguably the most important drug target class. However, they present increased complexity to model compared to soluble globular proteins. Despite breakthroughs, experimental X-ray crystal and cryo-EM structures are challenging to attain, meaning computational models of the receptor and ligand binding mode are sometimes necessary. This leads to uncertainty in understanding ligand-protein binding induced changes such as, water positioning and displacement, side chain positioning, hydrogen bond networks, and the overall structure of the hydration shell around the ligand and protein. In other words, the very elements that define structure activity relationships (SARs) and are crucial for accurate binding free energy calculations are typically more uncertain for GPCRs. In this work we use free energy perturbation (FEP) to predict the relative binding free energies for ligands of two different GPCRs. We pinpoint the key aspects for success such as the important role of key water molecules, amino acid ionization states, and the benefit of equilibration with specific ligands. Initial calculations following typical FEP setup and execution protocols delivered no correlation with experiment, but we show how results are improved in a logical and systematic way. This approach gave, in the best cases, a coefficient of determination (R2) compared with experiment in the range of 0.6-0.9 and mean unsigned errors compared to experiment of 0.6-0.7 kcal/mol. We anticipate that our findings will be applicable to other difficult-to-model protein ligand data sets and be of wide interest for the community to continue improving FE binding energy predictions.

Selective inhibition of intestinal guanosine 3',5'-cyclic monophosphate signaling by small-molecule protein kinase inhibitors.

The guanosine 3',5'-cyclic monophosphate (cGMP)-dependent protein kinase II (cGKII) serine/threonine kinase relays signaling through guanylyl cyclase C (GCC) to control intestinal fluid homeostasis. Here, we report the discovery of small-molecule inhibitors of cGKII. These inhibitors were imidazole-aminopyrimidines, which blocked recombinant human cGKII at submicromolar concentrations but exhibited comparatively little activity toward the phylogenetically related protein kinases cGKI and cAMP-dependent protein kinase (PKA). Whereas aminopyrimidyl motifs are common in protein kinase inhibitors, molecular modeling of these imidazole-aminopyrimidines in the ATP-binding pocket of cGKII indicated an unconventional binding mode that directs their amine substituent into a narrow pocket delineated by hydrophobic residues of the hinge and the αC-helix. Crucially, this set of residues included the Leu-530 gatekeeper, which is not conserved in cGKI and PKA. In intestinal organoids, these compounds blocked cGKII-dependent phosphorylation of the vasodilator-stimulated phosphoprotein (VASP). In mouse small intestinal tissue, cGKII inhibition significantly attenuated the anion secretory response provoked by the GCC-activating bacterial heat-stable toxin (STa), a frequent cause of infectious secretory diarrhea. In contrast, both PKA-dependent VASP phosphorylation and intestinal anion secretion were unaffected by treatment with these compounds, whereas experiments with T84 cells indicated that they weakly inhibit the activity of cAMP-hydrolyzing phosphodiesterases. As these protein kinase inhibitors are the first to display selective inhibition of cGKII, they may expedite research on cGMP signaling and may aid future development of therapeutics for managing diarrheal disease and other pathogenic syndromes that involve cGKII.

The size matters? A computational tool to design bivalent ligands.

Motivation: Bivalent ligands are increasingly important such as for targeting G protein-coupled receptor (GPCR) dimers or proteolysis targeting chimeras (PROTACs). They contain two pharmacophoric units that simultaneously bind in their corresponding binding sites, connected with a spacer chain. Here, we report a molecular modelling tool that links the pharmacophore units via the shortest pathway along the receptors van der Waals surface and then scores the solutions providing prioritization for the design of new bivalent ligands. Results: Bivalent ligands of known dimers of GPCRs, PROTACs and a model bivalent antibody/antigen system were analysed. The tool could rapidly assess the preferred linker length for the different systems and recapitulated the best reported results. In the case of GPCR dimers the results suggest that in some cases these ligands might bind to a secondary binding site at the extracellular entrance (vestibule or allosteric site) instead of the orthosteric binding site. Availability and implementation: Freely accessible from the Molecular Operating Environment svl exchange server (https://svl.chemcomp.com/). Supplementary information: Supplementary data are available at Bioinformatics online.

Limitations of Ligand-Only Approaches for Predicting the Reactivity of Covalent Inhibitors.

Covalent inhibition has undergone a resurgence and is an important modern-day drug design and chemical biology approach. To avoid off-target interactions and to fine-tune reactivity, the ability to accurately predict reactivity is vitally important for the design and development of safer and more effective covalent drugs. Several ligand-only metrics have been proposed that promise quick and simple ways of determining covalent reactivity. In particular, we examine proton affinity and reaction energies calculated with the density functional B3LYP-D3/6-311+G**//B3LYP-D3/6-31G* method to assess the reactivity of a series of α,ß-unsaturated carbonyl compounds that form covalent adducts with cysteine. We demonstrate that while these metrics correlate well with experiment for a diverse range of small reactive molecules these approaches fail for predicting the reactivity of drug-like compounds. We conclude that ligand-only metrics such as proton affinity and reaction energies do not capture determinants of reactivity in situ and fail to account for important factors such as conformation, solvation, and intermolecular interactions.

Mechanisms Underlying Allosteric Molecular Switches of Metabotropic Glutamate Receptor 5.

The metabotropic glutamate 5 (mGlu5) receptor is a class C G protein-coupled receptor (GPCR) that is implicated in several CNS disorders making it a popular drug discovery target. Years of research have revealed allosteric mGlu5 ligands showing an unexpected complete switch in functional activity despite only small changes in their chemical structure, resulting in positive allosteric modulators (PAM) or negative allosteric modulators (NAM) for the same scaffold. Up to now, the origins of this effect are not understood, causing difficulties in a drug discovery context. In this work, experimental data was gathered and analyzed alongside docking and Molecular Dynamics (MD) calculations for three sets of PAM and NAM pairs. The results consistently show the role of specific interactions formed between ligand substituents and amino acid side chains that block or promote local movements associated with receptor activation. The work provides an explanation for how such small structural changes lead to remarkable differences in functional activity. While this work can greatly help drug discovery programs avoid these switches, it also provides valuable insight into the mechanisms of class C GPCR allosteric activation. Furthermore, the approach shows the value of applying MD to understand functional activity in drug design programs, even for such close structural analogues.

Computational Drug Design Applied to the Study of Metabotropic Glutamate Receptors.

Metabotropic glutamate (mGlu) receptors are a family of eight GPCRs that are attractive drug discovery targets to modulate glutamate action and response. Here we review the application of computational methods to the study of this family of receptors. X-ray structures of the extracellular and 7-transmembrane domains have played an important role to enable structure-based modeling approaches, whilst we also discuss the successful application of ligand-based methods. We summarize the literature and highlight the areas where modeling and experiment have delivered important understanding for mGlu receptor drug discovery. Finally, we offer suggestions of future areas of opportunity for computational work.

Molecular blueprint of allosteric binding sites in a homologue of the agonist-binding domain of the α7 nicotinic acetylcholine receptor.

The α7 nicotinic acetylcholine receptor (nAChR) belongs to the family of pentameric ligand-gated ion channels and is involved in fast synaptic signaling. In this study, we take advantage of a recently identified chimera of the extracellular domain of the native α7 nicotinic acetylcholine receptor and acetylcholine binding protein, termed α7-AChBP. This chimeric receptor was used to conduct an innovative fragment-library screening in combination with X-ray crystallography to identify allosteric binding sites. One allosteric site is surface-exposed and is located near the N-terminal α-helix of the extracellular domain. Ligand binding at this site causes a conformational change of the α-helix as the fragment wedges between the α-helix and a loop homologous to the main immunogenic region of the muscle α1 subunit. A second site is located in the vestibule of the receptor, in a preexisting intrasubunit pocket opposite the agonist binding site and corresponds to a previously identified site involved in positive allosteric modulation of the bacterial homolog ELIC. A third site is located at a pocket right below the agonist binding site. Using electrophysiological recordings on the human α7 nAChR we demonstrate that the identified fragments, which bind at these sites, can modulate receptor activation. This work presents a structural framework for different allosteric binding sites in the α7 nAChR and paves the way for future development of novel allosteric modulators with therapeutic potential.

Identification of Allosteric Modulators of Metabotropic Glutamate 7 Receptor Using Proteochemometric Modeling.

Proteochemometric modeling (PCM) is a computational approach that can be considered an extension of quantitative structure-activity relationship (QSAR) modeling, where a single model incorporates information for a family of targets and all the associated ligands instead of modeling activity versus one target. This is especially useful for situations where bioactivity data exists for similar proteins but is scarce for the protein of interest. Here we demonstrate the application of PCM to identify allosteric modulators of metabotropic glutamate (mGlu) receptors. Given our long-running interest in modulating mGlu receptor function we compiled a matrix of compound-target bioactivity data. Some members of the mGlu family are well explored both internally and in the public domain, while there are much fewer examples of ligands for other targets such as the mGlu7 receptor. Using a PCM approach mGlu7 receptor hits were found. In comparison to conventional single target modeling the identified hits were more diverse, had a better confirmation rate, and provide starting points for further exploration. We conclude that the robust structure-activity relationship from well explored target family members translated to better quality hits for PCM compared to virtual screening (VS) based on a single target.

Application of Free Energy Perturbation for the Design of BACE1 Inhibitors.

Novel spiroaminodihydropyrroles probing for optimized interactions at the P3 pocket of ß-secretase 1 (BACE1) were designed with the use of free energy perturbation (FEP) calculations. The resulting molecules showed pIC50 potencies in enzymatic BACE1 inhibition assays ranging from approximately 5 to 7. Good correlation was observed between the predicted activity from the FEP calculations and experimental activity. Simulations run with a default 5 ns approach delivered a mean unsigned error (MUE) between prediction and experiment of 0.58 and 0.91 kcal/mol for retrospective and prospective applications, respectively. With longer simulations of 10 and 20 ns, the MUE was in both cases 0.57 kcal/mol for the retrospective application, and 0.69 and 0.59 kcal/mol for the prospective application. Other considerations that impact the quality of the calculations are discussed. This work provides an example of the value of FEP as a computational tool for drug discovery.

A Versatile Approach to CF3 -Containing 2-Pyrrolidones by Tandem Michael Addition-Cyclization: Exemplification in the Synthesis of Amidine Class BACE1 Inhibitors.

The synthesis of new fluorinated pyrrolidones starting from unprotected amino esters and amino nitriles through a Michael addition-lactamization sequence is described. The resulting CF3 -containing building blocks, bearing a quaternary stereogenic center adjacent to the fluorinated group, have been converted into amino pyrrolidines that display potent ß-secretase 1 (BACE1) inhibitory activity. This work constitutes an example of selective fluorination as a valid strategy for the modulation of physicochemical and biological properties of lead compounds in drug discovery.