MHC-Fine: Fine-tuned AlphaFold for precise MHC-peptide complex prediction.

The precise prediction of major histocompatibility complex (MHC)-peptide complex structures is pivotal for understanding cellular immune responses and advancing vaccine design. In this study, we enhanced AlphaFold's capabilities by fine-tuning it with a specialized dataset consisting of exclusively high-resolution class I MHC-peptide crystal structures. This tailored approach aimed to address the generalist nature of AlphaFold's original training, which, while broad-ranging, lacked the granularity necessary for the high-precision demands of class I MHC-peptide interaction prediction. A comparative analysis was conducted against the homology-modeling-based method Pandora as well as the AlphaFold multimer model. Our results demonstrate that our fine-tuned model outperforms others in terms of root-mean-square deviation (median value for Cα atoms for peptides is 0.66 Å) and also provides enhanced predicted local distance difference test scores, offering a more reliable assessment of the predicted structures. These advances have substantial implications for computational immunology, potentially accelerating the development of novel therapeutics and vaccines by providing a more precise computational lens through which to view MHC-peptide interactions.

Expanding FTMap for Fragment-Based Identification of Pharmacophore Regions in Ligand Binding Sites.

The knowledge of ligand binding hot spots and of the important interactions within such hot spots is crucial for the design of lead compounds in the early stages of structure-based drug discovery. The computational solvent mapping server FTMap can reliably identify binding hot spots as consensus clusters, free energy minima that bind a variety of organic probe molecules. However, in its current implementation, FTMap provides limited information on regions within the hot spots that tend to interact with specific pharmacophoric features of potential ligands. E-FTMap is a new server that expands on the original FTMap protocol. E-FTMap uses 119 organic probes, rather than the 16 in the original FTMap, to exhaustively map binding sites, and identifies pharmacophore features as atomic consensus sites where similar chemical groups bind. We validate E-FTMap against a set of 109 experimentally derived structures of fragment-lead pairs, finding that highly ranked pharmacophore features overlap with the corresponding atoms in both fragments and lead compounds. Additionally, comparisons of mapping results to ensembles of bound ligands reveal that pharmacophores generated with E-FTMap tend to sample highly conserved protein-ligand interactions. E-FTMap is available as a web server at https://eftmap.bu.edu.

Understanding the role of antibiotic-associated adverse events in influencing antibiotic decision-making.

Objective: To (1) understand the role of antibiotic-associated adverse events (ABX-AEs) on antibiotic decision-making, (2) understand clinician preferences for ABX-AE feedback, and (3) identify ABX-AEs of greatest clinical concern. Design: Focus groups. Setting: Academic medical center. Participants: Medical and surgical house staff, attending physicians, and advanced practice practitioners. Methods: Focus groups were conducted from May 2022 to December 2022. Participants discussed the role of ABX-AEs in antibiotic decision-making and feedback preferences and evaluated the prespecified categorization of ABX-AEs based on degree of clinical concern. Thematic analysis was conducted using inductive coding. Results: Four focus groups were conducted (n = 15). Six themes were identified. (1) ABX-AE risks during initial prescribing influence the antibiotic prescribed rather than the decision of whether to prescribe. (2) The occurrence of an ABX-AE leads to reassessment of the clinical indication for antibiotic therapy. (3) The impact of an ABX-AE on other management decisions is as important as the direct harm of the ABX-AE. (4) ABX-AEs may be overlooked because of limited feedback regarding the occurrence of ABX-AEs. (5) Clinicians are receptive to feedback regarding ABX-AEs but are concerned about it being punitive. (6) Feedback must be curated to prevent clinicians from being overwhelmed with data. Clinicians generally agreed with the prespecified categorizations of ABX-AEs by degree of clinical concern. Conclusions: The themes identified and assessment of ABX-AEs of greatest clinical concern may inform antibiotic stewardship initiatives that incorporate reporting of ABX-AEs as a strategy to reduce unnecessary antibiotic use.

Redox modulation of oxidatively-induced DNA damage by ascorbate enhances both in vitro and ex-vivo DNA damage formation and cell death in melanoma cells.

Elevated genomic instability in cancer cells suggests a possible model-scenario for their selective killing via the therapeutic delivery of well-defined levels of further DNA damage. To examine this scenario, this study investigated the potential for redox modulation of oxidatively-induced DNA damage by ascorbate in malignant melanoma (MM) cancer cells, to selectively enhance both DNA damage and MM cell killing. DNA damage was assessed by Comet and ɣH2AX assays, intracellular oxidising species by dichlorofluorescein fluorescence, a key antioxidant enzymatic defence by assessment of catalase activity and cell survival was determined by clonogenic assay. Comet revealed that MM cells had higher endogenous DNA damage levels than normal keratinocytes (HaCaT cells); this correlated MM cells having higher intracellular oxidising species and lower catalase activity, and ranked with MM cell melanin pigmentation. Comet also showed MM cells more sensitive towards the DNA damaging effects of exogenous H2O2, and that ascorbate further enhanced this H2O2-induced damage in MM cells; again, with MM cell sensitivity to induced damage ranking with degree of cell pigmentation. Furthermore, cell survival data indicated that ascorbate enhanced H2O2-induced clonogenic cell death selectively in MM cells whilst protecting HaCaT cells. Finally, we show that ascorbate serves to enhance the oxidising effects of the MM therapeutic drug Elesclomol in both established MM cells in vitro and primary cell cultures ex vivo. Together, these results suggest that ascorbate selectively enhances DNA damage and cell-killing in MM cells. This raises the option of incorporating ascorbate into clinical oxidative therapies to treat MM.

MHC-Fine: Fine-tuned AlphaFold for Precise MHC-Peptide Complex Prediction.

The precise prediction of Major Histocompatibility Complex (MHC)-peptide complex structures is pivotal for understanding cellular immune responses and advancing vaccine design. In this study, we enhanced AlphaFold's capabilities by fine-tuning it with a specialized dataset comprised by exclusively high-resolution MHC-peptide crystal structures. This tailored approach aimed to address the generalist nature of AlphaFold's original training, which, while broad-ranging, lacked the granularity necessary for the high-precision demands of MHC-peptide interaction prediction. A comparative analysis was conducted against the homology-modeling-based method Pandora [13], as well as the AlphaFold multimer model [8]. Our results demonstrate that our fine-tuned model outperforms both in terms of RMSD (median value is 0.65 Å) but also provides enhanced predicted lDDT scores, offering a more reliable assessment of the predicted structures. These advances have substantial implications for computational immunology, potentially accelerating the development of novel therapeutics and vaccines by providing a more precise computational lens through which to view MHC-peptide interactions.

An Electronic Algorithm to Identify Vancomycin-Associated Acute Kidney Injury.

Background: The burden of vancomycin-associated acute kidney injury (V-AKI) is unclear because it is not systematically monitored. The objective of this study was to develop and validate an electronic algorithm to identify cases of V-AKI and to determine its incidence. Methods: Adults and children admitted to 1 of 5 health system hospitals from January 2018 to December 2019 who received at least 1 dose of intravenous (IV) vancomycin were included. A subset of charts was reviewed using a V-AKI assessment framework to classify cases as unlikely, possible, or probable events. Based on review, an electronic algorithm was developed and then validated using another subset of charts. Percentage agreement and kappa coefficients were calculated. Sensitivity and specificity were determined at various cutoffs, using chart review as the reference standard. For courses ≥48 hours, the incidence of possible or probable V-AKI events was assessed. Results: The algorithm was developed using 494 cases and validated using 200 cases. The percentage agreement between the electronic algorithm and chart review was 92.5% and the weighted kappa was 0.95. The electronic algorithm was 89.7% sensitive and 98.2% specific in detecting possible or probable V-AKI events. For the 11 073 courses of ≥48 hours of vancomycin among 8963 patients, the incidence of possible or probable V-AKI events was 14.0%; the V-AKI incidence rate was 22.8 per 1000 days of IV vancomycin therapy. Conclusions: An electronic algorithm demonstrated substantial agreement with chart review and had excellent sensitivity and specificity in detecting possible or probable V-AKI events. The electronic algorithm may be useful for informing future interventions to reduce V-AKI.

The ClusPro AbEMap web server for the prediction of antibody epitopes.

Antibodies play an important role in the immune system by binding to molecules called antigens at their respective epitopes. These interfaces or epitopes are structural entities determined by the interactions between an antibody and an antigen, making them ideal systems to analyze by using docking programs. Since the advent of high-throughput antibody sequencing, the ability to perform epitope mapping using only the sequence of the antibody has become a high priority. ClusPro, a leading protein-protein docking server, together with its template-based modeling version, ClusPro-TBM, have been re-purposed to map epitopes for specific antibody-antigen interactions by using the Antibody Epitope Mapping server (AbEMap). ClusPro-AbEMap offers three different modes for users depending on the information available on the antibody as follows: (i) X-ray structure, (ii) computational/predicted model of the structure or (iii) only the amino acid sequence. The AbEMap server presents a likelihood score for each antigen residue of being part of the epitope. We provide detailed information on the server's capabilities for the three options and discuss how to obtain the best results. In light of the recent introduction of AlphaFold2 (AF2), we also show how one of the modes allows users to use their AF2-generated antibody models as input. The protocol describes the relative advantages of the server compared to other epitope-mapping tools, its limitations and potential areas of improvement. The server may take 45-90 min depending on the size of the proteins.

Comet Assay Profiling of FLASH-Induced Damage: Mechanistic Insights into the Effects of FLASH Irradiation.

Numerous studies have demonstrated the normal tissue-sparing effects of ultra-high dose rate 'FLASH' irradiation in vivo, with an associated reduction in damage burden being reported in vitro. Towards this, two key radiochemical mechanisms have been proposed: radical-radical recombination (RRR) and transient oxygen depletion (TOD), with both being proposed to lead to reduced levels of induced damage. Previously, we reported that FLASH induces lower levels of DNA strand break damage in whole-blood peripheral blood lymphocytes (WB-PBL) ex vivo, but our study failed to distinguish the mechanism(s) involved. A potential outcome of RRR is the formation of crosslink damage (particularly, if any organic radicals recombine), whilst a possible outcome of TOD is a more anoxic profile of induced damage resulting from FLASH. Therefore, the aim of the current study was to profile FLASH-induced damage via the Comet assay, assessing any DNA crosslink formation as a putative marker of RRR and/or anoxic DNA damage formation as an indicative marker of TOD, to determine the extent to which either mechanism contributes to the "FLASH effect". Following FLASH irradiation, we see no evidence of any crosslink formation; however, FLASH irradiation induces a more anoxic profile of induced damage, supporting the TOD mechanism. Furthermore, treatment of WB-PBLs pre-irradiation with BSO abrogates the reduced strand break damage burden mediated by FLASH exposures. In summary, we do not see any experimental evidence to support the RRR mechanism contributing to the reduced damage burden induced by FLASH. However, the observation of a greater anoxic profile of damage following FLASH irradiation, together with the BSO abrogation of the reduced strand break damage burden mediated by FLASH, lends further support to TOD being a driver of the reduced damage burden plus a change in the damage profile mediated by FLASH.

Streptomyces RNases - Function and impact on antibiotic synthesis.

Streptomyces are soil dwelling bacteria that are notable for their ability to sporulate and to produce antibiotics and other secondary metabolites. Antibiotic biosynthesis is controlled by a variety of complex regulatory networks, involving activators, repressors, signaling molecules and other regulatory elements. One group of enzymes that affects antibiotic synthesis in Streptomyces is the ribonucleases. In this review, the function of five ribonucleases, RNase E, RNase J, polynucleotide phosphorylase, RNase III and oligoribonuclease, and their impact on antibiotic production will be discussed. Mechanisms for the effects of RNase action on antibiotic synthesis are proposed.

Anemia prevalence time trends and disparities in the US population: examination of NHANES 1999-2020.

While a rising prevalence of anemia in the United States was reported in older studies, recent data are lacking. To estimate the prevalence and time trends of anemia in the United States and to examine how these estimates differ by gender, age, race, and household income to poverty threshold ratio (HIPR), we used the National Health and Nutrition Examination Surveys from 1999 to 2020. The presence of anemia was determined using the World Health Organization criteria. Survey-weighted raw and adjusted prevalence ratios (PRs) were determined using generalized linear models for the overall population and by gender, age, race, and HIPR. In addition, an interaction between gender and race was explored. Complete data on anemia, age, gender, and race were available on 87,554 participants (mean age = 34.6 years, women = 49.8%, Whites = 37.3%). Anemia prevalence increased from 4.03% during the 1999-2000 survey cycle to 6.49% during 2017-2020. In adjusted analyses, anemia prevalence was higher in >65 than in 26-45 years old (PR = 2.14, 95% confidence interval (CI) = 1.95, 2.35), in Blacks than Whites (PR = 3.97, 95% CI = 3.63, 4.35), in women than men (PR = 1.98, 95% CI = 1.83, 2.13), and in those with HIPR ≤ 1 than >4 (PR = 0.68, 95% CI = 0.61, 0.75). Gender modified the relationship between anemia and race; when compared to their male counterparts, Black, Hispanic, and other women had higher anemia prevalence than White women (all interaction p values <0.05). The anemia prevalence in the United States has risen from 1999 to 2020 and remains high among the elderly, minorities, and women. The difference in anemia prevalence between men and women is larger in non-Whites.

Calibration of the comet assay using ionising radiation.

Several trials have attempted to identify sources of inter-laboratory variability in comet assay results, aiming at achieving more equal responses. Ionising radiation induces a defined level of DNA single-strand breaks (per dose/base pairs) and is used as a reference when comparing comet results but relies on accurately determined radiation doses. In this ring test we studied the significance of dose calibrations and comet assay protocol differences, with the object of identifying causes of variability and how to deal with them. Eight participating laboratories, using either x-ray or gamma radiation units, measured dose rates using alanine pellet dosimeters that were subsequently sent to a specialised laboratory for analysis. We found substantial deviations between calibrated and nominal (uncalibrated) dose rates, with up to 46% difference comparing highest and lowest values. Three additional dosimetry systems were employed in some laboratories: thermoluminescence detectors and two aqueous chemical dosimeters. Fricke's and Benzoic Acid dosimetry solutions gave reliable quantitative dose estimations using local equipment. Mononuclear cells from fresh human blood or mammalian cell lines were irradiated locally with calibrated (alanine) radiation doses and analysed for DNA damage using a standardised comet assay protocol and a lab-specific protocol. The dose response of eight laboratories, calculated against calibrated radiation doses, was linear with slope variance CV= 29% with the lab-specific protocol, reduced to CV= 16% with the standard protocol. Variation between laboratories indicate post-irradiation repair differences. Intra-laboratory variation was very low judging from the dose response of 8 donors (CV=4%). Electrophoresis conditions were different in the lab-specific protocols explaining some dose response variations which were reduced by systematic corrections for electrophoresis conditions. The study shows that comet assay data obtained in different laboratories can be compared quantitatively using calibrated radiation doses and that systematic corrections for electrophoresis conditions are useful.

Mapping of antibody epitopes based on docking and homology modeling.

Antibodies are key proteins produced by the immune system to target pathogen proteins termed antigens via specific binding to surface regions called epitopes. Given an antigen and the sequence of an antibody the knowledge of the epitope is critical for the discovery and development of antibody based therapeutics. In this work, we present a computational protocol that uses template-based modeling and docking to predict epitope residues. This protocol is implemented in three major steps. First, a template-based modeling approach is used to build the antibody structures. We tested several options, including generation of models using AlphaFold2. Second, each antibody model is docked to the antigen using the fast Fourier transform (FFT) based docking program PIPER. Attention is given to optimally selecting the docking energy parameters depending on the input data. In particular, the van der Waals energy terms are reduced for modeled antibodies relative to x-ray structures. Finally, ranking of antigen surface residues is produced. The ranking relies on the docking results, that is, how often the residue appears in the docking poses' interface, and also on the energy favorability of the docking pose in question. The method, called PIPER-Map, has been tested on a widely used antibody-antigen docking benchmark. The results show that PIPER-Map improves upon the existing epitope prediction methods. An interesting observation is that epitope prediction accuracy starting from antibody sequence alone does not significantly differ from that of starting from unbound (i.e., separately crystallized) antibody structure.

Investigations of the novel checkpoint kinase 1 inhibitor SRA737 in non-small cell lung cancer and colorectal cancer cells of differing tumour protein 53 gene status.

Aim: In response to DNA damage the serine/threonine-specific protein kinase checkpoint kinase 1 (CHK1) is activated allowing cells to enter S phase (S) and G2 phase (G2) cell-cycle arrest. CHK1 inhibitors are expected to prevent cells from entering such arrest, thereby enhancing DNA damage-induced cytotoxicity. In contrast, normal cells with intact ataxia-telangiectasia mutated (ATM), CHK2 and tumour suppressor protein 53 (P53) signalling are still able to enter cell-cycle arrest using the functioning G1/S checkpoint, thereby being rescued from enhanced cytotoxicity. The main objective of this work is to investigate the in vitro effects of the novel CHK1 inhibitor SRA737 on pairs of non-small cell lung cancer (NSCLC) and colorectal cancer (CRC) cell lines, all with genetic aberrations rendering them susceptible to replication stress but of differing tumour protein 53 (TP53) gene status, focusing on DNA damage induction and the subsequent effects on cell proliferation and viability. Methods: NSCLC cell lines H23 [TP53 mutant (MUT)] and A549 [TP53 wild-type (WT)] and CRC cell lines HT29 (TP53 MUT) and HCT116 (TP53 WT) were incubated with differing micromolar concentrations of SRA737 for 24 h and then analysed using alkaline comet and phosphorylated H2A.X variant histone (γH2AX)-foci assays to assess mostly DNA single strand break and double strand break damage, respectively. Cell-counting/trypan blue staining was also performed to assess cell proliferation/viability. Results: Clear concentration-dependent increases in comet formation and γH2AX-foci/cell were noted for the TP53 MUT cells with no or lower increases being noted in the corresponding TP53 WT cells. Also, greater anti-proliferative and cell killing effects were noted in the TP53 MUT cells than in the TP53 WT cells. Conclusions: This study's data suggests that P53 status/functioning is a key factor in determining the sensitivity of NSCLC and CRC cancer cells towards CHK1 inhibition, even in circumstances conducive to high replicative stress.

API Development Increases Access to Shared Computing Resources at Boston University.

Within the last few decades, increases in computational resources have contributed enormously to the progress of science and engineering (S & E). To continue making rapid advancements, the S & E community must be able to access computing resources. One way to provide such resources is through High-Performance Computing (HPC) centers. Many academic research institutions offer their own HPC Centers but struggle to make the computing resources easily accessible and user-friendly. Here we present SHABU, a RESTful Web API framework that enables S & E communities to access resources from Boston University's Shared Computing Center (SCC). The SHABU requirements are derived from the use cases described in this work.

ATM Regulates Differentiation of Myofibroblastic Cancer-Associated Fibroblasts and Can Be Targeted to Overcome Immunotherapy Resistance.

Myofibroblastic cancer-associated fibroblast (myoCAF)-rich tumors generally contain few T cells and respond poorly to immune-checkpoint blockade. Although myoCAFs are associated with poor outcome in most solid tumors, the molecular mechanisms regulating myoCAF accumulation remain unclear, limiting the potential for therapeutic intervention. Here, we identify ataxia-telangiectasia mutated (ATM) as a central regulator of the myoCAF phenotype. Differentiating myofibroblasts in vitro and myoCAFs cultured ex vivo display activated ATM signaling, and targeting ATM genetically or pharmacologically could suppress and reverse differentiation. ATM activation was regulated by the reactive oxygen species-producing enzyme NOX4, both through DNA damage and increased oxidative stress. Targeting fibroblast ATM in vivo suppressed myoCAF-rich tumor growth, promoted intratumoral CD8 T-cell infiltration, and potentiated the response to anti-PD-1 blockade and antitumor vaccination. This work identifies a novel pathway regulating myoCAF differentiation and provides a rationale for using ATM inhibitors to overcome CAF-mediated immunotherapy resistance. SIGNIFICANCE: ATM signaling supports the differentiation of myoCAFs to suppress T-cell infiltration and antitumor immunity, supporting the potential clinical use of ATM inhibitors in combination with checkpoint inhibition in myoCAF-rich, immune-cold tumors.

To wait or not to wait: Optimal time interval between the first and second blood-culture sets to maximize blood-culture yield.

The optimal timing of blood culture (BCx) sets collection has not been evaluated with continuous BCx detection systems. The yield of BCx was similar between short intervals (median, 3 minutes) and longer intervals (median, 16 or 43 minutes) among 5,856 BCx, except for improved polymicrobial bacteremia detection with long-interval BCx.

Insights into the Bacterial Diversity and Detection of Opportunistic Pathogens in Mexican Chili Powder.

Chili powder is the most frequently consumed spice in Mexican diets. Thus, the dissemination of microorganisms associated with chili powder derived from Capsicum annuum L. is significant during microbial quality analysis, with special attention on detection of potential pathogens. The results presented here describe the initial characterization of bacterial community structure in commercial chili powder samples. Our results demonstrate that, within the domain Bacteria, the most abundant family was Bacillaceae, with a relative abundance of 99% in 71.4% of chili powder samples, while 28.6% of samples showed an average relative abundance of 60% for the Enterobacteriaceae family. Bacterial load for aerobic mesophilic bacteria (AMB) ranged from 104 to 106 cfu/g, while for sporulated mesophilic bacteria (SMB), the count ranged from 102 to 105 cfu/g. Bacillus cereus sensu lato (s.l.) was observed at ca. ˂600 cfu/g, while the count for Enterobacteriaceae ranged from 103 to 106 cfu/g, Escherichia coli and Salmonella were not detected. Fungal and yeast counts ranged from 102 to 105 cfu/g. Further analysis of the opportunistic pathogens isolated, such as B. cereus s.l. and Kosakonia cowanii, using antibiotic-resistance profiles and toxinogenic characteristics, revealed the presence of extended-spectrum ß-lactamases (ESBLs) and Metallo-ß-lactamases (MBLs) in these organisms. These results extend our knowledge of bacterial diversity and the presence of opportunistic pathogens associated with Mexican chili powder and highlight the potential health risks posed by its use through the spread of antibiotic-resistance and the production of various toxins. Our findings may be useful in developing procedures for microbial control during chili powder production.

Elucidation of protein function using computational docking and hotspot analysis by ClusPro and FTMap.

Starting with a crystal structure of a macromolecule, computational structural modeling can help to understand the associated biological processes, structure and function, as well as to reduce the number of further experiments required to characterize a given molecular entity. In the past decade, two classes of powerful automated tools for investigating the binding properties of proteins have been developed: the protein-protein docking program ClusPro and the FTMap and FTSite programs for protein hotspot identification. These methods have been widely used by the research community by means of publicly available online servers, and models built using these automated tools have been reported in a large number of publications. Importantly, additional experimental information can be leveraged to further improve the predictive power of these approaches. Here, an overview of the methods and their biological applications is provided together with a brief interpretation of the results.

FTMove: A Web Server for Detection and Analysis of Cryptic and Allosteric Binding Sites by Mapping Multiple Protein Structures.

Protein mapping distributes many copies of different molecular probes on the surface of a target protein in order to determine binding hot spots, regions that are highly preferable for ligand binding. While mapping of X-ray structures by the FTMap server is inherently static, this limitation can be overcome by the simultaneous analysis of multiple structures of the protein. FTMove is an automated web server that implements this approach. From the input of a target protein, by PDB code, the server identifies all structures of the protein available in the PDB, runs mapping on them, and combines the results to form binding hot spots and binding sites. The user may also upload their own protein structures, bypassing the PDB search for similar structures. Output of the server consists of the consensus binding sites and the individual mapping results for each structure - including the number of probes located in each binding site, for each structure. This level of detail allows the users to investigate how the strength of a binding site relates to the protein conformation, other binding sites, and the presence of ligands or mutations. In addition, the structures are clustered on the basis of their binding properties. The use of FTMove is demonstrated by application to 22 proteins with known allosteric binding sites; the orthosteric and allosteric binding sites were identified in all but one case, and the sites were typically ranked among the top five. The FTMove server is publicly available at https://ftmove.bu.edu.

FLASH irradiation induces lower levels of DNA damage ex vivo, an effect modulated by oxygen tension, dose, and dose rate.

OBJECTIVE: FLASH irradiation reportedly produces less normal tissue toxicity, while maintaining tumour response. To investigate oxygen's role in the 'FLASH effect', we assessed DNA damage levels following irradiation at different oxygen tensions, doses and dose rates. METHODS: Samples of whole blood were irradiated (20 Gy) at various oxygen tensions (0.25-21%) with 6 MeV electrons at dose rates of either 2 kGy/s (FLASH) or 0.1 Gy/s (CONV), and subsequently with various doses (0-40 Gy) and intermediate dose rates (0.3-1000 Gy/s). DNA damage of peripheral blood lymphocytes (PBL) were assessed by the alkaline comet assay. RESULTS: Following 20 Gy irradiation, lower levels of DNA damage were induced for FLASH, the difference being significant at 0.25% (p < 0.05) and 0.5% O2 (p < 0.01). The differential in DNA damage at 0.5% O2 was found to increase with total dose and dose rate, becoming significant for doses ≥20 Gy and dose rates ≥30 Gy/s. CONCLUSION: This study shows, using the alkaline comet assay, that lower levels of DNA damage are induced following FLASH irradiation, an effect that is modulated by the oxygen tension, and increases with the total dose and dose rate of irradiation, indicating that an oxygen related mechanism, e.g. transient radiation-induced oxygen depletion, may contribute to the tissue sparing effect of FLASH irradiation. ADVANCES IN KNOWLEDGE: This paper is first to directly show that FLASH-induced DNA damage is modulated by oxygen tension, total dose and dose rate, with FLASH inducing significantly lower levels of DNA damage for doses ≥20 Gy and dose rates ≥30 Gy/s, at 0.5% O2.