Forced Desorption of Bovine Serum Albumin and Lysozyme from

Jul 15, 2016 - Forced Desorption of Bovine Serum Albumin and Lysozyme from Graphite: Insights from Molecular Dynamics Simulation. Christian Mücksch ...
0 downloads 0 Views 5MB Size
Subscriber access provided by United Arab Emirates University | Libraries Deanship

Article

Forced Desorption of BSA and Lysozyme from Graphite: Insights from Molecular Dynamics Simulation Christian Mücksch, and Herbert M. Urbassek J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.6b05234 • Publication Date (Web): 15 Jul 2016 Downloaded from http://pubs.acs.org on July 17, 2016

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry B 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 22

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

Forced Desorption of BSA and Lysozyme from Graphite: Insights from Molecular Dynamics Simulation Christian M¨ucksch and Herbert M. Urbassek∗ Physics Department and Research Center OPTIMAS, University Kaiserslautern, Erwin-Schr¨odinger-Straße, D-67663 Kaiserslautern, Germany E-mail: [email protected],Phone:+49(0)6312053022



To whom correspondence should be addressed

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 We use molecular dynamics (MD) simulation to study the adsorption and desorption of two widely different proteins – BSA and lysozyme – on a graphite surface. The adsorption is modeled using accelerated MD to allow the proteins to find optimum conformations on the surface. Our results demonstrate that the “hard protein” lysozyme retains much of its secondary structure during adsorption, while BSA loses it almost completely. BSA has a considerably larger adsorption energy compared to lysozyme, which does not scale with chain length. Desorption simulations are carried out using classical steered MD. The BSA molecule becomes fully unzipped during pull-off, while several helices survive this process in lysozyme. The unzipping process shows up in the force-distance curve of BSA as a series of peaks while only a single – or few, depending on protein orientation – force peaks occur for lysozyme. The maximum desorption force is larger for BSA than for lysozyme, but only by a factor of about 2.3.

Introduction When a protein approaches and adheres to a surface several outcomes of this process called adsorption are possible. The adsorbed protein can be in various conformation states from native to unfolded, depending on the hydrophilicity of the surface, charge distribution, surface topography and protein structure. 1 This nonspecific process and the following protein-cell interactions determine the surfaces biocompatibility. Graphite which is widely used in its pyrolytic form 2 as an implant material due to its physical properties 3 displays a high biocompatibility. Therefore, we study the adsorption and forced desorption of two very different proteins, BSA and lysozyme, from a graphite surface. While the adsorption of proteins on inorganic surfaces has been studied extensively using simulation in the past 4–11 considerably less research has been done on the forced desorption induced by pulling the protein with an AFM tip off the surface. 7,8 The equivalent bioadhesion experiments using desorption can be performed with the use of scanning force spectroscopy (SFS). 12,13 2

ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22

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 problems that are posed when simulating the adsorption process are well known and have been described in several publications. 10,14,15 (i) Long timescale of simulation are needed in order to accommodate the protein in its optimum conformation on the surface; 11 (ii) large simulation box sizes are essential to comprise all the water. Problem (i) can be solved by using accelerated MD (aMD). 16 For problem (ii) sometimes so-called implicit water simulations are performed where the water molecules are not included in the simulation, but only the dielectric constant of water and mostly not even the viscous drag exerted by it is included in the equations of motion of the protein. 9 It is not clear whether this simplification is adequate. 17 These difficulties become even larger when performing pull-off simulations. (i) In experiment, pulling velocities are of the order of 10−6 m/s, such that an extension of BSA to its full estimated chain length of about 210 nm would require simulation times of around 0.2 s. In consequence, all simulations are usually performed with much larger velocities, around 1–10 m/s. We note that recently 18 a scheme was proposed which incorporates aMD into a pull-off simulation. We refrained from using this scheme here, since its optimum realization requires the use of several pulling velocities, which is – for the case of BSA – beyond our present computational capabilities. (ii) Since during desorption the protein is strongly lengthened, correspondingly long simulation boxes are needed. This prevented up to now the simulation of the desorption of large molecules such as BSA. We chose two widely different proteins, BSA and lysozyme, for our study. These differ in several respects: (i) size, which can be attributed either to the number of amino acids or to the radius of gyration describing the protein extension; (ii) unfolding tendency or hardness which is defined by the protein’s internal stability with BSA being described as a “soft protein” and lysozyme as a “hard protein”. 19 In the present paper we perform adsorption simulations of BSA and lysozyme for 40 ns using accelerated MD on a hydrophobic graphite surface. We consider this simulation time as adequate for describing an initial adsorption on an implant surface with respect to the

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

computational limitations that allow the following desorption simulations at pulling velocities of 10 m/s.

Simulation Methods

Figure 1: Orientation 1 of BSA above the graphite surface in 0.1 M NaCl. Nonpolar residues are shown in blue, neutral residues in green, and polar residues in red. Water molecules are not shown for clarity while sodium ions are shown in yellow and the chloride ions in cyan. For performing molecular dynamics simulations NAMD 2.9 20 was used together with the CHARMM27 force field 21 and the TIP3P water model. 22 Throughout all simulations, the temperature was set to 310 K and the pressure to 1 atm. Prior to both adsorption and desorption simulations, energy minimization and equilibration were carried out for 200 ps to stabilize the thermodynamic properties. A 2 fs time step was enabled by employing the SHAKE algorithm. 23 Periodic boundary conditions in combination with the particle-meshEwald (PME) routine 24 were used to handle the long-range electrostatics. The structures of BSA and lysozyme were taken from the protein data base. In the case of BSA we now used structure 4F5S 25 at pH 6.5 with a net charge of −17 e that was neutralized at 0.1 M NaCl. This setup above the graphite surface in a NaCl solution is shown in Fig. 1. In our previous work 26 we had to use a homology model since structural data was not available back then. Structural data 1AKI 27 for hen egg-white lysozyme at pH of 4.48 4

ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22

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

with a net charge of +8 e was also neutralized at 0.1 M NaCl. Each of the proteins was simulated in 2 different orientations where orientation 2 was obtained from orientation 1 by a 90◦ rotation around an axis parallel to the surface. The initial orientations for the two proteins are shown in Fig. 2. From the total residue count of the structures (583 amino acids for BSA and 129 amino acids for lysozyme) one can estimate the protein chain length of the BSA structure to 210 nm while the lysozyme structure reaches about 46 nm. All of the snapshots were rendered using VMD 28 and Tachyon. 29

Adsorption Since hydrocarbon contamination is believed to be responsible for making graphite more hydrophobic 30,31 but we do not know exactly what kind of hydrocarbons are present on the surface (compare 32 ), we model the surface by using an effective hydrophobicity that leads to a macroscopic contact angle of 80◦ . This angle was observed in previous contact-angle experiments after 3h of exposure to air, 32 which we consider a typical time scale in adsorption experiment preparations. The graphite surface is described by the Lennard–Jones parameters of Rmin /2 = 1.763388 ˚ A as provided by Werder et al. 33 and a scaled ǫ = −0.056485 kcal/mol to achieve the desired hydrophobicity. Both proteins start the adsorption phase with a minimum distance of 6 ˚ A to the graphite surface. In the case of BSA the simulation box had dimensions of 379 × 435 × 102 ˚ A with a total of about 5 · 105 water molecules. Since lysozyme is much smaller, it was placed inside a simulation box of 207 × 237 × 74 ˚ A with a total of about 1 · 105 water molecules. To account for long-term adsorption phenomena the implementation of dual accelerated MD 34–36 in NAMD 37 was used as we did previously. 16 The acceleration parameters were adapted for the different proteins since they contain different amounts of residues and water molecules but were not varied between the two orientations.

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

Figure 2: Initial orientations of BSA and lysozyme above the graphite surface. Orientation 2 was obtained from orientation 1 by a 90◦ rotation around an axis parallel to the surface. Nonpolar residues shown in blue, neutral residues in green and polar residues in red. The N-terminus is shown in yellow and the C-terminus in orange. Note that the scale is identical in all subfigures.

6

ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22

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

Desorption In order to fully desorb the proteins from their adsorbed state on graphite very long simulation boxes in pulling direction were needed with lateral dimensions reduced to the sizes of the adsorbed state. For BSA these account to 211 × 211 × 2254 ˚ A with a total of about 3.3 · 106 water molecules, and for lysozyme to 90 × 90 × 443 ˚ A with a total of about 1 · 105 water molecules. These box sizes were estimated from the respective chain lengths of the unfolded protein structures which were ≈ 210 nm for BSA and ≈ 46 nm for lysozyme. Steered molecular dynamics (SMD) simulations with a desorption velocity of 10 m/s and a spring constant of 0.03 N/m typical of scanning force spectroscopy experiments were employed. The fast pulling velocity is needed to desorb the proteins with a reasonable computational effort, although SMD simulations for BSA required up to 5 · 105 CPU hours for a single orientation to complete. The simulation times ranged between 6 ns for lysozyme and 25 ns for BSA until complete desorption. As an attachment for the virtual cantilever the top nitrogen atom was selected as the pull group.

Results and Discussion Adsorption The adsorption process of BSA and lysozyme on graphite is simulated using accelerated molecular dynamics for a total of 2 · 107 time steps. Although the two proteins show a different adsorption behavior both start by maximizing their surface contact which is seen in the attractive van-der-Waals interactions with graphite in the left part of Fig. 3. The radius of gyration with its component parallel to the surface shown in Fig. 4 gives us an idea how much the protein is able to spread out on the surface. The initial orientation plays an important role as it decides on the flexibility of the protein for spreading on the surface. Our simulations show that lysozyme as a “hard protein”, meaning a high internal

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

stability, 19 has little tendency to unfold and spread on the graphite surface which is also seen in the adsorption snapshots Fig. 5. Our results are in contrast to the work of Raffaini et al. 9 that demonstrated complete unfolding of lysozyme on graphite. The reason hereto might lie in the different surface parameters used, our accelerated explicit water model vs. the implicit inviscid water model used by Raffaini et al. in combination with minimization routines and mainly our reduced amount of orientations due to the computational demands of the pull-off simulations which resulted in different starting configurations for the adsorption process. BSA on the other hand shows quite a strong reduction in secondary structure (see Table 1) and a small but noticeable amount of spreading especially for orientation 1 (see Fig. 4). Also the adsorption snapshots in Fig. 5 in comparison with Fig. 2 show the reduced secondary structure and spreading across the graphite surface. There has recently been another molecular dynamics investigation on the adsorption of BSA on graphene by Vilhena et al. 38 showing minimal structural changes. A comparison with our previous work in an implicit inviscid water environment regarding the secondary structure content 26 shows that it is very similar to that studied here (see Table 1) using accelerated molecular dynamics. Therefore, we argue that the 150 ns simulations in explicit solvent performed by Vilhena et al. 38 might suffer from insufficient sampling. Even for a small protein like lysozyme it was shown that a simulation time of 300 ns using classical molecular dynamics can be insufficient to completely describe the adsorption process. 11 Furthermore, differences in the surface parameters might also cause an altered adsorption outcome. When comparing these two proteins regarding their adsorption phase it becomes obvious that BSA shows a higher affinity towards graphite than lysozyme, especially by the much increased amount of unfolding shown in Table 1 and Table 2 and the relative change in the protein’s radius of gyration, Fig. 4. It is known that protein stability in the vicinity of surfaces depends generally on several influence factors, among which are direct surface-protein interactions, water mediated interactions (depending on surface hydrophobicity) and entropic effects. 39 In our case, especially due to the different lengths and hydrophobic inner parts

8

ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22

of the two proteins considered here, competing factors like intrapeptide vs. protein-surface interactions affect the adsorption outcome but complicate a detailed analysis. Interestingly, regarding surface interactions, we find a direct correlation regarding the area on which the proteins spread and the van-der-Waals interactions which both differ by a factor of about 7 between BSA and lysozyme. 0

Orientation 1 Orientation 2

-200

E (kcal/mol)

-400

-600

-800

-1000 adsorption

desorption

-1200 0*100

4*106

8*106

12*106

0 20*106

16*106

5

MD timesteps

10

15

20

25

Time (ns)

0

Orientation 1 Orientation 2

-20 -40

E (kcal/mol)

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

-60 -80 -100 -120 -140 adsorption

desorption

-160 0*100

4*106

8*106

12*106

16*106

MD timesteps

0 2 20*106

4

6

Time (ns)

Figure 3: Van-der-Waals energies of BSA (top) and lysozyme (bottom) with the graphite surface.

9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

30

25

20

°

Rg,|| (A)

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 10 of 22

15

BSA Orientation 1 BSA Orientation 2 Lysozyme Orientation 1 Lysozyme Orientation 2

10

5 0 0*10

4*10

6

6

6

8*10 12*10 MD timesteps

16*10

6

20*10

6

Figure 4: Evolution of radius of gyration – component parallel to the graphite surface.

Table 1: Secondary structure of BSA before adsorption (‘initial’) and at the end of the adsorption (‘ads’) and desorption (‘des’) simulations at zero van-der-Waals interactions, 2 · 107 time steps, calculated with DSSP. 40 after 2 · 107 time steps

initial minimized structure α-helix 310 -helix π-helix β-sheet

69.5 % 1.9 % 0% 0%

Orientation 1 ads 23.5 % 1.7 % 0% 0.5 %

Orientation 2 ads 23.3 % 3.8 % 0.5 % 1.0 %

Orientation 1 des 0% 0% 0% 12.4 %

Orientation 2 des 1.3 % 0.1 % 0% 0%

Table 2: Secondary structure of lysozyme before adsorption (‘initial’) and at the end of the adsorption (‘ads’) and desorption (‘des’) simulations at zero van-derWaals interactions, 2 · 107 time steps, calculated with DSSP. 40 after 2 · 107 time steps

initial minimized structure α-helix 310 -helix π-helix β-sheet

28.4 % 6.5 % 3.1 % 10.2 %

Orientation 1 ads 29.3 % 1.1 % 3.6 % 7.1 %

10

Orientation 2 ads 26.7 % 2.3 % 0% 8.4 %

ACS Paragon Plus Environment

Orientation 1 des 10.7 % 4.5 % 0% 1.6 %

Orientation 2 des 19.4 % 1.9 % 0% 1.6 %

Page 11 of 22

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 5: Topview snapshots of BSA and lysozyme after 2 · 107 time steps of accelerated adsorption. The color coding is explained in Fig. 2. The scale is identical in all subfigures.

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

Desorption The adsorbed protein structures where then subjected to an external force by means of steered molecular dynamics for up to 25 ns in the case of BSA and up to 6 ns for lysozyme until complete desorption did occur. We define complete desorption by zero Lennard-Jones interactions between protein and surface. Desorption lengths and desorption forces are listed in Table 3. The force distance curves in Fig. 6 show a similar pattern for both proteins: the first stage consists of a steep rise in force until the hindering secondary structure is unfolded and the mainly frictional forces of the solvent due to the fast pulling velocity are overcome to allow for further desorption. The high desorption forces for both proteins show that adsorption to the graphite surface is indeed quite strong. Interestingly, in the first stage between 0.5 to 1 ns of desorption all orientations show an increased interaction with the surface (see Fig. 3) which is induced by the external force acting on the surface-distant site of the protein, i.e., at the pull-group that allows the protein to cross an energetic barrier that was previously hindering further adsorption. BSA and lysozyme exhibit quite a different desorption behavior. BSA has spread out much more during adsorption and also due to its bigger size has more surface contact and therefore a stronger binding to the surface. Furthermore, BSA being known as a “soft protein” with a low internal stability 19 shows an almost complete unfolding during desorption also seen in the right part of Table 1. The unzipping of the BSA molecule during desorption can also be seen in Fig. 7 which shows the protein structures at their maximum desorption force. At full desorption BSA is unfolded to nearly its entire polypeptide chain length also shown in Table 3. Lysozyme on the other hand keeps at least some of its structure during the desorption process. Here, desorption occurs faster due to the smaller overall size and the reduced surface interactions. The force-distance curves for lysozyme in Fig. 6 also show less features than for BSA which is mainly due to the much shorter pulling distance. Interestingly, lysozyme has about a 4.5-fold shorter chain length than BSA but needs only about a tenth of BSA’s 12

ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22

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

desorption length. The reduced adsorption and fast desorption process while maintaining at least some of its secondary structure fit to the picture of lysozyme being a “hard protein” with high internal stability. 19 The relatively high unfolding forces of lysozyme can be attributed to the viscous forces due to the high pulling velocity. Table 3: Desorption force, i.e. maximum force in the force-distance curves (Fig. 6), and desorption length i.e., distance between surface and pulled protein atom at instant of vanishing van-der-Waals interaction between protein and graphite surface.

desorption force desorption length

BSA Orientation 1 Orientation 2 1.390 nN 1.415 nN 172.3 nm 203.9 nm

lysozyme Orientation 1 Orientation 2 0.590 nN 0.606 nN 21.5 nm 17.4 nm

Conclusions We study the adsorption and forced desorption of two widely different proteins, BSA and lysozyme, from a graphite surface. We took care that all parameters in the simulation – such as water environment, neutralizing salt content, surface conditions, pulling velocities – were identical for both molecules. The adsorption is performed using aMD in an explicit solvent since protein adsorption is a process that requires advanced sampling due to metastable states occurring while the protein is adapting to the surface. That allows the protein to optimize its conformation on the surface and to reach a state of high adsorption energy; this is a necessary prerequisite to obtain a good starting conformation for the desorption simulation. Previous studies already showed the advantage of using aMD over classical MD in adsorption studies. 16 Our results after adsorption are in agreement with previous work using an implicit solvation. 26 Advanced sampling is here of key importance when comparing with recent work by Vilhena et al. 38 where for BSA only minimal structural changes could be observed. Comparing the adsorption of BSA and lysozyme, we find:

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1400 1200

Force (pN)

1000 800 600 400

Orientation 1 Orientation 2

200 0 0

50

100 150 pulling distance (nm)

700

200

Orientation 1 Orientation 2

600 500 Force (pN)

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 14 of 22

400 300 200 100 0 0

5

10

15 20 25 pulling distance (nm)

30

35

40

Figure 6: Force-distance curves for BSA (top) and lysozyme (bottom) during the desorption process. The point of complete desorption is marked by a circle which is defined by zero van-der-Waals interactions between protein and surface.

14

ACS Paragon Plus Environment

Page 15 of 22

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 7: Snapshots of BSA and lysozyme at maximum desorption force in their different starting orientations (‘Or’). Note that the scale has changed between the BSA and lysozyme subfigures.

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

• The secondary structure of BSA (only helices) is largely destroyed, while it survives to some extent in lysozyme. This feature justifies the distinction of “hard” vs. “soft” proteins introduced by Norde and Giacomelli. 19 • Analogously the radius of gyration of BSA increases considerably more for BSA than for lysozyme during adsorption. • As a consequence of the increased spreading of BSA on the surface, the adsorption energy of BSA is considerably (by a factor of around 7) larger than that of lysozyme. This demonstrates that the adsorption energy is not simply proportional to the chain length of the molecules but rather to the surface area over which the protein has spread. The desorption simulations were carried out using steered MD in its classical formulation 41,42 at a velocity of 10 m/s. We find that • BSA becomes almost fully unzipped during pull-off, while several helices survive this process in lysozyme due to its reduced surface interactions and higher internal stability. Again this justifies the notion of hardness of lysozyme as compared to BSA. • The (complete or incomplete) unzippping of the molecule also shows up in the forcedistance diagram of the desorption process. A consecutive series of force peaks accompanies the unzipping process in BSA, while only a single or few (depending on protein orientation) force peaks occur for lysozyme. • The maximum desorption force is larger for BSA than for lysoyzme, but only by a factor of about 2.3. This means that simple estimates of the desorption force based on chain length or adsorption area do not describe our results. Future studies will focus on the desorption process and investigate in particular its dependence on the pull-off velocity. By performing a sequence of simulations with varied velocity, aMD may be used to extrapolate the results to low velocities approaching those used in force spectroscopy experiments. 18 16

ACS Paragon Plus Environment

Page 16 of 22

Page 17 of 22

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

Acknowledgement We acknowledge financial support by Deutsche Forschungsgemeinschaft within project Ur 32/26-1. Furthermore we appreciate the computational resources provided by the compute cluster ’Elwetritsch’ of the University of Kaiserslautern.

References (1) Ratner, B. D.; Bryant, S. J. BIOMATERIALS: Where We Have Been and Where We Are Going. Annu. Rev. Biomed. Eng. 2004, 6, 41–75. (2) More, R. B.; Sines, G.; Ma, L.; Bokros, J. C. Pyrolytic Carbon. Encyclopedia of Biomaterials and Biomedical Engineering 2004, 1308–1319. (3) Winter, M.; Chuinard, C.; Cikes, A.; Pelegri, C.; Bronsard, N.; de Peretti, F. Surgical Management of Elbow Dislocation Associated with Non–Reparable Fractures of the Radial Head. Chirurgie de la Main 2009, 28, 158–167. (4) Raffaini, G.; Ganazzoli, F. Simulation Study of the Interaction of Some Albumin Subdomains with a Flat Graphite Surface. Langmuir 2003, 19, 3403–3412. (5) Agashe, M.; Raut, V.; Stuart, S. J.; Latour, R. A. Molecular Simulation To Characterize the Adsorption Behavior of a Fibrinogen γ–Chain Fragment. Langmuir 2005, 21, 1103– 1117. (6) Yin, J.; Zhao, Y.-P.; Zhu, R.-Z. Molecular Dynamics Simulation of Barnacle Cement. Materials Science and Engineering 2005, 409, 160–166. (7) Horinek, D.; Serr, A.; Geisler, M.; Pirzer, T.; Slotta, U.; Lud, S.; Garrido, J.; Scheibel, T.; Hugel, T.; Netz, R. Peptide Adsorption on a Hydrophobic Surface Results from an Interplay of Solvation, Surface, and Intrapeptide Forces. Proc. Natl. Acad. Sci. USA 2008, 105, 2842–2847.

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

Page 18 of 22

(8) 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. (9) Raffaini, G.; Ganazzoli, F. Protein Adsorption on a Hydrophbic Surface: A Molecular Dynamics Study of Lysozyme on Graphite. Langmuir 2010, 26, 5679–5689. (10) Cohavi, O.; Corni, S.; De Rienzo, F.; Di Felice, R.; Gottschalk, K. E.; Hoefling, M.; Kokh, D.; Molinari, E.; Schreiber, G.; Vaskevich, A. et al. Protein–Surface Interactions: Challenging Experiments and Computations. J. Mol. Recognit. 2010, 23, 259–262. (11) Wei, T.; Carignano, M. A.; Szleifer, I. Lysozyme Adsorption on Polyethylene Surfaces: Why Are Long Simulations Needed? Langmuir 2011, 27, 12074–12081. (12) R¨osch, C.; Huber, C.; M¨ uller, C.; Umanskaya, N.; Hannig, M.; Ziegler, C. Influence of Protein Immobilization on Protein–Protein Interaction Measured by Scanning Force Spectroscopy. physica status solidi (a) 2013, 210, 945–951. (13) M¨ uller, C.;

Wald, J.;

Hoth-Hannig, W.;

Umanskaya, N.;

Scholz, D.;

Han-

nig, C., Matthias M¨ uller; Wald, J.; Hoth-Hannig, W.; Umanskaya, N.; Scholz, D. et al. Protein Adhesion on Dental Surfaces–a Combined Surface Analytical Approach. Analytical and Bioanalytical Chemistry 2011, 400, 679–689. (14) Latour, R. A. Molecular Simulation of Protein–Surface Interactions: Benefits, Problems, Solutions, and Future Directions (Review). Biointerphases 2008, 3, FC2–FC12. (15) Latour, R. A. Perspectives on the Simulation of Protein–Surface Interactions using Empirical Force Field Methods. Colloids and Surfaces B: Biointerfaces 2014, 124, 25– 37. (16) M¨ ucksch, C.; Urbassek, H. M. Enhancing Protein Adsorption Simulations by Using Accelerated Molecular Dynamics. PLOS ONE 2013, 8, e64883.

18

ACS Paragon Plus Environment

Page 19 of 22

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

(17) Sun, Y.; Dominy, B. N.; Latour, R. A. Comparison of Solvation-Effect Methods for the Simulation of Peptide Interactions with a Hydrophobic Surface. Journal of Computational Chemistry 2007, 28, 1883–1892. (18) M¨ ucksch, C.; Urbassek, H. M. Accelerating Steered Molecular Dynamics: Toward Smaller Velocities in Forced Unfolding Simulations. J. Chem. Theory Comput. 2016, 12, 1380–1384. (19) Norde, W.; Giacomelli, C. E. Conformational Changes in Proteins at Interfaces: From Solution to the Interface, and back. Macromol. Symp. 1999, 145, 125–136. (20) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kal´e, L.; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781–1802. (21) MacKerell, J., A. D.; Bashford, D.; Bellott, M.; Dunbrack, J., R. L.; Evanseck, J. D.; Field, M. J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S. et al. All–Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586–3616. (22) 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. (23) Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. 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. (24) 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.

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

(25) Bujacz, A. Structures of Bovine, Equine and Leporine Serum Albumin. Acta Crystallographica Section D 2012, 68, 1278–1289. (26) M¨ ucksch, C.; Urbassek, H. M. Molecular Dynamics Simulation of Free and Forced BSA Adsorption on a Hydrophobic Graphite Surface. Langmuir 2011, 27, 12938–12943. (27) Carter, D.; He, J.; Ruble, J.; Wright, B. PDB–ID: 1AKI. The Structure of the Orthorhombic Form of Hen Egg-White Lysozyme at 1.5 Angstroms Resolution. 1997; http://www.rcsb.org/pdb/explore/explore.do?structureId=1AKI. (28) Humphrey, W.; Dalke, A.; Schulten, K. VMD – Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33–38. (29) Stone, J. An Efficient Library for Parallel Ray Tracing and Animation. M.Sc. thesis, Computer Science Department, University of Missouri–Rolla, 1998. (30) Li, Z.; Wang, Y.; Kozbial, A.; Shenoy, G.; Zhou, F.; McGinley, R.; Ireland, P.; Morganstein, B.; Kunkel, A.; Surwade, S. P. et al. Effect of Airborne Contaminants on the Wettability of Supported Graphene and Graphite. Nature Mater. 2013, 12, 925–931. (31) Kozbial, A.; Li, Z.; Sun, J.; Gong, X.; Zhou, F.; Wang, Y.; Xu, H.; Liu, H.; Li, L. Understanding the Intrinsic Water Wettability of Graphite. Carbon 2014, 74, 218–225. (32) M¨ ucksch, C.; R¨osch, C.; M¨ uller-Renno, C.; Ziegler, C.; Urbassek, H. M. Consequences of Hydrocarbon Contamination for Wettability and Protein Adsorption on Graphite Surfaces. The Journal of Physical Chemistry C 2015, 119, 12496–12501. (33) Werder, T.; Walther, J. H.; Jaffe, R. L.; Halicioglu, T.; Koumoutsakos, P. On the Water–Carbon Interaction for Use in Molecular Dynamics Simulations of Graphite and Carbon Nanotubes. J. Am. Chem. Soc. 2002, 107, 1345–1352. (34) Hamelberg, D.; Mongan, J.; McCammon, J. A. Accelerated Molecular Dynamics: A

20

ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22

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

Promising and Efficient Simulation Method for Biomolecules. J. Chem. Phys. 2004, 120, 11919–11929. (35) de Oliveira, C. A. F.; Hamelberg, D.; McCammon, J. A. On the Application of Accelerated Molecular Dynamics to Liquid Water Simulations. J. Phys. Chem. B 2006, 110, 22695–22701. (36) Wereszczynski, J.; McCammon, J. A. In Computational Drug Discovery and Design; Baron, R., Ed.; Methods in Molecular Biology; Springer Science+Business Media, 2012; Vol. 819; pp 515–524. (37) Wang, Y.; Harrison, C. B.; Schulten, K.; McCammon, J. A. Implementation of Accelerated Molecular Dynamics in NAMD. Comput. Sci. Discov. 2011, 4, 015002. (38) Vilhena, J. G.; Rubio-Pereda, P.; Vellosillo, P.; Serena, P. A.; P´erez, R. Albumin (BSA) Adsorption over Graphene in Aqueous Environment: Influence of Orientation, Adsorption Protocol, and Solvent Treatment. Langmuir 2016, 32, 1742–1755. (39) Zerze, G. H.; Mullen, R. G.; Levine, Z. A.; Shea, J.-E.; Mittal, J. To What Extent Does Surface Hydrophobicity Dictate Peptide Folding and Stability near Surfaces? Langmuir 2015, 31, 12223–12230. (40) Kabsch, W.; Sander, C. Dictionary of Protein Secondary Structure: Pattern Recognition of Hydrogen-Bonded and Geometrical Features. Biopolymers 1983, 22, 2577–2637. (41) Grubm¨ uller, H.; Heymann, B.; Tavan, P. Ligand Binding: Molecular Mechanics Calculation of the Streptavidin–Biotin Rupture Force. Science 1996, 271, 997–999. (42) Izrailev, S.; Stepaniants, S.; Balsera, M.; Oono, Y.; Schulten, K. Molecular Dynamics Study of Unbinding of the Avidin-Biotin Complex. Biophys. J. 1997, 72, 1568–1581.

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

Graphical TOC Entry

22

ACS Paragon Plus Environment

Page 22 of 22