Fabry disease (FD) is an X-linked lysosomal disease whose ultimate consequences are the accumulation of sphingolipids and subsequent inflammatory events, mainly at the endothelial level. The outcomes include different nervous system manifestations as well as multiple organ damage. Despite the availability of known biomarkers, early detection of FD remains a medical need. This study aimed to develop an in silico model based on machine learning to identify candidate vascular and nervous system proteins for early FD damage detection at the cellular level. A combined systems biology and machine learning approach was carried out considering molecular characteristics of FD to create a computational model of vascular and nervous system disease. A data science strategy was applied to identify risk classifiers by using 10 K-fold cross-validation. Further biological and clinical criteria were used to prioritize the most promising candidates, resulting in the identification of 36 biomarker candidates with classifier abilities, which are easily measurable in body fluids. Among them, we propose four candidates, CAMK2A, ILK, LMNA, and KHSRP, which have high classification capabilities according to our models (cross-validated accuracy ≥ 90%) and are related to the vascular and nervous systems. These biomarkers show promise as high-risk cellular and tissue damage indicators that are potentially applicable in clinical settings, although in vivo validation is still needed.

Type 1 Gaucher disease (GD1) is a rare, autosomal recessive disorder caused by glucocerebrosidase deficiency. Skeletal manifestations represent one of the most debilitating and potentially irreversible complications of GD1. Although imaging studies are the gold standard, early diagnostic/prognostic tools, such as molecular biomarkers, are needed for the rapid management of skeletal complications. This study aimed to identify potential protein biomarkers capable of predicting the early diagnosis of bone skeletal complications in GD1 patients using artificial intelligence. An in silico study was performed using the novel Therapeutic Performance Mapping System methodology to construct mathematical models of GD1-associated complications at the protein level. Pathophysiological characterization was performed before modeling, and a data science strategy was applied to the predicted protein activity for each protein in the models to identify classifiers. Statistical criteria were used to prioritize the most promising candidates, and 18 candidates were identified. Among them, PDGFB, IL1R2, PTH and CCL3 (MIP-1α) were highlighted due to their ease of measurement in blood. This study proposes a validated novel tool to discover new protein biomarkers to support clinician decision-making in an area where medical needs have not yet been met. However, confirming the results using in vitro and/or in vivo studies is necessary.