Next Article in Journal
Acetone Ingestion Mimics a Fasting State to Improve Glucose Tolerance in a Mouse Model of Gestational Hyperglycemia
Next Article in Special Issue
Residue Folding Degree—Relationship to Secondary Structure Categories and Use as Collective Variable
Previous Article in Journal
Investigation of the Characteristics and Antibacterial Activity of Polymer-Modified Copper Oxide Nanoparticles
Previous Article in Special Issue
Residual Helicity at the Active Site of the Histidine Phosphocarrier, HPr, Modulates Binding Affinity to Its Natural Partners
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein

1
Department of Structural Molecular Science, School of Physical Sciences, SOKENDAI, The Graduate University for Advanced Studies, Shonan Village, Hayama 240-0193, Kanagawa, Japan
2
Research Center of Integrative Molecular Systems, Institute for Molecular Science (IMS), National Institutes of Natural Sciences (NINS), Okazaki 444-8585, Aichi, Japan
3
Department of Chemistry, Graduate School of Science, Chiba University, Chiba 263-8522, Japan
4
Protein Design Group, Exploratory Research Center on Life and Living Systems (ExCELLS), National Institutes of Natural Sciences (NINS), Okazaki 444-8585, Aichi, Japan
5
Membrane Protein Research Center, Chiba University, Chiba 263-8522, Japan
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(23), 12906; https://doi.org/10.3390/ijms222312906
Submission received: 24 September 2021 / Revised: 19 November 2021 / Accepted: 23 November 2021 / Published: 29 November 2021

Abstract

:
G-protein coupled receptors (GPCRs) are known for their low stability and large conformational changes upon transitions between multiple states. A widely used method for stabilizing these receptors is to make chimeric receptors by fusing soluble proteins (i.e., fusion partner proteins) into the intracellular loop 3 (ICL3) connecting the transmembrane helices 5 and 6 (TM5 and TM6). However, this fusion approach requires experimental trial and error to identify appropriate soluble proteins, residue positions, and linker lengths for making the fusion. Moreover, this approach has not provided state-targeting stabilization of GPCRs. Here, to rationally stabilize a class A GPCR, adenosine A2A receptor (A2AR) in a target state, we carried out the custom-made de novo design of α-helical fusion partner proteins, which can fix the conformation of TM5 and TM6 to that in an inactive state of A2AR through straight helical connections without any kinks or intervening loops. The chimeric A2AR fused with one of the designs (FiX1) exhibited increased thermal stability. Moreover, compared with the wild type, the binding affinity of the chimera against the agonist NECA was significantly decreased, whereas that against the inverse agonist ZM241385 was similar, indicating that the inactive state was selectively stabilized. Our strategy contributes to the rational state-targeting stabilization of GPCRs.

1. Introduction

G-protein coupled receptors (GPCRs) are the largest membrane protein family in the human genome [1]. GPCRs mediate signals related to diverse physiological functions; therefore, these receptors are major drug targets [2]. GPCRs are in equilibrium between multiple states [3] and are known for their innate instability, which causes difficulties in sample preparation, functional assays, or structure determinations. Therefore, methods for stabilization have been developed, including alanine scanning [4,5], directed evolution [6], and rational mutation [7,8].
The fusion partner strategy has been widely used, in which the intracellular loop 3 (ICL3) connecting the transmembrane helices 5 and 6 (TM5 and TM6) is replaced with soluble protein domains, such as T4-lysozyme [9,10], apocytochrome b562RIL (BRIL) [11,12], rubredoxin [13], and glycogen synthetase [14]. However, these fusion partner proteins and residue positions for fusion have been identified through experimental trial and error. Moreover, stabilizing GPCRs in a specific state using the fusion approach, which is useful for screening out state-dependent ligands and antibodies [15,16], has not been achieved.
Recently, principles for designing protein structures from scratch have been developed, which made it possible to create a wide range of new protein structures with high thermal stability [17,18,19]. We previously developed a method to create a diverse set of all-α protein structures ranging from bundle-like topologies with parallel-aligned helices to complicated ones with irregularly arranged helices [20]. Using the developed method, we sought to rationally design fusion partner proteins customized for not only thermally stabilizing GPCRs but also stabilizing them in a target state compared to the other states.
In this study, we designed fusion partner proteins customized for stabilizing one of the class A GPCRs, adenosine A2A receptor (A2AR), in an inactive state [21,22]. A2AR plays important physiological roles, such as the modulation of motor, vascular control, and immunosuppression; therefore, A2AR is a drug target for various diseases, including Parkinson’s disease, heart failure, and cancer [23]. Class A GPCRs are the largest subfamily of GPCRs, and the receptors in the class have been suggested to undergo large conformational changes in TM6 associated with TM5 upon the state transitions (Figure 1A). Therefore, we stabilized A2AR in the inactive state by making a fusion with de novo designed proteins customized to fix the conformation of the two helices in the inactive state (Figure 1B).

2. Results

2.1. Computational Design of α-Helical Fusion Partner Proteins

We assumed that the TM5 and TM6 conformation could be fixed in a specific state through straight helical connections between a fusion partner protein and A2AR. Therefore, we sought to design α-helical protein structures de novo, of which the N- and C-terminal helices are, respectively, connected to TM5 and TM6 of an inactive state A2AR structure (PDB ID: 3PWH; this structure is bound to the inverse agonist ZM241385 [22]) without any kinks or intervening loops (Figure 2) (details are described in the Materials and Methods).
Using 1688 globular all-α backbone structures with six helices [20] whose N- and C-terminal helices are close to each other, we elongated the N- and C-terminal helices by seven residues, respectively, to fuse with TM5 and TM6 of the A2AR inactive structure. We then selected a set of 389 backbone structures whose terminal helices were elongated without steric clash. From the generated set, we selected backbone structures whose N- and C- terminal helices are well-superimposable to TM5 and TM6 in the inactive-state A2AR structure by calculating root mean square deviation (RMSD) values for the main chain of superimposed residues.
Next, we designed amino-acid sequences that stabilize each of the selected backbone structures by carrying out the cycles of amino acid sequence optimization and optimization of the entire structure [26]. Among the resulting designs with tight core packing [27] and high compatibility between the local sequence and structure [17], the designs whose N- and C-terminal helices were better superimposable to TM5 and TM6 were selected. Note that the designed structures that were inside of the predicted membrane region or had clashes with A2AR after the fusion with A2AR were discarded. (The positional information of the membrane was obtained from the Orientation of Proteins in Membrane (OPM) database (https://opm.phar.umich.edu, accessed on 23 September 2021 [28]).
Next, we selected the designed proteins that exhibited funnel-shaped energy landscapes in Rosetta ab initio folding simulations [17]. Among the selected designs, we further selected those whose N- and C- terminal helices exhibited low fluctuation in molecular dynamics (MD) simulations for experimental characterization. Finally, the designed proteins, FiX1 and FiX2 (FiX stands for a Fusion partner protein customized for inactivation and eXtra stabilization), were selected (one of the residues in FiX2 was mutated manually using Foldit [29]. See the Material and Methods).
Figure 2. Computational protocol for designing fusion partner proteins. As an initial set of backbone structures, 1688 globular all-α backbone structures with six helices, extracted from the previously created all-α backbone structure library [20], were used. (Step 1) For the initial set, the N- and C- terminal helices of the structures were elongated to fuse into A2AR through TM5 and TM6. Then, 389 backbone structures whose helices were elongated without steric clash were selected. (Step 2) Backbone structures whose elongated N- and C- terminal helices were well-superimposable (the main-chain root mean square deviation (RMSD) for the fused region, equal to or less than 0.65 Å) to TM5 and TM6 of the inactive state A2AR structure were selected. (Step 3) For each selected backbone structure, amino acid sequences that stabilize the backbone structure were designed. (Step 4) Six designed structures were selected on the basis of the following criteria: the designs are superimposable to A2AR with main-chain RMSD values less than and equal to 0.4 Å for the fused region, have tight core packing [27] and high sequence-structure compatibility [17] and can be fused into A2AR without steric crash with A2AR and the membrane. (Step 5) Two designed proteins that exhibited funnel-shaped energy landscapes in Rosetta folding simulations [17] and whose N- and C- terminal helices show low fluctuation in molecular dynamics (MD) simulations were selected. In the MD simulations, blue and red lines show initial and averaged main-chain RMSD values, respectively. Details are described in the Results and Materials and Methods.
Figure 2. Computational protocol for designing fusion partner proteins. As an initial set of backbone structures, 1688 globular all-α backbone structures with six helices, extracted from the previously created all-α backbone structure library [20], were used. (Step 1) For the initial set, the N- and C- terminal helices of the structures were elongated to fuse into A2AR through TM5 and TM6. Then, 389 backbone structures whose helices were elongated without steric clash were selected. (Step 2) Backbone structures whose elongated N- and C- terminal helices were well-superimposable (the main-chain root mean square deviation (RMSD) for the fused region, equal to or less than 0.65 Å) to TM5 and TM6 of the inactive state A2AR structure were selected. (Step 3) For each selected backbone structure, amino acid sequences that stabilize the backbone structure were designed. (Step 4) Six designed structures were selected on the basis of the following criteria: the designs are superimposable to A2AR with main-chain RMSD values less than and equal to 0.4 Å for the fused region, have tight core packing [27] and high sequence-structure compatibility [17] and can be fused into A2AR without steric crash with A2AR and the membrane. (Step 5) Two designed proteins that exhibited funnel-shaped energy landscapes in Rosetta folding simulations [17] and whose N- and C- terminal helices show low fluctuation in molecular dynamics (MD) simulations were selected. In the MD simulations, blue and red lines show initial and averaged main-chain RMSD values, respectively. Details are described in the Results and Materials and Methods.
Ijms 22 12906 g002

2.2. Experimental Characterization of FiX1 and FiX2

We experimentally characterized the de novo designed proteins, FiX1 and FiX2, without A2AR Table S2. These two designs were expressed in Escherichia coli and purified using a Ni-NTA column. Both of the designs were found to be well expressed and highly soluble, and are then characterized via circular dichroism (CD) spectroscopy, size-exclusion chromatography combined with multi-angle light scattering (SEC-MALS), and 1H–15N heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR) spectroscopy (Figure 3). Both designs showed CD spectra of all-α proteins from 25 to 98 °C, were monomeric in SEC-MALS, and showed well-dispersed sharp NMR peaks. These results indicate that the designs fold into unique α-helical structures as monomers with high thermal stability.

2.3. Experimental Characterization of A2AR Fused with FiX1 and FiX2

The genes encoding the wild-type A2AR (A2AR WT) and chimeras with fusion partner proteins, A2AR fused with FiX1 (A2AR–FiX1), FiX2 (A2AR–FiX2), and BRIL (A2AR–BRIL), were constructed and expressed in yeast as described previously (for each construct, a red fluorescent protein (RFP) was appended at the C-terminus) Table S3 [30,31]. It is known that the innate instability of GPCRs gives rise to low yields upon detergent extraction from the membrane. As one of the stability metrics, we studied the solubilization efficiency upon detergent extraction (the ratio of the fluorescent intensity after solubilization using n-decyl-D-maltopyranoside (DM) over that before solubilization).
A2AR–FiX1 showed significantly improved solubilization efficiency compared to the wild-type; the efficiency was greater than that of A2AR–FiX2 (Table 1). Therefore, we further studied the stability of A2AR–FiX1 by measuring the apparent melting temperatures in the clear-native polyacrylamide gel electrophoresis (CN–PAGE) method [31]. The melting temperature was found to be significantly increased (Figure 4B), which is consistent with the solubilization efficiency results. For comparison, the solubilization efficiency and melting temperature for A2AR–BRIL were also measured; these values were comparable to those of A2AR–FiX1 (see Discussion).
Next, we investigated the ligand-binding affinities of A2AR WT, A2AR–BRIL, and A2AR–FiX1 using the radioligands of an inverse agonist [3H]-ZM241385 and an agonist [3H]-NECA (see Methods). The binding affinities of A2AR WT and A2AR–BRIL against ZM241385 were similar to each other, and those against NECA were also similar to each other (the obtained equilibrium dissociation constants (Kd) of A2AR WT and A2AR–BRIL against ZM241385 were 10.5 nM and 15.7 nM, respectively; Kd values of A2AR WT and A2AR–BRIL against NECA were 161.5 nM and 193.1 nM, respectively) (Table 2 and Figure 4C). However, the binding affinity of A2AR–FiX1 against ZM241385 was similar to those of A2AR WT and A2AR–BRIL, but significant binding against NECA was not observed (Table 2 and Figure 4C). Moreover, we investigated the binding affinity of A2AR–FiX1 to NECA using the inhibition assay, which was approximately 100-times lower than that of A2AR WT (Figure 4D).
For A2AR, the correspondence between the ligand-binding states and the conformational states has been well studied: the binding of the inverse agonist ZM241385 shifts the conformational equilibrium to the inactive state, and that of the agonist NECA shifts the equilibrium to the active state [32,33]. Therefore, these results indicate that the shifting of the conformational equilibrium of A2AR toward the inactive state by fusion with FiX1 was successful. The experimental structure information of A2AR–FiX1 is required to further support our conclusion.

3. Discussion

We succeeded in rationally designing a fusion partner protein, which thermally stabilized one of the class A GPCRs, A2AR, and stabilized it in an inactive state. We carried out the custom-made de novo design of fusion partner proteins, of which the N- and C- terminal helices are, respectably, connectable to the TM5 and TM6 in the inactive state without kinks or intervening loops. The de novo designed fusion partner proteins FiX1 and FiX2 were found to fold as monomers with high thermal stability. The fusion of A2AR with FiX1 was found to be not only thermally stabilized but also stabilized in the inactive state, as we designed.
We expected that A2AR–FiX1 would be more stable than A2AR–BRIL, since the melting temperature of FiX1 is over 98 °C in circular dichroism, which is far more than that of BRIL, around 65 °C [34]. However, the apparent melting temperatures of A2AR–FiX1 and A2AR–BRIL were almost the same. This suggests that the overall stability of the A2AR chimera is determined by the transmembrane helices or the other loops rather than the ICL3 region; therefore, the stabilization by making fusion proteins may be saturated.
Previous fusion partner strategies have used naturally occurring proteins or their mutants as fusion partner proteins. However, the number of naturally occurring proteins is limited, and their structures have been optimized during evolution to express their functions. Among naturally occurring protein structures, the ones in which the terminal helices are in close distances would be readily found. However, for fusion partner proteins to be connected to TM5 and TM6 in a specific state using straight helices, the terminal helices must have a specific distance, angle, and helical cycle. It would be difficult to find naturally occurring proteins whose terminal helix geometries exactly match all the three conditions. Moreover, most naturally occurring proteins are not stable; therefore, their folding ability can be impaired by only a few mutations. In contrast, our de novo protein design approach allows us to create stable proteins with specific helix geometries and high stability, without experimental trial and error.
The next question is whether our developed method can stabilize A2AR in other states (i.e., active-like intermediate or active states). GPCRs in the basal condition favor the inactive state over the other states [3]. Therefore, stabilization of A2AR in the other states may be more difficult than stabilization in the inactive state. Nevertheless, our success in the stabilization in the inactive state indicates that our developed method has the potential for stabilization of A2AR in the other states and further for the state-targeting stabilization of other GPCRs.

4. Materials and Methods

4.1. Selection of Backbone Structure Models for Fusion Partners

1688 backbone structures were extracted from the previously created all-α backbone structure library with six helices [20] with the following restrictions: (1) the maximum consecutive buried residues (a residue with accessible surface area < 5.0 Å2 calculated using FreeSASA (https://freesasa.github.io, accessed on 23 September 2021) [35] with a probe of radius 3.0 Å is regarded as buried) in a structure is less than 4, (2) the distance between the N- and C-terminal Cα atoms is less than 12.0 Å corresponding to a rough distance between TM5 and TM6, and (3) lower radius of gyration.
Next, the N- and C-terminal helices of the extracted backbones were extended by appending seven helical residues using the RosettaRemodel protocol [36] in Rosetta software (https://www.rosettacommons.org/software, accessed on 23 September 2021). The calculations were attempted 100 times, and if 100 backbone structures were successfully generated, their averaged structure was used as a backbone structure in the following calculation; ultimately, 389 averaged backbone structures were obtained. Next, among these structures, we selected those whose terminal helices were fusible to TM5 and TM6 of the inactive state A2AR structure (PDB: 3PWH, obtained from the PDB OPM database [28]).
Note that the structure of 3PWH is that of a thermally stabilized mutant; however, we used the wild-type sequence with the mutation N154Q for preventing glycosylation in the experiments. To this end, main-chain RMSD values were calculated by superimposing all pairs of three consecutive residues in the N- and C- terminal helices of the backbone structure (the residues are selected from those of the residue number from 2 to 11, and those from 113 to 123, respectively) against all pairs of the three consecutive residues in TM5 and TM6 at the cytoplasmic side (the residues were selected from those of the residue number from 204 to 211, and those from 219 to 229, respectively) (See Table S4). We selected the backbone structures that were superimposable with a main-chain RMSD value less than or equal to 0.65 Å as fusible ones; 64 backbone structures were obtained.

4.2. Sequence Design for Further Backbone Selection

We further screened 64 backbone structures via the sequence design of each backbone structure, followed by entire structure optimization, using the FlxbbDesign protocol in Rosetta (for the score function, talaris2014 [37] was used). In the sequence design, amino acid residue types used for each residue position were restricted based on the buriedness: hydrophobic residues were used in the protein core, hydrophilic residues on the surface, and both hydrophobic and hydrophilic residues at the boundary. Cysteine was not used to prevent unintentional disulfide bond formation; histidine was not used because of its several protonation states; glycine was used for the first and last helix residues. After the sequence design, the designs whose backbone ABEGO torsion patterns (“A” indicates the alpha region of the Ramachandran plot; “B” the beta region; “G” and “E”, the positive phi region; and “O”, the cis peptide conformation [38]) were different from those of original backbone structures were discarded. We performed the design calculations 50 times independently, and then selected backbone structures from which almost all designs were successfully generated without a change in the backbone ABEGO torsion pattern. Finally, we obtained three backbone structures.

4.3. Sequence Design

For each of the selected three backbone structures, sequence designs were performed 10,000 times, using the design protocol described in above section, with additional restrictions for used amino acid residue types. (1) When the backbone dihedral angle was classified as G based on the ABEGO classification [38], the amino acid type of the residue was fixed to glycine; (2) serine and threonine on α-helices were not used, except for the first and last helix residues, because these residues have a tendency to bend α-helices [39]; (3) positively charged residues, lysine and arginine, were not used in the first three helix residues, based on a previous report [40].

4.4. Selection Criteria after Sequence Design

After the sequence design, the designs were selected by the following criteria. (1) the main-chain RMSD value between the N- and C-terminal helices and TM5 and TM6 of A2AR is less than or equal to 0.4 Å and (2) designs with tight-core packing calculated by Rosetta Holes [27] (more than 0 and less than 2.0) and Packstat (more than 0.6) in Rosetta software (https://www.rosettacommons.org/software, accessed on 23 September 2021). Then, designed structures that had clashes with A2AR in a fused structure were discarded (a clash was identified by the distances between the Cα atoms of a design and A2AR being less than 5.5 Å). In addition, designed structures that were to be inside of the membrane were also discarded (it is not allowed that even one of the atoms of a designed protein in a fused structure is in the membrane region; the membrane region was obtained from the Orientation of Proteins in Membrane (OPM) database [28]). Moreover, designs with high compatibility between the local sequence and structure were selected in the following manner. For each nine-residue frame of a designed protein, 200 nine-residue fragments were collected from a non-redundant set of X-ray structures based on the sequence similarity and secondary structure prediction. Then, for each frame, the RMSD of the local structure against each of the 200 fragments was calculated. Designs were ranked according to the summation of the log-ratio of the fragments, for which the RMSD was less than 1.5 Å, across all nine-residue frames, and six design sequence with high values were selected.

4.5. Rosetta Folding Simulation

Energy landscapes of the designed sequences were obtained from Rosetta folding simulations [17]. For each amino acid sequence of designed proteins, 10,000 predicted structure models were generated starting from a completely extended structure. Furthermore, 200 energy-minimized structure models were generated starting from each of the designed protein structures. The energy landscape of each designed structure was evaluated by the shape of the scatter plot of the Rosetta score of the generated models versus the corresponding RMSD values to the designed structure. We confirmed that the predicted energy landscapes for all the designs were funnel-like.

4.6. Molecular Dynamics (MD) Simulation

MD simulations were performed to select designed structures whose N- and C-terminal helices did not fluctuate significantly in the simulations. Main-chain RMSD values were calculated by superimposing three consecutive residues in the N- and C- terminal helices of each snapshot structure generated during an MD trajectory against three consecutive residues in TM5 and TM6. The positions for the three consecutive residues were those used in the RMSD-based screening calculation (Table S4). The designs with average RMSD values of more than 0.75 Å or unexpected hydrophilic interactions in the MD simulation trajectories were discarded.
The AMBER16 software suite [41] was used to perform the MD simulations. Hydrogen atoms were added using the LEaP module in AMBER16, after removing those from the design models. A box with a 12 Å buffer of water models around the protein model in each direction was created. TIP3P [42] and AMBER ff99SB force fields [43] were used as the water model and protein force field, respectively. Periodic boundary conditions were set at a cut-off distance of 10 Å. Long-range electrostatic interactions were treated using the particle mesh Ewald method.
At the beginning of the simulations, energy minimization of the solvent was performed with harmonic restriction for protein atoms, and subsequently energy minimization without restriction was performed. Next, the temperature of the system was gradually increased from 0 to 300 K in 100 ps in an NVT ensemble with a Langevin thermostat and harmonic positional restriction for the protein atoms. After the heating step, a 100 ns MD simulation was performed at 1 atm at 300 K in an NPT ensemble with isotropic position scaling, setting one step as 0.002 ps.

4.7. A Manual Mutation Using Foldit

The Tyr residue at the position 67 in FiX2 was manually mutated to Leu using Foldit [29] to optimize the core packing.

4.8. Experiments of De Novo Designed Fusion Partner Proteins: Protein Expression and Purification

Plasmids with FiX1 or FiX2 DNA sequences between the NdeI and XhoI restriction sites in pET21b vectors were purchased from FASMAC (Kanagawa, Japan). E. coli BL21 Star (DE3) competent cells were transformed with the plasmids and cultured in MJ9 minimal media containing 15N-labeled ammonium sulfate as a nitrogen source, and 15N-labeled FiX1 and FiX2 were expressed. After the cells were spun down, they were suspended in phosphate-buffered saline (PBS) buffer, 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, and 1.8 mM KH2PO4 at pH 7.4 with BugBuster (EMD Millipore Corp., Billerica, MA, USA), protease inhibitor, lysozyme, and deoxyribonuclease.
From the cell lysates, FiX1 and FiX2 samples with a His-tag at C-terminus were purified using a Ni-NTA column. The purified samples were dialyzed against PBS buffer at pH 7.4. The purity of the FiX1 and FiX2 samples was confirmed via SDS-PAGE (Figure S1) and mass spectrometry.

4.9. Experiments of De Novo Designed Fusion Partner Proteins: Circular Dichroism (CD)

CD spectra were measured using J-1500 KS (JASCO Corp., Tokyo, Japan). By heating the samples from 20 to 98 °C at a rate of 1 °C per min, far-UV CD spectra were measured from 260 to 200 nm at various temperatures of 20, 40, 60, 80, and 98 °C using 10 μM FiX1 and FiX2 samples in PBS buffer (pH 7.4) in a 1-mm path length cuvette.

4.10. Experiments of De Novo Designed Fusion Partner Proteins: Size Exclusion Chromatography Combined with Multi-Angle Light Scattering (SEC-MALS)

SEC-MALS measurements were performed using a miniDAWN TREOS static light scattering detector (Wyatt Technology Corp., Santa Barbara, California, USA) and a high-performance liquid chromatography (HPLC) system (1260 Infinity LC, Agilent Technologies, Santa Clara, CA, USA). Approximately 180 μM FiX1 and FiX2 samples in PBS buffer (pH 7.4) were injected into a Superdex 75 increase 10/300 GL column (GE Healthcare) equilibrated with PBS at a flow rate of 0.5 mL/min. Sample concentrations were evaluated based on the absorbance at 280 nm detected by using HPLC system. Static light scattering data at the angles of 43.6°, 90.0°, and 136.4° were obtained using a 659 nm laser. The data were analyzed using ASTRA software (https://store.wyatt.com/shop/viscostar/viscostar-iii/astra-software/, accessed on 23 September 2021) (version 6.1.2, Wyatt Technology Corp., Santa Barbara, California, USA) with a dn/dc value of 0.185 mL/g.

4.11. Experiments of De Novo Designed Fusion Partner Proteins: 2D 1H–15N HSQC Measurement

For 15N-labeled FiX1 and FiX2 samples of 400 to 600 μΜ in 90% H2O/10% D2O PBS buffer (pH 7.4), 2D HSQC NMR spectrum measurements were performed using a JNM-ECA 600 MHz spectrometer (JEOL, Tokyo, Japan). The obtained NMR spectra were analyzed using the Delta NMR software (https://nmrsupport.jeol.com/Software, accessed on 23 September 2021) (version 5.2.1, JEOL, Tokyo, Japan).

4.12. Experiments of A2AR-Designed Fusion Partner Proteins: DNA Construction

The coding sequence of the human adenosine A2A receptor (A2AR) from residues 1−316 was amplified by using the polymerase chain reaction (PCR) method, in which N154 was replaced by Q to eliminate N-linked glycosylation [44]. The DNA fragment was inserted into the plasmid pDDGFP-2 [45], including TagRFP-His8 at the C-terminus [31]. The intracellular loop 3 (ICL3) of the A2AR (denote A2AR WT) was replaced with FiX1 or FiX2. The residue numbers of the de novo designed fusion partner proteins refer to the original amino acid sequences on the pET21b vectors (Table S1).

4.13. Experiments of A2AR-Designed Fusion Partner Proteins: Solubilization Efficiency

Wild-type A2AR and its variants (A2AR–BRIL, A2AR–FiX1, and A2AR–FiX2) were expressed in Saccharomyces cerevisiae strain FGY217, and the membranes were prepared as described previously [31]. Briefly, membranes were resuspended in a solubilization buffer (50 mM Tris, 120 mM NaCl, 20% glycerol, and 1 μg/mL 6-amidinonaphthalen-2-yl 4-guanidinobenzoate bis (methanesulfonate) (AFP) (Alfresa Pharma Corp., Osaka, Japan); pH 8.0). The membrane suspension (5 mg/mL) was solubilized using n-decyl-β-D-maltopyranoside (DM) (final concentration, 1%) (Anatrace, Maumee, OH, USA) for 30 min at 4 °C. The red fluorescent protein (RFP) intensity was measured before and after solubilization at 595 nm (excitation at 535 nm) using a FilterMax F5 microplate reader (Molecular Devices, Sunnyvale, CA, USA). The solubilization ratio was evaluated as the RFP intensity of the unpurified A2AR–RFP divided by that of the whole membrane protein mixture soon after DM solubilization.

4.14. Experiments of A2AR-Designed Fusion Partner Proteins: Clear-Native PAGE

We evaluated the apparent melting temperatures of the wild-type A2AR and its variants (A2AR–BRIL, A2AR–FiX1) solubilized in 1% n-dodecyl β-D-maltopyranoside (DDM; Anatrace, Maumee, OH, USA) containing 0.2% cholesterol hemisuccinate (CHS; Sigma-Aldrich, Saint Louis, MO, USA), using clear native polyacrylamide gel electrophoresis (CN-PAGE) with modified Coomassie Brilliant Blue G-250 (mCBB) stain [31]. The samples were fused with RFP at the C-terminus, exhibiting fluorescence at 595 nm.
The samples were heated at each prescribed temperature for 5 min (the temperatures were prescribed in the range of 25–80 °C), and then immediately cooled on ice. CN-PAGE was performed using 10% Tris-glycine separation gel applied to the treated samples with CN-PAGE buffer (200 mM Tris-HCl, 20% glycerol, 1.0% mCBB, and 1.0% DDM; pH 8.6) at a ratio of 1:1. The samples on the gel were visualized (i.e., gel imaging was performed) using FUSION SOLO 7S (Vilber–Lourmat, Marne-la-Vallée, France) after a 5-s exposure to green light at 530 nm with a 655 nm cutoff filter.
The melting temperatures of the samples were determined from the fluorescence intensities of the monomeric bands on the CN-PAGE gel. The normalized fluorescence intensity was calculated by dividing the fluorescence intensity of the monomeric bands after heating by that before heating and is represented as a percentage. The obtained intensities of the monomeric bands were quantified using the ImageJ software (https://imagej.nih.gov/ij, accessed on 23 September 2021). The melting temperatures were calculated using GraphPad Prism 4.0 (GraphPad Software, San Diego, CA, USA) as previously described [31].

4.15. Experiments of A2AR-Designed Fusion Partner Proteins: Radioligand Binding Assay

Radioligand binding assays were performed using yeast cell membranes expressing the wild-type A2AR and its variants (A2AR–BRIL and A2AR–FiX1). The protein concentrations of the membranes were determined by using the bicinchoninic acid method (Thermo Fisher Scientific, Waltham, MA, USA) with bovine serum albumin as a standard. All experiments were performed in triplicate (independent expressions). For the saturation-binding assay, 10 µg of membranes were incubated (3 h on ice) with the inverse agonist [3H]-ZM241385 (American Radiolabeled Chemicals, Saint Louis, MO, USA) at concentrations ranging from 5 to 80 nM and the agonist [3H]-NECA (PerkinElmer, Waltham, MA, USA) at concentrations ranging from 25 to 600 nM. Non-specific binding was determined in the presence of 10 μM ZM241385 (Sigma-Aldrich, St. Louis, MO, USA) and 100 μM NECA (Sigma-Aldrich, St. Louis, MO, USA), respectively.
For the competition-binding assay, 10 µg of membranes were incubated with 20 nM of [3H]-ZM241385 and unlabeled NECA at concentrations ranging from 1 nM to 100 µM of [3H]-NECA for 3 h on ice. The unbound ligand was removed by rapid vacuum filtration over GF/F filters (GE Healthcare, Chicago, IL, USA). Filtration was performed using a MINI-VAC (Yamato Scientific, Tokyo, Japan) at room temperature. The filters were washed twice with solubilization buffer. After adding 5 mL of Filter-Count (PerkinElmer, Waltham, MA, USA), radioactivity was measured using an LSC-6100 liquid scintillation counter (Hitachi ALOKA Medical, Tokyo, Japan). The collected data were analyzed by using a nonlinear regression-fitting program in GraphPad Prism 8 (GraphPad Software, San Diego, CA, USA).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms222312906/s1. Figure S1, SDS-PAGE results for FiX1 and FiX2. Table S1, Sequence regions for A2AR and fusion partner proteins in A2AR chimeras. Table S2, Amino acid sequences of FiX1 and FiX2. Table S3, Amino acid sequences of A2AR fused with FiX1 and FiX2. Table S4, Residues used for superposition between A2AR and de novo designed fusion partner proteins.

Author Contributions

M.M., T.K., T.M. and N.K. designed the study. M.M. computationally designed fusion partner proteins and performed all experiments for the designed fusion partner proteins. K.S. and R.N. performed the radioligand binding assay and competition binding assays. K.K. and R.N. performed the clear-native PAGE experiments. M.M., K.S., K.K., T.K., T.M. and N.K. wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI Grants-in-Aid for Scientific Research 18H05420 to T.K. and N.K., and 18H05425 to T.M.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank the Functional Genomics Facility, NIBB Core Research Facilities, especially Yumiko Makino, for mass spectrometry analysis, and the Instrument Center, Okazaki, Japan, for HSQC spectra measurements. We thank Koya Sakuma for extracting the 1688 all-α backbone structure models, Nanao Suzuki for performing gene construction for A2AR-FiX1 and A2AR-FiX2 and solubilization efficiency experiments, and Shintaro Minami for providing the useful computational tools for protein design. We also thank Rie Koga for the constructive comments and discussions on the manuscript. The computations were performed using the Research Center for Computational Science (RCCS), Okazaki, Japan.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lagerström, M.C.; Schiöth, H.B. Structural diversity of G protein-coupled receptors and significance for drug discovery. Nat. Rev. Drug Discov. 2008, 7, 339–357. [Google Scholar] [CrossRef]
  2. Hauser, A.S.; Attwood, M.M.; Rask-Andersen, M.; Schiöth, H.B.; Gloriam, D.E. Trends in GPCR drug discovery: New agents, targets and indications. Nat. Rev. Drug Discov. 2017, 16, 829–842. [Google Scholar] [CrossRef]
  3. Weis, W.I.; Kobilka, B.K. The molecular basis of G protein–coupled receptor activation. Annu. Rev. Biochem. 2018, 87, 897–919. [Google Scholar] [CrossRef]
  4. 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 2008, 105, 10744–10749. [Google Scholar] [CrossRef] [Green Version]
  5. Serrano-Vega, M.J.; Magnani, F.; Shibata, Y.; Tate, C.G. Conformational thermostabilization of the β1-adrenergic receptor in a detergent-resistant form. Proc. Natl. Acad. Sci. USA 2008, 105, 877–882. [Google Scholar] [CrossRef] [Green Version]
  6. Sarkar, C.A.; Dodevski, I.; Kenig, M.; Dudli, S.; Mohr, A.; Hermans, E.; Plückthun, A. Directed evolution of a G protein-coupled receptor for expression, stability, and binding selectivity. Proc. Natl. Acad. Sci. USA 2008, 105, 14808–14813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Chen, K.Y.; Zhou, F.; Fryszczyn, B.G.; Barth, P. Naturally evolved G protein-coupled receptors adopt metastable conformations. Proc. Natl. Acad. Sci. USA 2012, 109, 13284–13289. [Google Scholar] [CrossRef] [Green Version]
  8. Popov, P.; Peng, Y.; Shen, L.; Stevens, R.C.; Cherezov, V.; Liu, Z.-J.; Katritch, V. Computational design of thermostabilizing point mutations for G protein-coupled receptors. eLife 2018, 7, e34729. [Google Scholar] [CrossRef] [PubMed]
  9. Cherezov, V.; Rosenbaum, D.M.; Hanson, M.A.; Rasmussen, S.G.F.; Thian, F.S.; Kobilka, T.S.; Choi, H.-J.; Kuhn, P.; Weis, W.I.; Kobilka, B.K.; et al. High-resolution crystal structure of an engineered human beta2-adrenergic G protein-coupled receptor. Science 2007, 318, 1258–1265. [Google Scholar] [CrossRef] [Green Version]
  10. Rosenbaum, D.M.; Cherezov, V.; Hanson, M.A.; Rasmussen, S.G.; Thian, F.S.; Kobilka, T.S.; Choi, H.J.; Yao, X.J.; Weis, W.I.; Stevens, R.C.; et al. GPCR engineering yields high-resolution structural insights into beta2-adrenergic receptor function. Science 2007, 318, 1266–1273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Chun, E.; Thompson, A.A.; Liu, W.; Roth, C.B.; Griffith, M.T.; Katritch, V.; Kunken, J.; Xu, F.; Cherezov, V.; Hanson, M.A.; et al. Fusion partner toolchest for the stabilization and crystallization of G protein-coupled receptors. Structure 2012, 20, 967–976. [Google Scholar] [CrossRef] [Green Version]
  12. Liu, W.; Chun, E.; Thompson, A.A.; Chubukov, P.; Xu, F.; Katritch, V.; Han, G.W.; Roth, C.B.; Heitman, L.H.; Ijzerman, A.P.; et al. Structural basis for allosteric regulation of GPCRs by sodium ions. Science 2012, 337, 232–236. [Google Scholar] [CrossRef] [Green Version]
  13. Tan, Q.; Zhu, Y.; Li, J.; Chen, Z.; Han, G.W.; Kufareva, I.; Li, T.; Ma, L.; Fenalti, G.; Li, J.; et al. Structure of the CCR5 chemokine receptor-HIV entry inhibitor maraviroc complex. Science 2013, 341, 1387–1390. [Google Scholar] [CrossRef] [Green Version]
  14. Yin, J.; Mobarec, J.C.; Kolb, P.; Rosenbaum, D.M. Crystal structure of the human OX2 orexin receptor bound to the insomnia drug suvorexant. Nature 2015, 519, 247–250. [Google Scholar] [CrossRef]
  15. Jo, M.; Jung, S.T. Engineering therapeutic antibodies targeting G-protein–coupled receptors. Exp. Mol. Med. 2016, 48, e207. [Google Scholar] [CrossRef]
  16. Yang, D.; Zhou, Q.; Labroska, V.; Qin, S.; Darbalaei, S.; Wu, Y.; Yuliantie, E.; Xie, L.; Tao, H.; Cheng, J.; et al. G protein-coupled receptors: Structure- and function-based drug discovery. Signal Transduct. Target. Ther. 2021, 6, 7. [Google Scholar] [CrossRef] [PubMed]
  17. Koga, N.; Tatsumi-Koga, R.; Liu, G.; Xiao, R.; Acton, T.B.; Montelione, G.T.; Baker, D. Principles for designing ideal protein structures. Nature 2012, 491, 222–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Huang, P.-S.; Boyken, S.E.; Baker, D. The coming of age of de novo protein design. Nature 2016, 537, 320–327. [Google Scholar] [CrossRef]
  19. Koga, R.; Koga, N. Consistency principle for protein design. Biophys. Phys. 2019, 16, 304–309. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Sakuma, K.; Kobayashi, N.; Sugiki, T.; Nagashima, T.; Fujiwara, T.; Suzuki, K.; Kobayashi, N.; Murata, T.; Kosugi, T.; Koga, R.; et al. Design of complicated all-α protein structures. bioRxiv 2021. [Google Scholar] [CrossRef]
  21. Jaakola, V.P.; Griffith, M.T.; Hanson, M.A.; Cherezov, V.; Chien, E.Y.; Lane, J.R.; Ijzerman, A.P.; Stevens, R.C. The 2.6 angstrom crystal structure of a human A2A adenosine receptor bound to an antagonist. Science 2008, 322, 1211–1217. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Doré, A.S.; Robertson, N.; Errey, J.C.; Ng, I.; Hollenstein, K.; Tehan, B.; Hurrell, E.; Bennett, K.; Congreve, M.; Magnani, F.; et al. Structure of the adenosine A2A receptor in complex with ZM241385 and the xanthines XAC and caffeine. Structure 2011, 19, 1283–1293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Borea, P.A.; Gessi, S.; Merighi, S.; Vincenzi, F.; Varani, K. Pharmacology of adenosine receptors: The state of the art. Physiol. Rev. 2018, 98, 1591–1625. [Google Scholar] [CrossRef]
  24. Lebon, G.; Warne, T.; Edwards, P.C.; Bennett, K.; Langmead, C.J.; Leslie, A.G.W.; Tate, C.G. Agonist-bound adenosine A2A receptor structures reveal common features of GPCR activation. Nature 2011, 474, 521–525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Carpenter, B.; Nehmé, R.; Warne, T.; Leslie, A.G.W.; Tate, C.G. Structure of the adenosine A2A receptor bound to an engineered G protein. Nature 2016, 536, 104–107. [Google Scholar] [CrossRef] [PubMed]
  26. Kuhlman, B.; Dantas, G.; Ireton, G.C.; Varani, G.; Stoddard, B.L.; Baker, D. Design of a novel globular protein fold with atomic-level accuracy. Science 2003, 302, 1364–1368. [Google Scholar] [CrossRef] [Green Version]
  27. Sheffler, W.; Baker, D. RosettaHoles2: A volumetric packing measure for protein structure refinement and validation. Protein Sci. 2010, 19, 1991–1995. [Google Scholar] [CrossRef] [Green Version]
  28. Lomize, M.A.; Pogozheva, I.D.; Joo, H.; Mosberg, H.I.; Lomize, A.L. OPM database and PPM web server: Resources for positioning of proteins in membranes. Nucleic Acids Res. 2012, 40, D370–D376. [Google Scholar] [CrossRef] [PubMed]
  29. Kleffner, R.; Flatten, J.; Leaver-Fay, A.; Baker, D.; Siegel, J.B.; Khatib, F.; Cooper, S. Foldit Standalone: A video game-derived protein structure manipulation interface using Rosetta. Bioinformatics 2017, 33, 2765–2767. [Google Scholar] [CrossRef]
  30. Shiroishi, M.; Tsujimoto, H.; Makyio, H.; Asada, H.; Yurugi-Kobayashi, T.; Shimamura, T.; Murata, T.; Nomura, N.; Haga, T.; Iwata, S.; et al. Platform for the rapid construction and evaluation of GPCRs for crystallography in Saccharomyces cerevisiae. Microb. Cell Factories 2012, 11, 78. [Google Scholar] [CrossRef] [Green Version]
  31. Suzuki, N.; Takamuku, Y.; Asakawa, T.; Inai, M.; Hino, T.; Iwata, S.; Kan, T.; Murata, T. An efficient screening method for purifying and crystallizing membrane proteins using modified clear-native PAGE. Anal. Biochem. 2018, 548, 7–14. [Google Scholar] [CrossRef]
  32. Ye, L.; Van Eps, N.; Zimmer, M.; Ernst, O.P.; Scott Prosser, R. Activation of the A2A adenosine G-protein-coupled receptor by conformational selection. Nature 2016, 533, 265–268. [Google Scholar] [CrossRef]
  33. Huang, S.K.; Pandey, A.; Tran, D.P.; Villanueva, N.L.; Kitao, A.; Sunahara, R.K.; Sljoka, A.; Prosser, R.S. Delineating the conformational landscape of the adenosine A2A receptor during G protein coupling. Cell 2021, 184, 1884–1894. e14. [Google Scholar] [CrossRef] [PubMed]
  34. Chu, R.; Takei, J.; Knowlton, J.R.; Andrykovitch, M.; Pei, W.; Kajava, A.V.; Steinbach, P.J.; Ji, X.; Bai, Y. Redesign of a four-helix bundle protein by phage display coupled with proteolysis and structural characterization by NMR and X-ray crystallography. J. Mol. Biol. 2002, 323, 253–262. [Google Scholar] [CrossRef]
  35. Mitternacht, S. FreeSASA: An open source C library for solvent accessible surface area calculations. F1000Research 2016, 5, 189. [Google Scholar] [CrossRef] [PubMed]
  36. Huang, P.-S.; Ban, Y.-E.A.; Richter, F.; Andre, I.; Vernon, R.; Schief, W.R.; Baker, D. RosettaRemodel: A Generalized Framework for Flexible Backbone Protein Design. PLoS ONE 2011, 6, e24109. [Google Scholar] [CrossRef] [Green Version]
  37. O’Meara, M.J.; Leaver-Fay, A.; Tyka, M.D.; Stein, A.; Houlihan, K.; DiMaio, F.; Bradley, P.; Kortemme, T.; Baker, D.; Snoeyink, J.; et al. Combined covalent-electrostatic model of hydrogen bonding improves structure prediction with rosetta. J. Chem. Theory Comput. 2015, 11, 609–622. [Google Scholar] [CrossRef] [Green Version]
  38. Lin, Y.R.; Koga, N.; Tatsumi-Koga, R.; Liu, G.; Clouser, A.F.; Montelione, G.T.; Baker, D. Control over overall shape and size in de novo designed proteins. Proc. Natl. Acad. Sci. USA 2015, 112, E5478–E5485. [Google Scholar] [CrossRef] [Green Version]
  39. Ballesteros, J.A.; Deupi, X.; Olivella, M.; Haaksma, E.E.J.; Pardo, L. Serine and threonine residues bend α-helices in the χ1=g− conformation. Biophys. J. 2000, 79, 2754–2760. [Google Scholar] [CrossRef] [Green Version]
  40. Baker, E.G.; Bartlett, G.J.; Crump, M.P.; Sessions, R.B.; Linden, N.; Faul, C.F.J.; Woolfson, D.N. Local and macroscopic electrostatic interactions in single α-helices. Nat. Chem. Biol. 2015, 11, 221–228. [Google Scholar] [CrossRef] [Green Version]
  41. Case, D.A.; Betz, R.M.; Cerutti, D.S.; Cheatham III, T.E.; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.; Goetz, A.W.; Homeyer, N.; et al. Amber 16; University of California: San Francisco, CA, USA, 2016. [Google Scholar]
  42. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  43. Hornak, V.; Abel, R.; Okur, A.; Strockbine, B.; Roitberg, A.; Simmerling, C. Comparison of multiple Amber force fields and development of improved protein backbone parameters. Proteins Struct. Funct. Bioinform. 2006, 65, 712–725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Hino, T.; Arakawa, T.; Iwanari, H.; Yurugi-Kobayashi, T.; Ikeda-Suno, C.; Nakada-Nakura, Y.; Kusano-Arai, O.; Weyand, S.; Shimamura, T.; Nomura, N.; et al. G-protein-coupled receptor inactivation by an allosteric inverse-agonist antibody. Nature 2012, 482, 237–240. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Newstead, S.; Kim, H.; von Heijne, G.; Iwata, S.; Drew, D. High-throughput fluorescent-based optimization of eukaryotic membrane protein overexpression and purification in Saccharomyces cerevisiae. Proc. Natl. Acad. Sci. USA 2007, 104, 13936–13941. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Strategy for state-targeting stabilization of GPCR, using de novo designed fusion partner proteins. (A) Multiple conformational states of A2AR. From left, the presented structures, respectively, correspond to the ones in an inactive, active-like intermediate, and activate states (PDB ID: 3PWH, 2YDV, and 5G53, respectively) [22,24,25]. The structure in the inactive state binds to the inverse agonist ZM241385, and the structures in the active-like intermediate, and active states bind to the agonist NECA. ZM241385 and NECA are shown in a sphere model. TM5 and TM6 are colored in blue for the inactive state, yellow for the active-like intermediate state, and pink for the active state. The loop connecting TM5 and TM6 is called ICL3. TM6 with TM5 exhibits large conformational changes upon the state transitions. (B) Our strategy for the state-targeting stabilization of A2AR. We fix the TM5 and TM6 conformation in a targeted state through the fusion strategy. To this end, we design α-helical proteins, which can be fused into A2AR in the targeted state through straight helical connections without kinks or intervening loops. In this work, we tested this idea by stabilizing the inactive state.
Figure 1. Strategy for state-targeting stabilization of GPCR, using de novo designed fusion partner proteins. (A) Multiple conformational states of A2AR. From left, the presented structures, respectively, correspond to the ones in an inactive, active-like intermediate, and activate states (PDB ID: 3PWH, 2YDV, and 5G53, respectively) [22,24,25]. The structure in the inactive state binds to the inverse agonist ZM241385, and the structures in the active-like intermediate, and active states bind to the agonist NECA. ZM241385 and NECA are shown in a sphere model. TM5 and TM6 are colored in blue for the inactive state, yellow for the active-like intermediate state, and pink for the active state. The loop connecting TM5 and TM6 is called ICL3. TM6 with TM5 exhibits large conformational changes upon the state transitions. (B) Our strategy for the state-targeting stabilization of A2AR. We fix the TM5 and TM6 conformation in a targeted state through the fusion strategy. To this end, we design α-helical proteins, which can be fused into A2AR in the targeted state through straight helical connections without kinks or intervening loops. In this work, we tested this idea by stabilizing the inactive state.
Ijms 22 12906 g001
Figure 3. Characterization of the de novo designed fusion partner proteins without A2AR. De novo designed fusion partner proteins, FiX1 (upper) and FiX2 (lower). (A) Computational models. (B) The energy landscape of each designed protein obtained from Rosetta folding simulations [17]. Each dot represents the lowest energy structure obtained from an independent trajectory starting from an extended chain (black) or the design model (red), the x-axis shows the Cα RMSD from the design model, and the y-axis shows the Rosetta full-atom energy. (C) Structural fluctuations of the N- and C-terminal helices in the MD simulation, starting from each design model. The main-chain RMSD of the N- and C- terminal helices of each snapshot structure during a MD trajectory against TM5 and TM6 of A2AR is shown along the time course. (Red and blue lines are the averaged and initial RMSD values, respectively.) (D) Far-ultraviolet CD spectra at various temperatures from 20 to 98 °C. (E) Two-dimensional 1H–15N HSQC spectra at 25 °C and 600 MHz (in parts per million, p.p.m). (F) Size-exclusion chromatograms combined with multi-angle light scattering (SEC-MALS) demonstrate that these designed proteins are monomeric in solution. M.W. stands for molecular weight.
Figure 3. Characterization of the de novo designed fusion partner proteins without A2AR. De novo designed fusion partner proteins, FiX1 (upper) and FiX2 (lower). (A) Computational models. (B) The energy landscape of each designed protein obtained from Rosetta folding simulations [17]. Each dot represents the lowest energy structure obtained from an independent trajectory starting from an extended chain (black) or the design model (red), the x-axis shows the Cα RMSD from the design model, and the y-axis shows the Rosetta full-atom energy. (C) Structural fluctuations of the N- and C-terminal helices in the MD simulation, starting from each design model. The main-chain RMSD of the N- and C- terminal helices of each snapshot structure during a MD trajectory against TM5 and TM6 of A2AR is shown along the time course. (Red and blue lines are the averaged and initial RMSD values, respectively.) (D) Far-ultraviolet CD spectra at various temperatures from 20 to 98 °C. (E) Two-dimensional 1H–15N HSQC spectra at 25 °C and 600 MHz (in parts per million, p.p.m). (F) Size-exclusion chromatograms combined with multi-angle light scattering (SEC-MALS) demonstrate that these designed proteins are monomeric in solution. M.W. stands for molecular weight.
Ijms 22 12906 g003
Figure 4. Experimental characteristics of A2AR fused with or without fusion partner proteins. (A) Crystal structures of A2AR WT (PDB ID: 3VG9) [22], A2AR–BRIL (PDB ID: 4EIY) [12], and the computational model of A2AR–FiX1. (B) (Top) Monomer bands in the clear-native PAGE for each A2AR sample heated at various temperatures. (Bottom) The fluorescence intensities of the gel images for each A2AR sample with temperature. The thermal transition from soluble to aggregated states was fitted (solid line) to obtain the midpoint temperature, Tm (aggr.) (C) Saturation binding curves of [3H]-ZM241385 (left) and [3H]-NECA (right) to A2AR WT (solid line), A2AR–BRIL (dashed line), and A2AR–FiX1 (dotted line). (D) Inhibition of [3H]-ZM241385 binding to A2AR WT (solid line) and A2AR–FiX1 (dotted line) by NECA. The binding of [3H]-ZM241385 in the absence of NECA was set to 100%. All measurements were carried out three times independently; dots show the average and whiskers show the s.e.m. for n = 3.
Figure 4. Experimental characteristics of A2AR fused with or without fusion partner proteins. (A) Crystal structures of A2AR WT (PDB ID: 3VG9) [22], A2AR–BRIL (PDB ID: 4EIY) [12], and the computational model of A2AR–FiX1. (B) (Top) Monomer bands in the clear-native PAGE for each A2AR sample heated at various temperatures. (Bottom) The fluorescence intensities of the gel images for each A2AR sample with temperature. The thermal transition from soluble to aggregated states was fitted (solid line) to obtain the midpoint temperature, Tm (aggr.) (C) Saturation binding curves of [3H]-ZM241385 (left) and [3H]-NECA (right) to A2AR WT (solid line), A2AR–BRIL (dashed line), and A2AR–FiX1 (dotted line). (D) Inhibition of [3H]-ZM241385 binding to A2AR WT (solid line) and A2AR–FiX1 (dotted line) by NECA. The binding of [3H]-ZM241385 in the absence of NECA was set to 100%. All measurements were carried out three times independently; dots show the average and whiskers show the s.e.m. for n = 3.
Ijms 22 12906 g004
Table 1. Solubilization efficiencies of A2AR fused with or without fusion partner proteins. Results are reported as mean ± standard deviation for n = 3 independent measurements.
Table 1. Solubilization efficiencies of A2AR fused with or without fusion partner proteins. Results are reported as mean ± standard deviation for n = 3 independent measurements.
ConstructSolubilization Efficiency (%)
A2AR WT24 ± 10
A2AR–BRIL57 ± 19
A2AR–FiX159 ± 21
A2AR–FiX242 ± 17
Table 2. Dissociation constants (Kd) of A2AR fused with or without fusion partner proteins by saturation binding assay. N.D. stands for not detected. Results are reported as the mean ± standard deviation for an n = 3 independent assay.
Table 2. Dissociation constants (Kd) of A2AR fused with or without fusion partner proteins by saturation binding assay. N.D. stands for not detected. Results are reported as the mean ± standard deviation for an n = 3 independent assay.
ConstructKd (nM)
ZM241385NECA
A2AR WT10.5 ± 0.3162 ± 44
A2AR–BRIL15.7 ± 0.6191 ± 22
A2AR–FiX17.3 ± 0.5N. D.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mitsumoto, M.; Sugaya, K.; Kazama, K.; Nakano, R.; Kosugi, T.; Murata, T.; Koga, N. State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein. Int. J. Mol. Sci. 2021, 22, 12906. https://doi.org/10.3390/ijms222312906

AMA Style

Mitsumoto M, Sugaya K, Kazama K, Nakano R, Kosugi T, Murata T, Koga N. State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein. International Journal of Molecular Sciences. 2021; 22(23):12906. https://doi.org/10.3390/ijms222312906

Chicago/Turabian Style

Mitsumoto, Masaya, Kanna Sugaya, Kazuki Kazama, Ryosuke Nakano, Takahiro Kosugi, Takeshi Murata, and Nobuyasu Koga. 2021. "State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein" International Journal of Molecular Sciences 22, no. 23: 12906. https://doi.org/10.3390/ijms222312906

APA Style

Mitsumoto, M., Sugaya, K., Kazama, K., Nakano, R., Kosugi, T., Murata, T., & Koga, N. (2021). State-Targeting Stabilization of Adenosine A2A Receptor by Fusing a Custom-Made De Novo Designed α-Helical Protein. International Journal of Molecular Sciences, 22(23), 12906. https://doi.org/10.3390/ijms222312906

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop