Collagen breaks at weak sacrificial bonds taming its mechanoradicals.

Collagen is a force-bearing, hierarchical structural protein important to all connective tissue. In tendon collagen, high load even below macroscopic failure level creates mechanoradicals by homolytic bond scission, similar to polymers. The location and type of initial rupture sites critically decide on both the mechanical and chemical impact of these micro-ruptures on the tissue, but are yet to be explored. We here use scale-bridging simulations supported by gel electrophoresis and mass spectrometry to determine breakage points in collagen. We find collagen crosslinks, as opposed to the backbone, to harbor the weakest bonds, with one particular bond in trivalent crosslinks as the most dominant rupture site. We identify this bond as sacrificial, rupturing prior to other bonds while maintaining the material's integrity. Also, collagen's weak bonds funnel ruptures such that the potentially harmful mechanoradicals are readily stabilized. Our results suggest this unique failure mode of collagen to be tailored towards combatting an early onset of macroscopic failure and material ageing.

ColBuilder: A server to build collagen fibril models.

Type I collagen is the main structural component of many tissues in the human body. It provides excellent mechanical properties to connective tissue and acts as a protein interaction hub. There is thus a wide interest in understanding the properties and diverse functions of type I collagen at the molecular level. A precondition is an atomistic collagen I structure as it occurs in native tissue. To this end, we built full-atom models of cross-linked collagen fibrils by integrating the low-resolution structure of collagen fibril available from x-ray fiber diffraction with high-resolution structures of short collagen-like peptides from x-ray crystallography and mass spectrometry data. We created a Web resource of collagen models for 20 different species with a large variety of cross-link types and localization within the fibril to facilitate structure-based analyses and simulations of type I collagen in health and disease. To easily enable simulations, we provide parameters of the modeled cross-links for an Amber force field. The repository of collagen models is available at https://colbuilder.h-its.org.

Mechanoradicals in tensed tendon collagen as a source of oxidative stress.

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.

Hybrid Kinetic Monte Carlo/Molecular Dynamics Simulations of Bond Scissions in Proteins.

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.