Region-specific role of water in collagen unwinding and assembly

Multimodal characterization of the collagen hydrogel structure and properties in response to physiologically relevant pH fluctuations

pH fluctuations within the extracellular matrix (ECM) and its principal constituent collagen, particularly in solid tumors and chronic wounds, may influence its structure and function. Whereas previous research examined the impact of pH on collagen fibrillogenesis, this study focuses on determining how pH fluctuations affect collagen hydrogels that mimic the physiological ECM. Utilizing a type I collagen hydrogel, we examined the influence of pH fluctuations on its structure, properties, and function while keeping the collagen hydrated. We show that collagen's secondary structure remains unaltered during pathologically relevant microenvironmental pH changes. By employing cryo scanning electron microscopy and artificial intelligence-assisted image analysis, we show that at physiological pH, collagen hydrogel presents densely packed, aligned, and elongated fibrils, which upon a decrease to pH 6.5, are transformed into shorter, sparser, and disoriented fibrils. The collagen possesses a higher storage modulus yet a lower permeability at pH 7 and 7.8 compared with pH 6.5 and 7.4. Exposing acidified collagen to a basic buffer reinstates its native structure and viscoelastic properties. Our study offers an innovative approach to analyze and characterize perturbations in hydrated collagen-based systems with potential implications for better understanding and combating disease progression. STATEMENT OF SIGNIFICANCE: As the main component of the extracellular matrix, collagen undergoes conformational changes associated with pH changes during disease. We analyze the impact of pH on pre-formed collagen fibers mimicking healthy tissues subjected to disease, and do not focus on the more studied fibrillogenesis process. Using cryogenic SEM, which allowed imaging close to the native state, we show that even minor fluctuations in the pH affect the collagen thickness, length, fiber alignment, and rheological properties. Following exposure to acidic pH, the collagen had short fibers, lacked orientation, and had low mechanical strength. This acidic collagen restored its original properties after returning to a neutral pH. These findings can help determine how pH changes can be modulated to restore healthy collagen properties.

Elucidating the role of water in collagen self-assembly by isotopically modulating collagen hydration

Water is known to play an important role in collagen self-assembly, but it is still largely unclear how water-collagen interactions influence the assembly process and determine the fibril network properties. Here, we use the H[Formula: see text]O/D[Formula: see text]O isotope effect on the hydrogen-bond strength in water to investigate the role of hydration in collagen self-assembly. We dissolve collagen in H[Formula: see text]O and D[Formula: see text]O and compare the growth kinetics and the structure of the collagen assemblies formed in these water isotopomers. Surprisingly, collagen assembly occurs ten times faster in D[Formula: see text]O than in H[Formula: see text]O, and collagen in D[Formula: see text]O self-assembles into much thinner fibrils, that form a more inhomogeneous and softer network, with a fourfold reduction in elastic modulus when compared to H[Formula: see text]O. Combining spectroscopic measurements with atomistic simulations, we show that collagen in D[Formula: see text]O is less hydrated than in H[Formula: see text]O. This partial dehydration lowers the enthalpic penalty for water removal and reorganization at the collagen-water interface, increasing the self-assembly rate and the number of nucleation centers, leading to thinner fibrils and a softer network. Coarse-grained simulations show that the acceleration in the initial nucleation rate can be reproduced by the enhancement of electrostatic interactions. These results show that water acts as a mediator between collagen monomers, by modulating their interactions so as to optimize the assembly process and, thus, the final network properties. We believe that isotopically modulating the hydration of proteins can be a valuable method to investigate the role of water in protein structural dynamics and protein self-assembly.

Collagen Structured Hydration

Collagen is a triple-helical protein unique to the extracellular matrix, conferring rigidity and stability to tissues such as bone and tendon. For the [(PPG)10]3 collagen-mimetic peptide at room temperature, our molecular dynamics simulations show that these properties result in a remarkably ordered first hydration layer of water molecules hydrogen bonded to the backbone carbonyl (bb-CO) oxygen atoms. This originates from the following observations. The radius of gyration attests that the PPG triplets are organized along a straight line, so that all triplets (excepting the ends) are equivalent. The solvent-accessible surface area (SASA) for the bb-CO oxygens shows a repetitive regularity for every triplet. This leads to water occupancy of the bb-CO sites following a similar regularity. In the crystal-phase X-ray data, as well as in our 100 K simulations, we observe a 0-2-1 water occupancy in the P-P-G triplet. Surprisingly, a similar (0-1.7-1) regularity is maintained in the liquid phase, in spite of the sub-nsec water exchange rates, because the bb-CO sites rarely remain vacant. The manifested ordered first-shell water molecules are expected to produce a cylindrical electrostatic potential around the peptide, to be investigated in future work.

1 H-13 C cross-polarization kinetics to probe hydration-dependent organic components of bone extracellular matrix

Bone is a living tissue made up of organic proteins, inorganic minerals, and water. The organic component of bone (mainly made up of Type-I collagen) provides flexibility and tensile strength. Solid-state nuclear magnetic resonance (ssNMR) is one of the few techniques that can provide atomic-level structural insights of such biomaterials in their native state. In the present article, we employed the variable contact time cross-polarization (¹ H-¹³ C CP) kinetics experiments to study the hydration-dependent atomic-level structural changes in the bone extracellular matrix (ECM). The natural abundant ¹³ C CP intensity of the bone ECM is measured by varying CP contact time and best fitted to the nonclassical kinetic model. Different relaxation parameters were measured by the best-fit equation corresponding to the different hydration conditions of the bone ECM. The associated changes in the measured parameters due to varying levels of hydration observed at different sites of collagen protein have provided its structural arrangements and interaction with water molecules in bone ECM. Overall, the present study reveals a better understanding of the kinetics of the organic part inside the bone ECM that will help in comprehending the disease-associated pathways.

Hydrated Collagen: Where Physical Chemistry, Medical Imaging, and Bioengineering Meet

It is well-known that collagen is the most abundant protein in the human body; however, what is not often appreciated is its fascinating physical chemistry and molecular physics. In this Perspective, we aim to expose some of the physicochemical phenomena associated with the hydration of collagen and to examine the role collagen's hydration water plays in determining its biological function as well as applications ranging from radiology to bioengineering. The main focus is on the Magic-Angle Effect, a phenomenon observed in Nuclear Magnetic Resonance (NMR) spectroscopy and Magnetic Resonance Imaging (MRI) of anisotropic collagenous tissues such as articular cartilage and tendon. While the effect has been known in NMR and MRI for decades, its exact molecular mechanism remains a topic of debate and continuing research in scientific literature. We survey some of the latest research aiming to develop a comprehensive molecular-level model of the Magic-Angle Effect. We also touch on other fields where understanding of collagen hydration is important, particularly nanomechanics and mechanobiology, biomaterials, and piezoelectric sensors.

Collagen type II-hyaluronan interactions - the effect of proline hydroxylation: a molecular dynamics study

Hyaluronan-collagen composites have been employed in numerous biomedical applications. Understanding the interactions between hyaluronan and collagen is particularly important in the context of joint cartilage function and the treatment of joint diseases. Many factors affect the affinity of collagen for hyaluronan. One of the important factors is the ratio of 3- or 4-hydroxy proline to proline residues. This article presents the results from molecular dynamics calculations of HA-collagen type II interactions with hyaluronan. The applied protocol employed docking and geometry optimization of complexes built using collagen structures with different numbers of hydroxyl groups attached to proline moieties. It was established that the hydroxyproline/proline ratio affects both structural and energetic features of the collagen-hyaluronan complex. Proline hydroxylation was found to significantly influence the number of all identified types of molecular forces, hydrophobic interactions, water bridges and hydrogen bonds, which can be formed between collagen and hyaluronan. Importantly, an increase in the hydroxyproline/proline ratio in the collagen chain increases the binding affinity for hyaluronan. This is illustrated by the linear correlation between the binding free energy and the hydroxylation degree. A comparison of the results obtained for 3 and 4 hydroxylation of proline indicates that the hydroxyl group attachment position plays a minor role in complex stabilization. However, a slightly stronger affinity was observed for 4 hydroxylation. In order to evaluate the effect of the aqueous environment on the collagen-hyaluronan complex stability, the enthalpic and entropic contributions to the free energy of solvation were analyzed.

Sequence-dependent mechanics of collagen reflect its structural and functional organization

Extracellular matrix mechanics influence diverse cellular functions, yet surprisingly little is known about the mechanical properties of their constituent collagen proteins. In particular, network-forming collagen IV, an integral component of basement membranes, has been far less studied than fibril-forming collagens. A key feature of collagen IV is the presence of interruptions in the triple-helix-defining (Gly-X-Y) sequence along its collagenous domain. Here, we used atomic force microscopy to determine the impact of sequence heterogeneity on the local flexibility of collagen IV and of the fibril-forming collagen III. Our extracted flexibility profile of collagen IV reveals that it possesses highly heterogeneous mechanics, ranging from semiflexible regions as found for fibril-forming collagens to a lengthy region of high flexibility toward its N-terminus. A simple model in which flexibility is dictated only by the presence of interruptions fit the extracted profile reasonably well, providing insight into the alignment of chains and demonstrating that interruptions, particularly when coinciding in multiple chains, significantly enhance local flexibility. To a lesser extent, sequence variations within the triple helix lead to variable flexibility, as seen along the continuously triple-helical collagen III. We found this fibril-forming collagen to possess a high-flexibility region around its matrix-metalloprotease binding site, suggesting a unique mechanical fingerprint of this region that is key for matrix remodeling. Surprisingly, proline content did not correlate with local flexibility in either collagen type. We also found that physiologically relevant changes in pH and chloride concentration did not alter the flexibility of collagen IV, indicating such environmental changes are unlikely to control its compaction during secretion. Although extracellular chloride ions play a role in triggering collagen IV network formation, they do not appear to modulate the structure of its collagenous domain.

Improvement of meat protein digestibility in infants and the elderly

Meat is a valuable protein source with a balanced composition of essential amino acids and various nutrients. This review aims to identify methods to improve digestion of meat proteins, as well as evaluate the digestive characteristics of infants and the elderly. Immature digestive conditions in infants, including a high gastric pH and low protease concentration, can hinder protein digestion, thus resulting in inhibited growth and development. Likewise, gastrointestinal (GI) tract aging and chronic health problems, including tooth loss and atrophic gastritis, can lead to reduction in protein digestion and absorption in the elderly compared with those in young adults. Moderate heating and several non-thermal technologies, such as aging, enzymatic hydrolysis, ultrasound, high-pressure processing, and pulsed electric field can alter protein structure and improve protein digestion in individuals with low digestive capacity.

Physical and chemical mechanisms of tissue optical clearing

Advanced optical methods combined with various probes pave the way toward molecular imaging within living cells. However, major challenges are associated with the need to enhance the imaging resolution even further to the subcellular level for the imaging of larger tissues, as well as for in vivo studies. High scattering and absorption of opaque tissues limit the penetration of light into deep tissues and thus the optical imaging depth. Tissue optical clearing technique provides an innovative way to perform deep-tissue imaging. Recently, various optical clearing methods have been developed, which provide tissue clearing based on similar physical principles via different chemical approaches. Here, we introduce the mechanisms of the current clearing methods from fundamental physical and chemical perspectives, including the main physical principle, refractive index matching via various chemical approaches, such as dissociation of collagen, delipidation, decalcification, dehydration, and hyperhydration, to reduce scattering, as well as decolorization to reduce absorption.

Effect of the Application of a Dehydrothermal Treatment on the Structure and the Mechanical Properties of Collagen Film

Dehydrothermal (DHT) treatment was used to improve the properties of collagen casings because of its non-cytotoxicity. Understanding the effects of DHT treatment on the structure and mechanical properties of collagen films is beneficial to developing satisfying collagen casings. Herein, DHT treatment with various temperatures (85-145 °C) and timescales (1-7 days) were investigated. It was clarified that the chemical crosslinking covalent bond between collagen molecules was formed after the DHT treatment. Crosslinking density increased with increasing DHT treatment temperatures, contributing to the increase of tensile strength up to over three times of that of the untreated collagen film. The increased crosslinking density was also found when increasing the DHT treatment time, and the maximum was obtained in 3 days. Further DHT treatment time did not change the crosslinking density. The damage in the triple helix structure and the self-assembly of collagen molecules were observed from IR and SAXS. The extent of denaturation increased with increasing DHT treatment temperature and time, although the effect of the DHT treatment time on the denaturation was more moderate. When the DHT treatment temperature was as high as 145 °C or the DHT treatment time exceeded 5 days, serious denaturation occurs, leading to the deterioration of mechanical properties.

Assessment of the denaturation of collagen protein concentrates using different techniques

The use of collagen and gelatin in the field of regenerative medicine is widely extended. However, most of the studies in this topic are focused on the scaffolds' properties, but only a few are related to the properties of the raw material used. The raw material analysis not only consists of a study of the composition, but also of the denaturation degree that can influence the processing and properties of the structure of the scaffold. Thus, the denaturation degree analysis of different collagen proteins was performed and assessed by the comparison of four different methods: differential scanning calorimetry (DSC), Fourier transform Infrared Spectroscopy (FTIR) and circular dichroism (CD) spectra and sulfhydryls content analysis. DSC measurements put forward a glass transition between 88°C and 95°C as well as from the FTIR measurements; the characteristic peaks for proteins are evidenced. However, from the sulfur content, only a small proportion of free sulfhydryls are present with respect to their total amount. In addition, CD spectra allow to estimate the secondary structure of the protein by the analysis of the α-helix and β-strand and also quantify the denaturation degree with the 'positive/negative ratio' (RPN) from the CD profiles, obtaining values in the range between 25% and 100%.

Structure and Dynamics of Collagen Hydration Water from Molecular Dynamics Simulations: Implications of Temperature and Pressure

Dynamics of water molecules in hydrated collagen plays an important role in determining the structural and functional properties of collagenous tissues. Experimental results suggest that collagen-bridging water molecules exhibit dynamic and thermodynamic properties of one-dimensional ice. However, molecular dynamics (MD) studies performed to date have failed to identify icelike water bridges. It has been hypothesized that this discrepancy is due to the experimental measurements and computational MD analysis having been performed on very different systems: complete tissues with large-scale collagen fiber assemblies and individual tropocollagen fragments, respectively. In this work, we explore ways of emulating a tissuelike macromolecular environment in MD simulations of hydrated collagen without increasing the size of the system to computationally prohibitive levels. We have investigated the effects of temperature and pressure on the dynamics of a small hydrated tropocollagen fragment. The occupancy and bond energies of interchain hydrogen bonds were relatively insensitive to temperature, suggesting that they play a key role in the stability of the collagen triple helix. The lifetimes of water bridges lengthened with decreasing temperature, but even at 280 K, no bridging water molecules exhibited icelike dynamics. We discuss the implications of these findings for the ability to emulate tissuelike conditions in hydrated collagen.

Environmentally Controlled Curvature of Single Collagen Proteins

The predominant structural protein in vertebrates is collagen, which plays a key role in extracellular matrix and connective tissue mechanics. Despite its prevalence and physical importance in biology, the mechanical properties of molecular collagen are far from established. The flexibility of its triple helix is unresolved, with descriptions from different experimental techniques ranging from flexible to semirigid. Furthermore, it is unknown how collagen type (homo- versus heterotrimeric) and source (tissue derived versus recombinant) influence flexibility. Using SmarTrace, a chain-tracing algorithm we devised, we performed statistical analysis of collagen conformations collected with atomic force microscopy to determine the protein's mechanical properties. Our results show that types I, II, and III collagens-the key fibrillar varieties-exhibit similar molecular flexibilities. However, collagen conformations are strongly modulated by salt, transitioning from compact to extended as KCl concentration increases in both neutral and acidic pH. Although analysis with a standard worm-like chain model suggests that the persistence length of collagen can attain a wide range of values within the literature range, closer inspection reveals that this modulation of collagen's conformational behavior is not due to changes in flexibility but rather arises from the induction of curvature (either intrinsic or induced by interactions with the mica surface). By modifying standard polymer theory to include innate curvature, we show that collagen behaves as an equilibrated curved worm-like chain in two dimensions. Analysis within the curved worm-like chain model shows that collagen's curvature depends strongly on pH and salt, whereas its persistence length does not. Thus, we find that triple-helical collagen is well described as semiflexible irrespective of source, type, pH, and salt environment. These results demonstrate that collagen is more flexible than its conventional description as a rigid rod, which may have implications for its cellular processing and secretion.

Optical clearing: impact of optical and dielectric properties of clearing solutions on pulmonary tissue mechanics

Optical clearing allows tissue visualization under preservation of organ integrity. Optical clearing of organs with a physiological change in three-dimensional geometry (such as the lung) would additionally allow visualization of macroscopic and microscopic tissue geometry. A prerequisite, however, is the preservation of the native tissue mechanics of the optically cleared lung tissue. We investigated the impact of optical and dielectric properties of clearing solutions on biomechanics and clearing potency in porcine tissue strips of healthy lungs. After fixation, bleaching, and rehydration, four methods of optical clearing were investigated using eight different protocols. The mechanical and optical properties of the cleared lung tissue strips were investigated by uniaxial tensile testing and by analyzing optical transparency and translucency for red, green, and blue light before, during, and after the biochemical optical clearing process. Fresh tissue strips were used as controls. Best balance between efficient clearing and preserved mechanics was found for clearing with a 1:1 mixture of dimethyl sulfoxide (DMSO) and aniline. Our findings show that 1) the degree of tissue transparency and translucency correlated with the refractive index of the clearing solution index ( r = 0.976, P = 0.0004; and r = 0.91, P = 0.0046, respectively), 2) tissue mechanics were affected by dehydration and the type of clearing solution, and 3) tissue biomechanics and geometry correlated with the dielectric constant of the clearing solution ( r = −0.98, P < 0.00001; and r = 0.69, P = 0.013, respectively). We show that the lower the dielectric constant of the clearing solutions, the larger the effect on tissue stiffness. This suggests that the dielectric constant is an important measure in determining the effect of a clearing solution on lung tissue biomechanics. Optimal tissue transparency requires complete tissue dehydration and a refractive index of 1.55 of the clearing solution.

NEW & NOTEWORTHY Investigating optical clearing in porcine lung tissue strips, we found that refractive index and dielectric constant of the clearing solution affected tissue clearing and biomechanics. By documenting the impact of the composition of the clearing solution on clearing potency and preservation of tissue mechanics, our results help to compose optimal clearing solutions. In addition, the results allow conclusions on the molecular interaction of solvents with collagen fibers in tissue, thereby consolidating existing theories about the functionality of collagen.

Collagen interactions: Drug design and delivery

Collagen is a major component in a wide range of drug delivery systems and biomaterial applications. Its basic physical and structural properties, together with its low immunogenicity and natural turnover, are keys to its biocompatibility and effectiveness. In addition to its material properties, the collagen triple-helix interacts with a large number of molecules that trigger biological events. Collagen interactions with cell surface receptors regulate many cellular processes, while interactions with other ECM components are critical for matrix structure and remodeling. Collagen also interacts with enzymes involved in its biosynthesis and degradation, including matrix metalloproteinases. Over the past decade, much information has been gained about the nature and specificity of collagen interactions with its partners. These studies have defined collagen sequences responsible for binding and the high-resolution structures of triple-helical peptides bound to its natural binding partners. Strategies to target collagen interactions are already being developed, including the use of monoclonal antibodies to interfere with collagen fibril formation and the use of triple-helical peptides to direct liposomes to melanoma cells. The molecular information about collagen interactions will further serve as a foundation for computational studies to design small molecules that can interfere with specific interactions or target tumor cells. Intelligent control of collagen biological interactions within a material context will expand the effectiveness of collagen-based drug delivery.

Elastic Energy Partitioning in DNA Deformation and Binding to Proteins

We study the elasticity of DNA based on local principal axes of bending identified from over 0.9-μs all-atom molecular dynamics simulations of DNA oligos. The calculated order parameters describe motion of DNA as an elastic rod. In 10 possible dinucleotide steps, bending about the two principal axes is anisotropic yet linearly elastic. Twist about the centroid axis is largely decoupled from bending, but DNA tends to overtwist for unbending beyond the typical range of thermal motion, which is consistent with experimentally observed twist-stretch coupling. The calculated elastic stiffness of dinucleotide steps yield sequence-dependent persistence lengths consistent with previous single-molecule experiments, which is further analyzed by performing coarse-grained simulations of DNA. Flexibility maps of oligos constructed from simulation also match with those from the precalculated stiffness of dinucleotide steps. These support the premise that base pair interaction at the dinucleotide-level is mainly responsible for the elasticity of DNA. Furthermore, we analyze 1381 crystal structures of protein-DNA complexes. In most structures, DNAs are mildly deformed and twist takes the highest portion of the total elastic energy. By contrast, in structures with the elastic energy per dinucleotide step greater than about 4.16 kBT (kBT: thermal energy), the major bending becomes dominant. The extensional energy of dinucleotide steps takes at most 35% of the total elastic energy except for structures containing highly deformed DNAs where linear elasticity breaks down. Such partitioning between different deformational modes provides quantitative insights into the conformational dynamics of DNA as well as its interaction with other molecules and surfaces.

Dynamic Water-Mediated Hydrogen Bonding in a Collagen Model Peptide

In the canonical (G-X-Y)(n) sequence of the fibrillar collagen triple helix, stabilizing direct interchain hydrogen bonding connects neighboring chains. Mutations of G can disrupt these interactions and are linked to connective tissue diseases. Here we integrate computational approaches with nuclear magnetic resonance (NMR) to obtain a dynamic view of hydrogen bonding distributions in the (POG)(4)(-)(POA)-(POG)(5) peptide, showing that the solution conformation, dynamics, and hydrogen bonding deviate from the reported X-ray crystal structure in many aspects. The simulations and NMR data provide clear evidence of inequivalent environments in the three chains. Molecular dynamics (MD) simulations indicate direct interchain hydrogen bonds in the leading chain, water bridges in the middle chain, and nonbridging waters in the trailing chain at the G → A substitution site. Theoretical calculations of NMR chemical shifts using a quantum fragmentation procedure can account for the unusual downfield NMR chemical shifts at the substitution sites and are used to assign the resonances to the individual chains. The NMR and MD data highlight the sensitivity of amide shifts to changes in the acceptor group from peptide carbonyls to water. The results are used to interpret solution NMR data for a variety of glycine substitutions and other sequence triplet interruptions to provide new connections between collagen sequences, their associated structures, dynamical behavior, and their ability to recognize collagen receptors.

Does L to D-amino acid substitution trigger helix→sheet conformations in collagen like peptides adsorbed to surfaces?

The present work reports on the structural order, self assembling behaviour and the role in adsorption to hydrophilic or hydrophobic solid surfaces of modified sequence from the triple helical peptide model of the collagenase cleavage site in type I collagen (Uniprot accession number P02452 residues from 935 to 970) using (D)Ala and (D)Ile substitutions as given in the models below: Model-1: GSOGADGPAGAOGTOGPQGIAGQRGVV GLOGQRGER. Model-2: GSOGADGP(D)AGAOGTOGPQGIAGQRGVVGLOGQRGER. Model-3: GSOGADGPAGAOGTOGPQG(D)IAGQRGVVGLOGQRGER. Collagenase is an important enzyme that plays an important role in degrading collagen in wound healing, cancer metastasis and even in embryonic development. However, the mechanism by which this degradation occurs is not completely understood. Our results show that adsorption of the peptides to the solid surfaces, specifically hydrophobic triggers a helix to beta transition with order increasing in peptide models 2 and 3. This restricts the collagenolytic behaviour of collagenase and may find application in design of peptides and peptidomimetics for enzyme-substrate interaction, specifically with reference to collagen and other extra cellular matrix proteins.

Chain registry and load-dependent conformational dynamics of collagen

Degradation of fibrillar collagen is critical for tissue maintenance. Yet, understanding collagen catabolism has been challenging partly due to a lack of atomistic picture for its load-dependent conformational dynamics, as both mechanical load and local unfolding of collagen affect its cleavage by matrix metalloproteinase (MMP). We use molecular dynamics simulation to find the most cleavage-prone arrangement of α chains in a collagen triple helix and find amino acids that modulate stability of the MMP cleavage domain depending on the chain registry within the molecule. The native-like state is mechanically inhomogeneous, where the cleavage site interfaces a stiff region and a locally unfolded and flexible region along the molecule. In contrast, a triple helix made of the stable glycine-proline-hydroxyproline motif is uniformly flexible and is dynamically stabilized by short-lived, low-occupancy hydrogen bonds. These results provide an atomistic basis for the mechanics, conformation, and stability of collagen that affect catabolism.

Multi-phonon scattering processes in one-dimensional anharmonic biological superlattices: understanding the dissipation of mechanical waves in mineralized tissues

The scattering of elastic waves in a one dimensional phononic (PnC) crystal composed of alternate collagen and hydroxy-apatite constituent layers is studied. These superlattices are metaphors for mineralized tissues present in bones and teeth. The collagen is treated as an open system elastic medium with water content which can vary depending on the level of stress applied. The open system nature of the collagen-water system leads to a non-linear stress-strain response. The finite difference time domain method is employed to investigate the propagation of non-linear mechanical waves through the superlattice. The spectral energy density method enables the calculation of the non-linear vibrational wave band structure. The non-linearity in the mechanical response of the collagen-water system enables a variety of multi-phonon scattering processes resulting in an increase in the number of channels for the dissipation of elastic waves and therefore for the dissipation of mechanical energy. These results provide an explanation for the relationship between bone fragility and decreased hydration.

The role of collagen charge clusters in the modulation of matrix metalloproteinase activity

Members of the matrix metalloproteinase (MMP) family selectively cleave collagens in vivo. Several substrate structural features that direct MMP collagenolysis have been identified. The present study evaluated the role of charged residue clusters in the regulation of MMP collagenolysis. A series of 10 triple-helical peptide (THP) substrates were constructed in which either Lys-Gly-Asp or Gly-Asp-Lys motifs replaced Gly-Pro-Hyp (where Hyp is 4-hydroxy-L-proline) repeats. The stabilities of THPs containing the two different motifs were analyzed, and kinetic parameters for substrate hydrolysis by six MMPs were determined. A general trend for virtually all enzymes was that, as Gly-Asp-Lys motifs were moved from the extreme N and C termini to the interior next to the cleavage site sequence, kcat/Km values increased. Additionally, all Gly-Asp-Lys THPs were as good or better substrates than the parent THP in which Gly-Asp-Lys was not present. In turn, the Lys-Gly-Asp THPs were also always better substrates than the parent THP, but the magnitude of the difference was considerably less compared with the Gly-Asp-Lys series. Of the MMPs tested, MMP-2 and MMP-9 most greatly favored the presence of charged residues with preference for the Gly-Asp-Lys series. Lys-Gly-(Asp/Glu) motifs are more commonly found near potential MMP cleavage sites than Gly-(Asp/Glu)-Lys motifs. As Lys-Gly-Asp is not as favored by MMPs as Gly-Asp-Lys, the Lys-Gly-Asp motif appears advantageous over the Gly-Asp-Lys motif by preventing unwanted MMP hydrolysis. More specifically, the lack of Gly-Asp-Lys clusters may diminish potential MMP-2 and MMP-9 collagenolytic activity. The present study indicates that MMPs have interactions spanning the P23-P23' subsites of collagenous substrates.

Structural basis for matrix metalloproteinase 1-catalyzed collagenolysis

The proteolysis of collagen triple-helical structure (collagenolysis) is a poorly understood yet critical physiological process. Presently, matrix metalloproteinase 1 (MMP-1) and collagen triple-helical peptide models have been utilized to characterize the events and calculate the energetics of collagenolysis via NMR spectroscopic analysis of 12 enzyme-substrate complexes. The triple-helix is bound initially by the MMP-1 hemopexin-like (HPX) domain via a four amino acid stretch (analogous to type I collagen residues 782-785). The triple-helix is then presented to the MMP-1 catalytic (CAT) domain in a distinct orientation. The HPX and CAT domains are rotated with respect to one another compared with the X-ray "closed" conformation of MMP-1. Back-rotation of the CAT and HPX domains to the X-ray closed conformation releases one chain out of the triple-helix, and this chain is properly positioned in the CAT domain active site for subsequent hydrolysis. The aforementioned steps provide a detailed, experimentally derived, and energetically favorable collagenolytic mechanism, as well as significant insight into the roles of distinct domains in extracellular protease function.

Epitaxially guided assembly of collagen layers on mica surfaces

Ordered assembly of collagen molecules on flat substrates has potential for various applications and serves as a model system for studying the assembly process. While previous studies demonstrated self-assembly of collagen on muscovite mica into highly ordered layers, the mechanism by which different conditions affect the resulting morphology remains to be elucidated. Using atomic force microscopy, we follow the assembly of collagen on muscovite mica at a concentration lower than the critical fibrillogenesis concentration in bulk. Initially, individual collagen molecules adsorb to mica and subsequently nucleate into fibrils possessing the 67 nm D-periodic bands. Emergence of fibrils aligned in parallel despite large interfibril distances agrees with an alignment mechanism guided by the underlying mica. The epitaxial growth was further confirmed by the formation of novel triangular networks of collagen fibrils on phlogopite mica, whose surface lattice is known to have a hexagonal symmetry, whereas the more widely used muscovite does not. Comparing collagen assembly on the two types of mica at different potassium concentrations revealed that potassium binds to the negatively charged mica surface and neutralizes it, thereby reducing the binding affinity of collagen and enhancing surface diffusion. These results suggest that collagen assembly on mica follows the surface adsorption, diffusion, nucleation, and growth pathway, where the growth direction is determined at the nucleation step. Comparison with other molecules that assemble similarly on mica supports generality of the proposed assembly mechanism, the knowledge of which will be useful for controlling the resulting surface morphologies.

Role of hydration in collagen recognition by bacterial adhesins

Protein-protein recognition regulates the vast majority of physiological or pathological processes. We investigated the role of hydration in collagen recognition by bacterial adhesin CNA by means of first principle molecular-dynamics samplings. Our characterization of the hydration properties of the isolated partners highlights dewetting-prone areas on the surface of CNA that closely match the key regions involved in hydrophobic intermolecular interactions upon complex formation, suggesting that the hydration state of the ligand-free CNA predisposes the protein to the collagen recognition. Moreover, hydration maps of the CNA-collagen complex reveal the presence of a number of structured water molecules that mediate intermolecular interactions at the interface between the two proteins. These hydration sites feature long residence times, significant binding free energies, and a geometrical distribution that closely resembles the hydration pattern of the isolated collagen triple helix. These findings are striking evidence that CNA recognizes the collagen triple helix as a hydrated molecule. For this structural motif, the exposure of several unsatisfied backbone carbonyl groups results in a strong interplay with the solvent, which is shown to also play a role in collagen recognition.

Role of hydration force in the self-assembly of collagens and amyloid steric zipper filaments

In protein self-assembly, types of surfaces determine the force between them. Yet the extent to which the surrounding water contributes to this force remains as a fundamental question. Here we study three self-assembling filament systems that respectively have hydrated (collagen), dry nonpolar, and dry polar (amyloid) interfaces. Using molecular dynamics simulations, we calculate and compare local hydration maps and hydration forces. We find that the primary hydration shells are formed all over the surface, regardless of the types of the underlying amino acids. The weakly oscillating hydration force arises from coalescence and depletion of hydration shells as two filaments approach, whereas local water diffusion, orientation, or hydrogen-bonding events have no direct effect. Hydration forces between hydrated, polar, and nonpolar interfaces differ in the amplitude and phase of the oscillation relative to the equilibrium surface separation. Therefore, water-mediated interactions between these protein surfaces, ranging in character from "hydrophobic" to "hydrophilic", have a common molecular origin based on the robustly formed hydration shells, which is likely applicable to a broad range of biomolecular assemblies whose interfacial geometry is similar in length scale to those of the present study.

A molecular dynamics study of the interprotein interactions in collagen fibrils

Molecular dynamics simulations of collagen are used to investigate at the atomistic level the nature of the interprotein interactions that are present within a collagen fibril, and which are responsible for the fibril's thermodynamic stability. Simulations both of a collagen fibril and of a fully solvated tropcollagen are compared in order to study the interactions that arise between the proteins upon the process of fibrillogenesis. The interactions studied include direct interprotein hydrogen bonds, water-mediated interprotein hydrogen bonds, and hydrophobic interactions. The simulations are used to quantify the number of interprotein interactions that form; to study which functional groups contribute most towards the interactions; and to study the spatial distribution of interprotein interactions throughout the fibril's D period. The processes of collagen fibrillogenesis and protein folding are then compared with each other, because these two physical processes share many similarities in concept, and the latter has been more widely studied. Molecular dynamics simulations of a bacteriophage T4 lysozyme protein, both in its native state and in and unfolded state, are used as an illustrative example of a typical protein folding process, for direct comparison with the collagen simulations.

The contribution of interchain salt bridges to triple-helical stability in collagen

Studies on collagen and collagen-like peptides suggest that triple-helical stability can vary along the amino acid chain. In this regard, it has been shown that lysine residues in the Y position and acidic residues in the X' position of (GPO)(3)GXYGX'Y'(GPO)(3) peptides lead to triple-helical structures with melting temperatures similar to (GPO)(8) (where O is hydroxyproline), which is generally regarded as the most stable collagen-like sequence of this length. This enhanced stability has been attributed to the formation of salt bridges between adjacent collagen chains. In this study, we explore the relationship between interchain salt bridge formation and triple-helical stability using detailed molecular simulations. Although our results confirm that salt bridges promote triple-helical stability, we find that not all salt bridges are created equal. In particular, lysine-glutamate salt bridges are most stabilizing when formed between residues in the middle strand (B) and the trailing strand (C), whereas lysine-aspartate salt bridges are most stabilizing when formed between residues in the leading (A) and middle (B) strand-the latter observation being consistent with recent NMR data on a heterotrimeric model peptide. Overall, we believe these data clarify the role of salt bridges in modulating triple-helical stability and can be used to guide the design of collagen-like peptides that have specific interchain interactions.

Molecular basis for optical clearing of collagenous tissues

Molecular interactions of optical clearing agents were investigated using a combination of molecular dynamics (MD) simulations and optical spectroscopy. For a series of sugar alcohols with low to high optical clearing potential, Raman spectroscopy and integrating sphere measurements were used to quantitatively characterize tissue water loss and reduction in light scattering following agent exposures. The rate of tissue water loss was found to correlate with agent optical clearing potential, but equivalent tissue optical clearing was measured in native and fixed tissue in vitro, given long-enough exposure times to the polyol series. MD simulations showed that the rate of tissue optical clearing correlated with the preferential formation of hydrogen bond bridges between agent and collagen. Hydrogen bond bridge formation disrupts the collagen hydration layer and facilitates replacement by a chemical agent to homogenize tissue refractive index. However, the reduction in tissue light scattering did not correlate with the agent index of refraction. Our results suggest that a necessary property of optical clearing agents is hyperosmolarity to tissue, but that the most effective agents with the highest rates of optical clearing are a subset with the highest collagen solubilities.

An efficient technique for adjusting and maintaining specific hydration levels in soft biological tissues in vitro

Elucidating how mechanics is affected by hydration in soft biological tissues is critical for understanding the potential effects of diseases where tissue extracellular matrix (ECM) is altered. The ability to control ECM water content is necessary for studying hydration-dependent tissue mechanics and for minimizing confounding effects caused by differences in tissue water content among specimens. In this paper, we describe an approach to adjust and maintain water content using a two-stage hydration technique, in order to overcome unique challenges faced in mechanical testing of biological tissues. Bovine aortic tissue was selected to demonstrate the approach. A liquid phase approach using PEG solutions allowed for efficient initial adjustment of tissue hydration. This was followed by a vapor phase approach using a humidity chamber for maintaining stable water content for a defined test duration of 45 min. Incubation in PEG solution brought bovine aortic tissue samples to equilibrium water content in approximately 6 h, much more efficiently than using a humidity chamber alone. Characteristic relationships between tissue water content and PEG concentration as well as relative humidity were obtained. It was found that PEG concentrations ranging from 0 to 40% had an inverse relationship with tissue water content ranging from 80 to 380%, which corresponded to relative humidities between 53 and 99%.

Stability of the hydration layer of tropocollagen: A QM study

Collagen is a triple helical protein, highly hydrated in nature. Bella and Berman (J Mol Biol 1996, 264, 734) have reported the structure of the first hydration layer. Water molecules form bridges of different length around the POG repeats and self assemble into left-handed helical water threads. To explore the stability of these specifically hydrated places, we have designed suitable QM models: each comprises a triple helix formed by 18 residues surrounded by 8 to 12 explicit waters. Two sets of amino acids were used, one standing for the core structural subunit of tropocollagen (POG-model) and one for its natural enzyme recognition sites (AAG-model). We have determined the stability order of the water binding places, the strongest being -8.1 kcal mol(-1), while the weakest -6.1 kcal mol(-1) per hydrogen bond. In X-ray structures, each triplet of tropocollagen is shielded by six to nine water molecules. Beside the mandatory six, the "surplus" three water molecules further strengthen the binding of all the others. However, the displacement of selected water molecules turns out to be energy neutral. These water binding places on the surface of the triple helix can provide explanation on how an almost liquid-like hydration environment exists between the closely packed tropocollagens (Henkelman et al., Magn Reson Med 1994, 32, 592). It seems that these water reservoirs or buffers can provide space for "hole conduction" of water molecules and thus contribute to the elasticity of collagen.

The relationship between water bridges and the polyproline II conformation: a large-scale analysis of molecular dynamics simulations and crystal structures

It has been suggested that denatured proteins are predisposed toward the left-handed polyproline II (P(II)) conformation. One possible source of P(II) stability in the denatured state is water bridges. Water bridges are networks of water molecules that link nearby hydrogen bond acceptors and/or donors on proteins. On the basis of the proposed behavior of P(II) and water bridges, the propensity of a residue to participate in water bridges should be correlated with its P(II) propensity. To test this hypothesis, we analyzed the following data sets: 2351 high-resolution crystal structures, and the native and denatured states of 188 different proteins from all-atom, explicit-solvent molecular dynamics (MD) simulations, which are part of our Dynameomics effort. We found that water bridges do not explain the high frequency of P(II) in denatured states; such bridges are less frequent around P(II) than around other conformations. Thus, this analysis casts doubt on water bridges as a dominant factor determining the residue-based P(II) propensities.

Thermodynamic selection of steric zipper patterns in the amyloid cross-beta spine

At the core of amyloid fibrils is the cross-beta spine, a long tape of beta-sheets formed by the constituent proteins. Recent high-resolution x-ray studies show that the unit of this filamentous structure is a beta-sheet bilayer with side chains within the bilayer forming a tightly interdigitating "steric zipper" interface. However, for a given peptide, different bilayer patterns are possible, and no quantitative explanation exists regarding which pattern is selected or under what condition there can be more than one pattern observed, exhibiting molecular polymorphism. We address the structural selection mechanism by performing molecular dynamics simulations to calculate the free energy of incorporating a peptide monomer into a beta-sheet bilayer. We test filaments formed by several types of peptides including GNNQQNY, NNQQ, VEALYL, KLVFFAE and STVIIE, and find that the patterns with the lowest binding free energy correspond to available atomistic structures with high accuracy. Molecular polymorphism, as exhibited by NNQQ, is likely because there are more than one most stable structures whose binding free energies differ by less than the thermal energy. Detailed analysis of individual energy terms reveals that these short peptides are not strained nor do they lose much conformational entropy upon incorporating into a beta-sheet bilayer. The selection of a bilayer pattern is determined mainly by the van der Waals and hydrophobic forces as a quantitative measure of shape complementarity among side chains between the beta-sheets. The requirement for self-complementary steric zipper formation supports that amyloid fibrils form more easily among similar or same sequences, and it also makes parallel beta-sheets generally preferred over anti-parallel ones. But the presence of charged side chains appears to kinetically drive anti-parallel beta-sheets to form at early stages of assembly, after which the bilayer formation is likely driven by energetics.

Role of side chains in collagen triple helix stabilization and partner recognition

Collagen is a widespread protein family involved in a variety of biological processes. The complexity of collagen and its fibrous nature prevent detailed investigations on the full-length protein. Reductionist approaches conducted by dissecting the protein complexity through the use of model peptides have proved to be quite effective. There are, however, several issues regarding structure-stability relationships, aggregation in higher-order assemblies, and partner recognition that are still extensively investigated. In this review, we discuss the role that side chains play in triple helix stabilization and in partner recognition. On the basis of recent literature data, we show that collagen triple helix stability is the result of the interplay of different factors. As a general trend, interactions established by amino/imino acid side chains within the triple helix scaffold effectively modulate the intrinsic residue propensity for this common structural motif. The use of peptide models has also highlighted the role that side chains play in collagen self-association and in its interactions with receptors. Valuable examples in these fields are illustrated. Finally, future actions required to obtain more detailed information on the structure and the function of this complex protein are also delineated.