Computational Macrocyclization: From de novo Macrocycle Generation to Binding Affinity Estimation

Uncovering the selectivity mechanism of phosphodiesterase 7A/8A inhibitors through computational studies

The expression of phosphodiesterase 7A (PDE7A) and phosphodiesterase 8A (PDE8) genes is integral to human signaling pathways, and the inhibition of PDE7A has been associated with the onset of various diseases, including effects on the immune system and nervous system. The development of PDE7 selective inhibitors can promote research on immune and nervous system diseases, such as multiple sclerosis, chronic inflammation, and autoimmune responses. PDE8A is expressed alongside PDE8B, and its inhibitory mechanism is still unclear. Studying the mechanisms of selective inhibitors against different PDE subtypes is crucial to prevent potential side effects, such as nausea and cardiac toxicity, and the sequence similarity of the two protein subtypes was 55.9%. Therefore, it is necessary to investigate the differences of both subtypes' ligand binding sites. Selective inhibitors of two proteins were chosen to summarize the reason for their selectivity through molecular docking, molecular dynamics simulation, alanine scanning mutagenesis, and MM-GBSA calculation. We found that Phe384PDE7A, Leu401PDE7A, Gln413PDE7A, Tyr419PDE7A, and Phe416PDE7A in the active site positively contribute to the selectivity towards PDE7A. Additionally, Asn729PDE8A, Phe767PDE8A, Gln778PDE8A, and Phe781PDE8A positively contribute to the selectivity towards PDE8A.

The maximal and current accuracy of rigorous protein-ligand binding free energy calculations

Computational techniques can speed up the identification of hits and accelerate the development of candidate molecules for drug discovery. Among techniques for predicting relative binding affinities, the most consistently accurate is free energy perturbation (FEP), a class of rigorous physics-based methods. However, uncertainty remains about how accurate FEP is and can ever be. Here, we present what we believe to be the largest publicly available dataset of proteins and congeneric series of small molecules, and assess the accuracy of the leading FEP workflow. To ascertain the limit of achievable accuracy, we also survey the reproducibility of experimental relative affinity measurements. We find a wide variability in experimental accuracy and a correspondence between binding and functional assays. When careful preparation of protein and ligand structures is undertaken, FEP can achieve accuracy comparable to experimental reproducibility. Throughout, we highlight reliable protocols that can help maximize the accuracy of FEP in prospective studies.

Macrocyclization of linear molecules by deep learning to facilitate macrocyclic drug candidates discovery

Interest in macrocycles as potential therapeutic agents has increased rapidly. Macrocyclization of bioactive acyclic molecules provides a potential avenue to yield novel chemical scaffolds, which can contribute to the improvement of the biological activity and physicochemical properties of these molecules. In this study, we propose a computational macrocyclization method based on Transformer architecture (which we name Macformer). Leveraging deep learning, Macformer explores the vast chemical space of macrocyclic analogues of a given acyclic molecule by adding diverse linkers compatible with the acyclic molecule. Macformer can efficiently learn the implicit relationships between acyclic and macrocyclic structures represented as SMILES strings and generate plenty of macrocycles with chemical diversity and structural novelty. In data augmentation scenarios using both internal ChEMBL and external ZINC test datasets, Macformer display excellent performance and generalisability. We showcase the utility of Macformer when combined with molecular docking simulations and wet lab based experimental validation, by applying it to the prospective design of macrocyclic JAK2 inhibitors.

Best practices for constructing, preparing, and evaluating protein-ligand binding affinity benchmarks [Article v0.1]

Free energy calculations are rapidly becoming indispensable in structure-enabled drug discovery programs. As new methods, force fields, and implementations are developed, assessing their expected accuracy on real-world systems (benchmarking) becomes critical to provide users with an assessment of the accuracy expected when these methods are applied within their domain of applicability, and developers with a way to assess the expected impact of new methodologies. These assessments require construction of a benchmark-a set of well-prepared, high quality systems with corresponding experimental measurements designed to ensure the resulting calculations provide a realistic assessment of expected performance when these methods are deployed within their domains of applicability. To date, the community has not yet adopted a common standardized benchmark, and existing benchmark reports suffer from a myriad of issues, including poor data quality, limited statistical power, and statistically deficient analyses, all of which can conspire to produce benchmarks that are poorly predictive of real-world performance. Here, we address these issues by presenting guidelines for (1) curating experimental data to develop meaningful benchmark sets, (2) preparing benchmark inputs according to best practices to facilitate widespread adoption, and (3) analysis of the resulting predictions to enable statistically meaningful comparisons among methods and force fields. We highlight challenges and open questions that remain to be solved in these areas, as well as recommendations for the collection of new datasets that might optimally serve to measure progress as methods become systematically more reliable. Finally, we provide a curated, versioned, open, standardized benchmark set adherent to these standards (PLBenchmarks) and an open source toolkit for implementing standardized best practices assessments (arsenic) for the community to use as a standardized assessment tool. While our main focus is free energy methods based on molecular simulations, these guidelines should prove useful for assessment of the rapidly growing field of machine learning methods for affinity prediction as well.

Machine learning-accelerated quantum mechanics-based atomistic simulations for industrial applications

Atomistic simulations have become an invaluable tool for industrial applications ranging from the optimization of protein-ligand interactions for drug discovery to the design of new materials for energy applications. Here we review recent advances in the use of machine learning (ML) methods for accelerated simulations based on a quantum mechanical (QM) description of the system. We show how recent progress in ML methods has dramatically extended the applicability range of conventional QM-based simulations, allowing to calculate industrially relevant properties with enhanced accuracy, at reduced computational cost, and for length and time scales that would have otherwise not been accessible. We illustrate the benefits of ML-accelerated atomistic simulations for industrial R&D processes by showcasing relevant applications from two very different areas, drug discovery (pharmaceuticals) and energy materials. Writing from the perspective of both a molecular and a materials modeling scientist, this review aims to provide a unified picture of the impact of ML-accelerated atomistic simulations on the pharmaceutical, chemical, and materials industries and gives an outlook on the exciting opportunities that could emerge in the future.

Structural basis for tailor-made selective PI3K α/β inhibitors: a computational perspective

PI3K α and β are Class IA PI3K isoforms that share a highly homologous ATP binding site, differing only in a few residues around the binding site.

Automated Design of Macrocycles for Therapeutic Applications: From Small Molecules to Peptides and Proteins

Macrocycles and cyclic peptides are increasingly attractive therapeutic modalities as they often have improved affinity, are able to bind to extended protein surfaces, and otherwise have favorable properties. Macrocyclization of a known binder may stabilize its bioactive conformation and improve its metabolic stability, cell permeability, and in certain cases oral bioavailability. Herein, we present implementation and application of an approach that automatically generates, evaluates, and proposes cyclizations utilizing a library of well-established chemical reactions and reagents. Using the three-dimensional (3D) conformation of the linear molecule in complex with a target protein as the starting point, this approach identifies attachment points, generates linkers, evaluates their geometric compatibility, and ranks the resulting molecules with respect to their predicted conformational stability and interactions with the target protein. As we show here with prospective and retrospective case studies, this procedure can be applied for the macrocyclization of small molecules and peptides and even PROteolysis TArgeting Chimeras (PROTACs) and proteins.

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 (R²) 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.

Improving the Accuracy of Protein Thermostability Predictions for Single Point Mutations

Accurately predicting the protein thermostability changes upon single point mutations in silico is a challenge that has implications for understanding diseases as well as industrial applications of protein engineering. Free energy perturbation (FEP) has been applied to predict the effect of single point mutations on protein stability for over 40 years and emerged as a potentially reliable prediction method with reasonable throughput. However, applications of FEP in protein stability calculations in industrial settings have been hindered by a number of limitations, including the inability to model mutations to and from prolines in which the bonded topology of the backbone is modified and the complexity in modeling charge-changing mutations. In this study, we have extended the FEP+ protocol to enable the accurate modeling of the effects on protein stability from proline mutations and from charge-changing mutations. We also evaluated the influence of the unfolded model in the stability calculations using increasingly longer peptides with native sequence and conformations. With the abovementioned improvements, the accuracy of FEP predictions of protein stability over a data set of 87 mutations on five different proteins has drastically improved compared with previous studies, with a mean unsigned error of 0.86 kcal/mol and root mean square error of 1.11 kcal/mol, comparable with the accuracy of previously published state-of-the-art small-molecule relative binding affinity calculations, which have been shown to be capable of driving discovery projects.

Adapting free energy perturbation simulations for large macrocyclic ligands: how to dissect contributions from direct binding and free ligand flexibility

Large and flexible ligands gain increasing interest in the development of bioactive agents. They challenge the applicability of computational ligand optimization strategies originally developed for small molecules. Free energy perturbation (FEP) is often used for predicting binding affinities of small molecule ligands, however, its use for more complex ligands remains limited. Herein, we report the structure-based design of peptide macrocycles targeting the protein binding site of human adaptor protein 14-3-3. We observe a surprisingly strong dependency of binding affinities on relatively small variations in substituent size. FEP was performed to rationalize observed trends. To account for insufficient convergence of FEP, restrained calculations were performed and complemented with extensive REST MD simulations of the free ligands. These calculations revealed that changes in affinity originate both from altered direct interactions and conformational changes of the free ligand. In addition, MD simulations provided the basis to rationalize unexpected trends in ligand lipophilicity. We also verified the anticipated interaction site and binding mode for one of the high affinity ligands by X-ray crystallography. The introduced fully-atomistic simulation protocol can be used to rationalize the development of structurally complex ligands which will support future ligand maturation efforts.

Best Practices for Alchemical Free Energy Calculations [Article v1.0]

Alchemical free energy calculations are a useful tool for predicting free energy differences associated with the transfer of molecules from one environment to another. The hallmark of these methods is the use of "bridging" potential energy functions representing alchemical intermediate states that cannot exist as real chemical species. The data collected from these bridging alchemical thermodynamic states allows the efficient computation of transfer free energies (or differences in transfer free energies) with orders of magnitude less simulation time than simulating the transfer process directly. While these methods are highly flexible, care must be taken in avoiding common pitfalls to ensure that computed free energy differences can be robust and reproducible for the chosen force field, and that appropriate corrections are included to permit direct comparison with experimental data. In this paper, we review current best practices for several popular application domains of alchemical free energy calculations performed with equilibrium simulations, in particular relative and absolute small molecule binding free energy calculations to biomolecular targets.

Large scale relative protein ligand binding affinities using non-equilibrium alchemy

Ligand binding affinity calculations based on molecular dynamics (MD) simulations and non-physical (alchemical) thermodynamic cycles have shown great promise for structure-based drug design. However, their broad uptake and impact is held back by the notoriously complex setup of the calculations. Only a few tools other than the free energy perturbation approach by Schrödinger Inc. (referred to as FEP+) currently enable end-to-end application. Here, we present for the first time an approach based on the open-source software pmx that allows to easily set up and run alchemical calculations for diverse sets of small molecules using the GROMACS MD engine. The method relies on theoretically rigorous non-equilibrium thermodynamic integration (TI) foundations, and its flexibility allows calculations with multiple force fields. In this study, results from the Amber and Charmm force fields were combined to yield a consensus outcome performing on par with the commercial FEP+ approach. A large dataset of 482 perturbations from 13 different protein-ligand datasets led to an average unsigned error (AUE) of 3.64 ± 0.14 kJ mol-1, equivalent to Schrödinger's FEP+ AUE of 3.66 ± 0.14 kJ mol-1. For the first time, a setup is presented for overall high precision and high accuracy relative protein-ligand alchemical free energy calculations based on open-source software.

Fast, efficient fragment-based coordinate generation for Open Babel

Rapidly predicting an accurate three dimensional geometry of a molecule is a crucial task for cheminformatics and across a wide range of molecular modeling. Consequently, developing a fast, accurate, and open implementation of structure prediction is necessary for reproducible cheminformatics research. We introduce a fragment-based coordinate generation implementation for Open Babel, a widely-used open source toolkit for cheminformatics. The new implementation improves speed and stereochemical accuracy, while retaining or improving accuracy of bond lengths, bond angles, and dihedral torsions. Input molecules are broken into fragments by cutting at rotatable bonds. The coordinates of fragments are set according to a fragment library, prepared from open crystallographic databases. Since the coordinates of multiple atoms are decided at once, coordinate prediction is accelerated over the previous rules-based implementation in Open Babel, as well as the widely-used distance geometry methods in RDKit. This new implementation will be beneficial for a wide range of applications, including computational property prediction in polymers, molecular materials and drug design.

Protein-Ligand Binding Free Energy Calculations with FEP

Accurate and reliable calculation of protein-ligand binding free energy is of central importance in computational biophysics and structure-based drug design. Among the various methods to calculate protein-ligand binding affinities, alchemical free energy perturbation (FEP) calculations performed by way of explicitly solvated molecular dynamics simulations (FEP/MD) provide a thermodynamically rigorous and complete description of the binding event and should in turn yield highly accurate predictions. Although the original theory of FEP was proposed more than 60 years ago, subsequent applications of FEP to compute protein-ligand binding free energies in the context of drug discovery projects over much of that time period was sporadic and generally unsuccessful. This was mainly due to the limited accuracy of the available force fields, inadequate sampling of the protein-ligand conformational space, complexity of simulation set up and analysis, and the large computational resources required to pursue such calculations. Over the past few years, there have been advances in computing power, classical force field accuracy, enhanced sampling algorithms, and simulation setup. This has led to newer FEP implementations such as the FEP+ technology developed by Schrödinger Inc., which has enabled accurate and reliable calculations of protein-ligand binding free energies and positioned free energy calculations to play a guiding role in small-molecule drug discovery. In this chapter, we outline the methodological advances in FEP+, including the OPLS3 force fields, the REST2 (Replica Exchange with Solute Tempering) enhanced sampling, the incorporation of REST2 sampling with conventional FEP (Free Energy Perturbation) through FEP/REST, and the advanced simulation setup and data analysis. The validation of FEP+ method in retrospective studies and the prospective applications in drug discovery projects are also discussed. We then present the recent extension of FEP+ method to handle challenging perturbations, including core-hopping transformations, macrocycle modifications, and reversible covalent inhibitor optimization. The limitations and pitfalls of the current FEP+ methodology and the best practices in real applications are also examined.

Modeling the value of predictive affinity scoring in preclinical drug discovery

Drug discovery is widely recognized to be a difficult and costly activity in large part due to the challenge of identifying chemical matter which simultaneously optimizes multiple properties, one of which is affinity for the primary biological target. Further, many of these properties are difficult to predict ahead of expensive and time-consuming compound synthesis and experimental testing. Here we highlight recent work to develop compound affinity prediction models, and extensively investigate the value such models may provide to preclinical drug discovery. We demonstrate that the ability of these models to improve the overall probability of success is crucially dependent on the shape of the error distribution, not just the root-mean-square error. In particular, while scoring more molecule ideas generally improves the probability of project success when the error distribution is Gaussian, fat-tail distributions such as a Cauchy distribution, can lead to a situation where scoring more ideas actually decreases the overall probability of success.

Predicting Binding Free Energies of PDE2 Inhibitors. The Difficulties of Protein Conformation

A congeneric series of 21 phosphodiesterase 2 (PDE2) inhibitors are reported. Crystal structures show how the molecules can occupy a 'top-pocket' of the active site. Molecules with small substituents do not enter the pocket, a critical leucine (Leu770) is closed and water molecules are present. Large substituents enter the pocket, opening the Leu770 conformation and displacing the waters. We also report an X-ray structure revealing a new conformation of the PDE2 active site domain. The relative binding affinities of these compounds were studied with free energy perturbation (FEP) methods and it represents an attractive real-world test case. In general, the calculations could predict the energy of small-to-small, or large-to-large molecule perturbations. However, accurately capturing the transition from small-to-large proved challenging. Only when using alternative protein conformations did results improve. The new X-ray structure, along with a modelled dimer, conferred stability to the catalytic domain during the FEP molecular dynamics (MD) simulations, increasing the convergence and thereby improving the prediction of ΔΔG of binding for some small-to-large transitions. In summary, we found the most significant improvement in results when using different protein structures, and this data set is useful for future free energy validation studies.

Computational Macrocyclization: From de novo Macrocycle Generation to Binding Affinity Estimation

Macrocycles play an increasing role in drug discovery, but their synthesis is often demanding. Computational tools that suggest macrocyclization based on a known binding mode and that estimate the binding affinity of these macrocycles could have a substantial impact on the medicinal chemistry design process. For both tasks, we established a workflow with high practical value. For five diverse pharmaceutical targets we show that the effect of macrocyclization on binding can be calculated robustly and accurately. Applying this method to macrocycles designed by LigMac, a search tool for de novo macrocyclization, our results suggest that we have a robust protocol in hand to design macrocycles and prioritize them prior to synthesis.