Alfonso Gautieri

Hierarchical Structure and Nanomechanics of Collagen Microfibrils from the Atomistic Scale Up

Collagen constitutes one-third of the human proteome, providing mechanical stability, elasticity, and strength to organisms and is the prime construction material in biology. Collagen is also the dominating material in the extracellular matrix and its stiffness controls cell differentiation, growth, and pathology. However, the origin of the unique mechanical properties of collagenous tissues, and in particular its stiffness, extensibility, and nonlinear mechanical response at large deformation, remains unknown. By using X-ray diffraction data of a collagen fibril (Orgel, J. P. R. O. et al. Proc. Natl. Acad. Sci. 2006, 103, 9001) here we present an experimentally validated model of the nanomechanics of a collagen microfibril that incorporates the full biochemical details of the amino acid sequence of constituting molecules and the nanoscale molecular arrangement. We demonstrate by direct mechanical testing that hydrated (wet) collagen microfibrils feature a Youngs modulus of ≈300 MPa at small, and ≈1.2 GPa at larger deformation in excess of 10% strain, which is in excellent agreement with experimental data. We find that dehydrated (dry) collagen microfibrils show a significantly increased Youngs modulus of ≈1.8-2.25 GPa, which is in agreement with experimental measurements and owing to tighter molecular packing. Our results show that the unique mechanical properties of collagen microfibrils arise due to their hierarchical structure at the nanoscale, where key deformation mechanisms are straightening of twisted triple-helical molecules at small strains, followed by axial stretching and eventual molecular uncoiling. The establishment of a model of hierarchical deformation mechanisms explains the striking difference of the elastic modulus of collagen fibrils compared with single molecules, which is found in the range of 4.8 ± 2 GPa, or ≈10-20 times greater. We find that collagen molecules alone are not capable of providing the broad range of mechanical functionality required for physiological function of collagenous tissues. Rather, the existence of an array of deformation mechanisms, derived from the hierarchical makeup of the material, is critical to the materials ability to confer key mechanical properties, specifically large extensibility, strain hardening, and toughness, despite the limitation that collagenous materials are constructed from only few distinct amino acids. The atomistic model of collagen microfibril mechanics now enables the bottom-up elucidation of structure-property relationships in a broader class of collagen materials (e.g., tendon, bone, cornea), including studies of genetic disease where the incorporation of biochemical details is essential. The availability of a molecular-based model of collagen tissues may eventually result in novel nanomedicine approaches to develop treatments for a broad class of collagen diseases and the design of de novo biomaterials for regenerative medicine.

Molecular mechanics of mineralized collagen fibrils in bone

Bone is a natural composite of collagen protein and the mineral hydroxyapatite. The structure of bone is known to be important to its load-bearing characteristics, but relatively little is known about this structure or the mechanism that govern deformation at the molecular scale. Here we perform full-atomistic calculations of the three-dimensional molecular structure of a mineralized collagen protein matrix to try to better understand its mechanical characteristics under tensile loading at various mineral densities. We find that as the mineral density increases, the tensile modulus of the network increases monotonically and well beyond that of pure collagen fibrils. Our results suggest that the mineral crystals within this network bears up to four times the stress of the collagen fibrils, whereas the collagen is predominantly responsible for the material’s deformation response. These findings reveal the mechanism by which bone is able to achieve superior energy dissipation and fracture resistance characteristics beyond its individual constituents.

Deformation rate controls elasticity and unfolding pathway of single tropocollagen molecules

Collagen is an important structural protein in vertebrates and is responsible for the integrity of many tissues like bone, teeth, cartilage and tendon. The mechanical properties of these tissues are primarily determined by their hierarchical arrangement and the role of the collagen matrix in their structures. Here we report a series of Steered Molecular Dynamics (SMD) simulations in explicit solvent, used to elucidate the influence of the pulling rate on the Youngs modulus of individual tropocollagen molecules. We stretch a collagen peptide model sequence [(Gly-Pro-Hyp)(10)](3) with pulling rates ranging from 0.01 to 100 m/s, reaching much smaller deformation rates than reported in earlier SMD studies. Our results clearly demonstrate a strong influence of the loading velocity on the observed mechanical properties. Most notably, we find that Youngs modulus converges to a constant value of approximately 4 GPa tangent modulus at 8% tensile strain when the initially crimped molecule is straightened out, for pulling rates below 0.5 m/s. This enables us for the first time to predict the elastic properties of a single tropocollagen molecule at physiologically and experimentally relevant pulling rates, directly from atomistic-level calculations. At deformation rates larger than 0.5 m/s, Youngs modulus increases continuously and approaches values in excess of 15 GPa for deformation rates larger than 100 m/s. The analyses of the molecular deformation mechanisms show that the tropocollagen molecule unfolds in distinctly different ways, depending on the loading rate, which explains the observation of different values of Youngs modulus at different loading rates. For low pulling rates, the triple helix first uncoils completely at 10%-20% strain, then undergoes some recoiling in the opposite direction, and finally straightens for strains larger than 30%. At intermediate rates, the molecule uncoils linearly with increasing strain up to 35% strain. Finally, at higher velocities the triple helix does not uncoil during stretching.

Thickness of Hydroxyapatite Nanocrystal Controls Mechanical Properties of the Collagen−Hydroxyapatite Interface

Collagen-hydroxyapatite interfaces compose an important building block of bone structures. While it is known that the nanoscale structure of this elementary building block can affect the mechanical properties of bone, a systematic understanding of the effect of the geometry on the mechanical properties of this interface between protein and mineral is lacking. Here we study the effect of geometry, different crystal surfaces, and hydration on the mechanical properties of collagen-hydroxyapatite interfaces from an atomistic perspective, and discuss underlying deformation mechanisms. We find that the presence of hydroxyapatite significantly enhances the tensile modulus and strength compared with a tropocollagen molecule alone. The stiffening effect is strongly dependent on the thickness of the mineral crystal until a plateau is reached at 2 nm crystal thickness. We observe no significant differences due to the mineral surface (Ca surface vs OH surface) or due to the presence of water. Our result shows that the hydroxyapatite crystal with its thickness confined to the nanometer size efficiently increases the tensile modulus and strength of the collagen-hydroxyapatite composite, agreeing well with experimental observations that consistently show the existence of extremely thin mineral flakes in various types of bones. We also show that the collagen-hydroxyapatite interface can be modeled with an elastic network model which, based on the results of atomistic simulations, provides a good estimate of the surface energy and other mechanical features.

Viscoelastic properties of model segments of collagen molecules.

Collagen is the prime construction material in vertebrate biology, determining the mechanical behavior of connective tissues such as tendon, bone and skin. Despite extensive efforts in the investigation of the origin of collagen unique mechanical properties, a deep understanding of the relationship between molecular structure and mechanical properties remains elusive, hindered by the complex hierarchical structure of collagen-based tissues. In particular, although extensive studies of viscoelastic properties have been pursued at the macroscopic (fiber/tissue) level, fewer investigations have been performed at the smaller scales, including in particular collagen molecules and fibrils. These scales are, however, important for a complete understanding of the role of collagen as an important constituent in the extracellular matrix. Here, using an atomistic modeling approach, we perform in silico creep tests of a collagen-like peptide, monitoring the strain-time response for different values of applied external load. The results show that individual collagen molecules exhibit a nonlinear viscoelastic behavior, with a Youngs modulus increasing from 6 to 16GPa (for strains up to 20%), a viscosity of 3.84.±0.38Pa·s, and a relaxation time in the range of 0.24-0.64ns. The single molecule viscosity, for the first time reported here, is several orders of magnitude lower than the viscosity found for larger-scale single collagen fibrils, suggesting that the viscous behavior of collagen fibrils and fibers involves additional mechanisms, such as molecular sliding between collagen molecules within the fibril or the effect of relaxation of larger volumes of solvent. Based on our molecular modeling results we propose a simple structural model that describes collagen tissue as a hierarchical structure, providing a bottom-up description of elastic and viscous properties form the properties of the tissue basic building blocks.

Single molecule effects of osteogenesis imperfecta mutations in tropocollagen protein domains

Osteogenesis imperfecta (OI) is a genetic disease characterized by fragile bones, skeletal deformities and, in severe cases, prenatal death that affects more than 1 in 10,000 individuals. Here we show by full atomistic simulation in explicit solvent that OI mutations have a significant influence on the mechanical properties of single tropocollagen molecules, and that the severity of different forms of OI is directly correlated with the reduction of the mechanical stiffness of individual tropocollagen molecules. The reduction of molecular stiffness provides insight into the molecular‐scale mechanisms of the disease. The analysis of the molecular mechanisms reveals that physical parameters of side‐chain volume and hydropathy index of the mutated residue control the loss of mechanical stiffness of individual tropocollagen molecules. We propose a model that enables us to predict the loss of stiffness based on these physical characteristics of mutations. This finding provides an atomistic‐level mechanistic understanding of the role of OI mutations in defining the properties of the basic protein constituents, which could eventually lead to new strategies for diagnosis and treatment the disease. The focus on material properties and their role in genetic diseases is an important, yet so far only little explored, aspect in studying the mechanisms that lead to pathological conditions. The consideration of how material properties change in diseases could lead to a new paradigm that may expand beyond the focus on biochemical readings alone and include a characterization of material properties in diagnosis and treatment, an effort referred to as materiomics.

Hydration and distance dependence of intermolecular shearing between collagen molecules in a model microfibril.

In vertebrates, collagen tissues are the main component responsible for force transmission. In spite of the physiological importance of these phenomena, force transmission mechanisms are still not fully understood, especially at smaller scales, including in particular collagen molecules and fibrils. Here we investigate the mechanism of molecular sliding between collagen molecules within a fibril, by shearing a central molecule in a hexagonally packed bundle mimicking the collagen microfibril environment, using varied lateral distance between the molecules in both dry and solvated conditions. In vacuum, the central molecule slides under a stick-slip mechanism that is due to the characteristic surface profile of collagen molecules, enhanced by the breaking and reformation of H-bonds between neighboring collagen molecules. This mechanism is consistently observed for varied lateral separations between molecules. The high shearing force (>7 nN) found for the experimentally observed intermolecular distance (≈1.1 nm) suggests that in dry samples the fibril elongation mechanism relies almost exclusively on molecular stretching, which may explain the higher stiffnesses found in dry fibrils. When hydrated, the slip-stick behavior is observed only below 1.3 nm of lateral distance, whereas above 1.3 nm the molecule shears smoothly, showing that the water layer has a strong lubricating effect. Moreover, the average force required to shear is approximately the same in solvated as in dry conditions (≈2.5 nN), which suggests that the role of water at the intermolecular level includes the transfer of load between molecules.

Mechanical properties of physiological and pathological models of collagen peptides investigated via steered molecular dynamics simulations

In this work we used molecular simulations to investigate the elastic properties of collagen single chain and triple helix with the aim of understanding its features starting from first principles. We analysed ideal collagen peptides, homotrimeric and heterotrimeric collagen type I and pathological models of collagen. Triple helices were found much more rigid than single chains, thus enlightening the important role of interchain stabilizing forces, like hydrophobic interaction and hydrogen bonds. We obtained Youngs moduli close to 4.5GPa for the ideal model of collagen and for the physiological heterotrimer, while the physiological homotrimer presented a Youngs modulus of 2.51GPa, that can be related to a mild form of Osteogenesis Imperfecta in which only the homotrimeric form of collagen type I is produced. Otherwise, the pathological model (presenting a glycine to alanine substitution) showed an elastic modulus of 4.32GPa, thus only slightly lower than the ideal model. This suggests that this mutation only slightly affects the mechanical properties of the collagen molecule, but possibly acts on an higher scale, such as the packing of collagen fibrils.

Nanomechanics of collagen microfibrils

Collagen constitutes one third of the human proteome, providing mechanical stability, elasticity and strength to organisms and is thus the prime construction material in biology. Collagen is also the dominating material in the extracellular matrix where its stiffness controls cell differentiation, growth and pathology. We use atomistic-based hierarchical multiscale modeling to describe this complex biological material from the bottom up. This includes the use and development of large-scale computational modeling tools to investigate several aspects related to collagen-based tissues, including source of visco-elasticity and deformation mechanisms at the nanoscale level. The key innovation of this research is that until now, collagen materials have primarily been described at macroscopic scales, without explicitly understanding the mechanical contributions at the molecular and fibrillar levels. The major impact of this research will be the development of fundamental models of collagenous tissues, important to the design of new scaffolding biomaterials for regenerative medicine as well as for the understanding of collagen-related diseases.

Role of Intrafibrillar Collagen Mineralization in Defining the Compressive Properties of Nascent Bone

Bone is the sole biological material found in the human body that is able to sustain compressive loads. However, although the structure of bone is well-known (it is a natural composite of collagen protein and hydroxyapatite mineral with a complex hierarchical organization), the details about the mechanisms that govern deformation at the molecular scale under compressive loading are still not completely understood. To investigate the molecular origins of bones unique compressive properties, we perform full atomistic simulations of the three-dimensional molecular structure of a mineralized collagen fibril, focusing on the role of intrafibrillar mineral densities in dictating the mechanical performance under compressive loading. We find that as the mineral density increases, the compressive modulus of the mineralized collagen increases monotonically and well beyond that of pure collagen fibrils. These findings reveal the mechanism by which bone is able to achieve superior load bearing characteristics beyond its individual constituents. Moreover, we find that intrafibrillar mineralization leads to compressive moduli that are one order of magnitude lower than the macroscale modulus of bone, indicating that extrafibrillar mineralization is mandatory for providing the load bearing properties of bone, consistent with recent experimental observations.