Nucleation of Biomimetic Hydroxyapatite Nanoparticles on the Surface

Jan 4, 2019 - Biomineralization is one of the most widespread phenomena in biological processes, which can facilitate to harden or stiffen existing ti...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF LOUISIANA

C: Physical Processes in Nanomaterials and Nanostructures

Nucleation of Biomimetic Hydroxyapatite Nanoparticles on the Surface of I-type Collagen: Molecular Dynamics Investigations Zhi-yu Xue, Mingli Yang, and Dingguo Xu J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b10342 • Publication Date (Web): 04 Jan 2019 Downloaded from http://pubs.acs.org on January 8, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Nucleation of Biomimetic Hydroxyapatite Nanoparticles on the Surface of I-type Collagen: Molecular Dynamics Investigations Zhiyu Xue1, Mingli Yang2,3 and Dingguo Xu1,2* 1MOE

Key Laboratory of Green Chemistry and Technology, College of Chemistry,

Sichuan University, Chengdu, Sichuan, P. R. China 2Genome

Research Center of Biomaterial, Sichuan University, Chengdu, Sichuan

610064, P. R. China 3Institute

of Atomic and Molecular Physics, MOE Key Laboratory of High Energy

Density Physics and Technology, Sichuan University, Chengdu, Sichuan, 610065, P. R. China *

To

whom

correspondence

should

be

addressed:

Tel:86-28-85406156

1

ACS Paragon Plus Environment

[email protected]

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Abstract Biomineralization is one of the most widespread phenomena in biological processes, which can facilitate to harden or stiffen existing tissues in living organism via forming minerals. Understanding the mechanism of biologically controlled or induced mineralization would be particularly useful in the development or optimization of biomaterials. The nucleation of minerals stays at the core position during mineralization. In this work, the molecular dynamics simulation approach was employed to tackle the nucleation process of hydroxyapatite in aqueous phase with the type-I collagen as the seed. In our simulation, we can offer strong evidence to confirm the formation of Ca-P clusters on the collagen surface, and predominantly around the charged residues, especially when the Glu-Arg pair occurring at the same region. Before the nucleation on the collagen surface, small Ca-P clusters could be generated initially in solution with characteristics of Ponser's Cluster. Moreover, the surface free energy calculations could confirm that the energy would continuously decrease with the increasing of cluster size. Our simulation also indicates that different collagen properties could lead to different morphology of final amorphous calcium phosphate clusters.

2

ACS Paragon Plus Environment

Page 2 of 40

Page 3 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

1. Introduction Biomineralization is one of the most widespread phenomena in biological processes, which can facilitate to harden or stiffen existing tissues in living organism via forming minerals1-6. Generally, more than 60 different minerals have been identified in organisms. Among them, minerals based on phosphate and carbonate salts of calcium are widely distributed in vertebrates, which are further in conjunction with some organic polymers like collagen and chitin to support bone structures. In fact, the bone tissue engineering involves intensive mineral growth during the bone repairing or regeneration7. Therefore, understanding mechanism of biologically controlled or induced mineralization would be particularly useful in the development or optimization of biomaterials. The nucleation is often thought to be the key step in biomineralization1,

8-10.

Currently, two fundamental nucleation theories in solution have been proposed, namely non-classical and classical nucleation theories.11-13 Quite interestingly, the recent in situ transmission electron microscopy (TEM) imaging14 of the calcium carbonate nucleation revealed that two possible nucleation pathways might exist together. The nucleation process or regulation of morphology of minerals with the help of versatile biological systems seems to be more complicated. It has been suggested that the nucleation and crystallization of hydroxyapatite in physiological environment can be initiated or controlled by proteins like collagens.15, 16 In particular, for the growth of calcium phosphate or hydroxyapatite (HAP) in bone or enamel, several proteins have been found to be important, e.g., type I collagen17, bone saloprotein18, osteopontin,19,

20

and ameloglenin.21 Therefore,

understanding the mechanism of nucleation of calcium phosphate on protein surface is not just of theoretical importance, but also has possible applications in bone repair 3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 40

engineering and the treatment of pathological mineral formation in biological environments. In this work, we will focus on the functional role of human type-I collagen protein in the nucleation of HA. Basically, the collagen is one of the major organic components of the bone extracellular matrix (ECM).22,

23

Their important

roles in bone repairing24-28 and osteoinduction29 have long been recognized. Among collagens, type-I collagen has the highest content in the human body,17, 30 which is also one of excellent biomaterials due to its low immunity and good biocompatibility. The type-I collagen is a 300 nm heterotrimer triple-helical molecule consisting of two α1 chains and one α2 chain.17, 31-33 There is another homotrimeric isoform of type I collagen, which consists of three same α1 chains.34, 35 Numerous experimental studies6, 9, 10, 15, 16, 24, 36, 37 and theoretical simulations38-44 have been devoted to understand HAP mineralization mechanisms induced or controlled by collagen proteins. Traditionally, the amorphous calcium phosphate (ACP), has been widely accepted as the precursor phase for final crystallization of hydroxyapatite.45 Moreover, with the help of biological additives,46-49 citrate50, 51 and supramoleclar template,52 the enamel- or bone-like apatites could be produced by nanodimensional hydroxyapatite and ACP. A recent cryogenic transmission electron microscopy (cryoTEM) study53 on the biomimetic surface induced calcium phosphate crystallization clearly confirms the existence of the ACP precursor phase. On the other hand, the functional roles of collagen-I in the nucleation, growth, formation of bone apatite have been studied to some extent.9, 15, 16 In particular, the nucleation sites on the collagen-I seems to be charged residues since they can interact with either Ca2+ or PO43- directly.9 Theoretically, the nucleation process of calcium phosphate regulated by collagen or model peptide has been studied via molecular dynamics (MD) simulations by several groups.38-40,

44

In particular, models using the collagen 4

ACS Paragon Plus Environment

Page 5 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

mimetic peptide (CMP), which just contains the simple repeat collagen sequence (Gly-Pro-Hyp)n, were applied to mimic the nucleation process in the presence of collagen38, 44. The abbreviation of Hyp means hydroxyproline residue. Of course, the amino acid sequences for those CMP models are much different from human type-I collagen, which is thought to be dominant in bone regeneration or repair engineering7. Therefore, simulations based on CMP models might not provide complete information of the nucleation occurring on the surface of the human type-I collagen. In addition, the type-I “super-twisted” collagen model, which is one of the collagen microfibers structure formed by an electronic

interactions of each single collagen

molecule, was employed by Zhou et al.43 However, some details of the nucleation have not be well addressed, e.g., nucleation sites, the morphology of clusters and the roles of specific charged residues. In this work, we will try to understand the nucleation mechanism of the phosphate salt of calcium with the heterotrimer human type-I collagen participating as the seed. In particular, extensive MD simulations will be applied to understand the formation and growth of ACP clusters. 2. Computational Details 2.1. Collagen Model In general, the type I collagen protein consists of three α-chains and each chains have around 300 Gly-Xaa-Yaa triplets.31 Gly is a fixed amino acid in the smallest unit and some charged residues but not proline or hydroxyproline could occupy the Xaa and Yaa positions. In vertebrate extracellular matrices, the X-ray structure (PDB code 3HR2) with the full length of collagen type I protein was reported.17 However, this structure only contains backbone α-carbon atoms. In order to construct the initial triple helix collagen model, the THeBuScr script54 was applied on the basis of the 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

primary sequence of human type-I collagen (PubMed entry number NP_000079 for α1 chain and NP_000080 for α2 chain). Considering the length of the whole collagen protein, expensive computational cost might restrict it to be treated using all atom molecular dynamics simulation. We then divided the whole sequences into thirteen even short segments. To avoid the possible boundary effects, followed by the strategy proposed by Chang et al.,34 the same six residues of GFPGPK were added to N termini and other same six residues of GEQGPA were added to cap the C termini of three chains, respectively. Once the triple helix structure was obtained, the side chains of the collagen model was complemented using Modeller (version 9.17).55 Subsequently, all thirteen short collagen models were fully relaxed by geometry optimization and molecular dynamics simulation (40 ns) using GROMACS version 4.6.7 package.56 Each model collagen was fully solvated in the pre-equilibrium rectangular TIP3P57 water box with the size of 290 Å × 80 Å × 90 Å. The non-bond interactions cutoff was set to be 12 Å. Sodium and chloride ions were added to neutralize the target system. The integration step size was set to be 2 fs, pressure at 1atm and temperature at 310K. The particle mesh Ewald (PME) algorithm58 was applied to treat the long range electrostatic interactions. In all calculations, the CHARMM2759 force fields were employed to describe the whole collagen protein. It should be pointed out that the collagen protein features a non-standard residue of hydroxyproline (Hyp). We then used the modified CHARMM27 fore fields which include the charges60, bond, angle and dihedral parameters of Hyp based on the standard AMBER general atom fore filed (GAFF).38 2.2. Nucleation MD In this work, we will try to understand the nucleation of the salt of calcium phosphate on the surface of collagen protein. Therefore, ions of calcium and 6

ACS Paragon Plus Environment

Page 6 of 40

Page 7 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

phosphate should be added into the solution properly. The normal concentration of calcium ions is 2.5-2.75 mM in our body fluid.40,

61

However, much higher

concentration (0.7 M) was employed by several groups to enhance sampling and accelerate the simulation.38-40,

62

One of example structures for collagen models

complexed with ions was displayed in Figure 1. In our simulation, the 0.7 M calcium ions were randomly added into the system, and the concentration of other ions like phosphate and hydroxyl ions followed the stoichiometry of the HAP crystal. 0.1M sodium chloride was added to neutralize the system, which mimics the human body physiological environment. The size of unit periodic box is around 290 Å × 80 Å × 90 Å for each short collagen model. This could result in a total atom number of the whole simulated systems is over 2 million. After we obtained the collagen models mixed with calcium and phosphate ions, the LAMMPS suite of program63 was employed to understand the nucleation process regulated by collagen protein. Molecular configurations were visualized using VMD 1.9.164. First of all, total of 10000 step conjugate gradient minimization was applied for geometry relaxation. MD simulations were carried out in the isothermal-isobaric ensembles (NPT) firstly (0.1 nanosecond) and followed by the 10 nanoseconds canonical (NVT) ensemble simulations using Nose-Hoover thermostats. The long-range electrostatic interactions were treated using the PME algorithm65 with precision of 10-4.

The temperature is set to be 310 K and the pressure is 1 atm,

respectively. SHAKE algorithm66,

67

was used to maintain all covalent bonds

involving hydrogen atoms of water molecules. Newton’s equations of atomic motion were integrated with 2 fs time step by the Verlet algorithm. In this work, the recent developed INTERFACE force field, namely IFF,68 was used to describe motions of calcium, phosphate and hydroxyl ions. 7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

3. Results 3.1 Precipitation of Clusters First of all, due to expensive computational cost of all atom MD simulation on the full size I-type collagen, we have simply divided it into 13 segments. To avoid the artificial effects, we then added the same terminal peptides to the terminus for each segment. During the initial stage of MD simulations of the collagen models, all 13 models can maintain their triple helix structure quite well. All nucleation dynamics will then carried out based on these models. Ten nanoseconds MD simulation is sufficient to observe the precipitation of the cluster constituting by calcium and phosphate ions. From Figure S2, after ten nanoseconds MD simulation, all 13 segments of collagen models can still keep their stable triple helix structures. Impressively, clusters containing ions of Ca2+, PO43- and OH- with different size accumulate and spread around the surface of each collagen segment. Since the ions in solution is in supersaturated status, the ions aggregated both on the collagen surface and in solution quickly. Notably, the size of clusters on the collagen surface is significantly larger than those in solution. Figure 2A depicts the average size comparison of the clusters in solution and on the collagen surface. Clearly, the precipitation or the size of calcium phosphate in solvent will reach nearly equilibrated after a quite short simulation time, whereas the clusters on the surface can continuously aggregate along the time course. Therefore, we might simply anticipate that the phosphate salt of calcium could finally precipitate, although we have no idea how long it will take. The collagen in solution can facilitate cluster aggregation and even biomineralization for calcium phosphate without question. On the other hand, if we carefully check the components of the cluster on the collagen surface, we can find 8

ACS Paragon Plus Environment

Page 8 of 40

Page 9 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

that the sodium and chloride ions never appear in the clusters. This also means the cluster has the formula of Cax(PO43-)y(OH-)z. Another interesting phenomenon could be observed from Figure S2 is an uneven distribution of Ca-P clusters on the collagen surface. Figure 2B plots the occurrence rate of ions aggregation on the collagen surface, which is defined as the ration of number of ions occurring in clusters and number of total ions in solution for each segment. The possibility of formation of the cluster in the 7th segment collagen model seems to be the largest. On the other hand, those clusters with different sizes are only located at some specific binding sites, which indicate the selectivity of nucleation site. More details will be discussed in below. 3.2 Properties of Clusters As we know, the calcium and phosphate ions will generate clusters in solution when approaching a level of supersaturation. Indeed, such phenomena have been observed for calcium carbonate mineralization in liquid phase.69, theoretical computation69 and experimental observations70,

71,

69

On the basis of

the onset of CaCO3

mineralization shows a clear non-classical nucleation process with distinct sign of amorphous calcium carbonate (ACC), which contradicts classical nucleation theory with unstable clusters. On the other hand, recent experiments also show that the ions can aggregate to generate ACP with the different morphology during the mineralization of hydroxyapatite assisted by collagen proteins.15, 16, 36 In this work, we have systematically investigated the nucleation process in the liquid phase with the presence of collagen protein. In order to understand the properties of the generated clusters along the MD simulations, we have to carefully investigate their size, purity and morphology, since larger and denser clusters should be the key step of nucleation and final crystallization.11 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The size of these clusters on the collagen surface area typically ranges from 100 to 835 atoms within 10 ns MD simulation. Meanwhile, their topologies are dramatically different. Figure 3A plots one example of Ca-P clusters generated on the collagen surface. Interestingly, this cluster even features a distinct porous structure. If we recall that the real bone is a porous-like material constituted by calcium and phosphate ions, current topology might be a spontaneous result regulated by proteins like collagen. On the other hand, we randomly selected three clusters on the collagen surface from the 2nd, 4th and 7th collagen segments to understand the evolution of Ca/P ratio along the simulation time. As shown in Figure 4A, although in the initial nucleation stage, the Ca/P ratio differs dramatically among three clusters, the ratio finally reaches nearly the same value around 1.3 ~ 1.5. Recently, Jahromi et al48 conducted one experiment to investigate the effects of charged amino acids in the crystallization of hydroxyapatite, in which the Ca/P ratio was found to be 1.35 after three days-aged precipitation in solution with the presence of either Glu or Arg. Impressively, the statistically averaged Ca/P ratio for all clusters on the collagen surface is calculated to be about 1.37 in our simulation, which agrees with the experimental data very well.48,

72

At the same time, we also summarized the size

evolution of these three clusters using the number of atoms along the simulation time in Figure 4B. As we can see, they all fall in the range of the critical nucleus size, which should contain about 100 to 1000 atoms in the cluster.11 On the other hand, the Ponser’s cluster73, which has the chemical formula of Ca9(PO4)6, was usually considered to be the fundamental unit of the ACP, and thus critical for the nucleation of apatite in solution. If we carefully check the stoichiometry for all of Ca-P clusters formed on the collagen surface, we can find they all can be expressed as the form of multiple Ponser’s clusters. For example, two clusters shown in Figure 3A and S3A, 10

ACS Paragon Plus Environment

Page 10 of 40

Page 11 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

respectively. Both can be grouped into Ponser's clusters, i.e., [Ponser]12 for Ca108(PO4)72(OH)10 from the seventh segment, while [Ponser]20~[Ponser]21 for Ca182(PO4)127(OH)9 from the second segment if we only consider calcium and phosphate ions. To this point, we could believe all clusters formed in this work should be considered to be the experimentally proposed ACP.45 Another interesting issue is the distribution of water molecule around the clusters.45 We can also notice some water molecules distribute inside the pore. One example was depicted in Figure 3B for all water molecules within the distance of 5 Å around the cluster. Figure S4 gives the mean square deviation (MSD) for three types of ions, OH-, Ca2+ and PO43-. Clearly, the movement speed in solution has the sequence of OH- >> Ca2+ > PO43-. This is quite reasonable considering the mass of the ions. Meanwhile, this can also provide a reasonable explanation that during the precipitation, hydroxyl ions are very difficult to be trapped by collagen surface residues, while phosphate groups are pretty easy to be absorbed to collagen surface charged residues. Moreover, once the calcium or phosphate ions are bound to the collagen surface residues, the interaction is quite stable no matter the size of the Ca-P clusters. Typically, the closest distance between Ca2+ and the carboxylate group of Asp or Glu has the range of [2.5, 3.0] Å, while [3.5, 4.0] Å between phosphate group and positive charged group of Arg or Lys. Such distances clearly indicate that the interaction between Ca-P clusters and collagen protein is pretty strong and stable. Further, this also means that the formed Ca-P clusters will stay on the collagen protein surface, but not move to the solvent. Another property to characterize the cluster is the pair distribution function (PDF) for the Ca-P, Ca-Ca and P-P, which could tell us the micro structure of the generated clusters. For comparison, we also plotted the PDFs for both TCP and HAP in Figure 5. 11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Obviously, the PDF for the clusters is disordered no matter at the short or long range zone. However, for the PDF of Ca-P, two clear peaks occur at 3.16 and 3.78 Å, which are pretty close to crystal properties of HAP. Thus, current simulation can show that the generated clusters have the tendency to further grow up to form an ordered crystal structure. Meanwhile, to further characterize the clusters, we then calculated the average distances between atoms in the clusters, which are summarized in Table 1. According to Table 1, the fundamental geometric parameters of the clusters are quite close to the HAP crystal structure. 3.3 Nucleation Process Analysis One of important issues has to be clearly addressed is that the nucleation process for the calcium phosphate assisted by collagen protein. Traditionally, the nucleation mechanism has two proposals in aqueous phase, namely classical and a two-step nucleation pathway. For the first pathway, the ions could directly aggregate together, and finally precipitate in solution. The second pathway indicates that ions could self-assemble in solution to form many metastable amorphous clusters or spherical nano particles, which further aggregate together to fulfil final precipitation11. Recent experimental investigations15,

16, 36

of biomineralization of calcium phosphate

regulated by collagen protein support that the nucleation process has the characteristics of the non-classical nucleation mechanism or two-step nucleation mechanism evidenced by the observation of ACP. In our simulation, we can also have the similar phenomena with experimental observances. To elucidate the performance of ions of calcium and phosphate in the presence of collagen, we then selected one segment model (the 7th segment) from our 12

ACS Paragon Plus Environment

Page 12 of 40

Page 13 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

simulation for further investigation. In Figure 6, we plotted the evolution along the simulation time course for all calcium and phosphate ions that belonging to the clusters formed on the collagen surface, which is extracted from the seventh collagen segment. Four snapshots at 0, 2, 4, and 10 ns were included, respectively. For clarity, those phosphate ions that finally aggregate together were plotted using the pyramid style, while calcium ions were plotted in ball style. Along the simulation, we can easily identify some calcium and phosphate ions are absorbed by several charged residues. Meanwhile, even in a relatively short time (2 ns), some ions could also aggregate together in solution with a self-assembly way. Quite impressively, nearly 50% clusters could be grouped into Clusters with different size were generated according to Figure 6B. We then statistically summarized the size and stoichiometry form for these clusters in Table S1, in which the Ponser’s cluster was identified according to Ca/P ratio ranging from 1.40 to 1.55. This clearly means the existence of ACP during the collagen regulated nucleation and even final biomineralization. Cluster surface energy is another sign to estimate the stability of generated ΔE

clusters. We calculated the surface energies (),which is defined as 𝛾 = A ,74 for those s clusters bearing characteristics of Ponser's cluster. The model clusters were extracted from the MD trajectory at 2 ns during the nucleation process for the 7th collagen segment. The E can be defined as the energy of a single cluster. The surface area (As) of the cluster was estimated using solvent accessible surface area (SASA) method75. Figure 7 depicts the surface energy. Clearly, the surface free energy shows an opposite tendency with cluster size. For comparison, the first five small clusters, 13

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

which do not belong to Ponser's cluster, were also included. As we can see, with the size increasing, the surface energies reduced accordingly. Such energetic profile also suggests that the formation and growth of the Ca-P clusters should be a thermodynamic favorable process. If we carefully check all thirteen MD trajectories, we could identify at least three ways for the aggregation and growth of ions to produce clusters. First, some calcium and phosphate ions could assemble around the charged residues to form a small cluster, which could further absorb other ions in solution and become bigger gradually. The second pathway is the ionic clusters will firstly form in solution, which could be trapped by other small clusters formed on the collagen surface. The third pathway is the clusters formed on the collagen surface capture other smaller clusters from nearby sites. The simultaneous existence of these three ways to form larger clusters means the nucleation assisted by collagen protein is complicated and also a competitive ionic adsorption process. 3.4 Analysis of Nucleation Sites As shown in Figure S2, we have noticed that uneven distribution of cluster formation or nucleation of calcium phosphate. More importantly, those sites for the nucleation seem to only occur around those charged residues like Asp, Glu, Lys and Arg. First of all, we counted the participation number of amino acids in the formation of Ca-P clusters during the MD simulation. To avoid artificial estimation, we then included all residues from all thirteen segments within 5 Å away from the Ca-P clusters in our analysis. Such results can be seen in Figure S5, in which we calculated possibilities for all types of amino acids and charged residues of all and separate 14

ACS Paragon Plus Environment

Page 14 of 40

Page 15 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

collagen segments. Notably, from Figure 8A, the possibilities evolution trends of the ions occurring close to both all residues and only charged residues are nearly the same. Meanwhile, after relatively short MD simulation, the occurrence possibilities around the Ca-P clusters can reach equilibrium. All these results strongly suggest that those charged residues could be the main factors to facilitate the formation and growth of Ca-P clusters. To further illuminate this result, we then summarized all nucleation sites using amino acid sequence analysis in Figure 9. In particular, the red colour represents negative charged amino acids, while the blue colour denotes positive charged amino acids. Since some of the proline residues in collagen protein are replaced with Hyp, we simply labelled them with an abbreviation "P", but coloured by green to distinguish from normal prolines (colored by black). We can see some basic structural characteristics for the type-I collagen molecule. For example, there are same two or three adjacent amino acids of the same charge (Asp/Glu or Arg/Lys) occurring simultaneously in each helix. At the same time, residues with the same charge (Asp/Glu or Arg/Lys) are often situated nearly at the same position of among different helices. Another structural property for type I collagen molecule is that one charged residue and its opposite charged residue often occur in close proximity.10 Such kind of structure could largely affect the ions absorption and further aggregation without question. As we have discussed above, during the precipitation of calcium phosphate, the protein does show its critical role on the basis of the comparison of cluster size in solution and on the collagen surface. According to Figure 9, we can easily find residues in the square circled regions are mainly charged residues (Arg, Lys, Glu and Asp), around which most of the Ca-P clusters were produced. This suggests that 15

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

residues with both negative and positive charged residues could largely promote the aggregation of ions around themselves. Another interesting phenomenon is nearly all the squared regions contain not only one type of charged residues. To illustrate this, we then plotted the occurrence frequency of Ca-P clusters in Figure 8B upon the group of two residues containing opposite charge, in which the possibility of the nions

accumulation of the clusters can be defined to be ( ntotal ). 𝑛ions denotes the number of Ca2+, PO43- and OH- in the Ca-P cluster around the specific amino acids pairs (Glu-Lys, Glu-Arg, Asp-Lys and Asp-Arg) and 𝑛ions represents the total number of Ca2+, PO43- and OH- in solution. For comparison, we also included two more pairs with same charge like Glu-Asp and Arg-Lys. Impressively, the occurrence rate of Glu-Arg pair seems to be dominant compared to all other cases. On the other hand, we could deduce that the Glu-Arg pair should be the major factor to facilitate the formation of Ca-P clusters. Our results can also show that the nucleation process controlled by type-I collagen molecule has certain directionality. Interestingly, recent studies by Jahromi et al48, 49, 76 strongly show that individual dissolved amino acids could inhibit the precipitation of HA via chelating interactions with either Ca2+ or PO43-. However, the combination of Glu and Arg show less effect in controlling HA nucleation but inhibit the HA crystallization. Clearly, in order to completely understand the functional roles of the biological additives in orientating the HA nucleation, a proper physiological environment might be highly necessary. The detailed interaction patterns between those Ca-P clusters and collagen surface charged residues deserve special investigations. Without loss of generality, we choose two Ca-P clusters formed on the collagen surface of the second and seventh collagen segments for further investigations, which have the stoichiometry of 16

ACS Paragon Plus Environment

Page 16 of 40

Page 17 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Ca108(PO4)72(OH)10 and Ca182(PO4)127(OH)9, respectively. The detailed interactions between two clusters and nearby residues were depicted in Figure 10 and S6. The calcium ions can interact with the carboxylate group of negatively charged amino acids with the non-bond distance around 2.7 Å, whereas the distances between the phosphate ions and positively charged guanidine group of Arg or amine group of Lys are typically less than 4.0 Å. Such short distances could strongly suggest strong electrostatic interactions between Ca-P clusters and collagen surface. In fact, the electrostatic interactions have been proposed to be the major driving force for the interactions between protein and Ca-P based biomaterials.77-79 On the other hand, we can also identify the calcium ions will always interact with only one oxygen atom of the carboxylate group. Similarly, only one of oxygen atoms of phosphate group will be in direct contact with Arg or Lys residues. Moreover, the average coordination number of P around Ca is calculated to be 5.48 for the cluster in Figure 10, while 5.22 for the cluster in Figure S6. Both are pretty close to 5.40 in the HAP crystal structure. This may indicate that the clusters are getting closer to the crystal structure when the type-I collagen serves as the seed for the nucleation of calcium phosphate. 4. Discussion In this work, we have systematically investigated the nucleation process of ions of Ca2+, PO43- and OH- with the help of human type-I collagen molecule. The original human type-I collagen molecule has the repetitive unit of Gly-Xaa-Yaa, in which X and Y replaced with some charged residues and hydroxyproline ocassionally. During the MD simulation, we can specifically observed Ca-P clusters with different size and shape formed on the collagen surface. The distance between the cluster and collagen

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ranges from 2.5 to 4.0 Å, which indicates strong interactions between collagen and clusters. Furhter, during the MD simulation, we can observe some small clusters formed in solution before they migrate to the collagen surface. The stoichiometry of these clusters could have the characteristics of Ponser’s cluster, which is considered to be the prerequisite fundamental unit of crystallization of hydroxyapatite. Therefore, our simulation can offer some solid evidences to support the formation of Ponser’s cluster during the crystallization of hydroxyapatite, although the time scale of our computation cannot reach the limit of crystallization. Another important issue is that the functional role of charged residues in the process of formation and growth of the Ca-P clusters. The nucleation sites with large size Ca-P clusters are mainly around the charged residues, especially the site containing Glu-Arg pair. In other words, the Glu-Arg pair should be considered to be the main factor to promote the growth of Ca-P clusters, and thus be critical in the crystallization of hydroxyapatite. Previously, the mechanisms of nucleation of HAP have been explored to some extents. However, we have to point out that the collagen models used in some previous work are nearly based on collagen mimetic peptide (CMP) model, which contains the repetitive units of Gly-Pro-Hyp. This collagen model can perfectly mimic the tropocollagen style for collagen molecule. However, it shows intrinsically difference from human type-I collagen molecule, e.g., existence of charged residues. Therefore, nucleation simulation of HAP based on this kind of collagen model might not directly provide insights into the functional role of charged residues in the 18

ACS Paragon Plus Environment

Page 18 of 40

Page 19 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

formation of nucleus or clusters of HAP. To compare the performance of ions of calcium and phosphates in the presence of CMP and type-I collagen model, we further carried out the nucleation simulation when CMP system in solution. The coordinate of CMP can be adopted from the protein data bank (PDB entry code 1CGD). In this work, 26 repeat units of Gly-Pro-Hyp were employed as the collagen model system. All of other MD setup protocols are the same as we did for type-I collagen. The size for the system is 290 Å × 80 Å × 90 Å. Total of 10 ns MD simulation was performed. Notably, we cannot observe any clusters formed on the CMP based collagen surface. As shown in Figure 11, although some clusters seem to be generated in solution, the shortest distance from the protein to the cluster is even longer than 8 Å, which means no direction contact between Ca-P clusters and protein. Since the CMP collagen model contains no charged amino acids, it cannot function as the seed during nucleation of hydroxyapatite. In other word, we can confirm the importance of charged residues once again. The collagen model employed in this work was simply divided into 13 segments. Such kind of treatment can save a lot of computational cost and allow all atom simulation to be feasible. However, our treatment of the collagen moleucle might not provide a complete picture for the HAP nucleation. It would be better to conduct this type of simulation for collagen molecule with full length. On the other hand, it has been well known that the hyman type-I collagen molecule is usually in a microfibril form in real case, and further forms a hierarchical strucutre via some organic crosslinks.31, 80 Of course, the size of such system is much over the limit of all-atom 19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

MD simulation could reach. Coarse grained modeling approach might be more approapriate for such large-scale systems. 5. Conclusions Biologically induced or controlled biomineralization represents one of widespread phenomena in biological systems. The mineralization of the salt containing ions of calcium and phosphate has been considered to be particular useful in the bone repairing or the treatment of pathological mineral formation in biological environments. Moreover, the nucleation process of salt of calcium and phosphate is the key step before mineralization. To understand the Ca-P clusters growth process in solution phase in the presence of the human type-I collagen molecule, extensive MD simulation were carried out in this work. In our simulation, we can find that the clusters containing calcium and phosphate ions aggreegate in solution phase and on the collagen surfaces, but with much different size and morphorlogy. Much larger clusters can be found occuring assisted by charged amino acids. Of course, unevenly distributions of Ca-P clusters also indicate that the aggregation of ions containing calcium and phosphate could be induced or controlled by protein surface properties. In particular, the size of the clusters around the pair of Glu-Arg is much larger than other position, which means that the growth of cluster or nucleation should have site specificity or directionality. During the formation of ACP on the collagen surface, we could identify some smaller clusters, which also have the chemical formula of Ponser’s cluster. That strongly supports a non-classical nucleation pathway for the HA nucleation with the collagen-I 20

ACS Paragon Plus Environment

Page 20 of 40

Page 21 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

as the seed. On the other hand, those clusters aggregated on the collagen surface shows dramtically different morphorlogies. Even we could observe the clusters at some collagen segments have porus structure, which also suggest that type-I collagen does have possibility to induce or orientate the growth of calcium phosphate clusters, although direct connection between the morphorlogy and protein cannot be established so far. We still believe that this could provide a chance to pursue a way to regulate or optimize the mineralization via protein engineering technologies.

Acknowledgement This work was supported by National Key Research and Development Program (No. 2016YFB0700801). Some of results described in this paper were obtained on the National Supercomputing Center of Guangzhou and Supercomputing Center of Sichuan University.

Supporting Information The backbone RMSDs of thirteen collagen proteins along the dynamics simulation time. The generated clusters on the collagen surface for all thirteen segments after 10 ns MD. The details of evolution of cluster (Ca108(PO4)72(OH)10) in the 7th collagen segment. Mean square displacement (MSD) values for ions of Ca2+, PO43- and OHalong the simulation time. The probability of the nucleation on the collagen protein surface during the 10 ns simulation. The summary of Ponser’s clusters occurred at 2

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ns in the seventh segment collagen model. This material is available free of charge via the internet at http://pubs.acs.org.

22

ACS Paragon Plus Environment

Page 22 of 40

Page 23 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

References 1. Yao, S. S.; Jin, B. A.; Liu, Z. M.; Shao, C. Y.; Zhao, R. B.; Wang, X. Y.; Tang, R. K., Biomineralization: From Material Tactics to Biological Strategy. Adv. Mater. 2017, 29, 19. 2. Espinosa, H. D.; Rim, J. E.; Barthelat, F.; Buehler, M. J., Merger of Structure and Material in Nacre and Bone - Perspectives on De Novo Biomimetic Materials. Prog. Mater. Sci. 2009, 54, 1059-1100. 3. Barthelat, F., Biomimetics for Next Generation Materials. Philos. Trans. R. Soc. A-Math. Phys. Eng. Sci. 2007, 365, 2907-2919. 4. Vert, M.; Doi, Y.; Hellwich, K. H.; Hess, M.; Hodge, P.; Kubisa, P.; Rinaudo, M.; Schue, F., Terminology for Biorelated Polymers and Applications (Iupac Recommendations 2012). Pure Appl. Chem. 2012, 84, 377-408. 5. Livingston, B. T.; Killian, C. E.; Wilt, F.; Cameron, A.; Landrum, M. J.; Ermolaeva, O.; Sapojnikov, V.; Maglott, D. R.; Buchanan, A. M.; Ettensohn, C. A., A Genorne-Wide Analysis of Biomineralization-Related Proteins in the Sea Urchin Strongylocentrotus Purpuratus. Dev. Biol. 2006, 300, 335-348. 6. Nudelman, F.; Sommerdijk, N., Biomineralization as an Inspiration for Materials Chemistry. Angew. Chem.-Int. Edit. 2012, 51, 6582-6596. 7. Beniash, E., Biominerals-Hierarchical Nanocomposites: The Example of Bone. Wiley Interdiscip. Rev.-Nanomed. Nanobiotechnol. 2011, 3, 47-69. 8. Sosso, G. C.; Chen, J.; Cox, S. J.; Fitzner, M.; Pedevilla, P.; Zen, A.; Michaelides, A., Crystal Nucleation in Liquids: Open Questions and Future Challenges in Molecular Dynamics Simulations. Chem. Rev. 2016, 116, 7078-7116. 9. Silver, F. H.; Landis, W. J., Deposition of Apatite in Mineralizing Vertebrate Extracellular Matrices: A Model of Possible Nucleation Sites on Type I Collagen. Connect. Tissue Res. 2011, 52, 242-254. 10. Landis, W. J.; Jacquet, R., Association of Calcium and Phosphate Ions with Collagen in the Mineralization of Vertebrate Tissues. Calcif. Tissue Int. 2013, 93, 329-337. 11. Erdemir, D.; Lee, A. Y.; Myerson, A. S., Nucleation of Crystals from Solution: Classical and Two-Step Models. Accounts Chem. Res. 2009, 42, 621-629. 12. Demichelis, R.; Raiteri, P.; Gale, J. D.; Quigley, D.; Gebauer, D., Stable Prenucleation Mineral Clusters Are Liquid-Like Ionic Polymers. Nat. Commun. 2011, 2, 8. 13. Karthika, S.; Radhakrishnan, T. K.; Kalaichelvi, P., A Review of Classical and Nonclassical Nucleation Theories. Cryst. Growth Des. 2016, 16, 6663-6681. 14. Nielsen, M. H.; Aloni, S.; De Yoreo, J. J., In Situ Tem Imaging of CaCO3 Nucleation Reveals Coexistence of Direct and Indirect Pathways. Science 2014, 345, 1158-1162. 15. Wang, Y.; Azais, T.; Robin, M.; Vallee, A.; Catania, C.; Legriel, P.; Pehau-Arnaudet, G.; Babonneau, F.; Giraud-Guille, M. M.; Nassif, N., The Predominant Role of Collagen in the Nucleation, Growth, Structure and Orientation of Bone Apatite. Nat. Mater. 2012, 11, 724-733. 16. Nudelman, F.; Pieterse, K.; George, A.; Bomans, P. H. H.; Friedrich, H.; Brylka, L. J.; Hilbers, P. A. J.; de With, G.; Sommerdijk, N., The Role of Collagen in Bone Apatite Formation in the Presence of Hydroxyapatite Nucleation Inhibitors. Nat. Mater. 2010, 9, 1004-1009. 17. Orgel, J.; Irving, T. C.; Miller, A.; Wess, T. J., Microfibrillar Structure of Type I Collagen in Situ. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 9001-9005. 18. Ganss, B.; Kim, R. H.; Sodek, J., Bone Sialoprotein. Crit. Rev. Oral Biol. Med. 1999, 10, 79-98. 19. Du, J. M.; Hou, S.; Zhong, C.; Lai, Z.; Yang, H.; Dai, J. X.; Zhang, D. P.; Wang, H.; Guo, Y. J.; Ding, J. P., Molecular Basis of Recognition of Human Osteopontin by 23C3, a Potential Therapeutic Antibody for Treatment of Rheumatoid Arthritis. J. Mol. Biol. 2008, 382, 835-842.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

20. Fisher, L. W.; Torchia, D. A.; Fohr, B.; Young, M. F.; Fedarko, N. S., Flexible Structures of Sibling Proteins, Bone Sialoprotein, and Osteopontin. Biochem. Biophys. Res. Commun. 2001, 280, 460-465. 21. Zhu, D. H.; Paine, M. L.; Luo, W.; Bringas, P.; Snead, M. L., Altering Biomineralization by Protein Design. J. Biol. Chem. 2006, 281, 21173-21182. 22. Yang, L.; Zhong, C., Advanced Engineering and Biomimetic Materials for Bone Repair and Regeneration. Front. Mater. Sci. 2013, 7, 313-334. 23. Liu, T.; Ding, X. B.; Lai, D. Z.; Chen, Y. W.; Zhang, R. D.; Chen, J. Y.; Feng, X. X.et al.; Zhao, R. B.; Chen, K.; Kong, X. D., Enhancing in Vitro Bioactivity and in Vivo Osteogenesis of Organic-Inorganic Nanofibrous Biocomposites with Novel Bioceramics. J. Mat. Chem. B 2014, 2, 6293-6305. 24. Landis, W. J.; Silver, F. H.; Freeman, J. W., Collagen as a Scaffold for Biomimetic Mineralization of Vertebrate Tissues. J. Mater. Chem. 2006, 16, 1495-1503. 25. Liu, Y.; Luo, D.; Kou, X. X.; Wang, X. D.; Tay, F. R.; Sha, Y. L.; Gan, Y. H.; Zhou, Y. H., Hierarchical Intrafibrillar Nanocarbonated Apatite Assembly Improves the Nanomechanics and Cytocompatibility of Mineralized Collagen. Adv. Funct. Mater. 2013, 23, 1404-1411. 26. Xing, R. T.; Liu, K.; Jiao, T. F.; Zhang, N.; Ma, K.; Zhang, R. Y.; Zou, Q. L.; Ma, G. H.; Yan, X. H., An Injectable Self-Assembling Collagen-Gold Hybrid Hydrogel for Combinatorial Antitumor Photothermal/Photodynamic Therapy. Adv. Mater. 2016, 28, 3669-3676. 27. Fitzpatrick, L. E.; McDevitt, T. C., Cell-Derived Matrices for Tissue Engineering and Regenerative Medicine Applications. Biomater. Sci. 2015, 3, 12-24. 28. D'Elia, N. L.; Mathieu, C.; Hoemann, C. D.; Laiuppa, J. A.; Santillan, G. E.; Messina, P. V., Bone-Repair Properties of Biodegradable Hydroxyapatite Nano-Rod Superstructures. Nanoscale 2015, 7, 18751-18762. 29. Lee, C. H.; Singla, A.; Lee, Y., Biomedical Applications of Collagen. Int. J. Pharm. 2001, 221, 1-22. 30. Hulmes, D. J. S., Building Collagen Molecules, Fibrils, and Suprafibrillar Structures. J. Struct. Biol. 2002, 137, 2-10. 31. Shoulders, M. D.; Raines, R. T., Collagen Structure and Stability. Annu. Rev. Biochem. 2009, 78, 929-958. 32. Mizuno, M.; Fujisawa, R.; Kuboki, Y., Type I Collagen-Induced Osteoblastic Differentiation of Bone-Marrow Cells Mediated by Collagen-α2β1 Integrin Interaction. J. Cell. Physiol. 2000, 184, 207-213. 33. Komsa-Penkova, R.; Koynova, R.; Kostov, G.; Tenchov, B., Discrete Reduction of Type I Collagen Thermal Stability Upon Oxidation. Biophys. Chem. 2000, 83, 185-195. 34. Chang, S. W.; Shefelbine, S. J.; Buehler, M. J., Structural and Mechanical Differences between Collagen Homo- and Heterotrimers: Relevance for the Molecular Origin of Brittle Bone Disease. Biophys. J. 2012, 102, 640-648. 35. Myllyharju, J.; Lamberg, A.; Notbohm, H.; Fietzek, P. P.; Pihlajaniemi, T.; Kivirikko, K. I., Expression of Wild-Type and Modified Pro Alpha Chains of Human Type I Procollagen in Insect Cells Leads to the Formation of Stable [α1(I)]2[α2(I)] Collagen Heterotrimers and [α1(I)]3 Homotrimers but Not [α2(I)]3 Homotrimers. J. Biol. Chem. 1997, 272, 21824-21830. 36. Niu, L. N.; Jee, S. E.; Jiao, K.; Tonggu, L.; Li, M.; Wang, L.; Yang, Y. D.; Bian, J. H.; Breschi, L.; Jang, S. S.;et al., Collagen Intrafibrillar Mineralization as a Result of the Balance between Osmotic Equilibrium and Electroneutrality. Nat. Mater. 2017, 16, 370-378. 37. Kim, D.; Lee, B.; Thomopoulos, S.; Jun, Y. S., In Situ Evaluation of Calcium Phosphate Nucleation Kinetics and Pathways During Intra- and Extrafibrillar Mineralization of Collagen Matrices. Cryst. Growth Des. 2016, 16, 5359-5366. 38. Almora-Barrios, N.; De Leeuw, N. H., Molecular Dynamics Simulation of the Early Stages of Nucleation of Hydroxyapatite at a Collagen Template. Cryst. Growth Des. 2012, 12, 756-763. 24

ACS Paragon Plus Environment

Page 24 of 40

Page 25 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

39. Ma, J., A Molecular Dynamics Study on the Nucleation of Calcium Phosphate Regulated by Collagen. J. Mater. Sci. 2014, 49, 3099-3106. 40. Yang, Y.; Cui, Q. A.; Sahai, N., How Does Bone Sialoprotein Promote the Nucleation of Hydroxyapatite? A Molecular Dynamics Study Using Model Peptides of Different Conformations. Langmuir 2010, 26, 9848-9859. 41. Yang, Y.; Mkhonto, D.; Cui, Q.; Sahai, N., Theoretical Study of Bone Sialoprotein in Bone Biomineralization. Cells Tissues Organs 2011, 194, 182-187. 42. Xu, Z. J.; Yang, Y.; Wang, Z. Q.; Mkhonto, D.; Shang, C.; Liu, Z. P.; Cui, Q.; Sahai, N., Small Molecule-Mediated Control of Hydroxyapatite Growth: Free Energy Calculations Benchmarked to Density Functional Theory. J. Comput. Chem. 2014, 35, 70-81. 43. Zhou, Z.; Qian, D.; Minary-Jolandan, M., Clustering of Hydroxyapatite on a Super-Twisted Collagen Microfibril under Mechanical Tension. J. Mat. Chem. B 2017, 5, 2235-2244. 44. Kawska, A.; Hochrein, O.; Brickmann, A.; Kniep, R.; Zahn, D., The Nucleation Mechanism of Fluorapatite-Collagen Composites: Ion Association and Motif Control by Collagen Proteins. Angew. Chem.-Int. Edit. 2008, 47, 4982-4985. 45. Dorozhkin, S. V., Amorphous Calcium (Ortho)Phosphates. Acta Biomater. 2010, 6, 4457-4475. 46. Tao, J. H.; Pan, H. H.; Zeng, Y. W.; Xu, X. R.; Tang, R. K., Roles of Amorphous Calcium Phosphate and Biological Additives in the Assembly of Hydroxyapatite Nanoparticles. J. Phys. Chem. B 2007, 111, 13410-13418. 47. Chu, X. B.; Jiang, W. G.; Zhang, Z. S.; Yan, Y.; Pan, H. H.; Xu, X. R.; Tang, R. K., Unique Roles of Acidic Amino Acids in Phase Transformation of Calcium Phosphates. J. Phys. Chem. B 2011, 115, 1151-1157. 48. Jahromi, M. T.; Yao, G.; Cerruti, M., The Importance of Amino Acid Interactions in the Crystallization of Hydroxyapatite. J. R. Soc. Interface 2013, 10, 14. 49. Jahromi, M. T.; Cerruti, M., Amino Acid/Ion Aggregate Formation and Their Role in Hydroxyapatite Precipitation. Cryst. Growth Des. 2015, 15, 1096-1104. 50. Chen, Y.; Gu, W. J.; Pan, H. H.; Jiang, S. Q.; Tang, R. K., Stabilizing Amorphous Calcium Phosphate Phase by Citrate Adsorption. Crystengcomm 2014, 16, 1864-1867. 51. Shao, C. Y.; Zhao, R. B.; Jiang, S. Q.; Yao, S. S.; Wu, Z. F.; Jin, B.; Yang, Y. L.; Pan, H. H.; Tang, R. K., Citrate Improves Collagen Mineralization Via Interface Wetting: A Physicochemical Understanding of Biomineralization Control. Adv. Mater. 2018, 30, 7. 52. Newcomb, C. J.; Bitton, R.; Velichko, Y. S.; Snead, M. L.; Stupp, S. I., The Role of Nanoscale Architecture in Supramolecular Templating of Biomimetic Hydroxyapatite Mineralization. Small 2012, 8, 2195-2202. 53. Dey, A.; Bomans, P. H. H.; Muller, F. A.; Will, J.; Frederik, P. M.; de With, G.; Sommerdijk, N., The Role of Prenucleation Clusters in Surface-Induced Calcium Phosphate Crystallization. Nat. Mater. 2010, 9, 1010-1014. 54. Rainey, J. K.; Goh, M. C., An Interactive Triple-Helical Collagen Builder. Bioinformatics 2004, 20, 2458-2459. 55. Fiser, A.; Sali, A., Modeller: Generation and Refinement of Homology-Based Protein Structure Models. In Macromolecular Crystallography, Pt D, Carter, C. W.; Sweet, R. M., Eds. Elsevier Academic Press Inc: San Diego, 2003; Vol. 374, pp 461-491. 56. Hess, B.; Kutzner, C.; van der Spoel, D.; Lindahl, E., Gromacs 4: Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. J. Chem. Theory Comput. 2008, 4, 435-447. 57. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L., Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926-935. 58. Darden, T.; York, D.; Pedersen, L., Particle Mesh Ewald: An N⋅log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98, 10089-10092. 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

59. Brooks, B. R.; Bruccoleri, R. E.; Olafson, B. D.; States, D. J.; Swaminathan, S.; Karplus, M., Charmm: A Program for Macromolecular Energy, Minimization,and Dynamics Calculations. J. Comput. Chem. 1983, 4, 187-217. 60. Park, S.; Radmer, R. J.; Klein, T. E.; Pande, V. S., A New Set of Molecular Mechanics Parameters for Hydroxyproline and Its Use in Molecular Dynamics Simulations of Collagen-Like Peptides. J. Comput. Chem. 2005, 26, 1612-1616. 61. Committee on Minerals and Toxic Substances in Diets and Water for Animals, National Research Council. National Academies Press: Washington, DC, 1980. 62. Lukasheva, N. V.; Tolmachev, D. A., Cellulose Nanofibrils and Mechanism of Their Mineralization in Biomimetic Synthesis of Hydroxyapatite/Native Bacterial Cellulose Nanocomposites: Molecular Dynamics Simulations. Langmuir 2016, 32, 125-134. 63. Plimpton, S., Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1-19. 64. Humphrey, W.; Dalke, A.; Schulten, K., Vmd: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33-38. 65. Hockney, R. W.; Eastwood, J. W., Computer Simulation Using Particles. CRC Press New York, 1988. 66. Ryckaert, J. P.; Ciccotti, G.; Berendsen, H. J. C., Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of N-Alkanes. J. Comput. Phys. 1977, 23, 327-341. 67. Andersen, H. C., Rattle: A “Velocity” Version of the Shake Algorithm for Molecular Dynamics Calculations. J. Comput. Phys. 1983, 52, 24-34. 68. Lin, T. J.; Heinz, H., Accurate Force Field Parameters and Ph Resolved Surface Models for Hydroxyapatite to Understand Structure, Mechanics, Hydration, and Biological Interfaces. J. Phys. Chem. C 2016, 120, 4975-4992. 69. Wallace, A. F.; Hedges, L. O.; Fernandez-Martinez, A.; Raiteri, P.; Gale, J. D.; Waychunas, G. A.; Whitelam, S.; Banfield, J. F.; De Yoreo, J. J., Microscopic Evidence for Liquid-Liquid Separation in Supersaturated Caco3 Solutions. Science 2013, 341, 885-889. 70. Pouget, E. M.; Bomans, P. H. H.; Goos, J.; Frederik, P. M.; de With, G.; Sommerdijk, N., The Initial Stages of Template-Controlled CaCO3 Formation Revealed by Cryo-Tem. Science 2009, 323, 1455-1458. 71. Rieger, J.; Frechen, T.; Cox, G.; Heckmann, W.; Schmidt, C.; Thieme, J., Precursor Structures in the Crystallization/Precipitation Processes of Caco3 and Control of Particle Formation by Polyelectrolytes. Faraday Discuss. 2007, 136, 265-277. 72. G.H.Nancollas; M.S.Mohan, The Growth of Hydroxyapatite Crystals. Arch. Oral Biol. 1970, 15, 731-745. 73. Yin, X. L.; Stott, M. J., Biological Calcium Phosphates and Posner's Cluster. J. Chem. Phys. 2003, 118, 3717-3723. 74. Woodruff, D. P., The Chemical Physics of Solid Surfaces. Elsevier: 2002; Vol. 10. 75. Eisenhaber, F.; Lijnzaad, P.; Argos, P.; Sander, C.; Scharf, M., The Double Cubic Lattice Method: Efficient Approaches to Numerical Integration of Surface Area and Volume and to Dot Surface Contouring of Molecular Assemblies. J. Comput. Chem. 1995, 16, 273-284. 76. Tavafoghi, M.; Cerruti, M., The Role of Amino Acids in Hydroxyapatite Mineralization. J. R. Soc. Interface 2016, 13, 12. 77. Dong, X. L.; Zhou, H. L.; Wu, T.; Wang, Q., Behavior Regulation of Adsorbed Proteins Via Hydroxyapatite Surface Texture Control. J. Phys. Chem. B 2008, 112, 4751-4759. 78. Zhou, H.; Wu, T.; Dong, X.; Wang, Q.; Shen, J. W., Adsorption Mechanism of Bmp-7 on Hydroxyapatite (001) Surfaces. Biochem. Biophys. Res. Commun. 2007, 361, 91-96. 79. dos Santos, E. A.; Farina, M.; Soares, G. A.; Anselme, K., Surface Energy of Hydroxyapatite and Beta-Tricalcium Phosphate Ceramics Driving Serum Protein Adsorption and Osteoblast Adhesion. J. Mater. Sci.-Mater. Med. 2008, 19, 2307-2316. 80. Kwansa, A. L.; De Vita, R.; Freeman, J. W., Tensile Mechanical Properties of Collagen Type I and Its Enzymatic Crosslinks. Biophys. Chem. 2016, 214, 1-10. 26

ACS Paragon Plus Environment

Page 26 of 40

Page 27 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

27

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. An example initial state of collagen model in aqueous phase in presence of

calcium, phosphate and hydroxyl ions.

28

ACS Paragon Plus Environment

Page 28 of 40

Page 29 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2. (A)The size comparison of two largest clusters in solution and on the collagen surface in the 7th collagen segment model, respectively. For convenience, we only count the number of calcium and phosphate ions occurring in the clusters. (B) The proportion of the occurrences of Ca-P clusters on the collagen surface.

29

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. (A) One typical cluster (Ca108(PO4)72(OH)10) formed on the collagen surface extracted from the last frame of MD trajectories. (B) The water distribution around the cluster with distance of 5 Å.

30

ACS Paragon Plus Environment

Page 30 of 40

Page 31 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4. (A) Ca/P ratio for selected clusters as a function of the simulation time. (B) Numbers of atoms for these selected clusters.

31

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 5. Pair distribution functions for the Ca-P, Ca-Ca and P-P of (A) HAP (B) -TCP (C) Cluster in the seventh collagen segment model.

32

ACS Paragon Plus Environment

Page 32 of 40

Page 33 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 6. The evolution along the dynamics time course (0, 2, 4 and 10 ns) for those calcium and phosphate ions that belonging to the clusters formed on the collagen surface, which is extracted from the seventh collagen segment. For clarity, those phosphate ions that finally aggregate together were plotted using the pyramid style, while calcium ions were plotted using ball style.

33

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. The surface energy of clusters in the 7th collagen segment after 2ns MD simulation.

34

ACS Paragon Plus Environment

Page 34 of 40

Page 35 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 8. (A) The probability of the nucleation on the collagen protein surface during the 10 ns simulation. (B) The comparison of occurrence ratio of clusters upon the charge residues pair during the 10 ns simulation.

35

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 9. Sequence summarization for the nucleation sites. The N-terminal was capped by GFPGPK, whereas the C-terminal part was capped by GEQGPA for three chains in all thirteen collagen segment models. All negatively charged residues (Asp and Glu) were colored in red, while positively charged residues were in blue color. The red square denotes the number of calcium, phosphate and hydroxyl ions in the cluster is over 300 , pink is between 200 and 300, purple square contains 100-200 ions and less than 100 ions for black square region.

36

ACS Paragon Plus Environment

Page 36 of 40

Page 37 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 10.

The details of nucleation site between the collagen protein and and one of clusters extracted from the second segment of collagen models, in which the cluster has the stoichiometry of Ca108(PO4)72(OH)10.

37

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 11. Snapshot for the CMP model after 10 ns nucleation MD simulation, (A) side view (B) zoom-in view from C-terminus.

38

ACS Paragon Plus Environment

Page 38 of 40

Page 39 of 40 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table1. Average particle-particle distances for clusters on the surface and for HAP crystal. Distances (Å)

Crystal Structure

Clusters

Ca-Ca

3.94

4.48 ± 0.39

Ca-P

3.08~3.67

3.68 ± 0.45

Ca-O(H)

2.38

3.19 ± 0.40

Ca-O(P)

2.34~2.81

3.08 ± 0.35

39

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TOC Graphic

40

ACS Paragon Plus Environment

Page 40 of 40