The complex pathology of Alzheimer's disease (AD) emphasises the need for comprehensive modelling of the disease, which may lead to the development of efficient treatment strategies. To address this challenge, we analysed transcriptome data of post-mortem human brain samples of healthy elders and individuals with late-onset AD from the Religious Orders Study and Rush Memory and Aging Project (ROSMAP) and Mayo Clinic (MayoRNAseq) studies in the AMP-AD consortium. In this context, we conducted several bioinformatics and systems medicine analyses including the construction of AD-specific co-expression networks and genome-scale metabolic modelling of the brain in AD patients to identify key genes, metabolites and pathways involved in the progression of AD. We identified AMIGO1 and GRPRASP2 as examples of commonly altered marker genes in AD patients. Moreover, we found alterations in energy metabolism, represented by reduced oxidative phosphorylation and ATPase activity, as well as the depletion of hexanoyl-CoA, pentanoyl-CoA, (2E)-hexenoyl-CoA and numerous other unsaturated fatty acids in the brain. We also observed that neuroprotective metabolites (e.g., vitamins, retinoids and unsaturated fatty acids) tend to be depleted in the AD brain, while neurotoxic metabolites (e.g., beta-alanine, bilirubin) were more abundant. In summary, we systematically revealed the key genes and pathways related to the progression of AD, gained insight into the crucial mechanisms of AD and identified some possible targets that could be used in the treatment of AD.