Conserved water-mediated hydrogen bond network between TM-I, -II, -VI, and -VII in 7TM receptor activation

Activation fingerprints and allosteric modulation at the free fatty acid receptor 1 (FFAR1) revealed by molecular dynamics simulation

The free fatty acid receptor 1 (FFAR1) is a transmembrane G-protein coupled receptor that mediates the metabolic and insulinotropic effects of endogenous free fatty acids in pancreatic cells while also exerting neuro-regulatory effects in the brain. The complexity of FFAR1 derives from its multiple binding sites and the absence of conventional activation motifs observed in class A GPCRs. This study uses molecular dynamics simulations to investigate the molecular mechanisms that underpin endogenous signaling and allosteric regulation in the FFAR1. We investigated and compared three ligand-bound states and the APO state. The ligand-bound simulations included FFAR1 in complex with γ-linolenic acid, FFAR1 in complex with γ-linolenic acid and TAK875, and a fully activated FFAR1 bundle complexed with docosahexaenoic acid and G-protein. The results highlight distinct protein contact fingerprints and dynamics in the ligand-bound states relative to the APO state. While ligand binding, in the absence of stabilizing G-protein, destabilizes the intracellular domain of the receptor, the second extracellular loop exhibits greater stability and salt bridge contact with the transmembrane domain. Notably, simulations of FFAR1 complexed with γ-linolenic acid, bound at the intracellular domain, revealed stable interactions between γ-linolenic acid and the receptor, as well as similar activation fingerprints when compared to FFAR1 in complex with docosahexaenoic acid and Gq. This suggests an effective allosteric regulation of the receptor following γ-linolenic acid binding to the intracellular domain. Finally, a set of hydrophobic amino acid residues at the intracellular and extracellular domains appears to function as potential rotameric switches, facilitating water-mediated receptor activation.

ART-SM: Boosting Fragment-Based Backmapping by Machine Learning

In sequential multiscale molecular dynamics simulations, which advantageously combine the increased sampling and dynamics at coarse-grained resolution with the higher accuracy of atomistic simulations, the resolution is altered over time. While coarse-graining is straightforward once the mapping between atomistic and coarse-grained resolution is defined, reintroducing the atomistic details is still a nontrivial process called backmapping. Here, we present ART-SM, a fragment-based backmapping framework that learns from atomistic simulation data to seamlessly switch from coarse-grained to atomistic resolution. ART-SM requires minimal user input and goes beyond state-of-the-art fragment-based approaches by selecting from multiple conformations per fragment via machine learning to simultaneously reflect the coarse-grained structure and the Boltzmann distribution. Additionally, we introduce a novel refinement step to connect individual fragments by optimizing specific bonds, angles, and dihedral angles in the backmapping process. We demonstrate that our algorithm accurately restores the atomistic bond length, angle, and dihedral angle distributions for various small and linear molecules from Martini coarse-grained beads and that the resulting high-resolution structures are representative of the input coarse-grained conformations. Moreover, the reconstruction of the TIP3P water model is fast and robust, and we demonstrate that ART-SM can be applied to larger linear molecules as well. To illustrate the efficiency of the local and autoregressive approach of ART-SM, we simulated a large realistic system containing the surfactants TAPB and SDS in solution using the Martini3 force field. The self-assembled micelles of various shapes were backmapped with ART-SM after training on only short atomistic simulations of a single water-solvated SDS or TAPB molecule. Together, these results indicate the potential for the method to be extended to more complex molecules such as lipids, proteins, macromolecules, and materials in the future.

Computational design of highly signalling-active membrane receptors through solvent-mediated allosteric networks

Protein catalysis and allostery require the atomic-level orchestration and motion of residues and ligand, solvent and protein effector molecules. However, the ability to design protein activity through precise protein-solvent cooperative interactions has not yet been demonstrated. Here we report the design of 14 membrane receptors that catalyse G protein nucleotide exchange through diverse engineered allosteric pathways mediated by cooperative networks of intraprotein, protein-ligand and -solvent molecule interactions. Consistent with predictions, the designed protein activities correlated well with the level of plasticity of the networks at flexible transmembrane helical interfaces. Several designs displayed considerably enhanced thermostability and activity compared with related natural receptors. The most stable and active variant crystallized in an unforeseen signalling-active conformation, in excellent agreement with the design models. The allosteric network topologies of the best designs bear limited similarity to those of natural receptors and reveal an allosteric interaction space larger than previously inferred from natural proteins. The approach should prove useful for engineering proteins with novel complex protein binding, catalytic and signalling activities.

Visualizing Dynamic Weak Interaction Networks of Fluorine Atoms within Proteins via NMR

The utilization of fluorine probes in protein research has advanced our understanding of the nature of macromolecules. The diverse activities of proteins are governed by intricate weak interaction networks, yet the experimental detection presents a challenge. Herein, we have developed an NMR analytical method based on quantum coherent operations to characterize the weak chemical bonds of fluorine atoms within proteins that are bound to metal fluorides or labeled with fluorinated amino acids. Our approach offers a fast-screening method for identifying a weak interaction network without requiring crystallization.

Implications of β-Arrestin biased signaling by angiotensin II type 1 receptor for cardiovascular drug discovery and therapeutics

Angiotensin II receptors, Type 1 (AT1R) and Type 2 (AT2R) are 7TM receptors that play critical roles in both the physiological and pathophysiological regulation of the cardiovascular system. While AT1R blockers (ARBs) have proven beneficial in managing cardiac, vascular and renal maladies they cannot completely halt and reverse the progression of pathologies. Numerous experimental and animal studies have demonstrated that β-arrestin biased AT1R-ligands (such as SII-AngII, S1I8, TRV023, and TRV027) offer cardiovascular benefits by blocking the G protein signaling while retaining the β-arrestin signaling. However, these ligands failed to show improvement in heart-failure outcome over the placebo in a phase IIb clinical trial. One major limitation of current β-arrestin biased AT1R-ligands is that they are peptides with short half-lives, limiting their long-term efficacy in patients. Additionally, β-arrestin biased AT1R-ligand peptides, may inadvertently block AT2R, a promiscuous receptor, potentially negating its beneficial effects in post-myocardial infarction (MI) patients. Therefore, developing a small molecule β-arrestin biased AT1R-ligand with a longer half-life and specificity to AT1R could be more effective in treating heart failure. This approach has the potential to revolutionize the treatment of cardiovascular diseases by offering more sustained and targeted therapeutic effects.

Molecular determinants of neuropeptide-mediated activation mechanisms in tachykinin NK1 and NK2 receptors

Substance P and neurokinin A are closely related neuropeptides belonging to the tachykinin family. Their receptors are neurokinin one receptor (NK1R) and neurokinin two receptor (NK2R), G protein-coupled receptors that transmit Gs and Gq-mediated downstream signaling. We investigate the importance of sequence differences at the bottom of the receptor orthosteric site for activity and selectivity, focusing on residues that closely interact with the C-terminal methionine of the peptide ligands. We identify a conserved serine (NK1R-S2977.45) and the position of the tryptophan residue within the canonical "toggle switch" motif, CWxP of TM6, neighboring a phenylalanine in NK1R (NK1R-F2646.51) and a tyrosine in NK2R (NK2R-Y2666.51), giving rise to distinct microenvironments for the neuropeptide C terminals. Mutating these residues results in dramatic activity changes in both NK1R and NK2R due to a close interaction between the ligand and toggle switch. Structural analysis of active and inactive NKR structures suggests only a minor change in sidechain rotation of toggle switch residues upon activation. However, extensive molecular dynamics simulations of receptor:neuropeptide:G protein complexes indicate that a major, concerted motion happens in the toggle switch tryptophan indole group and the sidechains of the microswitch motif Pro-Ile-Phe (PIF). This rotation establishes a tight hydrogen bond interaction from the tryptophan indole to the conserved serine (NK1R-S2977.45) and a mainchain carbonyl (NK1R-A2947.41) in the kink of TM7. This interaction facilitates communication with the NPxxY microswitch motif of TM7, resulting in stabilization of the G protein-binding region. NK1R-S2977.45 is consequently identified as a central hub for the activation of NKRs.

Illuminating the structural basis of human neurokinin 1 receptor (NK1R) antagonism through classical all-atoms molecular dynamics simulations

Advances in structural biology have bestowed insights into the pleiotropic effects of neurokinin 1 receptors (NK1R) in diverse patho-physiological processes, thereby highlighting the potential therapeutic value of antagonists directed against NK1R. Herein, we investigate the mode of antagonist recognition to discern the obscure atomic facets germane for the function and molecular determinants of NK1R. To commence discernment of potent antagonists and the conformational changes in NK1R, induced upon antagonist binding, state-of-the-art classical all-atoms molecular dynamics (MD) simulations in lipid mimetic bilayers have been utilized. MD simulations of structural ensembles reveals the involvement of TM5 and TM6 in tight anchoring of antagonists through a network of interhelical hydrogen-bonds, while, the extracellular loop 2 (ECL2) governs the overall size and nature of the pocket, thereby modulating NK1R. Consistent comparison between experiments and MD simulation results discerns the predominant role of TM3, TM4, and TM6 in lipid-NK1R interaction. Correlation between hydrophobic index and helicity of TM domains elucidates their importance in maintaining the structural stability in addition to regulating NK1R antagonism. Taken together, we anticipate that our computational study marks a comprehensive structural basis of NK1R antagonism in lipid bilayers, which may facilitate designing of new therapeutics against associated diseases targeting human neurokinin receptors.

Structural insights into angiotensin receptor signaling modulation by balanced and biased agonists

The peptide hormone angiotensin II regulates blood pressure mainly through the type 1 angiotensin II receptor AT1 R and its downstream signaling proteins Gq and β-arrestin. AT1 R blockers, clinically used as antihypertensive drugs, inhibit both signaling pathways, whereas AT1 R β-arrestin-biased agonists have shown great potential for the treatment of acute heart failure. Here, we present a cryo-electron microscopy (cryo-EM) structure of the human AT1 R in complex with a balanced agonist, Sar1 -AngII, and Gq protein at 2.9 Å resolution. This structure, together with extensive functional assays and computational modeling, reveals the molecular mechanisms for AT1 R signaling modulation and suggests that a major hydrogen bond network (MHN) inside the receptor serves as a key regulator of AT1 R signal transduction from the ligand-binding pocket to both Gq and β-arrestin pathways. Specifically, we found that the MHN mutations N1113.35 A and N2947.45 A induce biased signaling to Gq and β-arrestin, respectively. These insights should facilitate AT1 R structure-based drug discovery for the treatment of cardiovascular diseases.

Enhanced Thermostability of D-Psicose 3-Epimerase from Clostridium bolteae through Rational Design and Engineering of New Disulfide Bridges

D-psicose 3-epimerase (DPEase) catalyzes the isomerization of D-fructose to D-psicose (aka D-allulose, a low-calorie sweetener), but its industrial application has been restricted by the poor thermostability of the naturally available enzymes. Computational rational design of disulfide bridges was used to select potential sites in the protein structure of DPEase from Clostridium bolteae to engineer new disulfide bridges. Three mutants were engineered successfully with new disulfide bridges in different locations, increasing their optimum catalytic temperature from 55 to 65 °C, greatly improving their thermal stability and extending their half-lives (t_1/2) at 55 °C from 0.37 h to 4−4.5 h, thereby greatly enhancing their potential for industrial application. Molecular dynamics simulation and spatial configuration analysis revealed that introduction of a disulfide bridge modified the protein hydrogen–bond network, rigidified both the local and overall structures of the mutants and decreased the entropy of unfolded protein, thereby enhancing the thermostability of DPEase.

Novel Sequence Feature of SecA Translocase Protein Unique to the Thermophilic Bacteria: Bioinformatics Analyses to Investigate Their Potential Roles

SecA is an evolutionarily conserved protein that plays an indispensable role in the secretion of proteins across the bacterial cell membrane. Comparative analyses of SecA homologs have identified two large conserved signature inserts (CSIs) that are unique characteristics of thermophilic bacteria. A 50 aa conserved insert in SecA is exclusively present in the SecA homologs from the orders Thermotogales and Aquificales, while a 76 aa insert in SecA is specific for the order Thermales and Hydrogenibacillus schlegelii. Phylogenetic analyses on SecA sequences show that the shared presence of these CSIs in unrelated groups of thermophiles is not due to lateral gene transfers, but instead these large CSIs have likely originated independently in these lineages due to their advantageous function. Both of these CSIs are located in SecA protein in a surface exposed region within the ATPase domain. To gain insights into the functional significance of the 50 aa CSI in SecA, molecular dynamics (MD) simulations were performed at two different temperatures using ADP-bound SecA from Thermotoga maritima. These analyses have identified a conserved network of water molecules near the 50 aa insert in which the Glu185 residue from the CSI is found to play a key role towards stabilizing these interactions. The results provide evidence for the possible role of the 50 aa CSI in stabilizing the binding interaction of ADP/ATP, which is required for SecA function. Additionally, the surface-exposed CSIs in SecA, due to their potential to make novel protein-protein interactions, could also contribute to the thermostability of SecA from thermophilic bacteria.

The role of water molecules in phototransduction of retinal proteins and G protein-coupled receptors

G protein coupled receptors (GPCRs) are a key family of membrane proteins in all eukaryotes and also very important drug targets for medical intervention. The extensively studied visual pigment rhodopsin is a prime example of a family A GPCR. Its chromophore ligand retinal is covalently linked to a lysine in helix seven forming a protonated Schiff base. Interestingly, this is the same situation in other-non-GPCR-retinal proteins, like the prototype light-driven microbial proton pump bacteriorhodopsin, albeit there is no (or only a very remote) phylogenetical link. Close to the retinal ligand, several water molecules help to organise a functionally important hydrogen bond network that undergoes significant changes during photo-activation. Such water-mediated networks are also critical for ligand binding of other GPCRs and they are becoming increasingly important in drug discovery. GPCRs also contain a partially conserved water mediated hydrogen bond network that stabilises the ground state of the receptor, and rearrangement of this network leads to the stabilization of the active state. Some water molecules have a specific role in this process to appropriately orient specific residues relative to the Schiff base, and to modulate the fine structure of the transmembrane bundle, particularly near the intracellular G protein binding site. While the atomic details of these mechanisms are still missing, the recent developments in free electron lasers (FELs) are enabling us to begin to observe the changes in waters and relevant side chains shortly after photo activation at an unprecedented level of spatial and temporal resolution.

Crystal structures of the human neurokinin 1 receptor in complex with clinically used antagonists

Neurokinins (or tachykinins) are peptides that modulate a wide variety of human physiology through the neurokinin G protein-coupled receptor family, implicated in a diverse array of pathological processes. Here we report high-resolution crystal structures of the human NK₁ receptor (NK₁R) bound to two small-molecule antagonist therapeutics – aprepitant and netupitant and the progenitor antagonist CP-99,994. The structures reveal the detailed interactions between clinically approved antagonists and NK₁R, which induce a distinct receptor conformation resulting in an interhelical hydrogen-bond network that cross-links the extracellular ends of helices V and VI. Furthermore, the high-resolution details of NK₁R bound to netupitant establish a structural rationale for the lack of basal activity in NK₁R. Taken together, these co-structures provide a comprehensive structural basis of NK₁R antagonism and will facilitate the design of new therapeutics targeting the neurokinin receptor family.

Neurokinin receptors are G protein-coupled receptors. Here the authors present three crystal structures of the neurokinin 1 receptor (NK₁R) in complex with small-molecule antagonists including aprepitant and netupitant and observe that these clinically approved compounds induce a conformational change in the receptor.

Calcium-sensing receptor residues with loss- and gain-of-function mutations are located in regions of conformational change and cause signalling bias

The calcium-sensing receptor (CaSR) is a homodimeric G-protein-coupled receptor that signals via intracellular calcium (Ca2+i) mobilisation and phosphorylation of extracellular signal-regulated kinase 1/2 (ERK) to regulate extracellular calcium (Ca2+e) homeostasis. The central importance of the CaSR in Ca2+e homeostasis has been demonstrated by the identification of loss- or gain-of-function CaSR mutations that lead to familial hypocalciuric hypercalcaemia (FHH) or autosomal dominant hypocalcaemia (ADH), respectively. However, the mechanisms determining whether the CaSR signals via Ca2+i or ERK have not been established, and we hypothesised that some CaSR residues, which are the site of both loss- and gain-of-function mutations, may act as molecular switches to direct signalling through these pathways. An analysis of CaSR mutations identified in >300 hypercalcaemic and hypocalcaemic probands revealed five 'disease-switch' residues (Gln27, Asn178, Ser657, Ser820 and Thr828) that are affected by FHH and ADH mutations. Functional expression studies using HEK293 cells showed disease-switch residue mutations to commonly display signalling bias. For example, two FHH-associated mutations (p.Asn178Asp and p.Ser820Ala) impaired Ca2+i signalling without altering ERK phosphorylation. In contrast, an ADH-associated p.Ser657Cys mutation uncoupled signalling by leading to increased Ca2+i mobilization while decreasing ERK phosphorylation. Structural analysis of these five CaSR disease-switch residues together with four reported disease-switch residues revealed these residues to be located at conformationally active regions of the CaSR such as the extracellular dimer interface and transmembrane domain. Thus, our findings indicate that disease-switch residues are located at sites critical for CaSR activation and play a role in mediating signalling bias.

Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations

Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.

Evolutionary action and structural basis of the allosteric switch controlling β2AR functional selectivity

Functional selectivity of G-protein-coupled receptors is believed to originate from ligand-specific conformations that activate only subsets of signaling effectors. In this study, to identify molecular motifs playing important roles in transducing ligand binding into distinct signaling responses, we combined in silico evolutionary lineage analysis and structure-guided site-directed mutagenesis with large-scale functional signaling characterization and non-negative matrix factorization clustering of signaling profiles. Clustering based on the signaling profiles of 28 variants of the β₂-adrenergic receptor reveals three clearly distinct phenotypical clusters, showing selective impairments of either the Gi or βarrestin/endocytosis pathways with no effect on Gs activation. Robustness of the results is confirmed using simulation-based error propagation. The structural changes resulting from functionally biasing mutations centered around the DRY, NPxxY, and PIF motifs, selectively linking these micro-switches to unique signaling profiles. Our data identify different receptor regions that are important for the stabilization of distinct conformations underlying functional selectivity.

Ligand-induced biased signaling is thought to result in part from ligand-specific receptor conformations that cause the engagement of distinct effectors. Here the authors trace and evaluate the impact of mutations of the β2–adrenergic receptor on multiple signaling outputs to provide structural-level insight into the determinants of GPCR functional selectivity.

Chlorophyll-Derivative Modulation of Rhodopsin Signaling Properties through Evolutionarily Conserved Interaction Pathways

Retinal is the light-absorbing chromophore that is responsible for the activation of visual pigments and light-driven ion pumps. Evolutionary changes in the intermolecular interactions of the retinal with specific amino acids allow for adaptation of the spectral characteristics, referred to as spectral tuning. However, it has been proposed that a specific species of dragon fish has bypassed the adaptive evolutionary process of spectral tuning and replaced it with a single evolutionary event: photosensitization of rhodopsin by chlorophyll derivatives. Here, by using a combination of experimental measurements and computational modeling to probe retinal-receptor interactions in rhodopsin, we show how the binding of the chlorophyll derivative, chlorin-e6 (Ce6) in the intracellular domain (ICD) of the receptor allosterically excites G-protein coupled receptor class A (GPCR-A) conserved long-range correlated fluctuations that connect distant parts of the receptor. These long-range correlated motions are associated with regulating the dynamics and intermolecular interactions of specific amino acids in the retinal ligand-binding pocket that have been associated with shifts in the absorbance peak maximum (λ_max) and hence, spectral sensitivity of the visual system. Moreover, the binding of Ce6 affects the overall global properties of the receptor. Specifically, we find that Ce6-induced dynamics alter the thermal stability of rhodopsin by adjusting hydrogen-bonding interactions near the receptor active-site that consequently also influences the intrinsic conformational equilibrium of the receptor. Due to the conservation of the ICD residues amongst different receptors in this class and the fact that all GPCR-A receptors share a common mechanism of activation, it is possible that the allosteric associations excited in rhodopsin with Ce6 binding are a common feature in all class A GPCRs.

Enhancing Structure Prediction and Design of Soluble and Membrane Proteins with Explicit Solvent-Protein Interactions

Solvent molecules interact intimately with proteins and can profoundly regulate their structure and function. However, accurately and efficiently modeling protein solvation effects at the molecular level has been challenging. Here, we present a method that improves the atomic-level modeling of soluble and membrane protein structures and binding by efficiently predicting de novo protein-solvent molecule interactions. The method predicted with unprecedented accuracy buried water molecule positions, solvated protein conformations, and challenging mutational effects on protein binding. When applied to homology modeling, solvent-bound membrane protein structures, pockets, and cavities were recapitulated with near-atomic precision even from distant homologs. Blindly refined atomic-level structures of evolutionary distant G protein-coupled receptors imply strikingly different functional roles of buried solvent between receptor classes. The method should prove useful for refining low-resolution protein structures, accurately modeling drug-binding sites in structurally uncharacterized receptors, and designing solvent-mediated protein catalysis, recognition, ligand binding, and membrane protein signaling.

Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis

Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.

G Protein-Coupled Receptors Contain Two Conserved Packing Clusters

G protein-coupled receptors (GPCRs) have evolved a seven-transmembrane helix framework that is responsive to a wide range of extracellular signals. An analysis of the interior packing of family A GPCR crystal structures reveals two clusters of highly packed residues that facilitate tight transmembrane helix association. These clusters are centered on amino acid positions 2.47 and 4.53, which are highly conserved as alanine and serine, respectively. Ala2.47 mediates the interaction between helices H1 and H2, while Ser4.53 mediates the interaction between helices H3 and H4. The helical interfaces outside of these clusters are lined with residues that are more loosely packed, a structural feature that facilitates motion of helices H5, H6, and H7, which is required for receptor activation. Mutation of the conserved small side chain at position 4.53 within packing cluster 2 is shown to disrupt the structure of the visual receptor rhodopsin, whereas sites in packing cluster 1 (e.g., positions 1.46 and 2.47) are more tolerant to mutation but affect the overall stability of the protein. These findings reveal a common structural scaffold of GPCRs that is important for receptor folding and activation.

Evaluation of a melanocortin-4 receptor (MC4R) agonist (Setmelanotide) in MC4R deficiency

Objective

Pro-opiomelanocortin (POMC)-derived peptides act on neurons expressing the Melanocortin 4 receptor (MC4R) to reduce body weight. Setmelanotide is a highly potent MC4R agonist that leads to weight loss in diet-induced obese animals and in obese individuals with complete POMC deficiency. While POMC deficiency is very rare, 1–5% of severely obese individuals harbor heterozygous mutations in MC4R. We sought to assess the efficacy of Setmelanotide in human MC4R deficiency.

Methods

We studied the effects of Setmelanotide on mutant MC4Rs in cells and the weight loss response to Setmelanotide administration in rodent studies and a human clinical trial. We annotated the functional status of 369 published MC4R variants.

Results

In cells, we showed that Setmelanotide is significantly more potent at MC4R than the endogenous ligand alpha-melanocyte stimulating hormone and can disproportionally rescue signaling by a subset of severely impaired MC4R mutants. Wild-type rodents appear more sensitive to Setmelanotide when compared to MC4R heterozygous deficient mice, while MC4R knockout mice fail to respond. In a 28-day Phase 1b clinical trial, Setmelanotide led to weight loss in obese MC4R variant carriers. Patients with POMC defects upstream of MC4R show significantly more weight loss with Setmelanotide than MC4R deficient patients or obese controls.

Conclusions

Setmelanotide led to weight loss in obese people with MC4R deficiency; however, further studies are justified to establish whether Setmelanotide can elicit clinically meaningful weight loss in a subset of the MC4R deficient obese population.

Evolution of nonspectral rhodopsin function at high altitudes

High-altitude environments present a range of biochemical and physiological challenges for organisms through decreases in oxygen, pressure, and temperature relative to lowland habitats. Protein-level adaptations to hypoxic high-altitude conditions have been identified in multiple terrestrial endotherms; however, comparable adaptations in aquatic ectotherms, such as fishes, have not been as extensively characterized. In enzyme proteins, cold adaptation is attained through functional trade-offs between stability and activity, often mediated by substitutions outside the active site. Little is known whether signaling proteins [e.g., G protein-coupled receptors (GPCRs)] exhibit natural variation in response to cold temperatures. Rhodopsin (RH1), the temperature-sensitive visual pigment mediating dim-light vision, offers an opportunity to enhance our understanding of thermal adaptation in a model GPCR. Here, we investigate the evolution of rhodopsin function in an Andean mountain catfish system spanning a range of elevations. Using molecular evolutionary analyses and site-directed mutagenesis experiments, we provide evidence for cold adaptation in RH1. We find that unique amino acid substitutions occur at sites under positive selection in high-altitude catfishes, located at opposite ends of the RH1 intramolecular hydrogen-bonding network. Natural high-altitude variants introduced into these sites via mutagenesis have limited effects on spectral tuning, yet decrease the stability of dark-state and light-activated rhodopsin, accelerating the decay of ligand-bound forms. As found in cold-adapted enzymes, this phenotype likely compensates for a cold-induced decrease in kinetic rates-properties of rhodopsin that mediate rod sensitivity and visual performance. Our results support a role for natural variation in enhancing the performance of GPCRs in response to cold temperatures.

Vibrational resonance, allostery, and activation in rhodopsin-like G protein-coupled receptors

G protein-coupled receptors are a large family of membrane proteins activated by a variety of structurally diverse ligands making them highly adaptable signaling molecules. Despite recent advances in the structural biology of this protein family, the mechanism by which ligands induce allosteric changes in protein structure and dynamics for its signaling function remains a mystery. Here, we propose the use of terahertz spectroscopy combined with molecular dynamics simulation and protein evolutionary network modeling to address the mechanism of activation by directly probing the concerted fluctuations of retinal ligand and transmembrane helices in rhodopsin. This approach allows us to examine the role of conformational heterogeneity in the selection and stabilization of specific signaling pathways in the photo-activation of the receptor. We demonstrate that ligand-induced shifts in the conformational equilibrium prompt vibrational resonances in the protein structure that link the dynamics of conserved interactions with fluctuations of the active-state ligand. The connection of vibrational modes creates an allosteric association of coupled fluctuations that forms a coherent signaling pathway from the receptor ligand-binding pocket to the G-protein activation region. Our evolutionary analysis of rhodopsin-like GPCRs suggest that specific allosteric sites play a pivotal role in activating structural fluctuations that allosterically modulate functional signals.

Can Specific Protein-Lipid Interactions Stabilize an Active State of the Beta 2 Adrenergic Receptor?

G-protein-coupled receptors are eukaryotic membrane proteins with broad biological and pharmacological relevance. Like all membrane-embedded proteins, their location and orientation are influenced by lipids, which can also impact protein function via specific interactions. Extensive simulations totaling 0.25 ms reveal a process in which phospholipids from the membrane's cytosolic leaflet enter the empty G-protein binding site of an activated β2 adrenergic receptor and form salt-bridge interactions that inhibit ionic lock formation and prolong active-state residency. Simulations of the receptor embedded in an anionic membrane show increased lipid binding, providing a molecular mechanism for the experimental observation that anionic lipids can enhance receptor activity. Conservation of the arginine component of the ionic lock among Rhodopsin-like G-protein-coupled receptors suggests that intracellular lipid ingression between receptor helices H6 and H7 may be a general mechanism for active-state stabilization.

Structure-Activity Analysis of Biased Agonism at the Human Adenosine A3 Receptor

Biased agonism at G protein-coupled receptors (GPCRs) has significant implications for current drug discovery, but molecular determinants that govern ligand bias remain largely unknown. The adenosine A3 GPCR (A3AR) is a potential therapeutic target for various conditions, including cancer, inflammation, and ischemia, but for which biased agonism remains largely unexplored. We now report the generation of bias "fingerprints" for prototypical ribose containing A3AR agonists and rigidified (N)-methanocarba 5'-N-methyluronamide nucleoside derivatives with regard to their ability to mediate different signaling pathways. Relative to the reference prototypical agonist IB-MECA, (N)-methanocarba 5'-N-methyluronamide nucleoside derivatives with significant N(6) or C2 modifications, including elongated aryl-ethynyl groups, exhibited biased agonism. Significant positive correlation was observed between the C2 substituent length (in Å) and bias toward cell survival. Molecular modeling suggests that extended C2 substituents on (N)-methanocarba 5'-N-methyluronamide nucleosides promote a progressive outward shift of the A3AR transmembrane domain 2, which may contribute to the subset of A3AR conformations stabilized on biased agonist binding.

Biased Gs versus Gq proteins and β-arrestin signaling in the NK1 receptor determined by interactions in the water hydrogen bond network

X-ray structures, molecular dynamics simulations, and mutational analysis have previously indicated that an extended water hydrogen bond network between trans-membranes I-III, VI, and VII constitutes an allosteric interface essential for stabilizing different active and inactive helical constellations during the seven-trans-membrane receptor activation. The neurokinin-1 receptor signals efficiently through Gq, Gs, and β-arrestin when stimulated by substance P, but it lacks any sign of constitutive activity. In the water hydrogen bond network the neurokinin-1 has a unique Glu residue instead of the highly conserved AspII:10 (2.50). Here, we find that this GluII:10 occupies the space of a putative allosteric modulating Na(+) ion and makes direct inter-helical interactions in particular with SerIII:15 (3.39) and AsnVII:16 (7.49) of the NPXXY motif. Mutational changes in the interface between GluII:10 and AsnVII:16 created receptors that selectively signaled through the following: 1) Gq only; 2) β-arrestin only; and 3) Gq and β-arrestin but not through Gs. Interestingly, increased constitutive Gs but not Gq signaling was observed by Ala substitution of four out of the six core polar residues of the network, in particular SerIII:15. Three residues were essential for all three signaling pathways, i.e. the water-gating micro-switch residues TrpVI:13 (6.48) of the CWXP motif and TyrVII:20 (7.53) of the NPXXY motif plus the totally conserved AsnI:18 (1.50) stabilizing the kink in trans-membrane VII. It is concluded that the interface between position II:10 (2.50), III:15 (3.39), and VII:16 (7.49) in the center of the water hydrogen bond network constitutes a focal point for fine-tuning seven trans-membrane receptor conformations activating different signal transduction pathways.

Quantifying conformational changes in GPCRs: glimpse of a common functional mechanism

Background

G-protein-coupled receptors (GPCRs) are important drug targets and a better understanding of their molecular mechanisms would be desirable. The crystallization rate of GPCRs has accelerated in recent years as techniques have become more sophisticated, particularly with respect to Class A GPCRs interacting with G-proteins. These developments have made it possible for a quantitative analysis of GPCR geometrical features and binding-site conformations, including a statistical comparison between Class A GPCRs in active (agonist-bound) and inactive (antagonist-bound) states.

Results

Here we implement algorithms for the analysis of interhelical angles, distances, interactions and binding-site volumes in the transmembrane domains of 25 Class A GPCRs (7 active and 18 inactive). Two interhelical angles change in a statistically significant way between average inactive and active states: TM3-TM6 (by -9°) and TM6-TM7 (by +12°). A third interhelical angle: TM5-TM6 shows a trend, changing by -9°. In the transition from inactive to active states, average van der Waals interactions between TM3 and TM7 significantly increase as the average distance between them decreases by >2 Å. Average H-bonding between TM3 and TM6 decreases but is seemingly compensated by an increase in H-bonding between TM5 and TM6. In five Class A GPCRs, crystallized in both active and inactive states, increased H-bonding of agonists to TM6 and TM7, relative to antagonists, is observed. These protein-agonist interactions likely favour a change in the TM6-TM7 angle, which creates a narrowing in the binding pocket of activated receptors and an average ~200 ų reduction in volume.

Conclusions

In terms of similar conformational changes and agonist binding pattern, Class A GPCRs appear to share a common mechanism of activation, which can be exploited in future drug development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0567-3) contains supplementary material, which is available to authorized users.

Retinal ligand mobility explains internal hydration and reconciles active rhodopsin structures

Rhodopsin, the mammalian dim-light receptor, is one of the best-characterized G-protein-coupled receptors, a pharmaceutically important class of membrane proteins that has garnered a great deal of attention because of the recent availability of structural information. Yet the mechanism of rhodopsin activation is not fully understood. Here, we use microsecond-scale all-atom molecular dynamics simulations, validated by solid-state (2)H nuclear magnetic resonance spectroscopy, to understand the transition between the dark and metarhodopsin I (Meta I) states. Our analysis of these simulations reveals striking differences in ligand flexibility between the two states. Retinal is much more dynamic in Meta I, adopting an elongated conformation similar to that seen in the recent activelike crystal structures. Surprisingly, this elongation corresponds to both a dramatic influx of bulk water into the hydrophobic core of the protein and a concerted transition in the highly conserved Trp265(6.48) residue. In addition, enhanced ligand flexibility upon light activation provides an explanation for the different retinal orientations observed in X-ray crystal structures of active rhodopsin.

Probing the metabotropic glutamate receptor 5 (mGlu₅) positive allosteric modulator (PAM) binding pocket: discovery of point mutations that engender a "molecular switch" in PAM pharmacology

Positive allosteric modulation of metabotropic glutamate receptor subtype 5 (mGlu₅) is a promising novel approach for the treatment of schizophrenia and cognitive disorders. Allosteric binding sites are topographically distinct from the endogenous ligand (orthosteric) binding site, allowing for co-occupation of a single receptor with the endogenous ligand and an allosteric modulator. Negative allosteric modulators (NAMs) inhibit and positive allosteric modulators (PAMs) enhance the affinity and/or efficacy of the orthosteric agonist. The molecular determinants that govern mGlu₅ modulator affinity versus cooperativity are not well understood. Focusing on the modulators based on the acetylene scaffold, we sought to determine the molecular interactions that contribute to PAM versus NAM pharmacology. Generation of a comparative model of the transmembrane-spanning region of mGlu₅ served as a tool to predict and interpret the impact of mutations in this region. Application of an operational model of allosterism allowed for determination of PAM and NAM affinity estimates at receptor constructs that possessed no detectable radioligand binding as well as delineation of effects on affinity versus cooperativity. Novel mutations within the transmembrane domain (TM) regions were identified that had differential effects on acetylene PAMs versus 2-methyl-6-(phenylethynyl)-pyridine, a prototypical NAM. Three conserved amino acids (Y658, T780, and S808) and two nonconserved residues (P654 and A809) were identified as key determinants of PAM activity. Interestingly, we identified two point mutations in TMs 6 and 7 that, when mutated, engender a mode switch in the pharmacology of certain PAMs.

Allosteric transition: a comparison of two models

Introduction

Two recent models are in use for analysis of allosteric drug action at receptor sites remote from orthosteric binding sites. One is an allosteric two-state mechanical model derived in 2000 by David Hall. The other is an extended operational model developed in 2007 by Arthur Christopoulos’s group. The models are valid in pharmacology, enzymology, transportology as well as several other fields of biology involving allosteric concentration effects.

Results

I show here that Hall’s model for interactions between an orthoster, an alloster, and a receptive unit is the best choice of model both for simulation and analysis of allosteric concentration-responses at equilibrium or steady-state.

Conclusions

As detailed knowledge of receptors systems becomes available, systems with several pathways and states and/ or more than two binding sites should be analysed by extended forms of the Hall model rather than for instance a Hill type exponentiation of terms as introduced in non-mechanistic (operational) model approaches; yielding semi-quantitative estimates of actual system parameters based on Hill’s unlikely simultaneity model for G protein-coupled receptors.

Critical hydrogen bond formation for activation of the angiotensin II type 1 receptor

G protein-coupled receptors contain selectively important residues that play central roles in the conformational changes that occur during receptor activation. Asparagine 111 (N111(3.35)) is such a residue within the angiotensin II type 1 (AT(1)) receptor. Substitution of N111(3.35) for glycine leads to a constitutively active receptor, whereas substitution for tryptophan leads to an inactivable receptor. Here, we analyzed the AT(1) receptor and two mutants (N111G and N111W) by molecular dynamics simulations, which revealed a novel molecular switch involving the strictly conserved residue D74(2.50). Indeed, D74(2.50) forms a stable hydrogen bond (H-bond) with the residue in position 111(3.35) in the wild-type and the inactivable receptor. However, in the constitutively active mutant N111G-AT(1) receptor, residue D74 is reoriented to form a new H-bond with another strictly conserved residue, N46(1.50). When expressed in HEK293 cells, the mutant N46G-AT(1) receptor was poorly activable, although it retained a high binding affinity. Interestingly, the mutant N46G/N111G-AT(1) receptor was also inactivable. Molecular dynamics simulations also revealed the presence of a cluster of hydrophobic residues from transmembrane domains 2, 3, and 7 that appears to stabilize the inactive form of the receptor. Whereas this hydrophobic cluster and the H-bond between D74(2.50) and W111(3.35) are more stable in the inactivable N111W-AT(1) receptor, the mutant N111W/F77A-AT(1) receptor, designed to weaken the hydrophobic core, showed significant agonist-induced signaling. These results support the potential for the formation of an H-bond between residues D74(2.50) and N46(1.50) in the activation of the AT(1) receptor.

Hybrid molecular mechanics/coarse-grained simulations for structural prediction of G-protein coupled receptor/ligand complexes

Understanding how ligands bind to G-protein coupled receptors (GPCRs) provides insights into a myriad of cell processes and is crucial for drug development. Here we extend a hybrid molecular mechanics/coarse-grained (MM/CG) approach applied previously to enzymes to GPCR/ligand complexes. The accuracy of this method for structural predictions is established by comparison with recent atomistic molecular dynamics simulations on the human β2 adrenergic receptor, a member of the GPCRs superfamily. The results obtained with the MM/CG methodology show a good agreement with previous all-atom classical dynamics simulations, in particular in the structural description of the ligand binding site. This approach could be used for high-throughput predictions of ligand poses in a variety of GPCRs.

Modulation of a pre-existing conformational equilibrium tunes adenylate kinase activity

Structural plasticity is often required for distinct microscopic steps during enzymatic reaction cycles. Adenylate kinase from Escherichia coli (AK(eco)) populates two major conformations in solution; the open (inactive) and closed (active) state, and the overall turnover rate is inversely proportional to the lifetime of the active conformation. Therefore, structural plasticity is intimately coupled to enzymatic turnover in AK(eco). Here, we probe the open to closed conformational equilibrium in the absence of bound substrate with NMR spectroscopy and molecular dynamics simulations. The conformational equilibrium in absence of substrate and, in turn, the turnover number can be modulated with mutational- and osmolyte-driven perturbations. Removal of one hydrogen bond between the ATP and AMP binding subdomains results in a population shift toward the open conformation and a resulting increase of k(cat). Addition of the osmolyte TMAO to AK(eco) results in population shift toward the closed conformation and a significant reduction of k(cat). The Michaelis constants (K(M)) scale with the change in k(cat), which follows from the influence of the population of the closed conformation for substrate binding affinity. Hence, k(cat) and K(M) are mutually dependent, and in the case of AK(eco), any perturbation that modulates k(cat) is mirrored with a proportional response in K(M). Thus, our results demonstrate that the equilibrium constant of a pre-existing conformational equilibrium directly affects enzymatic catalysis. From an evolutionary perspective, our findings suggest that, for AK(eco), there exists ample flexibility to obtain a specificity constant (k(cat)/K(M)) that commensurate with the exerted cellular selective pressure.

The arginine of the DRY motif in transmembrane segment III functions as a balancing micro-switch in the activation of the β2-adrenergic receptor

Recent high resolution x-ray structures of the β2-adrenergic receptor confirmed a close salt-bridge interaction between the suspected micro-switch residue ArgIII:26 (Arg3.50) and the neighboring AspIII:25 (Asp3.49). However, neither the expected "ionic lock" interactions between ArgIII:26 and GluVI:-06 (Glu6.30) in the inactive conformation nor the interaction with TyrV:24 (Tyr5.58) in the active conformation were observed in the x-ray structures. Here we find through molecular dynamics simulations, after removal of the stabilizing T4 lysozyme, that the expected salt bridge between ArgIII:26 and GluVI:-06 does form relatively easily in the inactive receptor conformation. Moreover, mutational analysis of GluVI:-06 in TM-VI and the neighboring AspIII:25 in TM-III demonstrated that these two residues do function as locks for the inactive receptor conformation as we observed increased G(s) signaling, arrestin mobilization, and internalization upon alanine substitutions. Conversely, TyrV:24 appears to play a role in stabilizing the active receptor conformation as loss of function of G(s) signaling, arrestin mobilization, and receptor internalization was observed upon alanine substitution of TyrV:24. The loss of function of the TyrV:24 mutant could partly be rescued by alanine substitution of either AspIII:25 or GluVI:-06 in the double mutants. Surprisingly, removal of the side chain of the ArgIII:26 micro-switch itself had no effect on G(s) signaling and internalization and only reduced arrestin mobilization slightly. It is suggested that ArgIII:26 is equally important for stabilizing the inactive and the active conformation through interaction with key residues in TM-III, -V, and -VI, but that the ArgIII:26 micro-switch residue itself apparently is not essential for the actual G protein activation.

Evolutionary Analysis of Functional Divergence among Chemokine Receptors, Decoy Receptors, and Viral Receptors

Chemokine receptors (CKRs) function in the inflammatory response and in vertebrate homeostasis. Decoy and viral receptors are two types of CKR homologs with modified functions from those of the typical CKRs. The decoy receptors are able to bind ligands without signaling. On the other hand, the viral receptors show constitutive signaling without ligands. We examined the sites related to the functional difference. At first, the decoy and viral receptors were each classified into five groups, based on the molecular phylogenetic analysis. A multiple amino acid sequence alignment between each group and the CKRs was then constructed. The difference in the amino acid composition between the group and the CKRs was evaluated as the Kullback-Leibler (KL) information value at each alignment site. The KL information value is considered to reflect the difference in the functional constraints at the site. The sites with the top 5% of KL information values were selected and mapped on the structure of a CKR. The comparisons with decoy receptor groups revealed that the detected sites were biased on the intracellular side. In contrast, the sites detected from the comparisons with viral receptor groups were found on both the extracellular and intracellular sides. More sites were found in the ligand binding pocket in the analyses of the viral receptor groups, as compared to the decoy receptor groups. Some of the detected sites were located in the GPCR motifs. For example, the DRY motif of the decoy receptors was often degraded, although the motif of the viral receptors was basically conserved. The observations for the viral receptor groups suggested that the constraints in the pocket region are loose and that the sites on the intracellular side are different from those for the decoy receptors, which may be related to the constitutive signaling activity of the viral receptors.

Pharmacological modulation of chemokine receptor function

G protein-coupled chemokine receptors and their peptidergic ligands are interesting therapeutic targets due to their involvement in various immune-related diseases, including rheumatoid arthritis, multiple sclerosis, inflammatory bowel disease, chronic obstructive pulmonary disease, HIV-1 infection and cancer. To tackle these diseases, a lot of effort has been focused on discovery and development of small-molecule chemokine receptor antagonists. This has been rewarded by the market approval of two novel chemokine receptor inhibitors, AMD3100 (CXCR4) and Maraviroc (CCR5) for stem cell mobilization and treatment of HIV-1 infection respectively. The recent GPCR crystal structures together with mutagenesis and pharmacological studies have aided in understanding how small-molecule ligands interact with chemokine receptors. Many of these ligands display behaviour deviating from simple competition and do not interact with the chemokine binding site, providing evidence for an allosteric mode of action. This review aims to give an overview of the evidence supporting modulation of this intriguing receptor family by a range of ligands, including small molecules, peptides and antibodies. Moreover, the computer-assisted modelling of chemokine receptor-ligand interactions is discussed in view of GPCR crystal structures. Finally, the implications of concepts such as functional selectivity and chemokine receptor dimerization are considered.

Structural basis of activation of bitter taste receptor T2R1 and comparison with Class A G-protein-coupled receptors (GPCRs)

The human bitter taste receptors (T2Rs) are non-Class A members of the G-protein-coupled receptor (GPCR) superfamily, with very limited structural information. Amino acid sequence analysis reveals that most of the important motifs present in the transmembrane helices (TM1-TM7) of the well studied Class A GPCRs are absent in T2Rs, raising fundamental questions regarding the mechanisms of activation and how T2Rs recognize bitter ligands with diverse chemical structures. In this study, the bitter receptor T2R1 was used to systematically investigate the role of 15 transmembrane amino acids in T2Rs, including 13 highly conserved residues, by amino acid replacements guided by molecular modeling. Functional analysis of the mutants by calcium imaging analysis revealed that replacement of Asn-66(2.65) and the highly conserved Asn-24(1.50) resulted in greater than 90% loss of agonist-induced signaling. Our results show that Asn-24(1.50) plays a crucial role in receptor activation by mediating an hydrogen bond network connecting TM1-TM2-TM7, whereas Asn-66(2.65) is essential for binding to the agonist dextromethorphan. The interhelical hydrogen bond between Asn-24(1.50) and Arg-55(2.54) restrains T2R receptor activity because loss of this bond in I27A and R55A mutants results in hyperactive receptor. The conserved amino acids Leu-197(5.50), Ser-200(5.53), and Leu-201(5.54) form a putative LXXSL motif which performs predominantly a structural role by stabilizing the helical conformation of TM5 at the cytoplasmic end. This study provides for the first time mechanistic insights into the roles of the conserved transmembrane residues in T2Rs and allows comparison of the activation mechanisms of T2Rs with the Class A GPCRs.

From molecular details of the interplay between transmembrane helices of the thyrotropin receptor to general aspects of signal transduction in family a G-protein-coupled receptors (GPCRs)

Transmembrane helices (TMHs) 5 and 6 are known to be important for signal transduction by G-protein-coupled receptors (GPCRs). Our aim was to characterize the interface between TMH5 and TMH6 of the thyrotropin receptor (TSHR) to gain molecular insights into aspects of signal transduction and regulation. A proline at TMH5 position 5.50 is highly conserved in family A GPCRs and causes a twist in the helix structure. Mutation of the TSHR-specific alanine (Ala-593⁵·⁵⁰) at this position to proline resulted in a 20-fold reduction of cell surface expression. This indicates that TMH5 in the TSHR might have a conformation different from most other family A GPCRs by forming a regular α-helix. Furthermore, linking our own and previous data from directed mutagenesis with structural information led to suggestions of distinct pairs of interacting residues between TMH5 and TMH6 that are responsible for stabilizing either the basal or the active state. Our insights suggest that the inactive state conformation is constrained by a core set of polar interactions among TMHs 2, 3, 6, and 7 and in contrast that the active state conformation is stabilized mainly by non-polar interactions between TMHs 5 and 6. Our findings might be relevant for all family A GPCRs as supported by a statistical analysis of residue properties between the TMHs of a vast number of GPCR sequences.

Evidence from Chlamydomonas on the photoactivation of rhodopsins without isomerization of their chromophore

Attachment of retinal to opsin forms the chromophore N-retinylidene, which isomerizes during photoactivation of rhodopsins. To test whether isomerization is crucial, custom-tailored chromophores lacking the β-ionone ring and any isomerizable bonds were incorporated in vivo into the opsin of a blind mutant of the eukaryote Chlamydomonas reinhardtii. The analogs restored phototaxis with the anticipated action spectra, ruling out the need for isomerization in photoactivation. To further elucidate photoactivation, responses to chromophores formed from naphthalene aldehydes were studied. The resulting action spectral shifts suggest that charge separation within the excited chromophore leads to electric field-induced polarization of nearby amino acid residues and altered hydrogen bonding. This redistribution of charge facilitates the reported multiple bond rotations and protein rearrangements of rhodopsin activation. These results provide insight into the activation of rhodopsins and related GPCRs.

Predicting novel binding modes of agonists to β adrenergic receptors using all-atom molecular dynamics simulations

Understanding the binding mode of agonists to adrenergic receptors is crucial to enabling improved rational design of new therapeutic agents. However, so far the high conformational flexibility of G protein-coupled receptors has been an obstacle to obtaining structural information on agonist binding at atomic resolution. In this study, we report microsecond classical molecular dynamics simulations of β(1) and β(2) adrenergic receptors bound to the full agonist isoprenaline and in their unliganded form. These simulations show a novel agonist binding mode that differs from the one found for antagonists in the crystal structures and from the docking poses reported by in silico docking studies performed on rigid receptors. Internal water molecules contribute to the stabilization of novel interactions between ligand and receptor, both at the interface of helices V and VI with the catechol group of isoprenaline as well as at the interface of helices III and VII with the ethanolamine moiety of the ligand. Despite the fact that the characteristic N-C-C-OH motif is identical in the co-crystallized ligands and in the full agonist isoprenaline, the interaction network between this group and the anchor site formed by Asp(3.32) and Asn(7.39) is substantially different between agonists and inverse agonists/antagonists due to two water molecules that enter the cavity and contribute to the stabilization of a novel network of interactions. These new binding poses, together with observed conformational changes in the extracellular loops, suggest possible determinants of receptor specificity.

The minor binding pocket: a major player in 7TM receptor activation

From the deep part of the main ligand-binding crevice, a minor, often shallower pocket extends between the extracellular ends of transmembrane domains (TM)-I, II, III and VII of 7TM receptors. This minor binding pocket is defined by a highly conserved kink in TM-II that is induced by a proline residue located in one of two adjacent positions. Here we argue that this minor binding pocket is important for receptor activation. Functional coupling of the receptors seems to be mediated through the hydrogen bond network located between the intracellular segments of these TMs, with the allosteric interface between TM-II and TM-VII being of particular significance. Importantly, the minor binding pocket, especially the proline-kink in TM-II, is involved in G protein versus arrestin pathway-biased signaling, for example in the angiotensin AT1 system. Consequently, this pocket could be specifically targeted in the development of functionally biased drugs.