Solution- and adsorbed-state structural ensembles predicted for the statherin-hydroxyapatite system

We have developed a multiscale structure prediction technique to study solution- and adsorbed-state ensembles of biomineralization proteins. The algorithm employs a Metropolis Monte Carlo-plus-minimization strategy that varies all torsional and rigid-body protein degrees of freedom. We applied the technique to fold statherin, starting from a fully extended peptide chain in solution, in the presence of hydroxyapatite (HAp) (001), (010), and (100) monoclinic crystals. Blind (unbiased) predictions capture experimentally observed macroscopic and high-resolution structural features and show minimal statherin structural change upon adsorption. The dominant structural difference between solution and adsorbed states is an experimentally observed folding event in statherin's helical binding domain. Whereas predicted statherin conformers vary slightly at three different HAp crystal faces, geometric and chemical similarities of the surfaces allow structurally promiscuous binding. Finally, we compare blind predictions with those obtained from simulation biased to satisfy all previously published solid-state NMR (ssNMR) distance and angle measurements (acquired from HAp-adsorbed statherin). Atomic clashes in these structures suggest a plausible, alternative interpretation of some ssNMR measurements as intermolecular rather than intramolecular. This work demonstrates that a combination of ssNMR and structure prediction could effectively determine high-resolution protein structures at biomineral interfaces.