SEARCH · Search PubMed

Results for “Mutation Signature”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,045 records · Page 58Linked to original sources

PubMed2005

The THAP domain of THAP1 is a large C2CH module with zinc-dependent sequence-specific DNA-binding activity.

We have recently described an evolutionarily conserved protein motif, designated the THAP domain, which defines a previously uncharacterized family of cellular factors (THAP proteins). The THAP domain exhibits similarities to the site-specific DNA-binding domain of Drosophila P element transposase, including a putative metal-coordinating C2CH signature (CX(2-4)CX(35-53)CX(2)H). In this article, we report a comprehensive list of approximately 100 distinct THAP proteins in model animal organisms, including human nuclear proapoptotic factors THAP1 and DAP4/THAP0, transcriptional repressor THAP7, zebrafish orthologue of cell cycle regulator E2F6, and Caenorhabditis elegans chromatin-associated protein HIM-17 and cell-cycle regulators LIN-36 and LIN-15B. In addition, we demonstrate the biochemical function of the THAP domain as a zinc-dependent sequence-specific DNA-binding domain belonging to the zinc-finger superfamily. In vitro binding-site selection allowed us to identify an 11-nucleotide consensus DNA-binding sequence specifically recognized by the THAP domain of human THAP1. Mutations of single nucleotide positions in this sequence abrogated THAP-domain binding. Experiments with the zinc chelator 1,10-o-phenanthroline revealed that the THAP domain is a zinc-dependent DNA-binding domain. Site-directed mutagenesis of single cysteine or histidine residues supported a role for the C2CH motif in zinc coordination and DNA-binding activity. The four other conserved residues (P, W, F, and P), which define the THAP consensus sequence, were also found to be required for DNA binding. Together with previous genetic data obtained in C. elegans, our results suggest that cellular THAP proteins may function as zinc-dependent sequence-specific DNA-binding factors with roles in proliferation, apoptosis, cell cycle, chromosome segregation, chromatin modification, and transcriptional regulation.

Amino Acid Motifs↗

PubMed2026

Distinct molecular profiles of indeterminate and malignant thyroid nodules in patients under 21 years of age.

Although uncommon, thyroid nodules (TN) in pediatric and young adult patients carry higher malignancy risk and often present with a high burden of metastatic disease than adults. The molecular features underlying this distinct clinical behavior remain unclear. We analyzed Afirma Genomic Sequencing Classifier (GSC) data from 283,621 TN, comparing patients <21 and ≥21 years. Cytology (Bethesda), GSC benign (B) vs suspicious (S) calls, and Afirma Xpression Atlas (XA) variant/fusion profiles were evaluated in GSC-S and Bethesda V/VI samples. Genome-wide expression was used to derive pathway signatures and thyroid cancer-related scores: BRAF-RAS score (BRS), ERK, follicular and epithelial-to-mesenchymal transition (FMT, EMT) and thyroid differentiation scores (TDS). Among 2,397 patients <21 (median age 18.9; 81.4% female) and 281,224 adults ≥21 (median age 59.8; 77.1% female), <21 samples showed more Bethesda V/VI cytology (14.5% vs 5.0%; p<0.0001) and a lower GSC-B rate (43.5% vs 68.8%; p<0.0001). In GSC-S samples, total variant detection was higher in <21 (45.3% vs 37.4%), with enriched BRAF p.V600E, TSHR, and DICER1 variants, while HRAS variants were more common in adults (all p<0.01). Gene fusions involving RET, NTRK3 and ALK were enriched in <21 (14.5% vs 5.5%; p<0.0001). TERT promoter mutations were absent in <21 yrs GSC-S and Bethesda V/VI samples (vs 4.2% and 9.3% in adults). GSC-S <21 showed cell-cycle pathway enrichment. RET/NTRK/ALK-positive <21 demonstrated enrichment of angiogenesis and EMT pathways, higher ERK/EMT/FMT scores, and lower BRS/TDS scores vs genotyped-matched adults. These molecular differences provide mechanistic insight into the more invasive phenotype in pediatric and young adult TN.

BRAF↗

PubMed2005

Initial coupling of binding to gating mediated by conserved residues in the muscle nicotinic receptor.

We examined functional consequences of intrasubunit contacts in the nicotinic receptor alpha subunit using single channel kinetic analysis, site-directed mutagenesis, and structural modeling. At the periphery of the ACh binding site, our structural model shows that side chains of the conserved residues alphaK145, alphaD200, and alphaY190 converge to form putative electrostatic interactions. Structurally conservative mutations of each residue profoundly impair gating of the receptor channel, primarily by slowing the rate of channel opening. The combined mutations alphaD200N and alphaK145Q impair channel gating to the same extent as either single mutation, while alphaK145E counteracts the impaired gating due to alphaD200K, further suggesting electrostatic interaction between these residues. Interpreted in light of the crystal structure of acetylcholine binding protein (AChBP) with bound carbamylcholine (CCh), the results suggest in the absence of ACh, alphaK145 and alphaD200 form a salt bridge associated with the closed state of the channel. When ACh binds, alphaY190 moves toward the center of the binding cleft to stabilize the agonist, and its aromatic hydroxyl group approaches alphaK145, which in turn loosens its contact with alphaD200. The positional changes of alphaK145 and alphaD200 are proposed to initiate the cascade of perturbations that opens the receptor channel: the first perturbation is of beta-strand 7, which harbors alphaK145 and is part of the signature Cys-loop, and the second is of beta-strand 10, which harbors alphaD200 and connects to the M1 domain. Thus, interplay between these three conserved residues relays the initial conformational change from the ACh binding site toward the ion channel.

Acetylcholine↗

PubMed2006

How reliable are empirical genomic scans for selective sweeps?

The beneficial substitution of an allele shapes patterns of genetic variation at linked sites. Thus, in principle, adaptations can be mapped by looking for the signature of directional selection in polymorphism data. In practice, such efforts are hampered by the need for an accurate characterization of the demographic history of the species and of the effects of positive selection. In an attempt to circumvent these difficulties, researchers are increasingly taking a purely empirical approach, in which a large number of genomic regions are ordered by summaries of the polymorphism data, and loci with extreme values are considered to be likely targets of positive selection. We evaluated the reliability of the "empirical" approach, focusing on applications to human data and to maize. To do so, we considered a coalescent model of directional selection in a sensible demographic setting, allowing for selection on standing variation as well as on a new mutation. Our simulations suggest that while empirical approaches will identify several interesting candidates, they will also miss many--in some cases, most--loci of interest. The extent of the trade-off depends on the mode of positive selection and the demographic history of the population. Specifically, the false-discovery rate is higher when directional selection involves a recessive rather than a co-dominant allele, when it acts on a previously neutral rather than a new allele, and when the population has experienced a population bottleneck rather than maintained a constant size. One implication of these results is that, insofar as attributes of the beneficial mutation (e.g., the dominance coefficient) affect the power to detect targets of selection, genomic scans will yield an unrepresentative subset of loci that contribute to adaptations.

Africa South of the Sahara↗

PubMed1981

Central-pair microtubular complex of Chlamydomonas flagella: polypeptide composition as revealed by analysis of mutants.

Four mutants of Chlamydomonas reinhardtii representing independent gene loci have been shown to lack totally (pf-18, pf-19, and pf-15) or nearly totally (pf-20) the central microtubular pair complex in isolated axonemal preparations. Analysis of 35S-labeled axonemal proteins, using two methods of electrophoresis, reveals that all four mutants lack or are markedly deficient in 18 polypeptides, ranging in molecular weight from 360,000 to 20,000, that are regularly present in wild-type axonemes. Analyses of axonemal proteins labeled by cellular growth on 32P-labeled medium indicates that a subset of 8 of the 18 polypeptides are phosphorylated. Mutant and wild-type axonemes and flagella have been analyzed for their content of tubulin subunits using a high resolution two-dimensional electrophoresis system combined with agarose gel overlays containing either anti-alpha or anti-beta tubulin sera prepared from Chlamydomonas tubulins. The immunoprecipitates identify two major alpha tubulins, a major beta tubulin, and a minor component which is also precipitated by the anti-beta serum. None of these tubulins shows a specific defect in mutant axonemes, nor do the tubulin polypeptides show altered two-dimensional map positions in the mutant flagella. The 18 polypeptides provide a useful signature for identifying other mutants affecting the central-pair microtubular complex. Such mutants could be useful in defining the structural or functional role of these polypeptides in the central microtubules. Efforts to obtain additional central-pair mutants based on the motility phenotype of the four mutants analyzed here have yielded mutants which are allelic to three of the four mutants.

Chlamydomonas↗

PubMed2002

The signature of positive selection at randomly chosen loci.

In Drosophila and humans, there are accumulating examples of loci with a significant excess of high-frequency-derived alleles or high levels of linkage disequilibrium, relative to a neutral model of a random-mating population of constant size. These are features expected after a recent selective sweep. Their prevalence suggests that positive directional selection may be widespread in both species. However, as I show here, these features do not persist long after the sweep ends: The high-frequency alleles drift to fixation and no longer contribute to polymorphism, while linkage disequilibrium is broken down by recombination. As a result, loci chosen without independent evidence of recent selection are not expected to exhibit either of these features, even if they have been affected by numerous sweeps in their genealogical history. How then can we explain the patterns in the data? One possibility is population structure, with unequal sampling from different subpopulations. Alternatively, positive selection may not operate as is commonly modeled. In particular, the rate of fixation of advantageous mutations may have increased in the recent past.

Adaptation, Biological↗

PubMed2026

Human Monocytic Models Reveal Genotype-Dependent Inflammatory Programs in VEXAS Syndrome.

OBJECTIVES: VEXAS syndrome is a severe X-linked autoinflammatory disorder caused by somatic mutations in ubiquitin-like modifier activating enzyme 1 (UBA1), with clinical outcomes that vary by UBA1 genotype. We aimed to elucidate genotype-specific inflammatory programs and identify potential therapeutic targets. METHODS: We conducted longitudinal deep phenotyping, including whole-blood RNA sequencing (RNA-seq) and clinical activity assessment. Peripheral blood samples were analyzed by single-cell RNA-seq. Human monocytic cell lines harboring each major UBA1 mutation (p.Met41Val, p.Met41Thr, or p.Met41Leu) were generated and subjected to transcriptomic and functional analyses. RESULTS: Thirteen patients with VEXAS syndrome contributed a total of 79 RNA-seq samples. Among genes upregulated in VEXAS syndrome, RNASE1 showed the strongest correlation with longitudinal disease activity (r = 0.70, FDR < 0.05) and was upregulated in patients' monocytes. In UBA1-mutant monocytic cell lines, genotype-dependent ubiquitination defects were observed in a graded manner (p.Met41Val > p.Met41Thr > p.Met41Leu), even in the absence of exogenous stimuli. These defects were accompanied by unfolded protein response activation, increased pro-inflammatory cytokine production, progressive cell death, and RNASE1 upregulation, all following the same graded pattern, recapitulating patient genotype-phenotype associations. Transcriptomic analyses demonstrated enrichment of pro-inflammatory, interferon, and necroptosis signatures in more severe genotypes. Notably, inhibition of receptor-interacting protein kinase 3 (RIPK3) markedly attenuated all pathological features, including RNASE1 upregulation. CONCLUSIONS: Our UBA1-mutant monocytic cell-line models, representing three distinct genotypes, recapitulate genotype-dependent inflammatory phenotypes that can be modulated by RIPK3 inhibition, providing a translational platform for mechanistic investigation and precision therapy development in VEXAS syndrome.

Journal Article↗

PubMed1992

RNA binding assays for Tat-derived peptides: implications for specificity.

RNA recognition by the HIV Tat protein is mediated in part by an arginine- and lysine-rich basic subdomain implicated as a signature element in proteins that bind RNA. Relative RNA binding affinities for a 14-residue peptide derived from Tat that spans the basic region are determined using a competition protocol. Binding specificity is compared with complexation by a 38-residue model for the RNA binding domain of Tat using the same approach. Binding strength for the minimal (14 residue) peptide is correlated with that for the longer peptide: both peptides recognize a short, bulged duplex. However, the shorter peptide dissociates more rapidly from the wild-type site and discriminates less well between nonspecific (double-stranded RNA) and specific sites. Relative dissociation constants for 38-residue peptide determined from direct partition and competition assays differ; the former assay consistently predicts stronger discrimination against RNAs with mutations in the stems flanking the bulge. Differences between the two assays are reconciled in terms of contributions from labile binding which is unstable to native gel electrophoresis. Kinetic stability may constitute a major specificity determinant for basic subdomain-mediated recognition of RNA.

Base Sequence↗

PubMed2000

Catalysis of DNA cleavage and nucleoside triphosphate synthesis by NM23-H2/NDP kinase share an active site that implies a DNA repair function.

NM23/NDP kinases play an important role in development and cancer but their biological function is unknown, despite an intriguing collection of biochemical properties including nucleoside-diphosphate kinase (NDP kinase), DNA binding and transcription, a mutator function, and cleavage of unusually structured DNA by means of a covalent enzyme-DNA complex. To assess the role of the nuclease in human NM23-H2, we sought to identify the amino acid responsible for covalent catalysis. By sequencing a DNA-linked peptide and by site-directed mutagenesis, we identified lysine-12, a phylogenetically conserved residue, as the amino acid forming the covalent complex with DNA. In particular, the epsilon-amino group acts as the critical nucleophile, because substitution with glutamine but not arginine completely abrogated covalent adduct formation and DNA cleavage, whereas the DNA-binding properties remained intact. These findings and chemical modification data suggest that phosphodiester-bond cleavage occurs by a DNA glycosylase/lyase-like mechanism known as the signature of base excision DNA repair nucleases. Involvement of NM23/NDP kinase in a DNA repair pathway would be consistent with its role in normal and tumor cell development. Additionally, lysine-12, which is known in the x-ray crystallographic structure to lie in the catalytic pocket involved in the NDP kinase phosphorylation reaction, was found essential also for the NDP kinase activity of NM23-H2, suggesting that the two catalytic activities of NM23-H2 are fundamentally connected.

Binding Sites↗

PubMed1999

SST versus EST in gene recognition.

The expressed sequence tag (EST) data provide a powerful tool for identification of transcribed DNA sequences. However, as EST are relatively short, many exons are poorly covered by EST, thus reducing the utility of EST data. Recently, signature sequence tag (SST) fingerprints were proposed as an alternative to EST fingerprints. Given a fingerprint set of probes, SST of a clone is a subset of probes from the fingerprint set that hybridize with the clone. We demonstrate that besides being a powerful technique for screening cDNA libraries, SST technology provides for very accurate gene predictions. Even with a small fingerprint set (600-800 probes), SST-based gene recognition outperforms many conventional and EST-based methods. The increase in the size of the fingerprint set to 1500 probes provides almost perfect gene recognition. Even more importantly, SST-based gene predictions miss very few exons and, therefore, provide an opportunity to bypass the cDNA sequencing step on the way from finished genomic sequence to mutation detection in gene-hunting projects. Because SST data can be obtained in a highly parallel and inexpensive way, SST technology has a potential of complementing EST technology for gene hunting.

Computational Biology↗

PubMed2026

Virtual Tumors Enable Prediction of Personalized Therapeutic Combinations for Non-Small Cell Lung Cancer.

UNLABELLED: The disease burden from non-small cell lung cancer (NSCLC) adenocarcinoma is substantial, with a million new cases diagnosed globally each year and a 5-year survival rate of less than 20%. The lack of therapeutic options personalized to individual patients leads to high variation in survival. The combination of patient stratification with personalized treatment has the potential to improve outcomes; however, the variation in mutations found in patients with NSCLC adenocarcinoma makes experimentally determining treatment combinations time-consuming and expensive. In this study, we developed an interpretable mechanistic model to decipher complex signaling interplay and guide personalized therapy in NSCLC adenocarcinoma. This "virtual tumor" model encompassed key tumor-intrinsic oncogenic signaling pathways for efficiently predicting rational drug-drug and drug-radiotherapy combination therapies in NSCLC. Diverse genetic profiles were simulated for testing more than 10,000 therapeutic strategies to identify optimal approaches to overcome resistance mechanisms specific to genetic profiles and p53 status. The virtual tumor model reproduced drug additivity screens, predicted radiosensitizing genes validated in a CRISPR screen, and identified 53BP1 as a potential drug target that improved the therapeutic window during radiotherapy. A 19-gene signature derived from the virtual tumor framework stratified patients most likely to benefit from radiotherapy, which was validated using The Cancer Genome Atlas (TCGA) data. These results show the utility of virtual tumors to predict effective therapeutic combinations and present a computational resource for large-scale screening of personalized therapies to guide clinical decision-making in patients with NSCLC. SIGNIFICANCE: A computational framework that simulates thousands of personalized treatment strategies offers a scalable, cost-effective way to tailor therapies and improve outcomes for patients with genetically diverse NSCLC.

Humans↗

PubMed2002

Functional selectivity for glycerol of the nodulin 26 subfamily of plant membrane intrinsic proteins.

The nodulin-like intrinsic protein (NIP) subfamily of water and solute channels in plants is named for nodulin 26 of legume nodules. Two NIPs, soybean nodulin 26 and Lotus japonicus LIMP2, show a distinct functional profile with a low intrinsic osmotic water permeability (P(f)) and the ability to flux uncharged polyols such as glycerol. NIPs have a conserved signature sequence within the 'aromatic/arginine' region that forms the selectivity filter for major intrinsic proteins. This sequence is a hybrid of glyceroporin and aquaporin residues as well as exhibiting substitutions unique to the NIP subfamily. Site-directed mutagenesis of a conserved tryptophan in helix 2 of LIMP2 shows that this is a major determinant of glycerol selectivity.

Amino Acid Sequence↗

PubMed2004

N-linked glycosylation sites of the motor protein prestin: effects on membrane targeting and electrophysiological function.

Prestin is a motor protein of outer hair cells (OHC) that plays a crucial role in mammalian hearing. Prestin is a putative N-glycoprotein with three potential N-linked glycosylation sites. It is not known whether glycosylation affects the function and activity of prestin. Therefore, the effects of N-glycosylation were investigated by producing single-point (N163Q and N166Q) or double-point mutations (NN163/166QQ and NN163/166AA) at putative N-glycosylation sites. Further, treatment with tunicamycin or glycopeptidase-F was used to determine the consequences of removing N-linked glycosylation in wild-type prestin. We determined the effects of these manipulations on prestin's cell surface expression, molecular mass, glycosylation pattern, and electrophysiological properties in different cell-types. Data indicate that prestin is a glycoprotein with N-linked glycosylation sites at N163 and N166. N163 and N166 may have differential programs for synthesis and trimming of the glycans. The N166 site appears to have greater extent of glycosylation than its companion. N-linked glycosylation is not required for plasma membrane targeting of prestin. Both glycosylated and deglycosylated prestin demonstrate non-linear capacitance, a signature of prestin's motor function. Compared to glycosylated prestin, the fully de-glycosylated protein has altered electrophysiological function, with a change in membrane potential at most effective charge transfer to more depolarized values. These data suggest that glycosylation of prestin may quantitatively affect OHC electromotility.

Animals↗

PubMed2003

An extended LXXLL motif sequence determines the nuclear receptor binding specificity of TRAP220.

The interaction of coactivators with the ligand-binding domain of nuclear receptors (NRs) is mediated by amphipathic alpha-helices containing the signature motif LXXLL. TRAP220 contains two LXXLL motifs (LXM1 and LXM2) that are required for its interaction with NRs. Here we show that the nuclear receptor interaction domain (NID) of TRAP220 interacts weakly with Class I NRs. In contrast, SRC1 NID binds strongly to both Class I and Class II NRs. Interaction assays using nine amino acid LXXLL core motifs derived from SRC1 and TRAP220 revealed no discriminatory NR binding preferences. However, an extended LXM1 sequence containing amino acids -4 to +9, (where the first conserved leucine is +1) showed selective binding to thyroid hormone receptor and reduced binding to estrogen receptor. Replacement of either TRAP220 LXXLL motif with the corresponding 13 amino acids of SRC1 LXM2 strongly enhanced the interaction of the TRAP220 NID with the estrogen receptor. Mutational analysis revealed combinatorial effects of the LXM1 core and flanking sequences in the determination of the NR binding specificity of the TRAP220 NID. In contrast, a mutation that increased the spacing between TRAP220 LXM1 and LXM2 had little effect on the binding properties of this domain. Thus, a 13-amino acid sequence comprising an extended LXXLL motif acts as the key determinant of the NR binding specificity of TRAP220. Finally, we show that the NR binding specificity of full-length TRAP220 can be altered by swapping extended LXM sequences.

Amino Acid Motifs↗

PubMed2026

A pan-cancer multi-omic SuperLearner for regulated cell death survival topologies.

INTRODUCTION: Regulated cell death (RCD) pathways influence tumor progression and immune modulation. We previously constructed a signature database mapping 25 RCD forms across seven multi-omic layers and 33 tumor types (CancerRCDShiny). Despite their ability to identify risk populations, translating these signatures into personalized clinical workflows requires a shift from cohort stratification to individualized risk mapping by modeling patient risk (survival topologies) to capture the non-linear dynamics of RCD signatures. METHODS: We engineered a pan-cancer multi-omic SuperLearner pipeline across 33 cancer types. Phase I performed zero-leakage harmonization and groupwise imputation to prevent cross-cohort amalgamation. Phase II deployed Elastic Net-regularized Cox regression as a CANARY diagnostic to map proportional hazards failures. Strata with a 35% missingness barrier entered Phase III, deploying a Quadripartite ensemble: Random Survival Forests, XGBoost, Survival-Boruta, and Multi-Task Logistic Regression, fused within an Elastic Net Multi-View Meta-Learner (MVL), with post-hoc TreeSHAP and LIME interpretability. RESULTS: The CANARY diagnostic demonstrated the structural invalidity of pan-cancer geometric proportional hazards. Across 96 admissible strata, Phase III executed algorithmic displacement: continuous multi-omic topologies suppressed static genomic mutations and copy number variations (85.7% vs. 0.0% apex retention). The MVL stabilized predictions against extreme variance; LIME surrogate validations (R 2 < 0.10) confirmed the systematic failure of linear interpretative proxies. N-dimensional TreeSHAP interaction mapping exposed synergistic and antagonistic rescue trajectories defining individualized Survival Topologies, which were invisible to additive models. The architecture was deployed as CancerRCDPredictor, a digital molecular tumor board with integrated LLM capabilities. The MVL SuperLearner achieved a median C-index of 0.749 (IQR: 0.722-0.836) across 96 modelable strata, with 95% bootstrap confidence intervals confirming precision (median width: 0.052) and permutation significance in 93.8% of strata (p < 0.001). External CPTAC validation across ten cancer types demonstrated significant cross-cohort generalizability in clear cell renal carcinoma (KIRC; C-index 0.675, p = 0.017) and modest performance across the remaining adequately powered cancers (median 0.582), underscoring the need for larger multi-institutional validation cohorts. CONCLUSION: This pan-cancer multi-omic SuperLearner bypasses linear topological failures, advancing beyond generalized stratification to establish a deterministically mapped architecture for predicting RCD-related survival topologies. Through the CancerRCDPredictor interface, multi-omic insights translate into individualized survival topology exploration, providing a foundation for future precision oncology validation.

SuperLearner↗

PubMed2004

State of the science in brain tumor classification.

OBJECTIVES: To review the incidence of metastatic and primary brain tumors and the most widely used brain tumor classification systems, and to discuss discoveries advancing the understanding, classification, and grading of selected brain tumor histologies. DATA SOURCES: Journal articles, text books, epidemiologic and statistical reports. CONCLUSION: Recent advances in understanding the molecular biology of brain tumors have shown that molecular and genetic signatures may predict brain tumor behavior and may soon guide tumor classification, diagnosis, and tumor-specific treatment strategies. IMPLICATIONS FOR NURSING PRACTICE: Understanding recent advances in the molecular biology of brain tumors is important because these advances may soon guide treatment decisions. New tumor-specific therapeutic opportunities may improve outcomes as well as the care of persons with brain tumors.

Astrocytoma↗

PubMed2026

Multi-Omics Integration Identifies a Five-Gene Metabolic Signature With Experimental Validation in Clear Cell Renal Cell Carcinoma.

BACKGROUND: Clear cell renal cell carcinoma (ccRCC) is hallmarked by profound metabolic reprogramming; however, its intricate crosstalk with the tumor immune microenvironment (TIME) and its clinical ramifications remain inadequately elucidated. This study aims to systematically decipher the metabolic-immune interplay in ccRCC through multi-omics integration, with the goal of identifying robust prognostic biomarkers and actionable therapeutic vulnerabilities. AIMS: This study aims to systematically decipher the metabolic-immune interplay in clear cell renal cell carcinoma (ccRCC) through multi‑omics integration, and to identify robust prognostic biomarkers and actionable therapeutic vulnerabilities that can inform precision risk stratification and individualized treatment strategies. METHODS: We integrated bulk transcriptomic, genomic, and clinical data from multiple ccRCC cohorts. Differential expression and functional enrichment analyses were performed to characterize metabolic pathway alterations. Mendelian randomization (MR) was employed to infer causal relationships between metabolic disorders and ccRCC risk. A machine learning-based prognostic framework, incorporating SHAP (SHapley Additive exPlanations) for feature interpretability, was constructed and rigorously validated. TIME heterogeneity was dissected using deconvolution algorithms, while drug sensitivity, tumor mutation burden (TMB), and TIDE scores were utilized to assess therapeutic responses and immune evasion. Candidate gene function was evaluated through in vitro gain- and loss-of-function assays, with expression validated via TCGA, HPA, western blot, and qRT-PCR. RESULTS: Enrichment analysis identified coordinated dysregulation in lipid metabolism, energy homeostasis, and hypoxia response pathways. MR analysis confirmed lipid metabolism disorders as a causal risk factor for ccRCC. Our machine-learning model, centered on five core SHAP-identified features (SUCLA2, ACAT1, PC, SUCLG1, and HMGCS2), demonstrated superior predictive accuracy over conventional clinical staging. Immune profiling unveiled dichotomous TIME states: the low-risk group retained active immune surveillance, whereas the high-risk group was enriched with immunosuppressive subsets. Drug sensitivity screening pinpointed LY2109761 and carmustine as high-risk-specific candidate agents. Furthermore, TMB and TIDE analyses stratified high-risk patients displaying genomic instability and immune evasion phenotypes. Functionally, SUCLA2 knockdown significantly enhanced ccRCC cell proliferation and invasion, while its overexpression suppressed these malignant phenotypes, corroborating its tumor-suppressive role. Expression patterns of the hub genes were consistently validated across multi-level datasets and experimental assays. CONCLUSION: This study establishes a precision oncology framework for ccRCC by functionally linking metabolic biomarkers, immunophenotypes, and stratified therapeutic strategies. Importantly, we identify SUCLA2 as a potential functional tumor suppressor and a promising target for further mechanistic and translational investigation.

Humans↗

PubMed2004

A novel database of disulfide patterns and its application to the discovery of distantly related homologs.

Disulfide bonds are conserved strongly among proteins of related structure and function. Despite the explosive growth of protein sequence databases and the vast numbers of sequence search tools, no tool exists to draw relations between the disulfide patterns of homologous proteins. We present a comprehensive database of disulfide bonding patterns and a search method to find proteins with similar disulfide patterns. The disulfide database was constructed using disulfide annotations extracted from SwissProt, and was expanded significantly from 16,736 to 94,499 disulfide-containing domains by an inference method that combines SwissProt annotations with Pfam multiple alignments. To search the database, we define a disulfide description, called the disulfide signature, which encodes both spacings between cysteine residues and cysteine connectivity. A web tool was developed that allows users to search for related disulfide patterns and for subpatterns resulting from the removal of one or more disulfides from the pattern. We explore the possibility of using disulfide pattern conservation to identify protein homologs that are undetectable by PSI-BLAST. Examples include the homology between a sea anemone antihypertensive/antiviral protein and a sea anemone neurotoxin, and the homology between tick anticoagulant peptide and bovine trypsin inhibitor. In both examples, there is a clear structural similarity and a functional relationship. We used the database to find structural homologs for the Cripto CFC domain. The identification of a von Willebrand Factor C (VWFC)-like domain agrees with its functional role and explains mutation data. We believe that the rapid increase in structure determinations arising from structural genomics efforts and advances in mass spectrometry techniques will greatly increase the number of disulfide annotations. This information will become a valuable resource for structural and functional annotations of proteins. The availability of a searchable disulfide pattern database will thus provide a powerful new addition to existing homolog discovery methods.

Amino Acid Sequence↗