Computational methods for studying G protein-coupled receptors (GPCRs)

PACAP key interactions with PAC1, VPAC1, and VPAC2 identified by molecular dynamics simulations

The neuropeptide pituitary adenylate cyclase-activating polypeptide (PACAP) belongs to the glucagon/secretin family. PACAP interacts with the pituitary adenylate cyclase-activating polypeptide receptor type 1 (PAC1) and vasoactive intestinal peptide receptors 1 and 2 (VPAC1 and VPAC2), exhibiting functions in the immune, endocrine, and nervous systems. This peptide is upregulated in numerous instances of brain injury, acting as a neuroprotective agent. It can also suppress HIV-1 and SARS-CoV-2 viral replication in vitro. This work aimed to identify, in each peptide-receptor system, the most relevant residues for complex stability and interaction energy communication via Molecular Dynamics (MD), Free Energy calculations, and Protein-energy networks, thus revealing in detail the underlying mechanisms of activation of these receptors. Hydrogen bond formation, interaction energies, and computational alanine scanning between PACAP and its receptors showed that His1, Asp3, Arg12, Arg14, and Lys15 are crucial to the peptide's stability. Furthermore, several PACAP interactions with structurally conserved positions deemed necessary in GPCR B1 activation, including Arg2.60, Lys2.67, and Glu7.42, were significant for the peptide's stability within the receptors. According to the protein-energy network, the connection between Asp3 of PACAP and the receptors' conserved Arg2.60 represents a critical energy communication hub in all complexes. Additionally, the ECDs of the receptors were also found to function as energy communication hubs for PACAP. Although the overall binding mode of PACAP in the three receptors was found to be highly conserved, Arg12 and Tyr13 of PACAP were more prominent in complex with PAC1, while Ser2 of PACAP was with VPAC2. The detailed analyses performed in this work pave the way for using PACAP and its receptors as therapeutic targets.Communicated by Ramaswamy H. Sarma.

Dynamics of the Apo µ-Opioid Receptor in Complex with Gi Protein

Opioid receptors, particularly the µ-opioid receptor (μOR), play a pivotal role in mediating the analgesic and addictive effects of opioid drugs. G protein signaling is an important pathway of μOR function, usually associated with painkilling effects. However, the molecular mechanisms underlying the interaction between the μOR and G protein remain poorly understood. In this study, we employed classical all-atom molecular dynamics simulations to investigate the structural changes occurring with the μOR-G protein complex under two different conditions: with the G protein in the apo form (open) and with the GDP bound G protein (closed, holo form). The receptor was in the apo form and active conformation in both cases, and the simulation time comprised 1µs for each system. In order to assess the effect of the G protein coupling on the receptor activation state, three parameters were monitored: the correlation of the distance between TM3 and TM6 and the RMSD of the NPxxYA motif; the universal activation index (A100); and the χ2 dihedral distribution of residue W2936.48. When complexed with the open G protein, receptor conformations with intermediate activation state prevailed throughout the molecular dynamics, whereas in the condition with the closed G protein, mostly inactive conformations of the receptor were observed. The major effect of the G protein in the receptor conformation comes from a steric hindrance involving an intracellular loop of the receptor and a β-sheet region of the G protein. This suggests that G-protein precoupling is essential for receptor activation, but this fact is not sufficient for complete receptor activation.

In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice

The aim of this study is to evaluate the anticonvulsant potential of schisandrin B, a main ingredient of Schisandra chinensis extracts. Schisandrin B showed anticonvulsant activity in the zebrafish larva pentylenetetrazole acute seizure assay but did not alter seizure thresholds in the intravenous pentylenetetrazole test in mice. Schisandrin B crosses the blood-brain barrier, which we confirmed in our in silico and in vivo analyses; however, the low level of its unbound fraction in the mouse brain tissue may explain the observed lack of anticonvulsant activity. Molecular docking revealed that the anticonvulsant activity of the compound in larval zebrafish might have been due to its binding to a benzodiazepine site within the GABAA receptor and/or the inhibition of the glutamate NMDA receptor. Although schisandrin B showed a beneficial anticonvulsant effect, toxicological studies revealed that it caused serious developmental impairment in zebrafish larvae, underscoring its teratogenic properties. Further detailed studies are needed to precisely identify the properties, pharmacological effects, and safety of schisandrin B.

Synthesis, Structural and Behavioral Studies of Indole Derivatives D2AAK5, D2AAK6 and D2AAK7 as Serotonin 5-HT1A and 5-HT2A Receptor Ligands

Serotonin receptors are involved in a number of physiological functions and regulate aggression, anxiety, appetite, cognition, learning, memory, mood, nausea, sleep, and thermoregulation. Here we report synthesis and detailed structural and behavioral studies of three indole derivatives: D2AAK5, D2AAK6, and D2AAK7 as serotonin 5-HT_1A and 5-HT_2A receptor ligands. X-ray studies revealed that the D2AAK5 compound crystallizes in centrosymmetric triclinic space group with one molecule in the asymmetric unit. The main interaction between the ligands and the receptors is the salt bridge between the protonatable nitrogen atom of the ligands and the conserved Asp (3.32) of the receptors. The complexes were stable in the molecular dynamic simulations. MD revealed that the studied ligands are relatively stable in their binding sites, with the exception of D2AAK7 in the serotonin 5-HT_1A receptor. D2AAK7 exerts anxiolytic activity in the EPM test, while D2AAK5 has a beneficial effect on the memory processes in the PA test.

Synthesis, Docking Studies and Pharmacological Evaluation of Serotoninergic Ligands Containing a 5-Norbornene-2-Carboxamide Nucleus

A new series of 5-norbornene-2-carboxamide derivatives was prepared and their affinities to the 5-HT_1A, 5-HT_2A, and 5-HT_2C receptors were evaluated and compared to a previously synthesized series of derivatives characterized by exo-N-hydroxy-5-norbornene-2,3-dicarboximidenucleus, in order to identify selective ligands for the above-mentioned subtype receptors. Arylpiperazines represents one of the most important classes of 5-HT_1AR ligands, and recent research concerning new derivatives has been focused on the modification of one or more portions of such pharmacophore. The combination of structural elements (heterocyclic nucleus, propyl chain and 4-substituted piperazine), known to be critical to the affinity to 5-HT_1A receptors, and the proper selection of substituents led to compounds with high specificity and affinity towards serotoninergic receptors. The most active compounds were selected for further in vivo assays to determine their functional activity. Finally, to rationalize the obtained results, molecular docking studies were performed. The results of the pharmacological studies showed that Norbo-4 and Norbo-18 were the most active and promising derivatives for the serotonin receptor considered in this study.

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.

N-(3-{4-[3-(trifluoromethyl)phenyl]piperazin-1-yl}propyl)-1H-indazole-3-carboxamide (D2AAK3) as a potential antipsychotic: In vitro, in silico and in vivo evaluation of a multi-target ligand

Schizophrenia is a mental illness of not adequately understood causes that is not satisfactorily enough treated by current antipsychotics. In search for novel potential antipsychotics we performed structure-based virtual screening aimed to identify new dopamine D₂ receptor antagonists. We found compound D2AAK3 with affinity to dopamine D₂ receptor of 115 nM. D2AAK3 possesses additional nanomolar or low micromolar affinity to D₁, D₃, 5-HT_1A, 5-HT_2A and 5-HT₇ receptors, which makes it a good hit for further development as a multifunctional ligand. The compound has also some affinity to M₁ and H₁ receptors. We used homology modeling, molecular docking and molecular dynamics to study interactions of D2AAK3 with its molecular targets at the molecular level. In behavioral studies D2AAK3 decreases amphetamine-induced hyperactivity (when compared to the amphetamine-treated group) measured as spontaneous locomotor activity in mice. In addition, passive avoidance test demonstrated that D2AAK3 improves memory consolidation after acute treatment in mice. Elevated plus maze tests indicated that D2AAK3 induces anxiogenic activity 30 min after acute treatment, whereas this effect has no longer been observed 60 min after administration of the studied compound in mice.

Synthesis, docking studies, and pharmacological evaluation of 2-hydroxypropyl-4-arylpiperazine derivatives as serotoninergic ligands

A new series of norbornene and exo-N-hydroxy-7-oxabicyclo[2.2.1]hept-5-ene-2,3-dicarboximide derivatives was prepared, and their affinities to the 5-HT_1A , 5-HT_2A , and 5-HT_2C receptors were evaluated and compared with a previously synthesized series of derivatives characterized by the same nuclei, to identify selective ligands for the subtype receptors. Arylpiperazines represent one of the most important classes of 5-HT_1A R ligands, and the research of new derivatives has been focused on the modification of one or more portions of this pharmacophore. The combination of structural elements (heterocyclic nucleus, hydroxyalkyl chain, and 4-substituted piperazine), known to be critical for the affinity to 5-HT_1A receptors, and the proper selection of substituents resulted in compounds with high specificity and affinity toward serotoninergic receptors. The most active compounds were selected for further in vivo assays to determine their functional activity. Finally, to rationalize the obtained results, molecular docking studies were performed. The results of the pharmacological studies showed that 3e, 4j, and 4n were the most active and promising derivatives for the serotonin receptor considered in this study.

Illuminating the Path to Target GPCR Structures and Functions

G-Protein-coupled receptors (GPCRs) are ubiquitous within eukaryotes, responsible for a wide array of physiological and pathological processes. Indeed, the fact that they are the most drugged target in the human genome is indicative of their importance. Despite the clear interest in GPCRs, most information regarding their activity has been so far obtained by analyzing the response from a "bulk medium". As such, this Perspective summarizes some of the common methods for this indirect observation. Nonetheless, by inspecting approaches applying super-resolution imaging, we argue that imaging is perfectly situated to obtain more detailed structural and spatial information, assisting in the development of new GPCR-targeted drugs and clinical strategies. The benefits of direct optical visualization of GPCRs are analyzed in the context of potential future directions in the field.

Ligand with Two Modes of Interaction with the Dopamine D2 Receptor-An Induced-Fit Mechanism of Insurmountable Antagonism

A solid understanding of the mechanisms governing ligand binding is crucial for rational design of therapeutics targeting the dopamine D2 receptor (D2R). Here, we use G protein-coupled inward rectifier potassium (GIRK) channel activation in Xenopus oocytes to measure the kinetics of D2R antagonism by a series of aripiprazole analogues, as well as the recovery of dopamine (DA) responsivity upon washout. The aripiprazole analogues comprise an orthosteric and a secondary pharmacophore and differ by the length of the saturated carbon linker joining these two pharmacophores. Two compounds containing 3- and 5-carbon linkers allowed for a similar extent of recovery from antagonism in the presence of 1 or 100 μM DA (>25 and >90% of control, respectively), whereas recovery was less prominent (∼20%) upon washout of the 4-carbon linker compound, SV-III-130, both with 1 and 100 μM DA. Prolonging the coincubation time with SV-III-130 further diminished recovery. Curve-shift experiments were consistent with competition between SV-III-130 and DA. Two mutations in the secondary binding pocket (V91A and E95A) of D2R decreased antagonistic potency and increased recovery from SV-III-130 antagonism, whereas a third mutation (L94A) only increased recovery. Our results suggest that the secondary binding pocket influences recovery from inhibition by the studied aripiprazole analogues. We propose a mechanism, supported by in silico modeling, whereby SV-III-130 initially binds reversibly to the D2R, after which the drug-receptor complex undergoes a slow transition to a second ligand-bound state, which is dependent on secondary binding pocket integrity and irreversible during the time frame of our experiments.

Recent Insights from Molecular Dynamics Simulations for G Protein-Coupled Receptor Drug Discovery

G protein-coupled receptors (GPCRs) are critical drug targets. GPCRs convey signals from the extracellular to the intracellular environment through G proteins. Some ligands that bind to GPCRs activate different downstream signaling pathways. G protein activation, or -arrestin biased signaling, involves ligands binding to receptors and stabilizing conformations that trigger a specific pathway. -arrestin biased signaling has become a hot target for structure-based drug discovery. However, challenges include that there are few crystal structures available in the Protein Data Bank and that GPCRs are highly dynamic. Hence, molecular dynamics (MD) simulations are especially valuable for obtaining detailed mechanistic information, including identification of allosteric sites and understanding modulators' interactions with receptors and ligands. Here, we highlight recent MD simulation studies and enhanced sampling methods used to study biased G protein-coupled receptor signaling and their conformational dynamics as well as applications to drug discovery.

Importance of protein dynamics in the structure-based drug discovery of class A G protein-coupled receptors (GPCRs)

Demand for novel GPCR modulators is increasing as the association between the GPCR signaling pathway and numerous diseases such as cancers, psychological and metabolic disorders continues to be established. In silico structure-based drug design (SBDD) offers an outlet where researchers could exploit the accumulating structural information of GPCR to expedite the process of drug discovery. The coupling of structure-based approaches such as virtual screening and molecular docking with molecular dynamics and/or Monte Carlo simulation aids in reflecting the dynamics of proteins in nature into previously static docking studies, thus enhancing the accuracy of rationally designed ligands. This review will highlight recent computational strategies that incorporate protein flexibility into SBDD of GPCR-targeted ligands.

Distinct Dopamine D₂ Receptor Antagonists Differentially Impact D₂ Receptor Oligomerization

Dopamine D₂ receptors (D₂R) are known to form transient homodimer complexes, of which the increased formation has already been associated with development of schizophrenia. Pharmacological targeting and modulation of the equilibrium of these receptor homodimers might lead to a better understanding of the critical role played by these complexes in physiological and pathological conditions. Whereas agonist addition has shown to prolong the D₂R dimer lifetime and increase the level of dimer formation, the possible influence of D₂R antagonists on dimerization has remained rather unexplored. Here, using a live-cell reporter assay based on the functional complementation of a split Nanoluciferase, a panel of six D₂R antagonists were screened for their ability to modulate the level of D_2LR dimer formation. Incubation with the D₂R antagonist spiperone decreased the level of D_2LR dimer formation significantly by 40–60% in real-time and after long-term (≥16 h) incubations. The fact that dimer formation of the well-studied A_2a–D_2LR dimer was not altered following incubation with spiperone supports the specificity of this observation. Other D₂R antagonists, such as clozapine, risperidone, and droperidol did not significantly evoke this dissociation event. Furthermore, molecular modeling reveals that spiperone presents specific Tyr199^5.48 and Phe390^6.52 conformations, compared to clozapine, which may determine D₂R homodimerization.

Understanding G Protein-Coupled Receptor Allostery via Molecular Dynamics Simulations: Implications for Drug Discovery

Unraveling the mystery of protein allostery has been one of the greatest challenges in both structural and computational biology. However, recent advances in computational methods, particularly molecular dynamics (MD) simulations, have led to its utility as a powerful and popular tool for the study of protein allostery. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric hot spots and the determination of the mechanistic basis for allostery. These structural and dynamic studies can provide a foundation for a wide range of applications, including rational drug design and protein engineering. In our laboratory, the use of MD simulations and network analysis assisted in the elucidation of the allosteric hotspots and intracellular signal transduction of G protein-coupled receptors (GPCRs), primarily on one of the adenosine receptor subtypes, A_2A adenosine receptor (A_2AAR). In this chapter, we describe a method for calculating the map of allosteric signal flow in different GPCR conformational states and illustrate how these concepts have been utilized in understanding the mechanism of GPCR allostery. These structural studies will provide valuable insights into the allosteric and orthosteric modulations that would be of great help to design novel drugs targeting GPCRs in pathological states.

Computational Methods Used in Hit-to-Lead and Lead Optimization Stages of Structure-Based Drug Discovery

GPCR modeling approaches are widely used in the hit-to-lead (H2L) and lead optimization (LO) stages of drug discovery. The aims of these modeling approaches are to predict the 3D structures of the receptor-ligand complexes, to explore the key interactions between the receptor and the ligand and to utilize these insights in the design of new molecules with improved binding, selectivity or other pharmacological properties. In this book chapter, we present a brief survey of key computational approaches integrated with hierarchical GPCR modeling protocol (HGMP) used in hit-to-lead (H2L) and in lead optimization (LO) stages of structure-based drug discovery (SBDD). We outline the differences in modeling strategies used in H2L and LO of SBDD and illustrate how these tools have been applied in three drug discovery projects.

Opportunities and Challenges in the Discovery of Allosteric Modulators of GPCRs

From the pharmacological point of view, allosteric modulators may present numerous advantages over orthosteric ligands. Growing availability of novel tools and experimental data provides a tempting opportunity to apply computational methods to improve known modulators and design novel ones. However, recent progress in understanding of complexity of allostery increases awareness of problems involved in design of modulators with desired properties. Deeper insight into phenomena such as probe dependence, altering signaling bias with minor changes in ligand structure, as well as influence of subtle endogenous allosteric factors turns out to be fundamental. These effects make the design of a modulator with precise pharmacological outcome a very challenging task, and need to be taken into consideration throughout the design process. In this chapter, we focus on nuances of targeting GPCR allosteric sites in computational drug design efforts, in particular with application of docking, virtual screening, and molecular dynamics.

Challenges and Opportunities in Drug Discovery of Biased Ligands

The observation of biased agonism in G protein-coupled receptors (GPCRs) has provided new approaches for the development of more efficacious and safer drugs. However, in order to rationally design biased drugs, one must understand the molecular basis of this phenomenon. Computational approaches can help in exploring the conformational universe of GPCRs and detecting conformational states with relevance for distinct functional outcomes. This information is extremely valuable for the development of new therapeutic agents that promote desired conformational receptor states and responses while avoiding the ones leading to undesired side-effects.This book chapter intends to introduce the reader to powerful computational approaches for sampling the conformational space of these receptors, focusing first on molecular dynamics and the analysis of the produced data through methods such as dimensionality reduction, Markov State Models and adaptive sampling. Then, we show how to seek for compounds that target distinct conformational states via docking and virtual screening. In addition, we describe how to detect receptor-ligand interactions that drive signaling bias and comment current challenges and opportunities of presented methods.

Recent Advances in Structure-Based Drug Design Targeting Class A G Protein-Coupled Receptors Utilizing Crystal Structures and Computational Simulations

G protein-coupled receptors (GPCRs) represent the largest and most physiologically important integral membrane protein family, and these receptors respond to a wide variety of physiological and environmental stimuli. GPCRs are among the most critical therapeutic targets for numerous human diseases, and approximately one-third of the currently marketed drugs target this receptor family. The recent breakthroughs in GPCR structural biology have significantly contributed to our understanding of GPCR function, ligand binding, and pharmacological action as well as to the design of new drugs. This perspective highlights the latest advances in GPCR structures with a focus on the receptor-ligand interactions of each receptor family in class A nonrhodopsin GPCRs as well as the structural features for their activation, biased signaling, and allosteric mechanisms. The current state-of-the-art methodologies of structure-based drug design (SBDD) approaches in the GPCR research field are also discussed.

Signaling within Allosteric Machines: Signal Transmission Pathways Inside G Protein-Coupled Receptors

In recent years, our understanding of function of G protein-coupled receptors (GPCRs) has changed from a picture of simple signal relays, transmitting only a particular signal to a particular G protein heterotrimer, to versatile machines, capable of various responses to different stimuli and being modulated by various factors. Some recent reports provide not only the data on ligands/modulators and resultant signals induced by them, but also deeper insights into exact pathways of signal migration and mechanisms of signal transmission through receptor structure. Combination of these computational and experimental data sheds more light on underlying mechanisms of signal transmission and signaling bias in GPCRs. In this review we focus on available clues on allosteric pathways responsible for complex signal processing within GPCRs structures, with particular emphasis on linking compatible in silico- and in vitro-derived data on the most probable allosteric connections.

Molecular Docking, Molecular Dynamics, and Structure-Activity Relationship Explorations of 14-Oxygenated N-Methylmorphinan-6-ones as Potent μ-Opioid Receptor Agonists

Among opioids, morphinans are of major importance as the most effective analgesic drugs acting primarily via μ-opioid receptor (μ-OR) activation. Our long-standing efforts in the field of opioid analgesics from the class of morphinans led to N-methylmorphinan-6-ones differently substituted at positions 5 and 14 as μ-OR agonists inducing potent analgesia and fewer undesirable effects. Herein we present the first thorough molecular modeling study and structure-activity relationship (SAR) explorations aided by docking and molecular dynamics (MD) simulations of 14-oxygenated N-methylmorphinan-6-ones to gain insights into their mode of binding to the μ-OR and interaction mechanisms. The structure of activated μ-OR provides an essential model for how ligand/μ-OR binding is encoded within small chemical differences in otherwise structurally similar morphinans. We reveal important molecular interactions that these μ-agonists share and distinguish them. The molecular docking outcomes indicate the crucial role of the relative orientation of the ligand in the μ-OR binding site, influencing the propensity of critical non-covalent interactions that are required to facilitate ligand/μ-OR interactions and receptor activation. The MD simulations point out minor differences in the tendency to form hydrogen bonds by the 4,5α-epoxy group, along with the tendency to affect the 3-7 lock switch. The emerged SARs reveal the subtle interplay between the substituents at positions 5 and 14 in the morphinan scaffold by enabling the identification of key structural elements that determine the distinct pharmacological profiles. This study provides a significant structural basis for understanding ligand binding and μ-OR activation by the 14-oxygenated N-methylmorphinan-6-ones, which should be useful for guiding drug design.

Recent Advances and Applications of Molecular Docking to G Protein-Coupled Receptors

The growing number of studies on G protein-coupled receptors (GPCRs) family are a source of noticeable improvement in our understanding of the functioning of these proteins. GPCRs are responsible for a vast part of signaling in vertebrates and, as such, invariably remain in the spotlight of medicinal chemistry. A deeper insight into the underlying mechanisms of interesting phenomena observed in GPCRs, such as biased signaling or allosteric modulation, can be gained with experimental and computational studies. The latter play an important role in this process, since they allow for observations on scales inaccessible for most other methods. One of the key steps in such studies is proper computational reconstruction of actual ligand-receptor or protein-protein interactions, a process called molecular docking. A number of improvements and innovative applications of this method were documented recently. In this review, we focus particularly on innovations in docking to GPCRs.

Interplay of G Protein-Coupled Receptors with the Membrane: Insights from Supra-Atomic Coarse Grain Molecular Dynamics Simulations

G protein-coupled receptors (GPCRs) are central to many fundamental cellular signaling pathways. They transduce signals from the outside to the inside of cells in physiological processes ranging from vision to immune response. It is extremely challenging to look at them individually using conventional experimental techniques. Recently, a pseudo atomistic molecular model has emerged as a valuable tool to access information on GPCRs, more specifically on their interactions with their environment in their native cell membrane and the consequences on their supramolecular organization. This approach uses the Martini coarse grain (CG) model to describe the receptors, lipids, and solvent in molecular dynamics (MD) simulations and in enough detail to allow conserving the chemical specificity of the different molecules. The elimination of unnecessary degrees of freedom has opened up large-scale simulations of the lipid-mediated supramolecular organization of GPCRs. Here, after introducing the Martini CGMD method, we review these studies carried out on various members of the GPCR family, including rhodopsin (visual receptor), opioid receptors, adrenergic receptors, adenosine receptors, dopamine receptor, and sphingosine 1-phosphate receptor. These studies have brought to light an interesting set of novel biophysical principles. The insights range from revealing localized and heterogeneous deformations of the membrane bilayer at the surface of the protein, specific interactions of lipid molecules with individual GPCRs, to the effect of the membrane matrix on global GPCR self-assembly. The review ends with an overview of the lessons learned from the use of the CGMD method, the biophysical-chemical findings on lipid-protein interplay.

The dopamine D2 receptor dimer and its interaction with homobivalent antagonists: homology modeling, docking and molecular dynamics

In order to apply structure-based drug design techniques to G protein-coupled receptor complexes, it is essential to model their 3D structure and to identify regions that are suitable for selective drug binding. For this purpose, we have developed and tested a multi-component protocol to model the inactive conformation of the dopamine D2 receptor dimer, suitable for interaction with homobivalent antagonists. Our approach was based on protein-protein docking, applying the Rosetta software to obtain populations of dimers as present in membranes with all the main possible interfaces. Consensus scoring based on the values and frequencies of best interfaces regarding four scoring parameters, Rosetta interface score, interface area, free energy of binding and energy of hydrogen bond interactions indicated that the best scored dimer model possesses a TM4-TM5-TM7-TM1 interface, which is in agreement with experimental data. This model was used to study interactions of the previously published dopamine D2 receptor homobivalent antagonists based on clozapine,1,4-disubstituted aromatic piperidines/piperazines and arylamidoalkyl substituted phenylpiperazine pharmacophores. It was found that the homobivalent antagonists stabilize the receptor-inactive conformation by maintaining the ionic lock interaction, and change the dimer interface by disrupting a set of hydrogen bonds and maintaining water- and ligand-mediated hydrogen bonds in the extracellular and intracellular part of the interface. Graphical Abstract Structure of the final model of the dopamine D2 receptor homodimer, indicating the distancebetween Tyr37 and Tyr 5.42 in the apo form (left) and in the complex with the ligand (right).

Interplay between Two Allosteric Sites and Their Influence on Agonist Binding in Human μ Opioid Receptor

Allostery is a widespread mechanism that allows for precise protein tuning. Its underlying mechanisms are elusive, particularly when there are multiple allosteric sites at the protein. This concerns also G-protein-coupled receptors (GPCRs), which are targets for a vast part of currently used drugs. To address this issue, we performed molecular dynamics simulations of a GPCR-human μ opioid receptor (MOR) in a native-like environment, with full agonist (R)-methadone, Na(+) ions, and a positive modulator BMS986122 in various configurations. We found that MOR's seventh transmembrane helix (TM VII) is central for allosteric signal transmission, and modulators affect its bending and rotation. The PAM stabilizes favorable agonist interactions, while Na(+) tends to disrupt agonist binding. We identified two residues involved in allosteric signal transmission: Trp 7.35 at the top and Tyr 7.53 at the bottom of TM VII.