Cryo-EM structure of the dual incretin receptor agonist, peptide-19, in complex with the glucagon-like peptide-1 receptor

Evaluation of Structure Prediction and Molecular Docking Tools for Therapeutic Peptides in Clinical Use and Trials Targeting Coronary Artery Disease

This study evaluates the performance of various structure prediction tools and molecular docking platforms for therapeutic peptides targeting coronary artery disease (CAD). Structure prediction tools, including AlphaFold 3, I-TASSER 5.1, and PEP-FOLD 4, were employed to generate accurate peptide conformations. These methods, ranging from deep-learning-based (AlphaFold) to template-based (I-TASSER 5.1) and fragment-based (PEP-FOLD), were selected for their proven capabilities in predicting reliable structures. Molecular docking was conducted using four platforms (HADDOCK 2.4, HPEPDOCK 2.0, ClusPro 2.0, and HawDock 2.0) to assess binding affinities and interactions. A 100 ns molecular dynamics (MD) simulation was performed to evaluate the stability of the peptide-receptor complexes, along with Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) calculations to determine binding free energies. The results demonstrated that Apelin, a therapeutic peptide, exhibited superior binding affinities and stability across all platforms, making it a promising candidate for CAD therapy. Apelin's interactions with key receptors involved in cardiovascular health were notably stronger and more stable compared to the other peptides tested. These findings underscore the importance of integrating advanced computational tools for peptide design and evaluation, offering valuable insights for future therapeutic applications in CAD. Future work should focus on in vivo validation and combination therapies to fully explore the clinical potential of these therapeutic peptides.

Class B1 GPCRs: insights into multireceptor pharmacology for the treatment of metabolic disease

The secretin-like, class B1 subfamily of seven transmembrane-spanning G protein-coupled receptors (GPCRs) consists of 15 members that coordinate important physiological processes. These receptors bind peptide ligands and use a distinct mechanism of activation that is driven by evolutionarily conserved structural features. For the class B1 receptors, the C-terminus of the cognate ligand is initially recognized by the receptor via an N-terminal extracellular domain that forms a hydrophobic ligand-binding groove. This binding enables the N-terminus of the ligand to engage deep into a large volume, open transmembrane pocket of the receptor. Importantly, the phylogenetic basis of this ligand-receptor activation mechanism has provided opportunities to engineer analogs of several class B1 ligands for therapeutic use. Among the most accepted of these are drugs targeting the glucagon-like peptide-1 (GLP-1) receptor for the treatment of type 2 diabetes and obesity. Recently, multifunctional agonists possessing activity at the GLP-1 receptor and the glucose-dependent insulinotropic polypeptide (GIP) receptor, such as tirzepatide, and others that also contain glucagon receptor activity, have been developed. In this article, we review members of the class B1 GPCR family with focus on receptors for GLP-1, GIP, and glucagon, including their signal transduction and receptor trafficking characteristics. The metabolic importance of these receptors is also highlighted, along with the benefit of polypharmacologic ligands. Furthermore, key structural features and comparative analyses of high-resolution cryogenic electron microscopy structures for these receptors in active-state complexes with either native ligands or multifunctional agonists are provided, supporting the pharmacological basis of such therapeutic agents.

Insights into the mechanism of crotamine and potential targets involved in obesity-related metabolic pathways

Crotamine (Ctm) is a peptide isolated from Crotalus durissus terrificus venom. This molecule has been demonstrated to diminish body weight gain and enhance browning in adipose tissue, glucose tolerance, and insulin sensitivity; hence, it has been postulated as an anti-obesogenic peptide. However, the mechanism to elicit the anti-obesogenic effects has yet to be elucidated. Thus, we investigated the possible interaction of Ctm with receptors involved in obesity-related metabolic pathways through protein-protein docking and molecular dynamics refinement. To test the anti-obesogenic mechanism of Ctm, we selected and retrieved 18 targets involved in obesity-related drug discovery from Protein Data Bank. Then, we performed protein-protein dockings. The best three Ctm-target models were selected and refined by molecular dynamics simulations. Molecular docking demonstrated that Ctm was able to interact with 13 of the 18 targets tested. Having a better docking score with glucagon-like peptide-1 receptor (GLP-1R) (-1430.2 kcal/mol), DPP-IV (dipeptidyl peptidase-IV) (-1781.7 kcal/mol) and α-glucosidase (-1232.3 kcal/mol). These three models were refined by molecular dynamics. Ctm demonstrated a higher affinity for GLP-1R (ΔG: -41.886 ± 2.289 kcal/mol). However, Ctm interaction was more stable with DPP-IV (RMSD: 0.360 ± 0.015 nm, Radius of gyration: 2.781 ± 0.009 nm). Moreover, the number of interactions and the molecular mechanics energies of Ctm residues suggest that the interaction of Ctm with these receptors is mainly mediated by basic-hydrophobic dyads Y1-K2, W31-R32, and W33-R34. Together, all these results allow elucidating a possible molecular mechanism behind the previously described anti-obesogenic effects.

The G Protein-First Mechanism for Activation of the Class B Glucagon-like Peptide 1 Receptor Coupled to N-Terminal Domain-Mediated Conformational Progression

Recently, there has been a great deal of excitement about new glucagon-like peptide 1 receptor (GLP-1R) agonists (e.g., semaglutide and tirzepatide) that have received FDA approval for type 2 diabetes and obesity. Although effective, these drugs come with side effects that limit their use. While research efforts continue to focus intensively on long-lasting, orally administered GLP-1R medications with fewer side effects, a major impediment to developing improved GLP-1R medications is that the mechanism by which an agonist activates GLP-1R to imitate signaling is not known. Here we present and validate the G protein (GP)-first mechanism for the GLP-1R supported by extensive atomistic simulations. We propose that GLP-1R is preactivated through the formation of a GLP-1R-GP precoupled complex at the cell membrane prior to ligand binding. Despite a transmembrane helix 6 (TM6)-bentout conformation characteristic of activated GLP-1R, this precoupled complex remains unactivated until an agonist binds to elicit signaling. Notably, this new hypothesis offers a unified and predictive model for the activities of a series of full and partial agonists, including the peptides ExP5, GLP-1(7-36), and GLP-1(9-36). Most surprisingly, our simulations reveal an N-terminus domain (NTD)-swing/agonist-insertion mechanism wherein the long extracellular NTD of GLP-1R tightly holds the C-terminal half of the peptide agonist and progressively shifts the N-terminal head of the peptide to facilitate insertion into the orthosteric pocket. Our findings provide novel mechanistic insights into the activation and function of class B GPCRs and should provide a realistic basis for structure-based ligand design.

Promiscuous G-protein activation by the calcium-sensing receptor

The human calcium-sensing receptor (CaSR) detects fluctuations in the extracellular Ca2+ concentration and maintains Ca2+ homeostasis1,2. It also mediates diverse cellular processes not associated with Ca2+ balance3-5. The functional pleiotropy of CaSR arises in part from its ability to signal through several G-protein subtypes6. We determined structures of CaSR in complex with G proteins from three different subfamilies: Gq, Gi and Gs. We found that the homodimeric CaSR of each complex couples to a single G protein through a common mode. This involves the C-terminal helix of each Gα subunit binding to a shallow pocket that is formed in one CaSR subunit by all three intracellular loops (ICL1-ICL3), an extended transmembrane helix 3 and an ordered C-terminal region. G-protein binding expands the transmembrane dimer interface, which is further stabilized by phospholipid. The restraint imposed by the receptor dimer, in combination with ICL2, enables G-protein activation by facilitating conformational transition of Gα. We identified a single Gα residue that determines Gq and Gs versus Gi selectivity. The length and flexibility of ICL2 allows CaSR to bind all three Gα subtypes, thereby conferring capacity for promiscuous G-protein coupling.

An imbalanced GLP-1R/GIPR co-agonist peptide with a site-specific N-terminal PEGylation to maximize metabolic benefits

Glycemic and body weight control gained from GLP-1R agonists remains an unmet need for diabetes and obesity treatment, leading to the development of GLP-1R/GIPR co-agonists. An imbalance in GLP-1R/GIPR agonism may extensively maximize the glucose- and weight-lowering effects. Hence, we prepared a potent and imbalanced GLP-1R/GIPR co-agonist, and refined its action time through a site-specific N-terminal PEGylation strategy. The pharmacological efficacy of these resulting long-acting co-agonists was interrogated both in vitro and in vivo. The results showed that peptide 1 possessed potent and imbalanced receptor-stimulating potency favoring GIP activity, but its hypoglycemic action was disrupted probably resulting from its short half-life. After PEGylation to improve the pharmacokinetics, the pharmacological effects were amplified compared to native peptide 1. Among the resulting derivatives, D-5K exhibited significant glycemic, HbA1c, body-weight, and food-intake control, outperforming GLP-1R mono-agonists. Based on its excellent pharmacological profiles, D-5K may hold the great therapeutic potential for diabetes and obesity treatment.

Structure and Molecular Mechanism of Signaling for the Glucagon-like Peptide-1 Receptor Bound to Gs Protein and Exendin-P5 Biased Agonist

The glucagon-like peptide-1 receptor (GLP-1R) is a key regulator of blood glucose and a prime target for the treatment of type II diabetes and obesity with multiple public drugs. Here we present a comprehensive computational analysis of the interactions of the activated GLP-1R-Gs signaling complex with a G protein biased agonist, Exendin P5 (ExP5), which possesses a unique N-terminal sequence responsible for the signal bias. Using a refined all-atom model of the ExP5-GLP-1R-Gs complex in molecular dynamics (MD) simulations, we propose a novel mechanism of conformation transduction in which the unique interaction network of ExP5 N-terminus propagates the binding signal across an array of conserved residues at the transmembrane domain to enhance Gs protein coupling at the cytoplasmic end of the receptor. Our simulations reveal previously unobserved interactions important for activation by ExP5 toward GDP-GTP signaling, providing new insights into the mechanism of class B G protein-coupled receptor (GPCR) signaling. These findings offer a framework for the structure-based design of more effective therapeutics.

Selective polypeptide ligand binding to the extracellular surface of the transmembrane domains of the class B GPCRs GLP-1R and GCGR

Crystal and cryo-EM structures of the glucagon-like peptide-1 receptor (GLP-1R) and glucagon receptor (GCGR) bound with their peptide ligands have been obtained with full-length constructs, indicating that the extracellular domain (ECD) is indispensable for specific ligand binding. This article complements these data with studies of ligand recognition of the two receptors in solution. Paramagnetic NMR relaxation enhancement measurements using dual labeling with fluorine-19 probes on the receptor and nitroxide spin labels on the peptide ligands provided new insights. The glucagon-like peptide-1 (GLP-1) was found to interact with GLP-1R by selective binding to the extracellular surface. The ligand selectivity toward the extracellular surface of the receptor was preserved in the transmembrane domain (TMD) devoid of the ECD. The dual labeling approach further provided evidence of cross-reactivity of GLP-1R and GCGR with glucagon and GLP-1, respectively, which is of interest in the context of medical treatments using combinations of the two polypeptides.

Blood pH Analysis in Combination with Molecular Medical Tools in Relation to COVID-19 Symptoms

The global outbreak of SARS-CoV-2/COVID-19 provided the stage to accumulate an enormous biomedical data set and an opportunity as well as a challenge to test new concepts and strategies to combat the pandemic. New research and molecular medical protocols may be deployed in different scientific fields, e.g., glycobiology, nanopharmacology, or nanomedicine. We correlated clinical biomedical data derived from patients in intensive care units with structural biology and biophysical data from NMR and/or CAMM (computer-aided molecular modeling). Consequently, new diagnostic and therapeutic approaches against SARS-CoV-2 were evaluated. Specifically, we tested the suitability of incretin mimetics with one or two pH-sensitive amino acid residues as potential drugs to prevent or cure long-COVID symptoms. Blood pH values in correlation with temperature alterations in patient bodies were of clinical importance. The effects of biophysical parameters such as temperature and pH value variation in relation to physical-chemical membrane properties (e.g., glycosylation state, affinity of certain amino acid sequences to sialic acids as well as other carbohydrate residues and lipid structures) provided helpful hints in identifying a potential Achilles heel against long COVID. In silico CAMM methods and in vitro NMR experiments (including ³¹P NMR measurements) were applied to analyze the structural behavior of incretin mimetics and SARS-CoV fusion peptides interacting with dodecylphosphocholine (DPC) micelles. These supramolecular complexes were analyzed under physiological conditions by ¹H and ³¹P NMR techniques. We were able to observe characteristic interaction states of incretin mimetics, SARS-CoV fusion peptides and DPC membranes. Novel interaction profiles (indicated, e.g., by ³¹P NMR signal splitting) were detected. Furthermore, we evaluated GM1 gangliosides and sialic acid-coated silica nanoparticles in complex with DPC micelles in order to create a simple virus host cell membrane model. This is a first step in exploring the structure-function relationship between the SARS-CoV-2 spike protein and incretin mimetics with conserved pH-sensitive histidine residues in their carbohydrate recognition domains as found in galectins. The applied methods were effective in identifying peptide sequences as well as certain carbohydrate moieties with the potential to protect the blood-brain barrier (BBB). These clinically relevant observations on low blood pH values in fatal COVID-19 cases open routes for new therapeutic approaches, especially against long-COVID symptoms.

Molecular insights into peptide agonist engagement with the PTH receptor

The parathyroid hormone (PTH) 1 receptor (PTH1R) is a G protein-coupled receptor (GPCR) that regulates skeletal development and calcium homeostasis. Here, we describe cryo-EM structures of the PTH1R in complex with fragments of the two hormones, PTH and PTH-related protein, the drug abaloparatide, as well as the engineered tool compounds, long-acting PTH (LA-PTH) and the truncated peptide, M-PTH(1-14). We found that the critical N terminus of each agonist engages the transmembrane bundle in a topologically similar fashion, reflecting similarities in measures of Gαs activation. The full-length peptides induce subtly different extracellular domain (ECD) orientations relative to the transmembrane domain. In the structure bound to M-PTH, the ECD is unresolved, demonstrating that the ECD is highly dynamic when unconstrained by a peptide. High resolutions enabled identification of water molecules near peptide and G protein binding sites. Our results illuminate the action of orthosteric agonists of the PTH1R.

Structural Understanding of Peptide-Bound G Protein-Coupled Receptors: Peptide-Target Interactions

The activation of G protein-coupled receptors (GPCRs) is triggered by ligand binding to their orthosteric sites, which induces ligand-specific conformational changes. Agonists and antagonists bound to GPCR orthosteric sites provide detailed information on ligand-binding modes. Among these, peptide ligands play an instrumental role in GPCR pharmacology and have attracted increased attention as therapeutic drugs. The recent breakthrough in GPCR structural biology has resulted in the remarkable availability of peptide-bound GPCR complexes. Despite the several structural similarities shared by these receptors, they exhibit distinct features in terms of peptide recognition and receptor activation. From this perspective, we have summarized the current status of peptide-bound GPCR structural complexes, largely focusing on the interactions between the receptor and its peptide ligand at the orthosteric site. In-depth structural investigations have yielded valuable insights into the molecular mechanisms underlying peptide recognition. This study would contribute to the discovery of GPCR peptide drugs with improved therapeutic effects.

New insights into the structure and function of class B1 GPCRs

G protein-coupled receptors (GPCRs) are the largest family of cell surface receptors. Class B1 GPCRs constitute a subfamily of 15 receptors that characteristically contain large extracellular domains (ECDs) and respond to long polypeptide hormones. Class B1 GPCRs are critical regulators of homeostasis, and as such, many are important drug targets. While most transmembrane proteins, including GPCRs, are recalcitrant to crystallization, recent advances in electron cryo-microscopy (cryo-EM) have facilitated a rapid expansion of the structural understanding of membrane proteins. As a testament to this success, structures for all the class B1 receptors bound to G proteins have been determined by cryo-EM in the past five years. Further advances in cryo-EM have uncovered dynamics of these receptors, ligands, and signalling partners. Here, we examine the recent structural underpinnings of the class B1 GPCRs with an emphasis on structure-function relationships.

Understanding VPAC receptor family peptide binding and selectivity

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.

Membranes under the Magnetic Lens: A Dive into the Diverse World of Membrane Protein Structures Using Cryo-EM

Membrane proteins are highly diverse in both structure and function and can, therefore, present different challenges for structure determination. They are biologically important for cells and organisms as gatekeepers for information and molecule transfer across membranes, but each class of membrane proteins can present unique obstacles to structure determination. Historically, many membrane protein structures have been investigated using highly engineered constructs or using larger fusion proteins to improve solubility and/or increase particle size. Other strategies included the deconstruction of the full-length protein to target smaller soluble domains. These manipulations were often required for crystal formation to support X-ray crystallography or to circumvent lower resolution due to high noise and dynamic motions of protein subdomains. However, recent revolutions in membrane protein biochemistry and cryo-electron microscopy now provide an opportunity to solve high resolution structures of both large, >1 megadalton (MDa), and small, <100 kDa (kDa), drug targets in near-native conditions, routinely reaching resolutions around or below 3 Å. This review provides insights into how the recent advances in membrane biology and biochemistry, as well as technical advances in cryo-electron microscopy, help us to solve structures of a large variety of membrane protein groups, from small receptors to large transporters and more complex machineries.

Biased Agonism and Polymorphic Variation at the GLP-1 Receptor: Implications for the Development of Personalised Therapeutics

Glucagon-like peptide-1 receptor (GLP-1R) is a well-studied incretin hormone receptor and target of several therapeutic drugs for type 2 diabetes (T2D), obesity and, more recently, cardiovascular disease. Some signalling pathways downstream of GLP-1R may be responsible for drug adverse effects such as nausea, while others mediate therapeutic outcomes of incretin-based T2D therapeutics. Understanding the interplay between different factors that alter signalling, trafficking, and receptor activity, including biased agonism, single nucleotide polymorphisms and structural modifications is key to develop the next-generation of personalised GLP-1R agonists. However, these interactions remain poorly described, especially for novel therapeutics such as dual and tri-agonists that target more than one incretin receptor. Comparison of GLP-1R structures in complex with G proteins and different peptide and non-peptide agonists has revealed novel insights into important agonist-residue interactions and networks crucial for receptor activation, recruitment of G proteins and engagement of specific signalling pathways. Here, we review the latest knowledge on GLP-1R structure and activation, providing structural evidence for biased agonism and delineating important networks associated with this phenomenon. We survey current biased agonists and multi-agonists at different stages of development, highlighting possible challenges in their translational potential. Lastly, we discuss findings related to non-synonymous genomic variants of GLP1R and the functional importance of specific residues involved in GLP-1R function. We propose that studies of GLP-1R polymorphisms, and specifically their effect on receptor dynamics and pharmacology in response to biased agonists, could have a significant impact in delineating precision medicine approaches and development of novel therapeutics.

Secretin Amino-Terminal Structure-Activity Relationships and Complementary Mutagenesis at the Site of Docking to the Secretin Receptor

Class B1 G protein-coupled receptors are activated by peptides, with amino-terminal regions critical for biologic activity. Although high resolution structures exist, understanding of key features of the peptide activation domain that drive signaling is limited. In the secretin receptor (SecR) structure, interactions are observed between peptide residues His1 and Ser2 and seventh transmembrane segment (TM7) receptor residue E373. We interrogated these interactions using systematic structure-activity analysis of peptide and receptor. His1 was critical for binding and cAMP responses, but its orientation was not critical, and substitution could independently modify affinity and efficacy. Ser2 was also critical, with all substitutions reducing peptide affinity and functional responses proportionally. Mutation of E373 to conserved acidic Asp (E373D), uncharged polar Gln (E373Q), or charge-reversed basic Arg (E373R) did not alter receptor expression, with all exhibiting secretin-dependent cAMP accumulation. All position 373 mutants displayed reduced binding affinities and cAMP potencies for many peptide analogs, although relative effects of position 1 peptides were similar whereas position 2 peptides exhibited substantial differences. The peptide including basic Lys in position 2 was active at SecR having acidic Glu in position 373 and at E373D while exhibiting minimal activity at those receptors in which an acidic residue is absent in this position (E373Q and E373R). In contrast, the peptide including acidic Glu in position 2 was equipotent with secretin at E373R while being much less potent than secretin at wild-type SecR and E373D. These data support functional importance of a charge-charge interaction between the amino-terminal region of secretin and the top of TM7. SIGNIFICANCE STATEMENT: This work refines our molecular understanding of the activation mechanisms of class B1 G protein-coupled receptors. The amino-terminal region of secretin interacts with the seventh transmembrane segment of its receptor with structural specificity and with a charge-charge interaction helping to drive functional activation.

Implications of ligand-receptor binding kinetics on GLP-1R signalling

G protein-coupled receptors (GPCRs) are the largest class of membrane proteins and in recent years there has been a growing appreciation of the importance in understanding temporal aspects of GPCR behaviour, including the kinetics of ligand binding and downstream receptor mediated signalling. Class B1 GPCRs are activated by peptide agonists and are validated therapeutic targets for numerous diseases. However, the kinetics of ligand binding and how this is linked to downstream activation of signalling cascades is not routinely assessed in development of peptide agonists for this receptor class. The glucagon-like peptide-1 receptor (GLP-1R) is a prototypical class B1 GPCR and a validated target for treatment of global health burdens, including type 2 diabetes and obesity. In this study we examined the kinetics of different steps in GLP-1R activation and subsequent cAMP production mediated by a series of GLP-1R peptide agonists, including the ligand-receptor interaction, ligand-receptor-mediated G protein engagement and conformational change and cAMP production. Our results revealed GLP-1R peptide agonist dissociation kinetics (K_off), but not association kinetics (K_on), were positively correlated with the onset of receptor-G protein coupling/conformational change, onset of cAMP production and duration of cAMP signalling. Thus, this study advances the understanding of molecular events that couple GLP-1R ligand binding to intracellular signaling, with the findings likely to have implications for mechanistic understanding of agonist action at other related class B1 GPCRs.