s2fieldsofstudy: Biology Chemistry Medicine

text:
Pharmacological evaluation of enantiomerically separated positive allosteric modulators of cannabinoid 1 receptor, GAT591 and GAT593

Positive allosteric modulation of the type 1 cannabinoid receptor (CB1R) has substantial potential to treat both neurological and immune disorders. To date, a few studies have evaluated the structure-activity relationship (SAR) for CB1R positive allosteric modulators (PAMs). In this study, we separated the enantiomers of the previously characterized two potent CB1R ago-PAMs GAT591 and GAT593 to determine their biochemical activity at CB1R. Separating the enantiomers showed that the R-enantiomers (GAT1665 and GAT1667) displayed mixed allosteric agonist-PAM activity at CB1R while the S-enantiomers (GAT1664 and GAT1666) showed moderate activity. Furthermore, we observed that the R and S-enantiomers had distinct binding sites on CB1R, which led to their distinct behavior both in vitro and in vivo. The R-enantiomers (GAT1665 and GAT1667) produced ago-PAM effects in vitro, and PAM effects in the in vivo behavioral triad, indicating that the in vivo activity of these ligands may occur via PAM rather than agonist-based mechanisms. Overall, this study provides mechanistic insight into enantiospecific interaction of 2-phenylindole class of CB1R allosteric modulators, which have shown therapeutic potential in the treatment of pain, epilepsy, glaucoma, and Huntington’s disease.

Introduction
The endogenous cannabinoid system (ECS) is composed of endogenous ligands (e.g. anandamide [AEA] and 2arachidonoylglycerol [2-AG]), anabolic and catabolic enzymes, and receptors (the predominant receptors being the type 1 and 2 cannabinoid receptors [CB1R, CB2R]) (Di Marzo, 2018;Estrada and Contereras, 2020). The ECS is ubiquitous in the human body . Both CB1R and CB2R are class A G protein-coupled receptors (GPCRs) (Lutz, 2020). CB1R is most-abundant in the brain and central nervous system (CNS) whereas CB2R is most-abundant in cells of the immune system (Iannotti et al., 2016;Hryhorowicz et al., 2019). The primary role of CB1R is to regulate mood, diet, sleep and pain sensation whereas CB2R regulates immune responses (Iannotti et al., 2016). Both receptors are activated by the endogenous ligands AEA and 2-AG (Carrera et al., 2020). Due to the ubiquitous nature of CB1R and CB2R, the ECS is considered a potential target for a wide array of diseases, but CB1R in particular is considered a target for the treatment of pain and neurological disorders such as Parkinson's disease, Huntington's disease, and epilepsy (Perez-Olives et al., 2021).
The intoxicating effects of Cannabis represent a major limitation to its use as medicine. The intoxicating effects of Cannabis are thought to be due to Δ 9 -tetrahydrocannabinol (THC) binding to the orthosteric site of CB1R (Slivicki et al., 2020). It has previously been hypothesized that the intoxicating properties of CB1R activation could be avoided if a drug bound to the allosteric site of the receptor promoted endogenous ligand activation without direct activation of CB1R because the endogenous cannabinoids AEA and 2-AG are not known to produce intoxicating effects, tolerance, or dependence (Slivicki et al., 2020). Allosteric modulators may be positive allosteric modulators (PAMs), negative allosteric modulators (NAMs), or silent allosteric modulators (Leo and Abood, 2021). A PAM enhances the effect of the primary ligand, a NAM reduces the effect of the primary ligand, and a silent allosteric ligand does not affect the pharmacology of the orthosteric ligand (Mielnik et al., 2021). The first described CB1R allosteric modulator was the indole carboxamide Org27569, which later drove the development of indole sulfonamides as potent CB1R NAMs (Greig et al., 2016;Dopart et al., 2018).
In a previous paper, we reported on the pharmacology of racemic CB1R allosteric modulators GAT591 and GAT593 (Figure 1), wherein we observed these ligands have both allosteric-agonist and positive allosteric modulator (ago-PAM) properties (Garai et al., 2020). We found that GAT591 and GAT593 act at allosteric sites due to their ability to enhance the binding of [ 3 H]CP55,940 at CB1R and at the same time modulate the receptor in absence of an orthostreric ligand (Garai et al., 2020). Biased agonism describes the ability of a ligand to preferentially activate one signalling pathway compared to another; for example G protein-versus βarrestin-mediated signaling (Leo and Abood, 2021;Patel et al., 2021). Our previous work has shown that the racemic mixtures GAT591 and GAT593 displayed bias towards Gα i/o signaling as compared to βarrestin (Garai et al., 2020), which earlier work indicated was correlated with improved cell viability (Laprairie et al., 2016(Laprairie et al., , 2019. Additional studies from our group support the idea that G protein bias of PAMs is correlated with improved outcomes in rodent models of Huntington's disease, pain, and absence epilepsy (Laprairie et al., 2019;Slivicki et al., 2020;Roebuck et al., 2021); whereas βarrestin bias may be correlated with reduced cell viability and increases pathology in animal models of absence epilepsy and Huntington's disease (Laprairie et al., 2016;Roebuck et al., 2020). One of the challenges in developing ligands for GPCRs is biased agonism (Leo and Abood, 2021;Patel et al., 2021). To further understand molecular mechanism(s) of action and probe for potential enantio specific interaction, we separated the enantiomers, determined their absolute stereochemistry and did biochemical characterization using cAMP and βarrestin2 assays to determine their allostericagonist and PAM activity. Previously we observed the distinct pharmacology between the enantiomers of (±)-GAT211, where (R)-GAT228 was an allosteric-agonist and the opposite enantiomer, (S)-GAT229, showed PAM activity at CB1R within the assays and cell line used . In this study we observed that the enantiomers of GAT591 and GAT593 display unique in vitro and in vivo effects that are likely to be associated with their unique modes of binding to CB1R.
2 Experimental section 2.1 In vitro evaluation 2.1.1 Compounds used CP55,940 and THC were purchased from Cayman (Ann Arbor, MI) and Sigma-Aldrich (Mississauga, ON), respectively. All other tested compounds were obtained from Dr. Ganesh Thakur, Northeastern University. All compounds were initially dissolved in dimethylsulfoxide (DMSO) and diluted in a 10% DMSO solution in phosphate-buffered saline (PBS). Compounds were added directly to cell culture at the time and concentrations indicated at a final concentration of 0.1% DMSO.
Frontiers in Pharmacology frontiersin.org 02 2.1.3 HitHunter cAMP assay cAMP inhibition was performed in the presence of 10 µM forskolin (FSK) using the DiscoveRx HitHunter assay in hCB1R CHO-K1 cells. Cells (16,000 cells/well in 96 well plates) were incubated overnight in Opti-MEM (Invitrogen) containing 1% FBS at 37°C and 5% CO 2 . Opti-MEM media was removed and replaced with cell assay buffer (DiscoveRx) and then cells were simultaneously treated with 10 µM FSK and experimental compounds (0.10 nM-10 µM) for 90 min cAMP antibody solution and working detection solutions were added according to the manufacturer's protocol (DiscoveRx), and cells were incubated for 60 min at room temperature. cAMP solution A was added according to the manufacturer's protocol (DiscoveRx) and cells were incubated for an additional 60 min at room temperature before chemiluminescence was measured on a Cytation5 plate reader (top read, gain 200, integration time 10 s). Data are presented as percent maximal CP55,940-dependent inhibition of cAMP accumulation.
2.1.4 PathHunter βarrestin2 assay βarrestin2 recruitment was determined using the hCB1R CHO-K1 cell PathHunter assay (DiscoveRx). Cells (16,000 cells/well in 96 well plates) were incubated overnight in Opti-MEM (Invitrogen) containing 1% FBS at 37°C and 5% CO 2 . Cells were then simultaneously treated with experimental compounds (0.10 nM-10 µM) for 90 min. Detection solution was added to cells according to the manufacturer's protocol (DiscoveRx), and cells were incubated for 60 min at room temperature. Chemiluminescence was measured on a Cytation5 plate reader (top read, gain 200, integration time 10 s). Data are presented as percent maximal CP55,940dependent stimulation.

[ 3 H]CP55,940 radioligand displacement assay
Previous work from our group has shown that CB1R PAMs enhance orthosteric agonist binding (e.g. CP55,940) Garai et al., 2020). Radioligand binding assays were carried out with 1 nM [ 3 H]CP55,940 in Tris buffer (75 mM Tris-HCl, 12.5 mM MgCl 2 , 1 mM ethylenediaminetetraacetic acid (EDTA), 1% bovine serum albumin (BSA), pH 7.4) with a total assay volume of 200 µl. The assay began with the addition of transfected hCB1R CHO-K1 cell membranes (25 µg protein per well). The assays were left to equilibrate at room temperature for 2 h before vacuum filtration using a Millipore Sigma 12-well sampling manifold and filter paper that had been soaked in wash buffer. Each reaction well was washed 3 times with a 2 ml aliquot of Tris-binding buffer. The filters were removed then submerged in 5 ml of scintillation fluid (Ultima Gold F, PerkinElmer, Buckinghamshire, United Kingdom). Radioactivity was quantified by liquid scintillation spectrometry. Specific binding was defined as the difference between binding that occurred in FIGURE 1 GAT591 previously studied was separated out into its respective enantiomers (GAT1664 and GAT1665) and GAT593 was separated out into its respective enantiomers (GAT1666 and GAT1667). The difference between GAT1664 and GAT1665 vs. GAT1666 and GAT1667 is the position of the fluorine in the indole ring.

[ 35 S]GTPγS assay
This assay was conducted as described in previous reports (Garai et al., 2020). To summarize, the assay was completed in the presence of [35S]GTPγS (0.1 nM), GDP (30 µM), GTPγS (30 µM) using membranes derived from CHO-K1 cells (1 mg/ ml) overexpressing hCB1R. Assay buffer consisted of 50 mM Tris, 10 mM MgCl2, 100 mM NaCl, 0.2 mM EDTA and 1 mM dithiothreitol (DTT) at pH 7.4. Membranes were incubated at 30°C for 90 min in a total volume of 500 ml. Reactions were ended by the addition of ice-cold wash buffer (50 mM Tris and 1 mg/ml BSA, pH 7.4) followed by rapid filtration under vacuum through pre-soaked Whatman GF/B glass-fibre filters in a 24well sampling manifold (Brandel Cell Harvester; Brandel Inc., Gaithersburg, MD, United States). Reaction wells were washed six times with a 1.2 ml aliquots of Tris-binding buffer. Filters were subsequently oven-dried for 60 min and placed in 3 ml of scintillation fluid (Ultima Gold XR, PerkinElmer). Bound radioactivity was determined by liquid scintillation counting. Basal [35S]GTPγS binding was determined in the presence of 20 mM GDP without any compounds present. Non-specific binding was determined in the presence of 10 mM GTPγS.

In vivo evaluation 2.2.1 Triad assessment in mice
Male C57BL/6 mice (Charles River, Senneville, QC) between 4 and 7 months of age were used for these studies. Animals were group housed at the Laboratory Animal Services Unit (LASU) at the University of Saskatchewan (3-5 animals/cage) with a standard 12:12 light-dark cycle, ad libitum access to food and water, and environmental enrichment. Compounds administered intraperitoneally (i.p.) were prepared in vehicle [ethanol and cremophor in saline (1:1:8)]. Catalepsy was assessed in the bar holding assay 5 min after drug administration. Mice were placed so that their forepaws clasped a 0.7 cm ring clamp positioned 4.5 cm above the surface of the testing space . The length of time the ring was held was recorded in seconds. The trial ended if the mouse turned its head or body or made 3 consecutive escape attempts to a maximum of 60 s. Body temperature was measured by rectal thermometer 15 min after drug administration. Anti-nociceptive effects were measured in the warm water (52 ± 2°C) tail-flick test 20 min after drug administration. Response in this case was defined by the removal of the tail from the warm water, with a maximal response time of 20 s. Catalepsy and tail flick data are presented as percent maximum possible effect (MPE). Compounds were administered at the doses indicated. Experimenters were blinded to treatment for all behavioral assessments and analyses. Animals were purchased, rather than bred, to reduce animal numbers. In all cases, experiments were performed with the approval of the University Animal Care Committee (UACC) at the University of Saskatchewan and are in keeping with the guidelines of the Canadian Council on Animal Care (CCAC) and the ARRIVE guidelines (Kilkenny et al., 2010).

In silico evaluation 2.3.1 Ligand Preparation
A conformational analysis was performed on GAT1664, GAT1665, GAT1666, and GAT1667 using Spartan'18 V1.4.5 (Spartan, 2018). A conformational search of each compound was run at ground state with molecular mechanics force fields (MMFF) (Spartan, 2018). Semi empirical calculations utilizing the PM3 Hamiltonian were carried out to calculate the equilibrium geometry of each conformer. After duplicates were eliminated, HF/6-31G* was applied to the remaining rotamers to obtain the global minimum energy conformation (Spartan, 2018).

CB1R model
The CB1R model has been described in detail in Hurst et al. (2019), but in brief: The CB1R model used in this study was based on the CB1R activated state crystal structure (PDB-ID: 5XRA) (Hua et al., 2017) and the cryo-EM CB1/Gi bound structure (PDB-ID: 6N4B) (Krishna Kumar et al., 2019). Structures were prepared with the Protein Preparation protocol (Suite 2019-1, Schrödinger, Inc.), mutations in the 5XRA structure were returned to wild-type, and both models were inspected for close contacts and crystal packing effects. Modification of the 5XRA structure to create the CB1R model involved calculating low free energy conformations for TMH2/7 to accommodate mutation data for S7.39, F2.61, F2.64 (Kapur et al., 2007;Shim et al., 2011). In addition, because of loop compression in the 5XRA structure, the IC1 loop from the 6N4B Gi bound structure was used in the model. TMH4 in the 5XRA CB1R* structure has tight crystal contacts on its IC end with TMH1 of another bundle because of antiparallel packing. For this reason, the conformation of TMH4 in the 6N4B Gi bound structure was used in the model. Modeller was used to remodel the EC3 loop (D6.58 to T7.33) and allow K (373) to interact with D2.63, consistent with mutation results (Fiser and Sali, 2003;Marcu et al., 2013). Modeller was also used to extend and model the N-terminus to S (88) with inclusion of the C (98) to C (107) disulfide bridge Farrens, 2013, 2015). The receptor/ligand complex was energy minimized in Prime (Suite 2019-1, Schrödinger, Inc.). The Prime implicit membrane functionality was employed. Hydrophilic residues facing the binding crevice and within the low dielectric region of the implicit membrane were excluded from the low dielectric via exclusion spheres placed on each residue. The Generalized Born/Surface Area (GB/SA) continuum Frontiers in Pharmacology frontiersin.org 04 solvation model for water was used with the dielectric set to 80 outside of the implicit membrane region and 2 within. A truncated Newton conjugate gradient minimization was performed using the OPLS3e force field for one iteration up to a maximum of 1,000 steps and with a 0.1 kcal/mol gradient endpoint. Constraints of 1 kcal/mol were placed on the C-alpha atoms of residues R3.50, Y5.58, L6.33, and Y7.53 to prevent the intracellular opening present in the R* structures from closing during the minimization.

Binding Site Identification
To identify potential binding site(s) for GAT1664, GAT1665, GAT1666 and GAT1667 at CB1R, we used the Forced-Biased Metropolis Monte Carlo simulated annealing program (MMC) as the first step (Guarnieri and Mezei, 1996;Clark et al., 2006). The MMC method has been used successfully to identify water binding sites in proteins and on DNA to identify potent and novel p38 kinase inhibitors; to identify thermolysin and T4 and lysosome biding sites; and, to identify the binding site of the negative allosteric modulator, pregnenolone (Morales et al., 2016;Zhang et al., 2020). In MMC, the molecule of interest is first divided into smaller fragments. A series of grand canonical ensembles of a molecular fragment interacting with the protein in a large simulation box are created. The Chemical potential of the system is then annealed at descending chemical potential levels, with each new level starting from the last ensemble generated from the previous one. At each step, fragment ligand poses are sampled throughout the box and over the entire protein. Fragments are treated as rigid solvents that are inserted and deleted millions of times until the lowest energy configuration is found at the explored annealing level. As the chemical potential is annealed, the number of fragments in the box decreases because the method eliminates any fragment that has a poorer free energy of interaction than the annealing level. The output of the calculation is an ensemble of ligand poses at each chemical potential level in the annealing schedule. The method is repeated for each molecular fragment. Data analysis is performed using the GENS tool in the MMC program available from the Mihaly Mezei Laboratory (Mezei, 2011).
Of particular interest in MMC calculations are fragments that persist at particular sites on the protein throughout the annealing schedule, since these fragments clearly have high affinity for those sites. Receptor regions at which all molecular fragments of a studied ligand collect are identified from MMC output. This set of sites is then refined to include only those sites at which the order of the fragments reflects the structure of the entire molecule.
GAT1664 and GAT1665 were broken into fragments that included a 6-fluoro-indole ring, a 2-fluoro-phenyl ring and a 3fluoro-phenyl ring. GAT1666 and GAT1667 were broken into fragments that included 7-fluoro-indole ring, 2-fluoro-phenyl ring and 3-fluoro-phenyl ring. All fragments were prepared with partial charges using the Amber 2002 force field, a point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum MMC (Wang et al., 2018). Six MMC runs were performed in which our CB1R receptor model was immersed in a box filled with copies of one of these fragments. Analysis of the MMC runs for GAT1664 and GAT1666 revealed that while each fragment bound to multiple positions on the CB1R receptor, there was only one region in which all fragments clustered. This was at the extracellular end of TMH2/3, just beneath the EC1 loop. Y2.59, a polar residue that faces lipid was found to attract the nitro group fragment, while D2.63 consistently was found interacting with the fluorinated indole fragments.
Docking: Once the two sites were confirmed, the structures of GAT1664-GAT1667 were docked in the intracellular TMH1/2/ 4 site with H2.41, F4.46, and W4.50 used as direct interaction sites. GAT1664 and GAT1666 were docked in the TM2/3/ EC1 site with Y2.59 and D2.63 as direct interaction sites. In addition, after extensive previously published molecular dynamics calculations on unfluorinated versions of the GAT compounds (GAT1600-3), R (148) was modeled to interact directly with the carbonyl oxygens of the last turn of TMH1 and not the nitro group within the compounds . The receptor/ligand complexes were energy minimized in Prime (Suite 2019-1, Schrödinger, Inc.). The Prime implicit membrane functionality was employed. Hydrophilic residues facing the binding crevice and within the low dielectric region of the implicit membrane were excluded from the low dielectric via exclusion spheres placed on each residue. The Generalized Born/Surface Area (GB/SA) continuum solvation model for water was used with the dielectric set to 80 outside of the implicit membrane region and 2 within. A truncated Newton conjugate gradient minimization was performed using the OPLS3e force field for 1 iteration, up to a maximum of 1,000 steps and with a 0.1 kcal/mol gradient endpoint. Constraints of 1 kcal/mol placed on the c-alpha atoms of residues R3.50, Y5.58, L6.33, and Y7.53 were set to prevent the intracellular opening from closing during the minimization. The resulting docks were refined with the Induced Fit protocol (Suite 2019-1, Schrödinger, Inc.). The Glide box size was set to 12 Å3 centered on the ligand and the SP docking algorithm employed. Residues within 5 Å of the docked ligand were included in the Prime refinement stage, except in the case of the TM2/3 PAM site where S2.60 and K3.28 were excluded based on mutation data (Song and Bonner, 1996;Kapur et al., 2007). The implicit membrane previously used during the initial Prime minimization was employed here as well.

MMGBSA analysis
Each of the GAT1664-GAT1667 ligand-receptor complexes were evaluated via MMGBSA. This was the scoring function used to evaluate each complex and not Glide scores. This method is used to estimate the relative binding affinities for a list of ligands Frontiers in Pharmacology frontiersin.org (reported in kcal/mol). As the MMGBSA binding energies are approximate free energies of binding, a more negative value indicated stronger binding. For this computation, the VSGB solvation model was employed while the chosen force field OPLS4 was used.

Statistical analysis
Data related to EC 50 and E max in Figures 3, 4 were obtained with nonlinear regression models (3-parameter model, GraphPad Prism 9.0, San Diego, CA). Results were calculated as percent response relative to the reference agonist CP55,940 (Figure 3), percent [ 3 H]CP55,940 bound ( Figure 4A), or percent stimulation above baseline (Figures 4B,C). The 3-parameter model was used for [ 3 H]CP55,940 binding assays ( Figure 4A) because standard radioligand competition models that fit for K d do not appropriately account for allosteric interactions. Significance was determined by one-or two-way ANOVA followed by Tukey's or Bonferroni's post-hoc test as indicated. p < 0.05 was considered significant and p values were only employed where ANOVA was used. Data related to Figure 5 was analyzed in GraphPad Prism 9.0 was used to analyze in vivo data and p < 0.05 was considered to be statistically significant. Data are expressed as mean ± SEM. Group sizes for all experiments are described in figure legends.

Results and discussion
Our in vitro study explored CB1R-dependent modulation of cAMP inhibition and βarrestin2 recruitment in CHO-K1 cells stably-expressing hCB1R with or without CP55,940. Inhibition of cAMP accumulation is Gα i/o -protein mediated whereas recruitment of βarrestin2 is G-protein independent (Hryhorowicz et al., 2019). In a previous paper, GAT591 and GAT593 displayed mixed agonist and PAM (i.e. ago-PAM) activity (Garai et al., 2020). The compounds GAT591 and GAT593 contain a 50:50 racemic mixture of both R and S enantiomers. In this paper each enantiomer was investigated to determine their allosteric agonist and PAM activity, as we previously observed these properties to be enantiomerically distinct in the compounds GAT211, GAT228, and GAT229 within those experimental conditions . Separating GAT591 and GAT593 into their respective enantiomers gave GAT1664, GAT1665, GAT1666 and GAT1667 (Figure 1).

Synthesis, chiral separation and absolute stereochemistry determination
We synthesized GAT591 and GAT593 on a multigram scale using our previously published method (Garai et al., 2020). Both enantiomers of each compound were separated in high optical purity (>99%) using superfluid chiral high performance liquid chromatography (HPLC). The absolute stereochemistry of (+)-GAT1664 and (+)-GAT1666 was determined by singlecrystal X-ray diffraction technique (for details see supporting formation) and was found to be "S" for both of these enantiomers (Figure 2). Based on this study, we can predict the absolute stereochemistry of each opposite enantiomers, (-)-GAT1665 and (-)-GAT1667 as "R" (Figures 1, 2).

In vitro evaluation
Previously, the racemic compounds GAT591 and GAT593 were characterized for their PAM activity in the presence of 100 nM CP55,940 (Garai et al., 2020). Here, CP55,940 displayed EC 50 values of 15 nM and 310 nM in the cAMP inhibition and βarrestin2 recruitment assays respectively ( Figure 3). Therefore, 100 nM CP55,940 was used in PAM assays in the current study for consistency between studies and as an intermediate concentration.
All enantiomers displayed G protein-mediated agonist activity at hCB1R as shown by their ability to inhibit cAMP production in the absence of an orthosteric agonist ( Figure 3A). The ability of these compounds to augment orthosteric agonist signaling was tested in the presence of 100 nM CP55,940 because we have used this concentration of CP55,940 in previous studies (Garai et al., 2020) and because 100 nM CP55,940 produce an~65% response relative to the maximum response observed prior to fitting the data to a nonlinear regression (Figure 3). All 4 compounds increased CP55,940's activity above that observed at the 100 nM level, consistent with either positive allosteric modulation or additive non-competitive agonism (i.e. ago-PAMs) of hCB1R in the presence of the orthosteric agonist CP55,940 in the cAMP assay ( Figure 3C). As an agonist, GAT1666 displayed greater efficacy than GAT1664, and GAT1667 displayed the greatest efficacy at inhibiting cAMP as a PAM, although this was not statistically different from other compounds tested. All enantiomers displayed low nanomolar potency as agonists ( Figure 3A). Although not statistically significant differences, both Renantiomers (GAT1665 and GAT1667) displayed greater PAM potency and efficacy than the S-enantiomers (GAT1664 and GAT1666), ( Figure 3C). Initial studies with the parent compound scaffold, GAT211, found that PAM activity was associated with GAT229 whereas allosteric agonist activity was associated with GAT228 . Importantly, repeated testing of GAT229 and its structural analogs over the years has shown assay-specific variability and that these enantiomers can display some allosteric agonist activity (Garai et al., 2020), including in the present study. Pure PAM activity has been observed for GAT229 in the autaptic hippocampal neuron ex vivo model systems (Mitjavila et al., 2018). These observations are in keeping with our earlier findings with the parent racemic ligands because PAM activity is predominantly attributable to GAT229, which shares the same spatial orientation as GAT1665 and GAT1667, and agonist activity is predominantly attributable to GAT228, which shares the same spatial orientation as GAT1664 and GAT1666 .
Frontiers in Pharmacology frontiersin.org 07 For consistency with cAMP inhibition experiments and previous studies (Garai et al., 2020), compounds were tested for their ability to augment CP55,940-dependent βarrestin2 recruitment with 100 nM CP55,940. Both GAT1665 and GAT1667 recruited βarrestin2 as low-potency agonists and PAMs at hCB1R (Figures 3B,D). GAT1664 did not recruit βarrestin2 to hCB1R while GAT1666 displayed a weak ability to do so as an agonist (Figures 3B,D). This suggests that as agonists, GAT1665 and GAT1667 (the R-enantiomers) recruit βarrestin2.
The operational model described by Kenakin et al. (2012) cannot be used to describe bias for non-competitive interactions . Data are expressed as % radioligand bound (A) or % stimulation above baseline (i.e. vehicle) levels (B,C). Data were fit to a nonlinear regression (3-parameter model, GraphPad v. 9.0) to determine EC 50 , and E max . Data are mean ± SEM (E max ) or 95% CI (EC 50 ); n = 6 independent experiments performed in triplicate, except for CP55,940 where n = 12 independent experiments performed in triplicate as this compound was included in all assays as a reference compound. ***p < 0.001 compared to CP55,940,^p < 0.001 compared to GAT1664, $$$ p < 0.01 compared to GAT1665, ### p < 0.01 compared to GAT1666, as determined by one-way ANOVA followed by Tukey's post-hoc test.
Frontiers in Pharmacology frontiersin.org 08 accurately. Therefore, bias (ΔΔlogR) cannot be directly estimated between cAMP inhibition and βarrestin2 recruitment. Examining the potency and efficacy of compounds directly, however, we observed that as an agonist, GAT1664 displayed no activity βarrestin2 recruitment and this compound appears to heavily favour G protein-dependent inhibition of cAMP (Figure 3). GAT1666 displayed greater potency in the cAMP inhibition assay (Figure 3). GAT1665's potency and efficacy generally favored cAMP inhibition (Figure 3). GAT1667 did appear to favour inhibition of cAMP as an agonist ( Figure 3). As PAMs, GAT1664 and GAT1666 displayed no activity in the βarrestin2 recruitment assay (Figure 3). GAT1665 and GAT1667 favoured cAMP inhibition relative to βarrestin2 recruitment ( Figure 3F). The parent compound GAT229 shares the same spatial orientation as GAT1665 and GAT1667 and also displayed greater potency in cAMP inhibition assays when tested as a PAM . Therefore, these data confirm previous observations that this scaffold's
Frontiers in Pharmacology frontiersin.org 09 spatial orientation promotes G protein PAM activity and extends earlier reports to demonstrate increased potency and efficacy of fluorine-substituted ligands as compared to parent compounds .
Next, we sought to assess the ability of these compounds to modulate orthosteric agonist binding to hCB1R. All enantiomers studied augmented [ 3 H]CP55,940 binding, consistent with their proposed PAM activity, with the R-enantiomers GAT1665 and to a greater extent GAT1667 increasing [ 3 H]CP55,940 binding ( Figure 4A). GAT1667 displayed the greatest potency among the enantiomers tested (EC 50 = 28 [15-52] nM). Enantiomers were also tested for their ability to promote G protein coupling in the GTPγS assay ( Figures 4B,C). All of the enantiomers tested increased G protein coupling alone (i.e. in the absence of the orthosteric agonist CP55,940) with the R-enantiomers GAT1665 and GAT1667 being more potent than the S-enantiomers and GAT1667 in particular displaying the greatest potency (EC 50 = 8.0 [3.4-17] nM) ( Figure 4B). When 1 μM of each enantiomer was tested in the presence of CP55,940, each enantiomer increased G protein coupling which may be the result of allosteric non-competitive agonism or PAM activity ( Figure 4C). Therefore, all enantiomers tested were able to stimulate G protein coupling alone or in the presence of CP55,940, which is consistent with their activity as ago-PAMs. Among the enantiomers tested, GAT1667 consistently displayed the greatest potency and efficacy.

In vivo evaluation
All enantiomers were evaluated in male C57BL/6 mice using a triad of outcomes consisting of catalepsy, body temperature, and nociception. Previous studies have shown that the racemic GAT591 and GAT593 did not produce catalepsy or hypothermia at 0.1, 1, 3, or 10 mg/kg (i.p.) compared to vehicle (Garai et al., 2020). Both GAT591 and GAT593 however did produce a dose-  Frontiers in Pharmacology frontiersin.org 10 dependent anti-nociceptive effect in the tail flick assay that was significant relative to vehicle at 3 and 10 mg/kg (Garai et al., 2020). In this study we evaluated all separated optically pure enantiomers in a triad assay. When tested alone, GAT1664, GAT1665, GAT1666 and GAT1667 did not produce catalepsy, hypothermia, or anti-nociceptive effects at 0.1, 1, 3, or 10 mg/kg i.p. compared to vehicle (Figures 5A-C). We also tested these enantiomers as PAMs by co-administering a sub-threshold dose of 1 mg/kg THC with each GAT compound. THC was selected as a common CB1R partial agonist that could produce a moderate in vivo response to be modulated by our compounds of interest. It was found that GAT1665 and GAT1667 displayed significant cataleptic, hypothermic, and anti-nociceptive effects at 3 and 10 mg/kg i.p. compared to vehicle (Figures 5D-F). Therefore, GAT1664, GAT1665, GAT1666, and GAT1667 did not affect animal responses to catalepsy, body temperature, or nociception in the warm water tail flick assay in normal, otherwise healthy adult mice as agonists; but GAT1665 and GAT1667 were able to augment THC's effects in vivo when administered once. For comparison, the racemic parent compounds, GAT591 and GAT593, do evoke anti-nociceptive responses in vivo when administered alone Garai et al., 2020).  Figure 6. The conformation of GAT1664 before it binds to hCB1R (B1) and then the conformation of GAT1664.after it binds to hCB1R (B2). (C) GAT1666 main interactions and the amino acid residues of the binding pocket displayed in tan in Figure 6. The conformation of GAT1666 before it binds to hCB1R (D1) and then the conformation of GAT1664 after it binds to hCB1R (D2).

In silico studies
Compounds were docked to a modelled structure of hCB1R to determine their putative binding sites within the receptor. GAT1664 and GAT1666 displayed the greatest modelled affinity to an allosteric site on hCB1R proximal to the intracellular face of the receptor in tan (Figures 6A,B). An H-bond exists between the indole hydrogen in white and the nitrogen atom in blue (H-N distance 1.96 Å) on H154 2.41 ( Figure 6A). One of the fluorobenzyl groups on GAT has a face-to-face π-π interaction with F237 4.46 (4.03Å) while the other fluorobenzyl has a T-shaped π-π interaction with each ring on W241 4.50 (4.96 Å to the 6-membered ring; 5.09 Å to the 5-membered ring) ( Figure 6A). Upon binding to hCB1R the conformation of the ortho-fluorobenzyl substituent on GAT1664 rotates 90°. The conformational cost for this rotation is 2.35 kcal/mol (Table 1). GAT1666 has the same binding pose, thus the binding energy was measured to look at FIGURE 8 GAT1665 and GAT1667 bind to a putative PAM site. (A) GAT1665 main interactions and the amino acid residues of the binding pocket displayed in pink in Figure 6. The conformation of GAT1665 before it binds to hCB1R (B1) and then the conformation of GAT1665 after it binds to hCB1R (B2). (C) GAT1667 main interactions and the amino acid residues of the binding pocket displayed in pink in Figure 6. The conformation of GAT1667 before it binds to hCB1R (D1) and then the conformation of GAT1667 after it binds to hCB1R (D2).
Frontiers in Pharmacology frontiersin.org the relative affinity. GAT1666 is more stable than GAT1664 by 2.31 kcal/mol (ΔG GAT1666-ΔG GAT1664) which suggests that it has a higher affinity to the allosteric site of hCB1R shown in tan ( Figure 6). This can be attributed to changing the fluorine on the sixth position of the indole (GAT1664) to the seventh position on the indole (GAT1666). The binding of a compound to this site may allow it to act as an agonist by promoting the X1 = g + → trans conformation of residue F4.46, facilitating breaking of the R3.50/D6.30 ionic lock which has been shown to activate CB1R (Figures 6, 7) (Hurst et al., 2019). GAT1665 and GAT1667 displayed the greatest modelled affinity to a distinct allosteric site on hCB1R relative to GAT1664 and GAT1666 (Figures 6, 8). The NO 2 group on GAT1665 has a H-bond between 1 of the oxygens shown in red and a hydrogen atom on Y172 2.59 (O-H distance 2.17 Å) (Figure 8). Another H-bond exists between the indole hydrogen in white and an oxygen atom in red on D176 2.63 (H-O distance 2.29 Å) (Figure 8). A cation-π interaction exists between R182 EC1 and each ring on the indole on GAT1665 (six-membered ring centroid on indole to NH 2 + on R182 EC1 distance 3.59 Å; five-membered ring centroid on indole to NH 2 + on R182 EC1 distance 3.66 Å) (Figure 8). One of the fluorobenzyl groups on GAT has a face-to-face π-π interaction with Y172 2.59 (4.05 Å) (Figure 8). Upon binding to hCB1R the conformation of the NO 2 rotates in order to H-bond with Y172 2.59 (Table 1). The conformational cost for this rotation is 3.27 kcal/mol. If we compare the binding energy of GAT1665 and GAT1667 it is apparent that GAT1667 is more stable than GAT1665 bound to the PAM site by 4.78 kcal/mol (ΔG GAT1667-ΔG GAT1665). This can be attributed to changing the fluorine on the sixth position of the indole (GAT1665) to the seventh position on the indole (GAT1667). The previously characterized compound GAT229 also binds to this site and lacks intrinsic activity thus both these compounds alone cannot induce the active state of hCB1R and thus are PAMs. The R-enantiomers (GAT1665 and GAT1667) undergo a smaller conformational cost than the S-enantiomers (GAT1664 and GAT1666) upon binding to their respective sites on hCB1R (Table 1). The R-enantiomers (GAT1665 and GAT1667) bind to the ago-PAM site (ΔG = −60.77 and −64.55 kcal/mol) significantly weaker than the S-enantiomers (GAT1664 and GAT1666) bind to the PAM site (−69.68 and −71.99 kcal/mol). Of note, binding of enantiomers to either site was not absolute but was more energetically favourable for the S-enantiomers to bind near the intracellular surface and the R-enantiomers to bind near the first extracellular loop.

Conclusion
In summary, the purified enantiomers of GAT591 and GAT593 displayed unique pharmacology. The purified enantiomers GAT1664, GAT1665, GAT1666, and GAT1667 bound to hCB1R in the low nM range and were shown to enhance binding of [ 3 H]CP55,940 at hCB1R. Enhancing the binding of an orthosteric ligand is a characteristic of a PAM and these observations are consistent with all compounds operating in part via positive allostery and/or non-competitive allosteric agonism (Alaverdashvili and Laprairie, 2018;Kenakin and Strachan, 2018). Future explorations of these compounds will utilize full concentration-response curves of these compounds with and without CP55,940 or an endogenous orthosteric agonist (e.g. 2-AG). Such more in-depth experiments will allow for the use an operational model of allosterism that quantifies co-operativity effects on potency, efficacy, and bias as our group has done previously for both PAMs and NAMs (Laprairie et al., 2016;Laprairie et al., 2017); and determine whether the tested compounds display probe-dependence.
In addition to this, we demonstrated that these enantiomers activate hCB1R-dependent inhibition of cAMP and βarrestin2 recruitment. All of the GAT compounds tested here inhibited cAMP accumulation as allosteric agonists and PAMs (i.e. ago-PAMs), with GAT1665 and GAT1667 displaying the greatest potency and efficacy as CB1R PAMs in the presence of 100 nM CP55,940. In the βarrestin2 assay, GAT1665 and GAT1667 (R-enantiomers) acted as ago-PAM's while GAT1664 and GAT1666 (S-enantiomers) showed little to no activity as either an agonist or PAM. In the G protein coupling assay, the R-enantiomers GAT1665 and GAT1667 also showed higher potency and efficacy than S-enantiomers as agonists and PAMs. As for the in vivo studies, no compounds tested produced catalepsy, hypothermia at or antinociceptive effects at 0.1, 1, 3 or 10 mg/kg i.p. compared to vehicle when tested alone, but GAT1665 and GAT1667 did potentiate the effects of THC as ago-PAMs. These data are congruent with previous observations that CB1R ago-PAMs are inactive or minimally active in vivo in wild-type, otherwise healthy animals under non-pathological and acute conditions (Alaverdashvili and Laprairie, 2018). In planned future studies, these ago-PAMs will be further assessed in rodent models of absence epilepsy and pain to determine whether they are able to augment the endogenous cannabinoid signaling as we observed previously with related allosteric ligands (Slivicki et al., 2020;Roebuck et al., 2021). In addition, the magnitude of observed in vivo effects may have been effected by pharmacokinetics of our compounds; planned future studies will assess and optimize the pharmacokinetics of CB1R allosteric ligands. Our in silico data support in vitro binding and signaling data that indicate these compounds interact with an allosteric site on CB1R that is distinct from the primary orthosteric agonist binding site (Price et al., 2005;Vigolo et al., 2015;Hua et al., 2017;Krishna Kumar et al., 2019). GAT1664 and GAT1666 preferentially bind to a putative allosteric agonist site of hCB1R near the intracellular face of the receptor between transmembrane helices I and IV; while GAT1665 and GAT1667 preferentially bind to a putative PAM site of CB1R on the extracellular receptor surface. Enantiomers may display low affinity to multiple allosteric sites, accounting for the allosteric agonist activity of GAT1665 and GAT1667 (Price et al., 2005;Kapur et al., 2007;Hurst et al., 2019). However, future dynamic modeling and crystallization studies will be required to Frontiers in Pharmacology frontiersin.org fully understand how the putative allosteric sites contribute to partial agonism and positive allosteric modulation, as the current models represent static interactions. Collectively, these data support retained, potent, G protein-selective, ago-PAM activity at hCB1R for the GAT211 ligand scaffold and demonstrate enhanced potency and efficacy of these ligands relative to the parent compound Garai et al., 2020). Moreover, recent data from the racemic mixtures of these enantiomers indicates fluorine addition improves both metabolic stability and blood brain barrier penetrance; with GAT593-and therefore GAT1667-being superior to GAT591 (Garai et al., 2020). Therefore, studies are now underway to assess the in vivo efficacy of GAT1667 in the contexts of absence epilepsy and pain during both acute and chronic treatment paradigms.

Permission to reuse and copyright
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.

Ethics statement
The animal study was reviewed and approved by University Animal Care Committee Animal Research Ethics Board, University of Saskatchewan.


s2fieldsofstudy: Biology Medicine

text:
Glioma-associated mesenchymal stem cells-mediated PD-L1 expression is attenuated by Ad5-Ki67/IL-15 in GBM treatment

Background Glioblastoma (GBM) is a highly immunosuppressive and vascular malignant brain tumor. Current therapeutic strategies targeting tumor cells have limited efficacy because of the immunosuppressive microenvironment and vascularization. Glioma-associated mesenchymal stem cells (GA-MSCs) have been identified as important stromal components of the tumor microenvironment, owing to their contribution to tumor angiogenesis and their potential to drive glioma stem cells. However, there are no reports on the effect of oncolytic Ad5-Ki67/IL-15 on programmed death ligand 1 (PD-L1) expression and angiogenesis induced by GA-MSCs. Methods Flow cytometry was respectively performed to detect the PD-L1 of glioma cells and programmed death protein 1 (PD-1), CD3, CD4 and CD8 in lymphocytes, as well as distribution of the cell cycle. CCK-8 assay investigated the proliferation of glioma cells and GA-MSCs in vitro. Tumor-bearing nude mice were established with U87-Luc cells and treated with the viruses, and further the IVIS spectrum was utilized to obtain luciferase images. Finally, the expression of PD-L1 in tumor tissues was also investigated using western blotting. Results We found that GA-MSCs had potential to induce PD-L1 upregulation and involved in vascular mimicry in vitro. Importantly, Ad5-Ki67/IL-15 reduced PD-L1 expression of glioma cells and neovascularization by targeting GA-MSCs. Furthermore, despite the presence of GA-MSCs, the virus has the ability to generate potent antitumor efficacy in vitro and vivo. Conclusions These findings suggest the use of oncolytic Ad5-Ki67/IL-15 targeting GA-MSCs to treat GBM, indicating potential clinical applications. Supplementary Information The online version contains supplementary material available at 10.1186/s13287-022-02968-z.

Background
Glioblastoma (GBM) is the most aggressive and lethal brain tumor, and patients with GBM have extremely poor prognosis and an expectancy time of almost 14.6 months [1,2]. Current treatments for glioma include radiation and chemotherapy, for instance, temozolomide and bevacizumab and aim to kill tumor cells and inhibit tumor neovascularization [3][4][5][6]. However, there have been few studies of treatment strategies targeting surrounding stromal cells such as mesenchymal stem cells (MSCs). Glioma-associated MSCs (GA-MSCs) have been identified as stromal components that contribute to tumorigenesis and have the potential to drive glioma stem cells (GSCs), particularly in the unique microenvironment of human brain tumors [7,8]. GA-MSCs are important mesenchymal cells in glioma microenvironment. Furthermore, the fraction of GA-MSCs in human high-grade gliomas is inversely correlated with the overall survival time of patients [9]. In addition, the function of GA-MSCs in glioma progression is associated with CD90 expression [10]. However, the immunomodulation roles of GA-MSCs remain to be elucidated.
Recent studies reveal that bone marrow mesenchymal stem cells (BM-MSCs) have crucial roles in determining immunosuppressive characteristics [11][12][13][14]. BM-MSCs have been reported to secrete immunosuppressive factors including prostaglandin E2 (PGE2), transforming growth factor beta (TGFβ), and interleukin-10 (IL-10), as well as intracellular enzyme indoleamine 2,3-dioxygenase (IDO), resulting in suppression of proliferation, activation and differentiation of T cells [15]. In addition, these mesenchymal stem cells (MSCs) produce various chemokines and express adhesion molecules that are responsible for immune cell recruitment and maintaining close proximity with immune cells [16,17]. Furthermore, the suppressive effect of MSCs on T cells is thought to be mediated by the release of TGFβ and hepatocyte growth factor (HGF), leading to a reduction in cyclin D2 and an increase in p27kip1 levels in T cells, with consequent arrest of the cell cycle in the G1 phase [18,19]. Nevertheless, the roles of GA-MSCs in glioma treatment have not well expounded. Therefore, deeper understanding of the immunological roles of GA-MSCs will facilitate the development of new strategies for the treatment of gliomas.
Programmed death ligand 1 (PD-L1), which is present in the tumor microenvironment (TME), is an immune inhibitory receptor ligand that leads to immune cell dysfunction and apoptosis by binding to its receptor, programmed death protein 1 (PD-1), which works in braking the inflammatory response and conspiring tumor immune evasion [20,21]. PD-1, another immune checkpoint receptor, is expressed on activated immune cells, which usually upregulated in the TME [22,23]. In addition, GBM is a malignant brain disease that abundant with vascular structure. Tumor angiogenesis is one of the important factors that lead to rapid growth and progression of glioma. Thus, it is critical to develop strategies that will enhance immune response during treatment and expand the range of brain gliomas that can be effectively treated.
Oncolytic virus (OV) therapy is a promising therapeutic approach for solid tumors that involves selectively infecting and killing tumor cells [24,25]. Most previous studies of OVs have focused on tumor cells, but the research about OVs targeting GA-MSCs have not been elucidated. Therefore, we designed a novel oncolytic adenovirus Ad5-Ki67/IL-15 to explore its potential therapeutic applications for targeting MSCs. Our previous work demonstrated that Ad5-Ki67/IL-15 selectively killed tumor cells and exhibited potent antiangiogenic capacity via reduction of VEGF secretion [25].
Here, we reveal that GA-MSCs contribute to PD-L1 upregulation in GBM. Furthermore, Ad5-Ki67/IL-15 can reduce GA-MSC-mediated PD-L1 expression and angiogenesis in glioma. In GBM model, the virus can effectively inhibit tumor growth even though the presence of GA-MSCs. This research is innovative because the immunosuppressive effects of GA-MSCs have not previously been explored or tested. These results indicate a potential new approach for the treatment of glioma, using Ad5-Ki67/IL-15 to target not only tumor cells but also GA-MSCs.

Isolation and culture of GA-MSCs
Isolation of human GA-MSCs was performed as described previously [10]. Briefly, the GA-MSCs were separated from the fresh glioma specimens; the tissue was cut into pieces after washed twice in PBS. Afterward, the trypsin was added to the pieces for digestion. Then, the cell suspension was transferred to a 70-μm filter and washed twice again with PBS. GA-MSCs were cultured in DMEM containing 10% fetal bovine serum (FBS; Gibco, USA) and 100 U/ml penicillin/streptomycin (Gibco, USA) in a humidified atmosphere at 37 °C with 5% CO 2 . All experiments were performed according to the flow chart (Additional file 1).

Mouse lymphocytes
Isolation of mouse spleen lymphocytes was performed according to the instructions. Briefly, fresh spleen was obtained and transferred to a 70-μm filter (Corning, USA) on a 50-ml centrifuge tube, then gently ground, washed continuously with precooled PBS and centrifuged at 1800 rpm for 5 min. The cell pellets were re-suspended in washing solution and separated with lymphocyte fluid at a 1:1 ratio density gradient centrifugation at 2000 rpm for 20 min. The collected cells were washed using washing buffer and cleaning buffer before use in experiments.

Collection of conditioned media
GA-MSCs, U251 and BT-01 cells were seeded in six-well plates containing the viruses at a multiplicity of infection (MOI) of 40 for 72 h. The supernatants were collected and centrifuged at 2000 rpm for 10 min to remove cells and cellular debris.

Immunofluorescence
GA-MSCs were fixed with 4% paraformaldehyde, permeabilized with 0.5% Triton, and blocked with 3% bovine serum albumin for 1 h at room temperature. The cells were incubated with primary antibodies against CD44 and CD105 (1:100, Proteintech, China), then overnight at 4 °C. The secondary antibodies used were as follows: Cy3-conjugated antibodies (1:100, Boster, WuHan, China). Nuclei were stained for visualization using DAPI. Immunofluorescence microscopy was performed with an Olympus microscope.

Fluorescence microscopy
GA-MSCs were treated with Ad5-Ki67/GFP at an MOI of 40 and were observed under an Olympus microscope. Images were taken 48 h after infection.

Cell cycle
The GA-MSCs (1 × 10 5 cells/well) were seeded into 6-well plate with 10% FBS (Invitrogen, China) and 100 U/ ml penicillin/streptomycin (Gibco, USA) in a humidified atmosphere at 37 °C with 5% CO 2 . After the cell adherence, the culture condition was added with the virus Ad5-Ki67/IL-15 (MOI = 40), and incubated for 72 h. The untreated cell was as a control. The cells were fixed in 95% ethanol (4 °C) overnight. Then, the fixed cells were washed twice with precooled PBS and incubated in 2.5 mg/ml DNase-free RNase A (Gene-Protein Link) and PI (1×, Gene-Protein Link) for 30 min at 37 °C. The distribution of the cell cycle was detected using a flow cytometry.
The cells (3000 cells/well) were seeded into a 96-well plate and cultured overnight. Then, the medium was replaced with 100 µl of different media and cultured for 1, 2, 3 days. At every time point, 10 µl of CCK-8 was added to each well and incubated for 2 h at 37 °C with 5% CO 2 . Then, the absorbance of each well was measured at 450 nm using a microplate reader (PerkinElmer, USA). At least three wells were used for each sample in different media.

Tube formation assay
Tube formation assays were performed according to previous descriptions [25]. Briefly, growth-factor-reduced Matrigel (BD, USA) was pre-added to 96-well plates. Cells (2 × 10 4 cells/well) were seeded into wells and incubated at 37 °C and 5% CO 2 . After 6 h, the cells were labeled using Calcein AM (Tocris, USA), and tube formation was imaged with an Olympus microscope. Capillary-like tube formation of GA-MSCs was analyzed in three random fields of view per well using ImageJ software (NIH, USA). The tube segment lengths and number of tubes of GA-MSCs cultured with different conditions were also quantified.

Coculture assay
For the coculture experiments, glioma cells GL261, U251, U87 and BT-01 (1 × 10 5 cells/well) and lymphocytes (1 × 10 6 cells/well) were, respectively, placed into the bottom wells of transwell permeable 6-well plate supports (Corning, Corning, NY), then the GA-MSCs (1 × 10 5 cells/well) were plated in the apical chamber, and the virus Ad5-Ki67/IL-15 (4 × 10 6 VP) was added into the upper chamber, which was co-cultured for 72 h. Inserts had a pore size of 0.4 μm, permitting the free exchange of molecules but preventing cell migration or contact. The experiments included different groups: cultures of glioma cells with addition of GA-MSCs, cultures of glioma cells with GA-MSCs plus virus, glioma cells with single virus and only glioma cells. All cells were cultured in a supplemented MSCs basal medium. After 3 days in coculture, the cells were collected to perform flow cytometry analysis.

Flow cytometry
Flow cytometry analysis was performed using fluorochrome-conjugated antibodies. Cultured glioma cells were detached with 0.05% trypsin-EDTA (ATCC), and enzymatic action was stopped by adding 10% FBS in PBS (Gibco, USA). The cells were washed in PBS, and then, the pellets were re-suspended in fluorescent-activated cell sorting (FACS) buffer. The glioma cell suspensions were stained with APC-conjugated antibodies against human or mouse PD-L1 (1:20, Proteintech, USA), and the lymphocytes were, respectively, stained with mouse PD-1 (1:100, Proteintech, USA), CD3, CD4 and CD8 (1:100, BioLegend, USA). All cells were incubated in the dark at 4 °C for 20 min. Then, the cells were centrifuged and re-suspended in PBS. Data were acquired within 2 h after staining on a BD Accuri C6 Plus (BD Biosciences), and analyzed using FlowJo software (Tree Star, Ashland, OR, USA).

Animal experiments
All animal studies were approved by Experimental Animal Welfare and Ethics Committee of Beijing Tiantan Hospital affiliated with Capital Medical University and approved all the animal experiments. Male BALB/nu-c (6 weeks old; Beijing Vital River Laboratory Animal) were kept in the animal facilities at Beijing Neurosurgical Institute and maintained under specific-pathogenfree conditions. Inhalation of isoflurane was used to anesthetize the animals in all experiments. BALB/nu-c (n = 5/group) was subcutaneously inoculated with a PBS suspension containing 1 × 10 6 U87-Luc cells (100 ul) in the right underarm using 1 ml syringe. The GA-MSCs suspension (100 ul, 1 × 10 5 cells) or Ad5-Ki67/IL-15 (3 × 10 11 vp) plus GA-MSCs (1 × 10 5 cells) were injected into the tumor sites 14 days after U87-Luc cell implantation. Tumor volumes were calculated according to the following formula: width 2 × length × 0.5.

IVIS imaging
The IVIS spectrum was utilized to obtain luciferase images. In vivo images were obtained on day 0 (prior to virus injection) and on days 3 and 6 after virus and GA-MSCs injection. Each mouse was treated with 100 μl 15 mg/ml d-luciferin (PerkinElmer, Waltham, MA, USA) i.p. in PBS.

Western blotting
The animals were anesthetized and sacrificed after treatment for 10 days. The tumor specimens were harvested from the mice, then washed two times and cut into pieces. The total protein was extracted with RIPA lysis buffers were treated by ultrasound and supplemented with mammalian protease inhibitor (Sigma-Aldrich, St Louis, MO, USA). Protein samples (20 μg per lane) were loaded into gels for separation via sodium dodecyl sulfate polyacrylamide gel electrophoresis. After separation, the proteins were transferred to nitrocellulose membranes (Thermo, Waltham, MA, USA), which were then blocked for 60 min at room temperature in TBST (Tris-buffered saline, 0.1% Tween 20) containing 5% nonfat milk. Membranes were washed with TBST and probed with primary antibodies against PD-L1 (1:800, all from Santa Cruz Biotechnology, USA) at 4 °C overnight. Membranes were washed with TBST and incubated with secondary antibodies for 1 h at room temperature. Signals were detected with an ECL detection system.

Statistical analysis
Statistical analyses were performed using GraphPad Prism 6.0. Sample sizes for each experiment are indicated in the Results or the Material and Methods section. The specific statistical tests used were t test for single comparisons and analysis of variance followed by Tukey's test for multiple comparisons, and P < 0.05 was considered to indicate statistical significance. Numerical values are reported as mean ± standard deviation.

GA-MSCs contributed to mediate PD-L1 expression in glioma cells
GA-MSCs were isolated from human fresh glioma specimens and showed similar classical MSC characteristics with fibroblastic morphology in standard medium in vitro (Fig. 1A, B). Immunofluorescence demonstrated that the GA-MSCs expressed MSC markers CD44 and CD105 (Fig. 1C). To explore GA-MSC-mediated PD-L1 expression of glioma, mouse glioma cell GL261 were co-cultured with GA-MSCs for 72 h. Flow cytometry showed that PD-L1 expression in GL261 cells was significantly upregulated compared with the control (Fig. 2A 47.67 ± 1.818) were, respectively, cocultured with GA-MSCs for 72 h. We also found that GA-MSCs contributed to mediate PD-L1 upregulation in human glioma cells (Fig. 2C-H).

The effects of Ad5-Ki67/IL-15 on growth and angiogenesis of GA-MSCs
To address the effects of the virus on GA-MSCs, we performed infection and proliferation assays for the cells. We found that the Ad5-Ki67/GFP can infect GA-MSCs (Fig. 3A). Furthermore, the cells were cocultured with the virus Ad5-Ki67/IL-15 and detected the cell cycle distribution by flow cytometry. The proportion of cell arrest in S (23.46%) and G2 (6.06%) phase in Ad5-Ki67/ IL-15 was increased compared to the control (S: 17.97%, G2: 2.99%) (Fig. 3B, C). CCK-8 assay displayed that Ad5-Ki67/IL-15 has the potential to suppress their growth (Fig. 3D, E). Furthermore, the conditioned media from The conditioned medium of glioma cell-treated with the virus significantly suppressed the cell growth compared with the untreated. All data are presented as mean ± SD of three independent experiments. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 glioma cells treated with Ad5-Ki67/IL-15 significantly inhibited proliferation of GA-MSCs (Fig. 3F, G). Tube formation assays performed in vitro showed that tube formation capacity was significantly reduced at 6 h in GA-MSCs treated with the virus compared with control (Fig. 4A). The results were also quantified; the tube segment lengths (P < 0.01, 53.56 ± 4.023 vs. 17.83 ± 4.426) and number of tubes (P < 0.01, 57.5 ± 4.95 vs. 13.5 ± 2.121) of the GA-MSCs treated with Ad5-Ki67/ IL-15 were significantly shorter than those of control (Fig. 4B, C). Furthermore, GA-MSCs incubated with supernatants from Ad5-Ki67/IL15-treated U251 cells had significantly decreased angiogenic capacity at 6 h compared with those incubated with conditioned medium from U251 cells; the tube segment lengths (P < 0.05, 41.91 ± 2.291 vs. 24.55 ± 1.782) and number of tubes (P < 0.05, 45.5 ± 6.364 vs. 25 ± 1.414) of the GA-MSCs treated with Ad5-Ki67/IL-15-CM were significantly shorter than those of control ( Fig. 4D-F). This provides a new strategy that inhibited tumor angiogenesis to treat GBM using Ad5-Ki67/IL-15 targeting GA-MSCs in future.

Ad5-Ki67/IL-15 overcame GA-MSCs to enhance therapeutic efficacy and reduce PD-L1 expression in glioma cells
Prior studies demonstrated that GA-MSCs promoted glioma growth and progression [8,10]. To evaluate the effects of Ad5-Ki67/IL-15 on PD-L1 expression in glioma cells, different glioma cells GL261, U251, U87 and BT-01, respectively, cocultured with GA-MSCs with the virus for 24 h, 48 h and 72 h, we confirmed that the virus has the potential to overcome GA-MSCs to effectively kill glioma cells (Fig. 5A-D). In addition, we found Ad5-Ki67/ IL-15 could attenuate PD-L1 expression of mouse glioma cells induced by GA-MSCs in GBM treatment (Fig. 6A (Fig. 6C-H). In addition, we also revealed that PD-1 in mouse lymphocytes was increased compared with controls and the virus could reduce GA-MSC-induced PD-1 expression and T cell inhibition (Additional file 2: Fig. S1 and Additional file 3: Fig. S2). Therefore, although GA-MSCs promote PD-L1 upregulation in glioma, the virus can attenuate this property of GA-MSCs.

Ad5-Ki67/IL-15 induced antitumor efficacy despite the presence of GA-MSCs in GBM model
To investigate the therapeutic efficacy of the virus in vivo, we established a human GBM model. IVIS images of mice bearing U87-derived tumors showed that tumor size was significantly decreased in the Ad5-Ki67/IL-15 plus GA-MSCs group compared with the GA-MSCs alone, especially 6 days after the virus injection (Fig. 7A). Intracranial tumors of glioma-bearing mice were collected; we found significant reductions in tumor volumes after the addition of Ad5-Ki67/IL-15 treatment (Fig. 7B,  C). Moreover, western blotting revealed that PD-L1 expression in tumor was downregulated by Ad5-Ki67/ IL-15 (Fig. 7D). These results suggested that oncolytic Ad5-Ki67/IL-15 can induce significant antitumor even if the presence of GA-MSCs in glioma.

Discussion
GBM is highly immunosuppressive and vascular tumor, resulting in poor immunotherapy outcomes. In this study, we for the first time discover that GA-MSCs contribute to the formation of PD-L1 upregulation in glioma. Further studies reveal that oncolytic Ad5-Ki67/IL-15 effectively reduces neovascularization in GA-MSCs, indicating potential anti-tumor angiogenic capacity. Furthermore, the virus decreased PD-L1 expression of glioma and induced antitumor efficacy by targeting GA-MSCs.
To our knowledge, this is the first description of GA-MSC-mediated PD-L1 expression in brain glioma.
GA-MSCs, as resident stromal components in the unique microenvironment of human brain tumors, severely affect prognosis and survival time of GBM patients [9]. The fraction of MSCs in glioma tissues may be distinct in different patients, even within the same histopathological characteristics [9]. This signifies that therapeutic strategy targeting MSCs may vary from patient to patient. Therefore, deeper research about GA-MSCs in glioma will provide a reliable perspective for the development of individualized treatment in future. Previous studies about GA-MSCs mainly focused on the effects of tumor initiation, progression, driving GSCs proliferation and therapeutic resistance [8,26,27]. However, there are few studies on its immune regulation in glioma. An important finding of this study was that GA-MSCs can upregulate PD-L1 expression and can be downregulated by the virus. Moreover, the virus was still able to achieve the desired antitumor effect in a human GBM model containing GA-MSCs. The source of GA-MSCs is unknown, and the most probable explanation for their presence is that the tumor recruits them from normal tissues [26,28], as well as GSCs differentiation [26]. The involvement of MSCs in angiogenic activity has been well established [10,29]. Furthermore, GA-MSCs have the potential to differentiate into pericytes and maintain tumor vascular structure, and the GA-MSCs with CD90 low expression show a more active angiogenic capacity, which contribute to glioma progression [10,30]. Studies on GA-MSCs have mainly focused on their effects on tumor progression [8,26], whereas their immunological properties and treatments targeting GA-MSCs have not been reported. Our study found that GA-MSCs showed potent immunosuppressive characteristics, the increased numbers of both PD-1 + and PD-L1 + cells were found.
Recently, the effects of MSCs on T cells have been intensively studied. Immunosuppressive factors released by MSCs, such as TGF-β, IL-10 and HGF, which are thought to be main causes of the inhibition of T cell proliferation [14]. These factors lead to cyclin D2 downregulation and p27kip1 upregulation in T cells, thereby arresting proliferation in the G1 phase, and have been shown to induce immunosuppression within tumors [15,18,19]. MSCs also generate a variety of chemokines and adhesion molecules, which play important parts in lymphocyte recruitment, thus ensuring their immunosuppressive function. Increased expression of immunoregulatory and adhesion molecules is essential for effective T cell inhibition [17,31,32]. In addition, MSCs have the capacity to alter the activation and differentiation of T cells. There is evidence that MSCs decrease IFN-γ and IL-17 secretion by T cells but promote IL-10 production by antagonizing the differentiation of Th1 and Th17 cells, thereby inducing the generation of regulatory T cells [11,13,33]. Furthermore, our study demonstrated that GA-MSCs can upregulate PD-L1 expression in glioma, indicating a potential immunosuppressive property.
PD-L1 is expressed in various tissues and is often upregulated in malignant tumors; higher PD-L1 expression in GBM has been found to be correlated with poorer patient prognosis in some studies [34,35]. PD-L1 is an immune inhibitory receptor ligand that leads to immune cell dysfunction and apoptosis by binding to its receptor PD-1 [36]. This facilitates the immunosuppressive microenvironment and tumor progression. An increasing number of studies hint that PD-1/PD-L1 is a promising target to reverse the immune evasion of GBM [37][38][39]. In this All data are presented as mean ± SD of three independent experiments. *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 study, we demonstrate that GA-MSCs have the potential to increase PD-L1 of glioma and PD-1 in lymphocytes, indicating a potential immunosuppressive capacity. PD-L1 hinder T cell activation, especially during the effector phase, leads to T cell exhaustion during chronic antigen exposure, and inhibits NK cells and B cells; however, antagonists of PD-L1 or its ligands have been shown to reverse immune function at least in part and to enhance OV therapy [40,41]. Belcaid et al. found that Delta24-RGD oncolytic virus therapy overcome gliomainduced immune suppression and increase PD-1 + tumorinfiltrating CD8 + T cells [42]. Our investigation revealed that oncolytic Ad5-Ki67/IL-15 downregulate PD-L1 expression in glioma with the presence of GA-MSCs.
There are some limitations in our study. For example, the GA-MSCs that we used in the experiment were derived from patients diagnosed with glioma, but the fraction of the cells is not abundant in some glioma tissues, and that the heterogeneity of MSC will present great challenges for this research. Our data demonstrate that OVs improve GBM treatment by remodeling GA-MSCs-mediated PD-L1 expression and angiogenesis. However, we do not determine the interaction of GA-MSCs-mediated PD-L1 and immune cells and also not evaluated the specific manner by which OVs affect angiogenic capacity in this study.

Conclusions
In summary, our results confirmed that GA-MSCs participate in the upregulation of PD-L1 expression of glioma. Our findings further indicate that GA-MSCmediated PD-L1 can be attenuated by Ad5-Ki67/IL-15 in treatment, thereby enhancing the effectiveness of immunotherapy in GBM. This study indicates a potential new approach for GBM treatment using Ad5-Ki67/IL-15 to target not only tumor cells but also MSCs.