Isoexergonic Conformations of Surface-Bound Citrate Regulated Bioinspired Apatite Nanocrystal Growth

The superior biomechanical properties of bone and dentin are dictated, in part, by the unique plate-like morphology of hydroxyapatite (HAP) nanocrysals within a hierarchically assembled collagen matrix. Understanding the mechanism of crystal growth and thus morphology is important to the rational design of bioinspired apatite nanocrystals for orthopedic and dental applications. Citrate has long been proposed to modulate apatite crystal growth, but major questions exist regarding the HAP-bound citrate conformations and the identities of the interacting functional groups and HAP surface sites. Here, we conducted a comprehensive investigation of the mechanism from the angstrom to submicrometer scale by detailed correlation of the results of high-level metadynamics simulations, employing force-fields benchmarked to experiment and density functional theory calculations with the results of high resolution transmission electron microscopy, nuclear magnetic resonance spectroscopy, solution analysis, and thermogravimetric analysis. Crystal morphology changed from needle- to plate-like with increasing citrate concentration. Citrate adsorbed more strongly on the HAP (100) face than on the (001) face, thus resulting in preferential growth in the [001] direction and the plate-like morphology. Two very different bound conformations were obtained, involving interactions of either one or both terminal carboxyl groups with three or five surface calcium ions, respectively, and a hydrogen bond between the citrate hydroxyl and the HAP surface. Remarkably, despite fewer interaction sites in the single bound carboxyl conformation, the structures were isoexergonic, so both exist at equilibrium. Identification of the former conformation is significant because it allows a greater adsorption density than is traditionally assumed and can help explain concentration-dependence of citrate in modulating crystal morphology. These unique results were enabled first by the application of advanced metadynamics, a technique necessary for the accurate simulation of ionic materials but which is rarely employed in the biomaterials and biomineralization fields and second by the detailed correlation of computational, spectroscopic, and analytical results.