The triclustering method employed in this study integrates the δ-Trimax approach with the fuzzy cuckoo search (FCS), thereby leveraging the Lévy flight and Gaussian distribution to analyze gene expression data in three dimensions. In this framework, the initial triclusters produced by δ-Trimax are further optimized using FCS, where the Lévy flight enhances global exploration and the Gaussian distribution intensifies local exploitation, thus achieving a balanced search for optimal solutions. Each tricluster set is evaluated using the tricluster quality index (TQI) to ensure coherence across genes, conditions, and time points. The method was applied to gene expression datasets from primary fibroblast cells and heart disease samples. In the fibroblast dataset, the best tricluster set was obtained with δ = 0.015 and yielded the lowest average TQI value. For the heart disease dataset, the most optimal solution was achieved with δ = 0.026, which yielded the lowest average TQI, and the best tricluster showed large gene coverage across multiple time points. A functional analysis of the selected triclusters using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways uncovered significant enrichment in pathways such as the NF-κB signaling pathway (hsa04064), TGF-β signaling pathway (hsa04350), and calcium signaling pathway (hsa04020), all of which are mechanistically relevant to immune modulation, extracellular matrix organization, and cardiac muscle function. These findings highlight the utility of the proposed hybrid framework in uncovering biologically meaningful gene modules and provide valuable insights into the molecular mechanisms underlying fibrotic and cardiovascular diseases. © 2026 the Author(s),