Search

Glucagon-like peptide-1 receptor (GLP-1R) agonists (GLP-1RAs) are a highly effective therapy class for type 2 diabetes (T2D) and obesity, yet there are variable patient responses. Variation in the human Glp1r gene leading to altered receptor structure, signal transduction, and function might be directly linked to variable therapeutic responses in patients. A naturally occurring, low-frequency, gain-of-function missense variant, rs10305492 G>A (A316T), protects against T2D and cardiovascular disease. Here we employ CRISPR/Cas9 technology to generate a humanised knock-in mouse model bearing the homozygous Glp1r A316T substitution. Human Glp1r A316T/A316T mice displayed lower fasting blood glucose levels and improved glucose tolerance, as well as increased plasma insulin levels and improved insulin secretion compared to human Glp1r +/+ littermates, even under metabolic stress. They also exhibited alterations in islet cytoarchitecture and β-cell identity under a high-fat, high-sucrose (HFHS) diet. This was however associated with blunted responses to pharmacological GLP-1RAs in vivo . Further investigations in several rodent and human β-cell models demonstrated that the human Glp1r A316T variant exhibits characteristics of constitutive activation but dampened GLP-1RA responses. Our results are further supported by the cryo-EM analysis and molecular dynamics (MD) simulations of the GLP-1R A316T structure, collectively demonstrating that the A316T Glp1r variant governs basal receptor activity and pharmacological responses to GLP-1R-targeting anti-diabetic therapies, highlighting the importance of the molecular characterisation of human Glp1r variants to predict individual therapy responses.

Antisense oligonucleotides (ASOs) offer ground-breaking possibilities for selective pharmacological intervention for any gene product-related disease. Therapeutic ASOs contain extensive chemical modifications that improve stability to enzymatic cleavage and modulate binding affinity relative to natural RNA/DNA. Molecular dynamics (MD) simulation can provide valuable insights into how such modifications affect ASO conformational sampling and target binding. However, force field parameters for chemically modified nucleic acids (NAs) are still underdeveloped. To bridge this gap, we developed parameters to allow simulations of ASOs with the widely applied phosphorothioate (PS) backbone modification, and validated these in extensive all-atom MD simulations of relevant PS-modified NA systems representing B-DNA, RNA, and DNA/RNA hybrid duplex structures. Compared to the corresponding natural NAs, single PS substitutions had marginal effects on the ordered DNA/RNA duplex, whereas substantial effects of phosphorothioation were observed in single-stranded RNA and B-DNA, corroborated by the experimentally derived structure data. We find that PS-modified NAs shift between high and low twist states, which could affect target recognition and protein interactions for phosphorothioated oligonucleotides. Furthermore, conformational sampling was markedly altered in the PS-modified ssRNA system compared to that of the natural oligonucleotide, indicating sequence-dependent effects on conformational preference that may in turn influence duplex formation., (© The Author(s) 2024. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

HIV-1 Vpr is a multifunctional accessory protein consisting of 96 amino acids that play a critical role in viral pathogenesis. Among its diverse range of activities, Vpr can create a cation-selective ion channel within the plasma membrane. However, the oligomeric state of this channel has not yet been elucidated. In this study, we investigated the conformational dynamics of Vpr helices to model the ion channel topology. First, we employed a series of multiscale simulations to investigate the specific structure of monomeric Vpr in a membrane model. During the lipid bilayer self-assembly coarse grain simulation, the C-terminal helix (residues 56-77) effectively formed the transmembrane region, while the N-terminal helix exhibited an amphipathic nature by associating horizontally with a single leaflet. All-atom molecular dynamics (MD) simulations of full-length Vpr inside a phospholipid bilayer show that the C-terminal helix remains very stable inside the bilayer core in a vertical orientation. Subsequently, using the predicted C-terminal helix orientation and conformation, various oligomeric states (ranging from tetramer to heptamer) possibly forming the Vpr ion channel were built and further evaluated. Among these models, the pentameric form exhibited consistent stability in MD simulations and displayed a compatible conformation for a water-assisted ion transport mechanism. This study provides structural insights into the ion channel activity of the Vpr protein and the foundation for developing therapeutics against HIV-1 Vpr-related conditions.

In the decade since the discovery of androglobin, a multi-domain hemoglobin of metazoans associated with ciliogenesis and spermatogenesis, there has been little advance in the knowledge of the biochemical and structural properties of this unusual member of the hemoglobin superfamily. Using a method for aligning remote homologues, coupled with molecular modelling and molecular dynamics, we have identified a novel structural alignment to other hemoglobins. This has led to the first stable recombinant expression and characterization of the circularly permuted globin domain. Exceptional for eukaryotic globins is that a tyrosine takes the place of the highly conserved phenylalanine in the CD1 position, a critical point in stabilizing the heme. A disulfide bond, similar to that found in neuroglobin, forms a closed loop around the heme pocket, taking the place of androglobin's missing CD loop and further supporting the heme pocket structure. Highly unusual in the globin superfamily is that the heme iron binds nitric oxide as a five-coordinate complex similar to other heme proteins that have nitric oxide storage functions. With rapid autoxidation and high nitrite reductase activity, the globin appears to be more tailored toward nitric oxide homeostasis or buffering. The use of our multi-template profile alignment method to yield the first biochemical characterisation of the circularly permuted globin domain of androglobin expands our knowledge of the fundamental functioning of this elusive protein and provides a pathway to better define the link between the biochemical traits of androglobin with proposed physiological functions., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Conventional measurements of fasting and postprandial blood glucose levels investigated in genome-wide association studies (GWAS) cannot capture the effects of DNA variability on 'around the clock' glucoregulatory processes. Here we show that GWAS meta-analysis of glucose measurements under nonstandardized conditions (random glucose (RG)) in 476,326 individuals of diverse ancestries and without diabetes enables locus discovery and innovative pathophysiological observations. We discovered 120 RG loci represented by 150 distinct signals, including 13 with sex-dimorphic effects, two cross-ancestry and seven rare frequency signals. Of these, 44 loci are new for glycemic traits. Regulatory, glycosylation and metagenomic annotations highlight ileum and colon tissues, indicating an underappreciated role of the gastrointestinal tract in controlling blood glucose. Functional follow-up and molecular dynamics simulations of lower frequency coding variants in glucagon-like peptide-1 receptor (GLP1R), a type 2 diabetes treatment target, reveal that optimal selection of GLP-1R agonist therapy will benefit from tailored genetic stratification. We also provide evidence from Mendelian randomization that lung function is modulated by blood glucose and that pulmonary dysfunction is a diabetes complication. Our investigation yields new insights into the biology of glucose regulation, diabetes complications and pathways for treatment stratification., (© 2023. The Author(s).)

Nef is a small accessory protein pivotal in the HIV-1 viral replication cycle. It is a multifunctional protein and its interactions with kinases in host cells have been well characterized through many in vitro and structural studies. Nef forms a homodimer to activate the kinases and subsequently the phosphorylation pathways. The disruption of its homodimerization represents a valuable approach in the search for novel classes of antiretroviral. However, this research avenue is still underdeveloped as just a few Nef inhibitors have been reported so far, with very limited structural information about their mechanism of action. To address this issue, we have employed an in silico structure-based drug design strategy that combines de novo ligand design with molecular docking and extensive molecular dynamics simulations. Since the Nef pocket involved in homodimerization has high lipophilicity, the initial de novo-designed structures displayed poor drug-likeness and solubility. Taking information from the hydration sites within the homodimerization pocket, structural modifications in the initial lead compound have been introduced to improve the solubility and drug-likeness, without affecting the binding profile. We propose lead compounds that can be the starting point for further optimizations to deliver long-awaited, rationally designed Nef inhibitors., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

The term cancer refers to a plethora of diseases characterized by the development of abnormal cells that divide uncontrollably and can infiltrate further proximal or distal body tissues. Each type of cancer can be defined by aggressiveness, localization, metabolism, and response to available treatments. Among the most common hallmarks of cancer is a more acidic intracellular microenvironment. Offset pH values are due to an excess of lactate and an increased hypoxia-inducible factor (HIF) expression, which leads to a hypoxic state and a metabolic shift towards glycolysis to produce adenosine-5'-triphosphate (ATP) necessary for cellular metabolism. Warburg's hypothesis underpins this concept, making glycolysis and its central enzyme pyruvate kinase (hPKM2), an ideal target for drug development. Using molecular docking and extensive molecular dynamics (MD) simulations we investigated the binding mode of phosphoenolpyruvate (PEP) inside the hPKM2 active site, and then evaluated a set of known bio-isosteric inhibitors to understand the differences caused by their substitutions on their binding mode. Ultimately, we propose a new molecular entity to hamper hPKM2, unbalance cellular energy, and possibly trigger autophagic mechanisms., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

The spike protein is key to SARS-CoV-2 high infectivity because it facilitates the receptor binding domain (RBD) encounter with ACE2. As targeting subunit S1 has not yet delivered an ACE2-binding inhibitor, we have assessed the druggability of the conserved segment of the spike protein stalk within subunit S2 by means of an integrated computational approach that combines the molecular docking of an optimized library of fragments with high-throughput molecular dynamics simulations. The high propensity of the spike protein to mutate in key regions that are responsible for the recognition of the human angiotensin-converting enzyme 2 (hACE2) or for the recognition of antibodies, has made subunit S1 of the spike protein difficult to target. Despite the inherent flexibility of the stalk region, our results suggest two hidden interhelical binding sites, whose accessibility is only partially hampered by glycan residues.

The vasoactive intestinal peptide (VIP) and pituitary adenylate cyclase-activating polypeptide (PACAP) receptors are key regulators of neurological processes. Despite recent structural data, a comprehensive understanding of peptide binding and selectivity among different subfamily receptors is lacking. Here, we determine structures of active, Gs-coupled, VIP-VPAC1R, PACAP27-VPAC1R, and PACAP27-PAC1R complexes. Cryo-EM structural analyses and molecular dynamics simulations (MDSs) reveal fewer stable interactions between VPAC1R and VIP than for PACAP27, more extensive dynamics of VIP interaction with extracellular loop 3, and receptor-dependent differences in interactions of conserved N-terminal peptide residues with the receptor core. MD of VIP modelled into PAC1R predicts more transient VIP-PAC1R interactions in the receptor core, compared to VIP-VPAC1R, which may underlie the selectivity of VIP for VPAC1R over PAC1R. Collectively, our work improves molecular understanding of peptide engagement with the PAC1R and VPAC1R that may benefit the development of novel selective agonists., (© 2022. The Author(s).)

A series of benzyloxy and phenoxy derivatives of the adenosine receptor agonists N 6 -cyclopentyl adenosine (CPA) and N 6 -cyclopentyl 5'- N -ethylcarboxamidoadenosine (CP-NECA) were synthesized, and their potency and selectivity were assessed. We observed that the most potent were the compounds with a halogen in the meta position on the aromatic ring of the benzyloxy- or phenoxycyclopentyl substituent. In general, the NECA-based compounds displayed greater A 1 R selectivity than the adenosine-based compounds, with N 6 -2-(3-bromobenzyloxy)cyclopentyl-NECA and N 6 -2-(3-methoxyphenoxy)cyclopentyl-NECA showing ∼1500-fold improved A 1 R selectivity compared to NECA. In addition, we quantified the compounds' affinity and kinetics of binding at both human and rat A 1 R using a NanoBRET binding assay and found that the halogen substituent in the benzyloxy- or phenoxycyclopentyl moiety seems to confer high affinity for the A 1 R. Molecular modeling studies suggested a hydrophobic subpocket as contributing to the A 1 R selectivity displayed. We believe that the identified selective potent A 1 R agonists are valuable tool compounds for adenosine receptor research.

After the SARS-CoV-2 Wuhan variant that gave rise to the pandemic, other variants named Delta, Omicron, and Omicron-2 sequentially became prevalent, with mutations spread around the viral genome, including on the spike (S) protein; in order to understand the resultant in gains in infectivity, we interrogated in silico both the equilibrium binding and the binding pathway of the virus' receptor-binding domain (RBD) to the angiotensin-converting enzyme 2 (ACE2) receptor. We interrogated the molecular recognition between the RBD of different variants and ACE2 through supervised molecular dynamics (SuMD) and classic molecular dynamics (MD) simulations to address the effect of mutations on the possible S protein binding pathways. Our results indicate that compensation between binding pathway efficiency and stability of the complex exists for the Omicron BA.1 receptor binding domain, while Omicron BA.2's mutations putatively improved the dynamic recognition of the ACE2 receptor, suggesting an evolutionary advantage over the previous strains.

The SARS-CoV-2 virus is responsible for the COVID-19 pandemic the world experience since 2019. The protein responsible for the first steps of cell invasion, the spike protein, has probably received the most attention in light of its central role during infection. Computational approaches are among the tools employed by the scientific community in the enormous effort to study this new affliction. One of these methods, namely molecular dynamics (MD), has been used to characterize the function of the spike protein at the atomic level and unveil its structural features from a dynamic perspective. In this review, we focus on these main findings, including spike protein flexibility, rare S protein conformational changes, cryptic epitopes, the role of glycans, drug repurposing, and the effect of spike protein variants., (© 2022 The Authors. BioEssays published by Wiley Periodicals LLC.)

The development of therapeutic agonists for G protein-coupled receptors (GPCRs) is hampered by the propensity of GPCRs to couple to multiple intracellular signalling pathways. This promiscuous coupling leads to numerous downstream cellular effects, some of which are therapeutically undesirable. This is especially the case for adenosine A 1 receptors (A 1 Rs) whose clinical potential is undermined by the sedation and cardiorespiratory depression caused by conventional agonists. We have discovered that the A 1 R-selective agonist, benzyloxy-cyclopentyladenosine (BnOCPA), is a potent and powerful analgesic but does not cause sedation, bradycardia, hypotension or respiratory depression. This unprecedented discrimination between native A 1 Rs arises from BnOCPA's unique and exquisitely selective activation of Gob among the six Gαi/o subtypes, and in the absence of β-arrestin recruitment. BnOCPA thus demonstrates a highly-specific Gα-selective activation of the native A 1 R, sheds new light on GPCR signalling, and reveals new possibilities for the development of novel therapeutics based on the far-reaching concept of selective Gα agonism., (© 2022. The Author(s).)

Recent advances in G-protein-coupled receptor (GPCR) structural elucidation have strengthened previous hypotheses that multidimensional signal propagation mediated by these receptors depends, in part, on their conformational mobility; however, the relationship between receptor function and static structures is inherently uncertain. Here, we examine the contribution of peptide agonist conformational plasticity to activation of the glucagon-like peptide 1 receptor (GLP-1R), an important clinical target. We use variants of the peptides GLP-1 and exendin-4 (Ex4) to explore the interplay between helical propensity near the agonist N terminus and the ability to bind to and activate the receptor. Cryo-EM analysis of a complex involving an Ex4 analog, the GLP-1R and G s heterotrimer revealed two receptor conformers with distinct modes of peptide-receptor engagement. Our functional and structural data, along with molecular dynamics (MD) simulations, suggest that receptor conformational dynamics associated with flexibility of the peptide N-terminal activation domain may be a key determinant of agonist efficacy., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

The glucagon-like peptide-1 receptor (GLP-1R) has broad physiological roles and is a validated target for treatment of metabolic disorders. Despite recent advances in GLP-1R structure elucidation, detailed mechanistic understanding of how different peptides generate profound differences in G protein-mediated signalling is still lacking. Here we combine cryo-electron microscopy, molecular dynamics simulations, receptor mutagenesis and pharmacological assays, to interrogate the mechanism and consequences of GLP-1R binding to four peptide agonists; glucagon-like peptide-1, oxyntomodulin, exendin-4 and exendin-P5. These data reveal that distinctions in peptide N-terminal interactions and dynamics with the GLP-1R transmembrane domain are reciprocally associated with differences in the allosteric coupling to G proteins. In particular, transient interactions with residues at the base of the binding cavity correlate with enhanced kinetics for G protein activation, providing a rationale for differences in G protein-mediated signalling efficacy from distinct agonists., (© 2022. The Author(s).)

Glucose-dependent insulinotropic polypeptide (GIP) and its receptor (GIPR) are part of the incretin system that regulates glucose homeostasis. A series of GIPR residues putatively important for ligand binding and receptor activation were mutated and pharmacologically evaluated using GIPR selective agonists in cAMP accumulation, ERK1/2 phosphorylation (pERK1/2) and β-arrestin 2 recruitment assays. The impact of mutation on ligand efficacy was determined by operational modelling of experimental data for each mutant, with results mapped onto the full-length, active-state GIPR structure. Two interaction networks, comprising transmembrane helix (TM) 7, TM1 and TM2, and extracellular loop (ECL) 2, TM5 and ECL3 were revealed, respectively. Both networks were critical for Gα s -mediated cAMP accumulation and the recruitment of β-arrestin 2, however, cAMP response was more sensitive to alanine substitution, with most mutated residues displaying reduced signaling. Unlike the other two assays, activation of ERK1/2 was largely independent of the network involving ECL2, TM5 and ECL3, indicating that pERK1/2 is at least partially distinct from Gα s or β-arrestin pathways and this network is also crucial for potential biased agonism at GIPR. Collectively, our work advances understanding of the structure-function relationship of GIPR and provides a framework for the design and/or interpretation of GIP analogues with unique signaling profiles., (Copyright © 2021. Published by Elsevier Inc.)

Objective: Glucagon-like peptide-1 and glucagon receptor (GLP-1R/GCGR) co-agonism can maximise weight loss and improve glycaemic control in type 2 diabetes and obesity. In this study, we investigated the cellular and metabolic effects of modulating the balance between G protein and β-arrestin-2 recruitment at GLP-1R and GCGR using oxyntomodulin (OXM)-derived co-agonists. This strategy has been previously shown to improve the duration of action of GLP-1R mono-agonists by reducing target desensitisation and downregulation., Methods: Dipeptidyl dipeptidase-4 (DPP-4)-resistant OXM analogues were generated and assessed for a variety of cellular readouts. Molecular dynamic simulations were used to gain insights into the molecular interactions involved. In vivo studies were performed in mice to identify the effects on glucose homeostasis and weight loss., Results: Ligand-specific reductions in β-arrestin-2 recruitment were associated with slower GLP-1R internalisation and prolonged glucose-lowering action in vivo. The putative benefits of GCGR agonism were retained, with equivalent weight loss compared to the GLP-1R mono-agonist liraglutide despite a lesser degree of food intake suppression. The compounds tested showed only a minor degree of biased agonism between G protein and β-arrestin-2 recruitment at both receptors and were best classified as partial agonists for the two pathways measured., Conclusions: Diminishing β-arrestin-2 recruitment may be an effective way to increase the therapeutic efficacy of GLP-1R/GCGR co-agonists. These benefits can be achieved by partial rather than biased agonism., (Copyright © 2021 The Authors. Published by Elsevier GmbH.. All rights reserved.)

Class B1 G protein-coupled receptors (GPCRs) are important targets for many diseases, including cancer, diabetes, and heart disease. All the approved drugs for this receptor family are peptides that mimic the endogenous activating hormones. An understanding of how agonists bind and activate class B1 GPCRs is fundamental for the development of therapeutic small molecules. We combined supervised molecular dynamics (SuMD) and classic molecular dynamics (cMD) simulations to study the binding of the calcitonin gene-related peptide (CGRP) to the CGRP receptor (CGRPR). We also evaluated the association and dissociation of the antagonist telcagepant from the extracellular domain (ECD) of CGRPR and the water network perturbation upon binding. This study, which represents the first example of dynamic docking of a class B1 GPCR peptide, delivers insights on several aspects of ligand binding to CGRPR, expanding understanding of the role of the ECD and the receptor-activity modifying protein 1 (RAMP1) on agonist selectivity., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Deganutti, Atanasio, Rujan, Sexton, Wootten and Reynolds.)

Despite being a target for about one-third of approved drugs, G protein-coupled receptors (GPCRs) still represent a tremendous reservoir for therapeutic strategies against diseases. For example, several cardiovascular and central nervous system conditions could benefit from clinical agents that activate the adenosine 1 receptor (A 1 R); however, the pursuit of A 1 R agonists for clinical use is usually impeded by both on- and off-target side effects. One of the possible strategies to overcome this issue is the development of positive allosteric modulators (PAMs) capable of selectively enhancing the effect of a specific receptor subtype and triggering functional selectivity (a phenomenon also referred to as bias). Intriguingly, besides enforcing the effect of agonists upon binding to an allosteric site, most of the A 1 R PAMs display intrinsic partial agonism and orthosteric competition with antagonists. To rationalize this behavior, we simulated the binding of the prototypical PAMs PD81723 and VCP171, the full-agonist NECA, the antagonist 13B, and the bitopic agonist VCP746. We propose that a single PAM can bind several A 1 R sites rather than a unique allosteric pocket, reconciling the structure-activity relationship and the mutagenesis results.

The recent outbreak of the respiratory syndrome-related coronavirus (SARS-CoV-2) is stimulating an unprecedented scientific campaign to alleviate the burden of the coronavirus disease (COVID-19). One line of research has focused on targeting SARS-CoV-2 proteins fundamental for its replication by repurposing drugs approved for other diseases. The first interaction between the virus and the host cell is mediated by the spike protein on the virus surface and the human angiotensin-converting enzyme (ACE2). Small molecules able to bind the receptor-binding domain (RBD) of the spike protein and disrupt the binding to ACE2 would offer an important tool for slowing, or even preventing, the infection. Here, we screened 2421 approved small molecules in silico and validated the docking outcomes through extensive molecular dynamics simulations. Out of six drugs characterized as putative RBD binders, the cephalosporin antibiotic cefsulodin was further assessed for its effect on the binding between the RBD and ACE2, suggesting that it is important to consider the dynamic formation of the heterodimer between RBD and ACE2 when judging any potential candidate.

Despite being among the most characterized G protein-coupled receptors (GPCRs), adenosine receptors (ARs) have always been a difficult target in drug design. To date, no agonist other than the natural effector and the diagnostic regadenoson has been approved for human use. Recently, the structure of the adenosine A1 receptor (A 1 R) was determined in the active, G i protein complexed state; this has important repercussions for structure-based drug design. Here, we employed supervised molecular dynamics simulations and mutagenesis experiments to extend the structural knowledge of the binding of selective agonists to A 1 R. Our results identify new residues involved in the association and dissociation pathway, they suggest the binding mode of N 6 -cyclopentyladenosine (CPA) related ligands, and they highlight the dramatic effect that chemical modifications can have on the overall binding mechanism, paving the way for the rational development of a structure-kinetics relationship of A 1 R agonists., Competing Interests: The authors declare no competing financial interest., (© 2021 American Chemical Society.)

The free fatty acid receptor 1 (FFAR1, formerly GPR40), is a potential G protein-coupled receptor (GPCR) target for the treatment of type 2 diabetes mellitus (T2DM), as it enhances glucose-dependent insulin secretion upon activation by endogenous long-chain free fatty acids. The presence of two allosterically communicating binding sites and the lack of the conserved GPCR structural motifs challenge the general knowledge of its activation mechanism. To date, four X-ray crystal structures are available for computer-aided drug design. In this study, we employed molecular dynamics (MD) and supervised molecular dynamics (SuMD) to deliver insights into the (un)binding mechanism of the agonist MK-8666, and the allosteric communications between the two experimentally determined FFAR1 binding sites. We found that FFAR1 extracellular loop 2 (ECL2) mediates the binding of the partial agonist MK-8666. Moreover, simulations showed that the agonists MK-8666 and AP8 are reciprocally stabilized and that AP8 influences MK-8666 unbinding from FFAR1.

The class B secretin GPCR (SecR) has broad physiological effects, with target potential for treatment of metabolic and cardiovascular disease. Molecular understanding of SecR binding and activation is important for its therapeutic exploitation. We combined cryo-electron microscopy, molecular dynamics, and biochemical cross-linking to determine a 2.3 Å structure, and interrogate dynamics, of secretin bound to the SecR:Gs complex. SecR exhibited a unique organization of its extracellular domain (ECD) relative to its 7-transmembrane (TM) core, forming more extended interactions than other family members. Numerous polar interactions formed between secretin and the receptor extracellular loops (ECLs) and TM helices. Cysteine-cross-linking, cryo-electron microscopy multivariate analysis and molecular dynamics simulations revealed that interactions between peptide and receptor were dynamic, and suggested a model for initial peptide engagement where early interactions between the far N-terminus of the peptide and SecR ECL2 likely occur following initial binding of the peptide C-terminus to the ECD.

Adenosine receptors (ARs), like many otherGprotein-coupledreceptors (GPCRs), are targets of primary interest indrug design. However, one of the main limits for the development of drugs for this class of GPCRs is the complex selectivity profile usually displayed by ligands. Numerous efforts have been madefor clarifying the selectivity of ARs, leading to the development of many ligand-based models. The structure of the AR subtype A 1 (A 1 AR) has been recently solved,providing important structural insights. In the present work, we rationalized the selectivity profile of two selective A 1 AR and A 2A AR antagonists, investigating their recognition trajectories obtained by Supervised Molecular Dynamics from an unbound state and monitoring the role of the water molecules in the binding site., Competing Interests: The authors declare no conflict of interest.

The recent paradigm shift toward the use of the kinetics parameters in place of thermodynamic constants is leading the computational chemistry community to develop methods for studying the mechanisms of drug binding and unbinding. From this standpoint, molecular dynamics (MD) plays an important role in delivering insight at the molecular scale. However, a known limitation of MD is that the time scales are usually far from those involved in ligand-receptor unbinding events. Here, we show that the algorithm behind supervised MD (SuMD) can simulate the dissociation mechanism of druglike small molecules while avoiding the input of any energy bias to facilitate the transition. SuMD was tested on seven different intermolecular complexes, covering four G protein-coupled receptors: the A 2A and A 1 adenosine receptors, the orexin 2 and the muscarinic 2 receptors, and the soluble globular enzyme epoxide hydrolase. SuMD well-described the multistep nature of ligand-receptor dissociation, rationalized previous experimental data and produced valuable working hypotheses for structure-kinetics relationships.

Adrenomedullin (AM) and calcitonin gene-related peptide (CGRP) receptors are critically important for metabolism, vascular tone, and inflammatory response. AM receptors are also required for normal lymphatic and blood vascular development and angiogenesis. They play a pivotal role in embryo implantation and fertility and can provide protection against hypoxic and oxidative stress. CGRP and AM receptors are heterodimers of the calcitonin receptor-like receptor (CLR) and receptor activity-modifying protein 1 (RAMP1) (CGRPR), as well as RAMP2 or RAMP3 (AM 1 R and AM 2 R, respectively). However, the mechanistic basis for RAMP modulation of CLR phenotype is unclear. In this study, we report the cryo-EM structure of the AM 1 R in complex with AM and Gs at a global resolution of 3.0 Å, and structures of the AM 2 R in complex with either AM or intermedin/adrenomedullin 2 (AM2) and Gs at 2.4 and 2.3 Å, respectively. The structures reveal distinctions in the primary orientation of the extracellular domains (ECDs) relative to the receptor core and distinct positioning of extracellular loop 3 (ECL3) that are receptor-dependent. Analysis of dynamic data present in the cryo-EM micrographs revealed additional distinctions in the extent of mobility of the ECDs. Chimeric exchange of the linker region of the RAMPs connecting the TM helix and the ECD supports a role for this segment in controlling receptor phenotype. Moreover, a subset of the motions of the ECD appeared coordinated with motions of the G protein relative to the receptor core, suggesting that receptor ECD dynamics could influence G protein interactions. This work provides fundamental advances in our understanding of GPCR function and how this can be allosterically modulated by accessory proteins., Competing Interests: The authors declare no competing financial interest., (Copyright © 2020 American Chemical Society.)

A 2A and A 2B adenosine receptors (ARs) are closely related G protein-coupled receptor subtypes, which represent important (potential) drug targets. Despite their almost identical binding sites for adenosine, A 2A ARs are activated by low (nanomolar) adenosine concentrations, while A 2B ARs require micromolar concentrations. In the present study, we exchanged the extracellular loop 2 (ECL2) of the human A 2A AR for that of the A 2B AR. The resulting chimeric A 2A (ECL2-A 2B )AR was investigated in radioligand binding and cAMP accumulation assays in comparison to the wildtype A 2A AR. While the ribose-modified adenosine analog N-ethylcarboxamidoadenosine (NECA) and its 2-substituted derivative CGS-21680 did not exhibit significant changes, adenosine showed dramatically reduced potency and affinity for the A 2A (ECL2-A 2B )AR mutant displaying similarly low potency as for the wt A 2B AR. Supervised molecular dynamics simulation studies predicted a meta-binding site with high affinity for adenosine, but not for NECA, which may contribute to the observed effects., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Class B G-protein-coupled receptors are major targets for the treatment of chronic diseases, including diabetes and obesity 1 . Structures of active receptors reveal peptide agonists engage deep within the receptor core, leading to an outward movement of extracellular loop 3 and the tops of transmembrane helices 6 and 7, an inward movement of transmembrane helix 1, reorganization of extracellular loop 2 and outward movement of the intracellular side of transmembrane helix 6, resulting in G-protein interaction and activation 2-6 . Here we solved the structure of a non-peptide agonist, TT-OAD2, bound to the glucagon-like peptide-1 (GLP-1) receptor. Our structure identified an unpredicted non-peptide agonist-binding pocket in which reorganization of extracellular loop 3 and transmembrane helices 6 and 7 manifests independently of direct ligand interaction within the deep transmembrane domain pocket. TT-OAD2 exhibits biased agonism, and kinetics of G-protein activation and signalling that are distinct from peptide agonists. Within the structure, TT-OAD2 protrudes beyond the receptor core to interact with the lipid or detergent, providing an explanation for the distinct activation kinetics that may contribute to the clinical efficacy of this compound series. This work alters our understanding of the events that drive the activation of class B receptors.

One of the most intriguing findings highlighted from G protein-coupled receptor (GPCR) crystallography is the presence, in many members of class A, of a partially hydrated sodium ion in the middle of the seven transmembrane helices (7TM) bundle. In particular, the human adenosine A 2A receptor (A 2A AR) is the first GPCR in which a monovalent sodium ion was crystallized in a distal site from the canonical orthosteric one, corroborating, from a structural point of view, its role as a negative allosteric modulator. However, the molecular mechanism by which the sodium ion influences the recognition of the A 2A AR agonists is not yet fully understood. In this study, the supervised molecular dynamics (SuMD) technique was exploited to analyse the sodium ion recognition mechanism and how its presence influences the binding of the endogenous agonist adenosine. Due to a higher degree of flexibility of the receptor extracellular (EC) vestibule, we propose the sodium-bound A 2A AR as less efficient in stabilizing the adenosine during the different steps of binding.

Amylin is coexpressed with insulin in pancreatic islet β-cells and has potent effects on gastric emptying and food intake. The effect of amylin on satiation has been postulated to involve AMY 3 receptors (AMY 3 R) that are heteromers of the calcitonin receptor (CTR) and receptor activity-modifying protein 3 (RAMP3). Understanding the molecular control of signaling through the AMY 3 R is thus important for peptide drug targeting of this receptor. We have previously used alanine scanning mutagenesis to study the contribution of the extracellular surface of the CTR to binding and signaling initiated by calcitonin (CT) and related peptides (Dal Maso, E., et al . (2019) The molecular control of calcitonin receptor signaling. ACS Pharmacol. Transl. Sci. 2 , 31-51). That work revealed ligand- and pathway-specific effects of mutation, with extracellular loops (ECLs) 2 and 3 particularly important in the distinct propagation of signaling mediated by individual peptides. In the current study, we have used equivalent alanine scanning of ECL2 and ECL3 of the CTR in the context of coexpression with RAMP3 to form AMY 3 Rs, to examine functional affinity and efficacy of peptides in cAMP accumulation and extracellular signal-regulated kinase (ERK) phosphorylation (pERK). The effect of mutation was determined on representatives of the three major distinct classes of CT peptide, salmon CT (sCT), human CT (hCT), and porcine CT (pCT), as well as rat amylin (rAmy) or human α-CGRP (calcitonin gene-related peptide, hCGRP) whose potency is enhanced by RAMP interaction. We demonstrate that the dynamic nature of CTR ECL2 and ECL3 in propagation of signaling is fundamentally altered when complexed with RAMP3 to form the AMY 3 R, despite only having predicted direct interactions with ECL2. Moreover, the work shows that the role of these loops in receptor signaling is highly peptide dependent, illustrating that even subtle changes to peptide sequence may change signaling output downstream of the receptor., Competing Interests: The authors declare no competing financial interest., (Copyright © 2019 American Chemical Society.)

Molecular dynamics is a state of the art computational tool for the investigation of biophysics phenomenon at a molecular scale, as it enables the modeling of dynamic processes, such as conformational motions, molecular solvation and ligand binding. The recent advances in structural biology have led to a bloom in published G-protein-coupled receptor structures, representing a solid and valuable resource for molecular dynamics studies. During the last decade, indeed, a plethora of physiological and pharmacological facets of this membrane protein superfamily have been addressed by means of molecular dynamics simulations, including the activation mechanism, allosterism and, very recently, biased signaling. Here, we try to recapitulate some of the main contributions that molecular dynamics has recently produced in the field.

The allosteric modulation of G protein-coupled receptors (GPCRs) by sodium ions has received considerable attention as crystal structures of several receptors, in their inactive conformation, show a Na + ion bound to specific residues which, in the human A 2A adenosine receptor (hA 2A AR), are Ser91 3.39 , Trp246 6.48 , Asn280 7.45 , and Asn284 7.49 . A cluster of water molecules completes the coordination of the sodium ion in the putative allosteric site. It is absolutely consolidated that the progress made in the field of GPCRs structural determination has increased the adoption of docking-driven approaches for the identification or the optimization of novel potent and selective ligands. Despite the extensive use of docking protocols in virtual screening approaches, to date, almost any of these studies have been carried out without taking into account the presence of the sodium cation and its first solvation shell in the putative allosteric binding site. In this study, we have focused our attention on determining how the presence of sodium ion binding and additionally its first hydration sphere, in hA 2A AR could influence the ligand positioning accuracy during molecular docking simulations for most of the available resting and activated hA 2A AR crystal structures, using DockBench as a comparative benchmarking tool and implementing a new correlation coefficient (EM). This work provides indications on the evidence that the posing performance (accuracy and/or precision) of the docking protocols in reproducing the crystallographic poses of different hA 2A AR antagonists is generally increased in the presence of the sodium cation and its first solvation shell, in agreement with experimental observations. Consequently, the inclusion of sodium ion and its first solvation shell should be considered in order to facilitate the selection of new potential ligands in all molecular docking-based virtual screening protocols that aim to find novel GPCRs antagonists and inverse agonists.

Calcitonin gene-related peptide (CGRP) is a widely expressed neuropeptide that has a major role in sensory neurotransmission. The CGRP receptor is a heterodimer of the calcitonin receptor-like receptor (CLR) class B G-protein-coupled receptor and a type 1 transmembrane domain protein, receptor activity-modifying protein 1 (RAMP1). Here we report the structure of the human CGRP receptor in complex with CGRP and the G s -protein heterotrimer at 3.3 Å global resolution, determined by Volta phase-plate cryo-electron microscopy. The receptor activity-modifying protein transmembrane domain sits at the interface between transmembrane domains 3, 4 and 5 of CLR, and stabilizes CLR extracellular loop 2. RAMP1 makes only limited direct contact with CGRP, consistent with its function in allosteric modulation of CLR. Molecular dynamics simulations indicate that RAMP1 provides stability to the receptor complex, particularly in the positioning of the extracellular domain of CLR. This work provides insights into the control of G-protein-coupled receptor function.

One of the most largely accepted concepts in the G protein-coupled receptors (GPCRs) field is that the ligand, either agonist or antagonist, recognizes its receptor with a stoichiometry of 1 : 1. Recent experimental evidence, reporting ternary complexes formed by GPCR:orthosteric: allosteric ligands, has complicated the ligand-receptor 1 : 1 binding scenario. Molecular modeling simulations have been used to retrieve insights on the whole ligand-receptor recognition process, beyond information on the final bound state provided by experimental techniques. The simulation of adenosine binding pathways towards the A 2A adenosine receptor highlighted the presence of alternative binding sites (meta-binding sites) beside the canonical orthosteric one, mainly in proximity to the extracellular vestibule. In light of all these considerations, we investigated the possibility that a second molecule of adenosine engages its receptor when this is already in the holo form, generating a ternary complex with a stoichiometry of 2 : 1. Unexpectedly, supervised molecular dynamics (SuMD) simulations showed that the A 2A adenosine receptor could bind the second molecule of adenosine in one of the possible meta-binding sites as well as into its orthosteric site. The formation of this ternary complex, which favored the formation of the intracellular "ionic lock" between R102 (3.50) and E228 (6.30), could putatively be framed in the context of a negative allosteric regulation., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The pancreatic peptide hormone, amylin, plays a critical role in the control of appetite, and synergizes with other key metabolic hormones such as glucagon-like peptide 1 (GLP-1). There is opportunity to develop potent and long-acting analogues of amylin or hybrids between these and GLP-1 mimetics for treating obesity. To achieve this, interrogation of how the 37 amino acid amylin peptide engages with its complex receptor system is required. We synthesized an extensive library of peptides to profile the human amylin sequence, determining the role of its disulfide loop, amidated C -terminus and receptor "capture" and "activation" regions in receptor signaling. We profiled four signaling pathways with different ligands at multiple receptor subtypes, in addition to exploring selectivity determinants between related receptors. Distinct roles for peptide subregions in receptor binding and activation were identified, resulting in peptides with greater activity than the native sequence. Enhanced peptide activity was preserved in the brainstem, the major biological target for amylin. Interpretation of our data using full-length active receptor models supported by molecular dynamics, metadynamics, and supervised molecular dynamics simulations guided the synthesis of a potent dual agonist of GLP-1 and amylin receptors. The data offer new insights into the function of peptide amidation, how allostery drives peptide-receptor interactions, and provide a valuable resource for the development of novel amylin agonists for treating diabetes and obesity., Competing Interests: The authors declare no competing financial interest., (Copyright © 2018 American Chemical Society.)

Class B peptide hormone GPCRs are targets for the treatment of major chronic disease. Peptide ligands of these receptors display biased agonism and this may provide future therapeutic advantage. Recent active structures of the calcitonin (CT) and glucagon-like peptide-1 (GLP-1) receptors reveal distinct engagement of peptides with extracellular loops (ECLs) 2 and 3, and mutagenesis of the GLP-1R has implicated these loops in dynamics of receptor activation. In the current study, we have mutated ECLs 2 and 3 of the human CT receptor (CTR), to interrogate receptor expression, peptide affinity and efficacy. Integration of these data with insights from the CTR and GLP-1R active structures, revealed marked diversity in mechanisms of peptide engagement and receptor activation between the CTR and GLP-1R. While the CTR ECL2 played a key role in conformational propagation linked to Gs/cAMP signalling this was mechanistically distinct from that of GLP-1R ECL2. Moreover, ECL3 was a hotspot for distinct ligand- and pathway-specific effects, and this has implications for the future design of biased agonists of class B GPCRs., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Unquestionably, water appears to be an active player in the noncovalent protein-ligand binding process, as it can either bridge interactions between protein and ligand or can be replaced by the bound ligand. Accordingly, in the last decade, alternative computational methodologies have been sought with the aim of predicting the position and thermodynamic profile of water molecules (i.e., hydration sites) in the binding site using either the ligand-bound or ligand-free protein conformation. Herein, we present an alternative approach, named AquaMMapS, that provides a three-dimensional sampling of putative hydration sites. Interestingly, AquaMMapS can post-inspect molecular dynamics (MD) trajectories obtained from different MD engines using indifferently crystallographic or docking-driven structures as a starting point. Moreover, AquaMMapS is naturally integrated into supervised molecular dynamics (SuMD) simulations, presenting the possibility to inspect hydration sites during the ligand-protein association process. Finally, a penalty scoring method, named AquaMMapScoring(AMS), was developed to evaluate the number and nature of the water molecules displaced by a ligand approaching its binding site during the binding event, guiding a medicinal chemist to explore the most suitable regions of a ligand that can be decorated either with or without interfering with the interaction networks mediated by water molecules with specific recognition regions of the protein., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Ligand-protein binding kinetic rates are growing in importance as parameters to consider in drug discovery and lead optimization. In this study we analysed using surface plasmon resonance (SPR) the transition state (TS) properties of a set of six adenosine A 2A receptor inhibitors, belonging to both the xanthine and the triazolo-triazine scaffolds. SPR highlighted interesting differences among the ligands in the enthalpic and entropic components of the TS energy barriers for the binding and unbinding events. To better understand at a molecular level these differences, we developed suMetaD, a novel molecular dynamics (MD)-based approach combining supervised MD and metadynamics. This method allows simulation of the ligand unbinding and binding events. It also provides the system conformation corresponding to the highest energy barrier the ligand is required to overcome to reach the final state. For the six ligands evaluated in this study their TS thermodynamic properties were linked in particular to the role of water molecules in solvating/desolvating the pocket and the small molecules. suMetaD identified kinetic bottleneck conformations near the bound state position or in the vestibule area. In the first case the barrier is mainly enthalpic, requiring the breaking of strong interactions with the protein. In the vestibule TS location the kinetic bottleneck is instead mainly of entropic nature, linked to the solvent behaviour.

Adenosine deaminase converts adenosine into inosine. In contrast to adenosine, relatively little attention has been paid to the physiological roles of inosine. Nevertheless, recent studies have demonstrated that inosine has neuroprotective, cardioprotective, immunomodulatory, and antidepressive effects. Inosine was recently shown to be a less potent agonist than adenosine at the A 2A adenosine receptor. To better depict the differences in the mechanisms of receptor recognition between adenosine and inosine, we carried out supervised molecular dynamics (SuMD) simulations, and the results are analyzed herein., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Structure-driven fragment-based (SDFB) approaches have provided efficient methods for the identification of novel drug candidates. This strategy has been largely applied in discovering several pharmacological ligand classes, including enzyme inhibitors, receptor antagonists and, more recently, also allosteric (positive and negative) modulators. Recently, Siegal and collaborators reported an interesting study, performed on a detergent-solubilized StaR adenosine A 2A receptor, describing the existence of both fragment-like negative allosteric modulators (NAMs), and fragment-like positive allosteric modulators (PAMs). From this retrospective study, our results suggest that Supervised Molecular Dynamics (SuMD) simulations can support, on a reasonable time scale, the identification of fragment-like PAMs following their receptor recognition pathways and characterizing the possible allosteric binding sites.

Kinetic and thermodynamic ligand-protein binding parameters are gaining growing importance as key information to consider in drug discovery. The determination of the molecular structures, using particularly x-ray and NMR techniques, is crucial for understanding how a ligand recognizes its target in the final binding complex. However, for a better understanding of the recognition processes, experimental studies of ligand-protein interactions are needed. Even though several techniques can be used to investigate both thermodynamic and kinetic profiles for a ligand-protein complex, these procedures are very often laborious, time consuming and expensive. In the last 10 years, computational approaches have enormous potential in providing insights into each of the above effects and in parsing their contributions to the changes in both kinetic and thermodynamic binding parameters. The main purpose of this review is to summarize the state of the art of computational strategies for estimating the kinetic and thermodynamic parameters of a ligand-protein binding.

In this review, we present a survey of the recent advances carried out by our research groups in the field of ligand-GPCRs recognition process simulations recently implemented at the Molecular Modeling Section (MMS) of the University of Padova. We briefly describe a platform of tools we have tuned to aid the identification of novel GPCRs binders and the better understanding of their binding mechanisms, based on two extensively used computational techniques such as molecular docking and MD simulations. The developed methodologies encompass: (i) the selection of suitable protocols for docking studies, (ii) the exploration of the dynamical evolution of ligand-protein interaction networks, (iii) the detailed investigation of the role of water molecules upon ligand binding, and (iv) a glance at the way the ligand might go through prior reaching the binding site., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Molecular recognition is a crucial issue when aiming to interpret the mechanism of known active substances as well as to develop novel active candidates. Unfortunately, simulating the binding process is still a challenging task because it requires classical MD experiments in a long microsecond time scale that are affordable only with a high-level computational capacity. In order to overcome this limiting factor, we have recently implemented an alternative MD approach, named supervised molecular dynamics (SuMD), and successfully applied it to G protein-coupled receptors (GPCRs). SuMD enables the investigation of ligand-receptor binding events independently from the starting position, chemical structure of the ligand, and also from its receptor binding affinity. In this article, we present an extension of the SuMD application domain including different types of proteins in comparison with GPCRs. In particular, we have deeply analyzed the ligand-protein recognition pathways of six different case studies that we grouped into two different classes: globular and membrane proteins. Moreover, we introduce the SuMD-Analyzer tool that we have specifically implemented to help the user in the analysis of the SuMD trajectories. Finally, we emphasize the limit of the SuMD applicability domain as well as its strengths in analyzing the complexity of ligand-protein recognition pathways.

The search for G protein-coupled receptors (GPCRs) allosteric modulators represents an active research field in medicinal chemistry. Allosteric modulators usually exert their activity only in the presence of the orthosteric ligand by binding to protein sites topographically different from the orthosteric cleft. They therefore offer potentially therapeutic advantages by selectively influencing tissue responses only when the endogenous agonist is present. The prediction of putative allosteric site location, however, is a challenging task. In facts, they are usually located in regions showing more structural variation among the family members. In the present work, we applied the recently developed Supervised Molecular Dynamics (SuMD) methodology to interpret at the molecular level the positive allosteric modulation mediated by LUF6000 toward the human adenosine A3 receptor (hA3 AR). Our data suggest at least two possible mechanisms to explain the experimental data available. This study represent, to the best of our knowledge, the first case reported of an allosteric recognition mechanism depicted by means of molecular dynamics simulations., (Copyright © 2015 Elsevier Ltd. All rights reserved.)