Effect of Hydroxyapatite Surface on BMP-2 Biological Properties by

Subscriber access provided by UNIV OF LOUISIANA

B: Fluid Interfaces, Colloids, Polymers, Soft Matter, Surfactants, and Glassy Materials

The Effect of Hydroxyapatite Surface on BMP-2 Biological Properties by Docking and Molecular Simulation Approach Haojie Gu, Zhiyu Xue, Menghao Wang, Mingli Yang, Kefeng Wang, and Dingguo Xu J. Phys. Chem. B, Just Accepted Manuscript • Publication Date (Web): 26 Mar 2019 Downloaded from http://pubs.acs.org on March 27, 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 26 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

The Effect of Hydroxyapatite Surface on BMP-2 Biological Properties by Docking and Molecular Simulation Approach Haojie Gu†, Zhiyu Xue†, Menghao Wang#, Mingli Yang‡,£, Kefeng Wang§, ‡,* and Dingguo Xu†,‡* †College

of Chemistry, MOE Key Laboratory of Green Chemistry and Technology, and ‡Genome Research Center of Biomaterial, Sichuan University, Chengdu, Sichuan 610064, P. R. China.

§National

Engineering Research Center for Biomaterials, Sichuan University, Chengdu, Sichuan 610065, P. R. China.

#Key

Laboratory of Advanced Technologies of Materials, Ministry of Education, School of Materials Science and Engineering, Southwest Jiaotong University, Chengdu, Sichuan 610031, P. R. China. £Institute

of Atomic and Molecular Physics, MOE Key Laboratory of High Energy Density Physics and Technology, Sichuan University, Chengdu, Sichuan 610065, P. R. China.

*Corresponding authors: KF Wang, [email protected] DG Xu,

[email protected]

1

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

Abstract: The interactions between osteogenic proteins and biomaterials surface is crucial to the application of biomaterials, in which the conformational or orientational change of the adsorbed protein on the solid surfaces is one of the most important interactions other than the protein adsorption. Although some progress has been made in the mechanism of proteins adsorption on the surface of Hydroxyapatite (HAP) in recent years, there is still insufficient atomistic/molecular information about the conformation and orientation of protein upon adsorbing on solid surfaces. In this study, different orientations and conformations of Bone Morphological Protein- 2 (BMP-2) adsorbed on the surface of HAP were calculated through the protein-solid surface docking approach; the relationship between optimal adsorption and biological activity of BMP-2 was investigated by applying a in combination with Molecular Dynamic simulation (MD) and Steered Molecular Dynamic simulation (SMD). Two favorite adsorption conformers were screened out according to the docking results on the basis of orientations of BMP-2 with different epitopes. Subsequent MD and SMD results showed that the knuckle epitope of BMP-2 was easier to adsorb on the surface of HAP (100) than the wrist epitope companying with certain conformational changes. Such absorption mode led to the wrist epitope of BMP-2 being exposed to the environment. This can be specifically identified/interacted with type I receptors on the stem cell membrane, which further induces the differentiation of stem cells into osteoblasts. Current simulation provided a theoretical high-throughput screening method for the protein-biomaterials adsorption states. It can be extended to more researches about different proteins adsorption on the surface of different materials. The simulation results provided more information at the molecular and atomic level to further interpret the mechanism of osteoinductivity from the perspective of growth factor adsorption. Meanwhile, we believe it should be a meaningful attempt to screen biomaterials key factors by high-throughput method, which might become a promising way to develop or optimize new biomaterials.

2


Page 2 of 26

Page 3 of 26 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


1. Introduction The interaction of protein and biomaterials surfaces is the key to design biomaterials and bioinstruments. The behaviors of proteins or cells are just based on the complex phenomenon of biomolecular interaction with the surfaces. When biomaterials are implanted in to living body, proteins play an important role in the determination of biological properties. Proteins are not only the intermediary to study the interactions between cells and material surfaces, but also effectively enhance the proliferation and differentiation of cells by adsorption of target protein molecules onto the material surfaces. From the perspective of protein, biomaterials achieve its biological properties mainly through two aspects. One is the adsorption of certain efficient functional proteins.1-2 Another is the conformational change of the adsorbed target protein due to the specific biomaterials’ properties.3-5 The adsorption of protein not only ensures that the proteins have suitable orientations, but also promises that the proteins maintain their active conformations. It directly affects the sensitivity and effectiveness of the interaction between proteins and cells. Many works on protein adsorption behaviors have been carried out to better understand the mechanism of protein adsorptions and the relationship between biomaterials and its biological properties. It has been attracted lots of attentions in the process of bone tissue repairing, especially, in the early stage of bone mineralization. Born Morphogenetic Protein-2 (BMP-2), which belongs to the protein of the transforming growth factors- (TGF-), is one of the most important bone regeneration proteins and usually applied as the target protein of research. It can induce the formation and even take part in the growth, development and reconstruction of the bone and cartilage.6-9 BMP-2 binding to its SER-THR receptors triggers a signal cascade of osteoblast differentiation that induces the formation of bone and participates in the growth, development and reconstruction.10-15 A large number of experiments have shown that there is high-affinity interactions between BMP-2 and active type I receptor, and the affinity binding to type II receptor is relative lower.13, 15-16 Hydroxyapatite (HAP), Ca10(PO4)6(OH)2, one of the calcium phosphates, is the major mineral component of the human bone.17-19 It has been widely applied in bone regeneration for its good biocompatibility, bioactivity, even osteoinductivity which is proved by Zhang et al.20-23 Due to the advantages of HAP, many researchers have 3



Page 4 of 26

explored the interactions between its surfaces and proteins, especially with BMP-2.9, 24-29

For instance, Long et al. confirmed the helical structure of the salivary protein

terminal adsorbed on HAP surface by ssNMR.30 Dolatshahi-Pirouz et al studied the adsorption behaviors of fibronectin on surfaces of HAP.31 Zhang et al explored the adsorption of chitosan molecules on different crystal planes of HAP by MD simulations.32 These studies revealed that the conformation, binding sites and activity of protein could be affected by the surfaces. In addition, some studies showed that the osteoinductivity of HAP might be caused by the specific adsorption of BMP-2 and changes of the conformation of BMP-2.9, 28-29 A suitable conformation makes BMP-2 to interact naturally with its receptors on the cell membrane, thereby triggering a signal to induce cell differentiation. For now, there are lots of researches about the adsorption of BMP-2 but the study of conformational changes of BMP-2 caused by the surfaces of HAP is still rare.33-37 Importantly, monomer BMP-2 molecule, which is inactive in reality, was usually adopted in previous simulations.9, 25, 38-39 Thus, in order to better understand the bioactivity and osteoinduction mechanism of HAP through the view of BMP-2, further investigations are then required to reveal the effect of material surface properties on the conformation of BMP-2 homodimer. Although a few of experiments have been conducted to study the conformation changes of BMP-2 on different surfaces,32,

40-41

the information at the

atomistic/molecular level is still inadequate. Computer simulation has become an important means of materials genome research,42 which gives new protocols for studying the interactions between proteins and material surfaces from micro perspectives. It has been adopted to investigate the adsorption of proteins on material surfaces, for example, Wu and coworkers studied the micro-adsorption mechanism of the 10th type III module of fibronectin (FN-III10), BMP-2 and BMP-7 on the surface of HAP by using molecular dynamics (MD) and steered molecule dynamics (SMD).9, 43-44

Utesch et al simulated the initial adsorption process of BMP-2 on hydrophobic

graphite and hydrophilic titanium dioxide.26 Zhou et al studied the adsorption orientation of FN-III10 and FN-III7-10 adsorbed on HAP (001) surfaces in cavities of calcium ions with Parallel tempering Monte Carlo (PTMC) and MD simulation methods.45 Similar with above studies, most of researchers usually focused on the adsorption process and behaviors, but ignored the conformational changes of the adsorbed protein. Fortunately, some attempts have been reported to simulate the 4


Page 5 of 26 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


protein conformation changes upon adsorption. Liao et al investigated the adsorption behaviors and conformation changes of the basic fibroblast growth factor (bFGF) adsorbed on HAP (001) surface by the method of combination of replica-exchange molecular dynamics (REMD) and conventional MD (CMD).35 Cui et al established the most stable adsorption structure of peptide on HAP by RosettaSurface method combining with MD and predicted the relationship between structures and biological activity of HAP-binding peptide.46 Daggett et al employed MD with raising temperature to investigate the protein folding-unfolding pathway.47 Although above attempts expanded our knowledge about the protein conformation changes at certain extent, the computing methods still needed to be improved. MD with raising temperature has a problem of insufficient configuration of sampling. REMD method can achieve enhanced sampling for energetically favorable adsorption states, but it costs huge amount of calculation. On the other hand, RosettaSurface method can obtain representative samples with the lowest global energy and reasonable computational cost. Thus, it would be one of the practical methods for protein conformation computation. In this study, molecular docking and MD simulations were adopted to explore the conformational changes of BMP-2 homodimer on the HAP (100) surface which is one of the most important planes of HAP crystalline,40-41,

48-49

and the adsorption

behaviors of BMP-2 with energy favorable conformations. It would be expected to provide detailed information about the effects of HAP surface structure on the BMP-2 conformation changes and the corresponding adsorption characteristics. The results will be helpful to further reveal the bioactivity nature and the osteoinduction mechanism of Ca-P bioceramics at the atomic/molecular level. As the prospective study of biomaterials genome research, it also provides a new method for designing and screening of tissue inducible biomaterials such as calcium phosphates. 2. Models and methods 2.1 Models BMP-2. The configuration of BMP-2 homodimer was obtained from Protein Data Bank (PDB id 1REW), which contains two active regions, wrist epitope and knuckle epitope.10 The receptors and all crystallographic water molecules were removed with resulted BMP-2 model containing two homo chains, as shown in Figure 1A. 5



HAP (100) Surface. The initial molecule of HAP with space group P63/m was taken from the American Mineralogist Crystal Structure Database. The unit cell parameters are a = b = 9.424 Å and c =6.879 Å.50 The HAP (100) surface51 was prepared as the target model, as shown in Figure 1B. The material size was 18 Å x 130 Å x 130 Å in the x, y and z directions. Total number of atoms was estimated over 30000.

Figure 1. (A) The dimer crystal structure of BMP-2 in Protein Data Bank (PDB id 1REW). (B) HAP (100) surface structure: Side view and top view.

2.2 Methods Docking Method: RosettaSurface was adopted to implement the docking screening. It is a computational tool that simultaneously optimizes protein folding, side chain conformation, and rigid body orientation on solid surfaces. The method is based on the algorithm of all-atom Monte Carlo plus-minimization search and adds a weighted linear combination of van der Waals, electrostatic, hydrogen-bond and implicit solvent interaction on the basis of the all-atom score function of protein folding. It converges to a set of structures with low energy and favorable entropy. At present, this method has been applied in the interaction of peptides and HAP crystal face.52-55 Force Field: The BMH parameters derived from Hauptmann et al56 was usually applied in studying interactions of proteins and HAP surfaces.9, 38, 43-44, 46, 49 However, the BMH parameters have been considered to be particularly appropriate to describe the bulk properties of HAP.43-46 Recently, new parameters, namely INTERFACE force field (or IFF),57 were specifically designed to describe the interface properties of HAP, and thus have been applied to tackle the interactions between HAP surface and biological molecules.57-58 Previously, we even applied this force field to construct an ordered-to-disordered structure for HAP upon stimulated annealing simulation.59 Herein, the IFF parameters,57 which have been integrated in to CHARMM2760-61 6


Page 6 of 26

Page 7 of 26 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


force field, were employed in understanding interactions between BMP-2 and HAP (100). In addition, HAP crystals were semi-fixed rather than completely fixed through that no force was applied to the surface of about 12 Å. The intramolecular and intermolecular interactions were described by CHARMM27 all atom force field for BMP-260, 62 and TIP3P water model.63 2.3 Simulation Details The corresponding protein docking and conformational screening were carried out using the RosettaSurface52, 64-65 module implemented in the Rosetta 3.4 software65 with an implicit solvation model. Dictionary of Secondary Structure of Protein (DSSP)66 software was adopted to do the secondary structural analysis. All obtained complex models are further subjected to Cluster Analysis67. In order to have sufficient samples for screening, total of 1000 adsorption models of BMP-2 with different orientations on HAP (100) surface were obtained in this study. As we have mentioned above, the BMP-2 has two epitope sites, knuckle and wrist, with different physiological functions. Therefore, two different models, which are related with close interactions between two epitope sites and material surface, were screened out from these 1000 conformers. Subsequently, classical MD simulations were applied to get the target models fully relaxed. The secondary structure of BMP-2, adsorption characteristics/stability were then analyzed on the basis of the resulted trajectories. All MD simulations were performed by LAMMPS suite of program.68 The initial structures of HAP and BMP-2 were fully solvated in a pre-equilibrated TIP3P water box (100Å x 130Å x 130Å) using the GROMACS 4.6.7 genbox tool.69 Chloride and sodium ions with a concentration of 0.1 M were added to neutralize the whole system. Firstly, the energy minimization was performed by the method of the conjugated gradient.70 Subsequently, as long as 100 ns MD simulations were executed at 310 K with Nosè-Hoover temperature coupling under the NVT ensemble.71-72 The Verlet leapfrog algorithm was applied to deal with the equations of motions at 2 fs time step. The particle-particle-particle-mesh (PPPM) algorithm was adopted to treat the long-rang electrostatic interaction.73 Molecular configurations were visualized by VMD 1.9.1.74 In order to evaluate the binding affinities of BMP-2 molecule by HAP for different conformers, the potentials of mean force (PMFs) for desorption processes were further calculated by SMD.75-77 The BMP-2 molecules were slowly moved to 7



retain the quasi-static states.78-80 Typically, the last frame from the MD trajectory was extracted as the initial structure for SMD simulations. The constant velocity mode of pulling was adopted in current simulation. In the process of SMD, HAP atoms were fixed and a uniform pulling force was acted to atoms on the backbone of BMP-2 to make the BMP-2 move at a constant speed. Water molecules moved freely. The system was calculated under periodic boundary conditions. The spring constant was set at 20 kcal mol-1Å-2, and velocity (ν) was assigned to 10-5Å fs-1. 3. Results and Discussions 3.1 Docking Structure Evaluation RosettaSurface provides a suitable method to investigate and screen the most possible protein adsorption state on the material surface. In this work, total of 1000 docking structures were firstly generated by RosettaSurface module. Subsequently, these structures were divided into three groups according to the different adsorbed active regions of the BMP-2 as shown in Figure 2 and Supporting information Table S1. They were named as Model Ws, Model Ks and Others, respectively. Model Ws

Figure 2. The 2D pie chart of the proportion of various adsorption models in 1000 docking structures.

represented models with the adsorbed residues mainly located in wrist epitope of BMP-2. Similarly, the models with adsorbed residues mainly located in the knuckle epitope were called Model Ks. States not belonging to these two models were grouped into Others. As we have described above, the models related with absorption states of knuckle and wrist epitopes on the HAP surface deserve particular attentions. Therefore, two optimal models (shown in Figure 3A&B) were then screened out from these docked structures with the lowest 8


Page 8 of 26

Page 9 of 26 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. The initial models of MD simulation by Rosetta docking approach (A) Model W and (B) Model K. All residues belonging to wrist epitope are colored in yellow, while red color for residues of knuckle epitope are also labeled.

interaction energies and the greatest number of contacting residues belonging to either wrist or knuckle region. The screening results data were shown in Table S1 of Supporting Information. For the screened Model K, BMP-2 mainly contained 39.9% for -strand and -sheet, and 17.9% for the helix. There are 11 residues occurring in the adsorbed amino acids including Ala34, Pro35, Pro36, His39, Phe41, His54, Ile87, Ser88, Leu90, Asn95 and Glu109. These residues mostly belong to the knuckle epitope, namely the region of binding for the type II receptors. Although the secondary structure of BMP-2 in Model W was similar to Model K, the binding residues were found to be Asp25, Val26, Gly27, Trp31, Phe49, Pro50, His54, Val70, Ser72 and Asn95. These are mainly seated in wrist epitope, which can recognize type I receptors. Thus, the conformation and orientation of BMP-2 in the two optimal models were basically in line with the binding mode of the interaction with cells, in which they exerted different biological properties. Of course, only docking simulation cannot tell which one, Model W or Model K, is the closest to the practical adsorption state. Due to that the explicit water molecules and ions cannot calculated in RosettaSurface, it might result in simulations that were relative different from the realistic physiological environment. In order to reproduce the interactions between atoms, MD simulation is necessary to verify the stability and rationality of two docking optimal models.

9



3.2 MD Simulation of BMP2 Adsorption 100 ns MD simulations were performed in the explicit water model to observe the adsorption dynamic behaviors and the conformational changes of the two optimal docking structures. Figure 4 showed the changes of the model structures along the

Figure 4. (A) The structure change and active residues for the Model W during the MD simulation of 100 ns. (B) The structure change and active residues for the Model K during the MD simulation of 100 ns. Color designation is same as in Figure 3.

simulation time, and the relative position of BMP-2 active sites was also included. Meanwhile, it was interesting to note that the surface of HAP was partially dissolved by simulating with IFF force field parameters. PO43- and Ca2+ ions on the surface area became gradually disordered as shown in Figure 5A&B. It was found that the material surfaces of both models showed disordered layers with different degrees. The thickness of disordered layer was estimated to be about 15.2 Å for Model W, while about 17.6 Å for Model K, respectively. This result could be consistent with the

10


Page 10 of 26

Page 11 of 26 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. The disordered HAP (100) surfaces of IFF force field displayed using VMD quicksurf module. (A) Model W and (B) Model K. All Ca2+ are colored by cyan, the PO43- ion in red and OH- ions are labeled in pink.

surface characteristics of nanomaterials studied by Jäger81 and Bertinetti82 through UHRTEM. They revealed that solid materials had an amorphous and disordered layer of 1-2 nm thickness on the surface after high temperature sintering treatment. At the same time, this phenomenon was also revealed in the simulation by us.59 To this point, current results might ensure us to apply the IFF force field for further simulation of proteins interacted with the HAP surface. In addition, the Root Mean Square Deviation (RMSD) of BMP-2 backbone atoms revealed that all systems had reached equilibrium finally as shown in Figure 6. For Model K, its RMSD value was about 1 Å with pretty small fluctuation, which indicated that the conformation of BMP-2 was relatively stable. The RMSD of Model

Figure 6. The RMSD of the BMP-2 backbone atoms for Model W and Model K along the MD simulation time.

11



W was about 2 Å with significant fluctuation, but still in the equilibrium state. Much larger fluctuation than that of Model K manifests that BMP-2 skeleton movement of Model W should be more intense. Compared with the secondary structure of the initial docked BMP-2, the ratio of each secondary structure of the two models showed distinct fluctuations throughout the MD simulations, as shown in Figure 7. The β structures accounted for the majority

Figure 7. The percentage of second structure of protein during the MD simulation. (A) Model W and (B) Model K.

in the two models. In particular, the percentage of β structures in Model K showed a significant increasing trend (25.7% ~ 34.0%). However, for the structures of helix and turn, the change of its proportion showed a downward trend (13.1% to 9.7%, 10.2% to 3.40%). From the point of view of energy, the instable helix and turn structures was easy to change to a more stable β structures in the simulation process.83-84 It is worth noting that the knuckle epitope is mainly β-sheet structures, while the wrist epitope is mainly helixes including place-helix and pre-helix loop. The proportion of Model K changes in the whole secondary structure is about 4%. The proportion of secondary structure can affect the intermolecular recognition even the stiffness of the biomolecules.83 But at present, there are no more reports and data to interpret the effects of changes in the secondary structure proportion of BMP-2 on its biological activity. Further studies are needed in the future works. At the same time, residues around the adsorption sites of BMP-2 were also changing in accord with evolution time, as shown in Figure 8. Initially, the residues of BMP-2 responsible for adsorption in Model K are located between the β-7 (Ser88, Leu90, Leu92) and β-8 (Lys97, Leu100, Lys101, Asp105, Glu109) regions, which

12


Page 12 of 26

Page 13 of 26 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


belong to knuckle epitopes of BMP-2. On the other hand, the adsorption residues in Model W are mainly located in wrist epitope. With the simulation proceeding, the

Figure 8. The interacting residues of Model W on the HAP(100) surface during the MD simulation. (A) 500 ps, (B) 1 ns and (C) 10 ns. All name of the residues on the wrist epitope are colored by black, the name of residues on the knuckle epitope are colored by blue.

original residues (including Asp25, Gly27, Asp30, Trp31, Pro50, Val70,) responsible for adsorption moved away from the surface gradually. Meanwhile, those residues (including Glu94, Asn95, Glu96) in knuckle epitope area were getting close to the surface and contributed to new recognition modes. At 500ps, residues such as Asp25, Gly27, Asp30, Trp31, Pro50 and Val70 were mainly adsorbed on the crystal surface in the Model W, belonging to the wrist epitope. Then, the number of residues in the wrist epitope gradually decreased at 1 ns. At last, the interacting residues with the surface were located mainly in the knuckle epitopes, which were mainly Glu94, Asn94 and Glu96 at 10 ns, as shown in Figure 8. In fact, according to Figure 9A, we could find that the number of residues around the adsorption sites in wrist epitope

Figure 9. The number of interacting residues on the HAP(100) surface during the MD simulation. (A) Model W and (B) Model K. 13



Page 14 of 26

decreased continuously before stable adsorption. Overall, the number of the adsorption sites in knuckle epitope was more than that in wrist epitope in most of simulation time. Interestingly, for Model K, we could see from Figure 9B that the recognition of BMP-2 by HAP surface should predominantly attribute to residues belonging to knuckle epitope. To this point, we might simply postulate that the HAP (100) surface might prefer binding knuckle epitope of BMP-2. The residues responsible for recognition in Model K were mainly constituted by Asp, Lys and Glu, which could have strong electrostatic interactions with ions of PO43- and Ca2+ on HAP (100) surface. On the other hand, the interactions between protein and material surface in Model K mainly occured between Ca 2+ ion and residues of Asp and Glu. Indeed, the distances between Ca2+ and the carboxylate group of glutamic acids were calculated to be around 2.6 Å (Table 1) for two models, which Table 1. The distances between surface atoms and adsorbed residues during the MD simulation.

Distance(Å)

Model W

Model K

Glu94-Ca2+

_

2.63±0.11

Glu109-Ca2+

_

3.03±0.86

Asp53-Ca2+

_

3.04±0.81

Lys97-PO43-

_

3.11±0.71

Asp93-Ca2+

_

4.89±0.51

Glu96-Ca2+

2.65±0.14

_

Lys97-PO43-

3.20±0.89

_

clearly indicates strong adsorption states in our simulation. However, according to Table 1 and Figure 10, pretty weak -NH2-PO43- interactions can be found when the adsorption of BMP-2 on HAP (100) surface occurs. Current results are then not inconsistent with previous studies on the interactions between proteins and HAP surfaces.9, 43-44 In addition, due to the different surface properties and protein types, 14


Page 15 of 26 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 optimal adsorption orientation and preferred adsorption sites might be different. For example, Dong et al have found that it was impossible for Glu94 of

Figure 10. Snapshot representations for the complex structures of BMP-2 adsorbed on HAP (100) surface at 100 ns extracted from MD trajectories. Possible interactions between ions on HAP surface and charged residues from BMP-2 are also given for clarification. The main residues interacting with the surface are plotted in ball-and-sticks style. (A) Model W and (B) Model K. Distance units are Å.

15



BMP-2 to become a major adsorption residue due to its orientation. However, on the basis of our simulations, Glu94 was the main adsorption site in Model K, while Glu96 was the main adsorption site in Model W. It can be seen that the screening and optimization of the initial orientation and conformation play an important role in the simulation research of the interaction between proteins and material surfaces so that the results would be closer to reality.85-87 3.3 BMP-2 Desorption by Steered MD Simulation Although classical MD simulation could provide some insights into the binding performance of proteins on material surface, it could give contributions of individual residues to the binding affinity. Steered MD (SMD) approach is one of appropriate methods to solve this problem. In fact, the combination of SMD and MD have previously been applied to tackle this issue involved other BMPs on HAP.9, 44 Herein, the SMD was also employed to further reveal the adsorption mechanism of BMPs on HAP based on our newly constructed HAP/BMP-2 complex models. Surprisingly, the monomer model of BMP-2, which bears its biological activity in dimer form10-13, was widely adopted in previous studies. Lack of a complete picture of the protein, deviations should be expected between simulation results and practical experiments. Recently, Huang et al88 systematically studied the adsorption and desorption behaviors of BMP-2 on a nanostructured HAP surface by combining MD simulations and SMD simulations. The results showed that the nanostructured morphology of HAP surface played an important role in the dynamic behavior of BMP-2. However, their simulations provided no further binding information to the BMP-2 functional sites of knuckle and wrist epitopes, since the initial models in their work were prepared via an artificial positioning way. In this work, we will try to understand the behavior of those residues around the major binding epitopes during protein desorption process. On the basis of our MD results presented above, we can easily find out that the Model W is not stable. Number of residues in wrist epitope is even less than those in knuckle epitope during the MD simulation. Therefore, we then simply choose not to further investigate the desorption process of BMP-2 from HAP surface using the Model W as the initial structure. However, relatively stable material-protein complex model can be located in the MD simulation when knuckle epitope directly interacted with HAP (100) surface. Without loss of generality, we simply selected the last frame 16


Page 16 of 26

Page 17 of 26 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


structure from the MD trajectory in Model K as the initial model in our SMD simulation. The change of free energy in this process and the adsorption capacity of related sites were analyzed according to the potential of mean force (PMF). The results showed that the orientation of BMP-2 varies significantly in the period of 500 ps to 4000 ps, as shown in Figure 11A. The orientation gradually changed to upright from parallel to the HAP (100) surface at first. During the stage of 300-800 ps, two glutamic acids of Glu109 and Glu83 were desorbed from the surface. Due to the presence of two residues with negatively charged carboxylate group, the measured pulling force was quite large (Figure 11B). At the stage of 1500-2800 ps, the peak of

Figure 11. (A) Potential of mean force (PMF) and orientation evolution of BMP-2 during desorption process from the HAP (100) surface simulated using SMD approach for Model K. (B) Pulling forces with respect to Pulling Constant Velocity time upon SMD simulation for the Model K.

desorption appeared in Fig.11 revealed that Asp53 and Asp93 desorbed. After 3200 ps, the last adsorbed residue was Glu94, making the BMP-2 stick to the surface of HAP (Figure 11A&B). Surprisingly, this interaction was so strong that it could resist the external pulling force of SMD as long as 3ns. After 3 ns, the PMF and pulling force fluctuated around a fixed value, which indicated that BMP-2 was totally desorbed from the HAP surface and no longer had interactions with HAP. The SMD results showed that Glu94 had the strongest adsorption, which was consistent with the results of previous MD simulation86-87,

89.

Meanwhile, our simulation also suggested the

adsorption of BMP-2 on HAP was dominated by the acidic residues containing the 17



carboxylate group. It is worth noting that acidic residues account for about 20% in the knuckle epitope of BMP-2. Reasonably, the knuckle epitope had more stable and stronger adsorption on HAP (100) surface. Moreover, BMP-2 with knuckle epitope adsorbing was more likely to be the favorable orientation. It could allow BMP-2 wrist epitopes to be more exposed to environment, and thus have better chance to bind the type I receptors so that exert the bioactivity of osteoinduction. In summary, results presented above showed that the knuckle epitope region of BMP-2 could form much more stable interactions with the HAP (100) surface more stably, making the wrist epitope upward. This adsorption orientation was more conducive to that the wrist epitope exposed to specifically recognize type I receptors on the stem cell membranes, and exert biological activity (Figure 12). Thus, the biological function could be regulated through the surface modifications of biomaterials.

Figure 12. Graphical representation of the HAP (100) surface effects on the SMAD signaling pathway.

18


Page 18 of 26

Page 19 of 26 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


4. Conclusions Although high-throughput virtual screening method has been widely accepted in drug design,90-92 the applications of such a scenario are still rare46,

52-54

in material

design or investigation of interactions between solid surface and biological systems. In this work, we applied the screening technique, which is implemented in Rosetta software, to prepare the initial models for subsequent MD and SMD simulations. It is believed that such kind of automatic screening method could largely avoid the error brought by manually docking proteins on material surface. The application of RosettaSurface module can not only improve the efficiency and accuracy of materials genes screening, but also promote the development of materials genome on the theory and technology. Meanwhile, the usage of BMP-2 homodimer model and HAP IFF force field could provide more accurately addressing the functional roles of knuckle or wrist epitope with the help of implanted HAP bioceramics. The results showed that there was no obvious change in the conformation of the BMP-2 after docking, but the orientations changed a lot. The MD and SMD results also indicated that the HAP (100) surface prefer to adsorb the knuckle epitope instead of the wrist epitope. To this point, the wrist epitope or the active sites of type I could have much more opportunity to exposure to the environment and easily interact with type I receptors to exhibit its biological activity. All the calculations provided the molecular and atomic information to better understand the mechanism of osteoinduction from the perspective of protein adsorption. Of course, we have to point out only HAP (100) surface and BMP-2 homodimer as the target systems were applied in our work. Limited systems cannot provide complete information to address how HAP affects the absorbates or osteoinductive mechanisms. However, the applications of this high-throughput screening technique in the protein absorption on HAP surface provide a chance to extend to studies on more biomaterial forms and proteins. Furthermore, it might also provide materials genetic screening approach for the design and preparation of novel tissue induction materials.

Notes The authors declare no competing financial interest. 19



Acknowledgments This study was supported by the National Key Research and Development Program (no. 2016YFB0700801). Some of results described in this paper were obtained on the National Supercomputing Center of Guangzhou. Supporting Information The scoring energy of the 10 lowest docking structures. This material is available free of charge via the internet at http://pubs.acs.org/.

20


Page 20 of 26

Page 21 of 26 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


References (1) Fernández-Montes Moraleda, B.; Román, J. S.; Rodríguez-Lorenzo, L. M. Influence of Surface Features of Hydroxyapatite on the Adsorption of Proteins Relevant to Bone Regeneration. J. Biomed. Mater. Res. A 2013, 101A, 2332-2339. (2) Walkey, C. D.; Chan, W. C. W. Understanding and Controlling the Interaction of Nanomaterials with Proteins in a Physiological Environment. Chem. Soc. Rev. 2012, 41, 2780-2799. (3) Bhattacharjee, M.; Schultz-Thater, E.; Trella, E.; Miot, S.; Das, S.; Loparic, M.; Ray, A. R.; Martin, I.; Spagnoli, G. C.; Ghosh, S. The Role of 3D Structure and Protein Conformation on the Innate And adaptive Immune Responses to Silk-Based Biomaterials. Biomaterials 2013, 34, 8161-8171. (4) Molino, P. J.; Higgins, M. J.; Innis, P. C.; Kapsa, R. M. I.; Wallace, G. G. Fibronectin and Bovine Serum Albumin Adsorption and Conformational Dynamics on Inherently Conducting Polymers: A QCM-D Study. Langmuir 2012, 28, 8433-8445. (5) Rahmany, M. B.; Van Dyke, M. Biomimetic Approaches to Modulate Cellular Adhesion in Biomaterials: A Review. Acta Biomaterialia 2013, 9, 5431-5437. (6) Derner, R.; Anderson, A. C. The Bone Morphogenic Protein. Clin. Podiat. Med. Sur. 2005, 22, 607-618. (7) Na, K.; Kim, S. W.; Sun, B. K.; Woo, D. G.; Yang, H. N.; Chung, H. M.; Park, K. H. Osteogenic Differentiation of Rabbit Mesenchymal Stem Cells in Thermo-Reversible Hydrogel Constructs Containing Hydroxyapatite and Bone Morphogenic Protein-2 (BMP-2). Biomaterials 2007, 28, 2631-2637. (8) Kato, M.; Namikawa, T.; Terai, H.; Hoshino, M.; Miyamoto, S.; Takaoka, K. Ectopic Bone Formation in Mice Associated with a Lactic Acid/Dioxanone/Ethylene Glycol Copolymer-Tricalcium Phosphate Composite with Added Recombinant Human Bone Morphogenetic Protein-2. Biomaterials 2006, 27, 3927-3933. (9) Dong, X.; Wang, Q.; Wu, T.; Pan, H. Understanding Adsorption-Desorption Dynamics of Bmp-2 on Hydroxyapatite (001) Surface. Biophys. J. 2007, 93, 750-759. (10) Keller, S.; Nickel, J.; Zhang, J.; Sebald, W.; Mueller, T. D. Molecular Recognition of BMP-2 and BMP Receptor IA. Nat. Struct. Mol. Biol. 2004, 11, 481-488. (11) Allendorph, G. P.; Vale, W. W.; Choe, S. Structure of the Ternary Signaling Complex of a TGF- Superfamily Member. Proc. Natl. Acad. Sci. USA 2006, 103, 7643-7648. (12) Kirsch, T.; Sebald, W.; Dreyer, M. K. Crystal Structure of the BMP-2–BRIA Ectodomain Complex. Nat. Struct. Biol. 2000, 7, 492-496. (13) Kirsch, T.; Nickel, J.; Sebald, W. BMP-2 Antagonists Emerge from Alterations in the Low-Affinity Binding Epitope for Receptor BMPR-Ⅱ. EMBO J. 2000, 19, 3314-3324. (14) Kaneko, H.; Arakawa, T.; Mano, H.; Kaneda, T.; Ogasawara, A.; Nakagawa, M.; Toyama, Y.; Yabe, Y.; Kumegawa, M.; Hakeda, Y. Direct Stimulation of Osteoclastic Bone Resorption by Bone Morphogenetic Protein (BMP)-2 and Expression of BMP Receptors in Mature Osteoclasts. Bone 2000, 27, 479-486. (15) Weber, D.; Kotzsch, A.; Nickel, J.; Harth, S.; Seher, A.; Mueller, U.; Sebald, W.; Mueller, T. D. A Silent H-Bond Can Be Mutationally Activated for High-Affinity Interaction of BMP-2 and Activin Type ⅡB Receptor. BMC Struct. Biol. 2007, 7, 1-20. (16) Koenig, B. B.; Cook, J. S.; Wolsing, D. H.; Ting, J.; Tiesman, J. P.; Correa, P. E.; Olson, C. A.; Pecquet, A. L.; Ventura, F.; Grant, R. A et al. Characterization and Cloning of a Receptor for BMP-2 and BMP-4 from NIH 3T3 Cells. Mol. Cell. Biol. 1994, 14, 5961-5974. (17) Kim, S. R.; Lee, J. H.; Kim, Y. T.; Riu, D. H.; Jung, S. J.; Lee, Y. J.; Chung, S. C.; Kim, Y. H. Synthesis of Si, Mg Substituted Hydroxyapatites and Their Sintering Behaviors. Biomaterials 2003, 24, 1389-1398. (18) Barroug, A.; Glimcher, M. J. Hydroxyapatite Crystals as a Local Delivery System for Cisplatin: Adsorption and Release of Cisplatin in Vitro. J. Orthop. Res. 2002, 20, 274-280. 21



(19) Al Malyan, M.; Becchi, C.; Boncinelli, S.; Ashammakhi, N. Novel Drug Delivery Systems in Pain Therapy. Minerva. Anestesiol. 2007, 73, 173-179. (20) Yuan, H. P.; Li, Y. B.; de Bruijn, J. D.; de Groot, K.; Zhang, X. D. Tissue Responses of Calcium Phosphate Cement: A Study in Dogs. Biomaterials 2000, 21, 1283-1290. (21) Yuan, H. P.; Yang, Z. J.; Li, Y. B.; Zhang, X. D.; De Bruijn, J. D.; De Groot, K. Osteoinduction by Calcium Phosphate Biomaterials. J. Mater Sci-Mater. M. 1998, 9, 723-726. (22) Yuan, H. P.; Zou, P.; Yang, Z. J.; Zhang, X. D.; De Bruijn, J. D.; De Groot, K. Bone Morphogenetic Protein and Ceramic-Induced Osteogenesis. J. Mater. Sci- Mater. M. 1998, 9, 717-721. (23) Yang, Z. J.; Yuan, H. P.; Tong, W. D.; Zou, P.; Chen, W. Q.; Zhang, X. D. Osteogenesis in Extraskeletally Implanted Porous Calcium Phosphate Ceramics: Variability among Different Kinds of Animals. Biomaterials 1996, 17, 2131-2137. (24) K, Bessho.; David L, Carnes.; R, Cavin.; H, Chen.; Joo L, Ong. BMP Stimulation of Bone Response Adjacent to Titanium Implants in Vivo. Clin. Oral. Implan. Res. 2001, 10, 212-218. (25) Mücksch, C.; Urbassek, H. M. Adsorption of BMP-2 on a Hydrophobic Graphite Surface: A Molecular Dynamics Study. Chem. Phys. Lett. 2011, 510, 252-256. (26) Utesch, T.; Daminelli, G.; Mroginski, M. A. Molecular Dynamics Simulations of the Adsorption of Bone Morphogenetic Protein-2 on Surfaces with Medical Relevance. Langmuir 2011, 27, 13144-13153. (27) Begam, H.; Nandi, S. K.; Chanda, A.; Kundu, B. Effect of Bone Morphogenetic Protein on Zn-Hap and Zn-Hap/Collagen Composite: A Systematic in Vivo Study. Res. Vet. Sci. 2017, 115, 1-9. (28) Shi, F.; Liu, Y.; Zhi, W.; Xiao, D.; Li, H.; Duan, K.; Qu, S.; Weng, J. The Synergistic Effect of Micro/Nano-Structured and Cu2+-Doped Hydroxyapatite Particles to Promote Osteoblast Viability and Antibacterial Activity. Biomed. Mater. 2017, 12, 1-12. (29) Crouzier, T.; Sailhan, F.; Becquart, P.; Guillot, R.; Logeart-Avramoglou, D.; Picart, C. The Performance of BMP-2 Loaded TCP/HAP Porous Ceramics with a Polyelectrolyte Multilayer Film Coating. Biomaterials 2011, 32, 7543-7554. (30) Long, J. R.; Shaw, W. J.; Stayton, P. S.; Drobny, G. P. Structure and Dynamics of Hydrated Statherin on Hydroxyapatite as Determined by Solid-State Nmr. Biochemistry-US 2001, 40, 15451-15455. (31) Dolatshahi-Pirouz, A.; Jensen, T.; Foss, M.; Chevallier, J.; Besenbacher, F. Enhanced Surface Activation of Fibronectin Upon Adsorption on Hydroxyapatite. Langmuir 2009, 25, 2971-2978. (32) Zhang, H. P.; Lu, X.; Fang, L. M.; Qu, S. X.; Feng, B.; Weng, J. Atomic-Scale Interactions at the Interface of Biopolymer/Hydroxyapatite. Biomed. Mater. 2008, 3, 1-5. (33) Ishibe, T.; Goto, T.; Kodama, T.; Miyazaki, T.; Kobayashi, S.; Takahashi, T. Bone Formation on Apatite-Coated Titanium with Incorporated BMP-2/Heparin in Vivo. Oral Surg. Oral Med. O. 2009, 108, 867-875. (34) Boix, T.; Gómez-Morales, J.; Torrent-Burgués, J.; Monfort, A.; Puigdomènech, P.; Rodríguez-Clemente, R. Adsorption of Recombinant Human Bone Morphogenetic Protein RhBMP-2m onto Hydroxyapatite. J. Inorg. Biochem. 2005, 99, 1043-1050. (35) Liao, C.; Zhou, J. Replica-Exchange Molecular Dynamics Simulation of Basic Fibroblast Growth Factor Adsorption on Hydroxyapatite. J. Phys. Chem. B 2014, 118, 5843-5852. (36) Ono, I.; Tateshita, T.; Takita, H.; Kuboki, Y. Promotion of the Osteogenetic Activity of Recombinant Human Bone Morphogenetic Protein by Basic Fibroblast Growth Factor. J. Craniofac. surg. 1996, 7, 418-425. (37) Tsurushima, H.; Marushima, A.; Suzuki, K.; Oyane, A.; Sogo, Y.; Nakamura, K.; Matsumura, A.; Ito, A. Enhanced Bone Formation Using Hydroxyapatite Ceramic Coated with Fibroblast Growth Factor-2. Acta Biomater. 2010, 6, 2751-2759. 22


Page 22 of 26

Page 23 of 26 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


(38) Dong, X. L.; Qi, W.; Tao, W.; Ma, L. Y.; Fu, C. X. The Dynamic Behaviours of Protein Bmp-2 on Hydroxyapatite Nanoparticles. Mol. Simulat. 2011, 37, 1097-1104. (39) Mücksch, C.; Urbassek, H. M. Enhancing Protein Adsorption Simulations by Using Accelerated Molecular Dynamics. PLOS ONE 2013, 8, 1-7. (40) Almora-Barrios, N.; de Leeuw, N. H. A Density Functional Theory Study of the Interaction of Collagen Peptides with Hydroxyapatite Surfaces. Langmuir 2010, 26, 14535-14542. (41) 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. (42) Jain, A.; Ong, S. P.; Hautier, G.; Chen, W.; Richards, W. D.; Dacek, S.; Cholia, S.; Gunter, D.; Skinner, D.; Ceder, G.; Persson, K. A. Commentary: The Materials Project: A Materials Genome Approach to Accelerating Materials Innovation. APL Mater. 2013, 1, 011002. (43) Shen, J. W.; Wu, T.; Wang, Q.; Pan, H. H. Molecular Simulation of Protein Adsorption and Desorption on Hydroxyapatite Surfaces. Biomaterials 2008, 29, 513-532. (44) Zhou, H.; Wu, T.; Dong, X.; Wang, Q.; Shen, J. Adsorption Mechanism of BMP-7 on Hydroxyapatite (001) Surfaces. Biochem. Biophys. Res. Co. 2007, 361, 91-96. (45) Liao, C.; Xie, Y.; Zhou, J. Computer Simulations of Fibronectin Adsorption on Hydroxyapatite Surfaces. Rsc Adv. 2014, 4, 15759-15769. (46) Zhao, W.; Xu, Z.; Cui, Q.; Sahai, N. Predicting the Structure-Activity Relationship of Hydroxyapatite-Binding Peptides by Enhanced-Sampling Molecular Simulation. Langmuir 2016, 32, 7009-7022. (47) Daggett, V. Protein Folding-Simulation. Chem. Rev. 2006, 106, 1898-1916. (48) 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. 2013, 35, 70-81. (49) Pan, H. H.; Tao, J. H.; Xu, X. R.; Tang, R. K. Adsorption Processes of Gly and Glu Amino Acids on Hydroxyapatite Surfaces at the Atomic Level. Langmuir 2007, 23, 8972-8981. (50) Wilson, R. M.; Elliott, J. C.; Dowker, S. E. P. Rietveld Refinement of the Crystallographic Structure of Human Dental Enamel Apatites. Am. Mineral. 1999, 84, 1406-1414. (51) Chiatti, F.; Delle Piane, M.; Ugliengo, P.; Corno, M. Water at Hydroxyapatite Surfaces: The Effect of Coverage and Surface Termination as Investigated by All-Electron B3LYP-D* Simulations. Theor. Chem. Acc. 2016, 135, 54. (52) Gray, J. The Interaction of Proteins with Solid Surfaces. Curr. Opin. Struc. Biol. 2004, 14, 110-115. (53) Makrodimitris, K.; Masica, D. L.; Kim, E. T.; Gray, J. J. Structure Prediction of Protein-Solid Surface Interactions Reveals a Molecular Recognition Motif of Statherin for Hydroxyapatite. J. Am. Chem. Soc. 2007, 129, 13713-13722. (54) Masica, D. L.; Gray, J. J. Solution- and Adsorbed-State Structural Ensembles Predicted for the Statherin-Hydroxyapatite System. Biophys. J. 2009, 96, 3082-3091. (55) Addison, W. N.; Miller, S. J.; Ramaswamy, J.; Mansouri, A.; Kohn, D. H.; McKee, M. D. Phosphorylation-Dependent Mineral-Type Specificity for Apatite-Binding Peptide Sequences. Biomaterials 2010, 31, 9422-9430. (56) Hauptmann, S.; Dufner, H.; Brickmann, J.; M. Kast, S.; Berry, R. Potential Energy Function for Apatites. Phys. Chem. Chem. Phys. 2003, 5, 635-639. (57) Lin, T.; 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. (58) Heinz, H.; Lin, T. J.; Kishore Mishra, R.; Emami, F. S. Thermodynamically Consistent Force Fields for the Assembly of Inorganic, Organic, and Biological Nanostructures: The Interface Force Field. Langmuir 2013, 29, 1754-1765. 23



(59) Xie, Q.; Xue, Z.; Gu, H.; Hu, C.; Yang, M.; Wang, X.; Xu, D. Molecular Dynamics Exploration of Ordered-to-Disordered Surface Structures of Biomimetic Hydroxyapatite Nanoparticles. J. Phys. Chem. C 2018, 122, 6691-6703. (60)Mackerell, A.D.; Feig, M.; Brooks, C. L. Extending the Treatment of Backbone Energetics in Protein Force Fields: Limitations of Gas-Phase Quantum Mechanics in Reproducing Protein Conformational Distributions in Molecular Dynamics Simulations. J. Comput. Chem. 2004, 25, 1400-1415. (61) MacKerell, A.D.; D, Bashford ; M, Bellott ; R L, Dunbrack ; J D, Evanseck ; M J, Field ; S, Fischer ; J, Gao ; H, Guo ; S, Ha et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586-3616. (62) Bjelkmar, P.; Larsson, P.; Cuendet, M. A.; Hess, B.; Lindahl, E. Implementation of the Charmm Force Field in Gromacs: Analysis of Protein Stability Effects from Correction Maps, Virtual Interaction Sites, and Water Models. J. Chem. Theory Comput. 2010, 6, 459-466. (63) Jorgensen, W. L.; Chandrasekhar, J.; Madura, J.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926-935. (64) Pacella, M. S.; Koo, D. C. E.; Thottungal, R. A.; Gray, J. J. Using the Rosettasurface Algorithm to Predict Protein Structure at Mineral Surfaces. Method. Enzymol. 2013, 532, 343-366. (65) Rohl, C. A.; Strauss, C. E. M.; Misura, K. M. S.; Baker, D. Protein Structure Prediction Using Rosetta. Method. Enzymol. 2004, 383, 66-93. (66) Kabsch, W.; Sander, C. Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22, 2577-2637. (67) Jain, A. K.; Murty, M. N.; Flynn, P. J. Data Clustering: A Review. Acm. Comput. Surv. 1999, 31, 264-323. (68) Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1-19. (69) 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. (70) Kelly, D. W. The Self-Equilibration of Residuals and Complementary a Posteriori Error Estimates in the Finite Element Method. Int. J. Numer. Meth. Eng. 1984, 20, 1491-1506. (71) Nosé, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511-519. (72) Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A. 1985, 31, 1695-1697. (73) Hockney, R. W.; Eastwood, J. W. Computer Simulation Using Particles. Siam. Res. 1988, 25, 540-540. (74) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14, 33-38. (75) Park, S.; Schulten, K. Calculating Potentials of Mean Force from Steered Molecular Dynamics Simulations. J. Chem. Phys. 2004, 120, 5946-5961. (76) Jarzynski, C. Nonequilibrium Equality for Free Energy Differences. Phys. Rev. Lett. 1997, 78, 2690-2693. (77) Peter, D. Computational Molecular Dynamics: Challenges, Methods, Ideas: Proceedings of the 2nd International Symposium on Algorithms for Macromolecular Modellin; Springer: Berlin Heidelberg, 1999, 39-65. (78) Xu, Z. J.; Yang, X. N.; Yang, Z. On the Mechanism of Surfactant Adsorption on Solid Surfaces: Free-Energy Investigations. J. Phys. Chem. B 2008, 112, 13802-13811. (79) Park, S.; Khalili-Araghi, F.; Tajkhorshid, E.; Schulten, K. Free Energy Calculation from Steered Molecular Dynamics Simulations Using Jarzynski’s Equality. J. Chem. Phys. 2003, 119, 3559-3566. 24


Page 24 of 26

Page 25 of 26 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


(80) Trzesniak, D.; Kunz, A. E.; van Gunsteren, W. F. A Comparison of Methods to Compute the Potential of Mean Force. ChemPhysChem 2006, 8, 162-169. (81) Jäger, C.; Welzel, T.; Meyer-Zaika, W.; Epple, M. A Solid-State Nmr Investigation of the Structure of Nanocrystalline Hydroxyapatite. Magn. Reson. Chem. 2006, 44, 573-580. (82) Bertinetti, L.; Tampieri, A.; Landi, E.; Ducati, C.; Midgley, P. A.; Coluccia, S.; Martra, G. Surface Structure, Hydration, and Cationic Sites of Nanohydroxyapatite: UHR-TEM, IR, and Microgravimetric Studies. J. Phys. Chem. C 2007, 111, 4027-4035. (83) Forget, A.; Christensen, J.; Lüdeke, S.; Kohler, E.; Tobias, S.; Matloubi, M.; Thomann, R.; Shastri, V. P. Polysaccharide Hydrogels with Tunable Stiffness and Provasculogenic Properties Via -Helix to -Sheet Switch in Secondary Structure. Proc. Natl. Acad. Sci. 2013, 110, 12887. (84) Stanger, H. E.; Syud, F. A.; Espinosa, J. F.; Giriat, I.; Muir, T.; Gellman, S. H. Length-Dependent Stability and Strand Length Limits in Antiparallel -Sheet Secondary Structure. Proc. Natl. Acad. Sci. 2001, 98, 12015. (85) Price, P. A.; Poser, J. W.; Raman, N. Primary Structure of the Gamma-Carboxyglutamic Acid-Containing Protein from Bovine Bone. Proc. Natl Acad Sci 1976, 73, 3374. (86) Kawasaki, T.; Niikura, M.; Kobayashi, Y. Fundamental Study of Hydroxyapatite High-Performance Liquid Chromatography: III. Direct Experimental Confirmation of the Existence of Two Types of Adsorbing Surface on the Hydroxyapatite Crystal. J. Chromatogr. A 1990, 515, 125-148. (87) Kawasaki, T. Hydroxyapatite as a Liquid Chromatographic Packing. J. Chromatogr. A 1991, 544, 147-184. (88) Huang, B. L.; Lou, Y.; Li, T.; Lin, Z.; Sun, S. L.; Yuan, Y.; Liu, C. S.; Gu, Y. T. Molecular Dynamics Simulations of Adsorption and Desorption of Bone Morphogenetic Protein-2 on Textured Hydroxyapatite Surfaces. Acta Biomater. 2018, 80, 121-130. (89) Hunter, G. K.; O’Young, J.; Grohe, B.; Karttunen, M.; Goldberg, H. A. The Flexible Polyelectrolyte Hypothesis of Protein-Biomineral Interaction. Langmuir 2010, 26, 18639-18646. (90) Bleicher, K. H.; Böhm, H. J.; Müller, K.; Alanine, A. I. Hit and Lead Generation: Beyond High-Throughput Screening. Nat. Rev. Drug Discov. 2003, 2, 369. (91) Blundell, T. L.; Jhoti, H.; Abell, C. High-Throughput Crystallography for Lead Discovery in Drug Design. Nat. Rev. Drug Discov. 2002, 1, 45. (92) Bajorath, J. Integration of Virtual and High-Throughput Screening. Nat. Rev. Drug Discov. 2002, 1, 882.

25



TOC Graphic

26


Page 26 of 26

Effect of Hydroxyapatite Surface on BMP-2 Biological Properties by

Recommend Documents