CoChat
Sign Up
Cited by 14Receptor Mechanisms and SignalingProtein Structure and DynamicsChemical Synthesis and Analysis

The DOI matches its Crossref record; 3 of the 3 references it cites were found in a registry or an index.

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

Open access
Chen, K.-Y. M., Lai, J. K., Rudden, L. S. P., Wang, J., Russell, A. M., Conners, K., Rutter, M. E., Condon, B., Tung, F., Kodandapani, L., Chau, B., Zhao, X., Benach, J., Baker, K., Hembre, E. J., Barth, P.·2025·Nature Chemistry
DOI 10.1038/s41557-024-01719-2PMC PMC11882447

Published version of a preprint

Open this paper in CoChat — free accountOpen PDF

Abstract

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.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41557-024-01719-2

Cites

3 of 3 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1 Rational design of signalling activity through de novo solvent-mediated interaction networks. a, Dynamic remodelling of the solvent-mediated allosteric network during the activation of the μ-opioid receptor (mOR). The positions of water molecules and hydrogen bonds were extracted from the inactive (Protein Data Bank (PDB): 4DKL) and active (PDB: 5C1M) mOR structures. b, Schematic representation of a water-mediated hydrogen bond network bridging four helices designed using the SPaDES software. De novo amino acid sequences, conformations and water positions were searched concurrently for optimal water-mediated hydrogen bond connectivity at the helical interface. c, Snapshot of a low-energy water-mediated hydrogen bond network ensemble in the core of a GPCR structure. d–f, Computational design strategies for enhancing protein activity and conformational stability following three criteria: shifts in conformational energies (d), increased water-mediated dynamic contacts in the active state (e) or weakened ion-mediated locks in the inactive state (f).
Fig. 1 Rational design of signalling activity through de novo solvent-mediated interaction networks. a, Dynamic remodelling of the solvent-mediated allosteric network during the activation of the μ-opioid receptor (mOR). The positions of water molecules and hydrogen bonds were extracted from the inactive (Protein Data Bank (PDB): 4DKL) and active (PDB: 5C1M) mOR structures. b, Schematic representation of a water-mediated hydrogen bond network bridging four helices designed using the SPaDES software. De novo amino acid sequences, conformations and water positions were searched concurrently for optimal water-mediated hydrogen bond connectivity at the helical interface. c, Snapshot of a low-energy water-mediated hydrogen bond network ensemble in the core of a GPCR structure. d–f, Computational design strategies for enhancing protein activity and conformational stability following three criteria: shifts in conformational energies (d), increased water-mediated dynamic contacts in the active state (e) or weakened ion-mediated locks in the inactive state (f).CC BY · PubMed Central
Fig. 2 Computational-designed GPCRs with a wide range of water-mediated interaction networks at switchable TMH interfaces. a, Backbone representation of a GPCR scaffold highlighting static (dark grey) and switchable (light grey) TM helices. b, Relationship between the total number of buried water molecules mediating polar contacts and those at static–switchable TMH interfaces. Data for the native A2AR (light grey) are included for comparison. Designs with a higher and lower number of water-mediated static–switchable hydrogen bonds than A2AR are classified as Hyd_high and Hyd_low, respectively. c,d, Hydrogen bond network in the TMH core regions of the Hyd_low2 (c) and Hyd_high4 (d) structures. The positions of the residues highlighted as stick models are indicated by Ballesteros–Weinstein notation. e, Adenosine-induced activity of the designed receptor variants (as a percentage of adenosine-induced WT A2AR activity) as a function of the number of solvent-mediated interactions at static–switchable TMH interfaces. A baseline for the native A2AR (light grey) is included for comparison. The data are presented as the mean ± standard error of the mean (s.e.m.) for n = 3 independent experiments. Source data
Fig. 2 Computational-designed GPCRs with a wide range of water-mediated interaction networks at switchable TMH interfaces. a, Backbone representation of a GPCR scaffold highlighting static (dark grey) and switchable (light grey) TM helices. b, Relationship between the total number of buried water molecules mediating polar contacts and those at static–switchable TMH interfaces. Data for the native A2AR (light grey) are included for comparison. Designs with a higher and lower number of water-mediated static–switchable hydrogen bonds than A2AR are classified as Hyd_high and Hyd_low, respectively. c,d, Hydrogen bond network in the TMH core regions of the Hyd_low2 (c) and Hyd_high4 (d) structures. The positions of the residues highlighted as stick models are indicated by Ballesteros–Weinstein notation. e, Adenosine-induced activity of the designed receptor variants (as a percentage of adenosine-induced WT A2AR activity) as a function of the number of solvent-mediated interactions at static–switchable TMH interfaces. A baseline for the native A2AR (light grey) is included for comparison. The data are presented as the mean ± standard error of the mean (s.e.m.) for n = 3 independent experiments. Source dataCC BY · PubMed Central
Fig. 3 Designed receptor stability and function strongly depend on the density of solvent-mediated interactions at switchable TMH interfaces. a, Constitutive activity (as a percentage of maximal agonist-induced WT A2AR activity) as a function of the calculated conformational energy difference between the active and inactive states relative to A2AR. The reported energy difference was calculated from the average of the ten lowest-energy inactive- and active-state structures (Supplementary Table 3). The best linear fit to the data is shown as a dotted line (correlation coefficient r = 0.91). b, Measured constitutive activity of the designed receptor variants (as a percentage of adenosine-induced A2AR activity) as a function of the number of solvent-mediated interactions at static–switchable TMH interfaces. c,d, Apparent stability of Hyd_low (c) and Hyd_high (d) receptors compared with A2AR, measured as the functional half-life of adenosine-bound receptor variants incubated at 37 °C. The values for all of the designs are given in Supplementary Table 4. In b–d, the data are presented as the mean ± s.e.m. for n = 3 independent experiments. Source data
Fig. 3 Designed receptor stability and function strongly depend on the density of solvent-mediated interactions at switchable TMH interfaces. a, Constitutive activity (as a percentage of maximal agonist-induced WT A2AR activity) as a function of the calculated conformational energy difference between the active and inactive states relative to A2AR. The reported energy difference was calculated from the average of the ten lowest-energy inactive- and active-state structures (Supplementary Table 3). The best linear fit to the data is shown as a dotted line (correlation coefficient r = 0.91). b, Measured constitutive activity of the designed receptor variants (as a percentage of adenosine-induced A2AR activity) as a function of the number of solvent-mediated interactions at static–switchable TMH interfaces. c,d, Apparent stability of Hyd_low (c) and Hyd_high (d) receptors compared with A2AR, measured as the functional half-life of adenosine-bound receptor variants incubated at 37 °C. The values for all of the designs are given in Supplementary Table 4. In b–d, the data are presented as the mean ± s.e.m. for n = 3 independent experiments. Source dataCC BY · PubMed Central
Fig. 4 Polar interaction networks connecting static and switchable TMHs in the Hyd_low1 design and A2AR. a, Cartoon representation of the Hyd_low1 receptor active-state model, where the adenosine ligand and side chains in the vicinity of the designed mutations are highlighted as stick models. b,c, Enlarged views of the designed region in Hyd_low1 (b) and A2AR (c). Polar contacts are represented by red dotted lines. Residues are numbered using the Ballesteros–Weinstein notation.
Fig. 4 Polar interaction networks connecting static and switchable TMHs in the Hyd_low1 design and A2AR. a, Cartoon representation of the Hyd_low1 receptor active-state model, where the adenosine ligand and side chains in the vicinity of the designed mutations are highlighted as stick models. b,c, Enlarged views of the designed region in Hyd_low1 (b) and A2AR (c). Polar contacts are represented by red dotted lines. Residues are numbered using the Ballesteros–Weinstein notation.CC BY · PubMed Central
Fig. 5 The designed Hyd_high7 structure adopts an unforeseen activated conformation. a, Backbone representation of the designed Hyd_high7 structure bound to the ligand agonist CGS21680 and N-terminally fused to a T4Lcp variant. b, Cytoplasmic view of the following superimposed structures: antagonist-bound A2AR (PDB: 3EML, yellow), agonist-bound, computationally designed Hyd_high7 (green), Gt peptide-bound opsin active state (PDB: 3DQB, blue), and agonist- and mini-Gs-bound A2AR thermostabilized by alanine scanning (PDB: 5G53, magenta). The structures are ranked according to the magnitude of the conformational changes towards the fully active state (back arrow). c–e, Conformations of key conserved microswitches involved in receptor activation (in addition to the reference structures described in panel b, the following structure is also represented: agonist-bound A2AR thermostabilized by alanine scanning (PDB: 2YDO, grey)): the TMH3–TMH5–TMH6 interface around the W6.48 toggle switch in the receptor TM core (c), the TM3–TMH5–TMH6 interface around the intracellular DR3.50Y motif (d) and the TMH6–TMH7 interface around the intracellular N7.49PXXY7.53 motif (e). Conformational changes between inactive, partially active and fully active states are indicated by distinct black arrows. f–k, Comparison of the conformation of the designed (black) and neighbouring native (grey) residue microswitches in the agonist-bound designed Hyd_high7 structure (f) and the agonist- and Gs-bound Hyd_high7 design model (g) as well as the TMH5–TMH6 core interface (h), the TMH6–TMH7 intracellular interface with bound Gs (magenta) (i), the TMH1–TMH2–TMH3 core interface at the Na+ ion-conserved binding site (j) and the TMH3–TMH6–TMH7 intracellular interface (k).
Fig. 5 The designed Hyd_high7 structure adopts an unforeseen activated conformation. a, Backbone representation of the designed Hyd_high7 structure bound to the ligand agonist CGS21680 and N-terminally fused to a T4Lcp variant. b, Cytoplasmic view of the following superimposed structures: antagonist-bound A2AR (PDB: 3EML, yellow), agonist-bound, computationally designed Hyd_high7 (green), Gt peptide-bound opsin active state (PDB: 3DQB, blue), and agonist- and mini-Gs-bound A2AR thermostabilized by alanine scanning (PDB: 5G53, magenta). The structures are ranked according to the magnitude of the conformational changes towards the fully active state (back arrow). c–e, Conformations of key conserved microswitches involved in receptor activation (in addition to the reference structures described in panel b, the following structure is also represented: agonist-bound A2AR thermostabilized by alanine scanning (PDB: 2YDO, grey)): the TMH3–TMH5–TMH6 interface around the W6.48 toggle switch in the receptor TM core (c), the TM3–TMH5–TMH6 interface around the intracellular DR3.50Y motif (d) and the TMH6–TMH7 interface around the intracellular N7.49PXXY7.53 motif (e). Conformational changes between inactive, partially active and fully active states are indicated by distinct black arrows. f–k, Comparison of the conformation of the designed (black) and neighbouring native (grey) residue microswitches in the agonist-bound designed Hyd_high7 structure (f) and the agonist- and Gs-bound Hyd_high7 design model (g) as well as the TMH5–TMH6 core interface (h), the TMH6–TMH7 intracellular interface with bound Gs (magenta) (i), the TMH1–TMH2–TMH3 core interface at the Na+ ion-conserved binding site (j) and the TMH3–TMH6–TMH7 intracellular interface (k).CC BY · PubMed Central
Fig. 6 Hyd-high7 is built of a synthetic allosteric interaction network distant from those of natural GPCRs. a,b, Solvent-mediated polar interaction networks at static–switchable TMH interfaces in mOR (PDB: 5C1M) (a) and Hyd_high7 (b) active-state X-ray structures. Water molecules are represented as red spheres and water-mediated hydrogen bonds are indicated by yellow dotted lines. c,d, Interhelical (direct and water-mediated) polar interactions in the TMH core region extracted from the mOR (5c1m) (c) and Hyd_high7 (d) active-state X-ray structures viewed perpendicular to the membrane bilayer. Polar interactions are represented as dotted red lines on the backbone representation of the receptors. The designed mutations rewired the inter-TMH polar interaction network by eliminating native polar contacts between TMH3 and TMHs 5–7 while enhancing TMH7 connections with neighbouring TMHs 1, 2 and 6.
Fig. 6 Hyd-high7 is built of a synthetic allosteric interaction network distant from those of natural GPCRs. a,b, Solvent-mediated polar interaction networks at static–switchable TMH interfaces in mOR (PDB: 5C1M) (a) and Hyd_high7 (b) active-state X-ray structures. Water molecules are represented as red spheres and water-mediated hydrogen bonds are indicated by yellow dotted lines. c,d, Interhelical (direct and water-mediated) polar interactions in the TMH core region extracted from the mOR (5c1m) (c) and Hyd_high7 (d) active-state X-ray structures viewed perpendicular to the membrane bilayer. Polar interactions are represented as dotted red lines on the backbone representation of the receptors. The designed mutations rewired the inter-TMH polar interaction network by eliminating native polar contacts between TMH3 and TMHs 5–7 while enhancing TMH7 connections with neighbouring TMHs 1, 2 and 6.CC BY · PubMed Central
Extended Data Fig. 1 Polar interaction networks connecting static and switchable TMHs in the Hyd_high2 design and A2AR. a, Cartoon representation of the Hyd_high2 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high2 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 1 Polar interaction networks connecting static and switchable TMHs in the Hyd_high2 design and A2AR. a, Cartoon representation of the Hyd_high2 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high2 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 2 Polar interaction networks connecting static and switchable TMHs in the Hyd_low2 design and A2AR. a, Cartoon representation of the Hyd_low2 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_low2 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 2 Polar interaction networks connecting static and switchable TMHs in the Hyd_low2 design and A2AR. a, Cartoon representation of the Hyd_low2 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_low2 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 3 Polar interaction networks connecting static and switchable TMHs in the Hyd_high3 design and A2AR. a, Cartoon representation of the Hyd_high3 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high3 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 3 Polar interaction networks connecting static and switchable TMHs in the Hyd_high3 design and A2AR. a, Cartoon representation of the Hyd_high3 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high3 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 4 Polar interaction networks connecting static and switchable TMHs in the Hyd_low3 design and A2AR. a, Cartoon representation of the Hyd_low3 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_low3 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 4 Polar interaction networks connecting static and switchable TMHs in the Hyd_low3 design and A2AR. a, Cartoon representation of the Hyd_low3 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_low3 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 5 Polar interaction networks connecting static and switchable TMHs in the Hyd_high4 design and A2AR. a, Cartoon representation of the Hyd_high4 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high4 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 5 Polar interaction networks connecting static and switchable TMHs in the Hyd_high4 design and A2AR. a, Cartoon representation of the Hyd_high4 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high4 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 6 Polar interaction networks connecting static and switchable TMHs in the Hyd_high5 design and A2AR. a, Cartoon representation of the Hyd_high5 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high5 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 6 Polar interaction networks connecting static and switchable TMHs in the Hyd_high5 design and A2AR. a, Cartoon representation of the Hyd_high5 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high5 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 7 Polar interaction networks connecting static and switchable TMHs in the Hyd_high1 design and A2AR. a, Cartoon representation of the Hyd_high1 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high1 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.
Extended Data Fig. 7 Polar interaction networks connecting static and switchable TMHs in the Hyd_high1 design and A2AR. a, Cartoon representation of the Hyd_high1 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high1 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation.CC BY · PubMed Central
Extended Data Fig. 8 Polar interaction networks connecting static and switchable TMHs in the Hyd_high6 design and A2AR. a, Cartoon representation of the Hyd_high6 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high6 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation. Y315 and E316 belong to the G-protein Gαs.
Extended Data Fig. 8 Polar interaction networks connecting static and switchable TMHs in the Hyd_high6 design and A2AR. a, Cartoon representation of the Hyd_high6 receptor active state model where the adenosine ligand and side-chains in the vicinity of the designed mutations are highlighted in sticks. b,c, zoomed view on the designed region in Hyd_high6 (b) and A2AR (c). Polar contacts are represented with red dotted lines. Residues are numbered using the Ballesteros Weinstein notation. Y315 and E316 belong to the G-protein Gαs.CC BY · PubMed Central
Extended Data Fig. 9 Fluorescent Thermal Stability Assay of T4Lcp- Hyd_high constructs. Measurements were performed for N-terminus circular permutant T4-Lyzozyme - designed Hyd_high 9 (green) and 7 (red) crystallization constructs bound to the ligand agonist CGS21680. Source data
Extended Data Fig. 9 Fluorescent Thermal Stability Assay of T4Lcp- Hyd_high constructs. Measurements were performed for N-terminus circular permutant T4-Lyzozyme - designed Hyd_high 9 (green) and 7 (red) crystallization constructs bound to the ligand agonist CGS21680. Source dataCC BY · PubMed Central
Extended Data Fig. 10 Structure prediction of the Hyd_high7 design by Alpha-Fold2. a, Superposition between Hyd_high7 X-ray structure and Alpha-Fold2 model (AF2). b, AF2 model colored according to pLDDT score. c, pAE plot. d, e, Structural differences between the models generated by SPaDES or AF2 and the Hyd_high7 X-ray structure (provided as mean ± std calculated over 10 lowest energy models). d, Cα RMSD over the transmembrane region. e, RMSD of designed residue side-chains after local alignment of surrounding 6 residues. Source data
Extended Data Fig. 10 Structure prediction of the Hyd_high7 design by Alpha-Fold2. a, Superposition between Hyd_high7 X-ray structure and Alpha-Fold2 model (AF2). b, AF2 model colored according to pLDDT score. c, pAE plot. d, e, Structural differences between the models generated by SPaDES or AF2 and the Hyd_high7 X-ray structure (provided as mean ± std calculated over 10 lowest energy models). d, Cα RMSD over the transmembrane region. e, RMSD of designed residue side-chains after local alignment of surrounding 6 residues. Source dataCC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 3 references

Research works printed but not linked: 50

  1. Changeux, J. P. & Christopoulos, A. Allosteric modulation as a unifying mechanism for receptor function and regulation. Cell 166, 1084–1102 (2016). · DOI
  2. Yuan, S., Hu, Z., Filipek, S. & Vogel, H. W2466.48 opens a gate for a continuous intrinsic water pathway during activation of the adenosine A2A receptor. Angew. Chem. Int. Ed. 54, 556–559 (2015). · DOI
  3. Yuan, S., Filipek, S., Palczewski, K. & Vogel, H. Activation of G-protein-coupled receptors correlates with the formation of a continuous internal water pathway. Nat. Commun. 5, 4733 (2014). · DOI
  4. Setny, P. & Wiśniewska, M. D. Water-mediated conformational preselection mechanism in substrate binding cooperativity to protein kinase A. Proc. Natl Acad. Sci. USA 115, 3852–3857 (2018). · DOI
  5. Ye, L. et al. Mechanistic insights into allosteric regulation of the A2A adenosine G protein-coupled receptor by physiological cations. Nat. Commun. 9, 1372 (2018). · DOI
  6. Lu, P. et al. Accurate computational design of multipass transmembrane proteins. Science 359, 1042–1046 (2018). · DOI
  7. Boyken, S. E. et al. De novo design of protein homo-oligomers with modular hydrogen-bond network-mediated specificity. Science 352, 680–687 (2016). · DOI
  8. Jacobs, T. M. et al. Design of structurally distinct proteins using strategies inspired by evolution. Science 352, 687–690 (2016). · DOI
  9. Dou, J. et al. De novo design of a fluorescence-activating β-barrel. Nature 561, 485–491 (2018). · DOI · Paper page
  10. Angel, T. E., Chance, M. R. & Palczewski, K. Conserved waters mediate structural and functional activation of family A (rhodopsin-like) G protein-coupled receptors. Proc. Natl Acad. Sci. USA 106, 8555–8560 (2009). · DOI · Paper page
  11. Venkatakrishnan, A. J. et al. Diverse GPCRs exhibit conserved water networks for stabilization and activation. Proc. Natl Acad. Sci. USA 116, 3288–3293 (2019). · DOI
  12. Zarzycka, B., Zaidi, S. A., Roth, B. L. & Katritch, V. Harnessing ion-binding sites for GPCR pharmacology. Pharmacol. Rev. 71, 571–595 (2019). · DOI
  13. Feng, X., Ambia, J., Chen, K.-Y. M., Young, M. & Barth, P. Computational design of ligand-binding membrane receptors with high selectivity. Nat. Chem. Biol. 13, 715–723 (2017). · DOI
  14. Lai, J. K., Ambia, J., Wang, Y. & Barth, P. Enhancing structure prediction and design of soluble and membrane proteins with explicit solvent–protein interactions. Structure 25, 1758–1770.e8 (2017). · DOI
  15. Hauser, A. S. et al. GPCR activation mechanisms across classes and macro/microscales. Nat. Struct. Mol. Biol. 28, 879–888 (2021). · DOI
  16. Chen, K. M., Keri, D. & Barth, P. Computational design of G protein-coupled receptor allosteric signal transductions. Nat. Chem. Biol. 16, 77–86 (2020). · DOI
  17. Rasmussen, S. G. et al. Crystal structure of the β2 adrenergic receptor–Gs protein complex. Nature 477, 549–555 (2011). · DOI · Paper page
  18. Vickery, O. N. et al. Intracellular transfer of Na+ in an active-state G-protein-coupled receptor. Structure 26, 171–180.e2 (2018). · DOI
  19. Magnani, F., Shibata, Y., Serrano-Vega, M. J. & Tate, C. G. Co-evolving stability and conformational homogeneity of the human adenosine A2a receptor. Proc. Natl Acad. Sci. USA 105, 10744–10749 (2008). · DOI
  20. Carpenter, B., Nehme, R., Warne, T., Leslie, A. G. & Tate, C. G. Structure of the adenosine A2A receptor bound to an engineered G protein. Nature 536, 104–107 (2016). · DOI
  21. Ahuja, S. & Smith, S. O. Multiple switches in G protein-coupled receptor activation. Trends Pharmacol. Sci. 30, 494–502 (2009). · DOI
  22. Nygaard, R., Frimurer, T. M., Holst, B., Rosenkilde, M. M. & Schwartz, T. W. Ligand binding and micro-switches in 7TM receptor structures. Trends Pharmacol. Sci. 30, 249–259 (2009). · DOI
  23. Lebon, G. et al. Agonist-bound adenosine A2A receptor structures reveal common features of GPCR activation. Nature 474, 521–525 (2011). · DOI
  24. Huang, W. et al. Structural insights into micro-opioid receptor activation. Nature 524, 315–321 (2015). · DOI
  25. Venkatakrishnan, A. J. et al. Molecular signatures of G-protein-coupled receptors. Nature 494, 185–194 (2013). · DOI
  26. White, K. L. et al. Structural connection between activation microswitch and allosteric sodium site in GPCR signaling. Structure 26, 259–269.e5 (2018). · DOI
  27. Massink, A. et al. Sodium ion binding pocket mutations and adenosine A2A receptor function. Mol. Pharmacol. 87, 305–313 (2015). · DOI
  28. Nygaard, R., Valentin-Hansen, L., Mokrosinski, J., Frimurer, T. M. & Schwartz, T. W. Conserved water-mediated hydrogen bond network between TM-I, -II, -VI, and -VII in 7TM receptor activation. J. Biol. Chem. 285, 19625–19636 (2010). · DOI
  29. Mravic, M. et al. Packing of apolar side chains enables accurate design of highly stable membrane proteins. Science 363, 1418–1423 (2019). · DOI
  30. Joh, N. H. et al. De novo design of a transmembrane Zn2+-transporting four-helix bundle. Science 346, 1520–1524 (2014). · DOI
  31. Yin, J. et al. Structure of a D2 dopamine receptor–G-protein complex in a lipid membrane. Nature 584, 125–129 (2020). · DOI
  32. Ingwersen, J. et al. Dual roles of the adenosine A2a receptor in autoimmune neuroinflammation. J. Neuroinflammation 13, 48 (2016). · DOI
  33. Graber, S. G., Figler, R. A. & Garrison, J. C. Expression and purification of functional G protein alpha subunits using a baculovirus expression system. J. Biol. Chem. 267, 1271–1278 (1992). · DOI
  34. Xu, F. et al. Structure of an agonist-bound human A2A adenosine receptor. Science 332, 322–327 (2011). · DOI
  35. Chen, K. Y., Sun, J., Salvo, J. S., Baker, D. & Barth, P. High-resolution modeling of transmembrane helical protein structures from distant homologues. PLoS Comput. Biol. 10, e1003636 (2014). · DOI
  36. Jaakola, V. P. et al. The 2.6 angstrom crystal structure of a human A2A adenosine receptor bound to an antagonist. Science 322, 1211–1217 (2008). · DOI
  37. Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003). · DOI · Paper page
  38. Pándy-Szekeres, G. et al. GPCRdb in 2018: adding GPCR structure models and ligands. Nucleic Acids Res. 46, D440–D446. · DOI
  39. Evans, R. et al. Protein complex prediction with AlphaFold-Multimer. Preprint at bioRxiv https://doi.org/10.1101/2021.10.04.463034 (2022). · DOI · Paper page
  40. Paradis, J. S. et al. Computationally designed GPCR quaternary structures bias signaling pathway activation. Nat. Commun. 13, 6826 (2022). · DOI
  41. Matsumura, M. & Matthews, B. W. Control of enzyme activity by an engineered disulfide bond. Science 243, 792–794 (1989). · DOI
  42. Cellitti, J. et al. Exploring subdomain cooperativity in T4 lysozyme I: structural and energetic studies of a circular permutant and protein fragment. Protein Sci. 16, 842–851 (2007). · DOI
  43. Caffrey, M. & Cherezov, V. Crystallizing membrane proteins using lipidic mesophases. Nat. Protoc. 4, 706–731 (2009). · DOI
  44. Vonrhein, C. et al. Data processing and analysis with the autoPROC toolbox. Acta Crystallogr. D 67, 293–302 (2011). · DOI
  45. Kabsch, W. Integration, scaling, space-group assignment and post-refinement. Acta Crystallogr. D 66, 133–144 (2010). · DOI
  46. Evans, P. Scaling and assessment of data quality. Acta Crystallogr. D 62, 72–82 (2006). · DOI
  47. Winn, M. D. et al. Overview of the CCP4 suite and current developments. Acta Crystallogr. D 67, 235–242 (2011). · DOI · Paper page
  48. Segala, E. et al. Controlling the dissociation of ligands from the adenosine A2A receptor through modulation of salt bridge strength. J. Med. Chem. 59, 6470–6479 (2016). · DOI
  49. Murshudov, G. N. et al. REFMAC5 for the refinement of macromolecular crystal structures. Acta Crystallogr. D 67, 355–367 (2011). · DOI
  50. DiMaio, F. et al. Improved low-resolution crystallographic refinement with Phenix and Rosetta. Nat. Methods 10, 1102–1104 (2013). · DOI

Metadata

Authors
Chen, K.-Y. M., Lai, J. K., Rudden, L. S. P., Wang, J., Russell, A. M., Conners, K., Rutter, M. E., Condon, B., Tung, F., Kodandapani, L., Chau, B., Zhao, X., Benach, J., Baker, K., Hembre, E. J., & Barth, P.
Affiliations
  • École Polytechnique Fédérale de Lausanne (CH)
  • Baylor College of Medicine (US)
  • AbbVie (United States) (US)
  • Eli Lilly (United States) (US)
Corresponding author
Barth, P.
ORCID
Venue
Nature Chemistry
Publisher
Springer Science and Business Media LLC
Published
2025-01-23
Volume
17
Issue
3
Pages
429-438
ISSN
  • 1755-4330 (ISSN-L)
  • 1755-4349
DOI
10.1038/s41557-024-01719-2
PMC
PMC11882447
OpenAlex
W4406768741
Cited by
14
References
3
Citations by year
2026: 9 · 2025: 5
FWCI
4.21
Type
journal article
Language
en
Open access
Licence
Locations
  • Nature Chemistry (publishedVersion, cc-by)
  • Research Square (acceptedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • Infoscience (Ecole Polytechnique Fédérale de Lausanne) (submittedVersion)
Funders
  • U.S. Department of Health & Human Services | NIH | National Institute of General Medical Sciences
  • Eli Lilly and Company
  • Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Receptor Mechanisms and Signaling
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Protein Structure and Dynamics
  • Physical Sciences › Chemistry › Organic Chemistry › Chemical Synthesis and Analysis
Keywords
cell surface receptors, G protein-coupled receptors, allosteric networks, protein-solvent interactions, protein design, thermal stability, protein-ligand interactions
MeSH terms
Allosteric Regulation, Humans, Ligands, Models, Molecular, Receptors, Cell Surface, Solvents, Protein Engineering, Signal Transduction

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review

How to cite

Cite CoChat for the checks on this page, and the paper for its own findings.

This page
CoChat. (2026). Citation check: Computational design of highly signalling-active membrane receptors through solvent-mediated allosteric networks. https://app.cochat.ai/papers/10.1038/s41557-024-01719-2
The paper
Chen, K.-Y. M., Lai, J. K., Rudden, L. S. P., Wang, J., Russell, A. M., Conners, K., Rutter, M. E., Condon, B., Tung, F., Kodandapani, L., Chau, B., Zhao, X., Benach, J., Baker, K., Hembre, E. J., & Barth, P. (2025). Computational design of highly signalling-active membrane receptors through solvent-mediated allosteric networks. Nature Chemistry. https://doi.org/10.1038/s41557-024-01719-2