Investigation of ECD conformational transition mechanism of GLP-1R by molecular dynamics simulations and Markov state model

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.

Unraveling the Interplay of Extracellular Domain Conformational Changes and Parathyroid Hormone Type 1 Receptor Activation in Class B1 G Protein-Coupled Receptors: Integrating Enhanced Sampling Molecular Dynamics Simulations and Markov State Models

Parathyroid hormone (PTH) type 1 receptor (PTH1R), as a typical class B1 G protein-coupled receptor (GPCR), is responsible for regulating bone turnover and maintaining calcium homeostasis, and its dysregulation has been implicated in the development of several diseases. The extracellular domain (ECD) of PTH1R is crucial for the recognition and binding of ligands, and the receptor may exhibit an autoinhibited state with the closure of the ECD in the absence of ligands. However, the correlation between ECD conformations and PTH1R activation remains unclear. Thus, this study combines enhanced sampling molecular dynamics (MD) simulations and Markov state models (MSMs) to reveal the possible relevance between the ECD conformations and the activation of PTH1R. First, 22 intermediate structures are generated from the autoinhibited state to the active state and conducted for 10 independent 200 ns simulations each. Then, the MSM is constructed based on the cumulative 44 μs simulations with six identified microstates. Finally, the potential interplay between ECD conformational changes and PTH1R activation as well as cryptic allosteric pockets in the intermediate states during receptor activation is revealed. Overall, our findings reveal that the activation of PTH1R has a specific correlation with ECD conformational changes and provide essential insights for GPCR biology and developing novel allosteric modulators targeting cryptic sites.

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.

In-silico study of the interactions between acylated glucagon like-peptide-1 analogues and the native receptor

We have performed a series of multiple molecular dynamics (MD) simulations of glucagon-like peptide-1 (GLP-1) and acylated GLP-1 analogues in complex with the endogenous receptor (GLP-1R) to obtain a molecular understanding of how fatty acid (FA) chain structure, acylation position on the peptide, and presence of a linker affect the binding. MD simulations were analysed to extract heatmaps of receptor-peptide interaction patterns and to determine the free energy of binding using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach. The extracted free energies from MM-PBSA calculations are in qualitative agreement with experimentally determined potencies. Furthermore, the interaction patterns seen in the receptor-GLP-1 complex simulations resemble previously reported binding interactions validating the simulations. Analysing the receptor-GLP-1 analogue complex simulations, we found that the major differences between the systems stem from FA interactions and positioning of acylation in the peptide. Hydrophobic interactions between the FA chain and a hydrophobic patch on the extracellular domain contribute significantly to the binding affinity. Acylation on Lys26 resulted in noticeably more interactions between the FA chain and the extracellular domain hydrophobic patch than found for acylation on Lys34 and Lys38, respectively. The presence of a charged linker between the peptide and FA chain can potentially stabilise the complex by forming hydrogen bonds to arginine residues in the linker region between the extracellular domain and the transmembrane domain. A molecular understanding of the fatty acid structure and its effect on binding provides important insights into designing acylated agonists for GLP-1R.Communicated by Ramaswamy H. Sarma.

Molecular insights into the allosteric coupling mechanism between an agonist and two different transducers for μ-opioid receptors

G protein-coupled receptors (GPCRs) as the most important class of pharmacological targets regulate G-protein and β-arrestin-mediated signaling through allosteric interplay, which are responsible for different biochemical and physiological actions like therapeutic efficacy and side effects. However, the allosteric mechanism underlying preferentially recruiting one transducer versus the other has been poorly understood, limiting drug design. Motivated by this issue, we utilize accelerated molecular dynamics simulation coupled with potential of mean force (PMF), molecular mechanics Poisson Boltzmann surface area (MM/PBSA) and protein structure network (PSN) to study two ternary complex systems of a representative class A GPCR (μ-opioid receptor (μOR)) bound by an agonist and one specific transducer (G-protein and β-arrestin). The results show that no significant difference exists in the whole structure of μOR between two transducer couplings, but displays transducer-dependent changes in the intracellular binding region of μOR, where the β-arrestin coupling results in a narrower crevice with TM7 inward movement compared with the G-protein. In addition, both the G-protein and β-arrestin coupling can increase the binding affinity of the agonist to the receptor. However, the interactions between the agonist and μOR also exhibit transducer-specific changes, in particular for the interaction with ECL2 that plays an important role in recruiting β-arrestin. The allosteric network analysis further indicates that Y148^3.33, F152^3.37, F156^3.41, N191^4.49, T160^3.45, Y106^2.42, W293^6.48, F289^6.44, I248^5.54 and Y252^5.58 play important roles in equally activating G-protein and β-arrestin. In contrast, M161^3.46 and R165^3.50 devote important contributions to preferentially recruit G-protein while D164^3.49 and R179^ICL2 are revealed to be important for selectively activating β-arrestin. The observations provide useful information for understanding the biased activation mechanism.

Independent gradient model based on Hirshfeld partition: A new method for visual study of interactions in chemical systems

The powerful independent gradient model (IGM) method has been increasingly popular in visual analysis of intramolecular and intermolecular interactions in recent years. However, we frequently observed that there is an evident shortcoming of IGM map in graphically studying weak interactions, that is its isosurfaces are usually too bulgy; in these cases, not only the graphical effect is poor, but also the color on some areas on the isosurfaces is inappropriate and may lead to erroneous analysis conclusions. In addition, the IGM method was originally proposed based on promolecular density, which is quite crude and does not take actual electronic structure into account. In this article, we propose an improvement version of IGM, namely IGM based on Hirshfeld partition of molecular density (IGMH), which replaces the free-state atomic densities involved in the IGM method with the atomic densities derived by Hirshfeld partition of actual molecular electron density. This change makes IGM have more rigorous physical background. A large number of application examples in this article, including molecular and periodic systems, weak and chemical bond interactions, fully demonstrate the important value of IGMH in intuitively understanding interactions in chemical systems. Comparisons also showed that the IGMH usually has markedly better graphical effect than IGM and overcomes known problems in IGM. Currently IGMH analysis has been supported in our wavefunction analysis code Multiwfn (http://sobereva.com/multiwfn). We hope that IGMH will become a new useful method among chemists for exploring interactions in wide variety of chemical systems.

Prediction of Protein–Protein Interactions Between Alsin DH/PH and Rac1 and Resulting Protein Dynamics

Alsin is a protein of 1,657 amino acids known for its crucial role in vesicular trafficking in neurons thanks to its ability to interact with two guanosine triphosphatases, Rac1 and Rab5. Evidence suggests that Rac1 can bind Alsin central region, composed by a Dbl Homology (DH) domain followed by a Pleckstrin Homology (PH) domain, leading to Alsin relocalization. However, Alsin three-dimensional structure and its relationship with known biological functions of this protein are still unknown. In this work, a homology model of the Alsin DH/PH domain was developed and studied through molecular dynamics both in the presence and in the absence of its binding partner, Rac1. Due to different conformations of DH domain, the presence of Rac1 seems to stabilize an open state of the protein, while the absence of its binding partner results in closed conformations. Furthermore, Rac1 interaction was able to reduce the fluctuations in the second conserved region of DH motif, which may be involved in the formation of a homodimer. Moreover, the dynamics of DH/PH was described through a Markov State Model to study the pathways linking the open and closed states. In conclusion, this work provided an all-atom model for the DH/PH domain of Alsin protein; moreover, molecular dynamics investigations suggested underlying molecular mechanisms in the signal transduction between Rac1 and Alsin, providing the basis for a deeper understanding of the whole structure–function relationship for Alsin protein.

Acylation of the Incretin Peptide Exendin-4 Directly Impacts Glucagon-Like Peptide-1 Receptor Signaling and Trafficking

The glucagon-like peptide-1 receptor (GLP-1R) is a class B G protein-coupled receptor and mainstay therapeutic target for the treatment of type 2 diabetes and obesity. Recent reports have highlighted how biased agonism at the GLP-1R affects sustained glucose-stimulated insulin secretion through avoidance of desensitization and downregulation. A number of GLP-1R agonists (GLP-1RAs) feature a fatty acid moiety to prolong their pharmacokinetics via increased albumin binding, but the potential for these chemical changes to influence GLP-1R function has rarely been investigated beyond potency assessments for cAMP. Here, we directly compare the prototypical GLP-1RA exendin-4 with its C-terminally acylated analog, exendin-4-C16. We examine relative propensities of each ligand to recruit and activate G proteins and β-arrestins, endocytic and postendocytic trafficking profiles, and interactions with model and cellular membranes in HEK293 and HEK293T cells. Both ligands had similar cAMP potency, but exendin-4-C16 showed ∼2.5-fold bias toward G protein recruitment and a ∼60% reduction in β-arrestin-2 recruitment efficacy compared with exendin-4, as well as reduced GLP-1R endocytosis and preferential targeting toward recycling pathways. These effects were associated with reduced movement of the GLP-1R extracellular domain measured using a conformational biosensor approach and a ∼70% increase in insulin secretion in INS-1 832/3 cells. Interactions with plasma membrane lipids were enhanced by the acyl chain. Exendin-4-C16 showed extensive albumin binding and was highly effective for lowering of blood glucose in mice over at least 72 hours. Our study highlights the importance of a broad approach to the evaluation of GLP-1RA pharmacology. SIGNIFICANCE STATEMENT: Acylation is a common strategy to enhance the pharmacokinetics of peptide-based drugs. This work shows how acylation can also affect various other pharmacological parameters, including biased agonism, receptor trafficking, and interactions with the plasma membrane, which may be therapeutically important.

Understanding the Mechanism of Activation/Deactivation of GLP-1R via Accelerated Molecular Dynamics Simulation

Glucagon-like peptide-1 receptor (GLP-1R), as a member of the class B G protein-coupled receptors (GPCRs), plays a crucial role in regulating blood glucose level signal recognition through its activation. The conformation changes during the activation pathway are of particular importance for its function. To investigate the activation mechanism of GLP-1R, the crystal structures of active and inactive forms are chosen to perform a total of 2 μs of accelerated molecular dynamics (aMD) simulations and 400ns of conventional molecular dynamics (cMD) simulations. With the aid of structural analysis and potential of mean force (PMF) calculations, we reveal the role of different helices in the activation and deactivation process and obtain the intermediate states during activation and deactivation that are difficult to capture in experiments. Protein structure network (PSN) was utilised to clarify the allosteric communication pathways of activation and deactivation and reveal the mechanisms of its activation and deactivation. The results could advance our understanding of the activation mechanism of GLP-1R and the related drug design.

Differential GLP-1R Binding and Activation by Peptide and Non-peptide Agonists

Peptide drugs targeting class B1 G-protein-coupled receptors (GPCRs) can treat multiple diseases; however, there remains substantial interest in the development of orally delivered non-peptide drugs. Here, we reveal unexpected overlap between signaling and regulation of the glucagon-like peptide-1 (GLP-1) receptor by the non-peptide agonist PF 06882961 and GLP-1 that was not observed for another compound, CHU-128. Compounds from these patent series, including PF 06882961, are currently in clinical trials for treatment of type 2 diabetes. High-resolution cryoelectron microscopy (cryo-EM) structures reveal that the binding sites for PF 06882961 and GLP-1 substantially overlap, whereas CHU-128 adopts a unique binding mode with a more open receptor conformation at the extracellular face. Structural differences involving extensive water-mediated hydrogen bond networks could be correlated to functional data to understand how PF 06882961, but not CHU-128, can closely mimic the pharmacological properties of GLP-1. These findings will facilitate rational structure-based discovery of non-peptide agonists targeting class B GPCRs.

Zinc-mediated conformational preselection mechanism in the allosteric control of DNA binding to the zinc transcriptional regulator (ZitR)

The zinc transcriptional regulator (ZitR) functions as a metalloregulator that fine tunes transcriptional regulation through zinc-dependent DNA binding. However, the molecular mechanism of zinc-driven allosteric control of the DNA binding to ZitR remains elusive. Here, we performed enhanced sampling accelerated molecular dynamics simulations to figure out the mechanism, revealing the role of protein dynamics in the zinc-induced allosteric control of DNA binding to ZitR. The results suggest that zinc-free ZitR samples distinct conformational states, only a handful of which are compatible with DNA binding. Remarkably, zinc binding reduces the conformational plasticity of the DNA-binding domain of ZitR, promoting the population shift in the ZitR conformational ensemble towards the DNA binding-competent conformation. Further co-binding of DNA to the zinc–ZitR complex stabilizes this competent conformation. These findings suggest that ZitR–DNA interactions are allosterically regulated in a zinc-mediated conformational preselection manner, highlighting the importance of conformational dynamics in the regulation of transcription factor family.

Conformation Transition of Intracellular Part of Glucagon Receptor in Complex With Agonist Glucagon by Conventional and Accelerated Molecular Dynamics Simulations

The inactive conformations of glucagon receptor (GCGR) are widely reported by crystal structures that support the precision structure for drug discovery of type 2 diabetes. The previous study shows that the intracellular part is open in the glucagon-bound GCGR (glu-GCGR) and closed in the apo-GCGR by accelerated molecular dynamics (aMD) simulations. However, the crystal structure of GCGR in complex with partial agonist shows that the intracellular part is closed in the inactive conformation. To understand the differences between the studies of aMD simulations and crystal structure, the 2,500 ns conventional molecular dynamics (cMD) simulations are performed on the simulated model of glu-GCGR. The result shows that the transmembrane helices (TMH) 6 of glu-GCGR is outward ~4 Å to drive the intracellular part of glu-GCGR open until ~390 ns cMD simulations. The (TMH) 6 of glu-GCGR becomes closed after ~490 ns cMD simulations, which are consistent with the crystal structure of GCGR in complex with the partial agonist. To further elucidate the activation mechanism of GCGR deeply, the simulated models of glu-GCGR, apo-GCGR, and antagonist-bound GCGR (ant-GCGR) are constructed to perform 10 of parallel 300 ns aMD simulations, respectively. The results show that both of glu-GCGR and apo-GCGR can generate the open conformations of the intracellular part. But the glu-GCGR has the higher percentage of open conformations than apo-GCGR. The ant-GCGR is restricted to generate the open conformations of the intracellular part by antagonist MK-0893. It indicates that the glu-GCGR, apo-GCGR, and ant-GCGR can be distinguished by the aMD simulated method. Free energy landscape shows that the open conformations of the intracellular part of GCGR are in intermediate state. Our results show that aMD simulations enhance the space samplings of open conformations of GCGR via adding extra boost potential. It indicates that the aMD simulations are an effective way for drug discovery of GCGR.

New Insights into Beta-Cell GLP-1 Receptor and cAMP Signaling

Harnessing the translational potential of the GLP-1/GLP-1R system in pancreatic beta cells has led to the development of established GLP-1R-based therapies for the long-term preservation of beta cell function. In this review, we discuss recent advances in the current research on the GLP-1/GLP-1R system in beta cells, including the regulation of signaling by endocytic trafficking as well as the application of concepts such as signal bias, allosteric modulation, dual agonism, polymorphic receptor variants, spatial compartmentalization of cAMP signaling and new downstream signaling targets involved in the control of beta cell function.