Ranran
Tian
,
Mengbo
Luo
* and
Jingyuan
Li
*
Department of Physics, Zhejiang University, Hangzhou, 310027, China. E-mail: jingyuanli@zju.edu.cn; luomengbo@zju.edu.cn
First published on 27th November 2017
When nanoparticles enter blood or other biological fluids, they tend to contact with a variety of proteins which may hamper their application and even bring adverse impacts. Such nonspecific protein binding is usually weak and proteins reach dynamic equilibrium between adsorption and desorption. Studies on the spontaneous desorption of weakly bound proteins should not only be helpful to enrich our understanding about such nonspecific interactions but also shed light on the strategy to avoid nonspecific binding of proteins. Here, we use molecular dynamics simulations to investigate the interactions of human serum albumin with the self-assembled monolayer (SAM)-coated gold(111) surface, a typical facet of various gold nanoparticles. The response of the protein interfacial behavior to solution pH, especially the spontaneous desorption, is studied. When the solution pH is relatively low, the protein can adsorb to the SAM surface. Such adhesion is attributed to several salt bridges between acidic residues and SAM's protonated groups, and there are water molecules distributed under the adsorbed protein, i.e. interlayer water. Interestingly, the increase of the solution pH reduces the binding affinity of the protein, which engenders the lateral diffusion of the protein and the increase of interlayer water. When the solution pH is further increased, the enhanced lateral diffusion of the protein and the growth of interlayer water disrupt the formation of salt bridges between the protein and the SAM, and the protein progressively dissociates from the SAM. The spontaneous desorption process of the protein and the corresponding mechanisms illustrated in this work may be helpful to develop an antifouling surface and enhance the biosafety of nanomaterials.
Extensive experimental works have investigated the effects of environmental factors on the interactions between proteins and nanoparticles, especially the impacts of solution pH.26–33 In most cases, efforts have been focused on the change of protein adsorption amounts. Proteins tend to aggregate on the nanomaterial surface at pH close to the isoelectric point of proteins, leading to the increase of their adsorption.26–28 Besides, the solution pH also influences the zeta potential of nanomaterials, which affects the interaction with proteins and the resulting protein adsorption. For example, at acidic pH (e.g., pH = 1) the cerium oxide nanoparticles adsorb much BSA because of electrostatic attraction with proteins, the adsorption decreases as pH increases.29 In short, protein adsorption can be regulated by the solution pH. However, detailed information about the response of the interactions between proteins and nanoparticles, such as the interaction mode, kinetic processes of protein adsorption and desorption as well as the underlying mechanism, still remains elusive.
Molecular dynamics (MD) simulations have been widely used to study the interactions between proteins and nanomaterials, such as graphene,16,34–36 carbon nanotubes,13,37 gold nanorods,15 gold nanospheres38 and monolayer-protected gold nanoparticles.39 Useful information about the protein adsorption, such as the adsorption process, binding stability and the influencing factors (e.g. the size and surface properties of nanoparticles), has been illustrated in these studies. A comprehensive understanding about the interactions between proteins and nanoparticles requires full description of the protein interfacial behavior. However, there is still limited work about other aspects of the interfacial behavior of proteins, e.g., the evolution of the adsorbed protein and the process of spontaneous desorption.
In the present study, we use MD simulations to investigate the protein adsorption, as well as the desorption of weakly bound proteins on the self-assembled monolayer (SAM)-coated gold(111) surface, a typical facet of various gold nanomaterials, including nanorods and prisms.15,40–43 The most abundant plasma protein, human serum albumin (HSA), is selected as the model protein, and the interfacial behavior of the protein in response to the change in the solution pH is systematically investigated. When the solution pH is relatively low, the protein can stably bind to the SAM surface. Such binding is largely attributed to several salt bridges formed by acidic residues with SAM's protonated groups, and there are water molecules distributed under the adsorbed protein (i.e. interlayer water). As the solution pH increases and the protonation degree of the SAM reduces, the adsorption stability of the protein dramatically decreases. The protein undergoes substantial lateral diffusion and the number of interlayer water increases. When the solution pH is further increased, the enhanced lateral diffusion of the protein, the growth of interlayer water and the sparse population of protonated groups disrupt the formation of salt bridges between the protein and the SAM, and the protein progressively dissociates from the SAM surface. In short, our results illustrate the kinetic process and the underlying mechanism of protein desorption which may be helpful for the development of an antifouling surface by enhancing protein desorption.
The structure of human serum albumin (HSA) is obtained from the Protein Data Bank (PDB entry 1N5U).56 The corresponding protonation states of HSA residues at different pHs are determined by PDB2PQR 2.0.0 (http://nbcr-222.ucsd.edu/pdb2pqr_2.0.0).57 HSA is a heart-shaped protein composed of 583 residues with three domains (named I, II and III). A previous study has shown that HSA tends to adopt side-on orientations when it adsorbs on the hydrophilic surface.58 Therefore, there are two possible binding surfaces that are denoted by plane A and plane B (Fig. S1, ESI†).
First, SAMs with 50% and 25% protonation degree are prepared respectively (Fig. 1b and Fig. S2, ESI†). The SAMs and HSA are then equilibrated separately. Subsequently, HSA is placed at ∼5 Å above the SAMs with side-on orientations. The adsorption processes with both plane A and plane B orientating toward the SAMs are studied (Fig. 1c and d). Chlorine and sodium ions are added to maintain the system as neutral with an ionic strength of 0.2 M. After initial steepest descent minimization, 10 ns simulation is performed with all heavy atoms of both HSA and SAMs restrained to equilibrate the water molecules and ions. Next, 80 ns simulation is conducted to study the adsorption behavior of HSA on the 50% (i.e. adsorption-50%) and 25% (i.e. adsorption-25%) protonated SAMs.
To study the response of the adsorbed protein to the changes in the solution pH, we construct the other two systems on the basis of the final configuration of the system adsorption-50%, wherein the protonation degree of SAMs is switched to 25% (i.e. adjustment-25%) and 6.3% (i.e. adjustment-6.3%), respectively.59 The behaviors of adsorbed HSA in these two systems are investigated by 30 ns simulation.
All MD simulations are performed in the NPT ensemble (1 atm and 300 K). The pressure is controlled using a semi-isotropic Parrinello–Rahman barostat,60,61 and only the fluctuation in the z-axis is allowed. The temperature is maintained using a Nose–Hoover thermostat62 with a relaxation time of 0.5 ps. Parameters of SAMs are taken from the CHARMM General Force Field (CGENFF);63 and the CHARMM36 Force Field64 is adopted to describe protein molecules, sodium and chloride ions. The TIP3P model65 which is consistent with the chosen force fields is employed for water. The particle-mesh Ewald (PME) method66 is employed to calculate long-range electrostatic interactions, and a typical cutoff distance of 12 Å is applied to compute short-range electrostatic and van der Waals energies. Standard periodic boundary conditions are applied to all simulations, and the use of the LINCS algorithm allows a time step of 2 fs.67 All simulations are carried out using the Gromacs 4.6.5 package,68 and snapshots are rendered by the visual molecular dynamics (VMD) program.69
Fig. 2 shows that all three domains contribute to the adsorption of HSA with plane A, while only a fraction of the surface region directly contacts with the SAM (i.e. binding regions, Fig. 2b). And about 92.8% of the salt bridges are attributed to nine acidic residues, i.e. E167 and D173 (domain I), E277, E280, E297 and D301 (domain II), and E396, E442 and D563 (domain III), in these binding regions. In other words, the protein is essentially anchored by several salt bridges formed by these key residues. And there are water molecules distributed under the binding regions of the protein (denoted as interlayer water) after the adsorption is accomplished (Fig. 3a). Even though the protein adsorption does not fully exclude the hydration layer of SAMs, the approaching protein still needs to displace the hydration water. We compute the number of interlayer water (i.e. water molecules under the protein binding regions), and the number considerably decreases from ∼500 to ∼140 (Fig. 3b). Similar to many cases of the protein adsorbed onto hydrophilic surface,70–73 the interfacial water mediates the adsorption process of HSA onto the 10-amino-1-decanethiol (–S(CH2)10NH2)-protected gold nanoparticles.
The response of adsorbed HSA to the change of the solution pH is then investigated. As the solution pH increases from 10.2 to 10.7, the protonation degree decreases to 25% correspondingly. HSA remains bound to the SAM (Fig. 4a, adjustment-25%), while the interaction between the protein and SAM is substantially weakened. More interestingly, the protein undergoes considerable lateral diffusion on the SAM surface. The self-diffusion coefficient of HSA in both adsorption-50% and adjustment-25% is then determined by measuring its in-plane mean square displacement (based on the last 10 ns trajectory, Fig. 4b). The calculated self-diffusion coefficient of HSA in adjustment-25% is 2.6 × 10−7 cm2 s−1, which is comparable to the other weakly adsorbed proteins (e.g., Aβ trimer on the COOH-SAM, 2.3 × 10−7 cm2 s−1)74 and one order of magnitude larger than the self-diffusion coefficient of the protein in adsorption-50% (2.8 × 10−8 cm2 s−1). Such enhanced lateral diffusion of HSA should reciprocally contribute to the reduction of adsorption stability. On top of it, there are more water molecules distributed under the binding regions of the protein. The average number of interlayer water molecules increases to ∼370 (based on the last 10 ns trajectory, Fig. S4a, ESI†), much higher than the case of adsorption-50% (∼150). It should be noted that some water molecules are in the vicinity of salt bridges formed by acidic residues with protonated amine groups. Some of these water molecules form hydrogen bonds with the salt-bridge groups (Fig. S5, ESI†), and the increase of interlayer water should also disturb the interactions between HSA and the SAM. The average number of salt bridges (i.e. the last 10 ns of the simulation, Fig. S4b, ESI†) is only 1.9, which is one-sixth of the corresponding number in adsorption-50%. That is, the extent of salt bridge reduction is more dramatic than the reduction of the SAM's protonation degree. Besides, only two domains (domain II and III) participate in the protein adsorption, and the salt-bridge interactions are mainly attributed to three residues: D301 of domain II, E396 and D563 of domain III. The distance between the mass center of HSA and the SAM surface is about 3.1 nm, which is considerably higher than that in adsorption-50% (2.5 nm, see Fig. S6, ESI†). Taken together, the increase of the solution pH and the concomitantly decreased protonation degree of the SAM considerably affect the interfacial behavior of the adsorbed protein.
To further characterize the behaviors of HSA on the 25%-protonated SAM surface, we also assess its adsorption behavior (i.e. adsorption-25%). The protein can approach and bind to the SAM. The adsorbed protein similarly undergoes substantial lateral diffusion on the SAM surface, and the corresponding diffusion coefficient is 4.1 × 10−7 cm2 s−1 (based on the last 10 ns trajectory, Fig. S7, ESI†), comparable to the value of adjustment-25%. The average number of salt bridges (Fig. S8, ESI†) is 2.0, which is also close to the case of adjustment-25%. Besides, there are also two domains (domains II and III) involved in the adsorption of HSA (Fig. S9, ESI†), and water molecules are similarly distributed in the region between the protein and SAM. Moreover, the distance between the mass center of HSA and the SAM surface (∼3.0 nm, Fig. S10, ESI†) is still similar to the case of adjustment-25%. In short, the interfacial behavior of the protein in adsorption-25% converges with the adjusted behavior of the pre-adsorbed protein in adjustment-25% and is distinct to the case of adsorption-50%. Such an interfacial behavior can thus be denoted as the weakly-bound state. In addition, our results suggest that HSA can still adsorb on the surface as long as there are two salt bridges.
In addition, we examine the response of the adsorbed protein when the protonation degree of the SAM is further decreased to 6.3% (i.e. adjustment-6.3%). In this case, the adsorption stability of HSA continuously decreases and the protein essentially desorbs from the surface after t = 21 ns (Fig. 5): HSA no longer forms salt bridges with the SAM, and the average number of heavy atoms in contact with the SAM (within 0.5 nm from nitrogen atoms, the uppermost heavy atoms, of the SAM) becomes merely 0.4. Moreover, the protein progressively moves away from the SAM: the mass center distance, i.e. the distance between the mass center of HSA and the SAM surface, increases from 2.4 nm to 4.1 nm.
Even in the early stage of desorption (the first 4.5 ns), the protein undergoes considerable lateral diffusion with a diffusion coefficient of 1.4 × 10−7 cm2 s−1. Such diffusion indeed disrupts the interaction with the SAM. The number of salt bridge decreases to 0, and the heavy atom contact number decreases dramatically to about 10. Meanwhile, there are more water molecules distributed in the region between the protein and the SAM. The number of interlayer water increases from ∼140 to ∼240 at t = 4.5 ns (Fig. S11, ESI†). As mentioned above, the growth of interlayer water disturbs the interaction between the protein and SAM. Especially, the surrounding water molecules facilitate the breakage and prohibit the re-formation of salt-bridge interactions. For example, among the water molecules in the vicinity of the salt bridge between E277 and the protonated amine group, there are at least five water molecules that form hydrogen bonds with the salt-bridge groups (Fig. 6a). When the salt bridge breaks due to the lateral diffusion of the protein or the fluctuation of residue side chains, one water molecule nimbly inserts into the broken salt bridges, forming a water bridge between E227 and the amine group (Fig. 6b). The salt bridge is then replaced by the water bridge. In addition, the invasion of water molecules also enhances the hydration of salt-bridge groups, and the re-formation of the salt bridge is thus suppressed by the dehydration penalty. Similar phenomena are found during the breakage of salt bridges of E396 and D563 (Fig. 6c, d and Fig. S12, ESI†). In short, the invasion of interlayer water along with the lateral diffusion of protein disturbs the interactions between the protein and SAM.
After t = 4.5 ns, the lateral diffusion of the protein is further enhanced and the diffusion coefficient increases up to 5.3 × 10−7 cm2 s−1 (Fig. S13, ESI†), even larger than the case of the weakly-bound protein. Due to the limited number of protonated amine groups, the probability of acidic residues moving to the vicinity of these groups is very low. Besides, relatively rapid lateral diffusion of the protein further prevents the formation of the salt bridge with the SAM even when the acidic residues are close to the protonated amine groups. In addition, along with the dissociation of the protein from the SAM, the number of interlayer water further increases to ∼480 (Fig. S11, ESI†). An extended hydration layer forms in-between the protein and the SAM, and both protein's acidic residues and SAM's amine groups are well solvated. Therefore, the re-formation of salt bridges with the SAM needs to remove multiple water molecules and surmount considerable desolvation penalty. In short, the enhanced lateral diffusion of the protein, the existence of interlayer water and the sparse population of protonated amine groups substantially prevent the protein from forming the stable salt bridge with the SAM. After t = 21 ns, there are no salt-bridge interactions with the SAM. The protein diffuses away from the SAM surface even though there is still weak attraction with the surface.
During the desorption process, three domains sequentially depart from the SAM surface. The heavy atom contact number of each domain is calculated separately (Fig. 7). At t = 7 ns, the contact number of domain I decreases to 0, and this domain is then considered to be desorbed. The following dissociation of domain III occurs at t = 11 ns. In addition, domain II finally desorbs from the SAM surface at t = 21 ns, leading to the departure of the whole protein. It should be noted that the dissociation of each domain is triggered by the elimination of the salt bridges with the SAM. For example, domain III no longer forms the salt bridge with the SAM after t = 4.5 ns (Fig. S14, ESI†), the contact number of domain III progressively decreases to 0, and the domain desorbs at t = 11 ns. Similarly, there is no longer the salt bridge between domain II and the SAM after t = 21 ns (Fig. S14, ESI†). As the number of salt bridges decreases to ∼0, it is hard for the domain to contact with the SAM and thus the protein diffuses away from the surface. Hence, the breakage of salt bridges plays a critical role for the desorption of each domain.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7cp05515c |
This journal is © the Owner Societies 2018 |