Studying the drought-responsive transcriptome is of high interest as it can serve as a blueprint for stress adaptation strategies. Despite extensive studies in this area, there are still many details to be uncovered, such as the importance of each gene involved in the stress response as well as the relationship between these genes and the physiochemical processes governing stress tolerance. This study was designed to address such important details and to gain insights into molecular responses of barley (Hordeum vulgare L.) to drought stress. To that, we combined RNA-seq data analysis with field and greenhouse drought experiments in a systems biology approach. RNA-sequence analysis identified a total of 665 differentially expressed genes (DEGs) belonging to diverse functional categories. A gene network was derived from the DEGs, which comprised of a total of 131 nodes and 257 edges. Gene network topology analysis highlighted two programmed cell death (PCD) modulating genes, MC1 (metacaspase 1) and TSN1 (Tudor-SN 1), as important (hub) genes in the predicted network. Based on the field trial, a drought-tolerant and a drought-susceptible barley genotype was identified from eight tested cultivars. Identified genotypes exhibited different physiochemical characteristics, including proline content, chlorophyll concentration, percentage of electrolyte leakage and malondialdehyde content as well as expression profiles of MC1 and TSN1 genes. Machine learning and correspondence analysis revealed a significant relationship between drought tolerance and measured characteristics in the context of PCD. Our study provides new insights which bridge barley drought tolerance to PCD through MC1 and TSN1 pathway.