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

• Main Authors: Marcus Conrad, Christian A. Söldner, Yinglong Miao, Heinrich Sticht
• Format: Article
• Language: English
• Published: MDPI AG 2020-09-01
• Series: International Journal of Molecular Sciences
• Subjects: receptor–ligand interactions; G protein-coupled receptors (GPCRs); G_s protein; ternary complex; molecular dynamics simulations; metadynamics
• Online Access: https://www.mdpi.com/1422-0067/21/18/6693

The histamine H<inline-formula><math display="inline"><semantics><msub><mrow></mrow><mn>2</mn></msub></semantics></math></inline-formula> receptor (H₂R) plays an important role in the regulation of gastric acid secretion. Therefore, it is a main drug target for the treatment of gastroesophageal reflux or peptic ulcer disease. However, there is as of yet no 3D-structural information available hampering a mechanistic understanding of H₂R. Therefore, we created a model of the histamine-H₂R-G_s complex based on the structure of the ternary complex of the <inline-formula><math display="inline"><semantics><msub><mi mathvariant="sans-serif">β</mi><mn>2</mn></msub></semantics></math></inline-formula>-adrenoceptor and investigated the conformational stability of this active GPCR conformation. Since the physiologically relevant motions with respect to ligand binding and conformational changes of GPCRs can only partly be assessed on the timescale of conventional MD (cMD) simulations, we also applied metadynamics and Gaussian accelerated molecular dynamics (GaMD) simulations. A multiple walker metadynamics simulation in combination with cMD was applied for the determination of the histamine binding mode. The preferential binding pose detected is in good agreement with previous data from site directed mutagenesis and provides a basis for rational ligand design. Inspection of the H₂R-G_s interface reveals a network of polar interactions that may contribute to H₂R coupling selectivity. The cMD and GaMD simulations demonstrate that the active conformation is retained on a <inline-formula><math display="inline"><semantics><mi mathvariant="sans-serif">μ</mi></semantics></math></inline-formula>s-timescale in the ternary histamine-H₂R-G_s complex and in a truncated complex that contains only G_s helix <inline-formula><math display="inline"><semantics><mi mathvariant="sans-serif">α</mi></semantics></math></inline-formula>5 instead of the entire G protein. In contrast, histamine alone is unable to stabilize the active conformation, which is in line with previous studies of other GPCRs.