Research Papers

Ab Initio Investigation of Strain Dependent Atomistic Interactions at Two Tropocollagen-Hydroxyapatite Interfaces

[+] Author and Article Information
Vikas Tomar

e-mail: tomar@purdue.edu
School of Aeronautics and Astronautics,
Purdue University,
West Lafayette, IN 47907

1Corresponding author.

Contributed by the Materials Division of ASME for publication in the JOURNAL OF ENGINEERING MATERIALS AND TECHNOLOGY. Manuscript received July 30, 2012; final manuscript received January 29, 2013; published online March 28, 2013. Assoc. Editor: Xi Chen.

J. Eng. Mater. Technol 135(2), 021015 (Mar 28, 2013) (8 pages) Paper No: MATS-12-1185; doi: 10.1115/1.4023782 History: Received July 30, 2012; Revised January 29, 2013

Tropocollagen (TC) and hydroxyapatite (HAP) interfaces are one of the main load bearing entities in bone family of materials. Atomistic interactions in such interfaces occur in a variety of chemical environments under a range of biomechanical loading conditions. It is challenging to investigate such interactions using traditional analytical or using classical molecular simulation approaches owing to their limitations in predicting bond strength change as a function of change in chemical environment. In the present work, 3D ab initio molecular dynamics simulations are used to understand such atomistic interactions by analyzing tensile strain dependent deformation mechanism and strength of two structurally distinct idealized TC-HAP interfaces in hydrated as well as unhydrated environments. Analyses suggest that the presence of water molecules leads to modification of H-bond density at the interfaces that also depends upon the level of strain. TC molecules become stiffer in the presence of water due to the presence of H-bonds. Bond forming-and-breaking cycle change as a function of H-bond density lies at the heart of TC-HAP interfacial shear deformation. Consequently, interfaces with TC molecule placed flat on the HAP crystal surface experience significantly higher shear stress during deformation in comparison to the interfaces with TC molecule placed with their axes perpendicular to the HAP surface.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Currey, J. D., 1977, “Mechanical-Properties of Mother of Pearl in Tension,” Proc. Royal Soc. London, Series B, 196(1125), pp. 443–463. [CrossRef]
Rho, J. Y., Kuhn-Spearing, L., and Zioupos, P., 1998, “Mechanical Properties and the Hierarchical Structure of Bone,” Med. Eng. Phys., 20(2), pp. 92–102. [CrossRef] [PubMed]
Weiner, S., and Wagner, H. D., 1998, “The Material Bone: Structure Mechanical Function Relations,” Ann. Rev. Mater. Sci., 28, pp. 271–298. [CrossRef]
Ji, B., and Gao, H., 2004, “Mechanical Properties of Nanostructure of Biological Materials,” J. Mech. Phys. Solids, 52, pp. 1963–2000. [CrossRef]
Jager, I., and Fratzl, P., 2000, “Mineralized Collagen Fibrils: A Mechanical Model With a Staggered Arrangement of Mineral Particles,” Biophys. J., 79(4), pp. 1737–1746. [CrossRef] [PubMed]
Fratzl, P., Gupta, H. S., Paschalis, E. P., and Roschger, P., 2004, “Structure and Mechanical Quality of the Collagen-Mineral Nano-Composite in Bone,” J. Mater. Chem., 14(14), pp. 2115–2123. [CrossRef]
Gupta, H. S., Seto, J., Wagermaier, W., Zaslansky, P., Boesecke, P., and Fratzl, P., 2006, “Cooperative Deformation of Mineral and Collagen in Bone at the Nanoscale, Proc. Natl. Acad. Sci. U.S.A., 103, pp. 17741–17746. [CrossRef] [PubMed]
Gupta, H. S., Wagermaier, W., Zickler, G. A., Aroush, D. R.-B., Funari, S. S., Roschger, P., Wagner, H.-D., and Fratzl, P., 2005, “Nanoscale Deformation Mechanisms in Bone,” Nano Lett., 5(10), pp. 2108–2111. [CrossRef] [PubMed]
Sasaki, N., and Odajima, S., 1996, “Stress-Strain Curve and Young’s Modulus of a Collagen Molecule as Determined by the X-Ray Diffraction Technique,” J. Biomech., 29(5), pp. 655–658. [CrossRef] [PubMed]
Eppell, S. J., Smith, B. N., Kahn, H., and Ballarini, R., 2005, “Nano Measurements With Micro-Devices: Mechanical Properties of Hydrated Collagen Fibrils,” J. R. Soc., Interface, 3, pp. 117–121. [CrossRef]
Hodge, A. J., and Petruska, J. A., 1963, “Recent Studies With the Electron Microscope on Ordered Aggregates of the Tropocollagen Macromolecule,” Aspects of Protein Structure. Proceedings of a Symposium, G. N.Ramachandran, ed., Academic, New York, pp. 289–300.
Fantner, G. E., Hassenkam, T., Kindt, J. H., Weaver, J. C., Birkedal, H., Pechenik, L., Cutroni, J. A., Cidade, G. A. G., Stucky, G. D., Morse, D. E., and Hansma, P. K., 2005, “Sacrificial Bonds and Hidden Length Dissipate Energy as Mineralized Fibrils Separate During Bone Fracture,” Nature Mater., 4(8), pp. 612–616. [CrossRef]
Thurner, P. J., Erickson, B., Jungmann, R., Schriock, Z., Weaver, J. C., Fantner, G. E., Schitter, G., Morse, D. E., and Hansma, P. K., 2007, “High-Speed Photography of Compressed Human Trabecular Bone Correlates Whitening to Microscopic Damage,” Eng. Fract. Mech., 74(12), pp. 1928–1941. [CrossRef]
Gao, H., 2006, “Application of Fracture Mechanics Concepts to Hierarchical Biomechanics of Bone and Bone-Like Materials,” Int. J. Fracture, 138, pp. 101–137. [CrossRef]
Ji, B. H., 2008, “A Study of the Interface Strength Between Protein and Mineral in Biological Materials,” J. Biomech., 41(2), pp. 259–266. [CrossRef] [PubMed]
Lorenzo, A. C., and Caffarena, E. R. l., 2005, “Elastic Properties, Young’s Modulus Determination and Structural Stability of the Tropocollagen Molecule: A Computational Study by Steered Molecular Dynamics,” J. Biomech., 38, pp. 1527–1533. [CrossRef] [PubMed]
Buehler, M. J., 2006, “Atomistic and Continuum Modeling of Mechanical Properties of Collagen: Elasticity, Fracture, and Self-Assembly,” J. Mater, Res., 21(8), pp. 1947–1962. [CrossRef]
Buehler, M. J., 2007, “Entropic Elasticity Controls Nanomechanics of Single Tropocollagen Molecules,” Biophys. J., 93, pp. 37–43. [CrossRef] [PubMed]
Ghosh, P., Katti, D. R., and Katti, K. S., 2007, “Mineral Proximity Influences Mechanical Response of Proteins in Biological Mineral-Protein Hybrid Systems,” Biomacromolecules, 8(3), pp. 851–856. [CrossRef] [PubMed]
Bhowmik, R., Katti, K. S., Verma, D., and Katti, D. R., 2007, “Probing Molecular Interactions in Bone Biomaterials: Through Molecular Dynamics and Fourier Transform Infrared Spectroscopy,” Mater. Sci. Eng., 27, pp. 352–371. [CrossRef]
Israelowitz, M., Rizvi, S. W. H., Kramer, J., and Schroeder, H. P. v., 2005, “Computational Modeling of Type I Collagen Fibers to Determine the Extracellular Matrix Structure of Connective Tissues,” Protein Eng. Des. Sel., 18(7), pp. 329–335. [CrossRef] [PubMed]
Buehler, M. J., 2008, “Nanomechanics of Collagen Fibrils Under Varying Cross-Link Densities: Atomistic and Continuum Studies,” J. Mech. Behav. Biomed. Mater., 1, pp. 59–67. [CrossRef] [PubMed]
Handgraaf, J. W., and Zerbetto, F., 2006, “Molecular Dynamics Study of Onset of Water Gelation Around the Collagen Triple Helix,” Proteins: Struct., Funct., Bioinf., 64(3), pp. 711–718. [CrossRef]
Zhang, D., Chippada, U., and Jordan, K., 2007, “Effect of the Structural Water on the Mechanical Properties of Collagen-Like Microfibrils: A Molecular Dynamics Study,” Ann. Biomed. Eng., 35, pp. 1216–1230. [CrossRef] [PubMed]
Radmer, R. J., and Klein, T. E., 2006, “Triple Helical Structure and Stabilization of Collagen-Like Molecules With 4(R)-Hydroxyproline in the Xaa Position,” Biophys. J., 90, pp. 578–588. [CrossRef] [PubMed]
Dubey, D. K., and Tomar, V., 2009, “Understanding the Influence of Structural Hierarchy and Its Coupling With Chemical Environment on the Strength of Idealized Tropocollagen–Hydroxyapatite Biomaterials,” J. Mech. Phys. Solids, 57(10), pp. 1702–1717. [CrossRef]
Dubey, D. K., and Tomar, V., 2009, “Role of the Nanoscale Interfacial Arrangement in Mechanical Strength of Tropocollagen-Hydroxyapatite Based Hard Biomaterials,” Acta Biomater., 5(7), pp. 2704–2716. [CrossRef] [PubMed]
Dubey, D. K., and Tomar, V., 2009, “The Effect of Tensile and Compressive Loading on the Hierarchical Strength of Idealized Tropocollagen-Hydroxyapatite Biomaterials as a Function of the Chemical Environment,” J. Phys.: Condens. Matter, 21(20), p. 205103. [CrossRef] [PubMed]
Niebur, G., 2013, available at http://www.nd.edu/~gniebur
Almora-Barrios, N., and de Leeuw, N. H., 2010, “A Density Functional Theory Study of the Interaction of Collagen Peptides With Hydroxyapatite Surfaces,” Langmuir, 26(18), pp. 14535–14542. [CrossRef] [PubMed]
Ching, W. Y., Rulis, P., and Misra, A., 2009, “Ab Initio Elastic Properties and Tensile Strength of Crystalline Hydroxyapatite,” Acta Biomater., 5(8), pp. 3067–3075. [CrossRef] [PubMed]
Ren, F., Leng, Y., Xin, R., and Ge, X., 2010, “Synthesis, Characterization and Ab Initio Simulation of Magnesium-Substituted Hydroxyapatite,” Acta Biomater., 6(7), pp. 2787–2796. [CrossRef] [PubMed]
Chappell, H., Duer, M., Groom, N., Pickard, C., and Bristowe, P., 2008, “Probing the Surface Structure of Hydroxyapatite Using NMR Spectroscopy and First Principles Calculations,” Phys. Chem. Chem. Phys., 10(4), pp. 600–606. [CrossRef] [PubMed]
Prockop, D. J., and Kivirikko, K. I., 1995, “Collagens: Molecular Biology, Diseases, and Potentials for Therapy,” Ann. Rev. Biochem., 64, pp. 403–434. [CrossRef]
Makarov, V., Pettitt, B. M., and Feig, M., 2002, “Solvation and Hydration of Proteins and Nucleic Acids: A Theoretical View of Simulations and Experiment,” Acc. Chem. Res., 35, pp. 376–384. [CrossRef] [PubMed]
Voth, G. A., 2006, “Computer Simulation of Proton Solvation and Transport in Aqueous and Biomolecular Systems,” Acc. Chem. Res., 39(2), pp. 143–150. [CrossRef] [PubMed]
Valiev, M., Bylaska, E. J., Govind, N., Kowalski, K., Straatsma, T. P., Van Dam, H. J. J., Wang, D., Nieplocha, J., Apra, E., Windus, T. L., and de Jong, W., 2010, “NWChem: A Comprehensive and Scalable Open-Source Solution for Large Scale Molecular Simulations,” Comput. Phys. Commun., 181(9), pp. 1477–1489. [CrossRef]
Qin, Z., Gautieri, A., Nair, A. K., Inbar, H., and Buehler, M. J., 2012, “Thickness of Hydroxyapatite Nanocrystal Controls Mechanical Properties of the Collagen-Hydroxyapatite Interface,” Langmuir, 28(4), pp. 1982–1992. [CrossRef] [PubMed]
Dubey, D. K., and Tomar, V., 2009, “Role of Hydroxyapatite Crystal Shape in Nanoscale Mechanical Behavior of Model Tropocollagen-Hydroxyapatite Hard Biomaterials,” Mater. Sci. Eng., C, 29(7), pp. 2133–2140. [CrossRef]
Horstemeyer, M. F., Baskes, M. I., and Plimpton, S. J., 2001, “Length Scale and Time Scale Effects on the Plastic Flow of FCC Metals,” Acta Mater., 49(20), pp. 4363–4374. [CrossRef]


Grahic Jump Location
Fig. 1

A representation of the hierarchical structure in bone, (a) microstructure of trabecular bone resembling the honeycomb structure with strut like structures, (b) one such strut called trabecula, having aligned mineralized collagen fibers, (c) fracture surface of human bone showing mineralized collagen fibrils, and (d) a schematic of staggered and layered assembly of tropocollagen molecules (wavy structures) and hydroxyapatite platelets to form a mineralized collagen fibril. Solid and dotted ellipses show two possible types of interactions between TC and HAP surface, one with TC parallel to HAP [001] surface (Transverse Interface Model) and other with TC perpendicular to HAP surface (L). Image in (a) is borrowed from [29] and (b, c) from [12], respectively.

Grahic Jump Location
Fig. 2

Snapshots showing TC-HAP interfaces, (a) with longitudinal interface (i.e., TC longitudinal axis perpendicular to HAP [001] surface), and (b) transverse interface (i.e., TC longitudinal axis parallel to HAP [001] surface). In both, transverse and longitudinal cases, supercells with water and without water are shown. The direction of loading (x direction) applied to all supercell is also shown in the transverse hydrated case.

Grahic Jump Location
Fig. 3

Snapshot of longitudinal interface model showing TC and HAP constituents in detailed atomistic representations

Grahic Jump Location
Fig. 4

Plots showing von-Mises stress, virial stress, and maximum shear stress in the case of (a) transverse unhydrated interface, (b) transverse hydrated interface, (c) longitudinal unhydrated interface, and (d) longitudinal hydrated interface. Corresponding snapshots of all interfaces at 50% strain value are depicted on top of each plot.

Grahic Jump Location
Fig. 5

Plots showing number of hydrogen bonds associated with TC molecule and HAP surface as a function of strain in the case of (a) transverse interface and (b) longitudinal interface for hydrated and unhydrated cases. The shown hydrogen bonds do no include hydrogen bond contribution from water-water interaction. On an average, transverse interface involves more hydrogen bonds compared to longitudinal interface.

Grahic Jump Location
Fig. 6

Plots showing a comparison of atomic level forces acting on TC molecule as a function of deformation strain value and with or without presence of water for (a) transverse and (b) longitudinal interfaces. In plot (a) dependent axis is a broken axis and is plotted in this manner for the sake of visual clarity in variation of forces with respect to strain values.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In