Computational Modeling Analysis of Generation of Reactive Oxygen Species by Mitochondrial Assembled and Disintegrated...

Reactive oxygen species (ROS) function as critical mediators in a broad range of cellular signaling processes. The mitochondrial electron transport chain is one of the major contributors to ROS formation in most cells. Increasing evidence indicates that the respiratory Complex II (CII) can be the predominant ROS generator under certain conditions. A computational, mechanistic model of electron transfer and ROS formation in CII was developed in the present study to facilitate quantitative analysis of mitochondrial ROS production. The model was calibrated by fitting the computer simulated results to experimental data obtained on submitochondrial particles (SMP) prepared from bovine and rat heart mitochondria upon inhibition of the ubiquinone (Q)-binding site by atpenin A5 (AA5) and Complex III by myxothiazol, respectively. The model predicts that only reduced flavin adenine dinucleotide (FADH₂) in the unoccupied dicarboxylate state and flavin semiquinone radical (FADH^•) feature the experimentally observed bell-shaped dependence of the rate of ROS production on the succinate concentration upon inhibition of respiratory Complex III (CIII) or Q-binding site of CII, i.e., suppression of succinate-Q reductase (SQR) activity. The other redox centers of CII such as Fe-S clusters and Q-binding site have a hyperbolic dependence of ROS formation on the succinate concentration with very small maximal rate under any condition and cannot be considered as substantial ROS generators in CII. Computer simulation results show that CII disintegration (which results in dissociation of the hydrophilic SDHA/SDHB subunits from the inner membrane to the mitochondrial matrix) causes crucial changes in the kinetics of ROS production by CII that are qualitatively and quantitatively close to changes in the kinetics of ROS production by assembled CII upon inhibition of CIII or Q-binding site of CII. Thus, the main conclusions from the present computational modeling study are the following: (i) the impairment of the SQR activity of CII resulting from inhibition of CIII or Q-binding site of CII and (ii) CII disintegration causes a transition in the succinate-dependence of ROS production from a small-amplitude sigmoid (hyperbolic) shape, determined by Q-binding site or [3Fe-4S] cluster to a high-amplitude bell-shaped kinetics with a shift to small subsaturated concentrations of succinate, determined by the flavin site.