Computational chemistry for complex systems : open-shell molecules to conjugated organic materials

This thesis focuses on two different, but equally challenging, areas of computational chemistry: transition metal organic molecule interactions and parameterisation of organic conjugated polymers for molecular dynamics simulations. The metal-binding properties are important for understanding of biomolecular action of type 2 diabetes drug and development of novel protocols for redox calculations of copper systems. In this area the challenge is mainly related to the complex electronic structure of the open-shell transition metals. The main challenges for the parameterisation of conjugated polymers are due to the size of the studied systems, their conjugated nature and inclusion of environment. Metal-binding properties as well as electronic structures of copper complexes of type 2 diabetes drug metformin (Metf) and other similar, but often inactive, compounds were examined using DFT method. It was found that for neutral compounds it is not possible to explain the differences in their biological effects solely by examining the copper-binding properties. Further, the proposed mechanism potentially explaining the difference in the biomolecular mode of action involves a possible deprotonation of biguanide and Metf compounds under higher mitochondrial pH which would lead to formation of more stable copper complexes and potentially affecting the mitochondrial copper homeostasis. In addition, redox properties of copper-biguanide complexes could interfere with the sensitive redox chemistry or interact with important metalloproteins in the mitochondria. Understanding the copper-binding properties is also important for a systematic development and testing of computational protocols for calculations of reduction potentials of copper complexes. Copper macrocyclic complexes previously used as model systems for redox-active metalloenzymes and for which experimentally determined redox potentials are available were used as model systems. First adequacy of using single reference methods such as DFT was examined for these systems and then various DFT functionals and basis sets were tested in order to develop accurate redox potential protocol. It was shown that good relative cor-relations were obtained for several functionals while the best absolute agreement was obtained with either the M06/cc-pVTZ functional with the SMD or either M06L or TPSSTPSS functional with cc-pVTZ basis set and the PCM solvation model. Organic conjugated polymers have a great potential due to their application in organic optoelectronics. Various wavefunction and DFT methods are utilized in order to systematically develop parameterisation scheme that can be used to derive selected force-field parameters such as torsional potentials between monomer units that are critical for these systems and partial charges. Moreover, critical points of such a parameterisation are addressed in order to obtain accurate MD simulations that could provide valuable insight into material morphology and conformation that affect their optical properties and conductivity. It was shown that a two step approach of geometry optimisation with CAM-B3LYP/631G* and single point (SP) energy scan with CAM-B3LYP/cc-pVTZ is able to yield accurate dihedral potentials in agreement with the potentials calculated using higher level methods such as MP2 and CBS limit CCSD(T). Further, investigating partial charge distribution for increasing backbone length of fluorene and thiophene it has been found that it is possible to obtain a three residue model of converged charge distributions using the RESP scheme. The three partial charge residues can be then used to build and simulate much longer polymers without the need to re-parametrize charge distributions. In the case of side-chains, it was found that it is not possible to obtain converged charge sets for sidechain lengths of up to 10 carbons due to the strong asymmetry between the side-chain ends. Initial validation of derived force-field parameters performed by simulations of 32mers of fluorene with octyl side-chains (PF8) and thiophene with hexyl side-chains (P3HT) in chloroform and calculation of persistence lengths and end-to-end lengths showed close correspondence to experimentally obtained values

Sprouty2 loss‐induced IL6 drives castration‐resistant prostate cancer through scavenger receptor B1

Metastatic castration‐resistant prostate cancer (mCRPC) is a lethal form of treatment‐resistant prostate cancer and poses significant therapeutic challenges. Deregulated receptor tyrosine kinase (RTK) signalling mediated by loss of tumour suppressor Sprouty2 (SPRY2) is associated with treatment resistance. Using pre‐clinical human and murine mCRPC models, we show that SPRY2 deficiency leads to an androgen self‐sufficient form of CRPC. Mechanistically, HER2‐IL6 signalling axis enhances the expression of androgen biosynthetic enzyme HSD3B1 and increases SRB1‐mediated cholesterol uptake in SPRY2‐deficient tumours. Systemically, IL6 elevated the levels of circulating cholesterol by inducing host adipose lipolysis and hepatic cholesterol biosynthesis. SPRY2‐deficient CRPC is dependent on cholesterol bioavailability and SRB1‐mediated tumoral cholesterol uptake for androgen biosynthesis. Importantly, treatment with ITX5061, a clinically safe SRB1 antagonist, decreased treatment resistance. Our results indicate that cholesterol transport blockade may be effective against SPRY2‐deficient CRPC

Sleeping Beauty screen reveals Pparg activation in metastatic prostate cancer

Prostate cancer (CaP) is the most common adult male cancer in the developed world. The paucity of biomarkers to predict prostate tumor biology makes it important to identify key pathways that confer poor prognosis and guide potential targeted therapy. Using a murine forward mutagenesis screen in a Pten-null background, we identified peroxisome proliferator-activated receptor gamma (Pparg), encoding a ligand-activated transcription factor, as a promoter of metastatic CaP through activation of lipid signaling pathways, including up-regulation of lipid synthesis enzymes [fatty acid synthase (FASN), acetyl-CoA carboxylase (ACC), ATP citrate lyase (ACLY)]. Importantly, inhibition of PPARG suppressed tumor growth in vivo, with down-regulation of the lipid synthesis program. We show that elevated levels of PPARG strongly correlate with elevation of FASN in human CaP and that high levels of PPARG/FASN and PI3K/pAKT pathway activation confer a poor prognosis. These data suggest that CaP patients could be stratified in terms of PPARG/FASN and PTEN levels to identify patients with aggressive CaP who may respond favorably to PPARG/FASN inhibition

Germline biallelic mutation affecting the transcription factor Helios causes pleiotropic defects of immunity

Helios, a member of the Ikaros family of transcription factors, is predominantly expressed in developing thymocytes, activated T cells, and regulatory T cells (T-regs). Studies in mice have emphasized its role in maintenance of T-reg immunosuppressive functions by stabilizing Foxp3 expression and silencing the Il2 locus. However, its contribution to human immune homeostasis and the precise mechanisms by which Helios regulates other T cell subsets remain unresolved. Here, we investigated a patient with recurrent respiratory infections and hypogammaglobulinemia and identified a germline homozygous missense mutation in IKZF2 encoding Helios (p.Ile325Val). We found that Helios(I325V) retains DNA binding and dimerization properties but loses interaction with several partners, including epigenetic remodelers. Whereas patient T-regs showed increased IL-2 production, patient conventional T cells had decreased accessibility of the IL2 locus and consequently reduced IL-2 production. Reduced chromatin accessibility was not exclusive to the IL2 locus but involved a variety of genes associated with T cell activation. Single-cell RNA sequencing of peripheral blood mononuclear cells revealed gene expression signatures indicative of a shift toward a proinflammatory, effector-like status in patient CD8(+) T cells. Moreover, patient CD4(+) T cells exhibited a pronounced defect in proliferation with delayed expression of surface checkpoint inhibitors, suggesting an impaired onset of the T cell activation program. Collectively, we identified a previously uncharacterized, germline-encoded inborn error of immunity and uncovered a cell-specific defect in Helios-dependent epigenetic regulation. Binding of Helios with specific partners mediates this regulation, which is ultimately necessary for the transcriptional programs that enable T cell homeostasis in health and disease.Peer reviewe

18F-Fluciclovine PET metabolic imaging reveals prostate cancer tumour heterogeneity associated with disease resistance to androgen deprivation therapy

Background: Prostate cancer is highly prevalent worldwide. Androgen deprivation therapy (ADT) remains the treatment of choice for incurable prostate cancer, but majority of patients develop disease recurrence following ADT. There is therefore an urgent need for early detection of treatment resistance. Methods: Isogenic androgen-responsive (CWR22Res) and castration-resistant (22Rv1) human prostate cancer cells were implanted into the anterior lobes of the prostate in CD-1 Nu mice to generate prostate orthografts. Castrated mice bearing CWR22Res and 22Rv1 orthografts mimic clinical prostate cancer following acute and chronic ADT, respectively. 18F-Fluciclovine (1-amino-3-fluorocyclobutane-1-carboxylic acid) with a radiochemical purity of > 99% was produced on a FASTlab synthesiser. Ki67 staining in endpoint orthografts was studied. Western blot, quantitative RT-PCR and next-generation sequencing transcriptomic analyses were performed to assess the expression levels of amino acid transporters (including LAT1 and ASCT2, which have been implicated for Fluciclovine uptake). Longitudinal metabolic imaging with 18F-Fluciclovine-based positron emission tomography (PET) was performed to study tumour response following acute and chronic ADT. Results: Both immunohistochemistry analysis of endpoint prostate tumours and longitudinal 18F-Fluciclovine imaging revealed tumour heterogeneity, particularly following ADT, with in vivo 18F-Fluciclovine uptake correlating to viable cancer cells in both androgen-proficient and castrated environment. Highlighting tumour subpopulation following ADT, both SUVpeak and coefficient of variation (CoV) values of 18F-Fluciclovine uptake are consistent with tumour heterogeneity revealed by immunohistochemistry. We studied the expression of amino acid transporters (AATs) for 18F-Fluciclovine, namely LAT1 (SLC7A5 and SLC3A2) and ASCT2 (SLC1A5). SLC7A5 and SLC3A2 were expressed at relatively high levels in 22Rv1 castration-resistant orthografts following chronic ADT (modelling clinical castration-resistant disease), while SLC1A5 was preferentially expression in CWR22Res tumours following acute ADT. Additional AATs such as SLC43A2 (LAT4) were shown to be upregulated following chronic ADT by transcriptomic analysis; their role in Fluciclovine uptake warrants investigation. Conclusion: We studied in vivo 18F-Fluciclovine uptake in human prostate cancer orthograft models following acute and chronic ADT. 18F-Fluciclovine uptakes highlight tumour heterogeneity that may explain castration resistance and can be exploited as a clinical biomarker

BRF1 accelerates prostate tumourigenesis and perturbs immune infiltration

BRF1 is a rate-limiting factor for RNA Polymerase III-mediated transcription and is elevated in numerous cancers. Here, we report that elevated levels of BRF1 associate with poor prognosis in human prostate cancer. In vitro studies in human prostate cancer cell lines demonstrated that transient overexpression of BRF1 increased cell proliferation whereas the transient downregulation of BRF1 reduced proliferation and mediated cell cycle arrest. Consistent with our clinical observations, BRF1 overexpression in a Pten-deficient mouse (Pten BRF1 ) prostate cancer model accelerated prostate carcinogenesis and shortened survival. In Pten BRF1 tumours, immune and inflammatory processes were altered, with reduced tumoral infiltration of neutrophils and CD4 positive T cells, which can be explained by decreased levels of complement factor D (CFD) and C7 components of the complement cascade, an innate immune pathway that influences the adaptive immune response. We tested if the secretome was involved in BRF1-driven tumorigenesis. Unbiased proteomic analysis on BRF1-overexpresing PC3 cells confirmed reduced levels of CFD in the secretome, implicating the complement system in prostate carcinogenesis. We further identify that expression of C7 significantly correlates with expression of CD4 and has the potential to alter clinical outcome in human prostate cancer, where low levels of C7 associate with poorer prognosis

Landscape of Bone Marrow Metastasis in Human Neuroblastoma Unraveled by Transcriptomics and Deep Multiplex Imaging

Bone marrow metastasis frequently occurs in patients with solid cancers and most often leads to poor outcome. Yet, the composition of bone marrow metastases, including tumor and surrounding cells, has so far not been characterized. Herein, we aimed to investigate the diversity of tumor and surrounding cells, i.e., the microenvironment, in bone marrow metastases, using the childhood tumor neuroblastoma as a model. To this end, we screened genome-wide datasets to define a panel of cell-specific markers for multiplex microscopy of metastatic bone marrow samples, and developed DeepFLEX, a computational pipeline for subsequent image analysis. Thereby, we identified 35,000 single cells covering metastasized tumor cells, and various types of developing immune and bone marrow cells. In parallel, we analyzed the transcriptome, i.e., all genes that are expressed as mRNA, of 38 patients with and without bone marrow metastasis. We found vast tumor cell diversity and identified a marker protein, FAIM2, which can help to identify a broader range of tumor cell variants. In addition we showed that tumor cell metastasis in the bone marrow is associated with an immune response resembling inflammation, and the presence of cells that can repress an immune attack against cancer cells. Our study suggests that metastatic tumor cells are shaping the bone marrow microenvironment and builds the basis to further investigate its clinical relevance

Sprouty2 loss‐induced IL6 drives castration‐resistant prostate cancer through scavenger receptor B1

Metastatic castration‐resistant prostate cancer (mCRPC) is a lethal form of treatment‐resistant prostate cancer and poses significant therapeutic challenges. Deregulated receptor tyrosine kinase (RTK) signalling mediated by loss of tumour suppressor Sprouty2 (SPRY2) is associated with treatment resistance. Using pre‐clinical human and murine mCRPC models, we show that SPRY2 deficiency leads to an androgen self‐sufficient form of CRPC. Mechanistically, HER2‐IL6 signalling axis enhances the expression of androgen biosynthetic enzyme HSD3B1 and increases SRB1‐mediated cholesterol uptake in SPRY2‐deficient tumours. Systemically, IL6 elevated the levels of circulating cholesterol by inducing host adipose lipolysis and hepatic cholesterol biosynthesis. SPRY2‐deficient CRPC is dependent on cholesterol bioavailability and SRB1‐mediated tumoral cholesterol uptake for androgen biosynthesis. Importantly, treatment with ITX5061, a clinically safe SRB1 antagonist, decreased treatment resistance. Our results indicate that cholesterol transport blockade may be effective against SPRY2‐deficient CRPC

2,4-dienoyl-CoA reductase regulates lipid homeostasis in treatment-resistant prostate cancer

International audienceDespite the clinical success of Androgen Receptor (AR)-targeted therapies, reactivation of AR signalling remains the main driver of castration-resistant prostate cancer (CRPC) progression. In this study, we perform a comprehensive unbiased characterisation of LNCaP cells chronically exposed to multiple AR inhibitors (ARI). Combined proteomics and metabolomics analyses implicate an acquired metabolic phenotype common in ARI-resistant cells and associated with perturbed glucose and lipid metabolism. To exploit this phenotype, we delineate a subset of proteins consistently associated with ARI resistance and highlight mitochondrial 2,4-dienoyl-CoA reductase (DECR1), an auxiliary enzyme of beta-oxidation, as a clinically relevant biomarker for CRPC. Mechanistically, DECR1 participates in redox homeostasis by controlling the balance between saturated and unsaturated phospholipids. DECR1 knockout induces ER stress and sensitises CRPC cells to ferroptosis. In vivo, DECR1 deletion impairs lipid metabolism and reduces CRPC tumour growth, emphasizing the importance of DECR1 in the development of treatment resistance