Synthesis of fucosyllactose using α-L-fucosidases GH29 from infant gut microbial metagenome.

Fucosyl-oligosaccharides (FUS) provide many health benefits to breastfed infants, but they are almost completely absent from bovine milk, which is the basis of infant formula. Therefore, there is a growing interest in the development of enzymatic transfucosylation strategies for the production of FUS. In this work, the α-L-fucosidases Fuc2358 and Fuc5372, previously isolated from the intestinal bacterial metagenome of breastfed infants, were used to synthesize fucosyllactose (FL) by transfucosylation reactions using p-nitrophenyl-α-L-fucopyranoside (pNP-Fuc) as donor and lactose as acceptor. Fuc2358 efficiently synthesized the major fucosylated human milk oligosaccharide (HMO) 2'-fucosyllactose (2'FL) with a 35% yield. Fuc2358 also produced the non-HMO FL isomer 3'-fucosyllactose (3'FL) and traces of non-reducing 1-fucosyllactose (1FL). Fuc5372 showed a lower transfucosylation activity compared to Fuc2358, producing several FL isomers, including 2'FL, 3'FL, and 1FL, with a higher proportion of 3'FL. Site-directed mutagenesis using rational design was performed to increase FUS yields in both α-L-fucosidases, based on structural models and sequence identity analysis. Mutants Fuc2358-F184H, Fuc2358-K286R, and Fuc5372-R230K showed a significantly higher ratio between 2'FL yields and hydrolyzed pNP-Fuc than their respective wild-type enzymes after 4 h of transfucosylation. The results with the Fuc2358-F184W and Fuc5372-W151F mutants showed that the residues F184 of Fuc2358 and W151 of Fuc5372 could have an effect on transfucosylation regioselectivity. Interestingly, phenylalanine increases the selectivity for α-1,2 linkages and tryptophan for α-1,3 linkages. These results give insight into the functionality of the active site amino acids in the transfucosylation activity of the GH29 α-L-fucosidases Fuc2358 and Fuc5372. KEY POINTS: Two α-L-fucosidases from infant gut bacterial microbiomes can fucosylate glycans Transfucosylation efficacy improved by tailored point-mutations in the active site F184 of Fuc2358 and W151 of Fuc5372 seem to steer transglycosylation regioselectivity.

NMR-Based Mitochondria Metabolomic Profiling: A New Approach To Reveal Cancer-Associated Alterations.

Studying metabolism may assist in understanding the relationship between normal and dysfunctional mitochondrial activity and various diseases, such as neurodegenerative, cardiovascular, autoimmune, psychiatric, and cancer. Nuclear magnetic resonance-based metabolomics represents a powerful method to characterize the chemical content of complex samples and has been successfully applied to studying a range of conditions. However, an optimized methodology is lacking for analyzing isolated organelles, such as mitochondria. In this study, we report the development of a protocol to metabolically profile mitochondria from healthy, tumoral, and metastatic tissues. Encouragingly, this approach provided quantitative information about up to 45 metabolites in one comprehensive and robust analysis. Our results revealed significant differences between whole-cell and mitochondrial metabolites, which supports a more refined approach to metabolic analysis. We applied our optimized methodology to investigate aggressive and metastatic breast cancer in mouse tissues, discovering that lung mitochondria exhibit an altered metabolic fingerprint. Specific amino acids, organic acids, and lipids showed significant increases in levels when compared with mitochondria from healthy tissues. Our optimized methodology could promote a better understanding of the molecular mechanisms underlying breast cancer aggressiveness and mitochondrial-related diseases and support the optimization of new advanced therapies.

Folding of the nascent polypeptide chain of a histidine phosphocarrier protein in vitro.

The phosphotransferase system (PTS), a metabolic pathway formed by five proteins, modulates the use of sugars in bacteria. The second protein in the chain is the histidine phosphocarrier, HPr, with the binding site at His15. The HPr kinase/phosphorylase (HPrK/P), involved in the bacterial use of carbon sources, phosphorylates HPr at Ser46, and it binds at its binding site. The regulator of sigma D protein (Rsd) also binds to HPr at His15. We have designed fragments of HPr, growing from its N-terminus and containing the His15. In this work, we obtained three fragments, HPr38, HPr58 and HPr70, comprising the first thirty-eight, fifty-eight and seventy residues of HPr, respectively. All fragments were mainly disordered, with evidence of a weak native-like, helical population around the binding site, as shown by fluorescence, far-ultraviolet circular dichroism, size exclusion chromatography and nuclear magnetic resonance. Although HPr38, HPr58 and HPr70 were disordered, they could bind to: (i) the N-terminal domain of first protein of the PTS, EIN; (ii) Rsd; and, (iii) HPrK/P, as shown by fluorescence and biolayer interferometry (BLI). The association constants for each protein to any of the fragments were in the low micromolar range, within the same range than those measured in the binding of HPr to each protein. Then, although acquisition of stable, native-like secondary and tertiary structures occurred at the last residues of the polypeptide, the ability to bind protein partners happened much earlier in the growing chain. Binding was related to the presence of the native-like structure around His15.

The muscle-relaxing C-terminal peptide from troponin I populates a nascent helix, facilitating binding to tropomyosin with a potent therapeutic effect.

The conserved C-terminal end segment of troponin I (TnI) plays a critical role in regulating muscle relaxation. This function is retained in the isolated C-terminal 27 amino acid peptide (residues 184-210) of human cardiac TnI (HcTnI-C27): When added to skinned muscle fibers, HcTnI-C27 reduces the Ca2+-sensitivity of activated myofibrils and facilitates relaxation without decreasing the maximum force production. However, the underlying mechanism of HcTnI-C27 function is unknown. We studied the conformational preferences of HcTnI-C27 and a myopathic mutant, Arg192His, (HcTnI-C27-H). Both peptides were mainly disordered in aqueous solution with a nascent helix involving residues from Trp191 to Ile195, as shown by NMR analysis and molecular dynamics simulations. The population of nascent helix was smaller in HcTnI-C27-H than in HcTnI-C27, as shown by circular dichroism (CD) titrations. Fluorescence and isothermal titration calorimetry (ITC) showed that both peptides bound tropomyosin (αTm), with a detectably higher affinity (∼10 µM) of HcTnI-C27 than that of HcTnI-C27-H (∼15 µM), consistent with an impaired Ca2+-desensitization effect of the mutant peptide on skinned muscle strips. Upon binding to αTm, HcTnI-C27 acquired a weakly stable helix-like conformation involving residues near Trp191, as shown by transferred nuclear Overhauser effect spectroscopy and hydrogen/deuterium exchange experiments. With the potent Ca2+-desensitization effect of HcTnI-C27 on skinned cardiac muscle from a mouse model of hypertrophic cardiomyopathy, the data support that the C-terminal end domain of TnI can function as an isolated peptide with the intrinsic capacity of binding tropomyosin, providing a promising therapeutic approach to selectively improve diastolic function of the heart.

Secondary structure determination of conserved SARS-CoV-2 RNA elements by NMR spectroscopy.

The current pandemic situation caused by the Betacoronavirus SARS-CoV-2 (SCoV2) highlights the need for coordinated research to combat COVID-19. A particularly important aspect is the development of medication. In addition to viral proteins, structured RNA elements represent a potent alternative as drug targets. The search for drugs that target RNA requires their high-resolution structural characterization. Using nuclear magnetic resonance (NMR) spectroscopy, a worldwide consortium of NMR researchers aims to characterize potential RNA drug targets of SCoV2. Here, we report the characterization of 15 conserved RNA elements located at the 5' end, the ribosomal frameshift segment and the 3'-untranslated region (3'-UTR) of the SCoV2 genome, their large-scale production and NMR-based secondary structure determination. The NMR data are corroborated with secondary structure probing by DMS footprinting experiments. The close agreement of NMR secondary structure determination of isolated RNA elements with DMS footprinting and NMR performed on larger RNA regions shows that the secondary structure elements fold independently. The NMR data reported here provide the basis for NMR investigations of RNA function, RNA interactions with viral and host proteins and screening campaigns to identify potential RNA binders for pharmaceutical intervention.

Residual Helicity at the Active Site of the Histidine Phosphocarrier, HPr, Modulates Binding Affinity to Its Natural Partners.

The phosphoenolpyruvate-dependent phosphotransferase system (PTS) modulates the preferential use of sugars in bacteria. The first proteins in the cascade are common to all organisms (EI and HPr). The active site of HPr involves a histidine (His15) located immediately before the beginning of the first α-helix. The regulator of sigma D (Rsd) protein also binds to HPr. The region of HPr comprising residues Gly9-Ala30 (HPr9-30), involving the first α-helix (Ala16-Thr27) and the preceding active site loop, binds to both the N-terminal region of EI and intact Rsd. HPr9-30 is mainly disordered. We attempted to improve the affinity of HPr9-30 to both proteins by mutating its sequence to increase its helicity. We designed peptides that led to a marginally larger population in solution of the helical structure of HPr9-30. Molecular simulations also suggested a modest increment in the helical population of mutants, when compared to the wild-type. The mutants, however, were bound with a less favorable affinity than the wild-type to both the N-terminal of EI (EIN) or Rsd, as tested by isothermal titration calorimetry and fluorescence. Furthermore, mutants showed lower antibacterial properties against Staphylococcus aureus than the wild-type peptide. Therefore, we concluded that in HPr, a compromise between binding to its partners and residual structure at the active site must exist to carry out its function.

Mild Muscle Mitochondrial Fusion Distress Extends Drosophila Lifespan through an Early and Systemic Metabolome Reorganization.

In a global aging population, it is important to understand the factors affecting systemic aging and lifespan. Mitohormesis, an adaptive response caused by different insults affecting the mitochondrial network, triggers a response from the nuclear genome inducing several pathways that promote longevity and metabolic health. Understanding the role of mitochondrial function during the aging process could help biomarker identification and the development of novel strategies for healthy aging. Herein, we interfered the muscle expression of the Drosophila genes Marf and Opa1, two genes that encode for proteins promoting mitochondrial fusion, orthologues of human MFN2 and OPA1. Silencing of Marf and Opa1 in muscle increases lifespan, improves locomotor capacities in the long term, and maintains muscular integrity. A metabolomic analysis revealed that muscle down-regulation of Marf and Opa1 promotes a non-autonomous systemic metabolome reorganization, mainly affecting metabolites involved in the energetic homeostasis: carbohydrates, lipids and aminoacids. Interestingly, the differences are consistently more evident in younger flies, implying that there may exist an anticipative adaptation mediating the protective changes at the older age. We demonstrate that mild mitochondrial muscle disturbance plays an important role in Drosophila fitness and reveals metabolic connections between tissues. This study opens new avenues to explore the link of mitochondrial dynamics and inter-organ communication, as well as their relationship with muscle-related pathologies, or in which muscle aging is a risk factor for their appearance. Our results suggest that early intervention in muscle may prevent sarcopenia and promote healthy aging.

Characterization of triple-negative breast cancer preclinical models provides functional evidence of metastatic progression.

Triple-negative breast cancer (TNBC), an aggressive, metastatic and recurrent breast cancer (BC) subtype, currently suffers from a lack of adequately described spontaneously metastatic preclinical models that faithfully reproduce the clinical scenario. We describe two preclinical spontaneously metastatic TNBC orthotopic murine models for the development of advanced therapeutics: an immunodeficient human MDA-MB-231-Luc model and an immunocompetent mouse 4T1 model. Furthermore, we provide a broad range of multifactorial analysis for both models that could provide relevant information for the development of new therapies and diagnostic tools. Our comparisons uncovered differential growth rates, stromal arrangements and metabolic profiles in primary tumors, and the presence of cancer-associated adipocyte infiltration in the MDA-MB-231-Luc model. Histopathological studies highlighted the more rapid metastatic spread to the lungs in the 4T1 model following a lymphatic route, while we observed both homogeneous (MDA-MB-231-Luc) and heterogeneous (4T1) metastatic spread to axillary lymph nodes. We encountered unique metabolomic signatures in each model, including crucial amino acids and cell membrane components. Hematological analysis demonstrated severe leukemoid and lymphoid reactions in the 4T1 model with the partial reestablishment of immune responses in the immunocompromised MDA-MB-231-Luc model. Additionally, we discovered ß-immunoglobulinemia and increased basal levels of G-CSF correlating with a metastatic switch, with G-CSF also promoting extramedullary hematopoiesis (both models) and causing hepatosplenomegaly (4T1 model). Overall, we believe that the characterization of these preclinical models will foster the development of advanced therapeutic strategies for TNBC treatment, especially for the treatment of patients presenting both, primary tumors and metastatic spread.

A translational approach to assess the metabolomic impact of stabilized gold nanoparticles by NMR spectroscopy.

Gold nanoparticles have high potential in the biomedical area, especially in disease diagnosis and treatment. The application of these nanoparticles requires the presence of stabilizers to avoid their agglomeration. Nowadays, there is a lack of reliable methods for characterising the effect of stabilised nanoparticles on biological systems. To this end, in this study, we apply an experimental approach based on nuclear magnetic resonance spectroscopy to study the effect of gold nanoparticles, stabilised with cerium oxide or chitosan, on a human cancer cell model. The results showed that both systems have a significant effect, even at non-toxic levels, on the cellular antioxidant system. However, although particles functionalised with chitosan exerted a strong effect on the aerobic respiration, nanoparticles stabilised with cerium oxide had a higher impact on the mechanisms associated with anaerobic energy production. Therefore, even though both systems contained similar gold nanoparticles, the presence of different stabilizers strongly influenced their mode of action and potential applications in biomedicine.

Bioinformatics tools for the analysis of NMR metabolomics studies focused on the identification of clinically relevant biomarkers.

Metabolomics, a systems biology approach focused on the global study of the metabolome, offers a tremendous potential in the analysis of clinical samples. Among other applications, metabolomics enables mapping of biochemical alterations involved in the pathogenesis of diseases, and offers the opportunity to noninvasively identify diagnostic, prognostic and predictive biomarkers that could translate into early therapeutic interventions. Particularly, metabolomics by Nuclear Magnetic Resonance (NMR) has the ability to simultaneously detect and structurally characterize an abundance of metabolic components, even when their identities are unknown. Analysis of the data generated using this experimental approach requires the application of statistical and bioinformatics tools for the correct interpretation of the results. This review focuses on the different steps involved in the metabolomics characterization of biofluids for clinical applications, ranging from the design of the study to the biological interpretation of the results. Particular emphasis is devoted to the specific procedures required for the processing and interpretation of NMR data with a focus on the identification of clinically relevant biomarkers.

The PipX Protein, When Not Bound to Its Targets, Has Its Signaling C-Terminal Helix in a Flexed Conformation.

PipX, an 89-residue protein, acts as a coactivator of the global nitrogen regulator NtcA in cyanobacteria. NtcA-PipX interactions are regulated by 2-oxoglutarate (2-OG), an inverse indicator of the ammonia abundance, and by PII, a protein that binds to PipX at low 2-OG concentrations. The structure of PipX, when bound to NtcA or PII, consists of an N-terminal, five-stranded ß-sheet (conforming a Tudor-like domain), and two long α-helices. These helices adopt either a flexed conformation, where they are in close contact and in an antiparallel mutual orientation, also packing against the ß-sheet, or an open conformation (observed only in the PII-PipX complex) where the last α-helix moves apart from the rest of the protein. The aim of this work was to study the structure and dynamics of isolated PipX in solution by NMR. The backbone chemical shifts, the hydrogen-exchange, and the NOE patterns indicated that the isolated, monomeric PipX structure was formed by an N-terminal five-stranded ß-sheet and two C-terminal α-helices. Furthermore, the observed NOEs between the two helices, and of α-helix2 with ß-strand2 suggested that PipX adopted a flexed conformation. The ß-strands 1 and 5 were highly flexible, as shown by the lack of interstrand backbone-backbone NOEs; in addition, the 15N-dynamics indicated that the C terminus of ß-strand4 and the following ß-turn (Phe42-Thr47), and the C-cap of α-helix1 (Arg70-Asn71) were particularly mobile. These two regions could act as hinges, allowing PipX to interact with its partners, including PlmA in the newly recognized PII-PipX-PlmA ternary complex.

Direct Monitoring of Exogenous γ-Hydroxybutyric Acid in Body Fluids by NMR Spectroscopy.

γ-Hydroxybutyric acid (GHB) is a popular drug increasingly associated with cases of drug-facilitated sexual assault (DFSA). Currently, expanding procedures of analysis and having forensic evidence of GHB intake in a long term are mandatory. Up to now, most studies have been performed using GC/MS and LC-MS as analytical platforms, which involve significant manipulation of the sample and, often, indirect measurements. In this work, procedures used in NMR-based metabolomics were applied to a GHB clinical trial on urine and serum. Detection, identification, and quantification of the drug by NMR methods were surveyed, as well as the use of NMR-based metabolomics for the search of potential surrogate biomarkers of GHB consumption. Results demonstrated the suitability of NMR spectroscopy, as a robust nondestructive technique, to fast and directly monitor (detect, identify, and quantify) exogenous GHB in almost intact body fluids and its high potential in the search for metabolites associated with GHB intake.

A Drosophila model of GDAP1 function reveals the involvement of insulin signalling in the mitochondria-dependent neuromuscular degeneration.

Charcot-Marie-Tooth disease is a rare peripheral neuropathy for which there is no specific treatment. Some forms of Charcot-Marie-Tooth are due to mutations in the GDAP1 gene. A striking feature of mutations in GDAP1 is that they have a variable clinical manifestation, according to disease onset and progression, histology and mode of inheritance. Studies in cellular and animal models have revealed a role of GDAP1 in mitochondrial morphology and distribution, calcium homeostasis and oxidative stress. To get a better understanding of the disease mechanism we have generated models of over-expression and RNA interference of the Drosophila Gdap1 gene. In order to get an overview about the changes that Gdap1 mutations cause in our disease model, we have combined a comprehensive determination of the metabolic profile in the flies by nuclear magnetic resonance spectroscopy with gene expression analyses and biophysical tests. Our results revealed that both up- and down-regulation of Gdap1 results in an early systemic inactivation of the insulin pathway before the onset of neuromuscular degeneration, followed by an accumulation of carbohydrates and an increase in the ß-oxidation of lipids. Our findings are in line with emerging reports of energy metabolism impairments linked to different types of neural pathologies caused by defective mitochondrial function, which is not surprising given the central role of mitochondria in the control of energy metabolism. The relationship of mitochondrial dynamics with metabolism during neurodegeneration opens new avenues to understand the cause of the disease, and for the discovery of new biomarkers and treatments.

Human milk and mucosal lacto- and galacto-N-biose synthesis by transgalactosylation and their prebiotic potential in Lactobacillus species.

Lacto-N-biose (LNB) and galacto-N-biose (GNB) are major building blocks of free oligosaccharides and glycan moieties of glyco-complexes present in human milk and gastrointestinal mucosa. We have previously characterized the phospho-ß-galactosidase GnbG from Lactobacillus casei BL23 that is involved in the metabolism of LNB and GNB. GnbG has been used here in transglycosylation reactions, and it showed the production of LNB and GNB with N-acetylglucosamine and N-acetylgalactosamine as acceptors, respectively. The reaction kinetics demonstrated that GnbG can convert 69 ± 4 and 71 ± 1 % of o-nitrophenyl-ß-D-galactopyranoside into LNB and GNB, respectively. Those reactions were performed in a semi-preparative scale, and the synthesized disaccharides were purified. The maximum yield obtained for LNB was 10.7 ± 0.2 g/l and for GNB was 10.8 ± 0.3 g/l. NMR spectroscopy confirmed the molecular structures of both carbohydrates and the absence of reaction byproducts, which also supports that GnbG is specific for ß1,3-glycosidic linkages. The purified sugars were subsequently tested for their potential prebiotic properties using Lactobacillus species. The results showed that LNB and GNB were fermented by the tested strains of L. casei, Lactobacillus rhamnosus (except L. rhamnosus strain ATCC 53103), Lactobacillus zeae, Lactobacillus gasseri, and Lactobacillus johnsonii. DNA hybridization experiments suggested that the metabolism of those disaccharides in 9 out of 10 L. casei strains, all L. rhamnosus strains and all L. zeae strains tested relies upon a phospho-ß-galactosidase homologous to GnbG. The results presented here support the putative role of human milk oligosaccharides for selective enrichment of beneficial intestinal microbiota in breast-fed infants.

Three-dimensional solution structure, dynamics and binding of thioredoxin m from Pisum sativum.

Thioredoxins (TRXs) are ubiquitous small, globular proteins involved in cell redox processes. In this work, we report the solution structure of TRX m from Pisum sativum (pea), which has been determined on the basis of 1444 nuclear Overhauser effect- (NOE-) derived distance constraints. The average pairwise root-mean-square deviation (RMSD) for the 20 best structures for the backbone residues (Val7-Glu102) was 1.42 ± 0.15 Å, and 1.97 ± 0.15 Å when all heavy atoms were considered. The structure corresponds to the typical fold of TRXs, with a central five-stranded ß-sheet flanked by four α-helices. Some residues had an important exchange dynamic contribution: those around the active site; at the C terminus of ß-strand 3; and in the loop preceding α-helix 4. Smaller NOE values were observed at the N and C-terminal residues forming the elements of the secondary structure or, alternatively, in the residues belonging to the loops between those elements. A peptide derived from pea fructose-1,6-biphosphatase (FBPase), comprising the preceding region to the regulatory sequence of FBPase (residues Glu152 to Gln179), was bound to TRX m with an affinity in the low micromolar range, as measured by fluorescence and NMR titration experiments. Upon peptide addition, the intensities of the cross-peaks of all the residues of TRX m were affected, as shown by NMR. The value of the dissociation constant of the peptide from TRX m was larger than that of the intact FBPase, indicating that there are additional factors in other regions of the polypeptide chain of the latter protein affecting the binding to thioredoxin.

Optimised protocols for the metabolic profiling of S. cerevisiae by 1H-NMR and HRMAS spectroscopy.

An optimised extraction protocol for the analysis of Saccharomyces cerevisiae aqueous and organic metabolites by nuclear magnetic resonance spectroscopy that allows the identification and quantification of up to 50 different compounds is presented. The method was compared with other metabolic profiling protocols for S. cerevisiae, where generally different analytical techniques are applied for metabolite quantification. In addition, the analysis of intact S. cerevisiae cells by HRMAS was implemented for the first time as a complementary method. The optimised protocols were applied to study the metabolic effect of glucose and galactose on S. cerevisiae growth. Furthermore, the metabolic reaction of S. cerevisiae to osmotic stress has been studied.

Plasma metabolomics by nuclear magnetic resonance reveals biomarkers and metabolic pathways associated with the control of HIV-1 infection/progression.

How the human body reacts to the exposure of HIV-1 is an important research goal. Frequently, HIV exposure leads to infection, but some individuals show natural resistance to this infection; they are known as HIV-1-exposed but seronegative (HESN). Others, although infected but without antiretroviral therapy, control HIV-1 replication and progression to AIDS; they are named controllers, maintaining low viral levels and an adequate count of CD4+ T lymphocytes. Biological mechanisms explaining these phenomena are not precise. In this context, metabolomics emerges as a method to find metabolites in response to pathophysiological stimuli, which can help to establish mechanisms of natural resistance to HIV-1 infection and its progression. We conducted a cross-sectional study including 30 HESN, 14 HIV-1 progressors, 14 controllers and 30 healthy controls. Plasma samples (directly and deproteinized) were analyzed through Nuclear Magnetic Resonance (NMR) metabolomics to find biomarkers and altered metabolic pathways. The metabolic profile analysis of progressors, controllers and HESN demonstrated significant differences with healthy controls when a discriminant analysis (PLS-DA) was applied. In the discriminant models, 13 metabolites associated with HESN, 14 with progressors and 12 with controllers were identified, which presented statistically significant mean differences with healthy controls. In progressors, the metabolites were related to high energy expenditure (creatinine), mood disorders (tyrosine) and immune activation (lipoproteins), phenomena typical of the natural course of the infection. In controllers, they were related to an inflammation-modulating profile (glutamate and pyruvate) and a better adaptive immune system response (acetate) associated with resistance to progression. In the HESN group, with anti-inflammatory (lactate and phosphocholine) and virucidal (lactate) effects which constitute a protective profile in the sexual transmission of HIV. Concerning the significant metabolites of each group, we identified 24 genes involved in HIV-1 replication or virus proteins that were all altered in progressors but only partially in controllers and HESN. In summary, our results indicate that exposure to HIV-1 in HESN, as well as infection in progressors and controllers, affects the metabolism of individuals and that this affectation can be determined using NMR metabolomics.

The Footprint of Type 1 Diabetes on Red Blood Cells: A Metabolomic and Lipidomic Study.

The prevalence of diabetes type 1 (T1D) in the world populations is continuously growing. Although treatment methods are improving, the diagnostic is still symptom-based and sometimes far after onset of the disease. In this context, the aim of the study was the search of new biomarkers of the disease in red blood cells (RBCs), until now unexplored. The metabolomic and the lipidomic profile of RBCs from T1D patients and matched healthy controls was determined by NMR spectroscopy, and different multivariate discrimination models were built to select the metabolites and lipids that change most significantly. Relevant metabolites were further confirmed by univariate statistical analysis. Robust separation in the metabolomic and lipidomic profiles of RBCs from patients and controls was confirmed by orthogonal projection on latent structure discriminant analysis (OPLS-DA), random forest analysis, and significance analysis of metabolites (SAM). The main changes were detected in the levels of amino acids, organic acids, creatine and phosphocreatine, lipid change length, and choline derivatives, demonstrating changes in glycolysis, BCAA metabolism, and phospholipid metabolism. Our study proves that robust differences exist in the metabolic and lipidomic profile of RBCs from T1D patients, in comparison with matched healthy individuals. Some changes were similar to alterations found already in RBCs of T2D patients, but others seemed to be specific for type 1 diabetes. Thus, many of the metabolic differences found could be biomarker candidates for an earlier diagnosis or monitoring of patients with T1D.