Accurate Force Field Parameters and pH Resolved Surface Models for Hydroxyapatite to Understand Structure, Mechanics, Hydration, and Biological Interfaces

TL;DR

This paper presents a highly accurate, pH-resolved hydroxyapatite force field integrated into major simulation platforms, enabling detailed atomistic studies of mineral-biomolecule interfaces with experimental accuracy.

Contribution

It introduces a novel, highly accurate force field and surface models for hydroxyapatite that improve simulation realism and are compatible with multiple force fields, advancing the study of biological mineralization.

Findings

Reproduces lattice constants within 0.5% deviation

Achieves water immersion energy accuracy within 5%

Models surface interactions and dissolution with high fidelity

Abstract

Mineralization of bone and teeth involves interactions between biomolecules and hydroxyapatite. Associated complex interfaces and processes remain difficult to analyze at the 1 to 100 nm scale using current laboratory techniques, and prior models for atomistic simulations are limited in the representation of chemical bonding, surface chemistry, and interfacial interactions. This work introduces an accurate force field along with pH-resolved surface models for hydroxyapatite to represent chemical bonding, structural, surface, interfacial, and mechanical properties in quantitative agreement with experiment. The accuracy is orders of magnitude higher in comparison to earlier models to facilitate quantitative monitoring of inorganic-biological assembly. The force field is integrated into the CHARMM, AMBER, OPLS-AA, PCFF, and INTERFACE force fields to enable realistic simulations of apatite-biological systems of any composition and ionic strength. Specifically, the parameters reproduce lattice constants (<0.5% deviation), IR spectrum, cleavage energies, immersion energies in water (<5% deviation), and elastic constants (<10% deviation) of hydroxyapatite in comparison to experiment. Interactions between mineral, water, and organic compounds are represented by standard combination rules in the force field without additional adjustable parameters and shown to achieve quantitative accuracy. Surface models for common (001), (010), (020), (101) facets and nanocrystals are introduced as a function of pH on the basis of extensive experimental data. New insight into surface and immersion energies, the structure of aqueous interfaces, density profiles, and superficial dissolution is described. Mechanisms of specific binding of peptides, drugs, and mineralization can be analyzed and the force field is extensible to substituted and defective apatites as well as to other calcium phosphate phases.