==== Front Sci Rep Sci Rep Scientific Reports 2045-2322 Nature Publishing Group UK London 77659 10.1038/s41598-020-77659-x Article The G protein coupled receptor CXCR4 designed by the QTY code becomes more hydrophilic and retains cell signaling activity Tegler Lotta 12 Corin Karolina kcorin@gmail.com 136 Pick Horst 4 Brookes Jennifer 178 Skuhersky Michael 5 Vogel Horst 4 Zhang Shuguang shuguang@mit.edu 1 1 grid.116068.80000 0001 2341 2786Center for Bits and Atoms, Massachusetts Institute of Technology, Cambridge, MA 02139-4307 USA 2 grid.5640.70000 0001 2162 9922Molecular Biotechnology/IFM, Linköping University, 58183 Linköping, Sweden 3 grid.11951.3d0000 0004 1937 1135Biomedical Engineering Research Group, School of Electrical and Information Engineering, and Department of Molecular Medicine and Haematology, University of the Witwatersrand, Johannesburg, South Africa 4 grid.5333.60000000121839049Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland 5 grid.116068.80000 0001 2341 2786Synthetic Neurobiology Group, Media Lab, Massachusetts Institute of Technology, Cambridge, MA 02139-4307 USA 6 grid.19006.3e0000 0000 9632 6718Present Address: Department of Chemistry and Biochemistry, University of California, Los Angeles, CA 90095-1570 USA 7 grid.83440.3b0000000121901201Present Address: London Centre for Nanotechnology, University College London, 17-19 Gordon Street, London, WC1H 0AH UK 8 grid.83440.3b0000000121901201Present Address: Biophysics, Computational Physics, Quantum Physics, University College London, London, UK 7 12 2020 7 12 2020 2020 10 2137126 11 2019 13 11 2020 © The Author(s) 2020Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.G protein-coupled receptors (GPCRs) are vital for diverse biological functions, including vision, smell, and aging. They are involved in a wide range of diseases, and are among the most important targets of medicinal drugs. Tools that facilitate GPCR studies or GPCR-based technologies or therapies are thus critical to develop. Here we report using our QTY (glutamine, threonine, tyrosine) code to systematically replace 29 membrane-facing leucine, isoleucine, valine, and phenylalanine residues in the transmembrane α-helices of the GPCR CXCR4. This variant, CXCR4QTY29, became more hydrophilic, while retaining the ability to bind its ligand CXCL12. When transfected into HEK293 cells, it inserted into the cell membrane, and initiated cellular signaling. This QTY code has the potential to improve GPCR and membrane protein studies by making it possible to design functional hydrophilic receptors. This tool can be applied to diverse α-helical membrane proteins, and may aid in the development of other applications, including clinical therapies. Subject terms ProteinsG protein-coupled receptorsMembrane proteinshttp://dx.doi.org/10.13039/501100001337Claude Leon FoundationYang Trust FundWennerGren FoundationSwedish Chemical SocietyWellcome Trust UKOH2 LaboratoriesMIT Center for Bits and Atoms ConsortiumBay Valley Innovation Centerissue-copyright-statement© The Author(s) 2020 ==== Body Introduction G protein-coupled receptors (GPCRs) are arguably the most important class of proteins. They are the largest class of receptors, and are involved in numerous physiological functions including sight, smell, inflammation, taste, and the immune response1–5. Moreover, they are involved in high threat and prolific diseases like cancer, HIV infection, and diabetes, and are the targets of ~ 50% of pharmaceutical drugs3,6–11. The diverse and important functions of membrane proteins make them attractive templates for technological innovations. For example, photosystem 1 (PS1) can be integrated into chips to create photovoltaic cells12, and it has been hypothesized that olfactory receptors may be integrated into microfluidic devices to create smell sensors capable of detecting early-stage diseases and dangerous chemicals13,14. It may even be possible to design proteins de novo for novel applications. One problem currently hampering the development of membrane protein-based, and particularly GPCR-based, technologies is the need to stabilize them with detergents outside of the cell membrane. Although various surfactants have traditionally been used, they often do not adequately stabilize the proteins for sufficient periods of time. Moreover, the optimal detergent for a given receptor must be empirically determined, requiring a significant amount of time and resources. An alternate approach would be to design GPCRs to become hydrophilic, thereby eliminating the need for detergents. Not only would this aid in the development of GPCR and membrane protein-based technologies, but it would also further build on the basic principles for de novo protein design. With recent advancements in computational tools and software, researchers have started to design various proteins15–17. Several attempts with membrane proteins have focused on introducing hydrophilic amino acids like lysine or arginine to make them water-soluble. This goal has been accomplished with two small membrane proteins and one GPCR18–22. Although a computational approach was typically used18–20, one group substituted charged residues for hydrophobic ones21, and a second group “grafted” the hydrophilic face of a soluble protein onto the hydrophobic membrane-facing section of a membrane protein22. These strategies either require significant computational resources to design a water-soluble receptor, or may have limitations when applied to other membrane proteins. Moreover, detergents may still be needed for purification18. We previously reported an alternative approach for designing hydrophilic GPCRs, which we call the “QTY code”23. Using this code, hydrophobic amino acids are exchanged for hydrophilic residues with similar shapes and alpha-helix forming tendencies. Specifically, leucine (L) is replaced with glutamine (Q), valine (V) or isoleucine (I) are replaced with threonine (T), and phenylalanine (F) is replaced with tyrosine (Y) (Fig. 1). We showed that this code could be used to design hydrophilic chemokine receptors23. However, in order to make these receptors completely water-soluble, all hydrophobic transmembrane residues were replaced with non-charged hydrophilic ones. Such large-scale changes might not be tolerated as well by other GPCRs or helical membrane proteins. Moreover, the activity of these receptors in living cells was not determined.Figure 1 Schematic representations of CXCR4 and the proposed CXCR4QTY29 mutations. (A) Computational model of CXCR4QTY29. Hydrophobic residues are shown in yellow. The lipid-facing side faces the reader and the dimerization interface is on the other side. (B) Side view of CXCR4 (green) superimposed on the computational model of CXCR4QTY29 (blue). The helices nearly overlap, suggesting that the proposed mutations do not significantly alter the structure of CXCR4. (C) Extracellular view of the image in (B). The similarity in structures suggests that receptor binding and signaling may not be significantly affected. (D) Full-length sequence of the native CXCR4 aligned with CXCR4QTY29. The changed residues, which comprise only ~ 8% of the protein, are highlighted in red. A C-terminal GSSG linker and rho1D4 purification tag are underlined. The models in (A)–(C) were visualized within the simulation software YASARA38. Here we show that the QTY code can successfully be used to design a more hydrophilic form of CXCR4 that retains cellular signaling activity in HEK293 cells by replacing only 29 lipid-facing transmembrane L, I, V and F residues. CXCR4 was chosen because of its prominent role in 23 types of cancer metastasis and HIV infection24–26. Hence, a water-soluble variant could have potential therapeutic benefits, such as binding free CXCL12 to prevent cancer metastasis. Moreover, CXCR4’s ligand CXCL12 (SDF1α) and binding affinity are well characterized, and its crystal structure has been determined27, making it possible to make specific QTY replacements. CXCR4QTY29 was indeed more water-soluble than the native CXCR4 when expressed in a cell-free system. The secondary structures of the purified native and QTY variants were comparable, as were the ligand-binding affinities. Most importantly, CXCR4QTY29 was able to signal through the G protein Gαq in HEK293 cells, albeit at ~ 4× higher CXCL12 ligand concentrations. Our results demonstrate that the QTY code can be used to design biologically functional hydrophilic GPCRs. Furthermore, the QTY code can likely be applied to other helical membrane proteins, thus opening the possibility for a wide range of biotechnological or medical uses. Results Computational modeling of CXCR4QTY29 Since we currently have not yet obtained the crystal structure of CXCR4QTY29, we tested whether CXCR4QTY29 is stable by simulating it in an explicit water environment at 24.85 °C, pH7.4, and 0.9% NaCl for 1 µs. After an initial 0.3 µs using the AMBER14 force field software, the overall structure was already formed and seemed to be stable. Additional 0.7 µs simulations did not further stabilize the structure (Supplementary Information). After the simulations, CXCR4QTY29 was superimposed with the crystal structure of the natural receptor (Fig. 1). The comparison showed that CXCR4QTY29 and CXCR4 have very similar structures. This indicates that the mutations do not adversely affect the protein structure when compared to the native receptor. The structural conservation observed further suggests that the QTY substitutions would not significantly affect the fold of the α-helices, and perhaps would not drastically affect the interaction of the ligand in the binding domain. Expression and water-solubility of CXCR4QTY29 Escherichia coli-based cell free extracts were used to express native CXCR4 and CXCR4QTY29. The detergent Brij-35 was used during expression because it has been previously demonstrated to be the optimal detergent for cell-free expressed GPCRs28. Both receptors were expressed in the following concentrations of Brij-35 (w/v): 0%, 0.01%, and 0.2%. As expected, native CXCR4 was most soluble at a Brij-35 concentration of 0.2%, and least soluble without detergent (Fig. 2). While CXCR4QTY29 showed a similar pattern of water-solubility, it was more soluble than the native CXCR4 at 0% and 0.01% Brij-35 (p < 0.01). This suggests that, although CXCR4QTY29 is not completely water-soluble, it appears to be more hydrophilic and more water-soluble than native CXCR4.Figure 2 Relative insoluble fractions of CXCR4 and CXCR4QTY29 immediately after expression and 1 week later. (A) Insoluble fractions of CXCR4 and CXCR4QTY29 immediately after expression in 0%, 0.01%, and 0.2% Brij-35 (w/v). All data was normalized to CXCR4 in 0% Brij-35. Error bars show ± SEM. For 0% and 0.01% Brij-35, P < 0.01 (Student’s t-test, n = 12). (B) Insoluble fractions of CXCR4 and CXCR4QTY29 1 week after expression in 0%, 0.01%, and 0.2% Brij-35 (w/v). All data was normalized to CXCR4 in 0% Brij-35. Error bars show ± SEM. For 0% and 0.01% Brij-35, P < 0.01 (Student’s t-test, n = 12). We then assayed the ability of CXCR4 and CXCR4QTY29 to remain water-soluble over time since it is known that most membrane proteins are not very stable long-term even in detergents. For these experiments, the expressed receptors were stored at 4 ºC for 1 week, after which they were centrifuged to remove any insoluble protein fractions. This insoluble fraction was compared to the fraction that remained soluble. As expected, the highest insoluble fraction was seen when no detergent was used, the lowest was seen with 0.2% Brij-35, and an intermediate insoluble fraction was observed at 0.01% Brij-35 (Fig. 2). Interestingly, Brij-35 at a concentration of 0.2% was able to keep nearly 100% of both receptors soluble for 1 week (Fig. 2). At Brij-35 concentrations of 0.01% and 0%, CXCR4QTY29 was significantly more water-soluble than native CXCR4 (p < 0.01). This observation is similar to the results immediately after in vitro translation. This suggests that, while CXCR4QTY29 is not fully water-soluble, it is more soluble than the native receptor. Furthermore, it is able to remain more water-soluble for 1 week. Purification of CXCR4QTY29 Immunoaffinity chromatography was used to purify CXCR4 and CXCR4QTY29 expressed in 0.01% and 0.2% Brij-35. An antibody against the C-terminal rho-tag was used to detect both receptors on Western blots (Fig. 3A). CXCR4 ran at ~ 34 kDa, while CXCR4QTY29 ran at ~ 36 kDa. The predicted molecular weights are 40.9 kDa and 41 kDa, respectively. Although the observed molecular weights are smaller than the predicted molecular weights, the blots likely show full-length receptor expression: membrane proteins typically run smaller than their predicted molecular weights13,28, the relative weights of CXCR4 and CXCR4QTY29 are as expected, and detection was against the C-terminal rho-tag. Silver stains were used to assay the purity of the receptors (Fig. 3B). After immunoaffinity chromatography alone, the receptors were typically ~ 90% pure, which was sufficient for downstream analyses. Interestingly, while Western blotting and silver staining indicated that both receptors could dimerize, this tendency was reduced with CXCR4QTY29. This is not surprising since the CXCR4QTY29 variant most likely disrupted or changed the dimerization interface27.Figure 3 Western blot using a ρ1D4 monoclonal antibody, and silver stain of CXCR4 and CXCR4QTY29 in 0.01% and 0.2% Brij-35. (A) Western blot of CXCR4 and CXCR4QTY29. Lane 1: CXCR4 in 0.01% Brij-35; Lane 2: CXCR4QTY29 in 0.01% Brij-35; Lane 3: CXCR4 in 0.2% Brij-35; Lane 4: CXCR4QTY29 in 0.2% Brij-35. Although dimers are visible in both CXCR4 and CXCR4QTY29, a higher percentage of the expressed CXCR4 dimerizes. (B) Silver stain of purified CXCR4 and CXCR4QTY29. Lane 1: Ladder; Lane 2: CXCR4 in 0.01% Brij-35; Lane 3: CXCR4QTY29 in 0.01% Brij-35; Lane 4: CXCR4 in 0.2% Brij-35; Lane 5: CXCR4QTY29 in 0.2% Brij-35. Secondary structure analysis The secondary structures of CXCR4 and CXCR4QTY29 were analyzed using circular dichroism. The spectra of both receptors in the presence of 0.01%, and 0.2% FC14 were recorded. All samples exhibited characteristic alpha-helical spectra, with minima at ~ 208 nm and ~ 220 nm (Fig. 4). Because GPCRs have seven transmembrane alpha-helices, these results are expected thus suggesting proper folding. However, the spectra for CXCR4 and CXCR4QTY29 in 0.2% Brij-35 are different from the spectra in 0.01% Brij-35. Moreover, the spectra for receptors in 0.2% Brij-35 nearly overlap, while the spectra for receptors in 0.01% Brij-35 are similar but distinct. This suggests that there may be subtle differences in the secondary structures of CXCR4 and CXCR4QTY29 in 0.01% Brij-35, and that these structures are different from the secondary structures of both receptors in 0.2% Brij-35. However, it is also possible that any differences in the spectra are due to slight impurities. Taken together with the water-solubility data, the circular dichroism data suggests that the presence of detergent keeps the receptors properly folded, and that the structure of CXCR4QTY29 at 0.01% Brij-35 is slightly more similar to the structure in 0.2% Brij-35 than that of the native receptor.Figure 4 Circular dichroism spectra of CXCR4 and CXCR4QTY29 expressed in 0.01% or 0.2% Brij-35. All spectra have minima at ~ 208 nm and ~ 220 nm, which are characteristic for alpha-helical secondary structures. The spectra for receptors in 0.2% Brij-35 nearly overlap, suggesting nearly identical secondary structures. However, the spectra for the receptors in 0.01% Brij-35 are similar but distinct. Ligand-binding analysis We asked if CXCR4QTY29 interacts with the ligand CXCL12 (SDF1α) in the same fashion as the native CXCR4. Surface Plasmon Resonance (SPR) was used to carry out the ligand binding measurements. Both receptors were captured in different flow cells using the Rho-tag or 1D4 antibody immobilized on the surface, as previously described29. The same ligand preparations were first flowed over CXCR4, and were then flowed over CXCR4QTY29 for comparison. To obtain a reasonable fit to the SPR data, both CXCR4 and CXCR4QTY29 were analyzed with the “heterogenous ligand” method. The results suggest that the cell free protein expression system yields one protein fraction with lower activity, and one fraction with higher activity. Nevertheless, the more active fractions of both CXCR4 and CXCR4QTY29 have affinities that are similar to those previously reported for the natural CXCR4 isolated from cells. The KD values measured in this work are ~ 29 nM and ~ 56 nM, respectively, compared to reported CXCR4 affinities ranging from 5 to 200 nM29,30, depending on the preparation methods used (Fig. 5). Most importantly, the cell-free produced CXCR4QTY29 binds CXCL12 with a similar interaction profile as the cell-free produced wild type CXCR4, suggesting that the ligand CXCL12 interacts with both receptors in a similar manner.Figure 5 Surface Plasmon Resonance data of the interaction between SDF1α and (A) CXCR4 or (B), CXCR4QTY29. Colored lines represent the experimental data, while black lines represent fitted curves using the heterogeneous ligand-binding model in the BiaEvaluation program. Parameters from the fitting are shown in-line. The closeness of the fit, described by the statistical value χ2 (Chi2), was 1.78 for CXCR, and 13.4 for CXCR4QTY29, indicating that the model used is slightly more true to reality for CXCR4. The high-activity fractions of CXCR4QTY29 and wild-type CXCR4 have similar affinities and binding kinetics. The low-activity fraction of CXCR4QTY29 has a similar binding kinetic profile as CXCR4, but has a higher association rate constant yielding a lower affinity to the ligand. The active fractions of both receptors have comparable affinities to SDF1α. The KD values measured in this work are ~ 29 nM (CXCR4) and ~ 56 nM (CXCR4QTY29), which are comparable to reported CXCR4 affinities ranging from 5 to 200 nM29,30. CXCR4 and CXCR4QTY29 localization in HEK293 cells Since CXCR4QTY29 still has hydrophobic residues in its 7 transmembrane α-helices, we assayed whether it could still insert itself into the cell membrane. In order to visualize where the expressed receptors were localized, HEK293 cells from the same parent population were transfected in parallel with the genes that express either CXCR4 or CXCR4QTY29. Immunohistochemical staining was then performed to locate the expressed receptors. Both receptors trafficked to the cell membrane (Fig. 6). However, significantly more of the native CXCR4 was visualized in the membrane than CXCR4QTY29. Because CXCR4QTY29 is designed to be more water soluble, it is possible that less CXCR4QTY29 is able to remain tightly embedded in the lipid bilayer. It is also possible that it is less capable of effectively inserting into the cell membrane. This difference in cell-surface receptor concentration is likely to be responsible for the reduced cell-signaling response observed in CXCR4QTY29-transfected HEK293 cells (Fig. 7).Figure 6 Laser-scanning confocal micrographs showing fluorescently-labeled cell surface expression of wild type CXCR4 and CXCR4QTY29. Left hand: fluorescence images (Ex. 488 nm/Em. 505 nm, LP filter). Right hand: transmission images. (A) Negative control, showing stable HEK293 cells not induced for CXCR4 expression. (B) HEK293 cells expressing CXCR4 wildtype. (C) HEK293 cells expressing CXCR4QTY29. Equal amounts of CXCR4 and CXCR4QTY29 were transfected into HEK293 cells in parallel and under identical conditions. These images indicate that CXCR4 and CXCR4QTY29 traffic to the cell membrane, but that the amount of membrane-bound CXCR4QTY29 is less than that of CXCR4. The size bar is 10 mm. Figure 7 Intracellular calcium signaling in CXCR4- or CXCR4QTY29-transfected HEK293 cells in response to the ligand SDF-1α. CXCR4 or CXCR4QTY29 were co-transfected with the G-protein Gaq into HEK293 cells, and the cytosolic Ca++ concentration in response to SDF1α was monitored with FuraRed. (A)–(C) In CXCR4-transfected cells, the intracellular calcium response to application of 25 ng/ml (~ 3.2 nM) SDF1α and subsequent application of ATP (100 mM) as a viability control were recorded. (D)–(F) In CXCR4QTY29-transfected cells, the intracellular calcium response to application of 100 ng/ml (~ 13 nM) SDF1α and the subsequent application of ATP (100 mM) as a viability control were recorded. The fluorescence traces in (A) and (D) show the response of a single cell to both reagents. (B) and (E) show zoomed fields of transfected cells. (C) and (F) Show the frequency of individual Ca2+ responses in populations of CXCR4- and CXCR4QTY29-transfected HEK293 cells in response to SDF1α. Circles highlight active cells. In both cases, PBS or mock-transfected cells were used as negative controls, and the fluorescence in response to SDF1α was indistinguishable from baseline, while ATP could still elicit a normal response. Cell-signaling of CXCR4QTY29 in HEK293 cells The biologically functional activity of CXCR4QTY29 was assayed by co-transfecting it with the G-protein Gaq into HEK293 cells. For comparison, the native CXCR4 was co-transfected with Gaq into HEK293 cells in parallel and under identical experimental conditions. Both CXCR4 and CXCR4QTY29-transfected cells showed a transient increase in intracellular calcium upon application of the ligand CXCL12 (SDF1α) (Fig. 7). However, native CXCR4 seemed more readily able to initiate signal transduction than CXCR4QTY29: a strong response was seen from native CXCR4 at a CXCL12 (SDF1α) concentration of 25 ng/ml (~ 3.2 nM), while 100 ng/ml (~ 13 nM) CXCL12 typically needed to be applied to CXCR4QTY29 in order to observe a similar response. This concentration reflects a fourfold increase. Moreover, it appears that more cells expressing native CXCR4 signaled upon application of CXCL12 (SDF1α) than cells expressing the QTY variant. The decreased signaling response seen from CXCR4QTY29 likely reflects a lower concentration of biologically available receptors as demonstrated via membrane localization experiments (Fig. 6). However, it is also possible that the reduced signaling is also due to disruption of the dimer interface, or nuanced changes in the binding pocket. ATP was used as a positive control, while PBS and mock-transfected cells were used as negative controls (Fig. 7). After application of CXCL12 (SDF1α), the calcium levels were allowed to return to baseline levels, after which the cells were stimulated with ATP. In each case, a signaling response was observed. PBS was added to the cells either before or after CXCL12 (SDF1α) was applied. No change in calcium concentration was observed in both cases. Similarly, mock-transfected cells did not respond to either PBS or CXCL12 (SDF1α), but did respond to ATP. Taken together, this data clearly shows that CXCR4QTY29 is capable of signaling in response to the native ligand CXCL12 (SDF1α), albeit in a fourfold reduced signaling capability. This further suggests that the CXCR4QTY29 variant did not completely disturb the ligand-binding domain, and that CXCR4QTY29 is likely properly folded, integrated into the membrane though perhaps less efficiently, and is functional. Discussion In spite of the recent expansion in GPCR structures31, the number of known structures is still miniscule compared to the size and importance of this receptor family. Indeed, GPCRs comprise > 2% of the human genome and include 825 different types of receptors6,32. Moreover, they are the targets of ~ 50% of pharmaceutical drugs, which are used to treat diseases ranging from hypertension, congestive heart failure, coronary heart disease, asthma, pre-term labor, schizophrenia, Parkinson’s disease, HIV, cancer, and diabetes mellitus-induced renal damage33. Yet, the structures of only 88 receptors have been determined. This relatively small number of determined structures is solely due to the difficulty of expressing sufficient amounts of soluble, stable, and functional receptors for structure and function assays. It is necessary to determine the molecular structures of GPCRs in order to understand how they function at the molecular level, and to design specific drug targets, clinical therapies, and innovative protein-based technologies. Almost all of the crystallized GPCRs needed to be modified in order to enable sufficient expression and to facilitate crystallization. In spite of these modifications, detergents were still required, which have been shown to hamper crystallization34. Moreover, disordered regions still existed where the structure could not be determined, and the most successful strategy (insertion of T4L) can interfere with the ability of the modified GPCRs to signal. The QTY code we present here is a simple alternative method of obtaining GPCRs suitable for structure and function studies that offers advantages over previous approaches. We have shown that the QTY code can be used to design a more hydrophilic form of the receptor CXCR4 that is biologically active, and properly folded and functional after purification. Because only the helix portions facing the membrane are changed, the designed receptor is capable of signaling through G-proteins upon ligand-binding. And unlike most previous GPCR modifications and computer-modeling to make water-soluble variants, this method reduces the need for detergents. This opens the door for diverse applications including sensing devices, the development of decoy therapies, and beyond. We have previously shown that the QTY code can be used to generate completely water-soluble receptors. However, this required a significant number of modifications: those receptors had > 20% mutations, with ~ 50% residue replacements in the transmembrane regions23. Moreover, cell-signaling and trafficking assays were not performed, so it was unknown whether QTY-modified receptors could still insert into the membrane or remain biologically active. Here, we were able to design a more hydrophilic receptor that is capable of inserting into the cell membrane and signaling. This variant receptor had only 29 residues replaced (~ 8.2%), and these lipid-facing residues were targeted and rationally designed. Of significance, our data indicates that the mutant CXCR4QTY29 retains some cell surface expression in HEK293 cells. This suggests that CXCR4QTY29 folding, intracellular trafficking, and membrane insertion is preserved to some extent. Without proper protein folding, the receptor would accumulate in the endoplasmic reticulum, and receptor translocation to the surface plasma membrane would be arrested35. Furthermore, cell surface expression is a critical prerequisite for cell signaling. Our confocal images indicate that the ligand-binding site of CXCR4QTY29 is accessible on the surface of HEK293 cells as it binds an externally added, receptor-specific ligand. A qualitative comparison of HEK293 cells expressing wild-type CXCR4 or CXCR4QTY29 indicates visibly weaker expression of the QTY receptor. This implies that, while binding is preserved, the QTY modification has some influence over cell surface expression. This is in agreement with the SPR experiments reported in this study. The QTY code described here will likely facilitate more structure and function studies since more water-soluble receptors can be designed without the need for complex computational analyses or extensive library screens. This will not only aid in structural studies, but may also make expression easier. Moreover, it may aid in the design of other helical transmembrane proteins. The principal significance of the QTY code is that it is a simple and potentially general protein design approach that can be applied to other GPCRs without the need for extensive computational calculations. Moreover, this QTY code can likely be applied to other alpha-helical transmembrane proteins as well. And, because the resulting receptor is more water-soluble, it can not only be used for structure and function studies, but also has the potential to be used in novel technological developments requiring water-soluble receptors or novel protein designs. Methods CXCR4QTY29 gene and protein design As a GPCR, CXCR4 is characterized by 7 transmembrane alpha helices arranged in a barrel-like conformation. The inside of this barrel is more hydrophilic, while the outside contacting the lipid bilayer is hydrophobic. Loops connecting the helices are located either inside or outside the cell membrane. Only residues in the helical transmembrane segments facing the lipid bilayer were replaced with hydrophilic QTY residues (Fig. 1) Residues in the external loops that interact with the ligand were not changed. Seven amino acids have higher helical forming tendencies: leucine (L), glutamine (Q), threonine (T), phenylalanine (F), and tyrosine (Y). Q, T, and Y can form hydrogen bonds and are more hydrophilic than L, I, V, and F, but they bear no positive or negative charges. Moreover, the following pairs have similar shapes and sizes: L and Q, V or I and T, and F and Y. In our proposed scheme, Q replaced L, T replaced V or I, and Y replaced F (hence QTY). To design CXCR4QTY29, the crystal structure of CXCR4 (PDB 3ODU) was visualized using VMD software36 (http://www.ks.uiuc.edu/Research/vmd/), and the QTY method above was applied to residues facing the lipid bilayer. A total of 29 residues (8.2% of the total amino acid sequence) were changed: F36Y, I39T, F49Y, L50Q, I53T, V82T, L86Q, I89T, F93Y, V96T, V99T, I162T, L165Q, I169T, F172Y, I173T, V242T, I243T, L246Q, L253Q, I257T, F264Y, I286T, L290Q, F293Y, L297T, I300T, L301Q and F304Y. The protein sequence for the natural CXCR4 (P61073) was obtained from Uniprot (http://www.uniprot.org/uniprot/P61073), and the protein sequence for CXCR4QTY29 was designed as described above and ordered from Geneart (Thermo Fisher Scientific, Waltham, MA). Both sequences were optimized for expression in E. coli. A C-terminal rho tag was added for purification, and the constructs were cloned into the pIVex2.3 vector for cell-free expression, and the pcDNA4T/O vector for cell-signaling assays. The genes were amplified using Qiagen (Germantown, MD) plasmid mini and maxiprep kits. The full final sequence of CXCR4QTY29 is shown in Fig. 1. Computer simulations of CXCR4 and CXCR4QTY29 The structure for CXCR4 was obtained from the Protein Data Bank (code 3ODU, chain A). An initial structure for CXCR4QTY29 was obtained from the predicted sequence and the GOMoDo modeling server37. The simulations were done in an explicit water environment at 24.85 °C, pH7.4, and 0.9% NaCl for 1 µs using the full-atom AMBER14 N self-parameterizing force field within the simulation software YASARA38. The MD was simulated at a 5 fs timestep using a simulation cell 20 Å larger than each candidate, but with re-centering of the candidate at each timestep, and with generation and removal of water at the boundary to negate the possibility of non-physical boundary effects39. The CXCR4QTY29 model was then aligned with its detergent-encapsulated counterpart CXCR4 using MUSTANG48 and superimposed40. The computer used for the simulations was built with an Intel Core i7-6950X10-Core 3.0 GHz Processor, GIGABYTE GeForce GTX 1080 Video Card, and 16 GB of DDR4 2800 memory. Cell-free protein expression Cell-free expression has been previously described28. Briefly, native CXCR4 or CXCR4QTY29 DNA was added to a commercial E. coli cell-free reaction mix (ThermoFisher, Qiagen, or RiNA) according to the manufacturer’s instructions. The detergent Brij-35 was added at concentrations ranging from 0 to 0.2% (w/v). The completed reactions were centrifuged at 10,000 rpm for 10 min. The supernatant with the solubilized receptor was removed, and the pellet was resuspended in an equivalent volume of PBS. The relative quantities of solubilized and precipitated protein were assayed with a dot blot and analyzed using ImageJ41. Immunoaffinity purification CNBr-activated Sepharose 4B beads (GE Healthcare) coupled to the rho1D4 monoclonal antibody (Cell Essentials, Boston, MA) were used to purify CXCR4 and CXCR4QTY29 as previously described28. Briefly, solubilized protein was captured on the beads overnight at 4 °C. The beads were washed with wash buffer (PBS + 0.01% or 0.2% FC14) until all impurities had been removed. The receptor was eluted from the beads with elution buffer (PBS + 0.01% or 0.2%FC14 + 800 µM elution peptide). The eluted protein was concentrated using 50 kDa MWCO (0.2% detergent samples), 30 kDa MWCO (0.2% or 0.01% detergent samples), or 10 kDa MWCO (0.01% detergent samples) filter columns. To remove excess elution peptide, the samples were washed with ~ 67 × volumes of wash buffer in the centrifugation filter columns at 4 ºC, 1000–1200 g, for intervals of 1–2 min. The centrifugation filter columns were from Millipore (Burlington, MA), and the elution peptide Ac-TETSQVAPA-CONH2 was synthesized by CPC Scientific, Inc (Sunnyvale, CA). CXCR4 and CXCR4QTY29 detection and purity analysis Dot blots were used to analyze the relative solubility of different samples, while western blots and silver stains were used to detect the receptors and analyze their purity as previously described42,43. Samples were loaded into Novex 10% Bis–Tris Gels (Life Technologies, Waltham, MA) according to the manufacturer’s instructions, with the exception that the samples were incubated at room temperature for 10 min as boiling can cause receptor aggregation. After electrophoresis, the samples were transferred to a nitrocellulose membrane, blocked with milk (5% w/v non-fat dried milk in TBST) for 1 h, and incubated with the rho1D4 monoclonal antibody (1:3000 in TBST, 1 h at room temperature or overnight at 4 ºC). A goat anti-mouse HRP-conjugated secondary antibody (Pierce, Rockford, IL) and the ECL-Plus Kit (GE Healthcare, Pittsburgh, PA) were used to visualize the receptors. The SilverXpress kit (Life Technologies) was used according to the manufacturer’s instructions to do total protein stains. All images were captured on a Biorad GelDoc system. ImageJ was used to compare sample intensities and to analyze sample purity. Secondary structure analysis using circular dichroism Circular dichroism (CD) spectra were recorded at 15 ºC from 195 to 245 nm with a step size of 1 nm, and an averaging time of 4 s (Aviv Biomedical, model 410, or Jasco J-1500). The spectra were blanked to wash buffer, and represent an average of 3–5 repeat runs. A quartz sample cell with a path length of 1 mm was used. 300 µL of receptor sample were used for each experiment. Ligand-binding analysis SPR analyses were performed with Biacore T200 and Biacore 2000 instruments (GE Healthcare). Biacore CM4 series S chips were used, and the ID4 antibody was immobilized on the chips using standard EDC/NHS coupling chemistry. Immobilization was performed with 50 mg/ml 1D4 antibody in 10 mM Acetate buffer, pH 5.0 and contact times between 5 and 7 min to achieve immobilization levels of 8000–12,000 RU. The chemokine receptors CXCR4 and CXCR4QTY29 were captured to 2500–5000 RU from crude supernatants from cell free reactions, and the surface was washed with PBS supplemented with 0.2% FC-14 until the surface bleeding was less than 1 RU/min (40–180 min). All experiments were performed at 25 °C in 50 mM Hepes (pH 7.5), 5 mM MgCl2, 1 mM CaCl2, and 1 mg/ml BSA. The natural ligand SDf1α (7.8 kDa) was flown over the surfaces at concentrations of 0, 3.125, 6.25, 12.5, 25, 50, 100, 200 and 400 nM. The injections were made in triplicate with a blank injection before each injection due to slow dissociation. Analyses were made in BiaEvaluation software (GE Healthcare). Cell-signalling assays HEK293 cells derived from a single culture were seeded onto 0.18 mm thick glass slides at a density of 0.5 × 106 cells per mL and were grown overnight at 37 °C in the incubator to allow adhesion of the cells to the surface of the slides. In parallel and under identical conditions, the cells were transiently transfected with Gαq and either CXCR4 or CXCR4QTY29 DNA using Lipofectamine2000 (Invitrogen) according to the manufacturer’s instructions. Six hours after transfection, the cells were induced with doxycyclin. Eighteen to 20 h after induction, the cells were washed with PBS and loaded with 10 µM Fura-Red-AM (ThermoFisher Scientific, Waltham, MA) for 30 min in serum-free DMEM/F12 medium. This medium was then replaced with fresh DMEM supplemented with 10% serum, and the cells were incubated an additional 30 min before imaging. Calcium signaling in response to 25 ng (for CXCR4) or 100 ng (for CXCR4QTY29) SDF1a was visualized using confocal fluorescence microscopy (Zeiss LSM 510) with a water immersion objective (Zeiss Achroplan 63 × NA 1.2). The cells were excited at 488 nm (Ar + laser), and emission was monitored at 650 nm. Images were collected every 2 s for a total of 100 s. ATP was used as a positive control, and PBS and mock-transfected cells were used as negative controls. The cells were incubated at 37 ºC and 5% CO2 in DMEM supplemented with 10% FCS (Life Technologies) unless noted otherwise. Cell-staining assays HEK293 cells expressing CXCR4 or CXCR4QTY29 were grown on 0.18 mm thick cover glasses using a seeding density of 3 × 105 cells/ml in DMEM/F12 medium containing 10% FCS (GIBCO). Visualization of cell surface-expressed receptors was achieved by labeling the cells with 1 mM receptor-specific peptide 6 (Alexa 488) in PBS buffer for 10 min. Cells were washed with PBS buffer to remove unbound ligand. Images were recorded using a 488 nm Ar/Kr laser line and a LP505 mm filter on a Zeiss LSM510 laser-scanning confocal microscope with a 63× water objective (1.2 numerical aperture). Supplementary information Supplementary Information 1. Supplementary Information 2. Supplementary Information 3. Supplementary Figures. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. These authors contributed equally: Lotta Tegler and Karolina Corin. Supplementary information is available for this paper at 10.1038/s41598-020-77659-x. Acknowledgements This work was initially partially funded by the Yang Trust Fund. LTT gratefully acknowledges postdoctoral fellowships from the WennerGren Foundation and the Swedish Chemical Society. KC gratefully acknowledges postdoctoral fellowships from the Yang Trust Fund, Claude Leon Foundation, and the University of the Witwatersrand. JCB gratefully acknowledges the Wellcome Trust UK for the Sir Henry Wellcome Post-doctoral Research Fellowship. Author contributions The QTY code was conceived by S.Z.. Experiments were performed by L.T., K.C., H.P., and S.Z. Computational modeling was done by J.B. and M.S. H.V. provided instrumental contribution and scientific discussions. The manuscript was written by K.C., L.T., and H.P. All authors edited and revised the manuscript. Competing interests Massachusetts Institute of Technology (MIT) filed several patent applications for the QTY code, and OH2 Laboratories licensed the technology from MIT. This research was funded in part by OH2 Laboratories. S.Z. has a minor equity in OH2 Laboratories. S.Z., K.C., and L.T. are listed as authors on the patents. H.P., J.B., M.S., and H.V. declare no competing interests. ==== Refs References 1. Hill SJ Distribution, properties, and functional characteristics of three classes of Histamine receptor Pharmacol. Rev. 1990 42 1 45 83 2164693 2. Kedzierski RM Yanagisawa M Endothelin system: The double-edged sword in health and disease Annu. Rev. Pharmacol. Toxicol. 2001 41 851 876 10.1146/annurev.pharmtox.41.1.851 11264479 3. Ramachandran R Noorbakhsh F DeFea K Hollenberg MD Targeting proteinase-activated receptors: Therapeutic potential and challenges Nat. Rev. Drug Discov. 2012 11 1 69 86 10.1038/nrd3615 22212680 4. Kroeze WK Sheffler DJ Roth BL G-protein-coupled receptors at a glance J. Cell Sci. 2003 116 Pt 24 4867 4869 10.1242/jcs.00902 14625380 5. Latek D Modzelewska A Trzaskowski B Palczewski K Filipek S G protein-coupled receptors—Recent advances Acta Biochim. Pol. 2012 59 4 515 529 10.18388/abp.2012_2086 23251911 6. Dorsam RT Gutkind JS G-protein-coupled receptors and cancer Nat. Rev. Cancer. 2007 7 2 79 94 10.1038/nrc2069 17251915 7. Barton M Yanagisawa M Endothelin: 20 years from discovery to therapy Can. J. Physiol. Pharmacol. 2008 86 485 498 10.1139/Y08-059 18758495 8. Rosano L Spinella F Bagnato A Endothelin 1 in cancer: Biological implications and therapeutic opportunities Nat. Rev. Cancer. 2013 13 9 637 651 10.1038/nrc3546 23884378 9. Overington JP Al-lazikani B Hopkins AL How many drug targets are there? Nat. Rev. Drug Discov. 2006 5 12 993 996 10.1038/nrd2199 17139284 10. Sodhi A Montaner S Gutkind JS Viral hijacking of G-protein-coupled-receptor signalling networks Nat. Rev. Mol. Cell Biol. 2004 5 12 998 1012 10.1038/nrm1529 15573137 11. Oh DY Olefsky JM G protein-coupled receptors as targets for anti-diabetic therapeutics Nat. Rev. Drug Discov. 2016 15 3 161 173 10.1038/nrd.2015.4 26822831 12. Mershin A Self-assembled photosystem-I biophotovoltaics on nanostructured TiO(2) and ZnO Sci. Rep. 2012 10.1038/srep00234 22355747 13. Cook B Ernberg KE Chung H Zhang S Study of a synthetic human olfactory receptor 17–4: Expression and purification from an inducible mammalian cell line PLoS ONE 2008 3 8 e2920 10.1371/journal.pone.0002920 18682799 14. Corin, K. Olfactory-Related Receptors: Methods Towards Enabling Structural and Functional Studies (MIT PhD thesis, MIT Library, 2011). http://hdl.handle.net/1721.1/67203. 15. Fersht A Structure and Mechanism in Protein Science: A Guide to Enzyme Catalysis and Protein Folding 1998 New York W.H. Freeman 16. Ho SP DeGrado WF Design of a 4-helix bundle protein: Synthesis of peptides which self-associate into a helical protein J. Am. Chem. Soc. 1987 109 22 6751 6758 10.1021/ja00256a032 17. Regan L Protein design: Past, present, and future Biopolymers 2015 104 4 334 350 10.1002/bip.22639 25784145 18. Perez-Aguilar JM A computationally designed water-soluble variant of a G-protein-coupled receptor: The human mu opioid receptor PLoS ONE 2013 8 6 e66009 10.1371/journal.pone.0066009 23799068 19. Slovic AM Summa CM Lear JD DeGrado W Computational design of a water-soluble analog of phopholamban Protein Sci. 2003 12 2 337 348 10.1110/ps.0226603 12538897 20. Slovic AM Kono H Lear JD Saven JG DeGrado W Computational design of water-soluble analogues of the potassium channel KcsA Proc. Natl. Acad. Sci. 2004 101 7 1828 1833 10.1073/pnas.0306417101 14766985 21. Li H Cocco MJ Steitz TA Engelman DM Conversion of phospholamban into a soluble pentameric helical bundle Biochemistry 2001 40 22 6636 6645 10.1021/bi0026573 11380258 22. Frank S Toward a high-resolution structure of phopholamban: Design of soluble transmembrane domain mutants Biochemistry 2000 39 23 6825 6831 10.1021/bi0000972 10841762 23. Zhang S QTY code enables design of detergent-free chemokine receptors that retain ligand-binding activities Proc. Natl. Acad. Sci. USA 2018 115 37 E8652 E8659 10.1073/pnas.1811031115 30154163 24. Zhao H CXCR4 over-expression and survival in cancer: A system review and meta-analysis Oncotarget. 2015 6 7 5022 5040 10.18632/oncotarget.3217 25669980 25. Domanska UM A review on CXCR4/CXCL12 axis in oncology: No place to hide Eur. J. Cancer. 2013 49 1 219 230 10.1016/j.ejca.2012.05.005 22683307 26. Alkhatib G The biology of CCR5 and CXCR4 Curr. Opin. HIV AIDS. 2009 4 2 96 103 10.1097/COH.0b013e328324bbec 19339947 27. Wu B Structures of the CXCR4 chemokine GPCR with small-molecule and cyclic peptide antagonists Science 2010 330 6007 1066 1071 10.1126/science.1194396 20929726 28. Corin K A robust and rapid method of producing soluble, stable, and functional G-protein coupled receptors PLoS ONE 2011 6 10 e23036 10.1371/journal.pone.0023036 22039398 29. Navratilova I Dioszegi M Myszka DG Analyzing ligand and small molecule binding activity of solubilized GPCRs using biosensor technology Anal. Biochem. 2006 355 1 132 139 10.1016/j.ab.2006.04.021 16762304 30. Stenlund P Babcock GJ Sodroski J Myszka DG Capture and reconstitution of G protein-coupled receptors on a biosensor surface Anal. Biochem. 2003 316 2 243 250 10.1016/S0003-2697(03)00046-0 12711346 31. Membrane Proteins of Known Structure run by the Stephen White Laboratory at UCI (2018). http://blanco.biomol.uci.edu/mpstruc/ 32. Takeda S Kadowaki S Haga T Takaesu H Mitaku S Identification of G protein-coupled receptor genes from the human genome sequence FEBS Lett. 2002 520 1–3 97 101 10.1016/S0014-5793(02)02775-8 12044878 33. Insel PA Tang CM Hahntow I Michel MC Impact of GPCRs in clinical medicine: Genetic variants and drug targets Biochim. Biophys. Acta. 2007 1768 4 994 1005 10.1016/j.bbamem.2006.09.029 17081496 34. Prive CG Detergents for the stabilization and crystallization of membrane proteins Methods 2007 41 4 388 397 10.1016/j.ymeth.2007.01.007 17367711 35. Dong C Filipeanu CM Duvernay MT Wu G Regulation of G protein-coupled receptor export trafficking Biochim. Biophys. Acta 2007 1768 853 870 10.1016/j.bbamem.2006.09.008 17074298 36. Humphrey W Dalke A Schulten K VMD—visual molecular dynamics J. Mol. Gr. 1996 14 33 38 10.1016/0263-7855(96)00018-5 37. Sandal M GOMoDo: A GPCRs online modeling and docking webserver PLoS ONE 2013 8 e74092 10.1371/journal.pone.0074092 24058518 38. Salomon-Ferrer R Case DA Walker RC An overview of the Amber biomolecular simulation package WIREs Comput. Mol. Sci. 2013 3 198 210 10.1002/wcms.1121 39. Krieger E Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: Four approaches that performed well in CASP8 Proteins. 2009 77 Suppl. 9 114 122 10.1002/prot.22570 19768677 40. Konagurthu AS Whisstock JC Stuckey PJ Lesk AM MUSTANG: A multiple structural alignment algorithm Proteins. 2006 64 559 574 10.1002/prot.20921 16736488 41. Rasband, W.S. ImageJ, U.S. National Institute of Health, Bethesda, Maryland, USA, http://imagej.nih.gov.ij/, 1997–2016. 42. Corin K Designer lipid-like peptides: A class of detergents for studying functional olfactory receptors using commercial cell-free systems PLoS ONE 2011 6 11 e25067 10.1371/journal.pone.0025067 22132066 43. Corin K Structure and function analyses of the purified GPCR human vomeronasal type 1 receptor 1 Sci. Rep. 2011 10.1038/srep00172 22355687