Computational Fluorine Scanning Using Free-Energy Perturbation

Scaffold Hopping and Optimization of Small Molecule Soluble Adenyl Cyclase Inhibitors Led by Free Energy Perturbation

Free energy perturbation is a computational technique that can be used to predict how small changes to an inhibitor structure will affect the binding free energy to its target. In this paper, we describe the utility of free energy perturbation with FEP+ in the hit-to-lead stage of a drug discovery project targeting soluble adenyl cyclase. The project was structurally enabled by X-ray crystallography throughout. We employed free energy perturbation to first scaffold hop to a preferable chemotype and then optimize the binding affinity to sub-nanomolar levels while retaining druglike properties. The results illustrate that effective use of free energy perturbation can enable a drug discovery campaign to progress rapidly from hit to lead, facilitating proof-of-concept studies that enable target validation.

Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design

Bioisosteres are molecules that differ in substituents but still have very similar shapes. Bioisosteric replacements are ubiquitous in modern drug design, where they are used to alter metabolism, change bioavailability, or modify activity of the lead compound. Prediction of relative affinities of bioisosteres with computational methods is a long-standing task; however, the very shape closeness makes bioisosteric substitutions almost intractable for computational methods, which use standard force fields. Here, we design a quantum mechanical (QM)-cluster approach based on the GFN2-xTB semi-empirical quantum-chemical method and apply it to a set of H → F bioisosteric replacements. The proposed methodology enables advanced prediction of biological activity change upon bioisosteric substitution of -H with -F, with the standard deviation of 0.60 kcal/mol, surpassing the ChemPLP scoring function (0.83 kcal/mol), and making QM-based ΔΔG estimation comparable to ∼0.42 kcal/mol standard deviation of in vitro experiment. The speed of the method and lack of tunable parameters makes it affordable in current drug research.

TIES 2.0: A Dual-Topology Open Source Relative Binding Free Energy Builder with Web Portal

Relative binding free energy (RBFE) calculations are widely used to aid the process of drug discovery. TIES, Thermodynamic Integration with Enhanced Sampling, is a dual-topology approach to RBFE calculations with support for NAMD and OpenMM molecular dynamics engines. The software has been thoroughly validated on publicly available datasets. Here we describe the open source software along with a web portal (https://ccs-ties.org) that enables users to perform such calculations correctly and rapidly.

Design, Synthesis, and Pharmacological Evaluation of Second-Generation Soluble Adenylyl Cyclase (sAC, ADCY10) Inhibitors with Slow Dissociation Rates

Soluble adenylyl cyclase (sAC: ADCY10) is an enzyme involved in intracellular signaling. Inhibition of sAC has potential therapeutic utility in a number of areas. For example, sAC is integral to successful male fertility: sAC activation is required for sperm motility and ability to undergo the acrosome reaction, two processes central to oocyte fertilization. Pharmacologic evaluation of existing sAC inhibitors for utility as on-demand, nonhormonal male contraceptives suggested that both high intrinsic potency, fast on and slow dissociation rates are essential design elements for successful male contraceptive applications. During the course of the medicinal chemistry campaign described here, we identified sAC inhibitors that fulfill these criteria and are suitable for in vivo evaluation of diverse sAC pharmacology.

In Silico Positional Analogue Scanning with Amber GPU-TI

Positional analogue scanning (PAS) is an accepted strategy for multiparameter lead optimization (MPO) in drug discovery. Small structural changes as introduced by PAS can lead to 10-fold changes in binding potency in ∼10-20% of cases, a significant parameter shift irrespective of other MPO objectives. Sometimes performing a complete PAS is challenging due to resource and time constraints, building block availability, or difficulty in synthesis. Calculating relative binding free energies (RBFEs) for all positions can contribute to prioritizing the most promising analogues for synthesis. We tested a well-established RBFE calculation method, Amber GPU-TI, for 20 positional analogue scans in 14 test systems (cyclin-dependent kinase 8 (CDK8), hepatitis C virus nonstructural protein 5B (HCV NS5B), tankyrase, RAC-α serine/threonine-protein kinase (Akt), phosphodiesterase 1B (PDE1B), orexin/hypocretin receptor type 1 (OX1R), orexin/hypocretin receptor type 2 (OX2R), histone acetyltransferase K (lysine) acetyltransferase 6A (KAT6A), peroxisome proliferator-activated receptor γ (PPARγ), extracellular signal-regulated kinases (ERK1/2), coactivator-associated arginine methyltransferase 1 (PRMT4), α_vβ₆, bromodomain 1 (BD1), human immunodeficiency virus-1 (HIV-1) entry) involving nitrogen, methyl, halogen, methoxy, and hydroxyl scans with at least four analogues per set. Among the 66 analogue positions explored, we found that in 18 cases Amber GPU-TI calculations predicted a more than 10-fold change in potency. In all of these cases, the experimentally observed direction of potency changes agreed with the predictions. In 16 cases, more than 10-fold changes in experimental potency were observed. Again, in all of these cases, Amber GPU-TI predicted the direction of the potency changes correctly. In none of these cases would a decision made for or against synthesis based on a 10-fold change in potency have resulted in missing an important analogue. Therefore, in silico RBFE calculations using Amber GPU-TI can meaningfully contribute to the prioritization of positional analogues before synthesis.

Comparing the Performance of Different AMBER Protein Forcefields, Partial Charge Assignments, and Water Models for Absolute Binding Free Energy Calculations

Identifying chemical starting points is a vital first step in small molecule drug discovery and can take significant time and money. For this reason, computational approaches to virtual screening are of great interest as they can lower the cost and shorten timeframes. However, simple approaches such as molecular docking and pharmacophore screening are of limited accuracy and provide a low probability of success. Alchemical binding free energies represent a promising approach for virtual screening as they naturally incorporate the key effects of water molecules, protein flexibility, and binding entropy. However, the calculations are technically very challenging, with performance depending on the specific forcefield used. For this reason, it is important that the community has access to benchmark test sets to assess prediction accuracy. In this paper, we present an approach to alchemical binding free energies using OpenMM. We identify effective simulation parameters using an existing BRD4(1) test set and present two new benchmark sets (cMET and PDE2A) that can be used in the community for validation purposes. Our findings also highlight the effectiveness of some AMBER forcefields, in particular, AMBER ff15ipq.

Preliminary modelling as the first stage of targeted organic synthesis

The review aims to present a classification and applicability analysis of methods for preliminary molecular modelling for targeted organic, catalytic and biocatalytic synthesis. The following three main approaches are considered as a primary classification of the methods: modelling of the target – ligand coordination without structural information on both the target and the resulting complex; calculations based on experimentally obtained structural information about the target; and dynamic simulation of the target – ligand complex and the reaction mechanism with calculation of the free energy of the reaction. The review is meant for synthetic chemists to be used as a guide for building an algorithm for preliminary modelling and synthesis of structures with specified properties.

The bibliography includes 353 references.

Identification of Optimal Ligand Growth Vectors Using an Alchemical Free-Energy Method

In this work, a novel method to rationally design inhibitors with improved steric contacts and enhanced binding free energies is presented. This new method uses alchemical single step perturbation calculations to rapidly optimize the van der Waals interactions of a small molecule in a protein-ligand complex in order to maximize its binding affinity. The results of the optimizer are used to predict beneficial growth vectors on the ligand, and good agreement is found between the predictions from the optimizer and a more rigorous free energy calculation, with a Spearman's rank order correlation of 0.59. The advantage of the method presented here is the significant speed up of over 10-fold compared to traditional free energy calculations and sublinear scaling with the number of growth vectors assessed. Where experimental data were available, mutations from hydrogen to a methyl group at sites highlighted by the optimizer were calculated with MBAR, and the mean unsigned error between experimental and calculated values of the binding free energy was 0.83 kcal/mol.

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.

Advances in the calculation of binding free energies

In recent years, calculations of binding affinities from molecular simulations seem to have matured significantly. While the number of applications of such methods in drug design and biotechnology increases, the number of truly new methodological developments decreases. This review provides an overview of the current status of the field as reflected in recent publications. The focus is on the challenges that remain when using endstate, alchemical and pathway methods. For endstate methods this is the calculation of entropic contributions. For alchemical methods there are unsolved problems associated with the solvation of the active site, sampling slow degrees of freedom and when modifying the net charge. For pathway methods achieving sufficient sampling remains challenging. New trends are also highlighted, including the use of pathway methods for the quantification of protein-protein interactions.

Optimization of Protein-Ligand Electrostatic Interactions Using an Alchemical Free-Energy Method

We present an explicit solvent alchemical free-energy method for optimizing the partial charges of a ligand to maximize the binding affinity with a receptor. This methodology can be applied to known ligand–protein complexes to determine an optimized set of ligand partial atomic changes. Three protein–ligand complexes have been optimized in this work: FXa, P38, and the androgen receptor. The sets of optimized charges can be used to identify design principles for chemical changes to the ligands which improve the binding affinity for all three systems. In this work, beneficial chemical mutations are generated from these principles and the resulting molecules tested using free-energy perturbation calculations. We show that three quarters of our chemical changes are predicted to improve the binding affinity, with an average improvement for the beneficial mutations of approximately 1 kcal/mol. In the cases where experimental data are available, the agreement between prediction and experiment is also good. The results demonstrate that charge optimization in explicit solvent is a useful tool for predicting beneficial chemical changes such as pyridinations, fluorinations, and oxygen to sulfur mutations.