Exposure to PFAS chemicals induces sex-dependent alterations in key rate-limiting steps of lipid metabolism in liver steatosis.

Toxicants with the potential to bioaccumulate in humans and animals have long been a cause for concern, particularly due to their association with multiple diseases and organ injuries. Per- and polyfluoro alkyl substances (PFAS) and polycyclic aromatic hydrocarbons (PAH) are two such classes of chemicals that bioaccumulate and have been associated with steatosis in the liver. Although PFAS and PAH are classified as chemicals of concern, their molecular mechanisms of toxicity remain to be explored in detail. In this study, we aimed to identify potential mechanisms by which an acute exposure to PFAS and PAH chemicals can induce lipid accumulation and whether the responses depend on chemical class, dose, and sex. To this end, we analyzed mechanisms beginning with the binding of the chemical to a molecular initiating event (MIE) and the consequent transcriptomic alterations. We collated potential MIEs using predictions from our previously developed ToxProfiler tool and from published steatosis adverse outcome pathways. Most of the MIEs are transcription factors, and we collected their target genes by mining the TRRUST database. To analyze the effects of PFAS and PAH on the steatosis mechanisms, we performed a computational MIE-target gene analysis on high-throughput transcriptomic measurements of liver tissue from male and female rats exposed to either a PFAS or PAH. The results showed peroxisome proliferator-activated receptor (PPAR)-α targets to be the most dysregulated, with most of the genes being upregulated. Furthermore, PFAS exposure disrupted several lipid metabolism genes, including upregulation of fatty acid oxidation genes (Acadm, Acox1, Cpt2, Cyp4a1-3) and downregulation of lipid transport genes (Apoa1, Apoa5, Pltp). We also identified multiple genes with sex-specific behavior. Notably, the rate-limiting genes of gluconeogenesis (Pck1) and bile acid synthesis (Cyp7a1) were specifically downregulated in male rats compared to female rats, while the rate-limiting gene of lipid synthesis (Scd) showed a PFAS-specific upregulation. The results suggest that the PPAR signaling pathway plays a major role in PFAS-induced lipid accumulation in rats. Together, these results show that PFAS exposure induces a sex-specific multi-factorial mechanism involving rate-limiting genes of gluconeogenesis and bile acid synthesis that could lead to activation of an adverse outcome pathway for steatosis.

Probing Liver Injuries Induced by Thioacetamide in Human In Vitro Pooled Hepatocyte Experiments.

Animal studies are typically utilized to understand the complex mechanisms associated with toxicant-induced hepatotoxicity. Among the alternative approaches to animal studies, in vitro pooled human hepatocytes have the potential to capture population variability. Here, we examined the effect of the hepatotoxicant thioacetamide on pooled human hepatocytes, divided into five lots, obtained from forty diverse donors. For 24 h, pooled human hepatocytes were exposed to vehicle, 1.33 mM (low dose), and 12 mM (high dose) thioacetamide, followed by RNA-seq analysis. We assessed gene expression variability using heat maps, correlation plots, and statistical variance. We used KEGG pathways and co-expression modules to identify underlying physiological processes/pathways. The co-expression module analysis showed that the majority of the lots exhibited activation for the bile duct proliferation module. Despite lot-to-lot variability, we identified a set of common differentially expressed genes across the lots with similarities in their response to amino acid, lipid, and carbohydrate metabolism. We also examined efflux transporters and found larger lot-to-lot variability in their expression patterns, indicating a potential for alteration in toxicant bioavailability within the cells, which could in turn affect the gene expression patterns between the lots. Overall, our analysis highlights the challenges in using pooled hepatocytes to understand mechanisms of toxicity.

AI algorithm for personalized resource allocation and treatment of hemorrhage casualties.

A deep neural network-based artificial intelligence (AI) model was assessed for its utility in predicting vital signs of hemorrhage patients and optimizing the management of fluid resuscitation in mass casualties. With the use of a cardio-respiratory computational model to generate synthetic data of hemorrhage casualties, an application was created where a limited data stream (the initial 10 min of vital-sign monitoring) could be used to predict the outcomes of different fluid resuscitation allocations 60 min into the future. The predicted outcomes were then used to select the optimal resuscitation allocation for various simulated mass-casualty scenarios. This allowed the assessment of the potential benefits of using an allocation method based on personalized predictions of future vital signs versus a static population-based method that only uses currently available vital-sign information. The theoretical benefits of this approach included up to 46% additional casualties restored to healthy vital signs and a 119% increase in fluid-utilization efficiency. Although the study is not immune from limitations associated with synthetic data under specific assumptions, the work demonstrated the potential for incorporating neural network-based AI technologies in hemorrhage detection and treatment. The simulated injury and treatment scenarios used delineated possible benefits and opportunities available for using AI in pre-hospital trauma care. The greatest benefit of this technology lies in its ability to provide personalized interventions that optimize clinical outcomes under resource-limited conditions, such as in civilian or military mass-casualty events, involving moderate and severe hemorrhage.

Identifying metabolic adaptations characteristic of cardiotoxicity using paired transcriptomics and metabolomics data integrated with a computational model of heart metabolism.

Improvements in the diagnosis and treatment of cancer have revealed long-term side effects of chemotherapeutics, particularly cardiotoxicity. Here, we present paired transcriptomics and metabolomics data characterizing in vitro cardiotoxicity to three compounds: 5-fluorouracil, acetaminophen, and doxorubicin. Standard gene enrichment and metabolomics approaches identify some commonly affected pathways and metabolites but are not able to readily identify metabolic adaptations in response to cardiotoxicity. The paired data was integrated with a genome-scale metabolic network reconstruction of the heart to identify shifted metabolic functions, unique metabolic reactions, and changes in flux in metabolic reactions in response to these compounds. Using this approach, we confirm previously seen changes in the p53 pathway by doxorubicin and RNA synthesis by 5-fluorouracil, we find evidence for an increase in phospholipid metabolism in response to acetaminophen, and we see a shift in central carbon metabolism suggesting an increase in metabolic demand after treatment with doxorubicin and 5-fluorouracil.

High-Throughput Transcriptomics Differentiates Toxic versus Non-Toxic Chemical Exposures Using a Rat Liver Model.

To address the challenge of limited throughput with traditional toxicity testing, a newly developed high-throughput transcriptomics (HTT) platform, together with a 5-day in vivo rat model, offers an alternative approach to estimate chemical exposures and provide reasonable estimates of toxicological endpoints. This study contains an HTT analysis of 18 environmental chemicals with known liver toxicity. They were evaluated using male Sprague Dawley rats exposed to various concentrations daily for five consecutive days via oral gavage, with data collected on the sixth day. Here, we further explored the 5-day rat model to identify potential gene signatures that can differentiate between toxic and non-toxic liver responses and provide us with a potential histopathological endpoint of chemical exposure. We identified a distinct gene expression pattern that differentiated non-hepatotoxic compounds from hepatotoxic compounds in a dose-dependent manner, and an analysis of the significantly altered common genes indicated that toxic chemicals predominantly upregulated most of the genes and several pathways in amino acid and lipid metabolism. Finally, our liver injury module analysis revealed that several liver-toxic compounds showed similarities in the key injury phenotypes of cellular inflammation and proliferation, indicating potential molecular initiating processes that may lead to a specific end-stage liver disease.

A stochastic B cell affinity maturation model to characterize mechanisms of protection for tetravalent dengue vaccine constructs.

Dengue annually infects millions of people from a regionally and seasonally varying dengue virus population circulating as four distinct serotypes. Effective protection against dengue infection and disease requires tetravalent vaccine formulations to stimulate a balanced protective immune response to all four serotypes. However, this has been a challenge to achieve, and several clinical trials with different leading vaccine candidates have demonstrated unbalanced replication and interference of interindividual serotype components, leading to low efficacy and enhanced disease severity for dengue-naïve populations. Production of serotype-specific neutralizing antibodies is largely viewed as a correlate of protection against severe dengue disease. However, the underlying mechanisms that lead to these protective immune responses are not clearly elucidated. In this work, using a stochastic model of B cell affinity maturation, we tested different live-attenuated vaccine constructs with varied viral replication rates and contrasted the initiation and progress of adaptive immune responses during tetravalent vaccination and after dengue virus challenge. Comparison of our model simulations across different disease-severity levels suggested that individual production of high levels of serotype-specific antibodies together with a lower cross-reactive antibody are better correlates for protection. Furthermore, evolution of these serotype-specific antibodies was dependent on the percent of viral attenuation in the vaccine, and production of initial B cell and T cell populations pre- and post-secondary dengue infection was crucial in providing protective immunity for dengue-naïve populations. Furthermore, contrasting disease severity with respect to different dengue serotypes, our model simulations showed that tetravalent vaccines fare better against DENV-4 serotype when compared to other serotypes.

Using a Graph Convolutional Neural Network Model to Identify Bile Salt Export Pump Inhibitors.

The bile salt export pump (BSEP) is a key transporter involved in the efflux of bile salts from hepatocytes to bile canaliculi. Inhibition of BSEP leads to the accumulation of bile salts within the hepatocytes, leading to possible cholestasis and drug-induced liver injury. Screening for and identification of chemicals that inhibit this transporter aid in understanding the safety liabilities of these chemicals. Moreover, computational approaches to identify BSEP inhibitors provide an alternative to the more resource-intensive, gold standard experimental approaches. Here, we used publicly available data to develop predictive machine learning models for the identification of potential BSEP inhibitors. Specifically, we analyzed the utility of a graph convolutional neural network (GCNN)-based approach in combination with multitask learning to identify BSEP inhibitors. Our analyses showed that the developed GCNN model performed better than the variable-nearest neighbor and Bayesian machine learning approaches, with a cross-validation receiver operating characteristic area under the curve of 0.86. In addition, we compared GCNN-based single-task and multitask models and evaluated their utility in addressing data limitation challenges commonly observed in bioactivity modeling. We found that multitask models performed better than single-task models and can be utilized to identify active molecules for targets with limited data availability. Overall, our developed multitask GCNN-based BSEP model provides a useful tool for prioritizing hits during early drug discovery and in risk assessment of chemicals.

Assessing Kidney Injury Induced by Mercuric Chloride in Guinea Pigs with In Vivo and In Vitro Experiments.

Acute kidney injury, which is associated with high levels of morbidity and mortality, affects a significant number of individuals, and can be triggered by multiple factors, such as medications, exposure to toxic chemicals or other substances, disease, and trauma. Because the kidney is a critical organ, understanding and identifying early cellular or gene-level changes can provide a foundation for designing medical interventions. In our earlier work, we identified gene modules anchored to histopathology phenotypes associated with toxicant-induced liver and kidney injuries. Here, using in vivo and in vitro experiments, we assessed and validated these kidney injury-associated modules by analyzing gene expression data from the kidneys of male Hartley guinea pigs exposed to mercuric chloride. Using plasma creatinine levels and cell-viability assays as measures of the extent of renal dysfunction under in vivo and in vitro conditions, we performed an initial range-finding study to identify the appropriate doses and exposure times associated with mild and severe kidney injuries. We then monitored changes in kidney gene expression at the selected doses and time points post-toxicant exposure to characterize the mechanisms of kidney injury. Our injury module-based analysis revealed a dose-dependent activation of several phenotypic cellular processes associated with dilatation, necrosis, and fibrogenesis that were common across the experimental platforms and indicative of processes that initiate kidney damage. Furthermore, a comparison of activated injury modules between guinea pigs and rats indicated a strong correlation between the modules, highlighting their potential for cross-species translational studies.

Metabolic responses in blood-stage malaria parasites associated with increased and decreased sensitivity to PfATP4 inhibitors.

BACKGROUND: Spiroindolone and pyrazoleamide antimalarial compounds target Plasmodium falciparum P-type ATPase (PfATP4) and induce disruption of intracellular Na+ homeostasis. Recently, a PfATP4 mutation was discovered that confers resistance to a pyrazoleamide while increasing sensitivity to a spiroindolone. Transcriptomic and metabolic adaptations that underlie this seemingly contradictory response of P. falciparum to sublethal concentrations of each compound were examined to understand the different cellular accommodation to PfATP4 disruptions. METHODS: A genetically engineered P. falciparum Dd2 strain (Dd2A211V) carrying an Ala211Val (A211V) mutation in PfATP4 was used to identify metabolic adaptations associated with the mutation that results in decreased sensitivity to PA21A092 (a pyrazoleamide) and increased sensitivity to KAE609 (a spiroindolone). First, sublethal doses of PA21A092 and KAE609 causing substantial reduction (30-70%) in Dd2A211V parasite replication were identified. Then, at this sublethal dose of PA21A092 (or KAE609), metabolomic and transcriptomic data were collected during the first intraerythrocytic developmental cycle. Finally, the time-resolved data were integrated with a whole-genome metabolic network model of P. falciparum to characterize antimalarial-induced physiological adaptations. RESULTS: Sublethal treatment with PA21A092 caused significant (p < 0.001) alterations in the abundances of 91 Plasmodium gene transcripts, whereas only 21 transcripts were significantly altered due to sublethal treatment with KAE609. In the metabolomic data, a substantial alteration (≥ fourfold) in the abundances of carbohydrate metabolites in the presence of either compound was found. The estimated rates of macromolecule syntheses between the two antimalarial-treated conditions were also comparable, except for the rate of lipid synthesis. A closer examination of parasite metabolism in the presence of either compound indicated statistically significant differences in enzymatic activities associated with synthesis of phosphatidylcholine, phosphatidylserine, and phosphatidylinositol. CONCLUSION: The results of this study suggest that malaria parasites activate protein kinases via phospholipid-dependent signalling in response to the ionic perturbation induced by the Na+ homeostasis disruptor PA21A092. Therefore, targeted disruption of phospholipid signalling in PA21A092-resistant parasites could be a means to block the emergence of resistance to PA21A092.

Development and validation of a mathematical model to simulate human cardiovascular and respiratory responses to battlefield trauma.

Mathematical models of human cardiovascular and respiratory systems provide a viable alternative to generate synthetic data to train artificial intelligence (AI) clinical decision-support systems and assess closed-loop control technologies, for military medical applications. However, existing models are either complex, standalone systems that lack the interface to other applications or fail to capture the essential features of the physiological responses to the major causes of battlefield trauma (i.e., hemorrhage and airway compromise). To address these limitations, we developed the cardio-respiratory (CR) model by expanding and integrating two previously published models of the cardiovascular and respiratory systems. We compared the vital signs predicted by the CR model with those from three models, using experimental data from 27 subjects in five studies, involving hemorrhage, fluid resuscitation, and respiratory perturbations. Overall, the CR model yielded relatively small root mean square errors (RMSEs) for mean arterial pressure (MAP; 20.88 mm Hg), end-tidal CO2 (ETCO2 ; 3.50 mm Hg), O2 saturation (SpO2 ; 3.40%), and arterial O2 pressure (PaO2 ; 10.06 mm Hg), but a relatively large RMSE for heart rate (HR; 70.23 beats/min). In addition, the RMSEs for the CR model were 3% to 10% smaller than the three other models for HR, 11% to 15% for ETCO2 , 0% to 33% for SpO2 , and 10% to 64% for PaO2 , while they were similar for MAP. In conclusion, the CR model balances simplicity and accuracy, while qualitatively and quantitatively capturing human physiological responses to battlefield trauma, supporting its use to train and assess emerging AI and control systems.

Enabling data-limited chemical bioactivity predictions through deep neural network transfer learning.

The main limitation in developing deep neural network (DNN) models to predict bioactivity properties of chemicals is the lack of sufficient assay data to train the network's classification layers. Focusing on feedforward DNNs that use atom- and bond-based structural fingerprints as input, we examined whether layers of a fully trained DNN based on large amounts of data to predict one property could be used to develop DNNs to predict other related or unrelated properties based on limited amounts of data. Hence, we assessed if and under what conditions the dense layers of a pre-trained DNN could be transferred and used for the development of another DNN associated with limited training data. We carried out a quantitative study employing more than 400 pairs of assay datasets, where we used fully trained layers from a large dataset to augment the training of a small dataset. We found that the higher the correlation r between two assay datasets, the more efficient the transfer learning is in reducing prediction errors associated with the smaller dataset DNN predictions. The reduction in mean squared prediction errors ranged from 10 to 20% for every 0.1 increase in r2 between the datasets, with the bulk of the error reductions associated with transfers of the first dense layer. Transfer of other dense layers did not result in additional benefits, suggesting that deeper, dense layers conveyed more specialized and assay-specific information. Importantly, depending on the dataset correlation, training sample size could be reduced by up to tenfold without any loss of prediction accuracy.

Age-dependent antibody profiles to plasmodium antigens are differentially associated with two artemisinin combination therapy outcomes in high transmission setting.

The impact of pre-existing immunity on the efficacy of artemisinin combination therapy is largely unknown. We performed in-depth profiling of serological responses in a therapeutic efficacy study [comparing artesunate-mefloquine (ASMQ) and artemether-lumefantrine (AL)] using a proteomic microarray. Responses to over 200 Plasmodium antigens were significantly associated with ASMQ treatment outcome but not AL. We used machine learning to develop predictive models of treatment outcome based on the immunoprofile data. The models predict treatment outcome for ASMQ with high (72-85%) accuracy, but could not predict treatment outcome for AL. This divergent treatment outcome suggests that humoral immunity may synergize with the longer mefloquine half-life to provide a prophylactic effect at 28-42 days post-treatment, which was further supported by simulated pharmacokinetic profiling. Our computational approach and modeling revealed the synergistic effect of pre-existing immunity in patients with drug combination that has an extended efficacy on providing long term treatment efficacy of ASMQ.

Metabolic changes accompanying the loss of fumarate hydratase and malate-quinone oxidoreductase in the asexual blood stage of Plasmodium falciparum.

In the glucose-rich milieu of red blood cells, asexually replicating malarial parasites mainly rely on glycolysis for ATP production, with limited carbon flux through the mitochondrial tricarboxylic acid (TCA) cycle. By contrast, gametocytes and mosquito-stage parasites exhibit an increased dependence on the TCA cycle and oxidative phosphorylation for more economical energy generation. Prior genetic studies supported these stage-specific metabolic preferences by revealing that six of eight TCA cycle enzymes are completely dispensable during the asexual blood stages of Plasmodium falciparum, with only fumarate hydratase (FH) and malate-quinone oxidoreductase (MQO) being refractory to deletion. Several hypotheses have been put forth to explain the possible essentiality of FH and MQO, including their participation in a malate shuttle between the mitochondrial matrix and the cytosol. However, using newer genetic techniques like CRISPR and dimerizable Cre, we were able to generate deletion strains of FH and MQO in P. falciparum. We employed metabolomic analyses to characterize a double knockout mutant of FH and MQO (ΔFM) and identified changes in purine salvage and urea cycle metabolism that may help to limit fumarate accumulation. Correspondingly, we found that the ΔFM mutant was more sensitive to exogenous fumarate, which is known to cause toxicity by modifying and inactivating proteins and metabolites. Overall, our data indicate that P. falciparum is able to adequately compensate for the loss of FH and MQO, rendering them unsuitable targets for drug development.

Search for a Shared Genetic or Biochemical Basis for Biofilm Tolerance to Antibiotics across Bacterial Species.

Is there a universal genetically programmed defense providing tolerance to antibiotics when bacteria grow as biofilms? A comparison between biofilms of three different bacterial species by transcriptomic and metabolomic approaches uncovered no evidence of one. Single-species biofilms of three bacterial species (Pseudomonas aeruginosa, Staphylococcus aureus, and Acinetobacter baumannii) were grown in vitro for 3 days and then challenged with respective antibiotics (ciprofloxacin, daptomycin, and tigecycline) for an additional 24 h. All three microorganisms displayed reduced susceptibility in biofilms compared to planktonic cultures. Global transcriptomic profiling of gene expression comparing biofilm to planktonic and antibiotic-treated biofilm to untreated biofilm was performed. Extracellular metabolites were measured to characterize the utilization of carbon sources between biofilms, treated biofilms, and planktonic cells. While all three bacteria exhibited a species-specific signature of stationary phase, no conserved gene, gene set, or common functional pathway could be identified that changed consistently across the three microorganisms. Across the three species, glucose consumption was increased in biofilms compared to planktonic cells, and alanine and aspartic acid utilization were decreased in biofilms compared to planktonic cells. The reasons for these changes were not readily apparent in the transcriptomes. No common shift in the utilization pattern of carbon sources was discerned when comparing untreated to antibiotic-exposed biofilms. Overall, our measurements do not support the existence of a common genetic or biochemical basis for biofilm tolerance against antibiotics. Rather, there are likely myriad genes, proteins, and metabolic pathways that influence the physiological state of individual microorganisms in biofilms and contribute to antibiotic tolerance.

Metabolic adjustments of blood-stage Plasmodium falciparum in response to sublethal pyrazoleamide exposure.

Due to the recurring loss of antimalarial drugs to resistance, there is a need for novel targets, drugs, and combination therapies to ensure the availability of current and future countermeasures. Pyrazoleamides belong to a novel class of antimalarial drugs that disrupt sodium ion homeostasis, although the exact consequences of this disruption in Plasmodium falciparum remain under investigation. In vitro experiments demonstrated that parasites carrying mutations in the metabolic enzyme PfATP4 develop resistance to pyrazoleamide compounds. However, the underlying mechanisms that allow mutant parasites to evade pyrazoleamide treatment are unclear. Here, we first performed experiments to identify the sublethal dose of a pyrazoleamide compound (PA21A092) that caused a significant reduction in growth over one intraerythrocytic developmental cycle (IDC). At this drug concentration, we collected transcriptomic and metabolomic data at multiple time points during the IDC to quantify gene- and metabolite-level alterations in the treated parasites. To probe the effects of pyrazoleamide treatment on parasite metabolism, we coupled the time-resolved omics data with a metabolic network model of P. falciparum. We found that the drug-treated parasites adjusted carbohydrate metabolism to enhance synthesis of myoinositol-a precursor for phosphatidylinositol biosynthesis. This metabolic adaptation caused a decrease in metabolite flux through the pentose phosphate pathway, causing a decreased rate of RNA synthesis and an increase in oxidative stress. Our model analyses suggest that downstream consequences of enhanced myoinositol synthesis may underlie adjustments that could lead to resistance emergence in P. falciparum exposed to a sublethal dose of a pyrazoleamide drug.

Genomics and metabolomics of early-stage thioacetamide-induced liver injury: An interspecies study between guinea pig and rat.

To study the complex processes involved in liver injuries, researchers rely on animal investigations, using chemically or surgically induced liver injuries, to extrapolate findings and infer human health risks. However, this presents obvious challenges in performing a detailed comparison and validation between the highly controlled animal models and development of liver injuries in humans. Furthermore, it is not clear whether there are species-dependent and -independent molecular initiating events or processes that cause liver injury before they eventually lead to end-stage liver disease. Here, we present a side-by-side study of rats and guinea pigs using thioacetamide to examine the similarities between early molecular initiating events during an acute-phase liver injury. We exposed Sprague Dawley rats and Hartley guinea pigs to a single dose of 25 or 100 mg/kg thioacetamide and collected blood plasma for metabolomic analysis and liver tissue for RNA-sequencing. The subsequent toxicogenomic analysis identified consistent liver injury trends in both genomic and metabolomic data within 24 and 33 h after thioacetamide exposure in rats and guinea pigs, respectively. In particular, we found species similarities in the key injury phenotypes of inflammation and fibrogenesis in our gene module analysis for liver injury phenotypes. We identified expression of several common genes (e.g., SPP1, TNSF18, SERPINE1, CLDN4, TIMP1, CD44, and LGALS3), activation of injury-specific KEGG pathways, and alteration of plasma metabolites involved in amino acid and bile acid metabolism as some of the key molecular processes that changed early upon thioacetamide exposure and could play a major role in the initiation of acute liver injury.

Stochastic Model of the Adaptive Immune Response Predicts Disease Severity and Captures Enhanced Cross-Reactivity in Natural Dengue Infections.

The dengue virus circulates as four distinct serotypes, where a single serotype infection is typically asymptomatic and leads to acquired immunity against that serotype. However, the developed immunity to one serotype is thought to underlie the severe manifestation of the disease observed in subsequent infections from a different serotype. We developed a stochastic model of the adaptive immune response to dengue infections. We first delineated the mechanisms initiating and sustaining adaptive immune responses during primary infections. We then contrasted these immune responses during secondary infections of either a homotypic or heterotypic serotype to understand the role of pre-existing and reactivated immune pathways on disease severity. Comparison of non-symptomatic and severe cases from heterotypic infections demonstrated that overproduction of specific antibodies during primary infection induces an enhanced population of cross-reactive antibodies during secondary infection, ultimately leading to severe disease manifestations. In addition, the level of disease severity was found to correlate with immune response kinetics, which was dependent on beginning lymphocyte levels. Our results detail the contribution of specific lymphocytes and antibodies to immunity and memory recall that lead to either protective or pathological outcomes, allowing for the understanding and determination of mechanisms of protective immunity.

Inter-study and time-dependent variability of metabolite abundance in cultured red blood cells.

BACKGROUND: Cultured human red blood cells (RBCs) provide a powerful ex vivo assay platform to study blood-stage malaria infection and propagation. In recent years, high-resolution metabolomic methods have quantified hundreds of metabolites from parasite-infected RBC cultures under a variety of perturbations. In this context, the corresponding control samples of the uninfected culture systems can also be used to examine the effects of these perturbations on RBC metabolism itself and their dependence on blood donors (inter-study variations). METHODS: Time-course datasets from five independent studies were generated and analysed, maintaining uninfected RBCs (uRBC) at 2% haematocrit for 48 h under conditions originally designed for parasite cultures. Using identical experimental protocols, quadruplicate samples were collected at six time points, and global metabolomics were employed on the pellet fraction of the uRBC cultures. In total, ~ 500 metabolites were examined across each dataset to quantify inter-study variability in RBC metabolism, and metabolic network modelling augmented the analyses to characterize the metabolic state and fluxes of the RBCs. RESULTS: To minimize inter-study variations unrelated to RBC metabolism, an internal standard metabolite (phosphatidylethanolamine C18:0/20:4) was identified with minimal variation in abundance over time and across all the samples of each dataset to normalize the data. Although the bulk of the normalized data showed a high degree of inter-study consistency, changes and variations in metabolite levels from individual donors were noted. Thus, a total of 24 metabolites were associated with significant variation in the 48-h culture time window, with the largest variations involving metabolites in glycolysis and synthesis of glutathione. Metabolic network analysis was used to identify the production of superoxide radicals in cultured RBCs as countered by the activity of glutathione oxidoreductase and synthesis of reducing equivalents via the pentose phosphate pathway. Peptide degradation occurred at a rate that is comparable with central carbon fluxes, consistent with active degradation of methaemoglobin, processes also commonly associated with storage lesions in RBCs. CONCLUSIONS: The bulk of the data showed high inter-study consistency. The collected data, quantification of an expected abundance variation of RBC metabolites, and characterization of a subset of highly variable metabolites in the RBCs will help in identifying non-specific changes in metabolic abundances that may obscure accurate metabolomic profiling of Plasmodium falciparum and other blood-borne pathogens.

Immunoprofiling Correlates of Protection Against SHIV Infection in Adjuvanted HIV-1 Pox-Protein Vaccinated Rhesus Macaques.

Human immunodeficiency virus type 1 (HIV-1) infection remains a major public health threat due to its incurable nature and the lack of a highly efficacious vaccine. The RV144 vaccine trial is the only clinical study to date that demonstrated significant but modest decrease in HIV infection risk. To improve HIV-1 vaccine immunogenicity and efficacy, we recently evaluated pox-protein vaccination using a next generation liposome-based adjuvant, Army Liposomal Formulation adsorbed to aluminum (ALFA), in rhesus monkeys and observed 90% efficacy against limiting dose mucosal SHIV challenge in male animals. Here, we analyzed binding antibody responses, as assessed by Fc array profiling using a broad range of HIV-1 envelope antigens and Fc features, to explore the mechanisms of ALFA-mediated protection by employing machine learning and Cox proportional hazards regression analyses. We found that Fcγ receptor 2a-related binding antibody responses were augmented by ALFA relative to aluminium hydroxide, and these responses were associated with reduced risk of infection in male animals. Our results highlight the application of systems serology to provide mechanistic insights to vaccine-elicited protection and support evidence that antibody effector responses protect against HIV-1 infection.

Using the antibody-antigen binding interface to train image-based deep neural networks for antibody-epitope classification.

High-throughput B-cell sequencing has opened up new avenues for investigating complex mechanisms underlying our adaptive immune response. These technological advances drive data generation and the need to mine and analyze the information contained in these large datasets, in particular the identification of therapeutic antibodies (Abs) or those associated with disease exposure and protection. Here, we describe our efforts to use artificial intelligence (AI)-based image-analyses for prospective classification of Abs based solely on sequence information. We hypothesized that Abs recognizing the same part of an antigen share a limited set of features at the binding interface, and that the binding site regions of these Abs share share common structure and physicochemical property patterns that can serve as a "fingerprint" to recognize uncharacterized Abs. We combined large-scale sequence-based protein-structure predictions to generate ensembles of 3-D Ab models, reduced the Ab binding interface to a 2-D image (fingerprint), used pre-trained convolutional neural networks to extract features, and trained deep neural networks (DNNs) to classify Abs. We evaluated this approach using Ab sequences derived from human HIV and Ebola viral infections to differentiate between two Abs, Abs belonging to specific B-cell family lineages, and Abs with different epitope preferences. In addition, we explored a different type of DNN method to detect one class of Abs from a larger pool of Abs. Testing on Ab sets that had been kept aside during model training, we achieved average prediction accuracies ranging from 71-96% depending on the complexity of the classification task. The high level of accuracies reached during these classification tests suggests that the DNN models were able to learn a series of structural patterns shared by Abs belonging to the same class. The developed methodology provides a means to apply AI-based image recognition techniques to analyze high-throughput B-cell sequencing datasets (repertoires) for Ab classification.