Polymorphism of stable collagen fibrils

Development of a facile method to compute collagen network pathological anisotropy using AFM imaging

Type I collagen, a fundamental extracellular matrix (ECM) component, is pivotal in maintaining tissue integrity and strength. It is also the most prevalent fibrous biopolymer within the ECM, ubiquitous in mammalian organisms. This structural protein provides essential mechanical stability and resilience to various tissues, including tendons, ligaments, skin, bone, and dentin. Collagen has been structurally investigated for several decades, and variation to its ultrastructure by histology has been associated with several pathological conditions. The current study addresses a critical challenge in the field of collagen research by providing a novel method for studying collagen fibril morphology at the nanoscale. It offers a computational approach to quantifying collagen properties, enabling a deeper understanding of how collagen type I can be affected by pathological conditions. The application of Fast Fourier Transform (FFT) coupled with Atomic Force Microscope (AFM) imaging distinguishes not only healthy and diseased skin but also holds potential for automated diagnosis of connective tissue disorders (CTDs), contributing to both clinical diagnostics and fundamental research in this area. Here we studied the changes in the structural parameters of collagen fibrils in Ehlers Danlos Syndrome (EDS). We have used skin extracted from genetically mutant mice that exhibit EDS phenotype as our model system (Col1a1Jrt/+ mice). The collagen fibrils were analyzed by AFM based descriptive-structural parameters, coupled with a 2D Fast Fourier Transform(2D-FFT) approach that automated the analysis of AFM images. In addition, each sample was characterized based on its FFT and power spectral density. Our qualitative data showed morphological differences in collagen fibril clarity (clearness of the collagen fibril edge with their neighbouring fibri), D-banding, orientation, and linearity. We have also demonstrated that FFT could be a new tool for distinguishing healthy from tissues with CTDs by measuring the disorganization of fibrils in the matrix. We have also employed FFT to reveal the orientations of the collagen fibrils, providing clinically relevant phenotypic information on their organization and anisotropy. The result of this study can be used to develop a new automated tool for better diagnosis of CTDs.

Torsion and bistability of double-twist elastomers

We investigate the elastic properties of anisotropic elastomers with a double-twist director field, which is a model for collagen fibrils or blue phases. We observe a significant Poynting-like effect, coupling torsion (fibril twist) and extension. For freely-rotating boundary conditions, we identify a structural bistability at very small extensional strains which undergoes a saddle-node bifurcation at a critical strain - at approximately 1% strain for a parameterization appropriate for collagen fibrils. With clamped boundary conditions appropriate for many experimental setups, the bifurcation is not present. We expect significant helical shape effects when fixed torsion does not equal the equilibrium torsion of freely-rotating boundary conditions, due to residual torques.

Mineralized Collagen Fibrils: An Essential Component in Determining the Mechanical Behavior of Cortical Bone

Bone comprises mechanically different materials in a specific hierarchical structure. Mineralized collagen fibrils (MCFs), represented by tropocollagen molecules and hydroxyapatite nanocrystals, are the fundamental unit of bone. The mechanical characterization of MCFs provides the unique adaptive mechanical competence to bone to withstand mechanical load. The structural and mechanical role of MCFs is critical in the deformation mechanisms of bone and the marvelous strength and toughness possessed by bone. However, the role of MCFs in the mechanical behavior of bone across multiple length scales is not fully understood. In the present study, we shed light upon the latest progress regarding bone deformation at multiple hierarchical levels and emphasize the role of MCFs during bone deformation. We propose the concept of hierarchical deformation of bone to describe the interconnected deformation process across multiple length scales of bone under mechanical loading. Furthermore, how the deterioration of bone caused by aging and diseases impairs the hierarchical deformation process of the cortical bone is discussed. The present work expects to provide insights on the characterization of MCFs in the mechanical properties of bone and lays the framework for the understanding of the multiscale deformation mechanics of bone.

Interplay between nematic and cholesteric interactions in self-consistent field theory

Chirality is a design feature of a number of biomolecules (e.g., collagen). In these molecules, cholesteric (chiral-nematic) behavior emerges from a combination of the tendency for the biopolymers to align (nematic interactions) and for the alignment direction to change with position, rotating around an axis normal to the alignment direction. This paper presents self-consistent field theory (SCFT) of chiral-nematic polymers, which takes into account polymer flexibility and the orientational degrees of freedom of polymer segments. Using the resulting SCFT, we construct a phase diagram showing regions of stability for isotropic, nematic, and cholesteric phases. Furthermore, we find that nematic interactions can stabilize the cholesteric phase, pushing the isotropic-cholesteric phase transition to lower cholesteric interaction strength, until the isotropic-nematic-cholesteric triple point is reached.

Tropocollagen springs allow collagen fibrils to stretch elastically

The mechanical properties of connective tissues are tailored to their specific function, and changes can lead to dysfunction and pathology. In most mammalian tissues the mechanical environment is governed by the micro- and nano-scale structure of collagen and its interaction with other tissue components, however these hierarchical properties remain poorly understood. In this study we use the human cornea as a model system to characterise and quantify the dominant deformation mechanisms of connective tissue in response to cyclic loads of physiological magnitude. Synchronised biomechanical testing, x-ray scattering and 3D digital image correlation revealed the presence of two dominant mechanisms: collagen fibril elongation due to a largely elastic, spring-like straightening of tropocollagen supramolecular twist, and a more viscous straightening of fibril crimp that gradually increased over successive loading cycles. The distinct mechanical properties of the two mechanisms suggest they have separate roles in vivo. The elastic, spring-like mechanism is fast-acting and likely responds to stresses associated with the cardiac cycle, while the more viscous crimp mechanism will respond to slower processes, such as postural stresses. It is anticipated that these findings will have broad applicability to understanding the normal and pathological functioning of other connective tissues such as skin and blood vessels that exhibit both helical structures and crimp.

Statement of significance

The tropocollagen spring mechanism allows collagen fibrils from some tissues to elongate significantly under small loads, and its recent discovery has the potential to change our fundamental understanding of how tissue deforms. This time-resolved study quantifies the contribution of the spring mechanism to the local strain in stretched tissue and compares it to the contribution associated with the straightening of fibril waviness, the widely accepted primary low-load strain mechanism. The spring mechanism contributed more to the local tissue strain than fibril straightening, and was found to be elastic while fibril straightening was more viscous. The results suggest that the viscoelastic behaviour of a biomaterial is controlled, at least in part, by the relative amount of fibril-scale crimp and tropocollagen supramolecular twist.

Self-consistent field theory of chiral nematic worm-like chains

Many macromolecules of biological and technological interest are both chiral and semi-flexible. DNA and collagen are good examples. Such molecules often form chiral nematic (or cholesteric) phases, as is well-documented in collagen and chitin. This work presents a method for studying cholesteric phases in the highly successful self-consistent field theory of worm-like chains, offering a new way of studying many biologically relevant molecules. The method involves an effective Hamiltonian with a chiral term inspired by the Oseen-Frank (OF) model of liquid crystals. This method is then used to examine the formation of cholesteric phases in chiral-nematic worm-like chains as a function of polymer flexibility, as well as the optimal cholesteric pitch and distribution of polymer segment orientations. Our approach not only allows for the determination of the isotropic-cholesteric transition and segment distributions, beyond what the OF model promises, but also explicitly incorporates polymer flexibility into the study of the cholesteric phase, offering a more complete understanding of the behavior of semiflexible chiral-nematic polymers.

Collagen fibril assembly: New approaches to unanswered questions

Collagen fibrils are essential for metazoan life. They are the largest, most abundant, and most versatile protein polymers in animals, where they occur in the extracellular matrix to form the structural basis of tissues and organs. Collagen fibrils were first observed at the turn of the 20th century. During the last 40 years, the genes that encode the family of collagens have been identified, the structure of the collagen triple helix has been solved, the many enzymes involved in the post-translational modifications of collagens have been identified, mutations in the genes encoding collagen and collagen-associated proteins have been linked to heritable disorders, and changes in collagen levels have been associated with a wide range of diseases, including cancer. Yet despite extensive research, a full understanding of how cells assemble collagen fibrils remains elusive. Here, we review current models of collagen fibril self-assembly, and how cells might exert control over the self-assembly process to define the number, length and organisation of fibrils in tissues.

Equilibrium mechanisms of self-limiting assembly

Self-assembly is a ubiquitous process in synthetic and biological systems, broadly defined as the spontaneous organization of multiple subunits (e.g. macromolecules, particles) into ordered multi-unit structures. The vast majority of equilibrium assembly processes give rise to two states: one consisting of dispersed disassociated subunits, and the other, a bulk-condensed state of unlimited size. This review focuses on the more specialized class of self-limiting assembly, which describes equilibrium assembly processes resulting in finite-size structures. These systems pose a generic and basic question, how do thermodynamic processes involving non-covalent interactions between identical subunits "measure" and select the size of assembled structures? In this review, we begin with an introduction to the basic statistical mechanical framework for assembly thermodynamics, and use this to highlight the key physical ingredients that ensure equilibrium assembly will terminate at finite dimensions. Then, we introduce examples of self-limiting assembly systems, and classify them within this framework based on two broad categories: self-closing assemblies and open-boundary assemblies. These include well-known cases in biology and synthetic soft matter - micellization of amphiphiles and shell/tubule formation of tapered subunits - as well as less widely known classes of assemblies, such as short-range attractive/long-range repulsive systems and geometrically-frustrated assemblies. For each of these self-limiting mechanisms, we describe the physical mechanisms that select equilibrium assembly size, as well as potential limitations of finite-size selection. Finally, we discuss alternative mechanisms for finite-size assemblies, and draw contrasts with the size-control that these can achieve relative to self-limitation in equilibrium, single-species assemblies.

Non-equilibrium growth and twist of cross-linked collagen fibrils

The lysyl oxidase (LOX) enzyme that catalyses cross-link formation during the assembly of collagen fibrils in vivo is too large to diffuse within assembled fibrils, and so is incompatible with a fully equilibrium mechanism for fibril formation. We propose that enzymatic cross-links are formed at the fibril surface during the growth of collagen fibrils; as a consequence no significant reorientation of previously cross-linked collagen molecules occurs inside collagen fibrils during fibril growth in vivo. By imposing local equilibrium only at the fibril surface, we develop a coarse-grained quantitative model of in vivo fibril structure that incorporates a double-twist orientation of collagen molecules and a periodic D-band density modulation along the fibril axis. Radial growth is controlled by the concentration of available collagen molecules outside the fibril. In contrast with earlier equilibrium models of fibril structure, we find that all fibrils can exhibit a core-shell structure that is controlled only by the fibril radius. At small radii a core is developed with a linear double-twist structure as a function of radius. Within the core the double-twist structure is largely independent of the D-band. Within the shell at larger radii, the structure approaches a constant twist configuration that is strongly coupled with the D-band. We suggest a stable radius control mechanism that corneal fibrils can exploit near the edge of the linear core regime; while larger tendon fibrils use a cruder version of growth control that does not select a preferred radius.

Mechanical basis for fibrillar bundle morphology

Understanding the morphology of self-assembled fibrillar bundles and aggregates is relevant to a range of problems in molecular biology, supramolecular chemistry and materials science. Here, we propose a coarse-grained approach that averages over specific molecular details and yields an effective mechanical theory for the spatial complexity of self-assembling fibrillar structures that arises due to the competing effects of (the bending and twisting) elasticity of individual filaments and the adhesive interactions between them. We show that our theoretical framework accounting for this allows us to capture a number of diverse fibril morphologies observed in natural and synthetic systems, ranging from Filopodia to multi-walled carbon nanotubes, and leads to a phase diagram of possible fibril shapes. We also show how the extreme sensitivity of these morphologies can lead to spatially chaotic structures. Together, these results suggest a common mechanical basis for mesoscale fibril morphology as a function of the nanoscale mechanical properties of its filamentous constituents.

Modeling Fibrillogenesis of Collagen-Mimetic Molecules

One of the most robust examples of self-assembly in living organisms is the formation of collagen architectures. Collagen type I molecules are a crucial component of the extracellular matrix, where they self-assemble into fibrils of well-defined axial striped patterns. This striped fibrillar pattern is preserved across the animal kingdom and is important for the determination of cell phenotype, cell adhesion, and tissue regulation and signaling. The understanding of the physical processes that determine such a robust morphology of self-assembled collagen fibrils is currently almost completely missing. Here, we develop a minimal coarse-grained computational model to identify the physical principles of the assembly of collagen-mimetic molecules. We find that screened electrostatic interactions can drive the formation of collagen-like filaments of well-defined striped morphologies. The fibril axial pattern is determined solely by the distribution of charges on the molecule and is robust to the changes in protein concentration, monomer rigidity, and environmental conditions. We show that the striped fibrillar pattern cannot be easily predicted from the interactions between two monomers but is an emergent result of multibody interactions. Our results can help address collagen remodeling in diseases and aging and guide the design of collagen scaffolds for biotechnological applications.

Nematic ordering of worm-like polymers near an interface

The phase behavior of semi-flexible polymers is integral to various contexts, from materials science to biophysics, many of which utilize or require specific confinement geometries as well as the orientational behavior of the polymers. Inspired by collagen assembly, we study the orientational ordering of semi-flexible polymers, modeled as Maier-Saupe worm-like chains, using self-consistent field theory. We first examine the bulk behavior of these polymers, locating the isotropic-nematic transition and delineating the limit of stability of the isotropic and nematic phases. This effort explains how nematic ordering emerges from the isotropic phase and offers insight into how different (e.g., mono-domain vs multi-domain) nematic phases form. We then clarify the influence of planar confinement on the nematic ordering of the polymers. We find that while the presence of a single confining wall does not shift the location of nematic transition, it aligns the polymers in parallel to the wall; above the onset of nematic transition, this preference tends to propagate into the bulk phase. Introducing a second, perpendicular, wall leads to a nematic phase that is parallel to both walls, allowing the ordering direction to be uniquely set by the geometry of the experimental setup. The advantage of wall-confinement is that it can be used to propagate mono-domain nematic phases into the bulk phase.

Chiral and achiral mechanisms of self-limiting assembly of twisted bundles

A generalized theory of the self-limiting assembly of twisted bundles of filaments and columns is presented. Bundles and fibers form in a broad variety of supramolecular systems, from biological to synthetic materials. A widely-invoked mechanism to explain their finite diameter relies on chirality transfer from the molecular constituents to collective twist of the assembly, the effect of which frustrates the lateral assembly and can select equilibrium, finite diameters of bundles. In this article, the thermodynamics of twisted-bundle assembly is analyzed to understand if chirality transfer is necessary for self-limitation, or instead, if spontaneously-twisting, achiral bundles also exhibit self-limited assembly. A generalized description is invoked for the elastic costs imposed by twist for bundles of various states of intra-bundle order from nematic to crystalline, as well as a generic mechanism for generating twist, classified both by chirality but also the twist susceptibility of inter-filament alignment. The theory provides a comprehensive set of predictions for the equilibrium twist and size of bundles as a function of surface energy as well as chirality, twist susceptibility, and elasticity of bundles. Moreover, it shows that while spontaneous twist can lead to self-limitation, assembly of twisted achiral bundles can be distinguished qualitatively in terms of their range of equilibrium sizes and thermodynamic stability relative to bulk (untwisted) states.

Revealing the assembly of filamentous proteins with scanning transmission electron microscopy

Filamentous proteins are responsible for the superior mechanical strength of our cells and tissues. The remarkable mechanical properties of protein filaments are tied to their complex molecular packing structure. However, since these filaments have widths of several to tens of nanometers, it has remained challenging to quantitatively probe their molecular mass density and three-dimensional packing order. Scanning transmission electron microscopy (STEM) is a powerful tool to perform simultaneous mass and morphology measurements on filamentous proteins at high resolution, but its applicability has been greatly limited by the lack of automated image processing methods. Here, we demonstrate a semi-automated tracking algorithm that is capable of analyzing the molecular packing density of intra- and extracellular protein filaments over a broad mass range from STEM images. We prove the wide applicability of the technique by analyzing the mass densities of two cytoskeletal proteins (actin and microtubules) and of the main protein in the extracellular matrix, collagen. The high-throughput and spatial resolution of our approach allow us to quantify the internal packing of these filaments and their polymorphism by correlating mass and morphology information. Moreover, we are able to identify periodic mass variations in collagen fibrils that reveal details of their axially ordered longitudinal self-assembly. STEM-based mass mapping coupled with our tracking algorithm is therefore a powerful technique in the characterization of a wide range of biological and synthetic filaments.

Type-I collagen fibrils: From growth morphology to local order

The length of type-I collagen fibrils in solution increases through the development and progress of pointed tips appearing successively at the two ends of an axis-symmetric shaft with constant diameter. Those tips, respectively fine ([Formula: see text]) or coarse ([Formula: see text]) have opposite molecular orientations. The [Formula: see text]-pointed tips, the first to appear, are particularly remarkable as they all show, on most of their length, a common parabolic profile which stays constant during the growth. Assuming that the latter occurs by lateral accretion of individual molecules in staggered configuration, we propose to give account of this prominent morphological feature along a purely geometrical argument, the profile of a tip being linked to the shape of the trajectories followed all along the accretion process. Among several possible trajectories, Fermat spirals lead to a parabolic profile in perfect agreement with the one observed for [Formula: see text]-pointed tips. This is to be put in relation with the presence of such spirals in phyllotactic patterns which ensure the best packing efficiency in cases of axis-symmetry, which is indeed that of dense collagen fibrils. Moreover, those patterns are structured by concentric circles of dislocations, constitutive of the structure itself, whose behaviour might contribute to the mechanical properties of the fibrils.

Thermodynamic modelling of acidic collagenous solutions: from free energy contributions to phase diagrams

Tropocollagen is considered one of the main precursors in the fabrication of collagen-based biomaterials. Triple helix acidic solutions of collagen I have been shown experimentally to lead to chiral plywood architectures found in bone and "cornea" like tissues. As these plywoods are solid analogues of liquid crystal architectures, bio-inspired processing and fabrication platforms based on liquid crystal physics and thermodynamics will continue to play an essential role. For tissue engineering applications, it has been shown that dilute isotropic collagen solutions need to be flow processed first and then dehydrated. Thus, a complete fundamental understanding of the thermodynamics and free energy contributions in acidic collagen aqueous solutions is necessary to avoid expensive trial-and-error fabrication. To achieve this goal, we analyze the microscopic mechanisms of ordering and interactions in solutions of triple helix collagen, namely mixing, attraction, excluded-volume and chirality. To capture the mentioned physics, we then incorporate and integrate the Flory-Huggins, Maier-Saupe, Onsager and Frank theories. Nonetheless, they together are incapable of providing an acceptable mesophasic description in acidic collagenous solutions because tropocollagen biomacromolecules are positively charged. We then explore a simple and accurate electrostatic mean-field potential. Our results on collagen are in good agreement with experiments and include phase diagrams, phase transition thresholds, and critical isotropic/cholesteric order parameters. The present extended theory is shown to properly converge to classical liquid crystal models and is used to express the phenomenological Landau-de Gennes parameters with more fundamental quantities. This study provides a platform to derive accurate process models for the fabrication of collagen-based materials, considering and benefitting from the full range of underlying interactions.