Matteo Aldeghi

Accurate calculation of the absolute free energy of binding for drug molecules

Free energy calculations based on molecular dynamics and thermodynamic cycles accurately reproduce experimental affinities of diverse bromodomain inhibitors.

Selective targeting of the BRG/PB1 bromodomains impairs embryonic and trophoblast stem cell maintenance.

PFI-3, a novel inhibitor targeting the bromodomains of essential components of the BAF/PBAF complex, affects the differentiation of ESC and TSC. Mammalian SWI/SNF [also called Brg/Brahma-associated factors (BAFs)] are evolutionarily conserved chromatin-remodeling complexes regulating gene transcription programs during development and stem cell differentiation. BAF complexes contain an ATP (adenosine 5′-triphosphate)–driven remodeling enzyme (either BRG1 or BRM) and multiple protein interaction domains including bromodomains, an evolutionary conserved acetyl lysine–dependent protein interaction motif that recruits transcriptional regulators to acetylated chromatin. We report a potent and cell active protein interaction inhibitor, PFI-3, that selectively binds to essential BAF bromodomains. The high specificity of PFI-3 was achieved on the basis of a novel binding mode of a salicylic acid head group that led to the replacement of water molecules typically maintained in other bromodomain inhibitor complexes. We show that exposure of embryonic stem cells to PFI-3 led to deprivation of stemness and deregulated lineage specification. Furthermore, differentiation of trophoblast stem cells in the presence of PFI-3 was markedly enhanced. The data present a key function of BAF bromodomains in stem cell maintenance and differentiation, introducing a novel versatile chemical probe for studies on acetylation-dependent cellular processes controlled by BAF remodeling complexes.

The Fragment Molecular Orbital Method Reveals New Insight into the Chemical Nature of GPCR–Ligand Interactions

Our interpretation of ligand-protein interactions is often informed by high-resolution structures, which represent the cornerstone of structure-based drug design. However, visual inspection and molecular mechanics approaches cannot explain the full complexity of molecular interactions. Quantum Mechanics approaches are often too computationally expensive, but one method, Fragment Molecular Orbital (FMO), offers an excellent compromise and has the potential to reveal key interactions that would otherwise be hard to detect. To illustrate this, we have applied the FMO method to 18 Class A GPCR-ligand crystal structures, representing different branches of the GPCR genome. Our work reveals key interactions that are often omitted from structure-based descriptions, including hydrophobic interactions, nonclassical hydrogen bonds, and the involvement of backbone atoms. This approach provides a more comprehensive picture of receptor-ligand interactions than is currently used and should prove useful for evaluation of the chemical nature of ligand binding and to support structure-based drug design.

Predictions of Ligand Selectivity from Absolute Binding Free Energy Calculations.

Binding selectivity is a requirement for the development of a safe drug, and it is a critical property for chemical probes used in preclinical target validation. Engineering selectivity adds considerable complexity to the rational design of new drugs, as it involves the optimization of multiple binding affinities. Computationally, the prediction of binding selectivity is a challenge, and generally applicable methodologies are still not available to the computational and medicinal chemistry communities. Absolute binding free energy calculations based on alchemical pathways provide a rigorous framework for affinity predictions and could thus offer a general approach to the problem. We evaluated the performance of free energy calculations based on molecular dynamics for the prediction of selectivity by estimating the affinity profile of three bromodomain inhibitors across multiple bromodomain families, and by comparing the results to isothermal titration calorimetry data. Two case studies were considered. In the first one, the affinities of two similar ligands for seven bromodomains were calculated and returned excellent agreement with experiment (mean unsigned error of 0.81 kcal/mol and Pearson correlation of 0.75). In this test case, we also show how the preferred binding orientation of a ligand for different proteins can be estimated via free energy calculations. In the second case, the affinities of a broad-spectrum inhibitor for 22 bromodomains were calculated and returned a more modest accuracy (mean unsigned error of 1.76 kcal/mol and Pearson correlation of 0.48); however, the reparametrization of a sulfonamide moiety improved the agreement with experiment.

Statistical Analysis on the Performance of Molecular Mechanics Poisson–Boltzmann Surface Area versus Absolute Binding Free Energy Calculations: Bromodomains as a Case Study

Binding free energy calculations that make use of alchemical pathways are becoming increasingly feasible thanks to advances in hardware and algorithms. Although relative binding free energy (RBFE) calculations are starting to find widespread use, absolute binding free energy (ABFE) calculations are still being explored mainly in academic settings due to the high computational requirements and still uncertain predictive value. However, in some drug design scenarios, RBFE calculations are not applicable and ABFE calculations could provide an alternative. Computationally cheaper end-point calculations in implicit solvent, such as molecular mechanics Poisson–Boltzmann surface area (MMPBSA) calculations, could too be used if one is primarily interested in a relative ranking of affinities. Here, we compare MMPBSA calculations to previously performed absolute alchemical free energy calculations in their ability to correlate with experimental binding free energies for three sets of bromodomain–inhibitor pairs. Different MMPBSA approaches have been considered, including a standard single-trajectory protocol, a protocol that includes a binding entropy estimate, and protocols that take into account the ligand hydration shell. Despite the improvements observed with the latter two MMPBSA approaches, ABFE calculations were found to be overall superior in obtaining correlation with experimental affinities for the test cases considered. A difference in weighted average Pearson () and Spearman () correlations of 0.25 and 0.31 was observed when using a standard single-trajectory MMPBSA setup ( = 0.64 and = 0.66 for ABFE; = 0.39 and = 0.35 for MMPBSA). The best performing MMPBSA protocols returned weighted average Pearson and Spearman correlations that were about 0.1 inferior to ABFE calculations: = 0.55 and = 0.56 when including an entropy estimate, and = 0.53 and = 0.55 when including explicit water molecules. Overall, the study suggests that ABFE calculations are indeed the more accurate approach, yet there is also value in MMPBSA calculations considering the lower compute requirements, and if agreement to experimental affinities in absolute terms is not of interest. Moreover, for the specific protein–ligand systems considered in this study, we find that including an explicit ligand hydration shell or a binding entropy estimate in the MMPBSA calculations resulted in significant performance improvements at a negligible computational cost.

Using the fragment molecular orbital method to investigate agonist-orexin-2 receptor interactions.

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity and is essential for an efficient structure-based drug discovery (SBDD) process. Clearly, to begin SBDD, a structure is needed, and although there has been fantastic progress in solving G-protein-coupled receptor (GPCR) crystal structures, the process remains quite slow and is not currently feasible for every GPCR or GPCR–ligand complex. This situation significantly limits the ability of X-ray crystallography to impact the drug discovery process for GPCR targets in ‘real-time’ and hence there is still a need for other practical and cost-efficient alternatives. We present here an approach that integrates our previously described hierarchical GPCR modelling protocol (HGMP) and the fragment molecular orbital (FMO) quantum mechanics (QM) method to explore the interactions and selectivity of the human orexin-2 receptor (OX2R) and its recently discovered nonpeptidic agonists. HGMP generates a 3D model of GPCR structures and its complexes with small molecules by applying a set of computational methods. FMO allows ab initio approaches to be applied to systems that conventional QM methods would find challenging. The key advantage of FMO is that it can reveal information on the individual contribution and chemical nature of each residue and water molecule to the ligand binding that normally would be difficult to detect without QM. We illustrate how the combination of both techniques provides a practical and efficient approach that can be used to analyse the existing structure–function relationships (SAR) and to drive forward SBDD in a real-world example for which there is no crystal structure of the complex available.

Application of an Integrated GPCR SAR-Modeling Platform To Explain the Activation Selectivity of Human 5-HT2C over 5-HT2B.

Agonism of the 5-HT2C serotonin receptor has been associated with the treatment of a number of diseases including obesity, psychiatric disorders, sexual health, and urology. However, the development of effective 5-HT2C agonists has been hampered by the difficulty in obtaining selectivity over the closely related 5-HT2B receptor, agonism of which is associated with irreversible cardiac valvulopathy. Understanding how to design selective agonists requires exploration of the structural features governing the functional uniqueness of the target receptor relative to related off targets. X-ray crystallography, the major experimental source of structural information, is a slow and challenging process for integral membrane proteins, and so is currently not feasible for every GPCR or GPCR-ligand complex. Therefore, the integration of existing ligand SAR data with GPCR modeling can be a practical alternative to provide this essential structural insight. To demonstrate this, we integrated SAR data from 39 azepine series 5-HT2C agonists, comprising both selective and unselective examples, with our hierarchical GPCR modeling protocol (HGMP). Through this work we have been able to demonstrate how relatively small differences in the amino acid sequences of GPCRs can lead to significant differences in secondary structure and function, as supported by experimental data. In particular, this study suggests that conformational differences in the tilt of TM7 between 5-HT2B and 5-HT2C, which result from differences in interhelical interactions, may be the major source of selectivity in G-protein activation between these two receptors. Our approach also demonstrates how the use of GPCR models in conjunction with SAR data can be used to explain activity cliffs.

Absolute Alchemical Free Energy Calculations for Ligand Binding: A Beginner's Guide.

Many thermodynamic quantities can be extracted from computer simulations that generate an ensemble of microstates according to the principles of statistical mechanics. Among these quantities is the free energy of binding of a small molecule to a macromolecule, such as a protein. Here, we present an introductory overview of a protocol that allows for the estimation of ligand binding free energies via molecular dynamics simulations. While we focus on the binding of organic molecules to proteins, the approach is in principle transferable to any pair of molecules.

Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties

Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.

Exploring GPCR-Ligand Interactions with the Fragment Molecular Orbital (FMO) Method

The understanding of binding interactions between any protein and a small molecule plays a key role in the rationalization of affinity and selectivity. It is essential for an efficient structure-based drug design (SBDD) process. FMO enables ab initio approaches to be applied to systems that conventional quantum-mechanical (QM) methods would find challenging. The key advantage of the Fragment Molecular Orbital Method (FMO) is that it can reveal atomistic details about the individual contributions and chemical nature of each residue and water molecule toward ligand binding which would otherwise be difficult to detect without using QM methods. In this chapter, we demonstrate the typical use of FMO to analyze 19 crystal structures of β1 and β2 adrenergic receptors with their corresponding agonists and antagonists.