PPARγ protein expression in FGR placental derived EVs is impaired.
To explore the difference in human placental derived EVs between normal and FGR pregnancy, the EVs were firstly isolated from normal and FGR term placenta via sucrose gradient centrifugation. Electron microscopy showed that exosomes have a typical double-membrane structure (Fig. S1A). Western blot showed that EVs highly express their biomarkers, such as CD9, CD81, CD63 and HSP90β, while barely express their negative marker Tubulin (Fig. S1B). Moreover, we ascertained that the size of the extracted EVs from normal and FGR term placenta was primarily within the 80 to 130 nm range, as measured by laser scattering using the NanoSight tracking system (NTA, Fig. S1C). To identify the proteins enriched in human normal and FGR term placental derived EVs, we conducted proteomic analysis of EVs derived from human normal and FGR term placenta using GC-MS. Wayne's analysis revealed that 634 proteins were commonly expressed in both normal and FGR placental derived EVs (Fig. 1A). The GO-molecular function analysis of the commonly expressed proteins in normal and FGR term placental derived EVs indicated a predominant focus on protein binding and fatty acid metabolism (Fig. 1B). In our analysis of fatty acid metabolism-related proteins, we identified the PPARγ protein, which demonstrates a high expression abundance and plays a crucial role in adipogenesis (Table 1–2).
To confirm the expression of PPARγ in the human placenta, primary cytotrophoblasts (CTBs) were isolated from normal and FGR placenta tissues in late pregnancy using percoll gradient centrifugation. The purity of the extracted CTBs was validated through an immunofluorescence assay, demonstrating high levels of epithelial cell markers CK7 and E-cadherin (98.09 ± 0.56% and 97.89 ± 0.65% of the cell population, respectively), with minimal expression of the mesenchymal cell marker vimentin (1.85 ± 0.61%) (Fig. 1C). Subsequently, EVs were isolated from primary CTBs (Fig. S1D-F, supplement Table 5). Western blot results showed that the expression level of PPARγ protein in EVs of FGR derived CTBs and FGR CTBs was both significantly lower than that in the normal group (Fig. 1D and E), which is consistent with the lipogenesis disorder in FGR fetuses reported in previous study[5, 6]. Furthermore, a notable decrease in PPARγ protein expression was observed in EVs derived from FGR placental tissues compared to the normal group (Fig. 1F).
Table 5
Analysis of the free energy of binding between PPARγ and β’-COP
Energy component | Average | SD (Prop.) | SD | SEM (Prop.) | SEM |
ΔBOND | 0.00 | 18.51 | 0.00 | 0.19 | 0.00 |
ΔANGLE | 0.00 | 26.80 | 0.00 | 0.27 | 0.00 |
ΔDIHED | 0.00 | 14.86 | 0.00 | 0.15 | 0.00 |
ΔVDWAALS | -31.60 | 24.20 | 5.52 | 0.24 | 0.06 |
ΔEEL | 25.92 | 8.14 | 7.59 | 0.08 | 0.08 |
Δ1–4 VDW | 0.00 | 10.03 | 0.00 | 0.10 | 0.00 |
Δ1–4 EEL | 0.00 | 3.58 | 0.00 | 0.04 | 0.00 |
ΔEGB | -26.99 | 6.99 | 7.00 | 0.07 | 0.07 |
ΔESURF | -5.88 | 1.47 | 0.84 | 0.01 | 0.01 |
ΔGGAS | -5.69 | 25.53 | 10.10 | 0.26 | 0.10 |
ΔGSOLV | -32.87 | 7.15 | 6.53 | 0.07 | 0.07 |
ΔTOTAL | -38.55 | 26.51 | 5.77 | 0.27 | 0.06 |
SD: Standard Deviation; SEM: Standard Error of Mean; Prop.: Property; BOND: Bond Stretch Energy; ANGLE: Angle Bending Energy; DIHED: Dihedral Angle Energy; VDWAALS: van der Waals (vdW) Non-bonded Energy; EEL: Electrostatic (Coulombic) Energy between Charged Groups; VDW: van der Waals; EGB: Generalized Born Model for Solvation Energy; ESURF: Surface Area Term in Solvation Energy Calculation; GGAS: Gas Phase Energy Contribution to the Total Energy; GSOLV: Solvent Accessible Surface Area Energy; TOTAL: Total Energy Calculation. |
COPI coat complex subunit proteins were enriched in both normal and FGR trophoblast derived EVs.
As PPARγ is found within EVs, we investigated the mechanism by which PPARγ protein enters EVs in trophoblasts. Through proteomic analysis of EVs from normal and FGR term placenta, we identified proteins associated with intracellular protein transport, including COPB2 (β’-COP), ARCN1 (δ-COP), COPG1 (γ-COP), COPA (α-COP) and COPE (ε-COP), which are part of the COPI coat complex (Table 3–4). Western Blot analysis showed that both normal and FGR term placenta and placental-derived EVs express all COPI complex subunits (Fig. 1G-H). There was no significant difference in their expression levels between the two groups, suggesting that the transportation of PPARγ into EVs may be linked to the COPI coat complex.
β’-COP interacts directly with PPARγ in trophoblast starting from the early stage of EVs formation.
To investigate the role of the COPI complex in sorting PPARγ into EVs, a Co-immunoprecipitation (Co-IP) assay was conducted to confirm the interaction between COPI complex subunits and PPARγ. The interaction of each COPI complex subunit with PPARγ was individually examined in primary CTBs from both normal and FGR placenta. The results revealed a direct interaction between PPARγ and the β’-COP subunit in primary CTBs, with no involvement of other subunits (Fig. 2A-G). This association was further supported by reciprocal immunoprecipitation of COPI complex subunits with the PPARγ antibody (Fig. 2H). Immunofluorescence staining demonstrated co-localization of β’-COP and PPARγ in human placenta tissues and primary CTBs from both normal and FGR groups (Fig. 2I), providing spatial evidence of their binding. These findings indicate a specific interaction between β’-COP and PPARγ in human trophoblast cells.
Studies have shown that EV formation begins with the transition from early endosomes to late endosomes and multivesicular endosomes (LE/MVEs), culminating in EV production[27]. To investigate the stage at which β’-COP and PPARγ interact during EV formation, we performed an immunofluorescence assay to assess their co-localization with the early endosome marker EEA1 and the LE/MVE marker CD63 in HTR8/SVneo trophoblast cells. The assay confirmed the co-localization of PPARγ and β’-COP with EEA1 and CD63 (Fig. 2J-K). Furthermore, a Co-IP assay demonstrated a direct interaction between PPARγ and β’-COP in trophoblast cells (Fig. 2L). These results suggest that the interaction between PPARγ and β’-COP occurs early in EV formation. However, no co-localization was observed with LAMP1, a lysosomal marker, or LC3 puncta, an indicator of autophagosome, indicating that PPARγ and β’-COP specifically localize to early endosomes and MVEs rather than lysosomes or autophagosomes (Fig. 2M-N). In conclusion, β’-COP directly binds to PPARγ in trophoblast cells, influencing PPARγ localization to early endosomes and MVEs during EV formation.
β’-COP mediates ab initio loading of PPARγ into t-Exos.
Subsequently, we investigated the role of β’-COP in loading PPARγ into EVs. Through lentiviral transfection, we successfully generated a β’-COP knockout (β’-COP-KO) HTR8/SVneo trophoblast cell line (Fig. 3A-C). The EVs derived from HTR8/SVneo were also characterized (Fig. S1G-I, supplement Table 5). However, immunofluorescence staining revealed the lack of co-localization between PPARγ with EEA1and CD63 following β’-COP knockout (Fig. 3D-E). Western Blot analysis demonstrated a significant reduction in PPARγ protein expression in trophoblast EVs in the β’-COP-KO group compared to the control group (Fig. 3F). These findings suggest that β’-COP facilitates the transport of PPARγ into EVs during the early stages of EV formation.
β’-COP mediated transportation of PPARγ into EVs relies on intact COPI complex.
To investigate the relationship between β’-COP mediated loading of PPARγ into EVs and the COPI complex, we constructed δ-COP knockout (δ-COP-KO) and ε-COP knockout (ε-COP-KO) HTR8/SVneo trophoblast cell lines by lentiviral transfection (Fig. 3G-L). Co-IP assays indicated that the knockout of δ-COP or ε-COP could significantly interfere with the interaction of PPARγ to β’-COP (Fig. 3M-N). Importantly, the deficiency of δ-COP or ε-COP can also result in an inhibition in the EVs transportation of PPARγ in trophoblast which manifested as declined co-localization of PPARγ with EEA1 and CD63, as well as reduced PPARγ protein expression in δ-COP-KO and ε-COP-KO HTR8/SVneo trophoblast (Fig. 3O-S). These findings demonstrated the essential role of intact COPI complex in the loading of PPARγ proteins into EVs.
The binding of PPARγ to β’-COP depends on solvation free energy and van der Waals forces.
To investigate the binding modes between PPARγ and β’-COP, we conducted molecular dynamics simulations using Gromacs 2019.6 software. The root mean square deviation (RMSD) analysis revealed that equilibrium in the β’-COP-PPARγ binding was reached after 20 ns of simulation (Fig. 4A), with the RMSD value stabilizing around 1.0 nm. Further assessment of the stability of the complex was performed through radius of gyration (Rg) and solvent accessible surface area (SASA) analysis (Fig. 4B-C). The Rg fluctuated within 4.7 nm post 20 ns of simulation, indicating sustained stability of the β’-COP-PPARγ binding complex (Fig. 4B). Additionally, the SASA of the complex remained constant throughout the simulation (after 20 ns) (Fig. 4C), suggesting persistent binding of PPARγ to the β’-COP surface (Fig. 4D-E).
To investigate the binding force between β’-COP and PPARγ, we conducted simulations to analyze hydrogen bond distribution during the binding process (Fig. 4F). The results show that β’-COP-PPARγ typically forms 7–8 hydrogen bonds over the simulation period of 0-100ns, indicating stable binding of PPARγ to the β’-COP protein surface (Fig. 4F-H). Upon further analysis, the binding free energy (BFE) between PPARγ and β’-COP was determined to be -38.55 kcal/mol, primarily driven by solvation free energy (ΔGSOLV, -32.87 kcal/mol) and Vander Waals forces (ΔVDWAALS, -31.60 kcal/mol) (Table 5). The Coulomb force was found to be the main repulsive force (ΔEEL, 25.92 kcal/mol) (Fig. 4I, Table 5). Decomposition of the BFE of the β’-COP-PPARγ complex into residues highlighted key amino acids crucial for their interaction (Fig. 4J-K). Nine of these amino acids were identified as hot spots with a binding energy less than − 1 kcal/mol, namely GLU175, HIS178, ASP203, ARG204, ASP207, ASN238, and ARG244 for β’-COP, and LYS273, GYN281, and PHE285 for PPARγ (Fig. 4J-K). This analysis and subsequent bioinformatics investigation identified the binding domains of PPARγ and β’-COP to be located between amino acid residues 258–286 on PPARγ.
ASP258, LYS273, GLN281, GLY282, PHE285 and ARG286 of PPARγ protein were the key binding sites of β’-COP-PPARγ complex.
In order to identify the essential amino acid binding sites for the PPARγ-β’-COP complex within the 258–286 binding domains, an analysis was conducted to assess the binding energies contributed by amino acid sites within this domain. Through mutation of all amino acids in the binding domain to alanine, we observed a significant impact on the BFE at positions ASP258, ILE279, PHE280, GLY282, PHE285, and ARG286 (Table 6). Additionally, we investigated the effect of amino acid residue mutation on protein stability and BFE using ΔΔG, revealing that mutations at ASP258, LYS273, GLN281, GLY282, GLN284, PHE285, and ARG286 significantly influence the binding stability of the β’-COP-PPARγ complex (Table 7). Overall, our findings highlight the critical role of ASP258, LYS273, GLN281, GLY282, PHE285, and ARG286 in facilitating the transport of PPARγ into EVs mediated by β’-COP.
Table 6
Analysis of the BFE between PPARγ and β’-COP after amino acid sites mutation
Amino acid site | Amino acid | BFE before mutation | BFE after mutation |
258 | ASP | -0.752 | 2.796 |
273 | LYS | -3.436 | -1.673 |
274 | GLU | -0.390 | -0.334 |
275 | VAL | -0.288 | 0.0822 |
276 | ALA | -0.031 | 0 |
277 | ILE | -0.360 | 0.663 |
278 | ARG | -0.438 | -0.470 |
279 | ILE | -0.038 | 2.434 |
280 | PHE | -0.063 | 3.369 |
281 | CYN | -1.000 | -0.293 |
282 | GLY | -0.542 | 4.934 |
283 | GYS | -0.140 | -1.288 |
284 | GLN | -0.860 | -1.406 |
285 | PHE | -5.400 | 1.964 |
286 | ARG | -0.949 | 2.284 |
BFE: Binding free energy. |
Table 7
ΔΔG analysis of the BFE between PPARγ and β’-COP after amino acid sites mutation
Amino acid site | Amino acid | ΔΔG(kcal/mol) |
258 | ASP | 1.79 |
273 | LYS | 1 |
281 | CYN | 1.03 |
282 | GLY | 1.06–1.31 |
284 | GLN | 1.04 |
285 | PHE | 1.23–3.31 |
286 | ARG | 1.0-1.1 |
BFE: Binding free energy |
The mutation of ASP258, LYS273, GLN281, GLY282, PHE285 and ARG286 of PPARγ weaken the binding between PPARγ and β’-COP.
To elucidate the functions of the six amino acids in the β’-COP-PPARγ interaction, we employed bioinformatics to simulate the mutation of these six amino acids in the PPARγ protein to alanine. The stability of the system was assessed using RMSD, revealing initial stability from 0 to 30 ns, followed by significant fluctuations after 35 ns (Fig. 5A). These fluctuations were attributed to the impact of mutations on system stability, leading to partial dissociation and substantial RMSD fluctuations. Similar trends were observed for Rg and SASA, with Rg increasing significantly after 35 ns and SASA showing an increasing trend after 30 ns, indicating system instability post-mutation (Fig. 5B-C). Analysis of hydrogen bond distribution revealed a gradual decrease in the number of hydrogen bonds from 8–10 at the start of the simulation to approximately 1–2 by the end (Fig. 5D-H), aligning with the dissociation process observed in RMSD, Rg, and SASA.
We investigated the impact of a PPARγ mutation on BFE of PPARγ and β’-COP. The BFE was − 31.11 kcal/mol, primarily influenced by solvation free energy (ΔGSOLV, -42.63 kcal/mol) and van der Waals force (ΔVDWAALS, -26.42 kcal/mol), with Coulomb force (ΔEEL, 37.94 kcal/mol) as the main repulsive force (Table 8, Fig. 5I). The BFE initially increased within the first 20ns, followed by a subsequent decrease to a minimum of approximately − 8 kcal/mol at 42 ns. It then rose again, stabilizing at around − 32 kcal/mol (Fig. 5J), weaker than the BFE of the unmutated β’-COP-PPARγ complex (-38.55 kJ/mol). Key amino acids responsible for the bonding between the β’-COP-PPARγ complex were also altered and reduced (Fig. 5K-L). The mutation led to a significant decrease in the BFE between PPARγ and β’-COP, indicating a weakening of their interaction compared to the pre-mutation state.
Table 8
Analysis of the free energy of binding between PPARγ and β’-COP after mutation
Energy component | Average | SD | SEM |
VDWAALS | -26.42 | 10.2 | 0.14 |
EEL | 37.94 | 5.56 | 0.08 |
EGB | -38.14 | 5.01 | 0.07 |
ESURF | -4.49 | 1.82 | 0.03 |
GGAS | 11.52 | 13.43 | 0.19 |
GSOLV | -42.63 | 4.24 | 0.06 |
TOTAL | -31.11 | 11.78 | 0.17 |
SD: Standard Deviation; SEM: Standard Error of Mean; VDWAALS: van der Waals (vdW) Non-bonded Energy; EEL: Electrostatic (Coulombic) Energy between Charged Groups; EGB: Generalized Born Model for Solvation Energy; ESURF: Surface Area Term in Solvation Energy Calculation; GGAS: Gas Phase Energy Contribution to the Total Energy; GSOLV: Solvent Accessible Surface Area Energy; TOTAL: Total Energy Calculation. |
The mutation of ASP258, LYS273, GLN281, GLY282, PHE285 and ARG286 of PPARγ resulted in decreased expression of PPARγ in EVs.
To verify the participation of the six amino acids in β’-COP-mediated EVs transportation, we utilized overexpression lentivirus to introduce mutated PPARγ with these residues substituted with alanine to HTR8/SVneo trophoblast cell line. Concurrently, the fluorescent protein mCherry was incorporated upstream of the PPARγ protein, and a Histag label protein was attached downstream, resulting in a labeled PPARγ ASP258, LYS273, GLN281, GLY282, PHE285 and ARG286 mutated HTR8/SVneo trophoblast cell line (PPARγM-OE HTR8/SVneo) (Fig. 6A). The mCherry-PPARγ-histag overexpressed HTR8/SVneo (PPARγ-OE HTR8/SVneo) were used as a control group (Fig. 6A).
Confocal microscopy reveals the high transfection efficiency of the overexpressed letivirus in HTR8/SVneo trophoblast cells (Fig. 6B). Additionally, Western Blot and RT-qPCR analysis revealed a significantly overexpression of mCherry, PPARγ and histag both in PPARγM-OE and PPARγ-OE HTR8/SVneo cells (Fig. 6C-D). After mutation, the immunofluorescence assay revealed a significantly reduction in the co-localization of mCherry fluorescent protein with the EEA1 and CD63 in the PPARγM-OE HTR8/SVneo group (Fig. 6E-F). In contrast, the PPARγ-OE HTR8/SVneo group exhibited colocalization of mCherry protein with EEA1 and CD63 (Fig. 6E-F). Western Blot analysis also revealed a significantly reduction of PPARγ expression level in the EVs of PPARγM-OE HTR8/SVneo group compared to the PPARγ-OE HTR8/SVneo group (Fig. 6G). Taken together, these findings indicate that ASP258, LYS273, GLN281, GLY282, PHE285 and ARG286 of PPARγ are crucial amino acid sites for the β’-COP-mediated transportation of PPARγ into EVs.