As established nearly a century ago, mechanoradicals originate from homolytic bond scission in polymers. The existence, nature and biological relevance of mechanoradicals in proteins, instead, are unknown. We here show that mechanical stress on collagen produces radicals and subsequently reactive oxygen species, essential biological signaling molecules. Electron-paramagnetic resonance (EPR) spectroscopy of stretched rat tail tendon, atomistic molecular dynamics simulations and quantum-chemical calculations show that the radicals form by bond scission in the direct vicinity of crosslinks in collagen. Radicals migrate to adjacent clusters of aromatic residues and stabilize on oxidized tyrosyl radicals, giving rise to a distinct EPR spectrum consistent with a stable dihydroxyphenylalanine (DOPA) radical. The protein mechanoradicals, as a yet undiscovered source of oxidative stress, finally convert into hydrogen peroxide. Our study suggests collagen I to have evolved as a radical sponge against mechano-oxidative damage and proposes a mechanism for exercise-induced oxidative stress and redox-mediated pathophysiological processes.

No abstract

Proteins are exposed to various mechanical loads that can lead to covalent bond scissions even before macroscopic failure occurs. Knowledge of these molecular breakages is important to understand mechanical properties of the protein. In regular molecular dynamics (MD) simulations, covalent bonds are predefined, and reactions cannot occur. Furthermore, such events rarely take place on MD time scales. Existing approaches that tackle this limitation either rely on computationally expensive quantum calculations (e.g., QM/MM) or complex bond order formalisms in force fields (e.g., ReaxFF). To circumvent these limitations, we present a new reactive kinetic Monte Carlo/molecular dynamics (KIMMDY) scheme. Here, bond rupture rates are calculated based on the interatomic distances in the MD simulation and then serve as an input for a kinetic Monte Carlo step. This easily scalable hybrid approach drastically increases the accessible time scales. Using this new technique, we investigate bond ruptures in a multimillion atom system of tensed collagen, a structural protein found in skin, bones, and tendons. Our findings show a clear concentration of bond scissions near chemical cross-links in collagen. We also examine subsequent dynamic relaxation steps. Our method exhibits only a minor slowdown compared to classical MD and is straightforwardly applicable to other complex (bio)­materials under load and related chemistries.