Conformation-Dependent •OH/H2O2 Hydrogen Abstraction Reaction Cycles of Gly and Ala Residues: A Comparative Theoretical Study

CADMA-Chem: A Computational Protocol Based on Chemical Properties Aimed to Design Multifunctional Antioxidants

A computational protocol aimed to design new antioxidants with versatile behavior is presented. It is called Computer-Assisted Design of Multifunctional Antioxidants and is based on chemical properties (CADMA-Chem). The desired multi-functionality consists of in different methods of antioxidant protection combined with neuroprotection, although the protocol can also be used to pursue other health benefits. The dM38 melatonin derivative is used as a study case to illustrate the protocol in detail. This was found to be a highly promising candidate for the treatment of neurodegeneration, in particular Parkinson's and Alzheimer's diseases. This also has the desired properties of an oral-drug, which is significantly better than Trolox for scavenging free radicals, and has chelates redox metals, prevents the ●OH production, via Fenton-like reactions, repairs oxidative damage in biomolecules (lipids, proteins, and DNA), and acts as a polygenic neuroprotector by inhibiting catechol-O-methyl transferase (COMT), acetylcholinesterase (AChE) and monoamine oxidase B (MAOB). To the best of our best knowledge, CADMA-Chem is currently the only protocol that simultaneously involves the analyses of drug-like behavior, toxicity, manufacturability, versatile antioxidant protection, and receptor-ligand binding affinities. It is expected to provide a starting point that helps to accelerate the discovery of oral drugs with the potential to prevent, or slow down, multifactorial human health disorders.

Non-energetic, Low-Temperature Formation of Cα-Glycyl Radical, a Potential Interstellar Precursor of Natural Amino Acids

The reaction of H atoms with glycine was investigated at 3.1 K in para-H₂, a quantum-solid host. The reaction was followed by IR spectroscopy, with the spectral analysis aided by quantum chemical computations. Comparison of the experimental IR spectrum with computed anharmonic frequencies and intensities proved that, regardless of the reactant glycine conformation, C_α-glycyl radical is formed in an H-atom-abstraction process with great selectivity. The product of the second H-atom abstraction, iminoacetic acid, was also observed in a smaller amount. The C_α-glycyl radical is sensitive to UV light and decomposes to iminoacetic acid and H atom upon 280 nm radiation. Since the reactive radical center is located on the C_α-atom, it is suggested that natural α-amino acids can be formed from glycine via the C_α-glycyl radical by non-energetic mechanisms in the solid phase of the interstellar medium.

Insight into elemental mercury (Hg0) removal from flue gas using UV/H2O2 advanced oxidation processes

Elemental mercury (Hg⁰) emitted from coal-fired power plants and municipal solid waste (MSW) incinerators has caused great harm to the environment and human beings. The strong oxidized ^•OH radicals produced by UV/H₂O₂ advanced oxidation processes were studied to investigate the performance of Hg⁰ removal from simulated flue gases. The results showed that when H₂O₂ concentration was 1.0 mol/L and the solution pH value was 4.1, the UV/H₂O₂ system had the highest Hg⁰ removal efficiency. The optimal reaction temperature was approximately 50 °C and Hg⁰ removal was inhibited when the temperature was higher or lower. The yield of ^•OH radicals during UV/H₂O₂ reaction was studied by electron paramagnetic resonance (EPR) analysis. UV radiation was the determining factor to remove Hg⁰ in UV/H₂O₂ system due to ^•OH generation during H₂O₂ decomposition. SO₂ had little influence on Hg⁰ removal whereas NO had an inhibitory effect on Hg⁰ removal. The detailed findings for Hg⁰ removal reactions over UV/H₂O₂ make it an attractive method for mercury control from flue gases.

Comprehensive Investigation of the Antioxidant and Pro-oxidant Effects of Phenolic Compounds: A Double-Edged Sword in the Context of Oxidative Stress?

Oxidative stress (OS) is a health-threatening process that is involved, at least partially, in the development of several diseases. Although antioxidants can be used as a chemical defense against OS, they might also exhibit pro-oxidant effects, depending on environmental conditions. In this work, such a dual behavior was investigated for phenolic compounds (PhCs) within the framework of the density functional theory and based on kinetic data. Multiple reaction mechanisms were considered in both cases. The presence of redox metals, the pH, and the possibility that PhCs might be transformed into benzoquinones were identified as key aspects in the antioxidant versus pro-oxidant effects of these compounds. The main virtues of PhCs as antioxidants are their radical trapping activity, their regeneration under physiological conditions, and their behavior as OH-inactivating ligands. The main risks of PhCs as pro-oxidants are predicted to be the role of phenolate ions in the reduction of metal ions, which can promote Fenton-like reactions, and the formation of benzoquinones that might cause protein arylation at cysteine sites. Although the benefits seem to overcome the hazards, to properly design chemical strategies against OS using PhCs, it is highly recommended to carefully explore their duality in this context.

Oxidation of Acid, Base, and Amide Side-Chain Amino Acid Derivatives via Hydroxyl Radical

Hydroxyl radical (^•OH) is known to be highly reactive. Herein, we analyze the oxidation of acid (Asp and Glu), base (Arg and Lys), and amide (Asn and Gln) containing amino acid derivatives by the consecutive attack of two ^•OH. In this work, we study the reaction pathway by means of density functional theory. The oxidation mechanism is divided into two steps: (1) the first ^•OH can abstract a H atom or an electron, leading to a radical amino acid derivative, which is the intermediate of the reaction and (2) the second ^•OH can abstract another H atom or add itself to the formed radical, rendering the final oxidized products. The studied second attack of ^•OH is applicable to situations where high concentration of ^•OH is found, e.g., in vitro. Carbonyls are the best known oxidation products for these reactions. This work includes solvent dielectric and confirmation's effects of the reaction, showing that both are negligible. Overall, the most favored intermediates of the oxidation process at the side chain correspond to the secondary radicals stabilized by hyperconjugation. Intermediates show to be more stable in those cases where the spin density of the unpaired electron is lowered. Alcohols formed at the side chains are the most favored products, followed by the double-bond-containing ones. Interestingly, Arg and Lys side-chain scission leads to the most favored carbonyl-containing oxidation products, in line with experimental results.

Modelling the Effect of Conformation on Hydrogen-Atom Abstraction from Peptides

Computational quantum chemistry is used to examine the effect of conformation on the kinetics of hydrogen-atom abstraction by HO• from amides of glycine and proline as peptide models. In accord with previous findings, it is found that there are substantial variations possible in the conformations and the corresponding energies, with the captodative effect, hydrogen bonding, and solvation being some of the major features that contribute to the variations. The ‘minimum-energy-structure-pathway’ strategy that is often employed in theoretical studies of peptide chemistry with small models certainly provides valuable fundamental information. However, one may anticipate different reaction outcomes in structurally constrained systems due to modified reaction thermodynamics and kinetics, as demonstrated explicitly in the present study. Thus, using a ‘consistent-conformation-pathway’ approach may indeed be more informative in such circumstances, and in this regard theory provides information that would be difficult to obtain from experimental studies alone.

Protein Stability and Unfolding Following Glycine Radical Formation

Glycine (Gly) residues are particularly susceptible to hydrogen abstraction; which results in the formation of the capto-dative stabilized C_α-centered Gly radical (GLR) on the protein backbone. We examined the effect of GLR formation on the structure of the Trp cage; tryptophan zipper; and the villin headpiece; three fast-folding and stable miniproteins; using all-atom (OPLS-AA) molecular dynamics simulations. Radicalization changes the conformation of the GLR residue and affects both neighboring residues but did not affect the stability of the Trp zipper. The stability of helices away from the radical center in villin were also affected by radicalization; and GLR in place of Gly15 caused the Trp cage to unfold within 1 µs. These results provide new evidence on the destabilizing effects of protein oxidation by reactive oxygen species.

Chemical repair of protein carbon-centred radicals: long-distance dynamic factors

The thermodynamic and kinetic study of the repair reactions of three damaged aliphatic amino acids (alanine, valine, and leucine) with dihydrolipoic acid (DHLA) in a polar and a nonpolar solvent is presented in this work. Two simplified protein models were explored in the most common conformations (alpha helix and beta sheet). Calculations are performed at the M06-2X-SMD/6-31++G(d,p) level of theory. DHLA has shown to be an excellent antioxidant repair agent through hydrogen-transfer reaction involving the thiol groups, with rate constants close to diffusion control in most cases. The stability of the initial protein radical is not the most important factor determining the rate of the repair reaction because stabilizing intermolecular interactions involving the protein and the antioxidant can provide additional stability to some transition states accelerating the repair of sites that would otherwise not be so quickly repaired.

A computational study of radical initiated protein backbone homolytic dissociation on all natural amino acids

Hydroxyl radical (˙OH) is known to be one of the most reactive species. The attack of this radical onto the backbone of all natural amino acids is investigated.

The effect of oxidative stress on the bursopentin peptide structure: a theoretical study

Bursopentin (BP(5), H-Cys(1)-Lys(2)-Arg(3)-Val(4)-Tyr(5)-OH), found in the bursa Fabricius of the chicken, is a pentapeptide that protects the organism from oxidative stress by reducing the intracellular generation of reactive oxygen species. Hydrogen abstraction, a common oxidative reaction occurring in proteins, often results in the formation of d amino acid residues. To study the effect of this phenomenon on the structure of bursopentin, each of its residues were converted from the l configuration to the d configuration, and the structures of these peptide epimers were compared to that of the wild-type bursopentin. The conformations, secondary structures, compactness and hydrogen bonding of bursopentin were compared to its epimers using molecular dynamics simulations and first principles quantum chemical computations. It was discovered that the repulsion between the side chains of Lys(2) and Arg(3) influenced the conformation of the peptide regardless of the configuration of these residues. Epimerisation of the Val(4) and Tyr(5) caused a reduction in the compactness of bursopentin. In all cases, the occurrence of a turn structure was relatively high, especially when Arg(3) was in the d configuration. Thermodynamic analysis of the epimerisation process showed that the formation of d amino acid residues is favourable.

Reactivity and selectivity patterns in hydrogen atom transfer from amino acid C-H bonds to the cumyloxyl radical: polar effects as a rationale for the preferential reaction at proline residues

Absolute rate constants for hydrogen atom transfer (HAT) from the C-H bonds of N-Boc-protected amino acids to the cumyloxyl radical (CumO(•)) were measured by laser flash photolysis. With glycine, alanine, valine, norvaline, and tert-leucine, HAT occurs from the α-C-H bonds, and the stability of the α-carbon radical product plays a negligible role. With leucine, HAT from the α- and γ-C-H bonds was observed. The higher kH value measured for proline was explained in terms of polar effects, with HAT that predominantly occurs from the δ-C-H bonds, providing a rationale for the previous observation that proline residues represent favored HAT sites in the reactions of peptides and proteins with (•)OH. Preferential HAT from proline was also observed in the reactions of CumO(•) with the dipeptides N-BocProGlyOH and N-BocGlyGlyOH. The rate constants measured for CumO(•) were compared with the relative rates obtained previously for the corresponding reactions of different hydrogen-abstracting species. The behavior of CumO(•) falls between those observed for the highly reactive radicals Cl(•) and (•)OH and the significantly more stable Br(•). Taken together, these results provide a general framework for the description of the factors that govern reactivity and selectivity patterns in HAT reactions from amino acid C-H bonds.

Modelling the chemical repair of protein carbon-centered radicals formed via oxidative damage with dihydrolipoic acid

Dihydrolipoic acid repairs carbon-centred radicals at diffusion-controlled rates via HAT mechanism.

Reactions of the cumyloxyl radical with secondary amides. The influence of steric and stereoelectronic effects on the hydrogen atom transfer reactivity and selectivity

A time-resolved kinetic study of the hydrogen atom transfer (HAT) reactions from secondary alkanamides to the cumyloxyl radical was carried out in acetonitrile. HAT predominantly occurs from the N-alkyl α-C-H bonds, and a >60-fold decrease in kH was observed by increasing the steric hindrance of the acyl and N-alkyl groups. The role of steric and stereoelectronic effects on the reactivity and selectivity is discussed in the framework of HAT reactions from peptides.

Molecular ageing: free radical initiated epimerization of thymopentin--a case study

The epimerization of amino acid residues increases with age in living organisms. In the present study, the structural consequences and thermodynamic functions of the epimerization of thymopentin (TP-5), the active site of the thymic hormone thymopoietin, were studied using molecular dynamics and density functional theory methods. The results show that free radical-initiated D-amino acid formation is energetically favoured (-130 kJmol(-1)) for each residue and induces significant changes to the peptide structure. In comparison to the wild-type (each residue in the L-configuration), the radius of gyration of the D-Asp(3) epimer of the peptide decreased by 0.5 Å, and disrupted the intramolecular hydrogen bonding of the native peptide. Beyond establishing important structural, energetic and thermodynamic benchmarks and reference data for the structure of TP-5, these results disseminate the understanding of molecular ageing, the epimerization of amino acid residues.

Molecular mechanisms for the reaction between (˙)OH radicals and proline: insights on the role as reactive oxygen species scavenger in plant stress

The accumulation of proline (Pro) and overproduction of reactive oxygen species (ROS) by plants exposed to stress is well-documented. In vitro assays show that enzyme inactivation by hydroxyl radicals ((•)OH) can be avoided in the presence of Pro, suggesting this amino acid might act as a (•)OH scavenger. Although production of hydroxyproline (Hyp) has been hypothesized in connection with such antioxidant activity, no evidence on the detailed mechanism of scavenging has been reported. To elucidate whether and how Hyp might be produced, we used density functional theory calculations coupled to a polarizable continuum model to explore 27 reaction channels including H-abstraction by (•)OH and (•)OH/H2O addition. The structure and energetics of stable species and transition states for each reaction channel were characterized at the PCM-(U)M06/6-31G(d,p) level in aqueous solution. Evidence is found for a main pathway in which Pro scavenges (•)OH by successive H-abstractions (ΔG(‡,298) = 4.1 and 7.5 kcal mol(-1)) to yield 3,4-Δ-Pro. A companion pathway with low barriers yielding Δ(1)-pyrroline-5-carboxylate (P5C) is also supported, linking with 5-Hyp through hydration. However, this connection remains unlikely in stressed plants because P5C would be efficiently recycled to Pro (contributing to its accumulation) by P5C reductase, hypothesis coined here as the "Pro-Pro cycle".

Atropisomerism of the Asn α radicals revealed by Ramachandran surface topology

C radicals are typically trigonal planar and thus achiral, regardless of whether they originate from a chiral or an achiral C-atom (e.g., C-H + (•)OH → C• + H2O). Oxidative stress could initiate radical formation in proteins when, for example, the H-atom is abstracted from the Cα-carbon of an amino acid residue. Electronic structure calculations show that such a radical remains achiral when formed from the achiral Gly, or the chiral but small Ala residues. However, when longer side-chain containing proteogenic amino acid residues are studied (e.g., Asn), they provide radicals of axis chirality, which in turn leads to atropisomerism observed for the first time for peptides. The two enantiomeric extended backbone structures, •βL and •βD, interconvert via a pair of enantiotopic reaction paths, monitored on a 4D Ramachandran surface, with two distinct transition states of very different Gibbs-free energies: 37.4 and 67.7 kJ/mol, respectively. This discovery requires the reassessment of our understanding on radical formation and their conformational and stereochemical behavior. Furthermore, the atropisomerism of proteogenic amino acid residues should affect our understanding on radicals in biological systems and, thus, reframes the role of the D-residues as markers of molecular aging.

Computational study on the attack of ·OH radicals on aromatic amino acids

The attack of hydroxyl radicals on aromatic amino acid side chains, namely phenylalanine, tyrosine, and tryptophan, have been studied by using density functional theory. Two reaction mechanisms were considered: 1) Addition reactions onto the aromatic ring atoms and 2) hydrogen abstraction from all of the possible atoms on the side chains. The thermodynamics and kinetics of the attack of a maximum of two hydroxyl radicals were studied, considering the effect of different protein environments at two different dielectric values (4 and 80). The obtained theoretical results explain how the radical attacks take place and provide new insight into the reasons for the experimentally observed preferential mechanism. These results indicate that, even though the attack of the first (·)OH radical on an aliphatic C atom is energetically favored, the larger delocalization and concomitant stabilization that are obtained by attack on the aromatic side chain prevail. Thus, the obtained theoretical results are in agreement with the experimental evidence that the aromatic side chain is the main target for radical attack and show that the first (·)OH radical is added onto the aromatic ring, whereas a second radical abstracts a hydrogen atom from the same position to obtain the oxidized product. Moreover, the results indicate that the reaction can be favored in the buried region of the protein.

Reactions of the cumyloxyl and benzyloxyl radicals with strong hydrogen bond acceptors. Large enhancements in hydrogen abstraction reactivity determined by substrate/radical hydrogen bonding

A kinetic study on hydrogen abstraction from strong hydrogen bond acceptors such as DMSO, HMPA, and tributylphosphine oxide (TBPO) by the cumyloxyl (CumO(•)) and benzyloxyl (BnO(•)) radicals was carried out in acetonitrile. The reactions with CumO(•) were described in terms of a direct hydrogen abstraction mechanism, in line with the kinetic deuterium isotope effects, k(H)/k(D), of 2.0 and 3.1 measured for reaction of this radical with DMSO/DMSO-d(6) and HMPA/HMPA-d(18). Very large increases in reactivity were observed on going from CumO(•) to BnO(•), as evidenced by k(H)(BnO(•))/k(H)(CumO(•)) ratios of 86, 4.8 × 10(3), and 1.6 × 10(4) for the reactions with HMPA, TBPO, and DMSO, respectively. The k(H)/k(D) of 0.91 and 1.0 measured for the reactions of BnO(•) with DMSO/DMSO-d(6) and HMPA/HMPA-d(18), together with the k(H)(BnO(•))/k(H)(CumO(•)) ratios, were explained on the basis of the formation of a hydrogen-bonded prereaction complex between the benzyloxyl α-C-H and the oxygen atom of the substrates followed by hydrogen abstraction. This is supported by theoretical calculations that show the formation of relatively strong prereaction complexes. These observations confirm that in alkoxyl radical reactions specific hydrogen bond interactions can dramatically influence the hydrogen abstraction reactivity, pointing toward the important role played by structural and electronic effects.