Evaluation of Hydroxyl Radical Reactivity by Thioether Group Proximity in Model Peptide Backbone: Methionine versus S-Methyl-Cysteine.

Hydroxyl radicals (HO•) have long been regarded as a major source of cellular damage. The reaction of HO• with methionine residues (Met) in peptides and proteins is a complex multistep process. Although the reaction mechanism has been intensively studied, some essential parts remain unsolved. In the present study we examined the reaction of HO• generated by ionizing radiation in aqueous solutions under anoxic conditions with two compounds representing the simplest model peptide backbone CH3C(O)NHCHXC(O)NHCH3, where X = CH2CH2SCH3 or CH2SCH3, i.e., the Met derivative in comparison with the cysteine-methylated derivative. We performed the identification and quantification of transient species by pulse radiolysis and final products by LC-MS and high-resolution MS/MS after γ-radiolysis. The results allowed us to draw for each compound a mechanistic scheme. The fate of the initial one-electron oxidation at the sulfur atom depends on its distance from the peptide backbone and involves transient species of five-membered and/or six-membered ring formations with different heteroatoms present in the backbone as well as quite different rates of deprotonation in forming α-(alkylthio)alkyl radicals.

Photo- and Radiation-Induced One-Electron Oxidation of Methionine in Various Structural Environments Studied by Time-Resolved Techniques.

Oxidation of methionine (Met) is an important reaction that plays a key role in protein modifications during oxidative stress and aging. The first steps of Met oxidation involve the creation of very reactive and short-lived transients. Application of complementary time-resolved radiation and photochemical techniques (pulse radiolysis and laser flash photolysis together with time-resolved CIDNP and ESR techniques) allowed comparing in detail the one-electron oxidation mechanisms initiated either by ●OH radicals and other one-electron oxidants or the excited triplet state of the sensitizers e.g., 4-,3-carboxybenzophenones. The main purpose of this review is to present various factors that influence the character of the forming intermediates. They are divided into two parts: those inextricably related to the structures of molecules containing Met and those related to external factors. The former include (i) the protection of terminal amine and carboxyl groups, (ii) the location of Met in the peptide molecule, (iii) the character of neighboring amino acid other than Met, (iv) the character of the peptide chain (open vs cyclic), (v) the number of Met residues in peptide and protein, and (vi) the optical isomerism of Met residues. External factors include the type of the oxidant, pH, and concentration of Met-containing compounds in the reaction environment. Particular attention is given to the neighboring group participation, which is an essential parameter controlling one-electron oxidation of Met. Mechanistic aspects of oxidation processes by various one-electron oxidants in various structural and pH environments are summarized and discussed. The importance of these studies for understanding oxidation of Met in real biological systems is also addressed.

Modulation of the Directionality of Hole Transfer between the Base and the Sugar-Phosphate Backbone in DNA with the Number of Sulfur Atoms in the Phosphate Group.

This work shows that S atom substitution in phosphate controls the directionality of hole transfer processes between the base and sugar-phosphate backbone in DNA systems. The investigation combines synthesis, electron spin resonance (ESR) studies in supercooled homogeneous solution, pulse radiolysis in aqueous solution at ambient temperature, and density functional theory (DFT) calculations of in-house synthesized model compound dimethylphosphorothioate (DMTP(O-)═S) and nucleotide (5'-O-methoxyphosphorothioyl-2'-deoxyguanosine (G-P(O-)═S)). ESR investigations show that DMTP(O-)═S reacts with Cl2•- to form the σ2σ*1 adduct radical -P-S[Formula: see text]Cl, which subsequently reacts with DMTP(O-)═S to produce [-P-S[Formula: see text]S-P-]-. -P-S[Formula: see text]Cl in G-P(O-)═S undergoes hole transfer to Gua, forming the cation radical (G•+) via thermally activated hopping. However, pulse radiolysis measurements show that DMTP(O-)═S forms the thiyl radical (-P-S•) by one-electron oxidation, which did not produce [-P-S[Formula: see text]S-P-]-. Gua in G-P(O-)═S is oxidized unimolecularly by the -P-S• intermediate in the sub-picosecond range. DFT thermochemical calculations explain the differences in ESR and pulse radiolysis results obtained at different temperatures.

Pulse Radiolysis Studies of Temperature Dependent Electron Transfers among Redox Centers in ba3-Cytochrome c Oxidase from Thermus thermophilus: Comparison of A- and B-Type Enzymes.

The functioning of cytochrome c oxidases involves orchestration of long-range electron transfer (ET) events among the four redox active metal centers. We report the temperature dependence of electron transfer from the CuAr site to the low-spin heme-(a)bo site, i.e., CuAr + heme-a(b)o → CuAo + heme-a(b)r in three structurally characterized enzymes: A-type aa3 from Paracoccus denitrificans (PDB code 3HB3) and bovine heart tissue (PDB code 2ZXW), and the B-type ba3 from T. thermophilus (PDB codes 1EHK and 1XME). k,T data sets were obtained with the use of pulse radiolysis as described previously. Semiclassical Marcus theory revealed that λ varies from 0.74 to 1.1 eV, Hab, varies from ∼2 × 10-5 eV (0.16 cm-1) to ∼24 × 10-5 eV (1.9 cm-1), and ßD varies from 9.3 to 13.9. These parameters are consistent with diabatic electron tunneling. The II-Asp111Asn CuA mutation in cytochrome ba3 had no effect on the rate of this reaction whereas the II-Met160Leu CuA-mutation was slower by an amount corresponding to a decreased driving force of ∼0.06 eV. The structures support the presence of a common, electron-conducting "wire" between CuA and heme-a(b). The transfer of an electron from the low-spin heme to the high-spin heme, i.e., heme-a(b)r + heme-a3o → heme-a(b)o + heme-a3r, was not observed with the A-type enzymes in our experiments but was observed with the Thermus ba3; its Marcus parameters are λ = 1.5 eV, Hab = 26.6 × 10-5 eV (2.14 cm-1), and ßD = 9.35, consistent also with diabatic electron tunneling between the two hemes. The II-Glu15Ala mutation of the K-channel structure, ∼ 24 Å between its CA and Fe-a3, was found to completely block heme-br to heme-a3o electron transfer. A structural mechanism is suggested to explain these observations.

Solvated Electron in Acetonitrile: Radiation Yield, Absorption Spectrum, and Equilibrium between Cavity- and Solvent-Localized States.

The equilibrium between a solvent cavity-localized electron, ecav-, and a dimeric solvent anion, (CH3CN)2•-, which are the two lowest energy states of the solvated electron in acetonitrile, has been investigated by pulse radiolysis at 233-353 K. The enthalpy and entropy for the ecav- to (CH3CN)2•- conversion amount to -11.2 ± 0.3 kcal/mol and -39.3 ± 1.2 cal/(mol K), corresponding to a 0.44 ± 0.35 equilibrium constant at 25 °C. The radiation yield of the solvated electron has been quantified using a Co(II) macrocycle that scavenges electrons with a 1.55 × 1011 M-1 s-1 rate constant. The apparent yield increases without saturation over the attainable scavenger concentration range, reaching 2.8 per 100 eV; this value represents the lower limit for the acetonitrile ionization yield in pulse radiolysis. The apparent molar absorption coefficient of (20.8 ± 1.5) × 103 M-1 cm-1 at 1450 nm and 20 °C for the solvated electron and individual vis-near-infrared (NIR) absorption spectra of ecav- and (CH3CN)2•- are derived from the data. Variances with previous reports are thoroughly discussed. Collectively, these results resolve several controversies concerning the solvated electron properties in acetonitrile and furnish requisite data for quantitative pulse radiolysis investigations in this commonly used solvent.

Geant4-DNA Modeling of Water Radiolysis beyond the Microsecond: An On-Lattice Stochastic Approach.

In this work, we use the next sub-volume method (NSM) to investigate the possibility of using the compartment-based ("on-lattice") model to simulate water radiolysis. We, first, start with a brief description of the reaction-diffusion master equation (RDME) in a spatially discretized simulation volume ("mesh"), which is divided into sub-volumes (or "voxels"). We then discuss the choice of voxel size and merging technique of a given mesh, along with the evolution of the system using the hierarchical algorithm for the RDME ("hRDME"). Since the compartment-based model cannot describe high concentration species of early radiation-induced spurs, we propose a combination of the particle-based step-by-step ("SBS") Brownian dynamics model and the compartment-based model ("SBS-RDME model") for the simulation. We, finally, use the particle-based SBS Brownian dynamics model of Geant4-DNA as a reference to test the model implementation through several benchmarks. We find that the compartment-based model can efficiently simulate the system with a large number of species and for longer timescales, beyond the microsecond, with a reasonable computing time. Our aim in developing this model is to study the production and evolution of reactive oxygen species generated under irradiation with different dose rate conditions, such as in FLASH and conventional radiotherapy.

Dissociations of free radicals to generate protons, electrophiles or nucleophiles: role in DNA strand breaks.

The concept behind the research described in this article was that of marrying the 'soft' methods of radical generation with the effectiveness and flexibility of nucleophile/electrophile synthetic procedures. Classic studies with pulse radiolysis and laser flash photolysis had shown that free radicals could be more acidic than their closed shell counterparts. QM computations harmonised with this and helped to define which radical centres and which structural types were most effective. Radicals based on the sulfonic acid moiety and on the Meldrum's acid moiety (2,2-dimethyl-1,3-dioxane-4,6-dione) were found to be extreme examples in the superacid class. The ethyne unit could be used as a very effective spacer between the radical centre and the site of proton donation. The key factor in promoting acidity was understood to be the thermodynamic stabilisation of the conjugate anion-radicals released on deprotonation. Solvation played a key part in promoting this and theoretical microhydration studies provided notable support. A corollary was that heterolytic dissociations of free radicals to yield either electrophiles or nucleophiles were also enhanced relative to non-radical models. The most effective radical types for spontaneous releases of both these types of reagents were identified. Ethyne units were again effective as spacers. The enhancement of release of phosphate anions by adjacent radical centres was an important special case. Reactive oxygen species and also diradicals from endiyne antibiotics generate C4'-deoxyribose radicals from nucleotides. Radicals of these types spontaneously release phosphate and triphosphate and this is a contributor to DNA and RNA strand breaks.

Radiation- and Photo-Induced Oxidation Pathways of Methionine in Model Peptide Backbone under Anoxic Conditions.

Within the reactive oxygen species (ROS) generated by cellular metabolisms, hydroxyl radicals (HO•) play an important role, being the most aggressive towards biomolecules. The reactions of HO• with methionine residues (Met) in peptides and proteins have been intensively studied, but some fundamental aspects remain unsolved. In the present study we examined the biomimetic model made of Ac-Met-OMe, as the simplest model peptide backbone, and of HO• generated by ionizing radiation in aqueous solutions under anoxic conditions. We performed the identification and quantification of transient species by pulse radiolysis and of final products by LC-MS and high-resolution MS/MS after γ-radiolysis. By parallel photochemical experiments, using 3-carboxybenzophenone (CB) triplet with the model peptide, we compared the outcomes in terms of short-lived intermediates and stable product identification. The result is a detailed mechanistic scheme of Met oxidation by HO•, and by CB triplets allowed for assigning transient species to the pathways of products formation.

First-Principles Simulations of Biological Molecules Subjected to Ionizing Radiation.

Ionizing rays cause damage to genomes, proteins, and signaling pathways that normally regulate cell activity, with harmful consequences such as accelerated aging, tumors, and cancers but also with beneficial effects in the context of radiotherapies. While the great pace of research in the twentieth century led to the identification of the molecular mechanisms for chemical lesions on the building blocks of biomacromolecules, the last two decades have brought renewed questions, for example, regarding the formation of clustered damage or the rich chemistry involving the secondary electrons produced by radiolysis. Radiation chemistry is now meeting attosecond science, providing extraordinary opportunities to unravel the very first stages of biological matter radiolysis. This review provides an overview of the recent progress made in this direction, focusing mainly on the atto- to femto- to picosecond timescales. We review promising applications of time-dependent density functional theory in this context.

Impact of high-energy electron irradiation on mechanical, structural and chemical properties of agarose hydrogels.

Due to their excellent biocompatibility and biodegradability, natural hydrogels are highly demanded biomaterials for biomedical applications such as wound dressing, tissue engineering, drug delivery or three dimensional cell culture. Highly energetic electron irradiation up to 10 MeV is a powerful and fast tool to sterilize and tailor the material's properties. In this study, electron radiation treatment of agarose hydrogels was investigated to evaluate radiation effects on physical, structural and chemical properties. The viscoelastic behavior, surface hydrophilicity and swelling behavior in a range of typical sterilization doses of 0 kGy to 30 kGy was analyzed. The mechanical properties were determined by rheology measurements and decreased by more than 20% compared to the initial moduli. The number average molecular weight between crosslinks was estimated based on rubber elasticity theory to judge on the radiation degradation. In this dose range, the number average molecular weight between crosslinks increased by more than 6%. Chemical structure was investigated by FTIR spectroscopy to evaluate the radiation resistance of agarose hydrogels. With increasing electron dose, an increasing amount of carbonyl containing species was observed. In addition, irradiation was accompanied by formation of gas cavities in the hydrogels. The gas products were specified for CO2, CO and H2O. Based on the radiolytic products, a radiolysis mechanism was proposed. Electron beam treatment under high pressure conditions was found to reduce gas cavity formation in the hydrogels.

Fast autooxidation of a bis-histidyl-ligated globin from the anhydrobiotic tardigrade, Ramazzottius varieornatus, by molecular oxygen.

Tardigrades, a phylum of meiofaunal organisms, exhibit extraordinary tolerance to various environmental conditions, including extreme temperatures (-273 to 151°C) and exposure to ionizing radiation. Proteins from anhydrobiotic tardigrades with homology to known proteins from other organisms are new potential targets for structural genomics. Recently, we reported spectroscopic and structural characterization of a hexacoordinated haemoglobin (Kumaglobin [Kgb]) found in an anhydrobiotic tardigrade. In the absence of its exogenous ligand, Kgb displays hexacoordination with distal and proximal histidines. In this work, we analysed binding of the molecular oxygen ligand following reduction of haem in Kgb using a pulse radiolysis technique. Radiolytically generated hydrated electrons (eaq-) reduced the haem iron of Kgb within 20 µs. Subsequently, ferrous haem reacted with O2 to form a ferrous-dioxygen intermediate with a second-order rate constant of 3.0 × 106 M-1 s-1. The intermediate was rapidly (within 0.1 s) autooxidized to the ferric form. Redox potential measurements revealed an E'0 of -400 mV (vs. standard hydrogen electrode) in the ferric/ferrous couple. Our results suggest that Kgb may serve as a physiological generator of O2▪- via redox signalling and/or electron transfer.

OH radical reactions with the hydrophilic component of sphingolipids.

In this work, using the example of model compounds, we studied the reactions resulting from the interaction of OH radicals with the hydrophilic part of sphingolipids. We compared the stopped-flow EPR spectroscopy and pulse radiolysis with optical detection methods to characterize radical intermediates formed in the reaction of OH radicals with glycerol, serinol and N-boc-serinol. Quantum chemical calculations were also performed to help interpret the observed experimental data. It was shown that H-abstraction from the terminal carbon atom is the main process that is realized for all the studied compounds. The presence of the unsubstituted amino group (-NH2) is seen to completely change the reaction properties of serinol in comparison with those observed in glycerol and N-boc serinol.

Radiation-Induced Oxidation Reactions of 2-Selenouracil in Aqueous Solutions: Comparison with Sulfur Analog of Uracil.

One-electron oxidation of 2-selenouracil (2-SeU) by hydroxyl (●OH) and azide (●N3) radicals leads to various primary reactive intermediates. Their optical absorption spectra and kinetic characteristics were studied by pulse radiolysis with UV-vis spectrophotometric and conductivity detection and by the density functional theory (DFT) method. The transient absorption spectra recorded in the reactions of ●OH with 2-SeU are dominated by an absorption band with an λmax = 440 nm, the intensity of which depends on the concentration of 2-SeU and pH. Based on the combination of conductometric and DFT studies, the transient absorption band observed both at low and high concentrations of 2-SeU was assigned to the dimeric 2c-3e Se-Se-bonded radical in neutral form (2●). The dimeric radical (2●) is formed in the reaction of a selenyl-type radical (6●) with 2-SeU, and both radicals are in equilibrium with Keq = 1.3 × 104 M-1 at pH 4 (below the pKa of 2-SeU). Similar equilibrium with Keq = 4.4 × 103 M-1 was determined for pH 10 (above the pKa of 2-SeU), which admittedly involves the same radical (6●) but with a dimeric 2c-3e Se-Se bonded radical in anionic form (2●-). In turn, at the lowest concentration of 2-SeU (0.05 mM) and pH 10, the transient absorption spectrum is dominated by an absorption band with an λmax = 390 nm, which was assigned to the ●OH adduct to the double bond at C5 carbon atom (3●) based on DFT calculations. Similar spectral and kinetic features were also observed during the ●N3-induced oxidation of 2-SeU. In principle, our results mostly revealed similarities in one-electron oxidation pathways of 2-SeU and 2-thiouracil (2-TU). The major difference concerns the stability of dimeric radicals with a 2c-3e chalcogen-chalcogen bond in favor of 2-SeU.

Gold Nanoparticle Formation via X-ray Radiolysis Investigated with Time-Resolved X-ray Liquidography.

We report the generation of gold nanoparticles (AuNPs) from the aqueous solution of chloro(2,2',2″-terpyridine)gold(III) ion ([Au(tpy)Cl]2+) through X-ray radiolysis and optical excitation at a synchrotron. The original purpose of the experiment was to investigate the photoinduced structural changes of [Au(tpy)Cl]2+ upon 400 nm excitation using time-resolved X-ray liquidography (TRXL). Initially, the TRXL data did not show any signal that would suggest structural changes of the solute molecule, but after an induction time, the TRXL data started to show sharp peaks and valleys. In the early phase, AuNPs with two types of morphology, dendrites, and spheres, were formed by the reducing action of hydrated electrons generated by the X-ray radiolysis of water, thereby allowing the detection of TRXL data due to the laser-induced lattice expansion and relaxation of AuNPs. Along with the lattice expansion, the dendritic and spherical AuNPs were transformed into smaller, raspberry-shaped AuNPs of a relatively uniform size via ablation by the optical femtosecond laser pulse used for the TRXL experiment. Density functional theory calculations confirm that the reduction potential of the metal complex relative to the hydration potential of X-ray-generated electrons determines the facile AuNP formation observed for [Au(tpy)Cl]2+.

Ultrafast Dynamics of Water Radiolysis: Hydrated Electron Formation, Solvation, Recombination, and Scavenging.

The ultrafast formation, solvation, and geminate recombination of hydrated electrons upon vacuum ultraviolet photoexcitation of liquid water and the static and dynamic scavenging by NO3- are investigated using femtosecond time-resolved photoelectron spectroscopy. The solvation time constant for excess electrons is typical of that for liquid water but increases slightly with increasing excitation energy. The electron survival probability for geminate recombination is found to be much lower than the literature values owing to previously unobserved ultrafast geminate recombination in a period of 5 ps. NO3- induces the ultrafast (static) scavenging of photoexcited electronic states of liquid water and the dynamic scavenging of detached electrons with a reaction rate that is dependent on the excitation energy. The formation of hydrated electrons at 7.7 eV is ascribed to a H-atom-transfer process, but it is plausible that additional formation channels open at higher energies.

Modelling Spatial Scales of Dose Deposition and Radiolysis Products from Gold Nanoparticle Sensitisation of Proton Therapy in A Cell: From Intracellular Structures to Adjacent Cells.

Gold nanoparticle (GNP) enhanced proton therapy is a promising treatment concept offering increased therapeutic effect. It has been demonstrated in experiments which provided indications that reactive species play a major role. Simulations of the radiolysis yield from GNPs within a cell model were performed using the Geant4 toolkit. The effect of GNP cluster size, distribution and number, cell and nuclear membrane absorption and intercellular yields were evaluated. It was found that clusters distributed near the nucleus increased the nucleus yield by 91% while reducing the cytoplasm yield by 7% relative to a disperse distribution. Smaller cluster sizes increased the yield, 200 nm clusters had nucleus and cytoplasm yields 117% and 35% greater than 500 nm clusters. Nuclear membrane absorption reduced the cytoplasm and nucleus yields by 8% and 35% respectively to a permeable membrane. Intercellular enhancement was negligible. Smaller GNP clusters delivered near sub-cellular targets maximise radiosensitisation. Nuclear membrane absorption reduces the nucleus yield, but can damage the membrane providing another potential pathway for biological effect. The minimal effect on adjacent cells demonstrates that GNPs provide a targeted enhancement for proton therapy, only effecting cells with GNPs internalised. The provided quantitative data will aid further experiments and clinical trials.

Ultra-high dose-rate (FLASH) radiotherapy: Generation of early, transient, strongly acidic spikes in the irradiated tumor environment.

Monte Carlo simulations of γ/fast electron-radiolysis of water show that the in situ formation of H3O+ temporarily renders each "native" isolated spur/track region very acidic. For pulsed (FLASH) irradiation with high dose rate, this early time, transient "acid-spike" response is shown to extend evenly across the entire irradiated volume. Since pH controls many cellular processes, this study highlights the need to consider these spikes of acidity in understanding the fundamental mechanisms underlying FLASH radiotherapy.

The physical stage of radiolysis of solvated DNA by high-energy-transfer particles: insights from new first principles simulations.

The primary processes that occur following direct irradiation of bio-macromolecules by ionizing radiation determine the multiscale responses that lead to biomolecular lesions. The so-called physical stage loosely describes processes of energy deposition and molecular ionization/excitation but remains largely elusive. We propose a new approach based on first principles density functional theory to simulate energy deposition in large and heterogeneous biomolecules by high-energy-transfer particles. Unlike traditional Monte Carlo approaches, our methodology does not rely on pre-parametrized sets of cross-sections, but captures excitation, ionization and low energy electron emission at the heart of complex biostructures. It furthermore gives access to valuable insights on ultrafast charge and hole dynamics on the femtosecond time scale. With this new tool, we reveal the mechanisms of ionization by swift ions in microscopic DNA models and solvated DNA comprising almost 750 atoms treated at the DFT level of description. We reveal a so-called ebb-and-flow ionization mechanism in which polarization of the irradiated moieties appears as a key feature. We also investigate where secondary electrons produced by irradiation localize on chemical moieties composing DNA. We compare irradiation of solvated DNA by light (H+, and He2+) vs. heavier (C6+) ions, highlighting the much higher probability of double ionization with the latter. Our methodology constitutes a stepping stone towards a greater understanding of the chemical stage and more generally towards the multiscale modelling of radiation damage in biology using first principles.

Hydrated electrons induce the formation of interstrand cross-links in DNA modified by cisplatin adducts.

Double-stranded oligonucleotides containing cisplatin adducts, with and without a mismatched region, were exposed to hydrated electrons generated by gamma-rays. Gel electrophoresis analysis demonstrates the formation of cisplatin-interstrand crosslinks from the cisplatin-intrastrand species. The rate constant per base for the reaction between hydrated electrons and the double-stranded oligonucleotides with and without cisplatin containing a mismatched region was determined by pulse radiolysis to be 7 × 109 and 2 × 109 M-1 s-1, respectively. These results provide a better understanding of the radiosensitizing effect of cisplatin adducts in hypoxic tumors and of the formation of interstrand crosslinks, which are difficult for cells to repair.