G-Protein-Coupled Receptor-Membrane Interactions Depend on the Receptor Activation state

Abstract

G‐protein‐coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one‐third of currently marketed drugs. In particular, adenosine A1 receptor (A1AR) is an important therapeutic target for treating cardiac ischemia–reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A1AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid–protein interactions in order to understand the functional mechanism of GPCRs. Here, all‐atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G‐protein bound) A1AR, which was embedded in a 1‐palmitoyl‐2‐oleoyl‐glycero‐3‐phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A1AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A1AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.

In Silico Drug Design for Purinergic GPCRs: Overview on Molecular Dynamics Applied to Adenosine and P2Y Receptors

“…More recently, for the A 1 AR two inactive antagonist-bound structures [ 24 , 38 ] and a cryo-EM active agonist-bound structure [ 39 ] were released, so more reliable molecular modeling studies are now possible. Recently, GaMD investigated the flexibility profile of the active adenosine-bound and inactive PSB36-bound A 1 AR [ 97 ]. In both simulations, conformations of the ECL2 region and the intracellular end of TM6 fluctuated, with higher flexibility of these receptor portions and of the intracellular end of TM5 in the active state.…”

Mechanisms of γ-Secretase Activation and Substrate Processing

“…GaMD has successfully revealed physical pathways and mechanisms of protein folding and ligand binding, which are consistent with experiments and long-timescale conventional MD simulations. 27,[33][34] It has also been applied to characterize protein-protein, [35][36] protein-membrane, 37 and protein-nucleic acid [38][39] interactions. Therefore, GaMD was applied in this study for enhanced sampling of the γ -secretase complex, a well-known slow enzyme.…”

Agonist Binding and G Protein Coupling in Histamine H2 Receptor: A Molecular Dynamics Study

“…Moreover, because the boost potential exhibits a Gaussian distribution, biomolecular free energy profiles can be properly recovered through cumulant expansion to the second order [ 23 ]. GaMD simulations have successfully revealed mechanisms of protein folding and conformational changes [ 23 , 25 , 26 , 27 , 28 , 29 ], ligand binding [ 23 , 26 , 27 , 30 , 31 , 32 , 33 , 34 , 35 ], and protein-protein/membrane/nucleic acid interactions [ 32 , 36 , 37 , 38 , 39 , 40 ].…”