The pesticide chlorpyrifos promotes obesity by inhibiting diet-induced thermogenesis in brown adipose tissue

Wang, Bo; Tsakiridis, Evangelia E.; Zhang, Shuman; Llanos, Andrea; Desjardins, Eric M.; Yabut, Julian M.; Green, Alexander E.; Day, Emily A.; Smith, Brennan K.; Lally, James S. V.; Wu, Jianhan; Raphenya, Amogelang R.; Srinivasan, Krishna A.; McArthur, Andrew G.; Kajimura, Shingo; Patel, Jagdish Suresh; Wade, Michael G.; Morrison, Katherine M.; Holloway, Alison C.; Steinberg, Gregory R.

doi:10.1038/s41467-021-25384-y

Download PDF

Article
Open access
Published: 27 August 2021

The pesticide chlorpyrifos promotes obesity by inhibiting diet-induced thermogenesis in brown adipose tissue

Nature Communications volume 12, Article number: 5163 (2021) Cite this article

20k Accesses
43 Citations
498 Altmetric
Metrics details

Subjects

Obesity

Abstract

Obesity results from a caloric imbalance between energy intake, absorption and expenditure. In both rodents and humans, diet-induced thermogenesis contributes to energy expenditure and involves the activation of brown adipose tissue (BAT). We hypothesize that environmental toxicants commonly used as food additives or pesticides might reduce BAT thermogenesis through suppression of uncoupling protein 1 (UCP1) and this may contribute to the development of obesity. Using a step-wise screening approach, we discover that the organophosphate insecticide chlorpyrifos suppresses UCP1 and mitochondrial respiration in BAT at concentrations as low as 1 pM. In mice housed at thermoneutrality and fed a high-fat diet, chlorpyrifos impairs BAT mitochondrial function and diet-induced thermogenesis, promoting greater obesity, non-alcoholic fatty liver disease (NAFLD) and insulin resistance. This is associated with reductions in cAMP; activation of p38MAPK and AMPK; protein kinases critical for maintaining UCP1 and mitophagy, respectively in BAT. These data indicate that the commonly used pesticide chlorpyrifos, suppresses diet-induced thermogenesis and the activation of BAT, suggesting its use may contribute to the obesity epidemic.

Arsenite exposure suppresses adipogenesis, mitochondrial biogenesis and thermogenesis via autophagy inhibition in brown adipose tissue

Article Open access 08 October 2019

Aurantio-obtusin ameliorates obesity by activating PPARα-dependent mitochondrial thermogenesis in brown adipose tissues

Article 24 April 2023

Obesogens: a unifying theory for the global rise in obesity

Article Open access 11 January 2024

Introduction

Obesity is a major risk factor for type 2 diabetes (T2D), non-alcoholic fatty liver disease (NAFLD), and cardiovascular disease¹ that arises from a caloric surplus of as little as 10–30 kcal per day². And while increased consumption of energy-dense foods and reduced physical activity are commonly thought to be the major contributors to this caloric imbalance¹, diet-induced thermogenesis is a quantitatively important component of the energy balance equation³. In adult humans, recent studies have indicated that diet-induced thermogenesis requires the activation of brown adipose tissue (BAT)^4,5,6, however, the determinants regulating this process and why they may differ between individuals are not fully understood.x

A key protein regulating BAT thermogenesis is uncoupling protein 1 (UCP1). Genetic removal of Ucp1 in mice promotes obesity and insulin resistance when mice are fed a diet high in fat and housed under thermoneutral conditions (29–30 °C), indicating a vital role for this protein in diet-induced thermogenesis⁷. UCP1 is also important for diet-induced thermogenesis in humans, as a genetic loss of function mutation reduces postprandial thermogenesis in response to a high-fat meal⁸. While genetic polymorphisms in Ucp1 have been linked to obesity in some studies^9,10,11, this has not been observed in large genome-wide association studies¹², suggesting that other factors besides genomic alterations in Ucp1 may be contributing to reduced BAT activity in obese humans¹³.

Over the last two decades, a number of environmental toxicants have been linked to the development of obesity through their effects on gut microbiota¹⁴, energy intake, or adipogenesis^15,16,17,18 (for review see refs. ^19,20), however, only a few of these studies have examined a role for these agents to contribute to obesity by inhibiting energy expenditure^21,22,23 and/or BAT thermogenesis^24,25 (reviewed in refs. ^18,20). We hypothesized that environmental toxicants commonly present in food might reduce diet-induced thermogenesis through suppression of UCP1, the defining protein of human BAT thermogenesis^{26,27,28,29,30,31}. Through a step-wise screening of several common food contaminants or additives, we discovered that the organophosphate insecticide chlorpyrifos (CPF) potently suppressed the expression of UCP1 and mitochondrial respiration in brown adipocytes at concentrations as low as 1 pM. Chloropyrifos-induced suppression of diet-induced thermogenesis in BAT has also been observed in mice fed a diet high in fat and housed at thermoneutrality where it promoted greater obesity, NAFLD, and insulin resistance. Reductions in BAT thermogenesis by CPF were associated with reductions in cAMP and protein kinases critical for regulating UCP1 and mitophagy. These data indicate that the commonly used pesticide CPF, at very low concentrations, suppresses the activation of BAT, suggesting that its use may contribute to the obesity epidemic.

Results

We examined the effects of 34 chemicals—selected based on their widespread presence in food due to agriculture practices, food processing, and/or packaging (Table 1) on the expression of Ucp1 in immortalized brown adipocytes generated from Ucp1-luciferase reporter mice³². Of the 34 chemicals tested, only the pesticide CPF significantly decreased Ucp1 promoter activity and mRNA expression by at least 25% across six orders of magnitude of concentrations tested down to 1 pM (Fig. 1a, b, Supplementary Fig. 1a–c).

Table 1 List of tested chemicals.

Full size table

**Fig. 1: CPF inhibits the thermogenic gene program and mitochondrial respiration in cultured brown adipocytes.**

CPF is an organophosphate insecticide with broad-spectrum activity against many foliar and soil insects and is commonly used for pest control in a wide range of field crops and fruit³³. Previous studies in animals have found that high doses of CPF, which inhibits plasma butyryl cholinesterase and brain acetylcholinesterase activity, can promote obesity and glucose dysregulation through mechanisms that are poorly defined but may involve increases in energy intake and/or alterations in the gut microbiome^{34,35,36,37,38,39}. However, whether doses of CPF below those that inhibit acetylcholinesterase and closer to human-relevant exposures promote obesity in rodents or humans is currently unknown. We found that consistent with the acute suppression of Ucp1 promoter activity and Ucp1 mRNA, chronic CPF (6 days) also suppressed UCP1 protein in brown adipocytes at a concentration of 1 pM (Fig. 1c and Supplementary Fig. 1c). To further characterize the effects of CPF, we completed global RNA-sequencing after just 4 h of CPF treatment and found that RNA transcripts related to mitochondrial function/metabolism were most affected (Fig. 1d). Notably, in addition to Ucp1, RNA transcripts are critical for regulating fatty acid oxidation (Cpt1a, Cpt1b, Acat3), and Cytochrome C oxidase assembly factor (Cox16) (Fig. 1e) were reduced with CPF treatment. Consistent with previous studies in UCP1 null mice^40,41, compensatory upregulation of some mitochondrial transcripts was also observed (Fig. 1e). Importantly, consistent with alterations in transcriptional programs regulating mitochondrial function, chronic treatment with CPF-lowered respiration in the presence of the mitochondrial ATP synthase inhibitor oligomycin and the mitochondrial uncoupler FCCP; indicating impaired mitochondrial leak and maximal respiration, respectively (Fig. 1f). CPF also lowered cytochrome c oxidase activity (Fig. 1g) and mitochondrial membrane potential (Fig. 1h). These reductions in mitochondrial function were not due to impairments in BAT differentiation as Oil Red O staining and key transcription factors/markers of BAT differentiation (Prdm16, Ppargc1a, Pparg, Cidea) were not altered by CPF (Supplementary Fig. 1d, e). Impaired mitochondrial function was not observed following acute CPF exposure (4 h) (Supplementary Fig. 1f, g) as had previously been observed at high concentrations of CPF (40 µM)⁴². These data indicate that concentrations of CPF significantly below those causing toxicity in animal studies and more aligned with real-world human exposures⁴³ inhibit transcriptional programs crucial for regulating BAT mitochondrial function.

CPF is readily detected in a wide range of foods including fruits, vegetables, grains, beans, nuts, legumes, dairy, meat, fish, and eggs, at concentrations ranging from 0.002 to 436 μg/kg⁴⁴. We examined the effects of CPF in male C57BL6J mice that were fed a control diet (CD, 10 kcal% fat) or high-fat diet (HFD, 45 kcal% fat) supplemented with two doses of CPF (5 or 20 mg/kg diet, which equated to an exposure of ~0.5 or 2 mg/kg body weight in mice fed a control diet housed at room temperature (CD-RT) with slightly lower exposure in HFD fed mice housed at thermoneutrality (TN, Supplementary Fig. 2a). In contrast to an acute gavage of CPF (100 mg/kg) which resulted in dramatic inhibition of ChE activity in the serum, we found that CPF in the diet did not inhibit acetyl-cholinesterase activity in skeletal muscle (Supplementary Fig. 2b and c), which is consistent to a previous study³⁴.

Mice housed at room temperature (RT, 21–23 °C), common in most animal vivarium, are under cold stress and have increased sympathetic tone that upregulates basal metabolic rate through adaptive thermogenesis^7,45. This elevation in basal metabolic rate masks the effects of diet-induced thermogenesis mediated through β-adrenergic driven activation of BAT UCP1 compared to housing mice at thermoneutrality (TN, 29–30 °C) as evidenced by the development of obesity in mice lacking UCP1 at TN but not RT. Thus TN housing is a temperature range for mice that better mimics a condition in which humans reside⁴⁵ and as such recent studies have indicated that it more accurately models the development of metabolic diseases including obesity, insulin resistance, and NAFLD^46,47. Therefore, mice were housed at either RT or TN. To the best of our knowledge, this is the first time an environmental toxicant has been tested in rodents housed at TN.

On a control diet, CPF supplementation did not influence body weight (Supplementary Fig. 3a, b), fat mass (Supplementary Fig. 3c, d), glucose tolerance (Supplementary Fig. 3e, f), or insulin sensitivity (Supplementary Fig. 3g, h) in RT or TN conditions. However, when mice were housed at TN and fed an HFD, CPF treatment increased body mass (Fig. 2a), adiposity (Fig. 2b), and reduced glucose tolerance (Fig. 2c), and insulin sensitivity (Fig. 2d) compared to non-CPF treated mice. Importantly, these effects were blunted in mice housed at RT such that only the highest dose of CPF (2.0 mg/kg) promoted weight gain, glucose intolerance, and insulin resistance (Supplementary Fig. 4a–e). Consistent with increased adiposity under TN-HFD-fed conditions (Fig. 2b, e), CPF increased the size of WAT lipid droplets (Fig. 2f) and elevated serum free fatty acids (FFA, Fig. 2g) and triglycerides (TGs, Fig. 2h). Consistent with increased obesity, insulin resistance and FFAs HFD-fed mice treated with CPF and housed at TN had increased liver weights (Fig. 2i), liver TGs (Fig. 2j, k) and serum ALT (Fig. 2l) and AST (Fig. 2m) indicative of greater NAFLD than HFD controls. Importantly, these effects of CPF were blunted when fed an HFD at RT and were not observed when mice were fed with CD indicating these markers of NAFLD/liver toxicity were likely secondary to the increased weight rather than toxicity to CPF (Supplementary Figs. 5 and 6, respectively). These data indicate that low concentrations of CPF, which do not inhibit AchE activity, when combined with TN housing and an HFD promote obesity, insulin resistance, and NAFLD.

**Fig. 2: CPF promotes obesity and metabolic dysfunction in mice fed a high-fat diet and at thermoneutrality.**

To investigate the mechanisms by which CPF promoted HFD-induced weight gain, caloric intake and energy expenditure were assessed in metabolic cages. Consistent with previous studies mice housed at TN had lower caloric intake compared to mice housed at RT⁴⁸ and regardless of housing temperature, CPF increased caloric intake^36,49 (Fig. 3a). No differences in physical activity were observed between groups (Fig. 3b). As anticipated, TN lowered absolute (not corrected for body mass) resting oxygen consumption compared to RT housing (Fig. 3c). Importantly, CPF suppressed oxygen consumption by ~15% in mice housed at TN but not RT (Fig. 3c). These data indicate that CPF suppresses energy expenditure at TN but not RT. To further examine whether reductions in energy expenditure were important for the weight gain observed in HFD-fed mice housed at TN, in separate experiments mice were fed with CPF and allowed ad libitum access to food or were pair-fed to controls (Fig. 3d). Importantly, even in the pair-fed group CPF treated mice had a strong tendency for higher body weight and increased adiposity after 9 weeks of HFD feeding when the amount of food intake was matched with that of the control mice (Fig. 3e and f). Collectively these data suggest that CPF reduces energy expenditure, contributing to the increased weight gain observed at TN.

**Fig. 3: Dietary exposure to chlorpyrifos promotes weight gain in HFD-fed mice housed at thermoneutrality even in the absence of increased food intake.**

Our data demonstrating reductions in oxygen consumption at TN, but not at RT with CPF, are similar to observations in UCP1 null mice, which also only develop HFD-induced obesity when housed at TN; an effect attributed to impairments in BAT diet-induced thermogenesis⁷, We found that in HFD-fed mice housed at TN, CPF reduced Ucp1 mRNA and protein expression (Fig. 4a and b) and increased the size and number of lipid droplets (Fig. 4c) within the BAT. In addition to suppressing UCP1, transmission electron microscopy revealed an increased abundance of large mitochondria (Fig. 4d and e) with disrupted cristae (Fig. 4d and f). These effects of CPF on UCP1 and BAT morphology and mitochondria were not evident when mice were housed at RT (Supplementary Fig. 7). These data indicate that consistent with increased weight gain and metabolic dysfunction when mice are housed at TN, CPF reduces BAT UCP1 expression and mitochondrial morphology.

**Fig. 4: Dietary exposure to chlorpyrifos impairs mitochondrial function and thermogenesis in brown adipose tissue in vivo.**

To directly test the hypothesis that CPF was suppressing diet-induced thermogenesis we conducted a subsequent study in which mice were housed at TN and fed an HFD with or without CPF for 5 weeks (Fig. 5a), a time point preceding any differences in body mass or adiposity (Supplementary Fig. 8a, b). Mice were then placed in metabolic cages and fasted overnight before refeeding. CPF-treated mice had similar oxygen consumption (Fig. 5b) and energy expenditure (Fig. 5c) during fasting. However, CPF-treated mice had an ~10% reduction in oxygen consumption and energy expenditure during refeeding (Fig. 5b and c) indicative of impaired diet-induced thermogenesis. Importantly, this was not due to differences in food intake during the refeed (Fig. 5d) or activity levels (Fig. 5e). These data suggest that CPF reduces energy expenditure by suppressing diet-induced thermogenesis.

**Fig. 5: Chlorpyrifos inhibits diet-induced thermogenesis.**

In mice diet-induced thermogenesis is mediated through norepinephrine stimulated β₃-adrenergic receptor (β₃-AR) activation, which results in activation of adenylyl cyclase, and subsequent increases in cAMP and PKA activity^50,51. To examine whether this pathway may be altered with CPF, we measured cAMP in response to the pan beta-agonist isoproterenol and found that CPF reduced cAMP in isolated brown adipocytes (Fig. 6a). To examine potential mechanisms which might be contributing to this reduction in cAMP several experiments were completed. We first hypothesized that CPF may directly inhibit the β₃-AR. We built a homology model of the β₃-AR which was used as a starting structure to carry out classical molecular dynamic simulations (MD). Snapshots obtained from MD simulation were then used to dock CPF. The most populated docked pose shows that CPF is predicted to bind to the known pocket in the β₃-AR and predominantly forms hydrophobic contacts and hydrogen bonds with S200. (Supplementary Fig. 9a) Surprisingly, the docking score of the most populated CPF docking pose was higher (−7.4 kcal/mol, lower the better) than cyanopindolol (−8.8 kcal/mol), a potent β-adrenergic receptor antagonist, suggesting a weaker binding. To validate the modeling observations, studies using a radioligand competition assay of CPF or its metabolite CPF-oxon with cyanopindolol were conducted. However, in contrast to the docking studies, we found that neither CPF nor CPF-oxon competed with cyanopindolol for β₃-AR binding (Supplemental Fig. 9b). These data suggest that neither CPF nor its metabolite CPF-oxon directly bind the β₃-AR.

**Fig. 6: Chlorpyrifos inhibits diet-induced thermogenesis and β-adrenergic signaling in brown adipose tissue.**

We subsequently interrogated the hypothesis that CPF could be reducing cAMP by altering the expression of β-ARs (Adrb), adenylyl cyclase (Adcy), or phosphodiesterases (Pde) but found no change in isolated BAT cells (Supplementary Fig. 9c–e). In addition to direct binding or alterations in expression, cAMP production is dependent on membrane localization of adenylate cyclase therefore we subsequently examined whether this might be reduced by CPF treatment. Interestingly, CPF lead to a reduction of membrane-bound ADCY3 with a corresponding increase in cytoplasmic localization (Supplementary Fig. 9f). These data suggest that CPF may reduce cAMP by reducing membrane localization, however, the exact mechanisms that contribute to this effect require further investigation.

To examine the downstream consequences of this reduction in cAMP we measured the phosphorylation and activation of hormone sensitive lipase (HSL), p38MAPK, and AMP-activated protein kinase (AMPK), which are important for increasing lipolysis, UCP1 expression, and mitophagy respectively^50,51,52,53. Reductions in cAMP in isolated brown adipocytes (Fig. 6a) were associated with reduced activating phosphorylation of HSL (Ser660) and p38MAPK (Fig. 6b); effects similar to the pan-beta antagonist propranolol (Supplementary Fig. 10a). Reductions in the phosphorylation of PKA, HSL, and p38 MAPK were also observed in BAT of CPF-treated mice fed an HFD and housed at TN (Supplementary Fig. 10b–d). Concomitant with these reductions in cAMP the expression of Zfp516, a cAMP-responsive transcription factor regulating UCP1⁵⁴, was also reduced (Fig. 6c).

Treatment of mouse and human BAT cells with norepinephrine increases cAMP and AMPK activity⁵⁵. Similar to CPF-treated mice, mice lacking AMPK in BAT have an increased abundance of large mitochondria with disrupted cristae, an effect which is mediated through AMPK-induced increases in mitophagy⁵⁵. We found that consistent with mice genetically lacking AMPK in BAT, CPF-treated mice and BAT cells have reduced activating phosphorylation of AMPK and ULK1 and reduced accumulation of LC3B-II and increased p62, the latter being two markers of reduced autophagic flux (Fig. 6d and Supplementary Fig. 10e). There was also reduced phosphorylation of ACC, another downstream substrate of AMPK important for regulating fatty acid oxidation and fatty acid synthesis (Fig. 6e and Supplementary Fig. 10e). To test whether this inhibition of AMPK was important for mediating the inhibitory effects of CPF on mitochondrial function, brown adipocytes were treated with the highly selective AMPK agonist MK-8722⁵⁶ which eliminated the ability of CPF to suppress mitochondrial uncoupling by FCCP (Fig. 6f). These data indicate that CPF reduces BAT cAMP and downstream signaling leading to reductions in UCP1 and pathways critical for controlling mitochondrial function.

Discussion

A small imbalance between energy expenditure and energy intake/absorption increases adiposity and over time can lead to the development of obesity. Due to the lipophilicity of environmental toxicants, adipose tissue is the major site of environmental toxicant accumulation^57,58. In agreement with previous studies^35,36,37,38, we found that treatment of HFD-fed mice with 2.0 mg/kg body weight of CPF increased caloric intake and promoted obesity, insulin resistance, and NAFLD when mice were housed at RT. Interestingly no effect of this dose was observed on metabolic parameters when this dose of CPF was combined with the control low-fat diet indicating important interactions between CPF and dietary fat. Importantly, this dose of CPF is below that associated with the inhibition of acetyl-cholinesterase activity or toxicity and is more aligned with real-world exposures in humans^39,59.

When mice are housed at TN, but not RT, a low concentration of CPF (0.5 mg/kg body weight) promotes weight gain, insulin resistance, and NAFLD. To the best of our knowledge, this is the first toxicological study conducted in rodents that have been completed under TN conditions. The findings demonstrating effects of low concentrations of CPF at TN, but not RT, are consistent with recent studies that have highlighted that RT does not allow for accurate modeling of human disease due to increased thermogenic demands that mask biological processes that occur in humans who predominately reside at TN⁶⁰. A key reason for this is because at RT mice require adaptive thermogenesis to maintain their body temperature which in turn masks diet-induced thermogenesis. Therefore, the observation of increased weight gain at TN but not RT, despite increases in caloric consumption under both housing conditions, suggested that low dose CPF may be inhibiting diet thermogenesis and this could be critical for mediating the increased weight gain. Consistent with this hypothesis pair-feeding experiments indicated that CPF mice housed at TN who were calorically matched to the lower food intake of controls still gained more weight suggesting reductions in energy expenditure were important. To further examine whether diet-induced thermogenesis was critical for the weight gain mice were housed at TN and fed an HFD containing CPF for 4 weeks, and before any differences in adiposity emerged, were fasted and refed in metabolic cages. These experiments demonstrated that despite similar activity levels and food consumption, oxygen consumption was reduced in CPF-treated mice during refeeding. These findings indicating that CPF suppresses diet-induced thermogenesis indicate a mechanism by which CPF promotes obesity and metabolic dysfunction. They also suggest that since to date all environmental toxicants have been tested in rodents housed at RT, the effects of these chemicals on human health may have been underestimated and therefore additional toxicological testing at TN may be important.

BAT and UCP1 are important for diet-induced thermogenesis in both mice and humans⁵⁰. Increased weight gain and metabolic dysfunction in mice treated with low doses of CPF housed at TN were associated with a reduction in UCP1 and alterations in BAT mitochondria morphology. Interestingly, this effect was not observed when mice are housed at RT potentially because of the increased adrenergic drive which may have overridden the suppressive effects of CPF. These findings are in agreement with our observations in vitro that low pM concentrations of CPF suppress cAMP and downstream signaling pathways in BAT cells, effects associated with reductions in UCP1 expression, mitochondrial function, and oxygen consumption. In rodents, diet-induced thermogenesis is dependent on norepinephrine-induced activation of the β3-AR and UCP1 and in the current study, CPF exposure led to reductions in both cAMP and UCP1. We hypothesized that this may be due to binding of CPF to the β3-AR or reductions in β3-AR expression levels, however, we found no effect of CPF on either of these parameters. Previous studies in developing neurons found that CPF may reduce the activity of adenyl cyclase⁶¹ and consistent with this observation we observed a reduction in membrane localization. The exact mechanisms contributing to this effect are unknown however one possibility is that given that CPF and its active metabolite CPF oxon has been shown to disrupt microtubule networks this may be affecting intracellular localization⁶². Additional studies are needed to further investigate the mechanisms by which CPF reduces adrenergic signaling in BAT.

One of the limitations of our study is that we only examined the effects of CPF in male mice. It is known that CPF interacts with sex hormones⁶³ and steroid hormone receptors⁶⁴. In terms of the neurotoxic effects, males are more susceptible to organophosphates than females⁶⁵. However, as reviewed⁶⁶, women may be more vulnerable to environmental pollution than men mainly because of the physiological differences caused by sex steroids. With a larger relative fat mass, women have a larger distribution volume for lipophilic substances than men⁶⁶. Indeed, in aged human beings, concentrations of organochlorine pesticides are elevated in females⁶⁷. Females may also have higher levels of thermogenic BAT than males⁶⁸. Upon adrenergic activation, women have more potential to induce browning of perirenal adipose tissue than men due to sex-specific intrinsic factors⁶⁹ and have higher lipolytic activities⁷⁰. Interestingly, females exhibit higher DIT than males at 20 °C but lower DIT than males at TN⁷¹. Future studies examining the effects of CPF in female mice are needed.

In conclusion, our study modeled the effects of an environmental toxicant under TN housing conditions. These studies have revealed that in addition to promoting food intake, CPF suppresses diet-induced thermogenesis in BAT and this exacerbates the development of obesity, NAFLD, and insulin resistance even at low doses which have no effect in mice housed at RT. These studies suggest that the effects of environmental toxicants on the development of obesity may have been underestimated as all studies to date have been conducted in mice housed at RT. Future studies examining the mechanisms driving reductions in β-AR signaling and whether there are associations between BAT metabolic activity and CPF in humans will be important.

Methods

In-vitro analysis

Immortalized brown adipocyte cell line generation

Cell lines were created from either FVB/N mice (wildtype immortalized brown adipocytes) or UCP-1-Luc2-TdTomato reporter mouse (immortalized UCP-1 reporter brown adipocytes), also known as the “ThermoMouse”³² using reported techniques⁷². Briefly, BAT tissue was harvested from 4-day-old pups and subsequently digested with collagenase II, filtered, and plated. Isolated cells were then infected with a retrovirus created by transfecting Phoenix-ECO cells with a pBABE-SV40 plasmid from Addgene. Cells positive for SV40 were selected for 7 days using 2 µg/ml puromycin antibiotic and underwent a 14-day treatment with the antibiotic ciprofloxacin (10 µg/ml) in order to eliminate any mycoplasma contamination in the cells. UCP-1 Luc2-TdTomato (UCP-1 reporter line) was used for all luciferase assays to quantify UCP-1 promoter activity. For the cyclic AMP assay wildtype immortalized brown adipocytes were transfected with Promega pGloSensor-22F which can detect dynamic changes in cyclic AMP from concentration 0.003–10 µM in living cells.

Immortalized murine brown adipocyte culture

Immortalized brown adipocytes were cultured in high glucose DMEM (Gibco^TM) containing 10% fetal bovine serum and 1% antibiotic–antimycotic solution. For differentiation, confluent cells were placed in induction media (0.5 mM IBMX; 0.125 mM indomethacin; 0.2 μg ml⁻¹ dexamethasone; 1 nM T3, 20 nM insulin) for 2 days and subsequently maintained on differentiation media (1 nM T3, 20 nM insulin) for 7 days. For all cell experiments, cells were fully differentiated prior to compound or vehicle (≤0.1% dimethyl sulfoxide; DMSO) treatment.

Cyclic AMP assay

Cells were grown to confluence and fully differentiation as described above. Fully differentiated cells were seeded into 96-well assay plates at a 1:1 dilution based on well area. Following 48 h post-seeding, cells were treated with vehicle, CPF, or positive control propranolol for 7 days. Fresh media containing compound concentrations was added every 48 h. On day 7 a kinetic luciferase assay was performed. A 30 min basal read was first obtained, followed by a 100 nM isoproterenol injection and a 30 min stimulated kinetic read. Cyclic AMP activity is expressed as luminescence (RLU). Data are expressed as fold basal change (stimulated read (5–20 min)/basal read (5–30 min)).

MitoTracker staining

After washing with PBS, cells were treated with 200 nM MitoTracker Red and 100 nM MitoTracker Green staining solutions in DMEM (serum-free) for 25 min at 37 °C. Cells were then washed with PBS and stained with 1 µg/ml DAPI in DMEM. Images were taken using fluorescence microscopy and quantified using Image J.

Cytochrome c oxidase activity

Wildtype immortalized murine brown adipocytes were differentiated in 12-well plates. Cytochrome c oxidase activity was determined calorimetrically. The reaction mixture containing potassium phosphate buffer, reduced cytochrome c, and ddH₂O, and cell lysate. Absorbance was immediately read at 550 nm wavelength for 3 min using SpectraMax M5 plate reader. The reaction was performed at 37 °C and cytochrome c oxidase activity was calculated as Eq. (1):

$${{\mathrm {Activity}}}=\frac{\triangle {{\mathrm {OD}}}/{{\mathrm {reaction}}}\,{{\mathrm {time}}}}{{{\mathrm {protein}}}\,{{\mathrm {concentration}}}\times {{\mathrm {sample}}}\,{{\mathrm {volume}}}}\times \frac{{{\mathrm {reaction}}}\,{{\mathrm {volume}}}}{{{\mathrm {exctinction}}}\,{{\mathrm {coefficient}}}\times {{\mathrm {path}}}\,{{\mathrm {length}}}}$$

(1)

An extinction coefficient of 18.5 mmol⁻¹ cm⁻¹ (reduced cytochrome c), a path length of 0.625 cm, and a reaction volume of 0.0002 L were used. Protein was quantified during final calculations to determine the activity of enzyme per µg of protein. Each sample was run in triplicate on a 96-well plate.

In vitro respiration assay of brown adipocytes

The immortalized brown preadipocytes cell line were fully differentiated in a 75 cm² flask, then transferred to a 96-well flat-bottomed OxoPlate (OP96C, PreSens Precision Sensing) at a dilution of 1:8 by culture area. Before the assay, cell culture media was replaced with 200 μl respiration medium (minimal DMEM-D5030 supplemented with 10 mM HEPES, 25 mM glucose, 31 mM NaCl, 2 mM sodium pyruvate, 2 mL l-glutamine, and 4% BSA, pH = 7.4), covered by 75 μl heavy mineral oil and prewarmed at 37 °C in a SpectraMax M5 plate reader (Molecular Devices, San Jose, CA, USA) for 45 min. For calibration, each plate contained six wells of 0% O₂ standard (H₂O with 10 mg/mL sodium sulfite) or 100% O₂ standard (respiration media). The plate was kinetically read prior (basal) and following the addition of 2.5 μM oligomycin, 1 μM FCCP or 2 μM Rotenone +2 μM Antimycin A injections for 20 min. Oxygen levels were calculated according to the manufacturer’s manual.

RNA-sequencing analysis

TMCON brown pre-adipocytes were differentiated in 24-well plates and treated with 1 pM CPF on day 9 of differentiation. Cells were treated for 4 h in Krebs–Ringer bicarbonate buffer (Sigma-Aldrich) supplemented with 10 mM HEPES (Janssen Pharmaceuticals, 3a2440 Geel, Belgium), 1 mM CaCl₂, and 1% fatty-acid free Bovine Serum Albumin (BSA) (EquiTech Bio Inc., TX, USA) at pH 7.4. Cell lysates were collected by administering 175 µL per well of RNeasy Lysis Buffer (RLT) (Qiagen, Hilden, Germany) supplemented with 143 mM beta-mercaptoethanol (Sigma-Aldrich) and vortexing until cells lifted from plates. Once cells were lifted from plates, 175 µL of 70% ethanol was added to each well, and lysates were transferred to a High Pure Filter tube from Roche RNA Isolation Kit (Roche Diagnostics, Mannheim, Germany). RNA isolation was done as outlined by manufacturers' specifications.

Raw sequencing reads (~12.5 million paired-end per sample; 50 base pairs each) were assessed for quality (such as GC content, PHRED scores, synthetic aptamer content, and sequence length) using FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Trimmomatic⁷³ (default parameters) was used to trim low-quality reads, with the final set of reads aligned to the Ensembl Mus musculus (GRCm38.92) genome using HISAT2⁷⁴ (default parameters). The number of reads mapping to each gene was identified using htseq-count, with stranded set to “reverse”, minimum alignment quality set to 10, feature type set to “exon” and “simple” advanced options. Mouse GTF gene annotation reference files were obtained from Ensembl. DESeq2⁷⁵ (default parameters) was used to normalize these data and a custom python script (McArthur) was used to convert the DESeq2 output into a ranked list of genes. The ranked list was used as input for Gene Set Enrichment Analysis (GSEA)^76,77, which was conducted using Cytoscape⁷⁸. GSEA was performed using 1000 permutations, a false discovery rate cut-off of 25%, and a mouse gene set database from the Bader Lab at the University of Toronto (http://download.baderlab.org/EM_Genesets/current_release/Mouse/).

Oil red O staining and quantification

Cells were rinsed with PBS and fixed with 10% formalin for 40 min at RT. Fixed cells were incubated with Oil-Red O staining solution (Sigma Aldrich) for 10 min and then rinsed in 60% isopropanol to remove excessive Oil-Red O. Following imaging, Oil Red O stain was solubilized in isopropanol, and absorbance measures were taken at 510 nm utilizing Asys UVM340 microtiter plate reader.

Radioligand ADRB3-binding competition assays

Radioligand-binding competition activity was tested on the recombinant adrenergic β3 receptor for CPF and CPF oxon utilizing filtration binding assays by Epics and Euroscreen FAST (Belgium). Assays were performed in a 96-well plate (Master Block, Greiner, 786201), in duplicate. Binding buffer (25 mM Hepes pH 7.4, 1 mM EDTA, 0.5% BSA, 10 µg/ml saponin), membranes prepared from β3-receptor-expressing CHO-K1 cells (2 µg protein/well), radiotracer ¹²⁵I-Cyanopindolol (0.165 nM final assay concentration) and test compound were utilized. A total of 9 concentrations were assessed per test compound between 0.0001 and 10,000 nM in duplicate. The assay was performed using reference competitor ZD 7114 (Tocris Pharmaceuticals) a β3 agonist. Nonspecific binding was determined by co-incubation with 200-fold excess of a cold competitors. The samples were incubated in a final volume of 0.1 ml at 25 °C and for 30 min and then filtered over GF/B Unifilter plates presoaked for 2 h in 0.5% BSA. Filters are washed nine times with 0.5 ml of ice-cold washing buffer (25 mM Hepes pH 7.4, 0.5 M NaCl) and 50 µl of Microscint 20 (Packard) are added to each well. The plates were incubated for 15 min on an orbital shaker and then counted with a TopCountTM for 1 min/well. Dose–response data from test compounds were analyzed with XLfit (IDBS) software using nonlinear regression applied to a sigmoidal dose–response model.

In-vivo analysis

Housing and diet

All experiments were approved by the McMaster University Animal Ethics Committee and conducted under the Canadian guidelines for animal research. All mice used in this study were ordered at 7-week-old C57BL/6J males (JAX^®, The Jackson Laboratory). They were housed in specific pathogen-free microisolator cages located in a room maintained at a constant temperature of 22 (room temperature, RT) or 30 °C (thermoneutrality, TN) on a 12-h light–dark cycle with lights on at 7:00 a.m., the humidity was 40–60%. The mice were fed a control diet (CD, 10 kcal% fat; D12450H research diet) or a high-fat diet (HFD, 45 kcal% fat; D12451 research diet) for 14 weeks. To avoid AChE inhibition⁷⁸ and mimic the level of environmental CPF exposure, the low dose CPF diet was supplemented with 2 mg/kg/DT (per kg diet) CPF for the first week; 3 mg/DT diet CPF (Toronto Research Chemicals, C425300, ON) for the second and third weeks and 5 mg/kg/DT CPF for 11 weeks, to provide 0.5 mg/kg/BW (per kg body weight) CPF. High dose CPF was 4 times the low dose. Based on food intake the actual CPF intake of each group was lower than the designed doses (Supplementary Fig. 2a). Diets and water were provided to the animal’s ad libitum.

Body composition measures

Bodyweight (per mouse, n = 10) and body composition via Bruker’s minispec Whole Body Composition Analyzer software and (Bruker-LF90-MRI) were monitored weekly.

Metabolic testing

Glucose- and insulin-tolerance tests (GTT and ITT) were performed after 12 and 13 weeks of treatment respectively. Mice were fasted for 6 h and intraperitoneally injected with glucose (1 g/kg in saline) or insulin (0.7 U/kg for CD mice and 1 U/kg for HFD mice). Blood samples were collected by tail vein bleed and analyzed using a glucometer (Accu-Chek^®, Roche, 04680448003) immediately before and at 20, 40, 60, 90, and 120 min after injection of glucose or insulin.

Metabolic cages, Columbus Instruments Comprehensive Lab Animal Monitoring System (CLAMS), and Oxymax software (Columbus Instruments) were used to obtain VO₂, VCO₂, food intake, and activity data.

Diet-induced thermogenesis measures

Mice were fed with HFD (45 kcal% fat) supplemented with or without 20 mg/kg CPF for 5 weeks. Prior to 12 h of fasting, mice were placed in the CLAMS for 12 h for acclimatization and motoring the basal metabolic rate and activities from 7 p.m. to 7 a.m. Mice were then fasted from 7 a.m. to 7 p.m. and re-fed with 2 g of respective diet at 7 p.m. the following day.

Tissue processing, histological examination, and transmission electron microscopy (TEM)

Tissues for histology were fixed in 10% neutral buffered formalin for 12 h at 4 °C, embedded in paraffin, and sectioned (5 µm thickness). Following deparaffinization, tissue sections were used for Hematoxylin-eosin (H&E) staining or immunostaining. The size of adipocytes was measured (three images per mouse) using Image J (NIH), representative images were selected according to the average adipocyte size. For immunostaining, sections were heated in citrate buffer for 20 min, blocked with 5% goat serum in TBS containing 0.3% Triton X-100 for 2 h, then incubated sequentially with primary antibodies overnight and secondary antibody for 1 h. The color was developed using a VECTOR^® NovaRED^TM Peroxidase (HRP) Substrate Kit (#SK-4800, Vector Laboratories). Sections were then mounted with a mounting medium (Vector Lab, Burlingame, CA). For liver used for lipid staining, tissues were infiltrated with 30% sucrose for 24 h after formalin fixation, then embed with OTC and frozen in isopentane cooled in liquid nitrogen. Cryo-sections of 10-µm-thick were rinsed with PBS, then stained with Oil Red O for 10 min, washed with 60% isopropanol for three times then mounted with an aqueous mounting media for further imaging.

Cell fractions of BAT tissue were isolated using a Membrane and Cytosol Protein Extraction kit according to the manufacturer’s instructions (Beyotime, Cat. No. P0033, Shanghai, China).

For transmission electron microscopy (TEM), four interscapular BAT in each group were randomly selected and fixed in 2% glutaraldehyde (2% v/v) in 0.1 M sodium cacodylate buffer (pH 7.4) for 24 h. Tissues were sectioned on a Leica UCT ultramicrotome and picked up onto Cu grids. Sections were post-stained with uranyl acetate and lead citrate. The sectioning was performed by the electron microscopy group at McMaster University Medical Center. Mitochondria sizes was analyzed using Image J (NIH). The number of total mitochondria and mitochondria with disrupted cristae were counted. Cristae with any observable disorganization, vacuolization, or dissolution of cristae within mitochondria were categorized as disrupted. For each animal, more than 20 images were taken and analyzed. Researchers were blinded during image capturing and quantification. Representative images were selected according to the quantification data.

Pair-feeding experiments

Food intake of the control animals was measured every 3 days and the amount of food was provided to the pair-fed group in the following 3 days. For most of the time, the pair-fed animals treated with CPF ingested all the provided food. The difference in cumulative food intake of control and pair-fed mice was <0.1 g.

AChE activity assay

Acetylcholine esterase activity was measured in mouse skeletal muscle (hindlimb) and serum using Acetylcholinesterase Activity Assay Kit (Sigma-Aldrich Cat# MAK119). Skeletal muscle samples were homogenized in 0.1 M phosphate buffer (pH 7.5), centrifuged at 12,100×g for 5 min, and the supernatant was obtained. Diluted serum samples or supernatant from skeletal muscle samples was utilized as per the manufacturer’s instruction. AChE activity was then normalized to protein content in the homogenates determined using the PierceTM BCA Protein Assay Kit. The assay utilizes the Ellman method in which thiocholine reacts with 5,5-dithiobis (2-nitrobenzoic acid) to produce a colorimetric product indicative of AChE activity.

Immunoblotting analysis

Proteins were extracted from tissue or cells using cell lysis buffer (50 mM HEPES pH 7.4, 150 mM NaCl, 100 mM NaF, 10 Napyrophosphate, 5 EDTA mM, 250 mM sucrose, 1 mM DTT, and 1 mM Na-orthovanadate, 1% Triton X and Complete protease inhibitor cocktail (Roche)). Samples (1 μg/μL) were prepared in 4× SDS sample buffer and boiled at 95 °C for 5 min. Sample proteins were separated in a 10% or 12% SDS–PAGE gel and transferred to a PVDF membrane at 100 V for 90 min. Membranes were blocked with 5% skim milk in TBST (50 mM Tris–HCl, 150 mM NaCl, 0.05 % Tween 20) for 1 h. Subsequently, membranes were subject to primary antibody (1:1000 dilution) in TBST (with 5% BSA) overnight at 4 °C and then incubated with secondary antibody (1:1000 dilution) at room temperature for 1 h. Protein bands were imaged using electrochemiluminescence and analyzed using Image J software (National Institutes of Health, Bethesda, USA). Antibodies against pACC^S79 (#3661), ACC (#3662), pAMPKα^T172 (#2535), AMPKα (#2532), pULK1^S555 (#5869), ULK1 (#8054), p-HSL^S660 (#4126), HSL (#4107), LC3B (#2775), p38 (#8690), p-p38 (#9211), p62 (#5114), pPKA substrate (#9621) and p62 (#5114) were purchased from Cell Signaling (Danvers, MA). Antibody against UCP1 (#UCP11-A) was purchased from Alpha Diagnostic International. Antibody against β-tubulin (#32-2600) was purchased from Invitrogen.

Tissue and serum triglyceride determination

Approximately 50 mg of tissue was homogenized in chloroform:methanol (2:1) and a portion of the organic phase was dried down and suspended in isopropanol. Serum or tissue samples were assayed for triglyceride amount using a Triglyceride Colorimetric Assay Kit (Cayman Chemical, MI).

RNA isolation and quantitative real-time PCR (qRT-PCR)

Total RNA was isolated using TRIzol (Life Technologies, Grand Island, NY, USA) and purified using an RNeasy kit (QIAGEN) column. cDNA was synthesized using iScript™ cDNA Synthesis Kit (Bio-Rad, Hercules, CA). qRT-PCR was performed in a qPCR thermocycler (Corbett Rotor Gene 6000, MBI, QC, Canada) using TaqMan primers purchased from Invitrogen (Supplementary Table 1). Relative gene expression was calculated using the comparative Ct (2^−ΔΔCt) method, where values were normalized to a housekeeping gene (Ppia).

In-silico analysis

Docking of CPF to β₃-AR

To predict the binding pocket and binding mode of CPF in human beta-3 adrenergic receptors (β₃-AR), we employed an approach combining molecular dynamics (MD) simulations and docking calculations. Since the experimental 3-D structure of β₃-AR is unavailable, we used the multiple sequence alignment approach to build a homology model. a β₃-AR homology model was used as a starting structure to carry out classical MD simulation. The MD simulations were performed using GROMACS v5.1.2⁷⁹. AutoDock Vina⁸⁰ was then used to dock CPF to the 100 snapshots each of β₃AR extracted from the MD simulation to determine the pocket and high scoring pose.

β₃-AR multiple sequence alignment and model building

The amino acid sequence of β₃-AR was downloaded from UniProt (http://www.uniprot.org/) using accession number P13945. A template structure search was first carried out using SWISS-MODEL (https://swissmodel.expasy.org/). The template structures with a sequence identity > 50% and high GMQE score was selected for building the three-dimensional (3-D) structure of β₃-AR. The selected template structures were chain B of activated turkey β₁-AR (56.33%; PDB ID: 6H7J)⁸¹ and chain A of human β₂-AR (50%; PDB ID: 3NY8)⁸². MODELLER v9.20⁸³ was then used to perform multiple sequence alignment and generate the 3-D structure of β₃-AR. Prior to multiple sequence alignment, the G protein-like antibody, T4-lysozyme, seen in the template structure of human β₂-AR was removed. Five 3-D structures of β₃-AR were generated. Stereochemical checks were performed using the SWISS-MODEL structure assessment tool (https://swissmodel.expasy.org/) on all five structures and the best was chosen based on minimal stereochemical deviation and low Molprobity score. The final selected model had 96.74% of amino acid residues in the Ramachandran favored region with no outlier. Molprobity score, which combines the clash score, rotamer, and Ramachandran evaluations into a single score, normalized to be on the same scale as X-ray resolution for the model was 2.80. root mean square deviation (RMSD) between backbone atoms of final β₃-AR model and β₂-AR experimental structure was 3.6 Å. The final modeled structure of β₃-AR lacked N-terminal amino acid residues 1–32, intracellular loop 237–284, and C-terminal amino acids 361–408.

System setup and molecular dynamics simulation

The homology model of β₃-AR was used as a starting structure to carry out MD simulation. The structure file was first uploaded on the web server, Prediction of Proteins in Membranes (http://opm.phar.umich.edu/server.php). Membrane boundaries provided by this server along with the protein model were then uploaded onto the CHARMM-GUI server (http://charmm-gui.org/) for further processing. Protonation states of amino acid residues were assigned at the physiological pH of 7.4. The protein was embedded in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine. The replacement method was used to pack the receptor within the lipid bilayer. The lipid layer thickness was chosen to be 1.6 (80 lipids in a top leaflet and 78 lipids in the bottom leaflet). This system was then placed in a rectangular solvent box, and a 15 Å TIP3P water layer was added to solvate intra-and extra-cellular space. Charge neutrality of the system was achieved by adding Na+ and Cl− ions at a concentration of 0.15 mol/L to the water layers. 10 ns MD simulation of this system was then performed using the protocol described in Patel et al.⁸⁴. During the 10 ns production simulation, snapshots were saved every 100 ps giving 100 snapshots of β₃-AR to be used for docking calculations.

Docking of CPF to β₃-AR

Each of the 100 snapshots obtained during MD simulations was used to dock CPF using AutoDock Vina⁸⁰. The docking program allows the ligand to be completely flexible during the conformational search. However, only a few restricted side chains on the protein can be assigned as flexible. The use of 100 snapshots from MD simulation allows us to at least partly overcome this limitation.

CPF was sketched in 2-D using the ChemDraw tool of the ChemBio office package (http://www.cambridgesoft.com). The 2-D structure was converted to 3-D using the Chem3D tool. Ligand geometry was then optimized by performing an energy minimization run using MMFF94force-field. Lastly, the optimized 3-D structure of CPF was converted to the PDBQT format using AutoDockTools (http://mgltools.scripps.edu/) with the default identification of rotatable bonds and Gasteiger partial charges.

Ligand docking was then carried out by creating a grid box of size 28 Å × 28 Å × 28 Å, centered at a known ligand-binding site, with a grid spacing of 1 Å. This grid was large enough to allow exploration of all possible CPF-binding orientations surrounding known AR-binding pocket in an unbiased manner. The input exhaustiveness parameter for the docking was set to 500. The number of top docking orientations with high docking scores was fixed to 15. This docking protocol was similarly applied to all 100 snapshots of β₃-AR. CPF docking pose with the highest docking score was extracted and analyzed.

Statistical analysis

All data were found to be normally distributed. Results were analyzed using Student’s t-test or ANOVA where appropriate, using GraphPad Prism software. A repeated-measures ANOVA was used for all bodyweight plots, fed blood glucose, and GTT and ITT data. A Bonferroni post hoc test was used to test for significant differences revealed by the ANOVA. Significance was accepted at p ≤ 0.05.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

RNA-sequencing data are deposited in Gene Expression Omnibus (GEO) under accession number GSE178366. All data supporting the findings of this study are provided within the paper and its supplementary information. All additional information will be made available upon reasonable request to the corresponding author. Source data are provided with this paper.

References

Malik, V. S., Willett, W. C. & Hu, F. B. Global obesity: trends, risk factors and policy implications. Nat. Rev. Endocrinol. 9, 13–27 (2013).
Article PubMed Google Scholar
Hill, J. O. Can a small-changes approach help address the obesity epidemic? A report of the Joint Task Force of the American Society for Nutrition, Institute of Food Technologists, and International Food Information Council. Am. J. Clin. Nutr. 89, 477–484 (2009).
Article CAS PubMed Google Scholar
Rothwell, N. J. & Stock, M. J. A role for brown adipose tissue in diet-induced thermogenesis. Nature 281, 31 (1979).
Article ADS CAS PubMed Google Scholar
Vosselman, M. J. et al. Brown adipose tissue activity after a high-calorie meal in humans. Am. J. Clin. Nutr. 98, 57–64 (2013).
Article CAS PubMed Google Scholar
M, U. D. et al. Postprandial oxidative metabolism of human brown fat indicates thermogenesis. Cell Metab. 28, 207–216 e203 (2018).
Article CAS Google Scholar
Hibi, M. et al. Brown adipose tissue is involved in diet-induced thermogenesis and whole-body fat utilization in healthy humans. Int. J. Obes. 40, 1655–1661 (2016).
Article CAS Google Scholar
Feldmann, H. M., Golozoubova, V., Cannon, B. & Nedergaard, J. UCP1 ablation induces obesity and abolishes diet-induced thermogenesis in mice exempt from thermal stress by living at thermoneutrality. Cell Metab. 9, 203–209 (2009).
Article CAS PubMed Google Scholar
Nagai, N., Sakane, N., Ueno, L. M., Hamada, T. & Moritani, T. The -3826 A -> G variant of the uncoupling protein-1 gene diminishes postprandial thermogenesis after a high fat meal in healthy boys. J. Clin. Endocr. Metab. 88, 5661–5667 (2003).
Article CAS PubMed Google Scholar
Jia, J. J. et al. The polymorphisms of UCP1 genes associated with fat metabolism, obesity and diabetes. Mol. Biol. Rep. 37, 1513–1522 (2010).
Article CAS PubMed Google Scholar
Yoneshiro, T. et al. Impact of UCP1 and beta3AR gene polymorphisms on age-related changes in brown adipose tissue and adiposity in humans. Int. J. Obes. 37, 993–998 (2013).
Article CAS Google Scholar
Sramkova, D. et al. The UCP1 gene polymorphism A-3826G in relation to DM2 and body composition in Czech population. Exp. Clin. Endocrinol. Diabetes 115, 303–307 (2007).
Article CAS PubMed Google Scholar
Tang, L. et al. Meta-analyses between 18 candidate genetic markers and overweight/obesity. Diagn. Pathol. 9, 56 (2014).
Article PubMed PubMed Central CAS Google Scholar
Leitner, B. P. et al. Mapping of human brown adipose tissue in lean and obese young men. Proc. Natl Acad. Sci. USA 114, 8649–8654 (2017).
Article CAS PubMed PubMed Central Google Scholar
Liang, Y. et al. Organophosphorus pesticide chlorpyrifos intake promotes obesity and insulin resistance through impacting gut and gut microbiota. Microbiome 7, 19 (2019).
Article PubMed PubMed Central Google Scholar
Baillie-Hamilton, P. F. Chemical toxins: a hypothesis to explain the global obesity epidemic. J. Altern. Complement. Med. 8, 185–192 (2002).
Article PubMed Google Scholar
Singh, K. & Chan, H. M. Persistent organic pollutants and diabetes among Inuit in the Canadian Arctic. Environ. Int. 101,183–189 (2017).
Kassotis, C. D. & Stapleton, H. M. Endocrine-mediated mechanisms of metabolic disruption and new approaches to examine the public health threat. Front. Endocrinol. 10, 39 (2019).
Article Google Scholar
Di Gregorio, I. et al. Environmental pollutants effect on brown adipose tissue. Front. Physiol. 9, 1891 (2018).
Article PubMed Google Scholar
Heindel, J. J. History of the obesogen field: looking back to look forward. Front. Endocrinol. 10, 14 (2019).
Article Google Scholar
Gutgesell, R. M., Tsakiridis, E. E., Jamshed, S., Steinberg, G. R. & Holloway, A. C. Impact of pesticide exposure on adipose tissue development and function. Biochem. J. 477, 2639–2653 (2020).
Article CAS PubMed Google Scholar
Castriota, F. et al. Chronic arsenic exposure impairs adaptive thermogenesis in male C57BL/6J mice. Am. J. Physiol. Endocrinol. Metab. 318, E667–E677 (2020).
Article CAS PubMed PubMed Central Google Scholar
Bae, J. et al. Arsenite exposure suppresses adipogenesis, mitochondrial biogenesis and thermogenesis via autophagy inhibition in brown adipose tissue. Sci. Rep. 9, 14464 (2019).
Article ADS PubMed PubMed Central CAS Google Scholar
Zuo, Z. et al. Prolonged inorganic arsenic exposure via drinking water impairs brown adipose tissue function in mice. Sci. Total Environ. 668, 310–317 (2019).
Article ADS CAS PubMed Google Scholar
Xu, C.-X. et al. Aryl hydrocarbon receptor deficiency protects mice from diet-induced adiposity and metabolic disorders through increased energy expenditure. Int. J. Obes. 39, 1300 (2015).
Article CAS Google Scholar
La Merrill, M. et al. Perinatal exposure of mice to the pesticide DDT impairs energy expenditure and metabolism in adult female offspring. PLoS ONE 9, e103337 (2014).
Article ADS PubMed PubMed Central CAS Google Scholar
Virtanen, K. A. et al. Functional brown adipose tissue in healthy adults. N. Engl. J. Med. 360, 1518–1525 (2009).
Article CAS PubMed Google Scholar
Ouellet, V. et al. Outdoor temperature, age, sex, body mass index, and diabetic status determine the prevalence, mass, and glucose-uptake activity of 18F-FDG-detected BAT in humans. J. Clin. Endocrinol. Metab. 96, 192–199 (2011).
Article CAS PubMed Google Scholar
Saito, M. et al. High incidence of metabolically active brown adipose tissue in healthy adult humans: effects of cold exposure and adiposity. Diabetes 58, 1526–1531 (2009).
Article CAS PubMed PubMed Central Google Scholar
van der Lans, A. A. et al. Cold acclimation recruits human brown fat and increases nonshivering thermogenesis. J. Clin. Investig. 123, 3395–3403 (2013).
Article PubMed PubMed Central CAS Google Scholar
Harms, M. & Seale, P. Brown and beige fat: development, function and therapeutic potential. Nat. Med. 19, 1252–1263 (2013).
Article CAS PubMed Google Scholar
Bartelt, A. & Heeren, J. Adipose tissue browning and metabolic health. Nat. Rev. Endocrinol. 10, 24–36 (2014).
Article CAS PubMed Google Scholar
Galmozzi, A. et al. ThermoMouse: an in vivo model to identify modulators of UCP1 expression in brown adipose tissue. Cell Rep. 9, 1584–1593 (2014).
Article CAS PubMed PubMed Central Google Scholar
Solomon, K. R. et al. Properties and uses of chlorpyrifos in the United States. Rev. Environ. Contam. Toxicol. 231, 13–34 (2014).
CAS PubMed Google Scholar
Huff, R. A., Abu-Qare, A. W. & Abou-Donia, M. B. Effects of sub-chronic in vivo chlorpyrifos exposure on muscarinic receptors and adenylate cyclase of rat striatum. Arch. Toxicol. 75, 480–486 (2001).
Article CAS PubMed Google Scholar
Akhtar, N., Srivastava, M. K. & Raizada, R. B. Assessment of chlorpyrifos toxicity on certain organs in rat, Rattus norvegicus. J. Environ. Biol. 30, 1047–1053 (2009).
CAS PubMed Google Scholar
Peris-Sampedro, F. et al. Chronic exposure to chlorpyrifos triggered body weight increase and memory impairment depending on human apoE polymorphisms in a targeted replacement mouse model. Physiol. Behav. 144, 37–45 (2015).
Article CAS PubMed Google Scholar
Meggs, W. J. & Brewer, K. L. Weight gain associated with chronic exposure to chlorpyrifos in rats. J. Med. Toxicol. 3, 89–93 (2007).
Article PubMed PubMed Central Google Scholar
Fang, B., Li, J. W., Zhang, M., Ren, F. Z. & Pang, G. F. Chronic chlorpyrifos exposure elicits diet-specific effects on metabolism and the gut microbiome in rats. Food Chem. Toxicol. 111, 144–152 (2018).
Article CAS PubMed Google Scholar
Wang, L., Liu, Z., Zhang, J., Wu, Y. & Sun, H. Chlorpyrifos exposure in farmers and urban adults: metabolic characteristic, exposure estimation, and potential effect of oxidative damage. Environ. Res. 149, 164–170 (2016).
Article CAS PubMed Google Scholar
Kazak, L. et al. UCP1 deficiency causes brown fat respiratory chain depletion and sensitizes mitochondria to calcium overload-induced dysfunction. Proc. Natl Acad. Sci. USA 114, 7981–7986 (2017).
Article CAS PubMed PubMed Central Google Scholar
Bond, L. M., Burhans, M. S. & Ntambi, J. M. Uncoupling protein-1 deficiency promotes brown adipose tissue inflammation and ER stress. PLoS ONE 13, e0205726 (2018).
Article PubMed PubMed Central CAS Google Scholar
Xu, M. Y., Wang, P., Sun, Y. J., Yang, L. & Wu, Y. J. Joint toxicity of chlorpyrifos and cadmium on the oxidative stress and mitochondrial damage in neuronal cells. Food Chem. Toxicol. 103, 246–252 (2017).
Article CAS PubMed Google Scholar
Shih, D. M. et al. Mice lacking serum paraoxonase are susceptible to organophosphate toxicity and atherosclerosis. Nature 394, 284–287 (1998).
Article ADS CAS PubMed Google Scholar
Riederer, A. M. et al. Pyrethroid and Organophosphorus Pesticides in Composite Diet Samples from Atlanta, USA Adults. Envir Sci. Tech. 44, 483–490 https://doi.org/10.1021/es902479h (2010).
Article CAS Google Scholar
Fischer, A. W., Cannon, B. & Nedergaard, J. Optimal housing temperatures for mice to mimic the thermal environment of humans: an experimental study. Mol. Metab. 7, 161–170 (2018).
Article CAS PubMed Google Scholar
Giles, D. A. et al. Thermoneutral housing exacerbates nonalcoholic fatty liver disease in mice and allows for sex-independent disease modeling. Nat. Med. 23, 829–838 (2017).
Article CAS PubMed PubMed Central Google Scholar
Tian, X. Y. et al. Thermoneutral housing accelerates metabolic inflammation to potentiate atherosclerosis but not insulin resistance. Cell Metab. 23, 165–178 (2016).
Article CAS PubMed Google Scholar
Cui, X. et al. Thermoneutrality decreases thermogenic program and promotes adiposity in high-fat diet-fed mice. Physiol. Rep. 4, e12799 (2016).
Peris-Sampedro, F. et al. Adulthood dietary exposure to a common pesticide leads to an obese-like phenotype and a diabetic profile in apoE3 mice. Environ. Res. 142, 169–176 (2015).
Article CAS PubMed Google Scholar
Fischer, A. W., Schlein, C., Cannon, B., Heeren, J. & Nedergaard, J. Intact innervation is essential for diet-induced recruitment of brown adipose tissue. Am. J. Physiol.-Endocrinol. Metab. 316, E487–E503 (2018).
Article PubMed PubMed Central CAS Google Scholar
Bachman, E. S. et al. betaAR signaling required for diet-induced thermogenesis and obesity resistance. Science 297, 843–845 (2002).
Article ADS CAS PubMed Google Scholar
Steinberg, G. R. Cellular energy sensing and metabolism-implications for treating diabetes: the 2017 outstanding scientific achievement award lecture. Diabetes 67, 169–179 (2018).
Article CAS PubMed Google Scholar
Desjardins, E. M. & Steinberg, G. R. Emerging role of AMPK in brown and beige adipose tissue (BAT): implications for obesity, insulin resistance, and type 2 diabetes. Curr. Diab. Rep. 18, 80 (2018).
Article PubMed CAS Google Scholar
Dempersmier, J. et al. Cold-inducible Zfp516 activates UCP1 transcription to promote browning of white fat and development of brown fat. Mol. Cell 57, 235–246 (2015).
Article CAS PubMed PubMed Central Google Scholar
Mottillo, E. P. et al. Lack of adipocyte AMPK exacerbates insulin resistance and hepatic steatosis through brown and beige adipose tissue function. Cell Metab. 24, 118–129 (2016).
Article CAS PubMed PubMed Central Google Scholar
Myers, R. W. et al. Systemic pan-AMPK activator MK-8722 improves glucose homeostasis but induces cardiac hypertrophy. Science 357, 507–511 (2017).
Article ADS CAS PubMed Google Scholar
Jackson, E., Shoemaker, R., Larian, N. & Cassis, L. Adipose tissue as a site of toxin accumulation. Compr. Physiol. 7, 1085–1135 (2017).
Article PubMed PubMed Central Google Scholar
Tanvir, E. M. et al. A model of chlorpyrifos distribution and its biochemical effects on the liver and kidneys of rats. Hum. Exp. Toxicol. 35, 991–1004 (2016).
Article CAS PubMed Google Scholar
Burke, R. D. et al. Developmental neurotoxicity of the organophosphorus insecticide chlorpyrifos: from clinical findings to preclinical models and potential mechanisms. J. Neurochem. 142, 162–177 (2017).
Article CAS PubMed PubMed Central Google Scholar
Ganeshan, K. & Chawla, A. Warming the mouse to model human diseases. Nat. Rev. Endocrinol. 13, 458–465 (2017).
Article PubMed PubMed Central Google Scholar
Meyer, A. et al. Critical periods for chlorpyrifos-induced developmental neurotoxicity: alterations in adenylyl cyclase signaling in adult rat brain regions after gestational or neonatal exposure. Environ. Health Perspect. 112, 295–301 (2004).
Article CAS PubMed PubMed Central Google Scholar
Jiang, W. et al. Mice treated with chlorpyrifos or chlorpyrifos oxon have organophosphorylated tubulin in the brain and disrupted microtubule structures, suggesting a role for tubulin in neurotoxicity associated with exposure to organophosphorus agents. Toxicol. Sci. 115, 183–193 (2010).
Article CAS PubMed PubMed Central Google Scholar
Peiris, D. C. & Dhanushka, T. Low doses of chlorpyrifos interfere with spermatogenesis of rats through reduction of sex hormones. Environ. Sci. Pollut. Res. 24, 20859–20867 (2017).
Article CAS Google Scholar
Hazarika, J. et al. Endocrine disruption: molecular interactions of chlorpyrifos and its degradation products with estrogen receptor. Struct. Chem. 31, 2011–2021 (2020).
Article CAS Google Scholar
Comfort, N. & Re, D. B. Sex-specific neurotoxic effects of organophosphate pesticides across the life course. Curr. Environ. Health Rep. 4, 392–404 (2017).
Article CAS PubMed PubMed Central Google Scholar
Butter, M. E. Are women more vulnerable to environmental pollution? J. Hum. Ecol. 20, 221–226 (2006).
Article Google Scholar
Thomas, A. et al. Concentrations of organochlorine pesticides in pooled human serum by age and gender. Environ. Res. 154, 10–18 (2017).
Article CAS PubMed Google Scholar
Rodríguez, A. M., Quevedo‐Coli, S., Roca, P. & Palou, A. Sex‐dependent dietary obesity, induction of UCPs, and leptin expression in rat adipose tissues. Obes. Res. 9, 579–588 (2001).
Article PubMed Google Scholar
van den Beukel, J. C. et al. Women have more potential to induce browning of perirenal adipose tissue than men. Obesity 23, 1671–1679 (2015).
Article PubMed CAS Google Scholar
Schmidt, S. L. et al. Adrenergic control of lipolysis in women compared with men. J. Appl. Physiol. 117, 1008–1019 (2014).
Article CAS PubMed PubMed Central Google Scholar
Morris, E. M. et al. Difference in housing temperature-induced energy expenditure elicits sex-specific diet-induced metabolic adaptations in mice. Obesity 28, 1922–1931 (2020).
Article CAS PubMed Google Scholar
Uldry, M. et al. Complementary action of the PGC-1 coactivators in mitochondrial biogenesis and brown fat differentiation. Cell Metab. 3, 333–341 (2006).
Article CAS PubMed Google Scholar
Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).
Article CAS PubMed PubMed Central Google Scholar
Kim, D., Paggi, J. M., Park, C., Bennett, C. & Salzberg, S. L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915 (2019).
Article CAS PubMed PubMed Central Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Article PubMed PubMed Central CAS Google Scholar
Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Article ADS CAS PubMed PubMed Central Google Scholar
Mootha, V. K. et al. PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat. Genet. 34, 267–273 (2003).
Article CAS PubMed Google Scholar
Shannon, P. et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 13, 2498–2504 (2003).
Article CAS PubMed PubMed Central Google Scholar
Berendsen, H. et al. GROMACS: fast, flexible, and free. J. Comput. Chem. 26, 1701–1718 (2005).
Article PubMed CAS Google Scholar
Vina, A. Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading Trott, Oleg; Olson, Arthur J. J. Comput. Chem. 31, 455–461 (2010).
Google Scholar
Warne, T., Edwards, P. C., Dore, A. S., Leslie, A. G. W. & Tate, C. G. Molecular basis for high-affinity agonist binding in GPCRs. Science 364, 775–778 (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Wacker, D. et al. Conserved binding mode of human [beta subscript 2] adrenergic receptor inverse agonists and antagonist revealed by X-ray crystallography. J. Am. Chem. Soc. 132, 11443–11445 (2010).
Webb, B. & Sali, A. Protein Structure Modeling with MODELLER. Methods Mol Biol. 2199, 239–255 (2021).
Patel, J. S., Brown, C. J., Ytreberg, F. M. & Stenkamp, D. L. Predicting peak spectral sensitivities of vertebrate cone visual pigments using atomistic molecular simulations. PLoS Comput. Biol. 14, e1005974 (2018).
Article ADS PubMed PubMed Central CAS Google Scholar

Download references

Acknowledgements

The funding for this study was provided by the Canadian Institutes of Health Research (CIHR) Programmatic Grants in Environment, Genes and Chronic Diseases (to A.G.M., M.G.W., K.M.M., A.C.H. and G.R.S.—Grant #144625-1). G.R.S. is supported by a Canada Research Chair and a J. Bruce Duncan Chair in Metabolic Diseases. E.M.D. is a Vanier Scholar. J.S.P. was supported by the Center for Modeling Complex Interactions sponsored by the NIGMS under award number P20 GM104420. Computational resources were provided by the Institute for Bioinformatics and Evolutionary Studies Computational Resources Core sponsored by the National Institutes of Health (NIH P30 GM103324) and by the High Performance Computing Center at Idaho National Laboratory, which is supported by the Office of Nuclear Energy of the U.S. DOE and the Nuclear Science User Facilities under Contract No. DE-AC07-05ID14517.

Author information

Authors and Affiliations

Centre for Metabolism, Obesity and Diabetes Research, McMaster University, Hamilton, ON, Canada
Bo Wang, Evangelia E. Tsakiridis, Shuman Zhang, Andrea Llanos, Eric M. Desjardins, Julian M. Yabut, Alexander E. Green, Emily A. Day, Brennan K. Smith, James S. V. Lally, Jianhan Wu, Amogelang R. Raphenya, Krishna A. Srinivasan, Andrew G. McArthur, Katherine M. Morrison, Alison C. Holloway & Gregory R. Steinberg
Division of Endocrinology and Metabolism, Department of Medicine, McMaster University, Hamilton, ON, Canada
Bo Wang, Evangelia E. Tsakiridis, Shuman Zhang, Eric M. Desjardins, Julian M. Yabut, Alexander E. Green, Emily A. Day, Brennan K. Smith, James S. V. Lally, Jianhan Wu & Gregory R. Steinberg
State Key Laboratory of Animal Nutrition, College of Animal Science and Technology, China Agricultural University, Beijing, PR China
Bo Wang
Department of Obstetrics and Gynecology, McMaster University, Hamilton, ON, Canada
Andrea Llanos & Alison C. Holloway
Department of Biochemistry and Biomedical Sciences, McMaster University, Hamilton, ON, Canada
Amogelang R. Raphenya, Krishna A. Srinivasan, Andrew G. McArthur & Gregory R. Steinberg
Beth Israel Deaconess Medical Center and Harvard Medical School, Boston, MA, USA
Shingo Kajimura
Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, ID, USA
Jagdish Suresh Patel
Department of Biological Sciences, University of Idaho, Moscow, ID, USA
Jagdish Suresh Patel
Environmental Health Science & Research Bureau, Health Canada, Ottawa, ON, Canada
Michael G. Wade
Department of Pediatrics, McMaster University, Hamilton, ON, Canada
Katherine M. Morrison

Authors

Bo Wang
View author publications
You can also search for this author in PubMed Google Scholar
Evangelia E. Tsakiridis
View author publications
You can also search for this author in PubMed Google Scholar
Shuman Zhang
View author publications
You can also search for this author in PubMed Google Scholar
Andrea Llanos
View author publications
You can also search for this author in PubMed Google Scholar
Eric M. Desjardins
View author publications
You can also search for this author in PubMed Google Scholar
Julian M. Yabut
View author publications
You can also search for this author in PubMed Google Scholar
Alexander E. Green
View author publications
You can also search for this author in PubMed Google Scholar
Emily A. Day
View author publications
You can also search for this author in PubMed Google Scholar
Brennan K. Smith
View author publications
You can also search for this author in PubMed Google Scholar
James S. V. Lally
View author publications
You can also search for this author in PubMed Google Scholar
Jianhan Wu
View author publications
You can also search for this author in PubMed Google Scholar
Amogelang R. Raphenya
View author publications
You can also search for this author in PubMed Google Scholar
Krishna A. Srinivasan
View author publications
You can also search for this author in PubMed Google Scholar
Andrew G. McArthur
View author publications
You can also search for this author in PubMed Google Scholar
Shingo Kajimura
View author publications
You can also search for this author in PubMed Google Scholar
Jagdish Suresh Patel
View author publications
You can also search for this author in PubMed Google Scholar
Michael G. Wade
View author publications
You can also search for this author in PubMed Google Scholar
Katherine M. Morrison
View author publications
You can also search for this author in PubMed Google Scholar
Alison C. Holloway
View author publications
You can also search for this author in PubMed Google Scholar
Gregory R. Steinberg
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

B.W., S.Z., A.L., E.T., E.M.D., J.Y., A.G., E.D., B.K., J.L., J.P., J.W. performed the experiments. B.W., E.M.D., A.G., E.D., J.P., A.R., K.S., A.M., S.K., M.W., K.M., A.H. and G.R.S. provided intellectual input. B.W., S.Z., A.L. and G.R.S. wrote the manuscript with contributions from all authors. E.M.D., B.W. and E.T. created the graphical abstract.

Corresponding author

Correspondence to Gregory R. Steinberg.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks Robert Sargis and the other anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Reporting Summary

Source data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Wang, B., Tsakiridis, E.E., Zhang, S. et al. The pesticide chlorpyrifos promotes obesity by inhibiting diet-induced thermogenesis in brown adipose tissue. Nat Commun 12, 5163 (2021). https://doi.org/10.1038/s41467-021-25384-y

Download citation

Received: 11 August 2020
Accepted: 02 August 2021
Published: 27 August 2021
DOI: https://doi.org/10.1038/s41467-021-25384-y

This article is cited by

Association of urinary chlorpyrifos, paraquat, and cyproconazole levels with the severity of fatty liver based on MRI
- Peiqi Ma
- Hongliang Gao
- Bin Xue
BMC Public Health (2024)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Subjects

Abstract

Similar content being viewed by others

Introduction

Results

Discussion

Methods

In-vitro analysis

Immortalized brown adipocyte cell line generation

Immortalized murine brown adipocyte culture

Cyclic AMP assay

MitoTracker staining

Cytochrome c oxidase activity

In vitro respiration assay of brown adipocytes

RNA-sequencing analysis

Oil red O staining and quantification

Radioligand ADRB3-binding competition assays

In-vivo analysis

Housing and diet

Body composition measures

Metabolic testing

Diet-induced thermogenesis measures

Tissue processing, histological examination, and transmission electron microscopy (TEM)

Pair-feeding experiments

AChE activity assay

Immunoblotting analysis

Tissue and serum triglyceride determination

RNA isolation and quantitative real-time PCR (qRT-PCR)

In-silico analysis

Docking of CPF to β3-AR

β3-AR multiple sequence alignment and model building

System setup and molecular dynamics simulation

Docking of CPF to β3-AR

Statistical analysis

Reporting summary

Data availability

References

Acknowledgements

Author information

Authors and Affiliations

Contributions

Corresponding author

Ethics declarations

Competing interests

Additional information

Supplementary information

Source data

Rights and permissions

About this article

Cite this article

Share this article

This article is cited by

Comments

Search

Quick links

Docking of CPF to β₃-AR

β₃-AR multiple sequence alignment and model building

Docking of CPF to β₃-AR