Molecular Understanding of Laccase Adsorption on Charged Self

School of Materials and Chemical Engineering, Hainan University, Haikou 570228,. China, P. R. China. *. Corresponding author. Tel./fax: +86 20 8711406...
0 downloads 0 Views 3MB Size
Article Cite This: J. Phys. Chem. B 2017, 121, 10610-10617

pubs.acs.org/JPCB

Molecular Understanding of Laccase Adsorption on Charged SelfAssembled Monolayers Jie Liu,†,‡ Yun Xie,§ Chunwang Peng,† Gaobo Yu,†,∥ and Jian Zhou*,† †

School of Chemistry and Chemical Engineering, Guangdong Provincial Key Lab for Green Chemical Product Technology, South China University of Technology, Guangzhou 510640, P. R. China ‡ Key Laboratory for Green Chemical Process of Ministry of Education, School of Chemical Engineering and Pharmacy, Wuhan Institute of Technology, Wuhan 430073, P. R. China § Huizhou University, Huizhou 516007, P. R. China ∥ School of Materials and Chemical Engineering, Hainan University, Haikou 570228, P. R. China S Supporting Information *

ABSTRACT: Controlling the orientation of laccase on electrodes is crucial for the achievement of fast direct electron transfer. It is important to find a short pathway between the T1 copper site of laccase and a substrate during the laccase immobilization. In this work, we studied the adsorption orientation and conformation of Trametes versicolor laccase (TvL) on two kinds of charged selfassembled monolayers (SAMs), including NH2−SAM and COOH−SAM, by parallel tempering Monte Carlo and all-atom molecular dynamics simulations. TvL adsorbs on positively and negatively charged surface with “end-on” and “lying” orientation, respectively. On the positively charged surface, T1 copper site of TvL is closer to the surface. The orientation of TvL on positively charged surface is narrower than that on negatively charged surface. Thus, the positively charged surface is more conducive to the immobilization of TvL. The conformational changes of TvL on the charged surfaces are analyzed by RMSD, superimposed structures, dipole moment, gyration radius, and eccentricity. Results show that native structures of TvL are well preserved when it adsorbs on the charged surfaces. This work provides atomistic insight into the mechanism of TvL adsorption on charged surface and is helpful for the design and development of laccase-based electrodes.

1. INTRODUCTION

An efficient enzymatic electrode can be established by developing efficient electron transfer between the electrode and catalytic sites of the enzyme. To achieve the fast electron transfer, two schemes, including mediated electron transfer (MET)11 and direct electron transfer (DET),12 are usually adopted by many researchers. The difference between MET and DET is whether it needs redox mediators or not. Redox mediators usually make the design of enzyme electrode more complicated, which may induce higher overpotentials and reduce operational stability during the electrocatalytic process. Redox mediators usually are toxic.13 Compared with MET, DET avoids the use of redox mediators. To achieve DET, it is necessary to orient the redox center of enzyme toward the electrode. For laccase, it is crucial to shorten the distance between the T1 copper site of laccase and the substrate. Many researchers focused on the immobilization of laccases on different kinds of solid materials, such as carbon nanotubes,14−19 carbon electrode,20 titanium dioxide,21 nanofibrous membranes,8,22 mesopurous silica materials,23−26 gold electrodes/nanoparticles,13,27,28 and self-assembled monolayers

Laccase (EC 1.10.3.2), which is subordinate to a subclass of multicopper oxidases, is popular for its ability of catalyzing the oxidation of various aromatic substrates, and oxygen is usually required as a second substrate.1 Four copper ions within the laccase, classified as type-1 (T1), type-2 (T2), and type-3 (T3), are catalytic sites.2 Either T1 or T2 copper site contains a single copper ion while the T3 copper site consists of two copper ions. T1 copper site is the primary oxidation site which accepts the electron from the electrode. Then the electron is transferred to the T2/T3 trinuclear copper sites through an intermolecular electron-transfer mechanism and induces the reduction reaction of oxygen.3,4 With a high redox potential of the T1 copper site and catalytic efficiency of oxidizing pollutants, laccase has been widely used in biofuel cells, biosensors, and pollutant degradation.5−7 In practical applications, the reaction conditions (e.g., organic solution and pH) can induce the denaturation of free laccases, especially when mixed with reactants. It is also hard to recycle the free laccases from mixed reactants, which limits the industrialization of laccases. Through immobilization of laccase, the thermal stability, reusability, and suitable range of working pH as well as temperature are all enhanced.8−10 © 2017 American Chemical Society

Received: September 2, 2017 Revised: October 26, 2017 Published: November 2, 2017 10610

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B (SAMs).29−32 Arzola et al.30 investigated the adsorption of laccase on bare, hexadecanethiol and highly ordered pyrolytic graphite (HOPG) modified gold electrodes, and both electrochemical response and catalytic activity can be retained by SAM modified electrodes. Gutiérrez-Sánchez et al.13 covalently combined gold nanoparticles with laccases and assembled these compounds into porous graphite electrodes. They found that the DET reaction can be enhanced by controlling the orientation of laccase on functionalized gold nanoparticles, and the ordered adsorption of laccase on electrode with T1 copper site oriented close to the electrode was significant for fast DET occurring. Most researchers focused on the laccase loading capacity of the nanomaterials and the optimization of reaction conditions (e.g., pH, temperature, solvent, et al.).8,15,24,33−36 However, the charged property of surface which can affect the ordered orientation of adsorbed enzymes also plays a key role in the enzyme immobilization.13,37−39 Olejnik et al.4 found that the electrocatalytic current value of laccase adsorbed on positively charged surface was stronger than that on negatively charged surface. They mentioned that the orientation of laccase on electrode surfaces played a dominant role in the electrocatalytic current. Up to now, atomistic-level information about the orientation of laccase on charged surfaces has rarely been explored. It is difficult for experiments to directly detect the orientation of a single enzyme on a substrate. Molecular simulation provides an efficient avenue to investigate molecular level interactions between proteins and solid surfaces.40−46 The preferred orientation of proteins on surfaces could hardly be caught by single conventional molecular dynamics (MD) simulation. Thus, it is necessary to perform a Monte Carlo (MC) simulation first, which can catch the optimized configuration of proteins adsorbed on surfaces quickly and save computational resources, before long-time MD simulations.47 In our previous work, we developed a parallel tempering Monte Carlo (PTMC) simulation, which can fast and efficiently predict the preferred orientation of a protein on a charged surface48 and provide the initial configurations for the later long-time MD simulations.49−51 To the best of our knowledge, the accurate orientation of laccase on charged surfaces and its effect on DET between laccase and electrode have not been studied yet. In this work, the combination of PTMC and all-atom molecular dynamics (AAMD) simulations were performed to study the orientation and conformational changes of laccase on charged SAMs. The orientation of T1 copper site within the laccase was analyzed in detail. Electrocatalytic current differences of laccase adsorbed on charged surfaces, which were found in experiment,4 will be illustrated theoretically. This work can provide a better understanding of laccase immobilization on electrodes.

Figure 1. Crystal structure of TvL and its three domains. D1, D2, and D3 represent domain 1, domain 2, and domain 3, respectively. Green beads represent the copper clusters (i.e., T1 and T2/T3 copper sites).

three domains in TvL are shown in Figure 1. In both PTMC and AAMD simulations, the relative positions of copper atoms within TvL are kept as its crystal positions. As TvL is in its oxidized state, the coppers in TvL are assigned with the charge of +2 e. The lysine and arginine residues along with the Nterminal of TvL were protonated, whereas the glutamic acid and aspartic acid along with the C-terminal were deprotonated. The protonation state of histidine was calculated by the PDB2PQR and APBS algorithms52,53 with CHARMM force field. The net charge of TvL is −4 e. 2.1. PTMC. The PTMC48 simulations were performed to get the preliminary orientation of TvL on charged surfaces. In order to keep the basic structure of protein, united-residue model41 (i.e., each residue is reduced to a sphere centered at its α-carbon) was used. During PTMC simulations, the protein backbone was maintained rigid, and SAM surface was treated as a flat surface. TvL can freely rotate and translate around its center of mass during the whole PTMC simulation. To ensure enough sampling of PTMC algorithm, we plotted the interaction energies between TvL and charged surfaces over electric dipole orientation and hydrophobic dipole orientation in Figure S1. In our previous work,54 we proposed that the orientation was controlled by both electric and hydrophobic dipoles. Sometimes the electric dipole played a dominant role, and sometimes the other. As can be seen in Figure S1, for TvLs adsorption on both charged surfaces, electric dipole dominates the adsorption. Thus, in section 3, electric dipole is discussed in detail. The force field parameters were taken from our previous works.48 Both van der Waals (vdW, Uvdw) and electrostatic interactions (Uele) were considered. In total, 40 000 000 MC cycles were performed; the former 20 000 000 MC cycles were used for system equilibration while the remaining cycles were used for result analysis. We adopted five replicas, at 310, 350, 400, 500, and 700 K, to ensure sufficient energy overlap between neighboring replicas and to allow for the acceptance of configuration swaps. Each swap was carried out every 500 cycles. The final results were taken as the initial orientation for further AAMD simulations.49,51,55 2.2. AAMD. AAMD simulations were performed by GROMACS 4.5.4 with CHARMM force field.56 The carboxyl(COOH-) and amino- (NH2-) terminated SAMs on Au(111) with √3 × √3 structures are simulated as the model charged surfaces. There are totally 224 thiol chains and 2016 gold

2. SIMULATION DETAILS The high-resolution crystal structure of Trametes versicolor laccase (TvL) which is in its oxidized, copper-complete state2 (PDB-ID: 1GYC) was obtained from the RCSB (www.rcsb. org). TvL structure is a monomer and organized in three sequentially arranged domains (i.e., domain 1, domain 2, and domain 3, see Figure 1). The T1 copper site is located in domain 3. The trinuclear copper cluster, constructed by T2 and T3 copper sites, is embedded between domains 1 and 3, and both domains provide residues for the coordination with copper atoms.2 The secondary structures, copper clusters, and 10611

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B atoms. Fifteen thiol chains were deprotonated/protonated, representing COOH−SAM and NH2−SAM with the surface charge densities of −0.05 and +0.05 C·m−2, respectively. The surfaces have dimensions of 6.99 × 6.92 nm2. During AAMD simulations, the sulfur atoms of thiol chains along with gold atoms were fixed. The potential parameters for SAMs were from the CHARMM force field for lipids.56 Hydrogen atoms within TvL were added by GROMACS 4.5.4.57 TvL was initially put above the surface with the closest distance of 0.5 nm. Water molecules, which were described by TIP3P water,58 were added into the single protein box with the dimension of 6.99 × 6.92 × 6.00 nm3. Water molecules within 0.28 nm of TvL and SAMs were removed to avoid overlaps of coordinates.47 Counter ions, including Na+ and Cl−, were added to keep the systems neutral. TvL in the bulk solution was also simulated as a reference system. Each system was initially optimized by the steepest descent method. During this step, TvL was kept rigid as its crystal structure. Then, NVT-MD simulations were performed at 310 K, and the time step was 2 fs. The Nose-Hoover thermostat59,60 with a time constant of 0.5 ps was adopted to control the temperature. Bonds containing hydrogen atoms were constrained by the LINCS algorithm.61 A cutoff distance of 1.0 nm was used to calculate the LJ interactions, and the 3dc particle mesh Ewald method62 was adopted to evaluate the electrostatic interactions. For systems including flat surfaces, only x and y directions were considered in periodic boundary conditions, and two hard walls were added at both top and bottom of the box.63 In order to study the early adsorption of TvL on charged surfaces, a 100 ns AAMD simulation was performed. To confirm the convergence of simulations, the evolution of rootmean-square deviation (RMSD) versus simulation time is plotted in Figure S2. It can be found that RMSD reached equilibrium in all systems after 60 ns. Thus, 100 ns simulation was sufficient for the simulation of TvL adsorption. TvL in the bulk solution was simulated in the same simulation parameters as mentioned above, except for the periodic boundary conditions. As there is no flat surface in bulk system, the periodic boundary conditions with three-dimensional directions were used. The visual molecular dynamics (VMD) program64 was used for structure visualization.

Figure 2. Adsorption configuration of TvL on charged SAMs simulated by PTMC (shown in parts (a) and (b)) and AAMD (shown in parts (c) and (d)). The blue and red planes in (a) and (b) are presented for positively and negatively charged surface, respectively. The surfaces terminated with blue groups and red groups in (c) and (d) refer to NH2−SAM and COOH−SAM, respectively. The orange arrow indicates the direction of the dipole. Copper sites are represented in green beads. The value of d indicates the distance between T1 copper site and the SAM surface.

3.1. Orientation of T1 Copper Site. The adsorption configurations of TvL on both charged SAMs (i.e., COOH− SAM and NH2−SAM), simulated by PTMC and AAMD, are shown in Figure 2. TvL adsorbs on NH2−SAM with the “endon” orientation (i.e., the principal axis of protein relatively parallel to the surface normal), whereas on COOH−SAM, the “lying” orientation (i.e., the principal axis of protein is relatively perpendicular to the surface normal) is found. The orientation trends obtained from AAMD are consistent with those from PTMC, which indicates that the PTMC algorithm is suitable for predicting orientation of proteins. From Figure 2, it can be seen that the dipole direction of TvL is almost parallel to the surface normal when it adsorbs on NH2−SAM. Because TvL has a net charge of −4 e, the orientation of TvL on NH2−SAM is induced by the electric dipole, which is consistent with what we had found in our previous works.37,38,49,51,54,55,65,66 For TvL adsorption on negatively charged surface, as enzyme has the same charge sign with the surface, the orientation is controlled by the distribution of charged patches above the protein surface. The key factor which should be emphasized in DET-based laccase electrode is an ordered orientation of laccases on electrode with the T1 copper site close to the surface.4,13,17 The distances between T1 copper site and top atoms of surfaces are noted in Figure 2. It is obvious that when TvL adsorbs on NH2−SAM, the distance is shorter. Thus, controlling the orientation of laccase by a positively charged surface is more conducive for the DET between laccase and the electrode. Olejnik et al.4 found that the electrocatalytic current value of laccase on a negatively charged surface was significantly smaller than that on a positively charged surface, although the amount of adsorbed enzyme was similar. They deduced that it may be caused by the distance between the electrode and the T1 copper site.4 In this simulation study, we present more detailed orientation information on laccase on charged surfaces. It is consistent with experiment results that the positively charged surface is favorable for the DET-based laccase electrodes. We

3. RESULTS AND DISCUSSION In this work, the PTMC and AAMD methods were combined to investigate the orientations of TvL on charged surfaces. The preliminary orientations of TvL on charged surfaces were determined by PTMC. The distance between T1 copper site and charged surfaces, which is the key factor for DET occurring, was analyzed. Furthermore, the RMSD, superimposed structures, gyration radius, eccentricity, and dipole moment of TvL were calculated to study the conformational changes of TvL during the adsorption. Results were obtained by averaging over the last 40 ns AAMD trajectories. Simulation results are summarized in Tables 1−4 and Figures 2−5. Table 1. Interaction Energies of TvL Adsorbed on Charged SAMs by AAMD NH2−SAM COOH−SAM

Utot (kJ·mol−1)

Uele (kJ·mol−1)

Uvdw (kJ·mol−1)

−466.7 ± 51.8 −297.0 ± 64.3

−307.9 ± 34.8 −219.4 ± 64.5

−158.8 ± 19.4 −77.5 ± 18.0 10612

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B

That means a positively charged surface can promote the stability and conductivity of immobilized TvL. 3.3. Binding Site. To better understand the interactions between TvL and charged surfaces, charged residues of TvL in contact with the SAMs (i.e., binding sites) in the final configuration of AAMD simulations were picked out and are shown in Figure 4 and Table 2. Residues containing atoms

deduce that the weak electrocatalytic current value of laccase on the negatively charged surface is mainly caused by the orientation of T1 copper site, which is almost oriented toward the solution (shown in Figure 2). 3.2. Orientation Distribution and Interaction Energy. Orientation distribution and the interaction energy between proteins and surfaces can be used to quantitatively describe the stability of proteins adsorbed on surfaces. The interaction energies are averaged over the last 40 ns of the MD simulation by block averaging (as shown in Table 1). The adsorption behaviors of TvL on both charged surfaces are induced by the synergy of electrostatic and vdW interactions, while the electrostatic interactions play a more dominant role in the adsorption process. In our previous works, we also concluded that the electrostatic interactions usually dominated the adsorption behavior of charged proteins on charged surfaces.49,51,54,55,65 Due to the stronger electrostatic interaction energy between laccase and NH2−SAM, more side chains move close to the surface, and the vdW interaction energy enhances accordingly. The orientation distributions of TvL on both charged surfaces are shown in Figure 3. The orientations of TvL

Figure 3. Orientation distributions of TvL adsorbed on oppositely charged surfaces by PTMC and AAMD.

calculated by PTMC are consistent with the global-minimumenergy areas in Figure S1, which means the preferred orientations are found by a well sampling. Compared with the orientation distributions of TvL on both charged surfaces predicted by PTMC, the results from AAMD present slight shifts. In PTMC simulations, we ignore the side chains of protein residues and only consider the locations of α-C. Therefore, the dipole direction is slightly deviated. On the other hand, the protein is kept rigid in PTMC, whereas the flexible side chains are taken into account in AAMD simulations, which also contribute to the slight shift of orientation distributions. Moreover, as PTMC is based on the rigid united-residue protein model, the orientation distribution is narrower than that by AAMD. Nevertheless, it is clear that the trends of the orientation distributions of laccase adsorbed on both charged surfaces are basically the same. It can also be found from Figure 2 that the final orientation configurations of TvL on charged surfaces from PTMC and AAMD are quite similar. As can be seen from Figure 3, the orientation distribution of TvL on NH2−SAM is narrower than that on COOH−SAM. The interaction energy between TvL and NH2−SAM is also stronger. It is mainly caused by the net charge of TvL. As TvL carries net negative charges, it tends to adsorb on NH2−SAM.

Figure 4. Final configurations of TvL on NH2−SAM (a, c, e) and COOH−SAM (b, d, f) by AAMD. (a, b) are electrostatic maps of TvL. (c, d) and (e, f) are the side-view and top-view configurations, respectively. Charged residues less than 3.5 Å from the surface are represented in VDW mode. Positively and negatively charged residues are rendered in blue and red, respectively.

within 0.35 nm of the surface are considered as contacting residues.54 As can be seen from Figure 4 and Table 2, four negatively charged residues of TvL, including GLU496, ASP492, ASP364, and GLN499, strongly interact with the Table 2. Key Binding Sites of TvL on Oppositely Charged Surfaces negatively charged residues TvL on NH2− SAM TvL on COOH− SAM 10613

GLU496, ASP492, ASP364, GLN499 ASP42, ASP138, GLU142

positively charged residues NA LYS482, ARG43, LYS39, LYS40, LYS130, LYS194, ALA1

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B NH2−SAM surface while no positively charged residue appears near the surface. Meanwhile, for TvL on COOH−SAM, there are three negatively charged residues (i.e., ASP42, ASP138, and GLU142) and seven positively charged residues (i.e., LYS482, ARG43, LYS39, LYS40, LYS130, LYS194, and ALA1) directly adsorbed on the surface. Among these charged residues, the GLN499 and ALA1 are C-terminal and N-terminal residues. Interactions between TvL and charged SAM can be more visible by examining the distribution of charged residues above the TvL surface by electrostatic maps. From Figures 4a and b, we can see that the negatively charged region of TvL is close to the NH2−SAM while the positively charged region of TvL is close to the COOH−SAM. The number of contact residues of TvL on NH2−SAM is less than that on COOH−SAM. In section 3.1, we mentioned that TvL adsorbed on NH2−SAM with the “end-on” orientation, while it adsorbed on COOH− SAM with the “lying” orientation. In other words, TvL adsorbed on negatively charged surface with more contact residues and had more flat contact. However, as discussed above, we indicated that the adsorption of TvL on NH2−SAM was more stable than that on COOH−SAM. The main reason is that, on NH2−SAM, the dipole direction is parallel to the surface normal and the net charge of TvL is −4 e. Thus, TvL can stably adsorb on NH2−SAM with a narrower orientation distribution and stronger interaction energy. On COOH−SAM, the dipole direction is not antiparallel to the surface normal, as mentioned above. From Figure 4, it can be found that the adsorption of TvL on COOH−SAM is induced by the cooperation of electrostatic interactions and the geometric feature. In our previous work, we also found that the geometric feature and electric dipole of a protein play an important role in determination of orientation distribution.55 3.4. Conformational Changes. The conformational changes of TvL adsorbed on charged surfaces were analyzed by RMSD, superimposed structures, dipole moment, gyration radius, and eccentricity. 3.4.1. RMSD and Superimposed Structures. RMSD and superimposed structures are used to qualitatively analyze the conformational changes of TvL adsorbed on charged surfaces compared with that in bulk solution and crystal state. RMSD values calculated in this work are based on structures of simulated and crystal TvL; only backbone structures are considered as well. The RMSDs of the whole protein and three domains are calculated and shown in Table 3. The three

slightly smaller than that in bulk. In some previous simulation works, researchers also found the same phenomenon that the charged surface can preserve the native conformation of proteins.55,67,68 This is mainly due to the hydrophilic charged surface which limits the free movement of protein’s side chain. For three domains, the conformational change of D3 is relatively larger than those of the other two domains, whether TvL is adsorbed on charged surfaces or dissolved in bulk. Piontek et al.2 mentioned that the D3 domain had more helical contents than the other two domains. It is known that the hydrogen bonds in helix, coil, and turn structures are relatively less than those in sheet structures.55 Therefore, the secondary structure of D3 is more flexible than those of D1 and D2. Yet, this does not mean that TvL will lose its native conformation. The conformation of a protein in the bulk solution can represent the native structures.47,55 RMSDs of TvL on charged surfaces are similar to that in bulk. Thus, the native structure of TvL is well preserved when it is adsorbed on the studied charged surfaces. Furthermore, to provide a visual assessment of TvL structural changes after AAMD, we analyzed the superimposed structures of three domains of the adsorbed TvL on its crystal structure (shown in Figure 5). It reveals that the structural deviation

Figure 5. Simulated structures (red) of three domains (D1, D2, and D3) of TvL superimposed on its crystal structure (cyan).

mainly occurs in domain 3. Structural changes take place in the coil, turn, and helix structures. From Figure 5, it can also be seen that most secondary structures are preserved, which means the changes of simulated structure, with respect to the crystal structure, do not affect the biological activity of TvL. 3.4.2. Dipole Moment, Gyration Radius, and Eccentricity. The conformational changes of TvL adsorbed on charged surfaces are also quantitatively analyzed by dipole moment, gyration radius, and eccentricity. The dipole moment is usually used to represent the distribution of protein charges. From Table 4, it can be seen that the dipole moments of TvL on charged surfaces, 681 D for COOH−SAM and 639 D for NH2−SAM, are relatively larger than that of TvL in bulk (623D). This is mainly caused by the charged surface which can

Table 3. RMSD of TvL and Its Three Domains in the Bulk Solution and Adsorbed on SAM Surfaces by AAMD

in bulk on COOH− SAM on NH2− SAM

RMSD_total (nm)

RMSD_D1 (nm)

RMSD_D2 (nm)

RMSD_D3 (nm)

0.45 ± 0.02 0.40 ± 0.02

0.11 ± 0.01 0.12 ± 0.01

0.13 ± 0.02 0.11 ± 0.01

0.41 ± 0.02 0.32 ± 0.01

0.43 ± 0.04

0.18 ± 0.01

0.12 ± 0.01

0.27 ± 0.02

Table 4. Averaged Properties of TvL in Bulk Solution and Adsorbed on Charged SAMs by AAMD

domains (i.e., domain 1, domain 2, and domain 3) are named as D1, D2, and D3, respectively. When TvL adsorbs on positively and negatively charged surfaces, the RMSD values are 0.43 and 0.40 nm, respectively. The RMSD is 0.45 nm when TvL solved in the bulk solution. It indicates that the charged surfaces do not induce significant conformational change of TvL. Indeed, the conformational changes of TvL on charged surfaces are

crystal in bulk on COOH−SAM on NH2−SAM 10614

orientation

dipole (Debye)

Rg (nm)

eccentricity

− − −0.36 0.99

621 623 ± 86 681 ± 47 639 ± 40

2.15 2.32 ± 0.01 2.31 ± 0.01 2.33 ± 0.02

0.15 0.19 0.19 .

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B

charged surfaces. The conformational changes in domain 3 are relatively more distinct than in the other two domains. But this does not matter in the inactivation of TvL, as the overall conformations of TvL on charged surfaces are similar to that in bulk. This work shed lights on the orientation control of TvL on charged surfaces and provides the mechanism of laccase adsorption on charged surfaces by molecular simulations. The findings are helpful for the design and development of laccasebased electrodes.

induce side chains of residues bearing the same sign of charge to orient away from the surface. It can also be found in Table 4 that the dipole moment of TvL on COOH−SAM is larger than that on NH2−SAM. When adsorbed on NH2−SAM with the “end-on” orientation, the principal axis of TvL is nearly perpendicular to the SAM; the effect of NH2−SAM on the electric dipole moment of TvL is weaker than that of the COOH−SAM, which contributes to TvL adsorption with the “lying” orientation. Moreover, TvL carries a net charge of −4 e; the number of negatively charged residues on the protein surface is more than that of positively charged residues. Thus, when it adsorbs on COOH−SAM, the negatively charged residues on the protein surface tend to move away from the surface. That is another reason why the dipole moment of TvL on COOH−SAM is larger than that on NH2−SAM. In our previous work, we also found the same interesting phenomenon that the dipole moment of ribonuclease A, a positively charged protein, adsorbed on NH2−SAM is larger than that on COOH−SAM.55 It is worth to mention that the dipole moment changes are induced by the movement of side chains, because the backbone structure of the protein is well preserved as shown in RMSD results. Gyration radius and eccentricity are usually used together to analyze the overall morphology of a protein. The gyration radius, defined as the mass-weighted root-mean-square average distance of all atoms in a protein from its center of mass,47 demonstrates the overall size of a protein. The eccentricity represents the value of 1 − Iave/Imax, in which Imax and Iave are the maximal principal moment of inertia and average of three principal moments of inertia, respectively.49 It reveals the overall shape or ellipsoidal degree of a protein. From Table 4, we can see that the gyration radii of TvL on charged surfaces are nearly the same as those in bulk. The same result can also be found in the eccentricity analysis. This further demonstrates that the charged surface does not influence the overall morphology of TvL.



ASSOCIATED CONTENT

* Supporting Information S

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcb.7b08738. Potential energy surface and orientations in PTMC simulation and RMSD of TvL during MD simulation (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Fax: +86 20 87114069. Tel.: +86 20 87114069. ORCID

Jian Zhou: 0000-0002-3033-7785 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS Support from National Natural Science Foundation of China (Nos. 21376089, 21706197, 21606053, 21706045, 21776093, and 91334202), the National Key Basic Research Program of China (No. 2013CB733500), Guangdong Science Foundation (No. 2014A030312007), and the Fundamental Research Funds for the Central Universities (SCUT-2015ZP033) is gratefully acknowledged. An allocation time from the SCUTGrid at South China University of Technology is gratefully acknowledged.

4. CONCLUSION In this work, the orientation and the effective DET of TvL on charged surfaces (i.e., COOH−SAM and NH2−SAM) are studied by the combination of PTMC and AAMD methods. The results are consistent with previous experiment investigations which are summarized below. The DET can be achieved by orienting the T1 copper site within TvL toward the surface. Our results show that the NH2− SAM can make the T1 copper site closer to the surface. Furthermore, from the analysis of orientation distribution and interaction energy, it can be found that TvL can adsorb stably on the positively charged surface. It means that the positively charged electrode is beneficial for the DET between TvL and electrode. The adsorptions of TvL on both charged surfaces are induced by the synergy of electrostatic and vdW interactions, and the electrostatic interactions play the dominant role. The orientation of TvL on the positively charged surface (i.e., “endon” orientation) is controlled by the electric dipole of the protein. The orientation of TvL on the negatively charged surface (i.e., “lying” orientation) is caused by the combination of electrostatic interaction and the geometric feature. The conformational changes of TvL on charged surfaces are analyzed by RMSD, superimposed structures, dipole moment, gyration radius, and eccentricity. Results show that the native conformations of TvL are well preserved when it adsorbs on



REFERENCES

(1) Mayer, A. M.; Staples, R. C. Laccase: New Functions for an Old Enzyme. Phytochemistry 2002, 60, 551−565. (2) Piontek, K.; Antorini, M.; Choinowski, T. Crystal Structure of a Laccase from the Fungus Trametes Versicolor at 1.90-Angstrom Resolution Containing a Full Complement of Coppers. J. Biol. Chem. 2002, 277, 37663−37669. (3) Shervedani, R. K.; Amini, A. Direct Electrochemistry of Dopamine on Gold-Agaricus Bisporus Laccase Enzyme Electrode: Characterization and Quantitative Detection. Bioelectrochemistry 2012, 84, 25−31. (4) Olejnik, P.; Palys, B.; Kowalczyk, A.; Nowicka, A. M. Orientation of Laccase on Charged Surfaces. Mediatorless Oxygen Reduction on Amino- and Carboxyl-Ended Ethylphenyl Groups. J. Phys. Chem. C 2012, 116, 25911−25918. (5) Wang, F.; Hu, Y. R.; Guo, C.; Huang, W.; Liu, C. Z. Enhanced Phenol Degradation in Coking Wastewater by Immobilized Laccase on Magnetic Mesoporous Silica Nanoparticles in a Magnetically Stabilized Fluidized Bed. Bioresour. Technol. 2012, 110, 120−124. (6) Hublik, G.; Schinner, F. Characterization and Immobilization of the Laccase from Pleurotus Ostreatus and Its Use for the Continuous Elimination of Phenolic Pollutants. Enzyme Microb. Technol. 2000, 27, 330−336. (7) Wang, X.; Sjoberg-Eerola, P.; Immonen, K.; Bobacka, J.; Bergelin, M. Immobilization of Trametes Hirsuta Laccase into Poly(3,4-

10615

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B Ethylenedioxythiophene) and Polyaniline Polymer-Matrices. J. Power Sources 2011, 196, 4957−4964. (8) Xu, R.; Zhou, Q. J.; Li, F. T.; Zhang, B. R. Laccase Immobilization on Chitosan/Poly(Vinyl Alcohol) Composite Nanofibrous Membranes for 2,4-Dichlorophenol Removal. Chem. Eng. J. 2013, 222, 321−329. (9) Duran, N.; Rosa, M. A.; D’Annibale, A.; Gianfreda, L. Applications of Laccases and Tyrosinases (Phenoloxidases) Immobilized on Different Supports: A Review. Enzyme Microb. Technol. 2002, 31, 907−931. (10) Madhavi, V.; Lele, S. S. Laccase: Properties and Applications. Bioresources 2009, 4, 1694−1717. (11) Mao, F.; Mano, N.; Heller, A. Long Tethers Binding Redox Centers to Polymer Backbones Enhance Electron Transport in Enzyme ″Wiring″ Hydrogels. J. Am. Chem. Soc. 2003, 125, 4951− 4957. (12) Leger, C.; Bertrand, P. Direct Electrochemistry of Redox Enzymes as a Tool for Mechanistic Studies. Chem. Rev. 2008, 108, 2379−2438. (13) Gutierrez-Sanchez, C.; Pita, M.; Vaz-Dominguez, C.; Shleev, S.; De Lacey, A. L. Gold Nanoparticles as Electronic Bridges for LaccaseBased Biocathodes. J. Am. Chem. Soc. 2012, 134, 17212−17220. (14) Oliveira, T.; Fátima Barroso, M.; Morais, S.; de Lima-Neto, P.; Correia, A. N.; Oliveira, M.; Delerue-Matos, C. Biosensor Based on Multi-Walled Carbon Nanotubes Paste Electrode Modified with Laccase for Pirimicarb Pesticide Quantification. Talanta 2013, 106, 137−143. (15) Cardoso, F. P.; Aquino Neto, S. A.; Ciancaglini, P.; de Andrade, A. R. The Use of Pamam Dendrimers as a Platform for Laccase Immobilization: Kinetic Characterization of the Enzyme. Appl. Biochem. Biotechnol. 2012, 167, 1854−1864. (16) Ardhaoui, M.; Zheng, M. H.; Pulpytel, J.; Dowling, D.; Jolivalt, C.; Khonsari, F. A. Plasma Functionalized Carbon Electrode for Laccase-Catalyzed Oxygen Reduction by Direct Electron Transfer. Bioelectrochemistry 2013, 91, 52−61. (17) Trohalaki, S.; Pachter, R.; Luckarift, H. R.; Johnson, G. R. Immobilization of the Laccases from Trametes Versicolor and Streptomyces Coelicolor on Single-Wall Carbon Nanotube Electrodes: A Molecular Dynamics Study. Fuel Cells 2012, 12, 656−664. (18) Park, J. H.; Xue, H.; Jung, J. S.; Ryu, K. Immobilization of Laccase on Carbon Nanomaterials. Korean J. Chem. Eng. 2012, 29, 1409−1412. (19) Navaee, A.; Salimi, A. Graphene-Supported Pyrene-Functionalized Amino-Carbon Nanotube: A Novel Hybrid Architecture of Laccase Immobilization as Effective Bioelectrocatalyst for Oxygen Reduction Reaction. J. Mater. Chem. A 2015, 3, 7623−7630. (20) Ardhaoui, M.; Bhatt, S.; Zheng, M. H.; Dowling, D.; Jolivalt, C.; Khonsari, F. A. Biosensor Based on Laccase Immobilized on Plasma Polymerized Allylamine/Carbon Electrode. Mater. Sci. Eng., C 2013, 33, 3197−3205. (21) Wang, Q. Q.; Peng, L.; Li, G. H.; Zhang, P.; Li, D. W.; Huang, F. L.; Wei, Q. F. Activity of Laccase Immobilized on Tio2Montmorillonite Complexes. Int. J. Mol. Sci. 2013, 14, 12520−12532. (22) Zhu, J.; Sun, G. Lipase Immobilization on GlutaraldehydeActivated Nanofibrous Membranes for Improved Enzyme Stabilities and Activities. React. Funct. Polym. 2012, 72, 839−845. (23) Tavares, A. P. M.; Rodriguez, O.; Fernandez-Fernandez, M.; Dominguez, A.; Moldes, D.; Sanroman, M. A.; Macedo, E. A. Immobilization of Laccase on Modified Silica: Stabilization, Thermal Inactivation and Kinetic Behaviour in 1-Ethyl-3-Methylimidazolium Ethylsulfate Ionic Liquid. Bioresour. Technol. 2013, 131, 405−412. (24) Zhang, H.; Xun, E. N.; Wang, J. X.; Chen, G.; Cheng, T. X.; Wang, Z.; Ji, T. F.; Wang, L. Immobilization of Laccase for Oxidative Coupling of Trans-Resveratrol and Its Derivatives. Int. J. Mol. Sci. 2012, 13, 5998−6008. (25) Wang, F.; Huang, W.; Guo, C.; Liu, C. Z. Functionalized Magnetic Mesoporous Silica Nanoparticles: Fabrication, Laccase Adsorption Performance and Direct Laccase Capture from Trametes

Versicolor Fermentation Broth. Bioresour. Technol. 2012, 126, 117− 122. (26) Wang, F.; Guo, C.; Yang, L. R.; Liu, C. Z. Magnetic Mesoporous Silica Nanoparticles: Fabrication and Their Laccase Immobilization Performance. Bioresour. Technol. 2010, 101, 8931−8935. (27) Gupta, G.; Rajendran, V.; Atanassov, P. Bioelectrocatalysis of Oxygen Reduction Reaction by Laccase on Gold Electrodes. Electroanalysis 2004, 16, 1182−1185. (28) Brondani, D.; de Souza, B.; Souza, B. S.; Neves, A.; Vieira, I. C. PEI-Coated Gold Nanoparticles Decorated with Laccase: A New Platform for Direct Electrochemistry of Enzymes and Biosensing Applications. Biosens. Bioelectron. 2013, 42, 242−247. (29) Shervedani, R. K.; Foroushani, M. S. Comparative Electrochemical Behavior of Proteins; Cytochrome c, Agaricus Bisporus Laccase, and Glucose Oxidase, Immobilized onto Gold-Thiol SelfAssembled Monolayer via Electrostatic, Covalent, and Covalent Coordinate Bond Methods. Electrochim. Acta 2016, 187, 646−654. (30) Arzola, K. G.; Orive, A. G.; Arevalo, M. C.; Vazquez, L.; Creus, A. H.; Falcon, M. A. Adsorption of a Laccase from Fusarium Proliferatum on Au(111) and Hopg Electrodes: A Scanning Probe Microscopy and Electrochemical Approach. Int. J. Electrochem. Sci. 2012, 7, 1011−1026. (31) Michota-Kaminska, A.; Wrzosek, B.; Bukowska, J. Resonance Raman Evidence of Immobilization of Laccase on Self-Assembled Monolayers of Thiols on Ag and Au Surfaces. Appl. Spectrosc. 2006, 60, 752−757. (32) Gupta, G.; Rajendran, V.; Atanassov, P. Laccase Biosensor on Monolayer-Modified Gold Electrode. Electroanalysis 2003, 15, 1577− 1583. (33) Wang, Y.; Zhang, D.; He, F. R.; Chen, X. C. Immobilization of Laccase by Cu2+ Chelate Affinity Interaction on Surface Modified Magnetic Silica Particles and Its Use for the Removal of Pentachlorophenol. Chin. Chem. Lett. 2012, 23, 197−200. (34) Wang, Y.; Chen, X. C.; Liu, J.; He, F. R.; Wang, R. Immobilization of Laccase by Cu2+ Chelate Affinity Interaction on Surface-Modified Magnetic Silica Particles and Its Use for the Removal of 2,4-Dichlorophenol. Environ. Sci. Pollut. Res. 2013, 20, 6222−6231. (35) Beneyton, T.; El Harrak, A.; Griffiths, A. D.; Hellwig, P.; Taly, V. Immobilization of Cota, an Extremophilic Laccase from Bacillus Subtilis, on Glassy Carbon Electrodes for Biofuel Cell Applications. Electrochem. Commun. 2011, 13, 24−27. (36) Miyazawa, N.; Tanaka, M.; Hakamada, M.; Mabuchi, M. Molecular Dynamics Study of Laccase Immobilized on Self-Assembled Monolayer-Modified Au. J. Mater. Sci. 2017, 52, 12848−12853. (37) Peng, C.; Liu, J.; Xie, Y.; Zhou, J. Molecular Simulations of Cytochrome C Adsorption on Positively Charged Surfaces: The Influence of Anion Type and Concentration. Phys. Chem. Chem. Phys. 2016, 18, 9979−9989. (38) Liu, J.; Zhou, J. Hydrolysis-Controlled Protein Adsorption and Antifouling Behaviors of Mixed Charged Self-Assembled Monolayer: A Molecular Simulation Study. Acta Biomater. 2016, 40, 23−30. (39) Zhao, D.; Li, L.; He, D.; Zhou, J. Molecular Dynamics Simulations of Conformation Changes of Hiv-1 Regulatory Protein on Graphene. Appl. Surf. Sci. 2016, 377, 324−334. (40) Mijajlovic, M.; Penna, M. J.; Biggs, M. J. Free Energy of Adsorption for a Peptide at a Liquid/Solid Interface via Nonequilibrium Molecular Dynamics. Langmuir 2013, 29, 2919−2926. (41) Zhou, J.; Chen, S. F.; Jiang, S. Y. Orientation of Adsorbed Antibodies on Charged Surfaces by Computer Simulation Based on a United-Residue Model. Langmuir 2003, 19, 3472−3478. (42) Zhou, J.; Tsao, H. K.; Sheng, Y. J.; Jiang, S. Y. Monte Carlo Simulations of Antibody Adsorption and Orientation on Charged Surfaces. J. Chem. Phys. 2004, 121, 1050−1057. (43) Xie, Y.; Liu, M. F.; Zhou, J. Molecular Dynamics Simulations of Peptide Adsorption on Self-Assembled Monolayers. Appl. Surf. Sci. 2012, 258, 8153−8159. (44) Zhao, D.; Peng, C.; Zhou, J. Lipase Adsorption on Different Nanomaterials: A Multi-Scale Simulation Study. Phys. Chem. Chem. Phys. 2015, 17, 840−850. 10616

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617

Article

The Journal of Physical Chemistry B (45) Yu, G.; Liu, J.; Zhou, J. Mesoscopic Coarse-Grained Simulations of Hydrophobic Charge Induction Chromatography (Hcic) for Protein Purification. AIChE J. 2015, 61, 2035−2046. (46) Peng, C.; Liu, J.; Zhou, J. Molecular Simulations of Cytochrome C Adsorption on a Bare Gold Surface: Insights for the Hindrance of Electron Transfer. J. Phys. Chem. C 2015, 119, 20773−20781. (47) Zhou, J.; Zheng, J.; Jiang, S. Y. Molecular Simulation Studies of the Orientation and Conformation of Cytochrome C Adsorbed on Self-Assembled Monolayers. J. Phys. Chem. B 2004, 108, 17418− 17424. (48) Xie, Y.; Zhou, J.; Jiang, S. Y. Parallel Tempering Monte Carlo Simulations of Lysozyme Orientation on Charged Surfaces. J. Chem. Phys. 2010, 132, 065101. (49) Liu, J.; Liao, C.; Zhou, J. Multiscale Simulations of Protein G B1 Adsorbed on Charged Self-Assembled Monolayers. Langmuir 2013, 29, 11366−11374. (50) Liao, C. Y.; Xie, Y.; Zhou, J. Computer Simulations of Fibronectin Adsorption on Hydroxyapatite Surfaces. RSC Adv. 2014, 4, 15759−15769. (51) Yu, G.; Liu, J.; Zhou, J. Mesoscopic Coarse-Grained Simulations of Lysozyme Adsorption. J. Phys. Chem. B 2014, 118, 4451−4460. (52) Hudson, S.; Cooney, J.; Hodnett, B. K.; Magner, E. Chloroperoxidase on Periodic Mesoporous Organosilanes: Immobilization and Reuse. Chem. Mater. 2007, 19, 2049−2055. (53) Weber, E.; Sirim, D.; Schreiber, T.; Thomas, B.; Pleiss, J.; Hunger, M.; Glaser, R.; Urlacher, V. B. Immobilization of P450 Bm-3 Monooxygenase on Mesoporous Molecular Sieves with Different Pore Diameters. J. Mol. Catal. B: Enzym. 2010, 64, 29−37. (54) Peng, C.; Liu, J.; Zhao, D.; Zhou, J. Adsorption of Hydrophobin on Different Self-Assembled Monolayers: The Role of Hydrophobic Dipole and Electric Dipole. Langmuir 2014, 30, 11401−11411. (55) Liu, J.; Yu, G.; Zhou, J. Ribonuclease A Adsorption onto Charged Self-Assembled Monolayers: A Multiscale Simulation Study. Chem. Eng. Sci. 2015, 121, 331−339. (56) Patel, S.; Brooks, C. L. Charmm Fluctuating Charge Force Field for Proteins: I Parameterization and Application to Bulk Organic Liquid Simulations. J. Comput. Chem. 2004, 25, 1−15. (57) 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. (58) 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. (59) Hoover, W. G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A: At., Mol., Opt. Phys. 1985, 31, 1695−1697. (60) Nose, S. A Unified Formulation of the Constant Temperature Molecular Dynamics Methods. J. Chem. Phys. 1984, 81, 511−519. (61) Hess, B.; Bekker, H.; Berendsen, H. J. C.; Fraaije, J. Lincs: A Linear Constraint Solver for Molecular Simulations. J. Comput. Chem. 1997, 18, 1463−1472. (62) Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Hsing, L.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103, 8577−8593. (63) Yeh, I. C.; Berkowitz, M. L. Ewald Summation for Systems with Slab Geometry. J. Chem. Phys. 1999, 111, 3155−3162. (64) Humphrey, W.; Dalke, A.; Schulten, K. Vmd: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33−38. (65) Xie, Y.; Liao, C. Y.; Zhou, J. Effects of External Electric Fields on Lysozyme Adsorption by Molecular Dynamics Simulations. Biophys. Chem. 2013, 179, 26−34. (66) Yu, G.; Zhou, J. Understanding the Curvature Effect of Silica Nanoparticles on Lysozyme Adsorption Orientation and Conformation: A Mesoscopic Coarse-Grained Simulation Study. Phys. Chem. Chem. Phys. 2016, 18, 23500−23507. (67) Zhao, J.; Wang, Q.; Liang, G.; Zheng, J. Molecular Dynamics Simulations of Low-Ordered Alzheimer Beta-Amyloid Oligomers from Dimer to Hexamer on Self-Assembled Monolayers. Langmuir 2011, 27, 14876−14887.

(68) Yang, C.; Peng, C.; Zhao, D.; Liao, C.; Zhou, J.; Lu, X. Molecular Simulations of Myoglobin Adsorbed on Rutile (110) and (001) Surfaces. Fluid Phase Equilib. 2014, 362, 349−354.

10617

DOI: 10.1021/acs.jpcb.7b08738 J. Phys. Chem. B 2017, 121, 10610−10617