Pith. sign in

REVIEW 5 major objections 6 minor 47 references

Phenylalanine modification in plasma-driven biocatalysis revealed by solvent accessibility and reactive dynamics in combination with protein mass spectrometry

T0 review · 5 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read A solvent-accessible-surface screen predicts which amino acid residues plasma species modify, and mass spectrometry confirms three phenylalanine oxidations on CviUPO.

desk verdict Useful screening workflow with an honest limitations section, but the central validation claim is not supported: the MS search was seeded by the prediction, and the residue-level top hits from SASA were not the experimentally modified sites. read the letter →

arxiv 2506.20205 v1 pith:T22UNPH7 submitted 2025-06-25 physics.bio-ph

classification physics.bio-ph
keywords solventaccessiblesurfaceareaplasma-generatedspeciesreactivemoleculardynamicsphenylalanineoxidationunspecificperoxygenaseplasma-drivenbiocatalysisproteinmodificationscreeningmassspectrometry
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Plasma-driven biocatalysis uses a discharge to supply H2O2 to enzymes, but the same plasma generates radicals that can oxidize the enzyme. This paper argues that a solvent-accessible-surface-area (SASA) interaction screen, combined with reactive molecular dynamics, can identify amino acid residues most likely to be modified by plasma-generated species. Applied to the peroxygenase CviUPO, the screen flagged phenylalanine and basic residues as OH interaction partners, and mass spectrometry confirmed oxidation of Phe51, Phe67, and Phe208 after plasma treatment, together with two methionines that the screen had not singled out. The paper's conclusion is that this SASA-based screen is a fast method for directing targeted mass spectrometry searches and for locating regions of a protein surface that small reactive molecules will attack.

What carries the argument

The machinery is a two-stage computational pipeline. First, a SASA interaction analysis places a probe molecule at each solvent-accessible point about 1.4 Å from the protein surface, rotates it to the best orientation, and computes a local interaction energy using the ReaxFF reactive force field, a bond-order potential that allows bonds to break and form during the simulation; post-processing ranks residues by total and frequency-normalized (relative) interactions, using the thirty strongest energy points. This yields a three-dimensional interaction map for each plasma species. Second, reactive molecular dynamics validate the hot spots: short 75 fs simulations at the ten lowest-energy positions test whether a covalent bond forms, and longer high-concentration simulations with randomly placed species check whether the residue preferences survive diffusion and competition with solvent. The pipeline was run on three enzymes, CviUPO, AaeUPO, and GapA, and experimentally checked on CviUPO.

What would settle it

Plasma-treat CviUPO with and without a hydroxyl-radical scavenger and compare oxidation of Phe51, Phe67, and Phe208 by mass spectrometry; if oxidation persists when OH radicals are scavenged, the claim that OH drives these modifications is wrong. A complementary check would compute the OH–phenylalanine reaction energy with a high-level quantum-chemical method on a model compound and compare it with the ReaxFF value used in the screening.

Watch

Extended reading notes

Core claim

The central claim is that local interaction energies computed at solvent-accessible points on a protein surface identify the residues most likely to be chemically modified by reactive plasma species. The method places one model plasma species at a time on surface points 1.4 Å from the protein, relaxes its orientation, and computes a local interaction energy $E_{\mathrm{int}} = E_{\mathrm{SASA}} - E_{\mathrm{macro\,mol}} - E_{\mathrm{probe\,mol}}$ from the ReaxFF reactive force field. Hot spots cluster on specific residues; reactive molecular dynamics at the ten lowest-energy sites shows that H, O, and OH form bonds there, and high-concentration simulations reproduce the same residue preferences even with random starting positions. For CviUPO, the screen predicted strong interactions with surface-exposed phenylalanines, and dielectric-barrier-discharge plasma treatment followed by tryptic digestion and LC-MS/MS detected significant oxidation of Phe51, Phe67, and Phe208 in all replicates. The paper therefore claims that SASA-based screening can direct targeted mass spectrometry searches and, more generally, locate protein-surface regions of interest for small reactive molecules.

Load-bearing premise

The load-bearing premise is that the reactive force field faithfully represents the chemistry of plasma species such as hydroxyl radicals, and the paper's own supporting information concedes that distinguishing OH radicals from OH ions in ReaxFF 'remains highly uncertain'; if that chemistry is misrepresented, the predicted interaction profiles and modification sites may not reflect real plasma conditions.

Editorial extensions

If this is right

  • If the method is right, any enzyme can be screened before plasma experiments to obtain a short list of residues to search for modification in MS data, reducing the search space for time-consuming proteomics.
  • For CviUPO, the oxidation of Phe51 and Phe67, both near the substrate channel, flags surface-exposed aromatic residues as a potential cause of activity loss under plasma treatment.
  • The solvent simulations say most reactive species react with water before reaching the enzyme; this directly supports immobilization or geometry changes that lengthen the diffusion path as protection strategies.
  • At elevated temperatures, unfolding exposes buried residues and the heme cofactor to attack, so thermal control of the reaction environment should limit damage to catalytically essential regions.
  • Because the interaction energy definition is not specific to plasma species, the same SASA screen could be used to locate binding or modification sites for other small reactive molecules.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable extension beyond this paper: benchmark the SASA screen on enzymes with already published plasma- or oxidant-induced modification maps, asking whether the rank order of predicted hot spots matches experimentally observed modified residues.
  • Editorial inference: the screen's residue preferences may be driven mostly by surface exposure and electrostatics rather than by accurate radical chemistry; comparing SASA hot spots against a surface-exposure-only ranking would show how much the reactive dynamics add.
  • If Phe oxidation near the substrate channel indeed impairs catalysis, engineering solvent-exposed phenylalanines into less oxidizable residues could improve enzyme stability under plasma-driven conditions.
  • The current MS search covered only a small panel of modifications, so the actual damage may be broader; a wider search including nitration and carbonyl formation could test whether other SASA-flagged residues, such as lysine and arginine, are modified under conditions that favor reactive nitrogen species.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 6 minor

Summary. The paper proposes a computational workflow to predict which amino-acid residues of enzymes are modified by plasma-generated species (PGS) such as H2O2, OH, O2, O, NO, and H. The workflow combines a SASA-based interaction analysis, in which a single probe species is placed at solvent-accessible surface points and ReaxFF interaction energies are computed (Eq. 1), with reactive molecular dynamics (ReaxFF with the Monti et al. force field) as a short-time validation of bond formation at predicted positions. The method is applied to three enzymes, with CviUPO as the main experimental subject. For CviUPO, the analysis predicts strong OH interactions with Lys, Phe, and Arg; a DBD plasma-treated sample is then analyzed by LC-MS/MS, and the authors report significant oxidation of Phe51, Phe67, Phe208, Met42, and Met220, claiming that the predicted modification of phenylalanine was experimentally confirmed. The Conclusion and the Supplementary Information explicitly acknowledge that the Phe prediction prompted the inclusion of Phe modifications in the MS database search and that ReaxFF cannot reliably distinguish OH radicals from OH ions.

Significance. If validated, the proposed screening workflow would be practically useful for the plasma-biocatalysis community as a cheap computational prescreen for ranking amino-acid residues on a protein surface with respect to reactivity toward plasma-generated species, and for directing time-consuming mass spectrometry searches. The manuscript has notable strengths: the SASA-analysis package is released on GitHub, the simulation and mass spectrometry data are provided, and the SI candidly documents the ReaxFF species limitation. However, the current evidence does not establish the central claim as stated. The MS validation is partly circular because the search space was seeded by the prediction; Phe oxidation is a known general plasma effect; the agreement is at the level of residue type rather than residue identity; and the force-field ambiguity acknowledged in the SI affects the entire interaction-energy ranking. The paper currently supports a screening heuristic with a plausible consistency check, not an independently confirmed predictive method.

major comments (5)
  1. [Conclusion] The Conclusion states that 'the Phe modification by plasma treatment has thus been brought into focus by the SASA analysis, which prompted the inclusion of Phe modifications in the MS database search.' Because the prediction seeded the database search, the subsequent detection of Phe oxidation is not an independent confirmation of the method; it is a consistency check whose search space was conditioned on the prediction. This circularity is compounded by the fact that Phe hydroxylation is a well-documented plasma-induced modification of amino acids in the very literature cited by the authors (refs. 43 and 44), so its observation would be expected even if the SASA screen carried no predictive information. The Abstract's claim that predictions were 'confirmed experimentally by mass spectrometry' should be weakened to a consistency statement, or the authors should provide a discriminating validation, e.g., a prospective or pre-registered MS search over a fixed modification list in which the Phe prediction was made before the search was executed.
  2. [Results, SASA interaction analysis (CviUPO); Supplementary Table 1] The experimental match is at the amino-acid-type level, not the residue level. The ten strongest OH interaction points at 300 K listed in Supplementary Table 1 map to Leu225, Asp91, Arg84, Val153, and other non-Phe residues, and neither the short reactive-MD validation nor the high-concentration simulations (Fig. 5) report bond formation at Phe51, Phe67, or Phe208; they report Phe only as a residue type. The SASA evidence for Phe comes from the frequency-normalized counts in Fig. 2(e), i.e., Phe is over-represented in the pool of strong-interaction points, not from residue-specific interaction energies at the three positions later observed by MS. The paper should explicitly state that the method predicts which amino-acid types are preferential targets and that residue-level correspondence between prediction and experiment is not currently demonstrated.
  3. [SI, 'Explanation to the species nature in ReaxFF'] The SI states that distinguishing OH radicals from OH ions in ReaxFF 'remains highly uncertain' and that 'any conclusions drawn from these simulations should be interpreted with caution,' and the main text acknowledges that the modeled 'rmdPGS ... are not real plasma-generated species.' Since the SASA interaction energies (Eq. 1) and all bond-formation statistics in the reactive MD are computed with this force field, the mechanistic attribution of Phe oxidation to OH and O in the Conclusion is not supported by a species-resolved model. This caveat should be moved into the main text and the mechanistic claims reframed as applying to the modeled species only. A concrete strengthening would be to benchmark the force field against an experimentally known relative reactivity order, e.g., the high susceptibility of Met and Cys reported in refs. 43-45, before using it to rationalize Phe-selective modification.
  4. [Results, Mass spectrometry of plasma-treated CviUPO] The screen produced a false negative for a substantial fraction of the experimentally observed modifications: the SASA analysis reports no interactions with Met or Cys because these residues are rare and buried, yet Met42 and Met220 are oxidized in all three replicates. The paper's explanation that buried residues become accessible post hoc does not rescue the predictive claim, since the same argument could apply to many buried residues and is not testable from the reported data. The authors should quantify the screen's miss rate (two of the five observed modified residues belong to a residue type the screen excluded) and discuss whether the Met oxidation indicates that the static SASA screen misses modifications that require protein dynamics, radical diffusion into the protein interior, or species not represented in the model. The weak experimental signal for the strongly predicted Lys/Arg modifications (detected only in some replicates and at low relative intensity) further limits what the current experiment can establish about the method's selectivity.
  5. [Table 1] The mass spectrometry table reports modified Phe51, Phe67, and Phe208, with Phe51 and Phe67 located in the same 33-residue tryptic peptide, but no site-localization evidence (e.g., fragment-ion coverage or a localization score) is provided, and no error bars or significance tests are given for the relative intensities. Because the headline claim is residue-specific oxidation, the authors should report the localization confidence for each modified site (e.g., the number of site-determining fragment ions) and descriptive statistics over the three replicates, and define the 'relative intensity > 100' acceptance criterion.
minor comments (6)
  1. [Figure 5] The figure caption is internally inconsistent: it is titled 'Bond analysis for high concentrations of hydrogen' while the surrounding text and the axis description refer to OH, and the caption itself mentions 'an additional OH atom' while the panels are described in terms of hydrogen. Please reconcile the caption with the plotted species.
  2. [References] References 39 and 43 appear to cite the same article by Takai et al. (2014) with different journal names (Journal of Physics D: Applied Physics vs. Plasma Sources Science and Technology); please deduplicate and correct.
  3. [Computational details] In Eqs. (2) and (3), the quantity 'Interactions Total Number' is not defined; please specify whether it is the total number of SASA points, the number of points above the top-30 cutoff, or the summed interaction counts over the analyzed set.
  4. [Computational details] The default cutoff of the thirty strongest interaction energies is described only as 'determined empirically' with no supporting data; a brief sensitivity check (e.g., top 10 vs. top 50) would help the reader judge the robustness of the residue rankings and is recommended given the load this parameter carries.
  5. [Table 1] The acceptance criterion 'relative intensity > 100' is undocumented; please state how the relative intensities were computed (e.g., normalization to peptide or protein abundance) and why the threshold of 100 was chosen.
  6. [Throughout] Several typos and LaTeX artifacts remain in the text, including 'occurring in textitCviUPO', the broken 'N V Ten' ensemble label, and the phrase 'are display in Figure 5'; a careful proofread is needed.

Circularity Check

1 steps flagged · score 4.0 of 10

Partial circularity: the MS 'confirmation' of Phe modification is not independent because the SASA prediction was used to seed the MS database search, although unanticipated Met oxidation and weak Lys/Arg signals prevent a fully closed loop.

  1. self definitional [Conclusion (last paragraph); compare 'Mass spectrometry of plasma-treated CviUPO' section]
    "The Phe modification by plasma treatment has thus been brought into focus by the SASA analysis, which prompted the inclusion of Phe modifications in the MS database search."

    The paper's central validation loop is: SASA analysis predicts strong Phe/Lys/Arg interactions, the MS search then includes Phe, Arg, and Lys modifications 'because of their predicted strong interaction with ROS such as OH,' and the same experiment finds Phe+O and is presented as confirmation. This is not an independent test: the modification being tested was added to the search because the prediction said to look for it, so the positive Phe hit is an expected consequence of the search design rather than evidence that the SASA screen adds predictive value.

full rationale

The main circularity is the seeded MS validation: the Phe hypothesis was generated by the SASA analysis and then the MS search was explicitly designed to include Phe modifications, so detecting Phe+O does not independently validate the method. This is a genuine partial circularity. However, the loop is not closed: Met oxidation was found although SASA predicted no Met/Cys interactions, and predicted Lys/Arg modifications were largely not observed, giving the experiment independent content. The simulation-internal checks (short ReaxFF MD at the ten lowest SASA-energy points, and high-concentration MD comparisons) are consistency checks using the same force field; they do not independently validate the SASA approximation but are not circular in the harmful sense. No load-bearing self-citation is present: the ReaxFF force field is external (Monti et al. 2013), and the authors explicitly note that they did not develop it. The SI caveat that ReaxFF 'lacks the precision required to unambiguously differentiate OH radicals from OH ions' and that conclusions 'should be interpreted with caution' weakens the computational prediction but is an honest limitation, not circularity. The known background that plasma readily oxidizes phenylalanine (refs. 43, 44) further weakens the evidential value of the Phe hit, but that is a correctness/novelty concern, not a circularity one. Overall score 4 reflects partial circularity while acknowledging independent elements in the experimental outcome.

Assumptions & free parameters 1 free parameters · 3 assumptions · 0 invented entities

The central claim rests on the adequacy of ReaxFF for plasma species chemistry, on the assumption that vacuum SASA interaction energies predict modification propensity, and on the completeness of the targeted MS search. No new physical entities are introduced. The only hand-tuned numeric choice flagged is the 30-point cutoff in SASA post-processing.

free parameters (1)
  • SASA post-processing cutoff = 30 strongest interaction energies
    The post-processing functions only take into account the thirty strongest interaction energies, a cutoff determined empirically rather than from a principled criterion. This choice directly affects which residues are reported as top targets.
assumptions (3)
  • domain assumption ReaxFF with the Monti et al. 2013 biomolecule force field adequately represents the reactive chemistry of plasma-generated species.
    The entire simulation pipeline relies on this force field. The SI states that distinguishing OH radicals from OH ions is 'highly uncertain' and that conclusions should be interpreted with caution.
  • domain assumption SASA interaction energies computed in vacuum with a single static probe at solvent-accessible points correlate with covalent modification propensity in solution.
    The method places one probe at each SASA point and computes interaction energy in vacuum; the link between this energy and actual bond formation is assumed and tested only for a few positions in short MD runs.
  • domain assumption The targeted mass spectrometry search covers the modifications that matter for validation.
    The MS search included only a predefined set of modifications and achieved 70% sequence coverage, so some plasma-induced modifications could have been missed, which would bias the validation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Phenylalanine modification in plasma-driven biocatalysis revealed by solvent accessibility and reactive dynamics in combination with protein mass spectrometry." pith.science (2026). https://pith.science/paper/T22UNPH7

@misc{pith2026250620205,
  author       = {Pith},
  title        = {Pith review of: Phenylalanine modification in plasma-driven biocatalysis revealed by solvent accessibility and reactive dynamics in combination with protein mass spectrometry},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/T22UNPH7}},
  note         = {Machine review of arXiv:2506.20205}
}
read the original abstract

Biocatalysis is an emerging field that provides an environmentally friendly alternative to conventional catalysis, but still it faces some challenges. One of the major difficulties for biocatalysts that require reactive species like H2O2 as co-substrates lies in the concentration of these reactive species. On the one hand, they are used as reactants, but on the other hand, they inactivate the enzymes at high concentrations. When utilizing non-thermal plasma to deliver H2O2 for biocatalysis, it is essential to understand the potential interactions between plasma-generated species (PGS) and enzymes. This is particularly important because, alongside \ch{H2O2}, other reactive species such as hydroxyl radicals, atomic oxygen, superoxide, and nitric oxide are also produced. The investigation of the localized reactivity of the solvent accessible surface area (SASA) of an enzyme, with certain species, is an important tool for predicting these interactions. In combination with reactive molecular dynamics (MD) simulations this enabled us to identify amino acid residues that are likely targets for modifications by the PGS. A subset of the theoretical predictions made in the present study were confirmed experimentally by mass spectrometry, underlining the utility of the SASA and MD based screening approach to direct time-consuming experiments and assist their interpretation.

Figures

Figures reproduced from arXiv: 2506.20205 by the authors.

Figure 1
Figure 1. 3D structural representations of the enzymes [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Interaction profiles from SASA for the enzyme [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. 3D representations of CviUPO (left) and SASA interaction profile (right) for OH at different temperatures: (a) at 300 K, (b) 400 K, and (c) at 500 K. The SASA procedure and analysis are consistent with those, used to generate the previous figure. The 3D repre￾sentations were generated using PyMOL, 11 with input structures derived from the output of the equilibration simulations at the specified temperatures. For a d… view at source ↗
Figures from the paper (20 more)
Figure 4
Figure 4. Figure 4: Comparison of the SASA interaction map of OH at 300 K between (a) [PITH_FULL_IMAGE:figures/full_fig_p017_4.png]
Figure 5
Figure 5. Figure 5: Bond analysis for high concentrations of hydrogen with [PITH_FULL_IMAGE:figures/full_fig_p021_5.png]
Figure 1
Figure 1. Figure 1: Amino acid sequence of the Collariella virescens UPO used in the experimental study. The poly-His tail sequence is shown in italics and thrombin recognition sequence is underlined. The amino acid positions without the poly-His tail sequence (utilized in the text) are s…
Figure 2
Figure 2. Figure 2: Interaction profiles from SASA for the enzyme [PITH_FULL_IMAGE:figures/full_fig_p030_2.png]
Figure 3
Figure 3. Figure 3: SASA interaction profile for H and OH with [PITH_FULL_IMAGE:figures/full_fig_p031_3.png]
Figure 4
Figure 4. Figure 4: SASA interaction profile for H2O2 and NO with CviUPO at different temperatures: (a) at 300 K, (b) 400 K, and (c) at 500 K. 32 [PITH_FULL_IMAGE:figures/full_fig_p032_4.png]
Figure 5
Figure 5. Figure 5: SASA interaction profile for O2 and O with CviUPO at different temperatures: (a) at 300 K, (b) 400 K, and (c) at 500 K. 33 [PITH_FULL_IMAGE:figures/full_fig_p033_5.png]
Figure 6
Figure 6. Figure 6: Interaction profiles from SASA for the enzyme [PITH_FULL_IMAGE:figures/full_fig_p036_6.png]
Figure 7
Figure 7. Figure 7: SASA total interactions per residue at 300 K for [PITH_FULL_IMAGE:figures/full_fig_p037_7.png]
Figure 8
Figure 8. Figure 8: Interaction profiles from SASA for the enzyme GapA at 300 K. (a) shows the [PITH_FULL_IMAGE:figures/full_fig_p039_8.png]
Figure 9
Figure 9. Figure 9: SASA total interactions per residue at 300 K for GapA. (a) shows the interactions [PITH_FULL_IMAGE:figures/full_fig_p040_9.png]
Figure 10
Figure 10. Figure 10: Bond analysis for high concentrations of H with [PITH_FULL_IMAGE:figures/full_fig_p042_10.png]
Figure 11
Figure 11. Figure 11: Bond analysis for high concentrations of H [PITH_FULL_IMAGE:figures/full_fig_p043_11.png]
Figure 12
Figure 12. Figure 12: Bond analysis for high concentrations of NO with [PITH_FULL_IMAGE:figures/full_fig_p044_12.png]
Figure 13
Figure 13. Figure 13: Bond analysis for high concentrations of O [PITH_FULL_IMAGE:figures/full_fig_p045_13.png]
Figure 14
Figure 14. Figure 14: Bond analysis for high concentrations of oxygen with [PITH_FULL_IMAGE:figures/full_fig_p046_14.png]
Figure 15
Figure 15. Figure 15: Bond analysis for high concentrations of PGS with [PITH_FULL_IMAGE:figures/full_fig_p048_15.png]
Figure 16
Figure 16. Figure 16: Bond analysis for high concentrations of PGS with [PITH_FULL_IMAGE:figures/full_fig_p049_16.png]
Figure 17
Figure 17. Figure 17: Bond analysis for high concentrations of PGS with GapA at 300 K. All panels [PITH_FULL_IMAGE:figures/full_fig_p051_17.png]
Figure 18
Figure 18. Figure 18: Bond analysis for high concentrations of PGS with GapA at 300 K. All panels [PITH_FULL_IMAGE:figures/full_fig_p052_18.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

47 extracted references · 46 canonical work pages

  1. [1]

    Catalysis: Concepts and Green Applications, 2nd ed.; Wiley-VCH Verlag, 2017

    Rothenberg, G. Catalysis: Concepts and Green Applications, 2nd ed.; Wiley-VCH Verlag, 2017

  2. [2]

    M.; Tymoczko, J

    Berg, J. M.; Tymoczko, J. L.; Gatto, G. J.; Stryer, L. Biochemie; Springer Berlin Heidelberg, 2018

  3. [3]

    L.; Finnigan, W.; France, S

    Bell, E. L.; Finnigan, W.; France, S. P.; Green, A. P.; Hayes, M. A.; Hepworth, L. J.; Lovelock, S. L.; Niikura, H.; Osuna, S.; Romero, E.; Ryan, K. S.; Turner, N. J.; Flitsch, S. L. Biocatalysis. Nature Reviews Methods Primers 2021, 1, 1--21

  4. [4]

    C.; Baldenius, K.; Bornscheuer, U

    Wu, S.; Snajdrova, R.; Moore, J. C.; Baldenius, K.; Bornscheuer, U. T. Biocatalysis: Enzymatic Synthesis for Industrial Applications. Angewandte Chemie International Edition 2021, 60, 88--119

  5. [5]

    T.; Weiland, R

    Vasudevan, P. T.; Weiland, R. H. Deactivation of catalase by hydrogen peroxide. Biotechnology and Bioengineering 1990, 36, 783--789

  6. [6]

    Suicide Inactivation of Peroxidases and the Challenge of Engineering More Robust Enzymes

    Valderrama, B.; Ayala, M.; Vazquez-Duhalt, R. Suicide Inactivation of Peroxidases and the Challenge of Engineering More Robust Enzymes . Chemistry & Biology 2002, 9, 555--565

  7. [7]

    Exploring the catalase activity of unspecific peroxygenases and the mechanism of peroxide-dependent heme destruction

    Karich, A.; Scheibner, K.; Ullrich, R.; Hofrichter, M. Exploring the catalase activity of unspecific peroxygenases and the mechanism of peroxide-dependent heme destruction. Journal of Molecular Catalysis B: Enzymatic 2016, 134, 238--246

  8. [8]

    L.; Grunden, A

    Wapshott-Stehli, H. L.; Grunden, A. M. In situ H2O2 generation methods in the context of enzyme biocatalysis. Enzyme and Microbial Technology 2021, 145, 109744

Show all 47 references
  1. [9]

    F.; Hollmann, F.; Alcalde, M.; Kourist, R.; Bandow, J

    Yayci, A.; Álvaro Gómez Baraibar; Krewing, M.; Fueyo, E. F.; Hollmann, F.; Alcalde, M.; Kourist, R.; Bandow, J. E. Plasma-Driven in Situ Production of Hydrogen Peroxide for Biocatalysis. ChemSusChem 2020, 13, 2072--2079

  2. [10]

    Yayci, A.; Dirks, T.; Kogelheide, F.; Alcalde, M.; Hollmann, F.; Awakowicz, P.; Bandow, J. E. Microscale Atmospheric Pressure Plasma Jet as a Source for Plasma-Driven Biocatalysis. ChemCatChem 2020, 12, 5893--5897

  3. [11]

    I.; Walsh, J

    Dickenson, A.; Britun, N.; Nikiforov, A.; Leys, C.; Hasan, M. I.; Walsh, J. L. The generation and transport of reactive nitrogen species from a low temperature atmospheric pressure air plasma source. Physical Chemistry Chemical Physics 2018, 20, 28499--28510

  4. [12]

    2022; https://zenodo.org/records/6637571

    GROMACS 2022.2 Source code. 2022; https://zenodo.org/records/6637571

  5. [13]

    J.; Murtola, T.; Schulz, R.; P \' a ll, S.; Smith, J

    Abraham, M. J.; Murtola, T.; Schulz, R.; P \' a ll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS : High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1-2, 19--25

  6. [14]

    Berendsen, H. J. C.; van der Spoel, D.; van Drunen, R. GROMACS : A message-passing parallel molecular dynamics implementation. Computer Physics Communications 1995, 91, 43--56

  7. [15]

    GROMACS 4: 0.167em Algorithms for Highly Efficient , Load - Balanced , and Scalable Molecular Simulation

    Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E. GROMACS 4: 0.167em Algorithms for Highly Efficient , Load - Balanced , and Scalable Molecular Simulation . Journal of Chemical Theory and Computation 2008, 4, 435--447

  8. [16]

    GROMACS 3.0: a package for molecular simulation and trajectory analysis

    Lindahl, E.; Hess, B.; van der Spoel, D. GROMACS 3.0: a package for molecular simulation and trajectory analysis. Molecular modeling annual 2001, 7, 306--317

  9. [17]

    J.; Kutzner, C.; Hess, B.; Lindahl, E

    P \' a ll, S.; Abraham, M. J.; Kutzner, C.; Hess, B.; Lindahl, E. Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS . Solving Software Challenges for Exascale . Cham, 2015; pp 3--27

  10. [18]

    Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS

    P \' a ll, S.; Zhmurov, A.; Bauer, P.; Abraham, M.; Lundborg, M.; Gray, A.; Hess, B.; Lindahl, E. Heterogeneous parallelization and acceleration of molecular dynamics simulations in GROMACS . The Journal of Chemical Physics 2020, 153, 134110

  11. [19]

    R.; Smith, J

    Pronk, S.; P \' a ll, S.; Schulz, R.; Larsson, P.; Bjelkmar, P.; Apostolov, R.; Shirts, M. R.; Smith, J. C.; Kasson, P. M.; van der Spoel, D.; Hess, B.; Lindahl, E. GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit. Bioinformatics 2013...

  12. [20]

    E.; Berendsen, H

    Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A. E.; Berendsen, H. J. C. GROMACS : Fast , flexible, and free. Journal of Computational Chemistry 2005, 26, 1701--1718

  13. [21]

    A.; Armacost, K

    Soteras Guti \' e rrez, I.; Lin, F.-Y.; Vanommeslaeghe, K.; Lemkul, J. A.; Armacost, K. A.; Brooks, C. L.; MacKerell, A. D. Parametrization of halogen bonds in the CHARMM general force field: Improved treatment of ligand--protein interactions. Bioorganic & Medicinal Chemistry ...

  14. [22]

    Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; Mackerell Jr., A. D. CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force ...

  15. [23]

    Vanommeslaeghe, K.; MacKerell, A. D. J. Automation of the CHARMM General Force Field ( CGenFF ) I : Bond Perception and Atom Typing . Journal of Chemical Information and Modeling 2012, 52, 3144--3154

  16. [24]

    P.; MacKerell, A

    Vanommeslaeghe, K.; Raman, E. P.; MacKerell, A. D. J. Automation of the CHARMM General Force Field ( CGenFF ) II : Assignment of Bonded Parameters and Partial Atomic Charges . Journal of Chemical Information and Modeling 2012, 52, 3155--3168

  17. [25]

    Yu, W.; He, X.; Vanommeslaeghe, K.; MacKerell Jr., A. D. Extension of the CHARMM general force field to sulfonyl-containing compounds and its utility in biomolecular simulations. Journal of Computational Chemistry 2012, 33, 2451--2468

  18. [26]

    From Proteins to Perturbed Hamiltonians: A Suite of Tutorials for the GROMACS -2018 Molecular Simulation Package

    Lemkul, J. From Proteins to Perturbed Hamiltonians: A Suite of Tutorials for the GROMACS -2018 Molecular Simulation Package. Living Journal of Computational Molecular Science 2019, 1

  19. [27]

    2023; https://zenodo.org/records/10806852

    LAMMPS : Large -scale Atomic / Molecular Massively Parallel Simulator . 2023; https://zenodo.org/records/10806852

  20. [28]

    M.; Fogarty, J

    Aktulga, H. M.; Fogarty, J. C.; Pandit, S. A.; Grama, A. Y. Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques. Parallel Computing 2012, 38, 245--259

  21. [29]

    Fast Parallel Algorithms for Short - Range Molecular Dynamics

    Plimpton, S. Fast Parallel Algorithms for Short - Range Molecular Dynamics . Journal of Computational Physics 1995, 117, 1--19

  22. [30]

    P.; Aktulga, H

    Thompson, A. P.; Aktulga, H. M.; Berger, R.; Bolintineanu, D. S.; Brown, W. M.; Crozier, P. S.; in 't Veld, P. J.; Kohlmeyer, A.; Moore, S. G.; Nguyen, T. D.; Shan, R.; Stevens, M. J.; Tranchida, J.; Trott, C.; Plimpton, S. J. LAMMPS - a flexible simulation tool for particle-b...

  23. [31]

    Duin, A. C. V.; Dasgupta, S.; Lorant, F.; Goddard, W. A. ReaxFF: A Reactive Force Field for Hydrocarbons. Journal of Physical Chemistry A 2001, 105, 9396--9409

  24. [32]

    L.; Shin, Y

    Monti, S.; Corozzi, A.; Fristrup, P.; Joshi, K. L.; Shin, Y. K.; Oelschlaeger, P.; Duin, A. C. V.; Barone, V. Exploring the conformational and reactive dynamics of biomolecules in solution using an extended version of the glycine reactive force field. Physical Chemistry Chemic...

  25. [33]

    VMD -- V isual M olecular D ynamics

    Humphrey, W.; Dalke, A.; Schulten, K. VMD -- V isual M olecular D ynamics. Journal of Molecular Graphics 1996, 14, 33--38

  26. [34]

    P.; Wright, W

    Varshney, A.; Brooks, F. P.; Wright, W. V. Linearly Scalable Computation of Smooth Molecular Surfaces. IEEE Computer Graphics and Applications 1994, 14, 19--25

  27. [35]

    Visualization and analysis of atomistic simulation data with OVITO –the Open Visualization Tool

    Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO –the Open Visualization Tool . Modelling and Simulation in Materials Science and Engineering 2009, 18, 015012

  28. [36]

    Studier, F. W. Protein production by auto-induction in high-density shaking cultures. Protein Expression and Purification 2005, 41, 207--234

  29. [37]

    Reaction of OH with phenylalanine in neutral aqueous solution

    Solar, S. Reaction of OH with phenylalanine in neutral aqueous solution. Radiation Physics and Chemistry (1977) 1985, 26, 103--108

  30. [38]

    Chemical modification of amino acids by atmospheric-pressure cold plasma in aqueous solution

    Takai, E.; Kitamura, T.; Kuwabara, J.; Ikawa, S.; Yoshizawa, S.; Shiraki, K.; Kawasaki, H.; Arakawa, R.; Kitano, K. Chemical modification of amino acids by atmospheric-pressure cold plasma in aqueous solution. Journal of Physics D: Applied Physics 2014, 47, 285403

  31. [39]

    Reactive molecular dynamics simulations on interaction mechanisms of cold atmospheric plasmas and peptides

    Guo, J.-S.; Tian, S.-Q.; Zhang, Y.-T. Reactive molecular dynamics simulations on interaction mechanisms of cold atmospheric plasmas and peptides. Physics of Plasmas 2023, 30, 043512

  32. [40]

    Chemical fingerprints of cold physical plasmas – an experimental and computational study using cysteine as tracer compound

    Lackmann, J.-W.; Wende, K.; Verlackt, C.; Golda, J.; Volzke, J.; Kogelheide, F.; Held, J.; Bekeschus, S.; Bogaerts, A.; Schulz-von der Gathen, V.; Stapelmann, K. Chemical fingerprints of cold physical plasmas – an experimental and computational study using cysteine as tracer c...

  33. [41]

    Gonz \'a lez-Benjumea, A.; Carro, J.; Renau-M \'i nguez, C.; Linde, D.; Fern \'a ndez-Fueyo, E.; Guti \'e rrez, A.; Mart \'i nez, A. T. Fatty acid epoxidation by Collariella virescens peroxygenase and heme-channel variants. Catalysis Science & Technology 2020, 10, 717--725

  34. [42]

    Chemical modification of amino acids by atmospheric-pressure cold plasma in aqueous solution

    Takai, E.; Kitamura, T.; Kuwabara, J.; Ikawa, S.; Yoshizawa, S.; Shiraki, K.; Kawasaki, H.; Arakawa, R.; Kitano, K. Chemical modification of amino acids by atmospheric-pressure cold plasma in aqueous solution. Plasma Sources Science and Technology 2014, 47, 285403

  35. [43]

    E.-D.; Aadil, R

    Liu, Z.-W.; Zhou, Y.-X.; Wang, F.; Tan, Y.-C.; Cheng, J.-H.; Bekhit, A. E.-D.; Aadil, R. M.; Liu, X.-B. Oxidation induced by dielectric barrier discharge (DBD) plasma treatment reduces IgG/IgE binding capacity and improves the functionality of glycinin. Food chemistry 2021, 36...

  36. [44]

    Lackmann, J.-W.; Baldus, S.; Steinborn, E.; Edengeiser, E.; Kogelheide, F.; Langklotz, S.; Schneider, S.; Leichert, L. I. O.; Benedikt, J.; Awakowicz, P.; Bandow, J. E. A dielectric barrier discharge terminally inactivates RNase A by oxidizing sulfur-containing amino acids and...

  37. [45]

    Yayci, A.; Dirks, T.; Kogelheide, F.; Alcalde, M.; Hollmann, F.; Awakowicz, P.; Bandow, J. E. Protection strategies for biocatalytic proteins under plasma treatment. Journal of Physics D: Applied Physics 2020, 54, 035204

  38. [46]

    J.; Van Genechten, K.; Gasteiger, J

    Mortier, W. J.; Van Genechten, K.; Gasteiger, J. Electronegativity equalization: application and parametrization. Journal of the American Chemical Society 1985, 107, 829--835

  39. [47]

    Explanation to the species nature in ReaxFF

    Senftle, T. P.; Hong, S.; Islam, M. M.; Kylasa, S. B.; Zheng, Y.; Shin, Y. K.; Junkermeier, C.; Engel-Herbert, R.; Janik, M. J.; Aktulga, H. M.; Verstraelen, T.; Grama, A.; Duin, A. C. V. The ReaxFF reactive force-field: development, applications and future directions. npj Com...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.