Next Article in Journal
An Intensional Probability Theory: Investigating the Link between Classical and Quantum Probabilities
Next Article in Special Issue
Mathematical and Experimental Model of Neuronal Oscillator Based on Memristor-Based Nonlinearity
Previous Article in Journal
Z-Symmetric Manifolds Admitting Schouten Tensor
Previous Article in Special Issue
Modeling the Cognitive Activity of an Individual Based on the Mathematical Apparatus of Self-Oscillatory Quantum Mechanics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Kinetics of the Singlet Oxygen Effect in Aqueous Solutions of Proteins Exposed to Thermal and Laser Radiation

by
Alexey V. Shkirin
1,2,
Sergey N. Chirikov
2,
Nikolai V. Suyazov
1,
Veronika E. Reut
3,
Daria V. Grigorieva
3,
Irina V. Gorudko
3,
Vadim I. Bruskov
4 and
Sergey V. Gudkov
1,*
1
Prokhorov General Physics Institute of the Russian Academy of Sciences, Vavilova st. 38, Moscow 119991, Russia
2
Laser Physics Department, National Research Nuclear University MEPhI, Kashirskoe sh. 31, Moscow 115409, Russia
3
Physics Department, Belarusian State University, 4 Nezavisimosti Av., 220030 Minsk, Belarus
4
Institute of Theoretical and Experimental Biophysics of the Russian Academy of Sciences, Institutskaya st. 3, Pushchino 142290, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(22), 4295; https://doi.org/10.3390/math10224295
Submission received: 18 October 2022 / Revised: 12 November 2022 / Accepted: 14 November 2022 / Published: 16 November 2022

Abstract

:
A system of kinetic equations describing the changes in the concentration of reactive oxygen species (ROS) in aqueous solutions of proteins was obtained from the analysis of chemical reactions involving singlet oxygen. Applying the condition of the stationarity of the intermediate products to the system, we determined the functional dependence of the hydrogen peroxide concentration on the protein concentration under the action of thermal and laser radiation. An approximate analytical solution to the nonlinear system of differential equations that define the ROS concentration dynamics was found. For aqueous solutions of bovine serum albumin (BSA) and bovine gamma globulin (BGG), the orders and rate constants of the reactions describing the ROS conversions were determined by minimizing the sum of squared deviations of the functions found by solving both the static and dynamic problems from experimentally measured dependences. When solving the optimization problem, the Levenberg–Marquardt algorithm was used.

1. Introduction

Modeling the kinetics of biochemical processes can give an understanding of the main formation mechanisms and the roles of the various intermediate and final reaction products (see, for example, [1,2,3]). More specifically, it helps to assess the influence of external factors (temperature, type of solvent, concentration of the initial components, exposure to ionizing and non-ionizing radiation, etc.), to determine the orders and reaction rates, and to predict the results when changing the type of some of the starting products. In the present work, this kind of mathematical modeling was carried out for reactive oxygen species (ROS) in aqueous solutions of proteins. The term reactive oxygen species (ROS) describes a range of oxygen-derived reactive molecules and free radicals. ROS are constantly formed in aerobic cells as a result of normal metabolism [4] and when exposed to external environmental factors such as ionizing, ultraviolet, extremely high frequency (EHF) microwave radiation [5], plasma [6], xenobiotics [7], thermal exposure [8], and visible and infrared optical radiation [9]. Light exposure (usually by lasers) on biochemical systems leading to the formation of ROS is known as the light-oxygen effect [10,11,12].
The relevance of this topic is associated, in particular, with the fact that an increase in the intracellular concentration of ROS above a certain level causes “oxidative stress”, which is accompanied by damage to biological vital structures, leading to lipid peroxidation, oxidative modification of nucleic acids and proteins, and the disruption of the processes necessary for the life of cells [13,14]. However, it was found that ROS exert not only a damaging effect but also perform a signal-regulatory role in the body associated with a violation of redox homeostasis [15,16]. This leads to the activation of adaptive and protective processes in cells [17,18].
For example, in mammals, a small change in the content of hydrogen peroxide caused by various physical influences may be an important factor in the therapeutic effect and adaptation of the body to adverse environmental conditions [19,20,21]. A number of physiotherapeutic procedures employing sources of coherent and incoherent light exposure are widely used in medical practice; however, the biological mechanisms of their therapeutic effect are still poorly studied [22]. Under the influence of the electromagnetic radiation of the visible and near-infrared ranges on aqueous solutions saturated with atmospheric oxygen, the formation of ROS (singlet oxygen, superoxide anion radical, hydroxyl radicals, and hydrogen peroxide) is observed [23]. Proteins, due to their high reactivity to oxidation, are the most vulnerable target for ROS compared to lipids, DNA, and other components of biological systems. Oxidation with the help of ROS leads to the loss of the specific functions of proteins, which may be the reason for their destruction and pathological processes leading to aging and the acceleration of its effects [24].
The aim of this work was to construct a kinetic model of the formation of ROS in aqueous solutions of proteins upon irradiation with a heat source and a helium–neon laser at a wavelength of 632.8 nm and then to determine the orders and rate constants for the reactions of formation and decomposition of the ROS by minimizing the sum of squared deviations of the functions found within the kinetic model from the experimentally measured dependences of the hydrogen peroxide concentration on time and the concentration of the dissolved protein in aqueous solutions of BGG and BSA.
To solve the optimization problems when applying the kinetic model to the analysis of the experimental data, the Levenberg–Marquardt method [25,26,27] was used. The method is based on an iterative procedure for finding the minimum of the objective function in a direction that is a combination of the Newton–Gauss direction and the steepest descent. When the parameters that determine the sum of squared deviations are far from their optimal value, steps are taken practically in the direction of gradient descent, when the parameters are close to their optimal value, in the direction determined by the Newton–Gauss method. Thus, an increase in the step along the gentle slopes of the function and a decrease along the steep descents are achieved, which makes it possible to efficiently find minima even with a ravine relief of the function. The Levenberg–Marquardt method is characterized by increased stability and is less demanding on the initial approximation. The method has quadratic convergence, like Newton’s method, if the Jacobi matrix is nondegenerate and Lipschitz continuous. It has become a standard method not only for solving nonlinear least squares problems, but is also widely used in solving other problems, such as in artificial neural networks [28], robotics [29], speech recognition systems [30], optimization of technological processes [31], and algorithms for determining the coordinates and orientation of an object [32]. The method has a disadvantage in that at each step of the iterative process it is necessary, in addition to the gradient of the function, to calculate the Hessian and solve the problem of matrix inversion. This affects the speed of the algorithm when solving multi-parameter (with more than a hundred parameters) problems. The impossibility of approximating the Hessian in terms of the Jacobian in the case of large residuals also reduces the performance of the algorithm. Some modified versions of the Levenberg–Marquardt method have been developed: a modified method with cubic convergence [33]; a high-order Levenberg–Marquardt method, which has biquadratic convergence [34]; a fourth-order method, which is called the accelerated modified Levenberg–Marquardt method [35]; and a three-step Levenberg–Marquardt (TSLM) algorithm based on fuzzy logic theory (FLT) [36]. However, the modified methods require a large amount of computation, which reduces the advantage of the Levenberg–Marquardt algorithm. When solving problems with a small number of parameters, the method gives excellent results, which is confirmed by its approbation on a number of standard test functions.

2. Materials and Methods

2.1. Mathematical Algorithm

The approach to determining the orders and rate constants of the reactions is as follows. The experimental data were approximated by a model function f x , b depending on m parameters b = b 1 , , b m , which is a solution to the system of kinetic equations. The parameter vector was found by minimizing the functional:
Φ b = 1 2 i = 1 m y i f x i , b 2 ,
where x i , y i are the coordinates of the experimentally measured points. When searching for the minimum of the functional in accordance with the Levenberg–Marquardt method [25,26,27], the step at the k-th iteration is carried out in the direction of the vector h :
h = H + λ · diag H 1 · Φ b k ,
where H is the Hessian, which, for small values of the residual, can be expressed in terms of the Jacobian J :
H b = J b T · J b ,
where the criterion for the end of the iterative procedure was the increase in the norm of the residual vector by a value not exceeding 10 8 of its value at the next step.

2.2. Physical Experiment Protocol

The experimental data on the concentration of hydrogen peroxide in the protein solutions were obtained according to the following protocol. To prepare samples of the protein aqueous solutions, BSA and BGG (Sigma Aldrich, St. Louis, MO, USA) and distilled water (Merck Millipore, San Francisco, CA, USA) were used. Thermal and laser radiation exposure of the samples was carried out with the aid of a U-10 liquid ultrathermostat (Prüfgeräte-WerkMedigen, Berlin, Germany) and a helium–neon laser with a wavelength of 632.8 nm and a power of 3 mW, respectively. The concentration of hydrogen peroxide was measured by enhanced chemiluminescence using a setup (Figure 1) described in detail earlier [24].
For the comparison of the thermal and laser effects, two experimental protocols were used: the first was to study the dependence of the hydrogen peroxide concentration on the protein concentration, and the second was to study the time dependence of the hydrogen peroxide concentration. In the first case, aqueous solutions of BGG and BSA at each concentration in the range of 0.1–25 µM were exposed separately to moderate thermal irradiation at 40 °C for 2 h or to laser radiation for 30 min, and then, the hydrogen peroxide concentration was measured. In the second case, aqueous solutions of BGG and BSA with concentrations of 2 µM and 10 µM, respectively, were exposed separately to moderate thermal irradiation at 45 °C for 2 h or to laser radiation for 15 min, after which the time dependence of the hydrogen peroxide concentration was recorded for 6 h.

3. Results

3.1. ROS Reactions in Aqueous Protein Solutions

Under the influence of both thermal radiation and a helium–neon laser on aqueous solutions of proteins saturated with atmospheric oxygen, the formation of singlet oxygen occurs, which is accompanied by the generation of ROS [24]. In the absence of protein in an aqueous solution saturated with atmospheric oxygen, singlet oxygen forms under the influence of electromagnetic radiation:
O2 + 1O2,
where is the energy quantum of the electromagnetic radiation corresponding to the transition of oxygen from the basic triplet to one of the singlet states: 1Δg or 1Σg+, whose position corresponds to wavelengths of 1269 nm and 763 nm [37]. The state 1Σg+ is short-lived and quickly relaxes in 1Δg. It is important that a transition in collision complexes (dimoles) is possible, 3Σg + 3Σg + 1Δg + 1Δg, with light absorption at a wavelength of ~630 nm (bandwidth of 18 nm), which allows the use of helium–neon laser radiation with a wavelength of 632.8 nm to generate singlet oxygen [10].
When singlet oxygen interacts with hydroxyl and hydrogen ions, hydroxyl, hydroperoxide radical, and hydrogen peroxide are formed [24]:
OH + 1O2 + H+OH + HO2,
HO2 + HO2 → H2O2 + 1O2,
OH + OH → H2O2,
The reaction of OH with a hydrogen atom at the α-carbon atom of the polypeptide chain in the protein molecule P leads to the detachment of the hydrogen atom and the formation of the α-carbon radical [38]:
P + OH → P + H,
As a result of the reaction of the α-carbon radical with singlet oxygen, a peroxyl radical of the protein is formed [11]:
P + 1O2 → POO,
Next is the elimination of the hydroperoxide from the peroxyradical protein:
2 POO+H → 2 Pe + HO2,
and the additional formation of hydrogen peroxide due to the dismutation of hydroperoxide radicals according to reaction (6). Along with this, when interacting with singlet oxygen, a number of lateral amino acid residues are also modified [39].
If we summarize reactions (8) and (9), we obtain the formation of a protein peroxy radical as a result of the combined action of hydroxyl radicals and singlet oxygen on the protein:
P + OH + 1O2 → POO+ H,
Peroxyradical protein is a long-lived reactive protein species (LRPS). The further conversion of LRPS is accompanied by the release of hydroperoxide, the dismutation of two molecules of which leads to the formation of hydrogen peroxide:
2POO + P → 2Pe + H2O2 + 1O2,
where Pe is the final product of the conversion of LRPS from the original protein.
However, for the non-monotonic generation of hydrogen peroxide with the participation of LRPS, one more reaction of its decomposition by protein is additionally required:
2H2O2 + P → 2 H2O + 1O2,
The kinetic model formulated below predicts the existence of reaction (12), which is a non-catalytic form of the catalase activity of all proteins [40].

3.2. Kinetic Model

Taking into account the reactions involving singlet oxygen in protein aqueous solutions (5)–(13), we can describe the hydrogen peroxide generation observed in the experiments by a simplified minimal system of chemical equations:
P   +   1 O 2   k 1   P e + H O 2 ,
H O 2   +   H O 2   k 2   H 2 O 2   +   1 O 2 ,
2 H 2 O 2   +   P     k 3   2 H 2 O   +   1 O 2   +   P ,
2 H 2 O 2   k 4   2 H 2 O   +   O 2 ,
Reaction (14) describes the formation of the HO2 radical, which, in turn, leads to the formation of H2O2 in accordance with (15). In the proposed model, two channels of peroxide decomposition are considered: decomposition due to the catalase activity of the protein (16) and self-decomposition (17). For the decay reactions, a mismatch of the experimentally determined reaction orders with the formal stoichiometric coefficients is typical; these orders can be both integer and fractional [41]. Therefore, the kinetic equations for the concentrations of ROS (HO2 and H2O2), corresponding to the above system of chemical reactions, can be written in the general form:
d [ H O 2 ] 2 / d t = k 1 [ P ] [ 1 O 2 ] 2   k 2 [ H O 2 ] 2 ,
d H 2 O 2 / d t = k 2 · [ HO 2 ] 2 2 · k 3 · [ P ] γ H 2 O 2 β 2 · k 4 · H 2 O 2 α ,
where k i (i = 1 … 4) are the rate constants of the corresponding reactions, and α, β, γ are the orders of the hydrogen peroxide decomposition reactions. The values of α, β, γ are assumed to be determined from the experimental data.

3.3. Stationary Dependences of the Hydrogen Peroxide Concentration on the Protein Concentration in Aqueous Solutions of BSA and BGG

Applying the stationarity condition for the intermediate reaction products to the system (14)–(17), i.e., the equality of profit and loss for H2O2 and HO2, we obtain
for HO2:
k 1 · [ P ] · [ O 2 1 ] = 2 · k 2 [ H O 2 ] 2 ,
for H2O2:
k 2 [ H O 2 ] 2 = 2 k 3 [ P ] γ H 2 O 2 β + 2 k 4 H 2 O 2 α .
After substituting (20) into (21), we have:
k 1 [ P ] [ O 2 1 ] = 4 k 3 [ P ] γ H 2 O 2 β + 4 k 4 H 2 O 2 α
It is known [41,42] that the self-decomposition of hydrogen peroxide (17) has the order of 1 ≤ α ≤ 2. It is reasonable to assume that the decomposition of peroxide by protein (16) is described, similarly to self-decomposition, by the order of peroxide 1 ≤ β ≤ 2. The reaction order for the protein in (16) is not known to us in advance; however, it cannot significantly exceed the order of peroxide; so, it makes sense to search for it in the range 1 ≤ γ ≤ 3.
The formation of H2O2 was experimentally studied under the influence of moderate thermal radiation at 40 °C for 2 h on the BGG and BSA solutions and upon exposure to He–Ne laser radiation for 30 min. The measured dependences of the H2O2 concentration on the concentration of BGG and BSA are shown in Figure 2 and Figure 3, respectively.
In order to simplify the derivation of an explicit function for the H2O2 concentration from the Equation (22), we set α and β equal to 1, which is typical for the decomposition of peroxide under normal conditions. Then, the dependence of the peroxide concentration on the protein concentration takes the form:
  H 2 O 2 = k 1 P O 2   1 4 k 3 P γ + k 4 .
This dependence can be represented as the following function:
y x = x a x γ + b ,
a = 4 k 3 k 1 O 2   1 ;     b = 4 k 4 k 1 O 2   1 .
We searched for the coefficients ( a , b , γ) by minimizing the deviations of the function (24) from the experimental data (Figure 2 and Figure 3) using the Levenberg–Marquardt method [25,26,27]; the found values ( a , b , γ) for the BSA BGG aqueous solutions are given in Table 1.

3.4. Dynamics of Hydrogen Peroxide Formation in Aqueous Solutions of BSA and BGG

Assuming the reaction orders (α, β) are the same as in Section 3.2, we can write the system of kinetic Equations (18) and (19) in the following functional form ( z t ≡ [HO2], y t ≡ [H2O2]) with free coefficients (A, B, C):
d d t z t = A 2 B z t 2 ,     z 0 = z 0 ,
d d t y t = B z t 2 C y t ,     y 0 = y 0 .
The time dependence y t will have a maximum when the condition z 0 > A 2 B is met. In this case, the solution to system (26), (27) has the form:
z t = A 2 B c o t h 2 A B t + arccoth z 0 2 B A ,    
y t = A 2 e C t 0 t c o t h 2 A B t   +   a r c c o t h z 0 2 B A 2 e C τ d τ + y 0 e C t .
The search for free coefficients is greatly simplified if (29) is represented through known functions. Replacing c o t h x 2 in the integral (29) for x > 0 with an approximate expression 1 / x 2 + 1 e 3 x / 2 / 3 x e 53 x / 60 / 2 (the relative error is no more than 1.3 · 10 2 ), we obtain the following function:
+ A y ( t ) y 0 exp ( C t ) + ξ η 0 2 + C 2 ξ η 0 C exp ( C t ) + 1 2 Ei η 0 C ξ + 1 2 Ei C t + η 0 C ξ C exp C t η 0 C ξ ξ 2 + 900 η 0 C 795 ξ η 0 900 ξ exp C t 53 η 0 60 ( 60 C 53 ξ ) 2 + 1 3 exp C t 3 η 0 2 2 C 3 ξ + ( 15 ξ t + η 0 ( 60 C 53 ξ ) + 900 ξ exp 53 60 ξ t + η 0 ( 60 C 53 ξ ) 2 1 3 exp 3 2 ξ t + η 0 2 C 3 ξ ξ 2 t ξ η 0 + C 2 ξ ξ t + η 0 C
Here ξ = 2 A B , η 0 = arccoth z 0 2 B A , and Ei x = x e s s d s is the exponential integral. The experimental time dependences of the hydrogen peroxide concentration in the BGG and BSA solutions (concentration 2 μM and 10 μM, respectively) exposed to thermal and laser radiation, as well as the results of their approximation by the Levenberg–Marquardt method using the solution (30), are shown in Figure 4 and Figure 5.
The found parameters of the function (30) in the units indicated on the graphs (Figure 4 and Figure 5) are shown in Table 2.

3.5. Chemical Constants

The rate constants for the reactions (14)–(17) are related to the parameters of the functions (24) and (30), which describe the process of hydrogen peroxide generation, by the relations:
2 k 4 + k 3 P γ = C ,     k 3 = a b k 4 ,     k 2 = B = ξ 2 2 A ,     k 1 e f f = k 1 O 2   1 = 4 k 4 b .
Having the initially known protein concentrations in the kinetic experiments with BGG ([P] = 2 10 micromole/L) and BSA ([P] = 10 micromole/L), and taking ( a , b , γ) from Table 1, we made estimates of the reaction rate constants k 1 , k 2 , k 3 , and k 4 (Table 3), according to (31).
Note that the found values of k 4 are in an order of magnitude consistent with the interpolated data on the temperature decomposition of hydrogen peroxide [42].

4. Discussion

Considering the effects of singlet oxygen leading to the transformation of LPRS, we determined the minimal system of chemical reactions for ROS conversions in aqueous solutions of proteins and derived the corresponding kinetic equations. By comparison of the obtained kinetic model and the measured dependences of H2O2 concentration in BSA and BGG solutions subjected to thermal and laser radiation, we determined the reaction order of the protein-induced decomposition of hydrogen peroxide with respect to the protein and the kinetic constants describing the ROS conversions.
The calculation results (Table 3) show that, compared with laser exposure, under thermal exposure the rate k 1 e f f of the first reaction (14) is several times greater, and the rates k 2 , k 3 , and k 4 of the subsequent stages (15)–(17) significantly decrease. In contrast with our earlier work, the performed kinetics modeling has an important new feature such that the reaction order with respect to the protein is also included in the set of parameters to be determined using the Levenberg–Marquardt optimization algorithm and, as a result, its value is determined more accurately. Such refinement of the model gives more reliable values of the rates of the intermediate reactions k j than those obtained with a fixed preset value of the reaction order. The time dependence of the H2O2 and HO2 concentrations derived in an analytical form containing the kinetic parameters is applicable for describing the ROS dynamics in solutions of various proteins under both types of external exposure.
The similarity of the functional dependences of the H2O2 concentration on the protein content and time in solutions of the different proteins, both under the action of thermal radiation and the radiation of the He–Ne laser, allowed us to conclude that the fundamental mechanisms of the formation of H2O2 are identical in these cases. Thus, similar mathematical approaches can be applied to modeling the ROS conversions in various protein systems, including those that are substantially close to living systems.

5. Conclusions

Hydrogen peroxide is a harmful byproduct of many normal metabolic processes; to prevent damage to cells and tissues, it must be quickly converted into other, less dangerous substances. In this work, we have shown that under thermal/laser irradiation in the visible range, the concentration of protein dissolved in water can serve as a parameter for controlling the decomposition of hydrogen peroxide into the less reactive gaseous molecules of oxygen and water. We have proposed a kinetic model which includes the reaction equation for the catalysis of hydrogen peroxide decomposition due to protein macromolecules. This made it possible to explain the passage of the maximum with the subsequent tending to zero for the experimental dependence of the peroxide concentration on the protein concentration. The dynamics of this process are well described by a system of nonlinear differential equations; good agreement was obtained between the theoretical solution to this system and the experimental measurements.

Author Contributions

Conceptualization, S.V.G., V.I.B. and A.V.S.; methodology, A.V.S., S.N.C. and N.V.S.; software, S.N.C.; validation, V.I.B.; formal analysis, N.V.S.; investigation, V.E.R., D.V.G., I.V.G., A.V.S., S.N.C. and N.V.S.; resources, S.V.G.; data curation, S.N.C.; writing—original draft preparation, A.V.S. and S.N.C.; writing—review and editing, S.V.G.; visualization, S.N.C.; supervision, A.V.S.; project administration, S.V.G.; funding acquisition, S.V.G. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant of the Ministry of Science and Higher Education of the Russian Federation (075-15-2022-315) for the organization and development of a world-class research center ‘‘Photonics”.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bitsouni, V.; Gialelis, N.; Stratis, I.G. Rigorous Analysis of the Quasi-Steady-State Assumption in Enzyme Kinetics. Mathematics 2022, 10, 1086. [Google Scholar] [CrossRef]
  2. Nikolov, S.G.; Wolkenhauer, O.; Nenov, M.; Vera, J. Detection and Analysis of Critical Dynamic Properties of Oligodendrocyte Differentiation. Mathematics 2022, 10, 2928. [Google Scholar] [CrossRef]
  3. Sukhinov, A.; Belova, Y.; Nikitina, A.; Sidoryakina, V. Sufficient Conditions for the Existence and Uniqueness of the Solution of the Dynamics of Biogeochemical Cycles in Coastal Systems Problem. Mathematics 2022, 10, 2092. [Google Scholar] [CrossRef]
  4. Snezhkina, A.V.; Kudryavtseva, A.V.; Kardymon, O.L.; Savvateeva, M.V.; Melnikova, N.V.; Krasnov, G.S.; Dmitriev, A.A. ROS generation and antioxidant defense systems in normal and malignant cells. Oxid. Med. Cell. Longev. 2019, 2019, 6175804. [Google Scholar] [CrossRef] [Green Version]
  5. Gapeyev, A.B.; Lukyanova, N.A.; Gudkov, S.V. Hydrogen peroxide induced by modulated electromagnetic radiation protects the cells from DNA damage. Cent. Eur. J. Biol. 2014, 9, 915–921. [Google Scholar] [CrossRef]
  6. Piskarev, I. The role of light emission of spark discharge in air in the generation of active species. High Energy Chem. 2019, 53, 92–93. [Google Scholar] [CrossRef]
  7. Filev, A.D.; Shmarina, G.V.; Ershova, E.S.; Veiko, N.N.; Martynov, A.V.; Borzikova, M.A.; Poletkina, A.A.; Dolgikh, O.A.; Veiko, V.P.; Bekker, A.A.; et al. Oxidized cell-free DNA role in the antioxidant defense mechanisms under stress. Oxid. Med. Cell. Longev. 2019, 2019, 1245749. [Google Scholar] [CrossRef]
  8. Bruskov, V.I.; Malakhova, L.V.; Masalimov, Z.K.; Chernikov, A.V. Heat-induced formation of reactive oxygen species and 8-oxoguanine, a biomarker of damage to DNA. Nucleic Acids Res. 2002, 30, 1354–1363. [Google Scholar] [CrossRef] [Green Version]
  9. Gudkov, S.V.; Bruskov, V.I.; Astashev, M.E.; Chernikov, A.V.; Yaguzhinsky, L.S.; Zakharov, S.D. Oxygen-dependent auto-oscillations of water luminescence triggered by the 1264 nm radiation. J. Phys. Chem. B 2011, 115, 7693–7698. [Google Scholar] [CrossRef]
  10. Zakharov, S.D.; Ivanov, A.V. Light-oxygen effect in cells and its potential applications in tumour therapy. Quantum Electron. 1999, 29, 1031. [Google Scholar] [CrossRef]
  11. Zakharov, S.D.; Ivanov, A.V.; Wolf, E.; Danilov, V.P.; Murina, T.M.; Nguen, K.; Novikov, E.; Panasenko, N.; Perov, S.; Skopinov, S. Structural rearrangements in the aqueous phase of cell suspensions and protein solutions induced by a light-oxygen effect. Quantum Electron. 2003, 33, 149. [Google Scholar] [CrossRef]
  12. Chin, K.K.; Trevithick-Sutton, C.C.; McCallum, J.; Jockusch, S.; Turro, N.J.; Scaiano, J.; Foote, C.S.; Garcia-Garibay, M.A. Quantitative determination of singlet oxygen generated by excited state aromatic amino acids, proteins, and immunoglobulins. J. Am. Chem. Soc. 2008, 130, 6912–6913. [Google Scholar] [CrossRef] [PubMed]
  13. Forman, H.J.; Zhang, H. Targeting oxidative stress in disease: Promise and limitations of antioxidant therapy. Nat. Rev. Drug Discov. 2021, 20, 689–709. [Google Scholar] [CrossRef] [PubMed]
  14. Pisoschi, A.M.; Pop, A.; Iordache, F.; Stanca, L.; Predoi, G.; Serban, A.I. Oxidative stress mitigation by antioxidants-an overview on their chemistry and influences on health status. Eur. J. Med. Chem. 2021, 209, 112891. [Google Scholar] [CrossRef] [PubMed]
  15. Brillo, V.; Chieregato, L.; Leanza, L.; Muccioli, S.; Costa, R. Mitochondrial dynamics, ROS, and cell signaling: A blended overview. Life 2021, 11, 332. [Google Scholar] [CrossRef]
  16. Ochoa, C.D.; Wu, R.F.; Terada, L.S. ROS signaling and ER stress in cardiovascular disease. Mol. Aspects Med. 2018, 63, 18–29. [Google Scholar] [CrossRef]
  17. Sharapov, M.; Novoselov, V. Catalytic and signaling role of peroxiredoxins in carcinogenesis. Biochem. 2019, 84, 79–100. [Google Scholar] [CrossRef]
  18. Nakamura, H.; Takada, K. Reactive oxygen species in cancer: Current findings and future directions. Cancer Sci. 2021, 112, 3945. [Google Scholar] [CrossRef]
  19. Vlasova, I.I. Peroxidase activity of human hemoproteins: Keeping the fire under control. Molecules 2018, 23, 2561. [Google Scholar] [CrossRef] [Green Version]
  20. Borisova-Mubarakshina, M.M.; Vetoshkina, D.V.; Ivanov, B.N. Antioxidant and signaling functions of the plastoquinone pool in higher plants. Physiol. Plant. 2019, 166, 181–198. [Google Scholar] [CrossRef]
  21. Smirnova, O.A.; Bartosch, B.; Zakirova, N.F.; Kochetkov, S.N.; Ivanov, A.V. Polyamine metabolism and oxidative protein folding in the ER as ROS-producing systems neglected in virology. Int. J. Mol. Sci. 2018, 19, 1219. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Chernov, A.; Reshetnikov, D.; Yu, A.K.; Manokhin, A.; Gudkov, S. Influence of wideband visible light with an padding red component on the functional state of mice embryos and embryonic stem cells. J. Photochem. Photobiol. 2018, 188, 77–86. [Google Scholar] [CrossRef] [PubMed]
  23. Amaroli, A.; Pasquale, C.; Zekiy, A.; Utyuzh, A.; Benedicenti, S.; Signore, A.; Ravera, S. Photobiomodulation and oxidative stress: 980 nm diode laser light regulates mitochondrial activity and reactive oxygen species production. Oxid. Med. Cell. Longev. 2021, 2021, 6626286. [Google Scholar] [CrossRef] [PubMed]
  24. Ivanov, V.; Karp, O.; Bruskov, V.; Andreev, S.; Bunkin, N.; Gudkov, S. Formation of long-lived reactive products in blood serum under heat treatment and low-intensity laser irradiation, their role in hydrogen peroxide generation and DNA damage. Indian J. Biochem. Biophys. 2019, 56, 214–223. [Google Scholar] [CrossRef]
  25. Levenberg, K. A method for the solution of certain non-linear problems in least squares. Q. Appl. Math. 1944, 2, 164–168. [Google Scholar] [CrossRef] [Green Version]
  26. Marquardt, D.W. An algorithm for least-squares estimation of nonlinear parameters. SIAM J. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  27. Nocedal, J.; Wright, S. Numerical Optimization; Springer Science and Business Media: New York, NY, USA, 2006. [Google Scholar]
  28. Hua, L.; Bo, L.; Tong, L.; Wang, M.; Fu, H.; Guo, R. Angular Acceleration Sensor Fault Diagnosis Based on LM-BP Neural Network. In Proceedings of the 37th Chinese Control Conference, Wuhan, China, 25–27 July 2018; pp. 6028–6032. [Google Scholar]
  29. Luo, G.; Zou, L.; Wang, Z.; Lv, C.; Ou, J.; Huang, Y. A novel kinematic parameters calibration method for industrial robot based on Levenberg-Marquardt and Differential Evolution hybrid algorithm. Robot. Comput. Integr. Manuf. 2021, 71, 102165. [Google Scholar] [CrossRef]
  30. Kumar, S.S.; Sowmya, R.; Shankar, B.M.; Lingaraj, N.; Sivakumar, S. Analysis of Connected Word Recognition systems using Levenberg Mar-quardt Algorithm for cockpit control in unmanned aircrafts. Mater. Today Proc. 2021, 37, 1813–1819. [Google Scholar] [CrossRef]
  31. Mahmoudabadi, Z.S.; Rashidi, A.; Yousefi, M. Synthesis of 2D-Porous MoS2 as a Nanocatalyst for Oxidative Desulfuriza-tion of Sour Gas Condensate: Process Parameters Optimization Based on the Levenberg–Marquardt Algorithm. J. Environ. Chem. Eng. 2021, 9, 105200. [Google Scholar] [CrossRef]
  32. Kudinov, I.A.; Pavlov, O.V.; Kholopov, I.S. Implementation of an algorithm for determining the spatial coordinates and the angular orientation of an object based on reference marks, using information from a single camera. Comput. Opt. 2015, 39, 413–419. [Google Scholar] [CrossRef]
  33. Fan, J. The modified Levenberg–Marquardt method for nonlinear equations with cubic convergence. Math. Comp. 2012, 81, 447–466. [Google Scholar] [CrossRef] [Green Version]
  34. Yang, X. A higher-order Levenberg–Marquardt method for nonlinear equations. Appl. Math. Comput. 2013, 219, 10682–10694. [Google Scholar] [CrossRef]
  35. Chen, L. A high-order modified Levenberg–Marquardt method for systems of nonlinear equations with fourth-order convergence. Appl. Math. Comput. 2016, 285, 79–93. [Google Scholar] [CrossRef]
  36. Derakhshandeh, S.Y.; Pourbagher, R.; Kargar, A. A novel fuzzy logic Levenberg-Marquardt method to solve the ill-conditioned power flow problem. Int. J. Electr. Power Energy Syst. 2018, 99, 299–308. [Google Scholar] [CrossRef]
  37. Shinkarenko, N.V.; Aleskovskii, V. Singlet oxygen: Methods of preparation and detection. Russ. Chem. Rev. 1981, 50, 220. [Google Scholar] [CrossRef]
  38. Davies, M.J. The oxidative environment and protein damage. Biochim. Biophys. Acta Proteins Proteom. 2005, 1703, 93–109. [Google Scholar] [CrossRef] [PubMed]
  39. Davies, M.J. Reactive species formed on proteins exposed to singlet oxygen. Photochem. Photobiol. Sci. 2004, 3, 17–25. [Google Scholar] [CrossRef]
  40. Hawkins, C.L.; Davies, M.J. Generation and propagation of radical reactions on proteins. Biochim. Biophys. Acta Bioenerg. 2001, 1504, 196–219. [Google Scholar] [CrossRef] [Green Version]
  41. Panchenkov, G.M.; Lebedev, V.P. Chemical Kinetics and Catalysis; Khimiya: Moscow, Russia, 1985. [Google Scholar]
  42. Takagi, J.; Ishigure, K. Thermal decomposition of hydrogen peroxide and its effect on reactor water monitoring of boiling water reactors. Nucl. Sci. Eng. 1985, 89, 177–186. [Google Scholar] [CrossRef]
Figure 1. Scheme of the setup for measuring chemoluminescence of protein solutions: (1) Biotox-7A chemoluminometer; (2) focusing mirror; (3) cell with protein solution; (4) photo shutter; (5) optical filter; (6) photomultiplier tube; (7) personal computer.
Figure 1. Scheme of the setup for measuring chemoluminescence of protein solutions: (1) Biotox-7A chemoluminometer; (2) focusing mirror; (3) cell with protein solution; (4) photo shutter; (5) optical filter; (6) photomultiplier tube; (7) personal computer.
Mathematics 10 04295 g001
Figure 2. Variation of H2O2 concentration with the protein concentration in aqueous solutions exposed to moderate thermal radiation at 40 °C for 2 h: (a) BGG; (b) BSA. The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Figure 2. Variation of H2O2 concentration with the protein concentration in aqueous solutions exposed to moderate thermal radiation at 40 °C for 2 h: (a) BGG; (b) BSA. The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Mathematics 10 04295 g002
Figure 3. Variation of H2O2 concentration with the protein concentration in aqueous solutions exposed to He–Ne laser radiation (3 mW) for 30 min: (a) BGG; (b) BSA. The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Figure 3. Variation of H2O2 concentration with the protein concentration in aqueous solutions exposed to He–Ne laser radiation (3 mW) for 30 min: (a) BGG; (b) BSA. The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Mathematics 10 04295 g003
Figure 4. Dynamics of hydrogen peroxide generation in solutions of proteins exposed to heat (45 °C, 2 h) in phosphate buffer (10 mM Na2HPO4 and 150 mM NaCl, pH 7.4): (a) BGG (2 μM); (b) BSA (10 μM). The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Figure 4. Dynamics of hydrogen peroxide generation in solutions of proteins exposed to heat (45 °C, 2 h) in phosphate buffer (10 mM Na2HPO4 and 150 mM NaCl, pH 7.4): (a) BGG (2 μM); (b) BSA (10 μM). The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Mathematics 10 04295 g004
Figure 5. Dynamics of hydrogen peroxide generation in solutions of proteins exposed to laser radiation He–Ne laser (3 mW) for 15 min in phosphate buffer (10 mM Na2HPO4 and 150 mM NaCl, pH 7.4): (a) BGG (2 μM); (b) BSA (10 μM). The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Figure 5. Dynamics of hydrogen peroxide generation in solutions of proteins exposed to laser radiation He–Ne laser (3 mW) for 15 min in phosphate buffer (10 mM Na2HPO4 and 150 mM NaCl, pH 7.4): (a) BGG (2 μM); (b) BSA (10 μM). The circles are experimentally measured H2O2 concentrations. The solid lines display theoretical approximation of the experimental points by Levenberg–Marquardt method.
Mathematics 10 04295 g005
Table 1. Parameters of the model function (24).
Table 1. Parameters of the model function (24).
ProteinExposurea (L2·Micromole−2)bγ
BGGThermal (Figure 2a)5.1723.332.45
Laser (Figure 3a)36.45315.142.23
BSAThermal (Figure 2b)3.3276.521.81
Laser (Figure 3b)69.84125.781.43
Table 2. Parameters of the model function (30).
Table 2. Parameters of the model function (30).
ProteinExposure ξ
(h−1)
η 0 C
(h−1)
y 0
(nmole·L−1)
A
(nmole·L−1·h−1)
BGGThermal (Figure 4a)0.0660.1360.4260.0281.38·10−3
Laser (Figure 5a)0.0960.1480.9390.000234.86·10−4
BSAThermal (Figure 4b)0.0290.1630.3540.0331.34·10−3
Laser (Figure 5b)0.0440.2950.3960.000419.77·10−4
Table 3. Reaction rate constants.
Table 3. Reaction rate constants.
ProteinExposure k 1 ( e f f ) ( s 1 ) k2
( L mole 1 s 1 )
k3
( L 2 mole 2 s 1 )
k4
( s 1 )
BGGThermal4.6⋅10−64.4⋅1055.9⋅1062.7⋅10−5
Laser1.07⋅10−62.6⋅1069.8⋅1068.5⋅10−5
BSAThermal6.8⋅10−70.87⋅1055.6⋅1051.3⋅10−5
Laser1.1⋅10−72.8⋅1051.9⋅1063.5⋅10−6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Shkirin, A.V.; Chirikov, S.N.; Suyazov, N.V.; Reut, V.E.; Grigorieva, D.V.; Gorudko, I.V.; Bruskov, V.I.; Gudkov, S.V. Modeling the Kinetics of the Singlet Oxygen Effect in Aqueous Solutions of Proteins Exposed to Thermal and Laser Radiation. Mathematics 2022, 10, 4295. https://doi.org/10.3390/math10224295

AMA Style

Shkirin AV, Chirikov SN, Suyazov NV, Reut VE, Grigorieva DV, Gorudko IV, Bruskov VI, Gudkov SV. Modeling the Kinetics of the Singlet Oxygen Effect in Aqueous Solutions of Proteins Exposed to Thermal and Laser Radiation. Mathematics. 2022; 10(22):4295. https://doi.org/10.3390/math10224295

Chicago/Turabian Style

Shkirin, Alexey V., Sergey N. Chirikov, Nikolai V. Suyazov, Veronika E. Reut, Daria V. Grigorieva, Irina V. Gorudko, Vadim I. Bruskov, and Sergey V. Gudkov. 2022. "Modeling the Kinetics of the Singlet Oxygen Effect in Aqueous Solutions of Proteins Exposed to Thermal and Laser Radiation" Mathematics 10, no. 22: 4295. https://doi.org/10.3390/math10224295

APA Style

Shkirin, A. V., Chirikov, S. N., Suyazov, N. V., Reut, V. E., Grigorieva, D. V., Gorudko, I. V., Bruskov, V. I., & Gudkov, S. V. (2022). Modeling the Kinetics of the Singlet Oxygen Effect in Aqueous Solutions of Proteins Exposed to Thermal and Laser Radiation. Mathematics, 10(22), 4295. https://doi.org/10.3390/math10224295

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop