Department of Chemistry, Institute for Biophysical Dynamics, James Franck Institute and Computation Institute, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, USA.
J Am Chem Soc. 2013 Jun 12;135(23):8749-59. doi: 10.1021/ja404391q. Epub 2013 May 29.
G-protein-coupled receptors (GPCRs) are membrane proteins with critical functions in cellular signal transduction, representing a primary class of drug targets. Acting by direct binding, many drugs modulate GPCR activity and influence the signaling pathways associated with numerous diseases. However, complete details of ligand-dependent GPCR activation/deactivation are difficult to obtain from experiments. Therefore, it remains unclear how ligands modulate a GPCR's activity. To elucidate the ligand-dependent activation/deactivation mechanism of the human adenosine A2A receptor (AA2AR), a member of the class A GPCRs, we performed large-scale unbiased molecular dynamics and metadynamics simulations of the receptor embedded in a membrane. At the atomic level, we have observed distinct structural states that resemble the active and inactive states. In particular, we noted key structural elements changing in a highly concerted fashion during the conformational transitions, including six conformational states of a tryptophan (Trp246(6.48)). Our findings agree with a previously proposed view that, during activation, this tryptophan residue undergoes a rotameric transition that may be coupled to a series of coherent conformational changes, resulting in the opening of the G-protein binding site. Further, metadynamics simulations provide quantitative evidence for this mechanism, suggesting how ligand binding shifts the equilibrium between the active and inactive states. Our analysis also proposes that a few specific residues are associated with agonism/antagonism, affinity, and selectivity, and suggests that the ligand-binding pocket can be thought of as having three distinct regions, providing dynamic features for structure-based design. Additional simulations with AA2AR bound to a novel ligand are consistent with our proposed mechanism. Generally, our study provides insights into the ligand-dependent AA2AR activation/deactivation in addition to what has been found in crystal structures. These results should aid in the discovery of more effective and selective GPCR ligands.
G 蛋白偶联受体(GPCRs)是细胞信号转导中具有关键功能的膜蛋白,是药物靶点的主要类别之一。许多药物通过直接结合来调节 GPCR 活性并影响与许多疾病相关的信号通路,这些药物作用于 GPCR。然而,从实验中很难获得配体依赖性 GPCR 激活/失活的完整细节。因此,配体如何调节 GPCR 的活性仍不清楚。为了阐明人腺苷 A2A 受体(AA2AR)的配体依赖性激活/失活机制,AA2AR 是 A 类 GPCR 中的一员,我们对嵌入膜中的受体进行了大规模无偏的分子动力学和元动力学模拟。在原子水平上,我们观察到了与活性和非活性状态相似的独特结构状态。特别是,我们注意到在构象转变过程中,关键结构元素以高度协调的方式发生变化,包括色氨酸(Trp246(6.48))的六个构象状态。我们的发现与先前提出的观点一致,即在激活过程中,该色氨酸残基经历构象转变,可能与一系列相干构象变化相关联,导致 G 蛋白结合位点打开。此外,元动力学模拟为该机制提供了定量证据,表明配体结合如何使活性和非活性状态之间的平衡发生转移。我们的分析还提出了几个特定的残基与激动剂/拮抗剂、亲和力和选择性有关,并表明配体结合口袋可以被认为具有三个不同的区域,为基于结构的设计提供了动态特征。与新型配体结合的 AA2AR 的附加模拟与我们提出的机制一致。总体而言,我们的研究除了在晶体结构中发现的内容外,还提供了对配体依赖性 AA2AR 激活/失活的深入了解。这些结果应该有助于发现更有效和选择性的 GPCR 配体。