Mapping of allosteric druggable sites in activation-associated conformers of the M2 muscarinic receptor

Application of Mixed-Solvent Molecular Dynamics Simulations for Prediction of Allosteric Sites on G Protein-Coupled Receptors

The development of small molecule allosteric modulators acting at G protein-coupled receptors (GPCRs) is becoming increasingly attractive. Such compounds have advantages over traditional drugs acting at orthosteric sites on these receptors, in particular target specificity. However, the number and locations of druggable allosteric sites within most clinically relevant GPCRs are unknown. In the present study, we describe the development and application of a mixed-solvent molecular dynamics (MixMD)-based method for the identification of allosteric sites on GPCRs. The method employs small organic probes with druglike qualities to identify druggable hotspots in multiple replicate short-timescale simulations. As proof of principle, we first applied the method retrospectively to a test set of five GPCRs (cannabinoid receptor type 1, C-C chemokine receptor type 2, M2 muscarinic receptor, P2Y purinoceptor 1, and protease-activated receptor 2) with known allosteric sites in diverse locations. This resulted in the identification of the known allosteric sites on these receptors. We then applied the method to the μ-opioid receptor. Several allosteric modulators for this receptor are known, although the binding sites for these modulators are not known. The MixMD-based method revealed several potential allosteric sites on the mu-opioid receptor. Implementation of the MixMD-based method should aid future efforts in the structure-based drug design of drugs targeting allosteric sites on GPCRs. SIGNIFICANCE STATEMENT: Allosteric modulation of G protein-coupled receptors (GPCRs) has the potential to provide more selective drugs. However, there are limited structures of GPCRs bound to allosteric modulators, and obtaining such structures is problematic. Current computational methods utilize static structures and therefore may not identify hidden or cryptic sites. Here we describe the use of small organic probes and molecular dynamics to identify druggable allosteric hotspots on GPCRs. The results reinforce the importance of protein dynamics in allosteric site identification.

Multi-state Model-Based Identification of Cryptic Allosteric Sites on Human Serotonin Transporter

Serotonin transporter (SERT) plays a fundamental role in taking the synaptic cleft serotonin back to the presynaptic neuron. The discovery of allosteric SERT modulators represents the next-generation medication for psychiatric disorders such as depression. Here, based on the cryo-EM structures of ibogaine in complex with SERT in distinct conformations, the multiple functional structures of the transporter bound to serotonin, including outward-open (OO_holo), outward-occluded (OC_holo), and inward-open (IO_holo and IO_holo'), were carefully characterized by induced-fit docking Gaussian-accelerated molecular dynamics (IFD-GaMD) simulation and the free-energy landscape analysis. Further MM/GBSA binding free energy, per-residue contribution, and molecular interaction fingerprint calculations revealed the interaction variations of serotonin with SERT in functional structures, which confirmed the allostery of SERT during serotonin reuptake. Moreover, five unique cryptic allosteric sites, which are druggable and capable of targeting by small molecules, were identified on the characterized multistate structures. These results provide structural and energetic information for the molecular mechanism of serotonin reuptake and will provide opportunities for the development of novel therapeutics based on the identified new allosteric sites on SERT.

Allosteric Regulation of G-Protein-Coupled Receptors: From Diversity of Molecular Mechanisms to Multiple Allosteric Sites and Their Ligands

Allosteric regulation is critical for the functioning of G protein-coupled receptors (GPCRs) and their signaling pathways. Endogenous allosteric regulators of GPCRs are simple ions, various biomolecules, and protein components of GPCR signaling (G proteins and β-arrestins). The stability and functional activity of GPCR complexes is also due to multicenter allosteric interactions between protomers. The complexity of allosteric effects caused by numerous regulators differing in structure, availability, and mechanisms of action predetermines the multiplicity and different topology of allosteric sites in GPCRs. These sites can be localized in extracellular loops; inside the transmembrane tunnel and in its upper and lower vestibules; in cytoplasmic loops; and on the outer, membrane-contacting surface of the transmembrane domain. They are involved in the regulation of basal and orthosteric agonist-stimulated receptor activity, biased agonism, GPCR-complex formation, and endocytosis. They are targets for a large number of synthetic allosteric regulators and modulators, including those constructed using molecular docking. The review is devoted to the principles and mechanisms of GPCRs allosteric regulation, the multiplicity of allosteric sites and their topology, and the endogenous and synthetic allosteric regulators, including autoantibodies and pepducins. The allosteric regulation of chemokine receptors, proteinase-activated receptors, thyroid-stimulating and luteinizing hormone receptors, and beta-adrenergic receptors are described in more detail.

Unbiased molecular dynamics simulation of a first-in-class small molecule inhibitor binds to oncostatin M

Small molecule inhibitors (SMIs) targeting oncostatin M (OSM) signaling pathway represent new therapeutics to combat cancer, inflammatory bowel disease (IBD) and CNS disease. Recently, the first-in-class SMI named SMI-10B that target OSM and block its interaction with receptor (OSMR) were reported. However, the binding pocket and interaction mode of the compound on OSM remain poorly understood, which hampering the rational design of SMIs that target OSM. Here, using SMI-10B as a probe, the multiple pockets on OSM for small molecules binding were extensively explored by unbiased molecular dynamics (MD) simulations. Then, the near-native structure of the complex was identified by molecular mechanics generalized Born surface area (MM/GBSA) binding energy funnel. Moreover, the binding stabilities of the protein-ligand complexes in near- and non-native conformations were verified by additional independent MD runs and absolute free energy perturbation (FEP) calculation. In summary, the unique feature of SMI-10B spontaneously binds to OSM characterized here not only provide detailed information for understanding the molecular mechanism of SMI-10B binding to OSM, but also will facilitate the rational design of novel and more potent SMIs to block OSM signaling.

Sigmoid Accelerated Molecular Dynamics: An Efficient Enhanced Sampling Method for Biosystems

Gaussian accelerated molecular dynamics (GaMD) is recognized as a popular enhanced sampling method for tackling long-standing challenges in biomolecular simulations. Inspired by GaMD, Sigmoid accelerated molecular dynamics (SaMD) is proposed in this work by adding a Sigmoid boost potential to improve the balance between the highest acceleration and accurate reweighting. Compared with GaMD, SaMD extends the accessible time scale and improves the computational efficiency as tested in three tasks. In the alanine dipeptide task, SaMD can produce the free energy landscape with better accuracy and efficiency. In the chignolin folding task, the estimated Gibbs free energy difference can converge to the experimental value ∼30% faster. In the protein-ligand binding task, the bound conformations are closer to the crystal structure with a minimal ligand root-mean-square deviation of 1.7 Å. The binding of the ligand XK263 to the HIV protease is reproduced by SaMD in ∼60% less simulation time.

Conservation of Allosteric Ligand Binding Sites in G-Protein Coupled Receptors

Despite the growing number of G protein-coupled receptor (GPCR) structures, only 39 structures have been cocrystallized with allosteric inhibitors. These structures have been studied by protein mapping using the FTMap server, which determines the clustering of small organic probe molecules distributed on the protein surface. The method has found druggable sites overlapping with the cocrystallized allosteric ligands in 21 GPCR structures. Mapping of Alphafold2 generated models of these proteins confirms that the same sites can be identified without the presence of bound ligands. We then mapped the 394 GPCR X-ray structures available at the time of the analysis (September 2020). Results show that for each of the 21 structures with bound ligands there exist many other GPCRs that have a strong binding hot spot at the same location, suggesting potential allosteric sites in a large variety of GPCRs. These sites cluster at nine distinct locations, and each can be found in many different proteins. However, ligands binding at the same location generally show little or no similarity, and the amino acid residues interacting with these ligands also differ. Results confirm the possibility of specifically targeting these sites across GPCRs for allosteric modulation and help to identify the most likely binding sites among the limited number of potential locations. The FTMap server is available free of charge for academic and governmental use at https://ftmap.bu.edu/.

Mapping the binding sites of challenging drug targets

An increasing number of medically important proteins are challenging drug targets because their binding sites are too shallow or too polar, are cryptic and thus not detectable without a bound ligand or located in a protein-protein interface. While such proteins may not bind druglike small molecules with sufficiently high affinity, they are frequently druggable using novel therapeutic modalities. The need for such modalities can be determined by experimental or computational fragment based methods. Computational mapping by mixed solvent molecular dynamics simulations or the FTMap server can be used to determine binding hot spots. The strength and location of the hot spots provide very useful information for selecting potentially successful approaches to drug discovery.

Unmasking allosteric binding sites: Novel targets for GPCR drug discovery

INTRODUCTION

Unexpected non-apparent and hidden allosteric binding sites are non-classical and non-apparent allosteric centers in 3-D X-ray protein structures until orthosteric or allosteric ligands bind to them. The orthosteric center of one protomer that modulates binding centers of the other protomers within an oligomer is also an unexpected allosteric site. Furthermore, another partner protein can also produce these effects, acting as an unexpected allosteric modulator.

AREAS COVERED

This review summarizes both classical and non-classical allosterism. The authors focus on G protein-coupled receptor (GPCR) oligomers as a paradigm of allosteric molecules. Moreover, they show several examples of unexpected allosteric sites such as hidden allosteric sites in a protomer that appear after the interaction with other molecules and the allosterism exerted between orthosteric sites within GPCR oligomer, emphasizing on the allosteric modulations that can occur between binding sites.

EXPERT OPINION

The study of these new non-classical allosteric sites will expand the diversity of allosteric control on the function of orthosteric sites within proteins, whether GPCRs or other receptors, enzymes or transporters. Moreover, the design of new drugs targeting these hidden allosteric sites or already known orthosteric sites acting as allosteric sites in protein homo- or hetero-oligomers will increase the therapeutic potential of allosterism.

IgG1-b12-HIV-gp120 Interface in Solution: A Computational Study

The use of broadly neutralizing antibodies against human immunodeficiency virus type 1 (HIV-1) has been shown to be a promising therapeutic modality in the prevention of HIV infection. Understanding the b12-gp120 binding mechanism under physiological conditions may assist the development of more broadly effective antibodies. In this work, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein gp120 of HIV-1 and the IgG1-b12 mAb are studied. Accelerated molecular dynamics (aMD) and ab initio hybrid molecular dynamics have been combined to determine the most persistent interactions between the most populated conformations of the antibody-antigen complex under physiological conditions. The results show the most persistent receptor-binding mapping in the conformations of the antibody-antigen interface in solution. The binding-free-energy decomposition reveals a small enhancement in the contribution played by the CDR-H3 region to the b12-gp120 interface compared to the crystal structure.

Computational Methods for Structure-Based Drug Design Through System Biology

The advances in computational chemistry and biology, computer science, structural biology, and molecular biology go in parallel with the rapid progress in target-based systems. This technique has become a powerful tool in medicinal chemistry for the identification of hit molecules. The recent developments in target-based systems have played a major role in the creation of libraries of compounds, and it has also been widely applied for the design of molecular docking methods. The main advantage of this method is that it hits the fragment that has the strongest binding, has relatively small size, and leads to better compounds in terms of pharmacokinetic properties when compared with virtual screening (VS) and high-throughput screening (HTS) hits. De novo design is an essential aspect of target-based systems and requires the synthesis of chemical to allow the design of promising compound.

Hidden allosteric sites and De-Novo drug design

INTRODUCTION

Hidden allosteric sites are not visible in apo-crystal structures, but they may be visible in holo-structures when a certain ligand binds and maintains the ligand intended conformation. Several computational and experimental techniques have been used to investigate these hidden sites but identifying them remains a challenge.

AREAS COVERED

This review provides a summary of the many theoretical approaches for predicting hidden allosteric sites in disease-related proteins. Furthermore, promising cases have been thoroughly examined to reveal the hidden allosteric site and its modulator.

EXPERT OPINION

In the recent past, with the development in scientific techniques and bioinformatics tools, the number of drug targets for complex human diseases has significantly increased but unfortunately most of these targets are undruggable due to several reasons. Alternative strategies such as finding cryptic (hidden) allosteric sites are an attractive approach for exploitation of the discovery of new targets. These hidden sites are difficult to recognize compared to allosteric sites, mainly due to a lack of visibility in the crystal structure. In our opinion, after many years of development, MD simulations are finally becoming successful for obtaining a detailed molecular description of drug-target interaction.

Probe Confined Dynamic Mapping for G Protein-Coupled Receptor Allosteric Site Prediction

Targeting G protein-coupled receptors (GPCRs) through allosteric sites offers advantages over orthosteric sites in identifying drugs with increased selectivity and potentially reduced side effects. In this study, we developed a probe confined dynamic mapping protocol that allows the prediction of allosteric sites at both the GPCR extracellular and intracellular sides, as well as at the receptor-lipid interface. The applied harmonic wall potential enhanced sampling of probe molecules in a selected area of a GPCR while preventing membrane distortion in molecular dynamics simulations. The specific probes derived from GPCR allosteric ligand structures performed better in allosteric site mapping compared to commonly used cosolvents. The M₂ muscarinic, β₂ adrenergic, and P₂Y₁ purinergic receptors were selected for the protocol's retrospective validation. The protocol was next validated prospectively to locate the binding site of [5-fluoro-4-(hydroxymethyl)-2-methoxyphenyl]-(4-fluoro-1H-indol-1-yl)methanone at the D₂ dopamine receptor, and subsequent mutagenesis confirmed the prediction. The protocol provides fast and efficient prediction of key amino acid residues surrounding allosteric sites in membrane proteins and facilitates the structure-based design of allosteric modulators.

Unravelling the molecular interactions between the SARS-CoV-2 RBD spike protein and various specific monoclonal antibodies

Vaccination against SARS-CoV-2 just started in most of the countries. However, the development of specific vaccines against SARS-CoV-2 is not the only approach to control the virus and monoclonal antibodies (mAbs) start to merit special attention as a therapeutic option to treat COVID-19 disease. Here, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein of SARS-CoV-2 (S protein) with two mAbs (CR3022 and S309) and the ACE2 cell receptor are studied as the main representatives of three different epitopes on the RBD of S protein. The combined approach of 1 μs accelerated molecular dynamics (aMD) and ab-initio hybrid molecular dynamics is used to identify the most predominant interactions under physiological conditions. Results allow to determine the main receptor-binding mapping, hydrogen bonding network and salt bridges in the most populated antigen-antibody interface conformations. The deep knowledge on the protein-protein interactions involving mAbs and ACE2 receptor with the spike glycoprotein of SARS-CoV-2 increases background knowledge to speed up the development of new vaccines and therapeutic drugs.

Tuning Down the Pain - An Overview of Allosteric Modulation of Opioid Receptors: Mechanisms of Modulation, Allosteric Sites, Modulator Syntheses

Opioid signaling plays a central role in pain perception. As such, it remains the main target in the development of antinociceptive agents, despite serious side effects involved. In recent years, hopes for improved opioid painkillers are rising, together with our understanding of allosterism and biased signaling mechanisms. In this review, we focus on recently discovered allosteric modulators of opioid receptors, insights into phenomena underlying their action, as well as on how they extend our understanding of mechanisms of previously known compounds. A brief overlook of their synthesis is also presented.

Dynamical important residue network (DIRN): network inference via conformational change

MOTIVATION

Protein residue interaction network has emerged as a useful strategy to understand the complex relationship between protein structures and functions and how functions are regulated. In a residue interaction network, every residue is used to define a network node, adding noises in network post-analysis and increasing computational burden. In addition, dynamical information is often necessary in deciphering biological functions.

RESULTS

We developed a robust and efficient protein residue interaction network method, termed dynamical important residue network, by combining both structural and dynamical information. A major departure from previous approaches is our attempt to identify important residues most important for functional regulation before a network is constructed, leading to a much simpler network with the important residues as its nodes. The important residues are identified by monitoring structural data from ensemble molecular dynamics simulations of proteins in different functional states. Our tests show that the new method performs well with overall higher sensitivity than existing approaches in identifying important residues and interactions in tested proteins, so it can be used in studies of protein functions to provide useful hypotheses in identifying key residues and interactions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization

BACKGROUND

Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug.

METHODS

In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan.

RESULTS

A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given.

CONCLUSION

The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.

Analysis of tractable allosteric sites in G protein-coupled receptors

Allosteric modulation of G protein-coupled receptors represent a promising mechanism of pharmacological intervention. Dramatic developments witnessed in the structural biology of membrane proteins continue to reveal that the binding sites of allosteric modulators are widely distributed, including along protein surfaces. Here we restrict consideration to intrahelical and intracellular sites together with allosteric conformational locks, and show that the protein mapping tools FTMap and FTSite identify 83% and 88% of such experimentally confirmed allosteric sites within the three strongest sites found. The methods were also able to find partially hidden allosteric sites that were not fully formed in X-ray structures crystallized in the absence of allosteric ligands. These results confirm that the intrahelical sites capable of binding druglike allosteric modulators are among the strongest ligand recognition sites in a large fraction of GPCRs and suggest that both FTMap and FTSite are useful tools for identifying allosteric sites and to aid in the design of such compounds in a range of GPCR targets.

GPCR Allosteric Modulator Discovery

G protein-coupled receptors (GPCRs) influence virtually every aspect of human physiology; about one-third of all marketed drugs target members of this family. GPCR allosteric ligands hold the promise of improved subtype selectivity, spatiotemporal sensitivity, and possible biased property over typical orthosteric ligands. However, only a small number of GPCR allosteric ligands have been approved as drugs or in clinical trials since the discovery process is very challenging. The rapid development of GPCR structural biology leads to the discovery of several allosteric sites and sheds light on understanding the mechanism of GPCR allosteric ligands, which is critical for discovering novel therapeutics. This book chapter summarized different GPCR allosteric modulating mechanisms and discussed validated mechanisms based on allosteric modulator-GPCR complex structures.

Intracellular Receptor Modulation: Novel Approach to Target GPCRs

Recent crystal structures of multiple G protein-coupled receptors (GPCRs) have revealed a highly conserved intracellular pocket that can be used to modulate these receptors from the inside. This novel intracellular site partially overlaps with the G protein and β-arrestin binding site, providing a new manner of pharmacological intervention. Here we provide an update of the architecture and function of the intracellular region of GPCRs, until now portrayed as the signaling domain. We review the available evidence on the presence of intracellular binding sites among chemokine receptors and other class A GPCRs, as well as different strategies to target it, including small molecules, pepducins, and nanobodies. Finally, the potential advantages of intracellular (allosteric) ligands over orthosteric ligands are also discussed.

Ensemble docking to difficult targets in early-stage drug discovery: Methodology and application to fibroblast growth factor 23

Ensemble docking is now commonly used in early-stage in silico drug discovery and can be used to attack difficult problems such as finding lead compounds which can disrupt protein-protein interactions. We give an example of this methodology here, as applied to fibroblast growth factor 23 (FGF23), a protein hormone that is responsible for regulating phosphate homeostasis. The first small-molecule antagonists of FGF23 were recently discovered by combining ensemble docking with extensive experimental target validation data (Science Signaling, 9, 2016, ra113). Here, we provide a detailed account of how ensemble-based high-throughput virtual screening was used to identify the antagonist compounds discovered in reference (Science Signaling, 9, 2016, ra113). Moreover, we perform further calculations, redocking those antagonist compounds identified in reference (Science Signaling, 9, 2016, ra113) that performed well on drug-likeness filters, to predict possible binding regions. These predicted binding modes are rescored with the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) approach to calculate the most likely binding site. Our findings suggest that the antagonist compounds antagonize FGF23 through the disruption of protein-protein interactions between FGF23 and fibroblast growth factor receptor (FGFR).

In silico mapping of allosteric ligand binding sites in type-1 cannabinoid receptor

The recent resolution of the crystal structure of type-1 cannabinoid receptor (CB₁ ) and the discovery of novel modulators for this target open the way to the possibility of elucidating the structural requirements for CB₁ binding, and thereby facilitate a rational drug design. Compounds that target the orthosteric site of CB₁ in some cases have shown side effects. Allosteric modulators could potentially avoid these side effects by influencing binding and/or efficacy of orthosteric ligands. Here, we summarize and compare previous data on different putative allosteric binding sites observed in CB₁ homology models with an in silico docking study of the recently published crystal structure of the same receptor on endogenous and natural hydrophobic ligands that act as positive allosteric modulators and negative allosteric modulators of CB₁ . In particular, a lipid-exposed pocket targeted by most of the tested molecules is reported and discussed.

Mapping the allosteric sites of the A2A adenosine receptor

The A2A adenosine receptor (A2A AR) is a G protein-coupled receptor that is pharmacologically targeted for the treatment of inflammation, sepsis, cancer, neurodegeneration, and Parkinson's disease. Recently, we applied long-timescale molecular dynamics simulations on two ligand-free receptor conformations, starting from the agonist-bound (PDB ID: 3QAK) and antagonist-bound (PDB ID: 3EML) X-ray structures. This analysis revealed four distinct conformers of the A2A AR: the active, intermediate 1, intermediate 2, and inactive. In this study, we apply the fragment-based mapping algorithm, FTMap, on these receptor conformations to uncover five non-orthosteric sites on the A2A AR. Two sites that are identified in the active conformation are located in the intracellular region of the transmembrane helices (TM) 3/TM4 and the G protein-binding site in the intracellular region between TM2/TM3/TM6/TM7. Three sites are identified in the intermediate 1 and intermediate 2 conformations, annexing a site in the lipid interface of TM5/TM6. Five sites are identified in the inactive conformation, comprising a site in the intracellular region of TM1/TM7 and in the extracellular region of TM3/TM4 of the A2A AR. We postulate that these sites on the A2A AR be screened for allosteric modulators for the treatment of inflammatory and neurological diseases.

Molecular Dynamics Simulation and Prediction of Druggable Binding Sites

Binding site identification and druggability evaluation are two essential steps in structure-based drug design. A druggable binding site tends to have high binding affinity to drug-like molecules. Predicting such sites can have a significant impact on a drug design campaign. This chapter focuses on summarizing the different methods that are used to predict druggable binding sites. The chapter also discusses the importance of including protein flexibility in the search process and the use of molecular dynamics simulations to address this aspect. Case studies from the literature are also summarized and discussed. We hope that this chapter would provide an overview on the different methods employed in binding site identification evaluation.

Identification of a Different Agonist-Binding Site and Activation Mechanism of the Human P2Y1 Receptor

The human P2Y1 receptor (P2Y1R) is a purinergic G-protein-coupled receptor (GPCR) that functions as a receptor for adenosine 5'-diphosphate (ADP). An antagonist of P2Y1R might potentially have antithrombotic effects, whereas agonists might serve as antidiabetic agents. On the basis of the antagonist-bound MRS2500-P2Y1R crystal structure, we constructed computational models of apo-P2Y1R and the agonist-receptor complex 2MeSADP-P2Y1R. We then performed conventional molecular dynamics (cMD) and accelerated molecular dynamics (aMD) simulations to study the conformational dynamics after binding with agonist/antagonist as well as the P2Y1R activation mechanism. We identified a new agonist-binding site of P2Y1R that is consistent with previous mutagenesis data. This new site is deeper than those of the agonist ADP in the recently simulated ADP-P2Y1R structure and the antagonist MRS2500 in the MRS2500-P2Y1R crystal structure. During P2Y1R activation, the cytoplasmic end of helix VI shifts outward 9.1 Å, the Ser1463.47-Tyr2375.58 hydrogen bond breaks, a Tyr2375.58-Val2626.37 hydrogen bond forms, and the conformation of the χ1 rotamer of Phe2696.44 changes from parallel to perpendicular to helix VI. The apo-P2Y1R system and the MRS2500-P2Y1R system remain inactive. The newly identified agonist binding site and activation mechanism revealed in this study may aid in the design of P2Y1R antagonists/agonists as antithrombotic/antidiabetic agents, respectively.

A computationally identified compound antagonizes excess FGF-23 signaling in renal tubules and a mouse model of hypophosphatemia

Fibroblast growth factor-23 (FGF-23) interacts with a binary receptor complex composed of α-Klotho (α-KL) and FGF receptors (FGFRs) to regulate phosphate and vitamin D metabolism in the kidney. Excess FGF-23 production, which causes hypophosphatemia, is genetically inherited or occurs with chronic kidney disease. Among other symptoms, hypophosphatemia causes vitamin D deficiency and the bone-softening disorder rickets. Current therapeutics that target the receptor complex have limited utility clinically. Using a computationally driven, structure-based, ensemble docking and virtual high-throughput screening approach, we identified four novel compounds predicted to selectively inhibit FGF-23-induced activation of the FGFR/α-KL complex. Additional modeling and functional analysis found that Zinc13407541 bound to FGF-23 and disrupted its interaction with the FGFR1/α-KL complex; experiments in a heterologous cell expression system showed that Zinc13407541 selectivity inhibited α-KL-dependent FGF-23 signaling. Zinc13407541 also inhibited FGF-23 signaling in isolated renal tubules ex vivo and partially reversed the hypophosphatemic effects of excess FGF-23 in a mouse model. These chemical probes provide a platform to develop lead compounds to treat disorders caused by excess FGF-23.

Exploring Flexibility of Progesterone Receptor Ligand Binding Domain Using Molecular Dynamics

Progesterone receptor (PR), a member of nuclear receptor (NR) superfamily, plays a vital role for female reproductive tissue development, differentiation and maintenance. PR ligand, such as progesterone, induces conformation changes in PR ligand binding domain (LBD), thus mediates subsequent gene regulation cascades. PR LBD may adopt different conformations upon an agonist or an antagonist binding. These different conformations would trigger distinct transcription events. Therefore, the dynamics of PR LBD would be of general interest to biologists for a deep understanding of its structure-function relationship. However, no apo-form (non-ligand bound) of PR LBD model has been proposed either by experiments or computational methods so far. In this study, we explored the structural dynamics of PR LBD using molecular dynamics simulations and advanced sampling tools in both ligand-bound and the apo-forms. Resolved by the simulation study, helix 11, helix 12 and loop 895–908 (the loop between these two helices) are quite flexible in antagonistic conformation. Several residues, such as Arg899 and Glu723, could form salt-bridging interaction between helix 11 and helix 3, and are important for the PR LBD dynamics. And we also propose that helix 12 in apo-form PR LBD, not like other NR LBDs, such as human estrogen receptor α (ERα) LBD, may not adopt a totally extended conformation. With the aid of umbrella sampling and metadynamics simulations, several stable conformations of apo-form PR LBD have been sampled, which may work as critical structural models for further large scale virtual screening study to discover novel PR ligands for therapeutic application.

Leveraging allostery to improve G protein-coupled receptor (GPCR)-directed therapeutics: cannabinoid receptor 1 as discovery target

INTRODUCTION

Allosteric modulators of G-protein coupled receptors (GPCRs) hold the promise of improved pharmacology and safety over typical orthosteric GPCR ligands. These features are particularly relevant to the cannabinoid receptor 1 (CB1R) GPCR, since typical orthosteric CB1R ligands are associated with adverse events that limit their translational potential. Areas covered: The contextual basis for applying allostery to CB1R is considered from pharmacological, drug-discovery, and medicinal standpoints. Rational design of small-molecule CB1R allosteric modulators as potential pharmacotherapeutics would be greatly facilitated by direct experimental characterization of structure-function correlates underlying the biological activity of chemically-diverse CB1R allosteric modulators, CB1R allosteric ligand-binding binding pockets, and amino acid contact residues critical to allosteric ligand engagement and activity. In these regards, designer covalent probes exhibiting well-characterized molecular pharmacology as CB1R allosteric modulators are emerging as valuable molecular reporters enabling experimental interrogation of CB1R allosteric site(s) and informing the design of new CB1R agents as drugs. Expert opinion: Synthesis and pharmacological profiling of CB1R allosteric ligands will continue to provide valuable insights into CB1R structure-function correlates. The resulting data should expand the repertoire of novel agents capable of exerting therapeutic benefit by modulating CB1R-dependent signaling.

Accelerated structure-based design of chemically diverse allosteric modulators of a muscarinic G protein-coupled receptor

Design of ligands that provide receptor selectivity has emerged as a new paradigm for drug discovery of G protein-coupled receptors, and may, for certain families of receptors, only be achieved via identification of chemically diverse allosteric modulators. Here, the extracellular vestibule of the M2 muscarinic acetylcholine receptor (mAChR) is targeted for structure-based design of allosteric modulators. Accelerated molecular dynamics (aMD) simulations were performed to construct structural ensembles that account for the receptor flexibility. Compounds obtained from the National Cancer Institute (NCI) were docked to the receptor ensembles. Retrospective docking of known ligands showed that combining aMD simulations with Glide induced fit docking (IFD) provided much-improved enrichment factors, compared with the Glide virtual screening workflow. Glide IFD was thus applied in receptor ensemble docking, and 38 top-ranked NCI compounds were selected for experimental testing. In [(3)H]N-methylscopolamine radioligand dissociation assays, approximately half of the 38 lead compounds altered the radioligand dissociation rate, a hallmark of allosteric behavior. In further competition binding experiments, we identified 12 compounds with affinity of ≤30 μM. With final functional experiments on six selected compounds, we confirmed four of them as new negative allosteric modulators (NAMs) and one as positive allosteric modulator of agonist-mediated response at the M2 mAChR. Two of the NAMs showed subtype selectivity without significant effect at the M1 and M3 mAChRs. This study demonstrates an unprecedented successful structure-based approach to identify chemically diverse and selective GPCR allosteric modulators with outstanding potential for further structure-activity relationship studies.

Allosteric communication pipelines in G-protein-coupled receptors

The binding of ligands to G-protein-coupled receptors (GPCRs) in the extracellular region transmits the signal to the intracellular region to initiate coupling to effector proteins. The mechanism of this allosteric communication remains largely unexplored. Knowledge of the residues involved in the pipeline of the allosteric communication from the extracellular to the intracellular region will provide means to (a) design ligands with bias in potency towards one signaling pathway over others, and (b) design allosteric modulators that show subtype selectivity in GPCRs. In this review we describe the current state of the computational methods that provide insights into the allosteric communication in GPCRs and elucidate how this information can be used to design allosteric modulators.

G-protein coupled receptors: advances in simulation and drug discovery

G-protein coupled receptors (GPCRs), the largest family of human membrane proteins, mediate cellular signaling and represent primary targets of about one third of currently marketed drugs. GPCRs undergo highly dynamic structural transitions during signal transduction, from binding of extracellular ligands to coupling with intracellular effector proteins. Molecular dynamics (MD) simulations have been utilized to investigate GPCR signaling mechanisms (such as pathways of ligand binding and receptor activation/deactivation) and to design novel small-molecule drug candidates. Future research directions point towards modeling cooperative binding of multiple orthosteric and allosteric ligands to GPCRs, GPCR oligomerization and interactions of GPCRs with different intracellular signaling proteins. Through methodological and supercomputing advances, MD simulations will continue to provide important insights into GPCR signaling mechanisms and further facilitate structure-based drug design.

Activation and conformational dynamics of a class B G-protein-coupled glucagon receptor

The human glucagon receptor (GCGR) is a class B G-protein-coupled receptor (GPCR). The GCGR can be activated by glucagon and regulates the release of glucose. The GCGR has been proposed to be an important drug target for type 2 diabetes. Based on the structural model of a full-length glucagon-bound GCGR (glu-GCGR), we performed accelerated molecular dynamics (aMD) simulations, potential of mean force (PMF) calculations, cross-correlation analysis and community network analysis to study the activation mechanism and the conformational dynamics during the activation process. The PMF map depicts three different conformational states of the GCGR: the inactive, intermediate and active states. The activation of the GCGR is characterized by the outward movement of the intracellular side of helix VI. In the active state of the GCGR, the Arg173(2.46)-Ser350(6.41) and Glu245(3.50)-Thr351(6.42) hydrogen bonds break, and the χ1 rotamer of Phe322(5.54) changes from perpendicular to parallel to helix VI. The binding of the agonist glucagon decreases the correlated motions of the extracellular loops (ELCs) and the helices around the glucagon-binding site. During the activation of the GCGR, the connections between the intracellular sides of helices become weaker, and the connections between glucagon and ECLs and the extracellular sides of helices become stronger. These facilitate G-protein coupling on the intracellular side and glucagon binding on the extracellular side, and stabilize the GCGR in the active state. We expect that this study can provide useful information on the activation mechanism of the GCGR and facilitate the future design of GCGR inhibitors.

Geomfinder: a multi-feature identifier of similar three-dimensional protein patterns: a ligand-independent approach

Background

Since the structure of proteins is more conserved than the sequence, the identification of conserved three-dimensional (3D) patterns among a set of proteins, can be important for protein function prediction, protein clustering, drug discovery and the establishment of evolutionary relationships. Thus, several computational applications to identify, describe and compare 3D patterns (or motifs) have been developed. Often, these tools consider a 3D pattern as that described by the residues surrounding co-crystallized/docked ligands available from X-ray crystal structures or homology models. Nevertheless, many of the protein structures stored in public databases do not provide information about the location and characteristics of ligand binding sites and/or other important 3D patterns such as allosteric sites, enzyme-cofactor interaction motifs, etc. This makes necessary the development of new ligand-independent methods to search and compare 3D patterns in all available protein structures.

Results

Here we introduce Geomfinder, an intuitive, flexible, alignment-free and ligand-independent web server for detailed estimation of similarities between all pairs of 3D patterns detected in any two given protein structures. We used around 1100 protein structures to form pairs of proteins which were assessed with Geomfinder. In these analyses each protein was considered in only one pair (e.g. in a subset of 100 different proteins, 50 pairs of proteins can be defined). Thus: (a) Geomfinder detected identical pairs of 3D patterns in a series of monoamine oxidase-B structures, which corresponded to the effectively similar ligand binding sites at these proteins; (b) we identified structural similarities among pairs of protein structures which are targets of compounds such as acarbose, benzamidine, adenosine triphosphate and pyridoxal phosphate; these similar 3D patterns are not detected using sequence-based methods; (c) the detailed evaluation of three specific cases showed the versatility of Geomfinder, which was able to discriminate between similar and different 3D patterns related to binding sites of common substrates in a range of diverse proteins.

Conclusions

Geomfinder allows detecting similar 3D patterns between any two pair of protein structures, regardless of the divergency among their amino acids sequences. Although the software is not intended for simultaneous multiple comparisons in a large number of proteins, it can be particularly useful in cases such as the structure-based design of multitarget drugs, where a detailed analysis of 3D patterns similarities between a few selected protein targets is essential.

Electronic supplementary material

The online version of this article (doi:10.1186/s13321-016-0131-9) contains supplementary material, which is available to authorized users.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

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.

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation

Accelerated molecular dynamics (aMD) simulations greatly improve the efficiency of conventional molecular dynamics (cMD) for sampling biomolecular conformations, but they require proper reweighting for free energy calculation. In this work, we systematically compare the accuracy of different reweighting algorithms including the exponential average, Maclaurin series, and cumulant expansion on three model systems: alanine dipeptide, chignolin, and Trp-cage. Exponential average reweighting can recover the original free energy profiles easily only when the distribution of the boost potential is narrow (e.g., the range ≤20k_BT) as found in dihedral-boost aMD simulation of alanine dipeptide. In dual-boost aMD simulations of the studied systems, exponential average generally leads to high energetic fluctuations, largely due to the fact that the Boltzmann reweighting factors are dominated by a very few high boost potential frames. In comparison, reweighting based on Maclaurin series expansion (equivalent to cumulant expansion on the first order) greatly suppresses the energetic noise but often gives incorrect energy minimum positions and significant errors at the energy barriers (∼2–3k_BT). Finally, reweighting using cumulant expansion to the second order is able to recover the most accurate free energy profiles within statistical errors of ∼k_BT, particularly when the distribution of the boost potential exhibits low anharmonicity (i.e., near-Gaussian distribution), and should be of wide applicability. A toolkit of Python scripts for aMD reweighting “PyReweighting” is distributed free of charge at http://mccammon.ucsd.edu/computing/amdReweighting/.