Recent advances in computational studies of GPCR-G protein interactions

Real-Time Determination of Intracellular cAMP Reveals Functional Coupling of Gs Protein to the Melatonin MT1 Receptor

Melatonin is a neuroendocrine hormone that regulates the circadian rhythm and many other physiological processes. Its functions are primarily exerted through two subtypes of human melatonin receptors, termed melatonin type-1 (MT1) and type-2 (MT2) receptors. Both MT1 and MT2 receptors are generally classified as Gi-coupled receptors owing to their well-recognized ability to inhibit cAMP accumulation in cells. However, it remains an enigma as to why melatonin stimulates cAMP production in a number of cell types that express the MT1 receptor. To address if MT1 can dually couple to Gs and Gi proteins, we employed a highly sensitive luminescent biosensor (GloSensorTM) to monitor the real-time changes in the intracellular cAMP level in intact live HEK293 cells that express MT1 and/or MT2. Our results demonstrate that the activation of MT1, but not MT2, leads to a robust enhancement on the forskolin-stimulated cAMP formation. In contrast, the activation of either MT1 or MT2 inhibited cAMP synthesis driven by the activation of the Gs-coupled β2-adrenergic receptor, which is consistent with a typical Gi-mediated response. The co-expression of MT1 with Gs enabled melatonin itself to stimulate cAMP production, indicating a productive coupling between MT1 and Gs. The possible existence of a MT1-Gs complex was supported through molecular modeling as the predicted complex exhibited structural and thermodynamic characteristics that are comparable to that of MT1-Gi. Taken together, our data reveal that MT1, but not MT2, can dually couple to Gs and Gi proteins, thereby enabling the bi-directional regulation of adenylyl cyclase to differentially modulate cAMP levels in cells that express different complements of MT1, MT2, and G proteins.

GPCR-IPL score: multilevel featurization of GPCR-ligand interaction patterns and prediction of ligand functions from selectivity to biased activation

G-protein-coupled receptors (GPCRs) mediate diverse cell signaling cascades after recognizing extracellular ligands. Despite the successful history of known GPCR drugs, a lack of mechanistic insight into GPCR challenges both the deorphanization of some GPCRs and optimization of the structure-activity relationship of their ligands. Notably, replacing a small substituent on a GPCR ligand can significantly alter extracellular GPCR-ligand interaction patterns and motion of transmembrane helices in turn to occur post-binding events of the ligand. In this study, we designed 3D multilevel features to describe the extracellular interaction patterns. Subsequently, these 3D features were utilized to predict the post-binding events that result from conformational dynamics from the extracellular to intracellular areas. To understand the adaptability of GPCR ligands, we collected the conformational information of flexible residues during binding and performed molecular featurization on a broad range of GPCR-ligand complexes. As a result, we developed GPCR-ligand interaction patterns, binding pockets, and ligand features as score (GPCR-IPL score) for predicting the functional selectivity of GPCR ligands (agonism versus antagonism), using the multilevel features of (1) zoomed-out 'residue level' (for flexible transmembrane helices of GPCRs), (2) zoomed-in 'pocket level' (for sophisticated mode of action) and (3) 'atom level' (for the conformational adaptability of GPCR ligands). GPCR-IPL score demonstrated reliable performance, achieving area under the receiver operating characteristic of 0.938 and area under the precision-recall curve of 0.907 (available in gpcr-ipl-score.onrender.com). Furthermore, we used the molecular features to predict the biased activation of downstream signaling (Gi/o, Gq/11, Gs and β-arrestin) as well as the functional selectivity. The resulting models are interpreted and applied to out-of-set validation with three scenarios including the identification of a new MRGPRX antagonist.

A multiscale predictive digital twin for neurocardiac modulation

Cardiac function is tightly regulated by the autonomic nervous system (ANS). Activation of the sympathetic nervous system increases cardiac output by increasing heart rate and stroke volume, while parasympathetic nerve stimulation instantly slows heart rate. Importantly, imbalance in autonomic control of the heart has been implicated in the development of arrhythmias and heart failure. Understanding of the mechanisms and effects of autonomic stimulation is a major challenge because synapses in different regions of the heart result in multiple changes to heart function. For example, nerve synapses on the sinoatrial node (SAN) impact pacemaking, while synapses on contractile cells alter contraction and arrhythmia vulnerability. Here, we present a multiscale neurocardiac modelling and simulator tool that predicts the effect of efferent stimulation of the sympathetic and parasympathetic branches of the ANS on the cardiac SAN and ventricular myocardium. The model includes a layered representation of the ANS and reproduces firing properties measured experimentally. Model parameters are derived from experiments and atomistic simulations. The model is a first prototype of a digital twin that is applied to make predictions across all system scales, from subcellular signalling to pacemaker frequency to tissue level responses. We predict conditions under which autonomic imbalance induces proarrhythmia and can be modified to prevent or inhibit arrhythmia. In summary, the multiscale model constitutes a predictive digital twin framework to test and guide high-throughput prediction of novel neuromodulatory therapy. KEY POINTS: A multi-layered model representation of the autonomic nervous system that includes sympathetic and parasympathetic branches, each with sparse random intralayer connectivity, synaptic dynamics and conductance based integrate-and-fire neurons generates firing patterns in close agreement with experiment. A key feature of the neurocardiac computational model is the connection between the autonomic nervous system and both pacemaker and contractile cells, where modification to pacemaker frequency drives initiation of electrical signals in the contractile cells. We utilized atomic-scale molecular dynamics simulations to predict the association and dissociation rates of noradrenaline with the β-adrenergic receptor. Multiscale predictions demonstrate how autonomic imbalance may increase proclivity to arrhythmias or be used to terminate arrhythmias. The model serves as a first step towards a digital twin for predicting neuromodulation to prevent or reduce disease.

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-β2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and β2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Decoding Partner Specificity of Opioid Receptor Family

This paper describes an exciting big data analysis compiled in a freely available database, which can be applied to characterize the coupling of different G-Protein coupled receptors (GPCRs) families with their intracellular partners. Opioid receptor (OR) family was used as case study in order to gain further insights into the physiological properties of these important drug targets, known to be associated with the opioid crisis, a huge socio-economic issue directly related to drug abuse. An extensive characterization of all members of the ORs family (μ (MOR), δ (DOR), κ (KOR), nociceptin (NOP)) and their corresponding binding partners (ARRs: Arr2, Arr3; G-protein: Gi1, Gi2, Gi3, Go, Gob, Gz, Gq, G11, G14, G15, G12, Gssh, Gslo) was carried out. A multi-step approach including models' construction (multiple sequence alignment, homology modeling), complex assembling (protein complex refinement with HADDOCK and complex equilibration), and protein-protein interface characterization (including both structural and dynamics analysis) were performed. Our database can be easily applied to several GPCR sub-families, to determine the key structural and dynamical determinants involved in GPCR coupling selectivity.

Specific Engineered G Protein Coupling to Histamine Receptors Revealed from Cellular Assay Experiments and Accelerated Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are targets of extracellular stimuli and hence occupy a key position in drug discovery. By specific and not yet fully elucidated coupling profiles with α subunits of distinct G protein families, they regulate cellular responses. The histamine H2 and H4 receptors (H2R and H4R) are prominent members of Gs- and Gi-coupled GPCRs. Nevertheless, promiscuous G protein and selective Gi signaling have been reported for the H2R and H4R, respectively, the molecular mechanism of which remained unclear. Using a combination of cellular experimental assays and Gaussian accelerated molecular dynamics (GaMD) simulations, we investigated the coupling profiles of the H2R and H4R to engineered mini-G proteins (mG). We obtained coupling profiles of the mGs, mGsi, or mGsq proteins to the H2R and H4R from the mini-G protein recruitment assays using HEK293T cells. Compared to H2R-mGs expressing cells, histamine responses were weaker (pEC50, Emax) for H2R-mGsi and -mGsq. By contrast, the H4R selectively bound to mGsi. Similarly, in all-atom GaMD simulations, we observed a preferential binding of H2R to mGs and H4R to mGsi revealed by the structural flexibility and free energy landscapes of the complexes. Although the mG α5 helices were consistently located within the HR binding cavity, alternative binding orientations were detected in the complexes. Due to the specific residue interactions, all mG α5 helices of the H2R complexes adopted the Gs-like orientation toward the receptor transmembrane (TM) 6 domain, whereas in H4R complexes, only mGsi was in the Gi-like orientation toward TM2, which was in agreement with Gs- and Gi-coupled GPCRs structures resolved by X-ray/cryo-EM. These cellular and molecular insights support (patho)physiological profiles of the histamine receptors, especially the hitherto little studied H2R function in the brain, as well as of the pharmacological potential of H4R selective drugs.

Strategies towards Targeting Gαi/s Proteins: Scanning of Protein-Protein Interaction Sites To Overcome Inaccessibility

Heterotrimeric G proteins are classified into four subfamilies and play a key role in signal transduction. They transmit extracellular signals to intracellular effectors subsequent to the activation of G protein-coupled receptors (GPCRs), which are targeted by over 30 % of FDA-approved drugs. However, addressing G proteins as drug targets represents a compelling alternative, for example, when G proteins act independently of the corresponding GPCRs, or in cases of complex multifunctional diseases, when a large number of different GPCRs are involved. In contrast to Gαq, efforts to target Gαi/s by suitable chemical compounds has not been successful so far. Here, a comprehensive analysis was conducted examining the most important interface regions of Gαi/s with its upstream and downstream interaction partners. By assigning the existing compounds and the performed approaches to the respective interfaces, the druggability of the individual interfaces was ranked to provide perspectives for selective targeting of Gαi/s in the future.

Prediction and targeting of GPCR oligomer interfaces

GPCR oligomerization has emerged as a hot topic in the GPCR field in the last years. Receptors that are part of these oligomers can influence each other's function, although it is not yet entirely understood how these interactions work. The existence of such a highly complex network of interactions between GPCRs generates the possibility of alternative targets for new therapeutic approaches. However, challenges still exist in the characterization of these complexes, especially at the interface level. Different experimental approaches, such as FRET or BRET, are usually combined to study GPCR oligomer interactions. Computational methods have been applied as a useful tool for retrieving information from GPCR sequences and the few X-ray-resolved oligomeric structures that are accessible, as well as for predicting new and trustworthy GPCR oligomeric interfaces. Machine-learning (ML) approaches have recently helped with some hindrances of other methods. By joining and evaluating multiple structure-, sequence- and co-evolution-based features on the same algorithm, it is possible to dilute the issues of particular structures and residues that arise from the experimental methodology into all-encompassing algorithms capable of accurately predict GPCR-GPCR interfaces. All these methods used as a single or a combined approach provide useful information about GPCR oligomerization and its role in GPCR function and dynamics. Altogether, we present experimental, computational and machine-learning methods used to study oligomers interfaces, as well as strategies that have been used to target these dynamic complexes.

Mechanistic Insights into Specific G Protein Interactions with Adenosine Receptors

Coupling between G-protein-coupled receptors (GPCRs) and the G proteins is a key step in cellular signaling. Despite extensive experimental and computational studies, the mechanism of specific GPCR-G protein coupling remains poorly understood. This has greatly hindered effective drug design of GPCRs that are primary targets of ∼¹/₃ of currently marketed drugs. Here, we have employed all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method to decipher the mechanism of the GPCR-G protein interactions. Adenosine receptors (ARs) were used as model systems based on very recently determined cryo-EM structures of the A₁AR and A_2AAR coupled with the G_i and G_s proteins, respectively. Changing the G_i protein to the G_s led to increased fluctuations in the A₁AR and agonist adenosine (ADO), while agonist 5'-N-ethylcarboxamidoadenosine (NECA) binding in the A_2AAR could be still stabilized upon changing the G_s protein to the G_i. Free energy calculations identified one stable low-energy conformation for each of the A₁AR-G_i and A_2AAR-G_s complexes as in the cryo-EM structures, similarly for the A_2AAR-G_i complex. In contrast, the ADO agonist and G_s protein sampled multiple conformations in the A₁AR-G_s system. GaMD simulations thus indicated that the A₁AR preferred to couple with the G_i protein to the G_s, while the A_2AAR could couple with both the G_s and G_i proteins, being highly consistent with experimental findings of the ARs. More importantly, detailed analysis of the atomic simulations showed that the specific AR-G protein coupling resulted from remarkably complementary residue interactions at the protein interface, involving mainly the receptor transmembrane 6 helix and the G_α α5 helix and α4-β6 loop. In summary, the GaMD simulations have provided unprecedented insights into the dynamic mechanism of specific GPCR-G protein interactions at an atomistic level.