Body Composition in Advanced Non-Small Cell Lung Cancer Treated With Immunotherapy.

Importance: The association between body composition (BC) and cancer outcomes is complex and incompletely understood. Previous research in non-small-cell lung cancer (NSCLC) has been limited to small, single-institution studies and yielded promising, albeit heterogeneous, results. Objectives: To evaluate the association of BC with oncologic outcomes in patients receiving immunotherapy for advanced or metastatic NSCLC. Design, Setting, and Participants: This comprehensive multicohort analysis included clinical data from cohorts receiving treatment at the Dana-Farber Brigham Cancer Center (DFBCC) who received immunotherapy given alone or in combination with chemotherapy and prospectively collected data from the phase 1/2 Study 1108 and the chemotherapy arm of the phase 3 MYSTIC trial. Baseline and follow-up computed tomography (CT) scans were collected and analyzed using deep neural networks for automatic L3 slice selection and body compartment segmentation (skeletal muscle [SM], subcutaneous adipose tissue [SAT], and visceral adipose tissue). Outcomes were compared based on baseline BC measures or their change at the first follow-up scan. The data were analyzed between July 2022 and April 2023. Main Outcomes and Measures: Hazard ratios (HRs) for the association of BC measurements with overall survival (OS) and progression-free survival (PFS). Results: A total of 1791 patients (878 women [49%]) with NSCLC were analyzed, of whom 487 (27.2%) received chemoimmunotherapy at DFBCC (DFBCC-CIO), 825 (46.1%) received ICI monotherapy at DFBCC (DFBCC-IO), 222 (12.4%) were treated with durvalumab monotherapy on Study 1108, and 257 (14.3%) were treated with chemotherapy on MYSTIC; median (IQR) ages were 65 (58-74), 66 (57-71), 65 (26-87), and 63 (30-84) years, respectively. A loss in SM mass, as indicated by a change in the L3 SM area, was associated with worse oncologic outcome across patient groups (HR, 0.59 [95% CI, 0.43-0.81] and 0.61 [95% CI, 0.47-0.79] for OS and PFS, respectively, in DFBCC-CIO; HR, 0.74 [95% CI, 0.60-0.91] for OS in DFBCC-IO; HR, 0.46 [95% CI, 0.33-0.64] and 0.47 [95% CI, 0.34-0.64] for OS and PFS, respectively, in Study 1108; HR, 0.76 [95% CI, 0.61-0.96] for PFS in the MYSTIC trial). This association was most prominent among male patients, with a nonsignificant association among female patients in the MYSTIC trial and DFBCC-CIO cohorts on Kaplan-Meier analysis. An increase of more than 5% in SAT density, as quantified by the average CT attenuation in Hounsfield units of the SAT compartment, was associated with poorer OS in 3 patient cohorts (HR, 0.61 [95% CI, 0.43-0.86] for DFBCC-CIO; HR, 0.62 [95% CI, 0.49-0.79] for DFBCC-IO; and HR, 0.56 [95% CI, 0.40-0.77] for Study 1108). The change in SAT density was also associated with PFS for DFBCC-CIO (HR, 0.73; 95% CI, 0.54-0.97). This was primarily observed in female patients on Kaplan-Meier analysis. Conclusions and Relevance: The results of this multicohort study suggest that loss in SM mass during systemic therapy for NSCLC is a marker of poor outcomes, especially in male patients. SAT density changes are also associated with prognosis, particularly in female patients. Automated CT-derived BC measurements should be considered in determining NSCLC prognosis.

Autoencoder-based multimodal prediction of non-small cell lung cancer survival.

The ability to accurately predict non-small cell lung cancer (NSCLC) patient survival is crucial for informing physician decision-making, and the increasing availability of multi-omics data offers the promise of enhancing prognosis predictions. We present a multimodal integration approach that leverages microRNA, mRNA, DNA methylation, long non-coding RNA (lncRNA) and clinical data to predict NSCLC survival and identify patient subtypes, utilizing denoising autoencoders for data compression and integration. Survival performance for patients with lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC) was compared across modality combinations and data integration methods. Using The Cancer Genome Atlas data, our results demonstrate that survival prediction models combining multiple modalities outperform single modality models. The highest performance was achieved with a combination of only two modalities, lncRNA and clinical, at concordance indices (C-indices) of 0.69 ± 0.03 for LUAD and 0.62 ± 0.03 for LUSC. Models utilizing all five modalities achieved mean C-indices of 0.67 ± 0.04 and 0.63 ± 0.02 for LUAD and LUSC, respectively, while the best individual modality performance reached C-indices of 0.64 ± 0.03 for LUAD and 0.59 ± 0.03 for LUSC. Analysis of biological differences revealed two distinct survival subtypes with over 900 differentially expressed transcripts.

Heritable transcriptional defects from aberrations of nuclear architecture.

Transcriptional heterogeneity due to plasticity of the epigenetic state of chromatin contributes to tumour evolution, metastasis and drug resistance1-3. However, the mechanisms that cause this epigenetic variation are incompletely understood. Here we identify micronuclei and chromosome bridges, aberrations in the nucleus common in cancer4,5, as sources of heritable transcriptional suppression. Using a combination of approaches, including long-term live-cell imaging and same-cell single-cell RNA sequencing (Look-Seq2), we identified reductions in gene expression in chromosomes from micronuclei. With heterogeneous penetrance, these changes in gene expression can be heritable even after the chromosome from the micronucleus has been re-incorporated into a normal daughter cell nucleus. Concomitantly, micronuclear chromosomes acquire aberrant epigenetic chromatin marks. These defects may persist as variably reduced chromatin accessibility and reduced gene expression after clonal expansion from single cells. Persistent transcriptional repression is strongly associated with, and may be explained by, markedly long-lived DNA damage. Epigenetic alterations in transcription may therefore be inherently coupled to chromosomal instability and aberrations in nuclear architecture.

Factors promoting nuclear envelope assembly independent of the canonical ESCRT pathway.

The nuclear envelope (NE) undergoes dynamic remodeling to maintain NE integrity, a process involving the inner nuclear membrane protein LEM2 recruiting CHMP7/Cmp7 and then ESCRT-III. However, prior work has hinted at CHMP7/ESCRT-independent mechanisms. To identify such mechanisms, we studied NE assembly in Schizosaccharomyces japonicus, a fission yeast that undergoes partial mitotic NE breakdown and reassembly. S. japonicus cells lacking Cmp7 have compromised NE sealing after mitosis but are viable. A genetic screen identified mutations that promote NE integrity in cmp7Δ cells. Unexpectedly, loss of Lem2 or its interacting partner Nur1 suppressed cmp7Δ defects. In the absence of Cmp7, Lem2 formed aggregates that appear to interfere with ESCRT-independent NE sealing. A gain-of-function mutation implicated a membrane and ESCRT-III regulator, Alx1, in this alternate pathway. Additional results suggest a potentially general role for unsaturated fatty acids in NE integrity. These findings establish the existence of mechanisms for NE sealing independent of the canonical ESCRT pathway.

Adipocytes sensitize melanoma cells to environmental TGF-ß cues by repressing the expression of miR-211.

Transforming growth factor-ß (TGF-ß) superfamily members are critical signals in tissue homeostasis and pathogenesis. Melanoma grows in the epidermis and invades the dermis before metastasizing. This disease progression is accompanied by increased sensitivity to microenvironmental TGF-ß. Here, we found that skin fat cells (adipocytes) promoted metastatic initiation by sensitizing melanoma cells to TGF-ß. Analysis of melanoma clinical samples revealed that adipocytes, usually located in the deeper hypodermis layer, were present in the upper dermis layer within proximity to in situ melanoma cells, an observation that correlated with disease aggressiveness. In a coculture system, adipocytes secreted the cytokines IL-6 and TNF-α, which induced a proliferative-to-invasive phenotypic switch in melanoma cells by repressing the expression of the microRNA miR-211. In a xenograft model, miR-211 exhibited a dual role in melanoma progression, promoting cell proliferation while inhibiting metastatic spread. Bioinformatics and molecular analyses indicated that miR-211 directly targeted and repressed the translation of TGFBR1 mRNA, which encodes the type I TGF-ß receptor. Hence, through this axis of cytokine-mediated repression of miR-211, adipocytes increased the abundance of the TGF-ß receptor in melanoma cells, thereby enhancing cellular responsiveness to TGF-ß ligands. The induction of TGF-ß signaling, in turn, resulted in a proliferative-to-invasive phenotypic switch in cultured melanoma cells. Pharmacological inhibition of TGF-ß prevented these effects. Our findings further reveal a molecular link between fat cells and metastatic progression in melanoma that might be therapeutically targeted in patients.

Breaching Self-Tolerance to Alu Duplex RNA Underlies MDA5-Mediated Inflammation.

Aberrant activation of innate immune receptors can cause a spectrum of immune disorders, such as Aicardi-Goutières syndrome (AGS). One such receptor is MDA5, a viral dsRNA sensor that induces antiviral immune response. Using a newly developed RNase-protection/RNA-seq approach, we demonstrate here that constitutive activation of MDA5 in AGS results from the loss of tolerance to cellular dsRNAs formed by Alu retroelements. While wild-type MDA5 cannot efficiently recognize Alu-dsRNAs because of its limited filament formation on imperfect duplexes, AGS variants of MDA5 display reduced sensitivity to duplex structural irregularities, assembling signaling-competent filaments on Alu-dsRNAs. Moreover, we identified an unexpected role of an RNA-rich cellular environment in suppressing aberrant MDA5 oligomerization, highlighting context dependence of self versus non-self discrimination. Overall, our work demonstrates that the increased efficiency of MDA5 in recognizing dsRNA comes at a cost of self-recognition and implicates a unique role of Alu-dsRNAs as virus-like elements that shape the primate immune system.

Cross-linking reveals laminin coiled-coil architecture.

Laminin, an ∼800-kDa heterotrimeric protein, is a major functional component of the extracellular matrix, contributing to tissue development and maintenance. The unique architecture of laminin is not currently amenable to determination at high resolution, as its flexible and narrow segments complicate both crystallization and single-particle reconstruction by electron microscopy. Therefore, we used cross-linking and MS, evaluated using computational methods, to address key questions regarding laminin quaternary structure. This approach was particularly well suited to the ∼750-Šcoiled coil that mediates trimer assembly, and our results support revision of the subunit order typically presented in laminin schematics. Furthermore, information on the subunit register in the coiled coil and cross-links to downstream domains provide insights into the self-assembly required for interaction with other extracellular matrix and cell surface proteins.

Codon-level information improves predictions of inter-residue contacts in proteins by correlated mutation analysis.

Methods for analysing correlated mutations in proteins are becoming an increasingly powerful tool for predicting contacts within and between proteins. Nevertheless, limitations remain due to the requirement for large multiple sequence alignments (MSA) and the fact that, in general, only the relatively small number of top-ranking predictions are reliable. To date, methods for analysing correlated mutations have relied exclusively on amino acid MSAs as inputs. Here, we describe a new approach for analysing correlated mutations that is based on combined analysis of amino acid and codon MSAs. We show that a direct contact is more likely to be present when the correlation between the positions is strong at the amino acid level but weak at the codon level. The performance of different methods for analysing correlated mutations in predicting contacts is shown to be enhanced significantly when amino acid and codon data are combined.

N-terminal domains in two-domain proteins are biased to be shorter and predicted to fold faster than their C-terminal counterparts.

Computational analysis of proteomes in all kingdoms of life reveals a strong tendency for N-terminal domains in two-domain proteins to have shorter sequences than their neighboring C-terminal domains. Given that folding rates are affected by chain length, we asked whether the tendency for N-terminal domains to be shorter than their neighboring C-terminal domains reflects selection for faster-folding N-terminal domains. Calculations of absolute contact order, another predictor of folding rate, provide additional evidence that N-terminal domains tend to fold faster than their neighboring C-terminal domains. A possible explanation for this bias, which is more pronounced in prokaryotes than in eukaryotes, is that faster folding of N-terminal domains reduces the risk for protein aggregation during folding by preventing formation of nonnative interdomain interactions. This explanation is supported by our finding that two-domain proteins with a shorter N-terminal domain are much more abundant than those with a shorter C-terminal domain.

Different mechanistic requirements for prokaryotic and eukaryotic chaperonins: a lattice study.

MOTIVATION: The folding of many proteins in vivo and in vitro is assisted by molecular chaperones. A well-characterized molecular chaperone system is the chaperonin GroEL/GroES from Escherichia coli which has a homolog found in the eukaryotic cytosol called CCT. All chaperonins have a ring structure with a cavity in which the substrate protein folds. An interesting difference between prokaryotic and eukaryotic chaperonins is in the nature of the ATP-mediated conformational changes that their ring structures undergo during their reaction cycle. Prokaryotic chaperonins are known to exhibit a highly cooperative concerted change of their cavity surface while in eukaryotic chaperonins the change is sequential. Approximately 70% of proteins in eukaryotic cells are multi-domain whereas in prokaryotes single-domain proteins are more common. Thus, it was suggested that the different modes of action of prokaryotic and eukaryotic chaperonins can be explained by the need of eukaryotic chaperonins to facilitate folding of multi-domain proteins. RESULTS: Using a 2D square lattice model, we generated two large populations of single-domain and double-domain substrate proteins. Chaperonins were modeled as static structures with a cavity wall with which the substrate protein interacts. We simulated both concerted and sequential changes of the cavity surfaces and demonstrated that folding of single-domain proteins benefits from concerted but not sequential changes whereas double-domain proteins benefit also from sequential changes. Thus, our results support the suggestion that the different modes of allosteric switching of prokaryotic and eukaryotic chaperonin rings have functional implications as it enables eukaryotic chaperonins to better assist multi-domain protein folding.

A tale of two tails: why are terminal residues of proteins exposed?

MOTIVATION: It is widely known that terminal residues of proteins (i.e. the N- and C-termini) are predominantly located on the surface of proteins and exposed to the solvent. However, there is no good explanation as to the forces driving this phenomenon. The common explanation that terminal residues are charged, and charged residues prefer to be on the surface, cannot explain the magnitude of the phenomenon. Here, we survey a large number of proteins from the PDB in order to explore, quantitatively, this phenomenon, and then we use a lattice model to study the mechanisms involved. RESULTS: The location of the termini was examined for 425 small monomeric proteins (50-200 amino acids) and it was found that the average solvent accessibility of termini residues is 87.1% compared with 49.2% of charged residues and 35.9% of all residues. Using a cutoff of 50% of the maximal possible exposure, 80.3% of the N-terminal and 86.1% of the C-terminal residues are exposed compared to 32% for all residues. In addition, terminal residues are much more distant from the center of mass of their proteins than other residues. Using a 2D lattice, a large population of model proteins was studied on three levels: structural selection of compact structures, thermodynamic selection of conformations with a pronounced energy gap and kinetic selection of fast folding proteins using Monte-Carlo simulations. Progressively, each selection raises the proportion of proteins with termini on the surface, resulting in similar proportions to those observed for real proteins.