Structure Prediction of Protein−Solid Surface Interactions Reveals a Molecular Recognition Motif of Statherin for Hydroxyapatite

Interaction of Statherin-Derived Peptide with the Surface of Hydroxyapatite: Perspectives Based on Molecular Dynamics Simulations

INTRODUCTION

Statherin-derived peptide (StatpSpS) has shown promise against erosive tooth wear. To elucidate its interaction with the hydroxyapatite (HAP) surface, the mechanism related to adsorption of this peptide with HAP was investigated through nanosecond-long all-atom molecular dynamics simulations.

METHODS

StatpSpS was positioned parallel to the HAP surface in 2 orientations: 1 - neutral and negative residues facing the surface and 2 - positive residues facing the surface. A system containing StatpSpS without HAP was also simulated as control. In the case of systems with HAP, both partially restrained surface and unrestrained surface were constructed. Structural analysis, interaction pattern, and binding-free energy were calculated.

RESULTS

In the peptide system without the HAP, there were some conformational changes during the simulation. In the presence of the surface, only moderate changes were observed. Many residues exhibited short and stable distances to the surface, indicating strong interaction. Specially, the residues ASP1 and SER2 have an important role to anchor the peptide to the surface, with positively charged residues, mainly arginine, playing a major role in the further stabilization of the peptide in an extended conformation, with close contacts to the HAP surface.

CONCLUSION

The interaction between StatpSpS and HAP is strong, and the negative charged residues are important to the anchoring of the peptide in the surface, but after the initial placement the peptide rearranges itself to maximize the interactions between positive charged residues.

Acquired enamel pellicle and biofilm engineering with a combination of acid-resistant proteins (CaneCPI-5, StN15, and Hemoglobin) for enhanced protection against dental caries - in vivo and in vitro investigations

OBJECTIVE

This study was designed in two-legs. In the in vivo, we explored the potential of a rinse solution containing a combination (Comb) of 0.1 mg/mL CaneCPI-5 (sugarcane-derive cystatin), 1.88 × 10- 5M StN15 (statherin-derived peptide) and 1.0 mg/mL hemoglobin (Hb) to change the protein profile of the acquired enamel pellicle(AEP) and the microbiome of the enamel biofilm. The in vitro, was designed to reveal the effects of Comb on the viability and bacterial composition of the microcosm biofilm, as well as on enamel demineralization.

MATERIALS AND METHODS

In vivo study, 10 participants rinsed (10mL,1 min) with either deionized water (H2O-control) or Comb. AEP and biofilm were collected after 2 and 3 h, respectively, after rinsing. AEP samples underwent proteomics analysis, while biofilm microbiome was assessed via 16 S-rRNA Next Generation Sequencing(NGS). In vitro study, a microcosm biofilm protocol was employed. Ninety-six enamel specimens were treated with: 1)Phosphate-Buffered Solution-PBS(negative-control), 2)0.12%Chlorhexidine, 3)500ppmNaF and 4)Comb. Resazurin, colony-forming-units(CFU) and Transversal Microradiography(TMR) were performed.

RESULTS

The proteomic results revealed higher quantity of proteins in the Comb compared to control associated with immune system response and oral microbial adhesion. Microbiome showed a significant increase in bacteria linked to a healthy microbiota, in the Comb group. In the in vitro study, Comb group was only efficient in reducing mineral-loss and lesion-depth compared to the PBS.

CONCLUSIONS

The AEP modification altered the subsequent layers, affecting the initial process of bacterial adhesion of pathogenic and commensal bacteria, as well as enamel demineralization.

CLINICAL RELEVANCE

Comb group shows promise in shaping oral health by potentially introducing innovative approaches to prevent enamel demineralization and deter tooth decay.

Nucleation Domains in Biomineralization: Biomolecular Sequence and Conformational Features

Biomolecules play a vital role in the regulation of biomineralization. However, the characteristics of practical nucleation domains are still sketchy. Herein, the effects of the representative biomolecular sequence and conformations on calcium phosphate (Ca-P) nucleation and mineralization are investigated. The results of computer simulations and experiments prove that the line in the arrangement of dual acidic/essential amino acids with a single interval (Bc (Basic) -N (Neutral) -Bc-N-Ac (Acidic)- NN-Ac-N) is most conducive to the nucleation. 2α-helix conformation can best induce Ca-P ion cluster formation and nucleation. "Ac- × × × -Bc" sequences with α-helix are found to be the features of efficient nucleation domains, in which process, molecular recognition plays a non-negligible role. It further indicates that the sequence determines the potential of nucleation/mineralization of biomolecules, and conformation determines the ability of that during functional execution. The findings will guide the synthesis of biomimetic mineralized materials with improved performance for bone repair.

Biomedical applications of solid-binding peptides and proteins

Over the past decades, solid-binding peptides (SBPs) have found multiple applications in materials science. In non-covalent surface modification strategies, solid-binding peptides are a simple and versatile tool for the immobilization of biomolecules on a vast variety of solid surfaces. Especially in physiological environments, SBPs can increase the biocompatibility of hybrid materials and offer tunable properties for the display of biomolecules with minimal impact on their functionality. All these features make SBPs attractive for the manufacturing of bioinspired materials in diagnostic and therapeutic applications. In particular, biomedical applications such as drug delivery, biosensing, and regenerative therapies have benefited from the introduction of SBPs. Here, we review recent literature on the use of solid-binding peptides and solid-binding proteins in biomedical applications. We focus on applications where modulating the interactions between solid materials and biomolecules is crucial. In this review, we describe solid-binding peptides and proteins, providing background on sequence design and binding mechanism. We then discuss their application on materials relevant for biomedicine (calcium phosphates, silicates, ice crystals, metals, plastics, and graphene). Although the limited characterization of SBPs still represents a challenge for their design and widespread application, our review shows that SBP-mediated bioconjugation can be easily introduced into complex designs and on nanomaterials with very different surface chemistries.

Gels containing statherin-derived peptide protect against enamel and dentin erosive tooth wear in vitro

The effect of gels containing a statherin-derived peptide (Stn) on the protection against enamel and dentin erosive tooth wear (ETW) in vitro was evaluated. Bovine enamel and dentin specimens were divided into 2 groups (n = 15 and 18/group for enamel and dentin, respectively) that were treated with Chitosan or Carboxymethylcellulose (CMC) gels containing Stn15pSpS at 1.88 × 10-5 M or 3.76 × 10-5 M. Chitosan or CMC gels without active ingredients served as negative controls, while chitosan gel containing 1.23% F (as NaF) and acidulated phosphate fluoride gel (1.23% F) served as positive controls. The gels were applied on the specimens for 4 min. Stimulated saliva was collected from 3 donors and used to form a 2-h acquired pellicle on the specimens. Then, the specimens were submitted to an erosive pH cycling protocol 4 times/day for 7 days (0.01 M HCl pH 2.0/45 s, artificial saliva/2 h, and artificial saliva overnight). The gels were applied again during pH cycling, 2 times/day for 4 min after the first and last erosive challenges. Enamel and dentin loss (μm) were assessed by contact profilometry. Scanning electron microscopy (SEM) was analyzed using a cold field emission. Data were analyzed by two-way ANOVA (for chitosan and CMC gels, separately) and Tukey's multiple comparison test. SEM images showed changes to enamel topography after application oft the gels containing Stn or F. Regarding CMC-based gels, for enamel, none of the treatments significantly reduced ETW in comparison with placebo; for dentin, however, gels containing Stn, regardless the concentration, significantly reduced the ETW. Moreover, Chitosan-based gels, regardless the Stn concentration, were able to protect enamel and dentin against ETW. Gels containing Stn might be a new approach to protect against ETW.

On the stability of peptide secondary structures on the TiO2 (101) anatase surface: a computational insight

The biological activity of proteins is partly due to their secondary structures and conformational states. Peptide chains are rather flexible so that finding ways inducing protein folding in a well-defined state is of great importance. Among the different constraint techniques, the interaction of proteins with inorganic surfaces is a fruitful strategy to stabilize selected folded states. Surface-induced peptide folding can have potential applications in different biomedicine areas, but it can also be of fundamental interest in prebiotic chemistry since the biological activity of a peptide can turn-on when folded in a given state. In this work, periodic quantum mechanical simulations (including implicit solvation effects) at the PBE-D2* level have been carried out to study the adsorption and the stability of the secondary structures (α-helix and β-sheet) of polypeptides with different chemical composition (i.e., polyglycine, polyalanine, polyglutamic acid, polylysine, and polyarginine) on the TiO₂ (101) anatase surface. The computational cost is reduced by applying periodic boundary conditions to both the surface and the peptides, thus obtaining full periodic polypeptide/TiO₂ surface systems. At variance with polyglycine, the interaction of the other polypeptides with the surface takes place with the lateral chain functionalities, leaving the secondary structures almost undistorted. Results indicate that the preferred conformation upon adsorption is the α-helix over the β-sheet, with the exception of the polyglutamic acid. According to the calculated adsorption energies, the affinity trend of the polypeptides with the (101) anatase surface is: polyarginine ≈ polylysine > polyglutamic acid > polyglycine ≈ polyalanine, both when adsorbed in gas phase and in presence of the implicit water solvent, which is very similar to the trend for the single amino acids. A set of implications related to the areas of surface-induced peptide folding, biomedicine and prebiotic chemistry are finally discussed.

Engineered Salivary Peptides Reduce Enamel Demineralization Provoked by Cariogenic S. mutans Biofilm

Engineering of the acquired enamel pellicle using salivary peptides has been shown to be a promising anticaries strategy. However, the mechanisms by which these peptides protect teeth against tooth decay are not fully understood. In this study, we evaluated the effect of the engineered salivary peptides DR9-DR9 and DR9-RR14 on enamel demineralization in two experimental conditions: (1) adsorbed onto the enamel surface forming the AEP, and (2) forming the AEP combined with their use to treat the biofilms 2×/day, using a validated cariogenic Streptococcus mutans in vitro biofilm model. Biofilms were grown for 144 h on enamel slabs and then collected to determine the bacterial viability (CFU/biofilm) and biofilm mass (mg protein/biofilm), and to extract cellular/extracellular proteins, which were characterized by mass spectrometry. The culture medium was changed 2×/day to fresh medium, and pH (indicator of biofilm acidogenicity) and calcium concentration (indicator of demineralization) was determined in used medium. DR9-RR14 peptide significantly reduced enamel demineralization (p < 0.0001) in both experimental conditions. However, this peptide did not have a significant effect on biofilm biomass (p > 0.05) nor did it modulate the expression of cellular and extracellular bacterial proteins involved in biofilm cariogenicity. These findings suggest that DR9-RR14 may control caries development mainly by a physicochemical mechanism.

Modulation of Streptococcus mutans Adherence to Hydroxyapatite by Engineered Salivary Peptides

Since the modification of the proteinaceous components of the Acquired Enamel Pellicle (AEP) could influence the adhesion of Streptococcus mutans, the most cariogenic bacteria, to dental surfaces, we assessed if engineered salivary peptides would affect the adherence and modulate the bacterial proteome upon adherence. Single-component AEPs were formed onto hydroxyapatite (HAp) discs by incubating them with statherin, histatin-3, DR9, DR9-DR9, DR9-RR14, RR14, and parotid saliva. Then, the discs were inoculated with S. mutans UA159 and the bacteria were allowed to adhere for 2 h, 4 h, and 8 h (n = 12/treatment/time point). The number of bacteria adhered to the HAp discs was determined at each time point and analyzed by two-way ANOVA and Bonferroni tests. Cell-wall proteins were extracted from adhered, planktonic, and inoculum (baseline) bacteria and proteome profiles were obtained after a bottom-up proteomics approach. The number of adhered bacteria significantly increased over time, being the mean values obtained at 8 h, from highest to lowest, as follows: DR9-RR14 > statherin > RR14 = DR9-DR9 > DR9 = histatin3 > saliva (p < 0.05). Treatments modulated the bacterial proteome upon adherence. The findings suggested a potential use of our engineered peptide DR9-DR9 to control S. mutans biofilm development by reducing bacterial colonization.

The restricted adhesion of bone marrow mesenchymal stem cells by stepped structures on surfaces of hydroxyapatite

Currently, many researches have developed several strategies to design the surface structures of hydroxyapatite (HA), and have proved that the surface structures are pivotal in guiding the adhesion of bone marrow mesenchymal stem cells (BMSCs) as well as subsequent cellular behaviours. Most of these strategies, such as altering roughness and constructing surface patterning of HA, involve the construction of geometric topographies at the micro/nanoscale. However, besides geometric topographies, crystal defects are also important characteristics of surface structures and would alter many local physicochemical properties, which is critical for contact between cells and bioceramic surfaces. For the practical applications of crystal defects, a major hindrance is that crystal defects are usually unstable and easily eliminated during crystallization, which limits the large-scale fabrication of materials with crystal defects. In this work, given that stepped structures contain massive stable crystal defects on their step edges and kinks, we proposed a feasible and efficient method to fabricate HA dishes with stepped structures on their surfaces. First, plate-like HA mesocrystals were prepared from CaHPO₄via topotactic transformation, and were shaped into HA dishes by vacuum-filtration. Then, a sintering process was applied to facilitate the formation of stepped structures on the surfaces. We demonstrated that the generation of stepped structures could restrict the adhesion of BMSCs and showed the restriction effect is highly correlated with the density of exposed stepped structures. This phenomenon is interesting and the construction of a cell adhesion model is robust and easy, the underlying mechanisms of which deserve further exploration. Furthermore, constructing stepped structures on surfaces may be a new useful strategy to regulate cell adhesion and could also cooperate with other methods that do not need change in the surface crystal structure.

Stepped structures largely exposed on surfaces of HA significantly restrict the adhesion of bone marrow mesenchymal stem cells.

Advanced non-fluoride approaches to dental enamel remineralization: The next level in enamel repair management

In modern dentistry, a minimally invasive management of early caries lesions or early-stage erosive tooth wear (ETW) with synthetic remineralization systems has become indispensable. In addition to fluoride, which is still the non-plus-ultra in these early caries/ETW treatments, a number of new developments are in the test phase or have already been commercialized. Some of these systems claim that they are comparable or even superior to fluoride in terms of their ability to remineralize enamel. Besides, their use can help avoid some of the risks associated with fluoride and support treatments of patients with a high risk of caries. Two individual non-fluoride systems can be distinguished; intrinsic and extrinsic remineralization approaches. Intrinsic (protein/peptide) systems adsorb to hydroxyapatite crystals/organics located within enamel prisms and accumulate endogenous calcium and phosphate ions from saliva, which ultimately leads to the re-growth of enamel crystals. Extrinsic remineralization systems function on the basis of the external (non-saliva) supply of calcium and phosphate to the crystals to be re-grown. This article, following an introduction into enamel (re)mineralization and fluoride-assisted remineralization, discusses the requirements for non-fluoride remineralization systems, particularly their mechanisms and challenges, and summarizes the findings that underpin the most promising advances in enamel remineralization therapy.

Molecular Dynamics Simulations of Phosphorylated Intrinsically Disordered Proteins: A Force Field Comparison

Phosphorylation is a common post-translational modification among intrinsically disordered proteins and regions, which helps regulate function by changing the protein conformations, dynamics, and interactions with binding partners. To fully comprehend the effects of phosphorylation, computer simulations are a helpful tool, although they are dependent on the accuracy of the force field used. Here, we compared the conformational ensembles produced by Amber ff99SB-ILDN+TIP4P-D and CHARMM36m, for four phosphorylated disordered peptides ranging in length from 14-43 residues. CHARMM36m consistently produced more compact conformations with a higher content of bends, mainly due to more stable salt bridges. Based on comparisons with experimental size estimates for the shortest and longest peptide, CHARMM36m appeared to overestimate the compactness. The difference between the force fields was largest for the peptide showing the greatest separation between positively charged and phosphorylated residues, in line with the importance of charge distribution. For this peptide, the conformational ensemble did not change significantly upon increasing the ionic strength from 0 mM to 150 mM, despite a reduction of the salt-bridging probability in the CHARMM36m simulations, implying that salt concentration has negligible effects in this study.

Structure and Dynamics Perturbations in Ubiquitin Adsorbed or Entrapped in Silica Materials Are Related to Disparate Surface Chemistries Resolved by Solid-State NMR Spectroscopy

Protein immobilization on material surfaces is emerging as a powerful tool in the design of devices and active materials for biomedical and pharmaceutical applications as well as for catalysis. Preservation of the protein's biological functionality is crucial to the design process and is dependent on the ability to maintain its structural and dynamical integrity while removed from the natural surroundings. The scientific techniques to validate the structure of immobilized proteins are scarce and usually provide limited information as a result of poor resolution. In this work, we benchmarked the ability of standard solid-state NMR techniques to resolve the effects of binding to dissimilar silica materials on a model protein. In particular, the interactions between ubiquitin and the surfaces of MCM41, SBA15, and silica formed in situ were tested for their influence on the structure and dynamics of the protein. It is shown that the protein's globular fold in the free state is only slightly perturbed in the three silica materials. Local motions on a residue level that are quenched by immobilization or, conversely, that arise from the process are also detailed. NMR measurements show that these perturbations are unique to each silica material and can serve as reporters of the characteristic surface chemistry.

Developments and Ongoing Challenges for Analysis of Surface-Bound Proteins

Proteins at surfaces and interfaces play important roles in the function and performance of materials in applications ranging from diagnostic assays to biomedical devices. To improve the performance of these materials, detailed molecular structure (conformation and orientation) along with the identity and concentrations of the surface-bound proteins on those materials must be determined. This article describes radiolabeling, surface plasmon resonance, quartz crystal microbalance with dissipation, X-ray photoelectron spectroscopy, secondary ion mass spectrometry, sum frequency generation spectroscopy, and computational techniques along with the information each technique provides for characterizing protein films. A multitechnique approach using both experimental and computation methods is required for these investigations. Although it is now possible to gain much insight into the structure of surface-bound proteins, it is still not possible to obtain the same level of structural detail about proteins on surfaces as can be obtained about proteins in crystals and solutions, especially for large, complex proteins. However, recent results have shown it is possible to obtain detailed structural information (e.g., backbone and side chain orientation) about small peptides (5-20 amino sequences) on surfaces. Current studies are extending these investigations to small proteins such as protein G B1 (∼6 kDa). Approaches for furthering the capabilities for characterizing the molecular structure of surface-bound proteins are proposed.

A whole area scanning-enabled direct-counting strategy for studying blocking efficiency in mitigating protein-solid surface binding

The study of protein-solid surface binding as well as blocking efficiency of blocking agents plays an important role in the development of high-performance immunoassays. Although conventional colorimetric based assays are widely employed to monitor protein non-specific binding on the surface of microplate wells and evaluate the performance of blocking agents, there is still a great need to develop new methods to achieve the same goal from a new perspective. In this study, an innovative whole area scanning (WAS)-enabled direct-counting strategy was developed and validated through studying the blocking efficiency of different blocking agents on the non-specific binding of streptavidin-alkaline phosphatase conjugate (Strep-ALP, a model protein) to the surface of 96-well microplates. After non-specific binding of Strep-ALP in wells with or without blocking agents' treatment and loading of ELF™ 97 phosphate (ELFP), ALP in Strep-ALP conjugates converts ELFP to water-insoluble ELF™ 97 alcohol (ELFA), which precipitates locally, self-assembles into large needle structures, and glows green fluorescence upon excitation. After quenching the reaction, WAS of the whole wells allows us to directly count the number of individual fluorescent precipitates, which can be used to calculate and compare the blocking efficiency of three commonly used blocking agents (BSA, casein, and dry milk) based on mitigating the non-specific binding of Strep-ALP. WAS-enabled counting of individual needle-type precipitates opens a new avenue to investigate protein-solid surface binding as well as the efficiency of blocking agents with high sensitivity.

A review on hydroxyapatite coatings for the biomedical applications: experimental and theoretical perspectives

The production of hydroxyapatite (HAP) composite coatings has continuously been investigated for bone tissue applications during the last few decades due to their significant bioactivity and osteoconductivity.

Molecular docking and molecular dynamics simulation studies on the adsorption/desorption behavior of bone morphogenetic protein-7 on the β-tricalcium phosphate surface

The adsorption/desorption behavior, and conformational and orientational changes of proteins on the surface of biomaterials are significant parameters for understanding how biomaterials perform their biological functions. In this study, for the first time, the interactions between BMP-7 and β-TCP (001) surface models with different ion-rich terminations (Ca-rich and P-rich) were investigated by molecular dynamics simulation (MD) and steered molecular dynamics simulation (SMD). The results indicated that BMP-7 preferentially interacts with both Ca-rich and P-rich β-TCP (001) surfaces at its wrist epitope residues with certain conformational changes, which led to more exposure of BMP-7 knuckle epitope residues to the environment and facilitation for binding to the type II receptor. Compared to the P-rich surface, it is speculated that the Ca-rich surface was more conducive to BMP-7 signal transduction since the upright orientation of the protein adsorption would lead to smaller hindrance for receptor binding. This study provided more atomistic and molecular information for better understanding the process of Ca-P surfaces affecting BMP-7 biological properties and further interpreted the osteoinductive mechanism from the perspective of growth factor adsorption. Moreover, the docking screening method adopted in this study is of guiding significance to the design and development of bioactive materials.

Solid-State NMR and MD Study of the Structure of the Statherin Mutant SNa15 on Mineral Surfaces

Elucidation of the structure and interactions of proteins at native mineral interfaces is key to understanding how biological systems regulate the formation of hard tissue structures. In addition, understanding how these same proteins interact with non-native mineral surfaces has important implications for the design of medical and dental implants, chromatographic supports, diagnostic tools, and a host of other applications. Here, we combine solid-state NMR spectroscopy, isotherm measurements, and molecular dynamics simulations to study how SNa15, a peptide derived from the hydroxyapatite (HAP) recognition domain of the biomineralization protein statherin, interacts with HAP, silica (SiO2), and titania (TiO2) mineral surfaces. Adsorption isotherms are used to characterize the binding affinity of SNa15 to HAP, SiO2, and TiO2. We also apply 1D 13C CP MAS, 1D 15N CP MAS, and 2D 13C-13C DARR experiments to SNa15 samples with uniformly 13C- and 15N-enriched residues to determine backbone and side-chain chemical shifts. Different computational tools, namely TALOS-N and molecular dynamics simulations, are used to deduce secondary structure from backbone and side-chain chemical shift data. Our results show that SNa15 adopts an α-helical conformation when adsorbed to HAP and TiO2, but the helix largely unravels upon adsorption to SiO2. Interactions with HAP are mediated in general by acidic and some basic amino acids, although the specific amino acids involved in direct surface interaction vary with surface. The integrated experimental and computational approach used in this study is able to provide high-resolution insights into adsorption of proteins on interfaces.

Computer simulations of the adsorption of an N-terminal peptide of statherin, SN15, and its mutants on hydroxyapatite surfaces

Salt-bridge adsorption of the SN15 peptide and its mutants on the HAP(001) surface.

Association between protein profile and postmortem interval in human bone remains

Proteomic techniques in bones forensic samples are increasingly, being applied. The main aim of forensic sciences is the estimation of postmortem interval. Most current techniques are useful for the first post-mortem stages. However, in the case of osseous remains, these techniques may be difficult to use due to the high level of decomposition of the sample. Our objective was to attempt to know whether there is a protein profile in human bone remains that would enable a late postmortem. interval ranging from 5 to 20 years postmortem to be estimated. A total of 40 femur bones from 40 different cadavers (data range 5-20 years) were use. Of the 275 total proteins, we excluded the circulating ones (n = 227), leaving a total of 48 proteins (29 structural and 19 functional) were found. A multiple correspondence analysis was applied on the 48 proteins. Finally selecting 32 proteins that allowed us to discriminate between the. two groups of postmortem interval. Analysis of the protein profile present in bone permits an approximation of the date of death within the studied interval, and could be used to complement other tests for estimating the postmortem interval.

Modulation of calcium oxalate dihydrate growth by phosphorylated osteopontin peptides

Osteopontin (OPN) is a significant component of kidney stone matrix and a key modulator of stone formation. Here, we investigated the effects of different phosphorylated states of a synthesized peptide of OPN (the ASARM peptide; acidic, serine- and aspartate-rich motif) on calcium oxalate dihydrate (COD) crystals, a major mineral phase of kidney stones. In vitro, phosphorylated OPN-ASARM peptides strongly inhibited COD crystal growth in solution as compared to the nonphosphorylated state, with increasing inhibitory potency correlating with the degree of peptide phosphorylation. Scanning electron microscopy revealed that the inhibition from the phosphopeptides resulted in distinctive, rosette-like crystal aggregates called spherulites. The OPN-ASARM peptides preferentially bound and specifically inhibited the {1 1 0} crystallographic faces of COD, as identified by combining atomic force microscopy and computational simulation approaches. These {1 1 0} surfaces of COD have high lattice calcium occupancy (exposure), providing preferential binding sites for the highly acidic peptides; binding and inhibition by OPN-ASARM peptides at the {1 1 0} faces led to crystal aggregation and intergrowth. The crystal spherulite formations obtained in vitro when using the most phosphorylated form of OPN-ASARM peptide at a high concentration, resembled crystal morphologies observed in vivo in a rat model of urolithiasis, in which crystal deposits in the kidney contain abundant OPN as revealed by immunogold labeling. A mechanistic model for spherulite formation is proposed based on the symmetry and crystallographic structure of COD, where the phosphate groups of OPN-ASARM bind to calcium atoms at [1 1 1] step risers on the COD {1 1 0} surface, inducing the periodic emergence of new COD crystals to form spherulites.

Controlling Redox Enzyme Orientation at Planar Electrodes

Redox enzymes, which catalyze reactions involving electron transfers in living organisms, are very promising components of biotechnological devices, and can be envisioned for sensing applications as well as for energy conversion. In this context, one of the most significant challenges is to achieve efficient direct electron transfer by tunneling between enzymes and conductive surfaces. Based on various examples of bioelectrochemical studies described in the recent literature, this review discusses the issue of enzyme immobilization at planar electrode interfaces. The fundamental importance of controlling enzyme orientation, how to obtain such orientation, and how it can be verified experimentally or by modeling are the three main directions explored. Since redox enzymes are sizable proteins with anisotropic properties, achieving their functional immobilization requires a specific and controlled orientation on the electrode surface. All the factors influenced by this orientation are described, ranging from electronic conductivity to efficiency of substrate supply. The specificities of the enzymatic molecule, surface properties, and dipole moment, which in turn influence the orientation, are introduced. Various ways of ensuring functional immobilization through tuning of both the enzyme and the electrode surface are then described. Finally, the review deals with analytical techniques that have enabled characterization and quantification of successful achievement of the desired orientation. The rich contributions of electrochemistry, spectroscopy (especially infrared spectroscopy), modeling, and microscopy are featured, along with their limitations.

A computational method for selecting short peptide sequences for inorganic material binding

Discovering or designing biofunctionalized materials with improved quality highly depends on the ability to manipulate and control the peptide-inorganic interaction. Various peptides can be used as assemblers, synthesizers, and linkers in the material syntheses. In another context, specific and selective material-binding peptides can be used as recognition blocks in mining applications. In this study, we propose a new in silico method to select short 4-mer peptides with high affinity and selectivity for a given target material. This method is illustrated with the calcite (104) surface as an example, which has been experimentally validated. A calcite binding peptide can play an important role in our understanding of biomineralization. A practical aspect of calcite is a need for it to be selectively depressed in mining sites.

Predicting the orientation of protein G B1 on hydrophobic surfaces using Monte Carlo simulations

A Monte Carlo algorithm was developed to predict the most likely orientations of protein G B1, an immunoglobulin G (IgG) antibody-binding domain of protein G, adsorbed onto a hydrophobic surface. At each Monte Carlo step, the protein was rotated and translated as a rigid body. The assumption about rigidity was supported by quartz crystal microbalance with dissipation monitoring experiments, which indicated that protein G B1 adsorbed on a polystyrene surface with its native structure conserved and showed that its IgG antibody-binding activity was retained. The Monte Carlo simulations predicted that protein G B1 is likely adsorbed onto a hydrophobic surface in two different orientations, characterized as two mutually exclusive sets of amino acids contacting the surface. This was consistent with sum frequency generation (SFG) vibrational spectroscopy results. In fact, theoretical SFG spectra calculated from an equal combination of the two predicted orientations exhibited reasonable agreement with measured spectra of protein G B1 on polystyrene surfaces. Also, in explicit solvent molecular dynamics simulations, protein G B1 maintained its predicted orientation in three out of four runs. This work shows that using a Monte Carlo approach can provide an accurate estimate of a protein orientation on a hydrophobic surface, which complements experimental surface analysis techniques and provides an initial system to study the interaction between a protein and a surface in molecular dynamics simulations.

Effects of the dielectric properties of the ceramic-solvent interface on the binding of proteins to oxide ceramics: a non-local electrostatic approach

The rapid development of nanoscience and nanotechnology has raised many fundamental questions that significantly impede progress in these fields. In particular, understanding the physicochemical processes at the interface in aqueous solvents requires the development and application of efficient and accurate methods. In the present work we evaluate the electrostatic contribution to the energy of model protein-ceramic complex formation in an aqueous solvent. We apply a non-local (NL) electrostatic approach that accounts for the effects of the short-range structure of the solvent on the electrostatic interactions of the interfacial systems. In this approach the aqueous solvent is considered as a non-ionic liquid, with the rigid and strongly correlated dipoles of the water molecules. We have found that an ordered interfacial aqueous solvent layer at the protein- and ceramic-solvent interfaces reduces the charging energy of both the ceramic and the protein in the solvent, and significantly increases the electrostatic contribution to their association into a complex. This contribution in the presented NL approach was found to be significantly shifted with respect to the classical model at any dielectric constant value of the ceramics. This implies a significant increase of the adsorption energy in the protein-ceramic complex formation for any ceramic material. We show that for several biocompatible ceramics (for example HfO2, ZrO2, and Ta2O5) the above effect predicts electrostatically induced protein-ceramic complex formation. However, in the framework of the classical continuum electrostatic model (the aqueous solvent as a uniform dielectric medium with a high dielectric constant ∼80) the above ceramics cannot be considered as suitable for electrostatically induced complex formation. Our results also show that the protein-ceramic electrostatic interactions can be strong enough to compensate for the unfavorable desolvation effect in the process of protein-ceramic complex formation.

Predicting the Structure-Activity Relationship of Hydroxyapatite-Binding Peptides by Enhanced-Sampling Molecular Simulation

Understanding the molecular structural and energetic basis of the interactions between peptides and inorganic surfaces is critical to their applications in tissue engineering and biomimetic material synthesis. Despite recent experimental progresses in the identification and functionalization of hydroxyapatite (HAP)-binding peptides, the molecular mechanisms of their interactions with HAP surfaces are yet to be explored. In particular, the traditional method of molecular dynamics (MD) simulation suffers from insufficient sampling at the peptide-inorganic interface that renders the molecular-level observation dubious. Here we demonstrate that an integrated approach combining bioinformatics, MD, and metadynamics provides a powerful tool for investigating the structure-activity relationship of HAP-binding peptides. Four low charge density peptides, previously identified by phage display, have been considered. As revealed by bioinformatics and MD, the binding conformation of the peptides is controlled by both the sequence and the amino acid composition. It was found that formation of hydrogen bonds between lysine residue and phosphate ions on the surface dictates the binding of positively charged peptide to HAP. The binding affinities of the peptides to the surface are estimated by free energy calculation using parallel-tempering metadynamics, and the results compare favorably to measurements reported in previous experimental studies. The calculation suggests that the charge density of the peptide primarily controls the binding affinity to the surface, while the backbone secondary structure that may restrain side chain orientation toward the surface plays a minor role. We also report that the application of enhanced-sampling metadynamics effects a major advantage over the steered MD method by significantly improving the reliability of binding free energy calculation. In general, our novel integration of diverse sampling techniques should contribute to the rational design of surface-recognition peptides in biomedical applications.

ProtPOS: a python package for the prediction of protein preferred orientation on a surface

Summary: Atomistic molecular dynamics simulation is a promising technique to investigate the energetics and dynamics in the protein–surface adsorption process which is of high relevance to modern biotechnological applications. To increase the chance of success in simulating the adsorption process, favorable orientations of the protein at the surface must be determined. Here, we present ProtPOS which is a lightweight and easy-to-use python package that can predict low-energy protein orientations on a surface of interest. It combines a fast conformational sampling algorithm with the energy calculation of GROMACS. The advantage of ProtPOS is it allows users to select any force fields suitable for the system at hand and provide structural output readily available for further simulation studies.

Availability and Implementation: ProtPOS is freely available for academic and non-profit uses at http://cbbio.cis.umac.mo/software/protpos

Supplementary information:Supplementary data are available at Bioinformatics online.

Contact: shirleysiu@umac.mo

Prediction and clarification of structures of (bio)molecules on surfaces

The design of future materials for biotechnological applications via deposition of molecules on surfaces will require not only exquisite control of the deposition procedure, but of equal importance will be our ability to predict the shapes and stability of individual molecules on various surfaces. Furthermore, one will need to be able to predict the structure patterns generated during the self-organization of whole layers of (bio)molecules on the surface. In this review, we present an overview over the current state of the art regarding the prediction and clarification of structures of biomolecules on surfaces using theoretical and computational methods.

Simulations of inorganic-bioorganic interfaces to discover new materials: insights, comparisons to experiment, challenges, and opportunities

Natural and man-made materials often rely on functional interfaces between inorganic and organic compounds. Examples include skeletal tissues and biominerals, drug delivery systems, catalysts, sensors, separation media, energy conversion devices, and polymer nanocomposites. Current laboratory techniques are limited to monitor and manipulate assembly on the 1 to 100 nm scale, time-consuming, and costly. Computational methods have become increasingly reliable to understand materials assembly and performance. This review explores the merit of simulations in comparison to experiment at the 1 to 100 nm scale, including connections to smaller length scales of quantum mechanics and larger length scales of coarse-grain models. First, current simulation methods, advances in the understanding of chemical bonding, in the development of force fields, and in the development of chemically realistic models are described. Then, the recognition mechanisms of biomolecules on nanostructured metals, semimetals, oxides, phosphates, carbonates, sulfides, and other inorganic materials are explained, including extensive comparisons between modeling and laboratory measurements. Depending on the substrate, the role of soft epitaxial binding mechanisms, ion pairing, hydrogen bonds, hydrophobic interactions, and conformation effects is described. Applications of the knowledge from simulation to predict binding of ligands and drug molecules to the inorganic surfaces, crystal growth and shape development, catalyst performance, as well as electrical properties at interfaces are examined. The quality of estimates from molecular dynamics and Monte Carlo simulations is validated in comparison to measurements and design rules described where available. The review further describes applications of simulation methods to polymer composite materials, surface modification of nanofillers, and interfacial interactions in building materials. The complexity of functional multiphase materials creates opportunities to further develop accurate force fields, including reactive force fields, and chemically realistic surface models, to enable materials discovery at a million times lower computational cost compared to quantum mechanical methods. The impact of modeling and simulation could further be increased by the advancement of a uniform simulation platform for organic and inorganic compounds across the periodic table and new simulation methods to evaluate system performance in silico.

Orientational switching of protein conformation as a function of nanoparticle curvature and their geometrical fitting

Among the various surface properties, nanoparticle curvature has a direct effect on the inner root of protein nanoparticle interaction. However, the orientation of adsorbed proteins onto the nanoparticle surface and its binding mechanism still remains elusive because of the lack of in-depth knowledge at the molecular level. Here, we demonstrate detail molecular insights of the orientational switching of several serum proteins as a function of nanoparticle curvature using theoretical simulation along with some experimental results. With the variation of binding stability, four distinctly different classes of orientation were observed for human serum albumin, whereas only two unique classes of conformations were observed for ubiquitin, insulin, and haemoglobin. As a general observation, our data suggested that orientations were exclusively dependent on the specific protein structure and the geometrical fitting onto the nanoparticle surface.

Controlling Nanostructures by Templated Templates: Inheriting Molecular Orientation in Binary Heterostructures

Precise preparation strategies are required to fabricate molecular nanostructures of specific arrangement. In bottom-up approaches, where nanostructures are gradually formed by piecing together individual parts to the final structure, the self-ordering mechanisms of the involved structures are utilized. In order to achieve the desired structures regarding morphology, grain size, and orientation of the individual moieties, templates can be applied, which influence the formation process of subsequent structures. However, this strategy is of limited use for complex architectures because the templates only influence the structure formation at the interface between the template and the first compound. Here, we discuss the implementation of so-called templated templates and analyze to what extent orientations of the initial layers are inherited in the top layers of another compound to enable structural control in binary heterostructures. For that purpose, we prepared crystalline templates of the organic semiconductors pentacene and perfluoropentacene in different exclusive orientations. We observe that for templates of both individual materials the molecular orientation is inherited in the top layers of the respective counterpart. This behavior is also observed for various other molecules, indicating the robustness of this approach.

A comprehensive molecular dynamics approach to protein retention modeling in ion exchange chromatography

In downstream processing, the underlying adsorption mechanism of biomolecules to adsorbent material are still subject of extensive research. One approach to more mechanistic understanding is simulating this adsorption process and hereby the possibility to identify the parameters with strongest impact. So far this method was applied with all-atom molecular dynamics simulations of two model proteins on one cation exchanger. In this work we developed a molecular dynamics tool to simulate protein-adsorber interaction for various proteins on an anion exchanger and ran gradient elution experiments to relate the simulation results to experimental data. We were able to show that simulation results yield similar results as experimental data regarding retention behavior as well as binding orientation. We could identify arginines in case of cation exchangers and aspartic acids in case of anion exchangers as major contributors to binding.

A study of phenylalanine side-chain dynamics in surface-adsorbed peptides using solid-state deuterium NMR and rotamer library statistics

Extracellular matrix proteins adsorbed onto mineral surfaces exist in a unique environment where the structure and dynamics of the protein can be altered profoundly. To further elucidate how the mineral surface impacts molecular properties, we perform a comparative study of the dynamics of nonpolar side chains within the mineral-recognition domain of the biomineralization protein salivary statherin adsorbed onto its native hydroxyapatite (HAP) mineral surface versus the dynamics displayed by the native protein in the hydrated solid state. Specifically, the dynamics of phenylalanine side chains (viz., F7 and F14) located in the surface-adsorbed 15-amino acid HAP-recognition fragment (SN15: DpSpSEEKFLRRIGRFG) are studied using deuterium magic angle spinning ((2)H MAS) line shape and spin-lattice relaxation measurements. (2)H NMR MAS spectra and T1 relaxation times obtained from the deuterated phenylalanine side chains in free and HAP-adsorbed SN15 are fitted to models where the side chains are assumed to exchange between rotameric states and where the exchange rates and a priori rotameric state populations are varied iteratively. In condensed proteins, phenylalanine side-chain dynamics are dominated by 180° flips of the phenyl ring, i.e., the "π flip". However, for both F7 and F14, the number of exchanging side-chain rotameric states increases in the HAP-bound complex relative to the unbound solid sample, indicating that increased dynamic freedom accompanies introduction of the protein into the biofilm state. The observed rotameric exchange dynamics in the HAP-bound complex are on the order of 5-6 × 10(6) s(-1), as determined from the deuterium MAS line shapes. The dynamics in the HAP-bound complex are also shown to have some solution-like behavioral characteristics, with some interesting deviations from rotameric library statistics.

Solid-state NMR studies of biomineralization peptides and proteins

Nature has evolved sophisticated strategies for engineering hard tissues through the interaction of proteins, and ultimately cells, with inorganic mineral phases. This process, called biomineralization, is how living organisms transform inorganic materials such as hydroxyapatite, calcite, and silica into highly intricate and organized structures. The remarkable material properties of shell, bone, and teeth come from the activities of proteins that function at the organic-inorganic interface. A better understanding of the biomolecular mechanisms used to promote or retard the formation of mineral-based structures could provide important design principles for the development of calcification inhibitors and promoters in orthopedics, cardiology, urology, and dentistry. With the knowledge of the structural basis for control of hard tissue growth by proteins, scientists could potentially develop materials using biomimetic principles with applications in catalysis, biosensors, electronic devices, and chromatographic separations, to name a few. Additionally, biomineralization also has potential applications in electronics, catalysis, magnetism, sensory devices, and mechanical design. Where man-made hard materials require the use of extreme temperatures, high pressure, and pH, biological organisms can accomplish these feats at ambient temperature and at physiological pH. Despite the fact that many researchers want to identify and control the structure of proteins at material and biomineral interfaces, there is a decided lack of molecular-level structure information available for proteins at biomaterial interfaces in general. In particular, this holds for mammalian proteins that directly control calcification processes in hard tissue. The most fundamental questions regarding the secondary and tertiary structures of proteins adsorbed to material surfaces, how proteins catalyze the formation of biomineral composites, or how proteins interact at biomaterial interfaces remain unanswered. This is largely due to a lack of methods capable of providing high-resolution structural information for proteins adsorbed to material surfaces under physiologically relevant conditions. In this Account, we highlight recent work that is providing insight into the structure and crystal recognition mechanisms of a salivary protein model system, as well as the structure and interactions of a peptide that catalyzes the formation of biosilica composites. To develop a better understanding of the structure and interactions of proteins in biomaterials, we have used solid-state NMR techniques to determine the molecular structure and dynamics of proteins and peptides adsorbed onto inorganic crystal surfaces and embedded within biomineral composites. This work adds to the understanding of the structure and crystal recognition mechanisms of an acidic human salivary phosphoprotein, statherin.

Acidulated phosphate fluoride application changes the protein composition of human acquired enamel pellicle

We evaluated, by proteomic analysis, whether the chemical changes provoked on enamel by acidulated phosphate fluoride (APF) application alter the protein composition of acquired enamel pellicle. Enamel slabs, pretreated with distilled water (negative control), phosphoric acid (active control) or APF solution, were immersed in human saliva for pellicle formation. The adsorbed proteins were extracted and analyzed by liquid chromatography-mass spectrometry/mass spectrometry. Fifty-six proteins were identified, 12 exclusive to APF and 11 to phosphoric acid. APF decreased the concentration of histatin-1, but increased the concentration of S100-A9, which is confirmed by immunoblotting. The findings suggest that APF application changes the acquired enamel pellicle composition.

Rapid calculation of protein pKa values using Rosetta

We developed a Rosetta-based Monte Carlo method to calculate the pK(a) values of protein residues that commonly exhibit variable protonation states (Asp, Glu, Lys, His, and Tyr). We tested the technique by calculating pK(a) values for 264 residues from 34 proteins. The standard Rosetta score function, which is independent of any environmental conditions, failed to capture pK(a) shifts. After incorporating a Coulomb electrostatic potential and optimizing the solvation reference energies for pK(a) calculations, we employed a method that allowed side-chain flexibility and achieved a root mean-square deviation (RMSD) of 0.83 from experimental values (0.68 after discounting 11 predictions with an error over 2 pH units). Additional degrees of side-chain conformational freedom for the proximal residues facilitated the capture of charge-charge interactions in a few cases, resulting in an overall RMSD of 0.85 pH units. The addition of backbone flexibility increased the overall RMSD to 0.93 pH units but improved relative pK(a) predictions for proximal catalytic residues. The method also captures large pK(a) shifts of lysine and some glutamate point mutations in staphylococcal nuclease. Thus, a simple and fast method based on the Rosetta score function and limited conformational sampling produces pK(a) values that will be useful when rapid estimation is essential, such as in docking, design, and folding.

Using the RosettaSurface algorithm to predict protein structure at mineral surfaces

Determination of protein structure on mineral surfaces is necessary to understand biomineralization processes toward better treatment of biomineralization diseases and design of novel protein-synthesized materials. To date, limited atomic-resolution data have hindered experimental structure determination for proteins on mineral surfaces. Molecular simulation represents a complementary approach. In this chapter, we review RosettaSurface, a computational structure prediction-based algorithm designed to broadly sample conformational space to identify low-energy structures. We summarize the computational approaches, the published applications, and the new releases of the code in the Rosetta 3 framework. In addition, we provide a protocol capture to demonstrate the practical steps to employ RosettaSurface. As an example, we provide input files and output data analysis for a previously unstudied mineralization protein, osteocalcin. Finally, we summarize ongoing challenges in energy function optimization and conformational searching and suggest that the fusion between experiment and calculation is the best route forward.

Docking of ubiquitin to gold nanoparticles

Protein-nanoparticle associations have important applications in nanoscience and nanotechnology such as targeted drug delivery and theranostics. However, the mechanisms by which proteins recognize nanoparticles and the determinants of specificity are still poorly understood at the microscopic level. Gold is a promising material in nanoparticles for nanobiotechnology applications because of the ease of its functionalization and its tunable optical properties. Ubiquitin is a small, cysteine-free protein (ubiquitous in eukaryotes) whose binding to gold nanoparticles has been characterized recently by nuclear magnetic resonance (NMR). To reveal the molecular basis of these protein-nanoparticle interactions, we performed simulations at multiple levels (ab initio quantum mechanics, classical molecular dynamics and Brownian dynamics) and compared the results with experimental data (circular dichroism and NMR). The results provide a model of the ensemble of structures constituting the ubiquitin-gold surface complex, and insights into the driving forces for the binding of ubiquitin to gold nanoparticles, the role of nanoparticle surfactants (citrate) in the association process, and the origin of the perturbations in the NMR chemical shifts.

Interactions binding mineral and organic phases in nanocomposites based on bacterial cellulose and calcium phosphates

The interactions responsible for the adhesion of calcium phosphate (CP) nanocrystals and bacterial cellulose (BC) nanofibrils in the composite material obtained by mixing aqueous suspensions of presynthesized CP and BC and the dependence of these interactions on the CP morphology and chemical structure have been elucidated by molecular mechanics calculations of the CP-BC interfacial structures. The interactions between the superficial CP and BC crystal layers have been simulated. Two crystalline CP structures (i.e., hydroxyapatite (HAP) and whitlockite) with two morphologies (plate-shaped and rod-shaped) were considered. Electrostatics has been found to be the major contributor to the adhesion of the CP crystallites and BC nanofibers, and the formation of interfacial hydrogen bonds makes a minor contribution to the interaction energy. It has also been found that, in general, the energy gain resulting from whitlockite-BC binding is greater than that for HAP-BC binding, and the binding of the rod-shaped crystallites of whitlockite with BC is the most profitable. The energy loss and entropy gain upon replacement of the BC-water and CP-water contacts by the BC-CP contacts have been estimated.

Does adsorption at hydroxyapatite surfaces induce peptide folding? Insights from large-scale B3LYP calculations

Large-scale periodic quantum mechanical calculations (509 atoms, 7852 atomic orbitals) based on the hybrid B3LYP functional focused on the peptide folding induced by the adsorption on the (001) and (010) hydroxyapatite (HA) surfaces give interesting insights on the role of specific interactions between surface sites and the peptide, which stabilize the helix conformation over the "native" random coil ones for in silico designed model peptides. The two peptides were derived from the 12-Gly oligomer, with one (P1, C-tGGKGGGGGGEGGN-t) and two (P2, C-tGGKGGKEGGEGGN-t) glutamic acid (E) and lysine (K) residue mutations. The most stable gas-phase "native" conformation for both peptides resulted in a random coil (RC) structure, with the helix (H) conformation being ≈100 kJ mol(-1) higher in free energy. The two peptide conformations interact with the HA (001) and (010) surfaces by C═O groups via Ca(2+) ions, by hydrogen bond between NH(2) groups and the basic PO(4)(3-) groups and by a relevant fraction due to dispersion forces. Peptide adsorption was studied on the dry (001) surface, the wet one envisaging 2 H(2)O per surface Ca(2+) and, on the latter, also considering the adsorption of microsolvated peptides with 4 H(2)O molecules located at sites responsible of the interaction with the surface. The P1 mutant does prefer to be adsorbed as a random coil by ≈160 kJ/mol, whereas the reverse is computed for P2, preferring the helix conformation by ≈50 kJ/mol. Adsorption as helix of both P1 and P2 mutants brings about proton transfer toward the HA surfaces with a large charge transfer component to the interaction energy.

Direct observation of phenylalanine orientations in statherin bound to hydroxyapatite surfaces

Extracellular biomineralization proteins such as salivary statherin control the growth of hydroxyapatite (HAP), the principal component of teeth and bones. Despite the important role that statherin plays in the regulation of hard tissue formation in humans, the surface recognition mechanisms involved are poorly understood. The protein-surface interaction likely involves very specific contacts between the surface atoms and the key protein side chains. This study demonstrates for the first time the power of combining near-edge X-ray absorption fine structure (NEXAFS) spectroscopy with element labeling to quantify the orientation of individual side chains. In this work, the 15 amino acid N-terminal binding domain of statherin has been adsorbed onto HAP surfaces, and the orientations of phenylalanine rings F7 and F14 have been determined using NEXAFS analysis and fluorine labels at individual phenylalanine sites. The NEXAFS-derived phenylalanine tilt angles have been verified with sum frequency generation spectroscopy.

Ab initio modelling of protein-biomaterial interactions: influence of amino acid polar side chains on adsorption at hydroxyapatite surfaces

The adsorption from the gas phase of five different amino acids (AAs), namely Gly, Ser, Lys, Gln and Glu, on three surface models of hexagonal hydroxyapatite (HA) has been studied at B3LYP level with Gaussian type basis set within a periodic approach. The AA adsorption was simulated on the (001) and (010) stoichiometric surfaces, the latter both in its pristine and water-reacted form. Low/high AA coverage has been studied by doubling the HA unit cell size. The AAs have been docked to the HA surfaces following the electrostatic complementarity between the electrostatic potentials of AA and the bare HA. Gly adsorbs as a zwitterion at the (001) surface, whereas at the (010) ones, the proton of the COOH group is transferred to the surface resulting in an HA(+)/Gly(-) ion pair. For the other AAs, the common COOH-CH-NH(2) moiety behaves like in Gly, while the specific side-chain functionalities adsorb at the HA surfaces by maximizing electrostatic and H-bond interactions. The interactions between the side chains and the HA surface impart a higher stability compared with the Gly case, with Glu being the strongest adsorbate owing to its high Ca affinity and H-bond donor propensity. For AAs of large size, the adsorption is more favourable in conditions of low coverage as repulsion between adjacent AAs is avoided. For all considered AAs, the strongest interaction is always established on the (010) faces rather than on the (001) one, therefore suggesting an easier growth along the c-direction of HA crystals from AA solutions.