Environmental Effect on Surface Immobilized Biological Molecules

Sep 29, 2014 - ... Structures Formed by Covalently Bound, Short, β-Stranded Peptides on Self-Assembled Monolayers. Jason W. Dugger and Lauren J. Webb...
0 downloads 0 Views 574KB Size
Article pubs.acs.org/JPCB

Environmental Effect on Surface Immobilized Biological Molecules Zunliang Wang,†,‡ Xiaofeng Han,† Nongyue He,*,† Zhan Chen,*,‡ and Charles L. Brooks, III*,‡ †

State Key Laboratory of Bioelectronics, School of Biological Science and Medical Engineering, Southeast University, Si Pai Lou 2, Nanjing 210096, China ‡ Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States ABSTRACT: Our recent sum frequency generation (SFG) vibrational spectroscopic experiment (J. Phys. Chem. B 2014, 118, 2904−2912) showed that immobilized antimicrobial peptide cecropin P1 (cCP1) on a self-assembled monolayer (SAM) surface via N-terminus exhibited significantly different conformational and/or orientational behaviors when exposed to pure water vs a 50% (v/v) 2,2,2-trifluoroethanol (TFE)/water mixture. Meanwhile, our recent molecular dynamics (MD) simulations (J. Phys. Chem. B 2014, 118, 5670−5680) further revealed that the immobilized cCP1 via N-terminus in pure water largely adopts an overall bent structure lying down on the SAM surface, consistent with the SFG observation. Here, MD simulations were performed on the immobilized cCP1 on a SAM surface via N-terminus while in contact with a 50% (v/v) TFE/water mixture to further investigate the effects of environment (water vs TFE/water mixture) on the interfacial structure and orientation of immobilized peptide. The simulation results demonstrated that the immobilized cCP1 on the SAM surface via the N-terminus with two different starting states with different orientations and conformations, when exposed to a 50% (v/v) TFE/water mixture, was eventually able to maintain a linear α-helical structure, standing upright on the SAM surface. Taken with the corresponding SFG observation, our simulation results indicate that the conformational behavior of the immobilized peptide is mediated by the local hydrophobic environments resulting from the TFE aggregation around the peptide. Such knowledge can be used to regulate the surface conformation and functionality of immobilized peptides via changing surrounding chemical environments (e.g., TFE cosolvent), which is important for the microbial detection and killing based on surface-immobilized antimicrobial peptides.

1. INTRODUCTION Natural antimicrobial peptides are important in the immune system, protecting the host from microbial invasion.1−4 They can selectively bind to and disrupt the microbial cell membranes to kill microbials. Due to their strong affinity for microbial membranes, surface-immobilized antimicrobial peptide arrays can be used for rapid detection of microbial targets including bacterial and viral pathogens.5−18 For such detection arrays, stable surface structure and optimized orientation of immobilized peptides are crucial for their surface activity.19−21 Understanding such peptide structure and orientation aids in the rational design and development of peptide arrays with improved performance in microbial capture and death. Besides the physicochemical properties of the substrate surface, the conformational/orientational behavior of immobilized peptides on surfaces must be also mediated by the complex interactions between the amino acids in the peptide and the surrounding environment. Effects of chemical environment on peptide conformation and stability have been extensively investigated, and 2,2,2-trifluoroethanol (TFE) is often used as a cosolvent in peptide and protein conformational studies, due to its ability to promote and reinforce the α-helical secondary structure.22−25 © XXXX American Chemical Society

Many experimental methods have been used to characterize the surface structure of immobilized peptides, including scanning tunneling microscopy (STM),26 electrical impedance spectroscopy (EIS),27 X-ray photoelectron spectroscopy (XPS),27 Fourier transform infrared (FTIR) spectroscopy,27 circular dichroism (CD) spectroscopy,4 and sum frequency generation (SFG) vibrational spectroscopy.28−61 Among these techniques, SFG spectroscopy stands out as a powerful characterization technique with submonolayer sensitivity in probing surface-immobilized peptides at the solid/liquid interface at the molecular level and in situ. We recently studied the surface immobilized antimicrobial peptide cecropin P1 (CP1, sequence SWLSKTAKKLENSAKKRISEGIAIAIQGGPR) using SFG.62 It was found that surface-immobilized CP1 via an N-terminally modified cysteine (cCP1) exhibited a much weaker SFG amide I signal in pure water but a very strong α-helical SFG amide I signal in the 50% (v/v) 2,2,2trifluoroethanol (TFE)/water mixture.62 The observed difference between the SFG spectra in the two different solvent Received: August 24, 2014 Revised: September 27, 2014

A

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

investigate the TFE-induced interfacial conformational transition of the immobilized peptide. The immobilized peptide cCP1 monomer consists of 31 amino acid residues, with an additional cysteine residue at the N-terminus used to covalently bind the peptide to the maleimide group on the SAM. The cCP1 monomer has positively and negatively charged N-terminal NH3+ and Cterminal COO− groups, respectively. In the peptide, Lys (K) and Arg (R) residues are positively charged, Glu (E) is negatively charged, and other residues are neutral, giving each of the cCP1 molecule a total charge of +5. In Simulation I, the initial tilt angle of the α-helical peptide vs the SAM surface normal was ∼10°. The SAM was modeled with a monolayer containing 457 silane-EG 4 -maleimide [(OH) 3 Si(CH 2 ) 2 NH-CONH-EG 4 (CH2)2N(CO)2C2H2] chains on a silica substrate at a density of 4.5 chains per nm2. This is consistent with the experimental condition of 4−5 chains per nm2 for a full coverage of SAM on silica.78−80 The initial SAM structure was constructed using the CHARMM program.81−83 The maleimide group of the EG4SAM was modeled using the CHARMM27-compatible force field parameters developed by Imhof and co-workers.84 The force field parameters for each EG4 unit were directly adopted from the CHARMM27 force field.85,86 All the silicon atoms and the hydroxylated α-quartz (011) slab used for the silica substrate were modeled using the CHARMM27-compatible silica force field developed by Lopes et al.87 To build the solvent box of the TFE/water mixture, we started from a pre-equilibrated TIP3P88 water box with a side length of 19 Å produced by the CHARMM program.81−83 To obtain the 50% v/v concentration, we randomly added TFE molecules into the water box, where the water molecules overlapping within 2.0 Å around the TFE molecules were removed. The final molecular number ratio of water/TFE was 4:1 (216 water and 54 TFE molecules), producing a density of 1.20 g/mL corresponding to the reported experimental data.89 The new solvent box containing the 50% (v/v) TFE/water mixture had a dimension of 23 × 23 × 23 Å3. First, the solvent box was equilibrated for 3 ns in the NVT ensemble at 300 K, and then continued for 3 ns in the NPT ensemble. For all simulations in this study, the TFE molecules were modeled using the CHARMM-27 compatible force field developed by MacKerell group.90 Both simulation systems were solvated in the equilibrated 50% (v/v) TFE/water solvent box noted above. Each solvated system was then neutralized by adding 5 Cl− ions. For Simulation I, there were 2580 TFE and 15 566 water molecules. For Simulation II, there were 2685 TFE and 16 174 water molecules. Each simulation system was then subjected to an energy minimization procedure to remove unfavorable contacts. This procedure first involved solvent minimization using 100 steps of the steepest descent (SD)81 and 1000 steps of adopted basis Newton−Raphson (ABNR)81 minimization, with the silicon atoms fixed and the non-hydrogen atoms of the solutes (peptide and SAMs) restrained with a harmonic force constant of 40 (kcal/mol)/Å2. This was followed by five similar rounds of the SD-ABNR minimizations, during which the harmonic restraints on the heavy atoms of the solutes were then gradually reduced to zero. 2.2. Simulation Protocol. After energy minimization, all simulation systems were then gradually heated from 0 to 300 K in increments of 25 K for every 5 ps, while harmonically restraining the peptide backbone atoms and the SAM to their

environments indicates that in pure water, the immobilized cCP1 via N-terminus is likely to adopt a random coil or a lyingdown α-helical structure on the self-assembled monolayer (SAM) surface, whereas after the addition of 50% TFE to the water solvent, the immobilized peptide may refold to an α-helix with an upright orientation on the SAM surface. Although this SFG observation can capture the TFE-induced difference in the conformational/orientational behavior of surface-immobilized peptides, the details of the structural and orientational changes as well as dynamics of the immobilized peptides has not been completely understood yet. Molecular dynamics (MD) simulations can provide an important complement and more detailed information and have been widely used in understanding complex surface conformational behavior of proteins and peptides.63−72 Using MD simulation methods, we and other research groups have investigated the stabilizing effect of TFE on the α-helical structure of peptide.73−76 Almost all of these studies mainly focused on the peptide conformational behavior in bulk solutions. The effect of TFE on the structure and orientation of the peptides immobilized on surfaces is still unclear. Our recent MD simulation study on the surface-immobilized peptide in pure water demonstrated that the N-terminally immobilized cCP1 adopted an overall bent structure, lying down on a SAM surface.77 Such structural and orientational behavior is well correlated with the SFG observation that very weak SFG signals were detected from the immobilized cCP1 in pure water.62 The strong hydrophobic interaction of the free Cterminal amino acid residues with water is the dominant driving force leading to the surface parallel orientation of immobilized cCP1 and its conformational transition from an surface normal linear structure to a lying-down bent structure. Here, we further use all-atom MD simulations to investigate the effect of TFE on the interfacial structure and orientation of surface-immobilized cCP1. The structural and orientational dynamics of the immobilized peptide were derived from the MD simulations and compared to the available SFG experimental data. In particular, the local concentration of the TFE molecules surrounding each peptide residue was examined. To further understand the interaction of TFE with the surface-immobilized peptide, we analyzed the TFE structural arrangement around the peptide using the radial distribution functions. The underlined molecular mechanism of the TFE-induced structural and orientational changes of surface-immobilized peptide were also elucidated by such analysis.

2. METHOD 2.1. System Setup. Current MD simulations of surfaceimmobilized peptide cCP1 were performed in a 50% (v/v) TFE/water mixture at 300 K. In this work, two initial conformational/orientational states were considered: an αhelical structure for the entire peptide standing upright on the SAM surface (Simulation I), and a bent structure lying down on the SAM surface (Simulation II). These two states were directly adopted from the starting (at 0 ns) and the final (at 60 ns) snapshot structures of the previous MD simulations carried out in pure water, respectively.77 The MD simulation (Simulation I) that started with the α-helical peptide standing upright on the SAM surface was used to evaluate the effect of TFE on the structural and orientational stability of the surface-immobilized peptide. The MD simulation (Simulation II) that started with the bent structure lying down on the SAM surface was used to B

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

initial positions. The heating process allowed the initial relaxation of water and TFE molecules around the peptide and the SAM surface. Following this, a full MD run of 50 ns was performed for each simulation system without any restraints on the peptide. All MD simulations in this work were performed with the NAMD2.8,91 using the CHARMM27 force field.85,86 Periodic boundary conditions92 were applied for both simulations. For Simulation I, the simulation box had a dimension of 94 × 94 × 129 Å3, containing 108 178 atoms. For Simulation II, the simulation box had a dimension of 96 × 96 × 130 Å3, containing 110 858 atoms. All of the MD trajectories were generated using the NVT ensemble at 300 K. During both MD simulations, all silicon atoms were fixed, and all covalent bonds involving hydrogen atoms were constrained using the RATTLE method.93 The velocity Verlet method94 was used to integrate Newton’s equations with a time step of 2.0 fs. The short-range VDW interactions were calculated by the switch function with a twin range cutoff between 12 and 14 Å. The long-range electrostatic interactions were calculated by the force-shifting function at a cutoff distance of 14 Å. 2.3. Data Analysis. In both simulations, the coordinates and velocities of the MD trajectories were saved every 2 ps for the analysis of the MD trajectory. The root-mean-square deviations (RMSDs) of the peptide backbone atoms enable the structural changes to be examined throughout the MD simulation. The RMSD in this work was calculated by

LTC(d) =

g (r ) =

3. RESULTS AND DISCUSSION 3.1. Effect of TFE on the Stability of the SurfaceImmobilized α-Helical Peptide in Simulation I. To investigate the stabilizing effect of TFE on the conformation and orientation of surface immobilized α-helical cCP1 in a 50% (v/v) TFE/water mixture, the root mean-square deviations (RMSDs) of the backbone atoms of cCP1 in Simulation I were calculated as a function of simulation time (Figure 1a). It is clearly shown that in the TFE/water mixture the immobilized cCP1 largely remains in its starting α-helical conformation because of the very small RMSD change around 2.5 Å during the entire 50 ns MD simulation. While in water, as demonstrated by our previous MD simulations, the RMSD value of cCP1 exhibits a large increase suggesting a more substantial structural change in the peptide because of its strong C-terminal hydrophobic fluctuation.77 Moreover, to observe the local motion of each amino acid residue, the root-mean-square fluctuation (RMSF) of the backbone Cα atom of every peptide residue was calculated as a function of time throughout the 50 ns simulation (Figure 1b). It is evident that almost every amino acid residue in cCP1 experiences a small structural fluctuation with a RMSF value below 2.0 Å (except for the two residues at the C-terminus which have much larger values), indicating that the α-helical conformation of the immobilized cCP1 in TFE/ water mixture is well maintained despite the C-terminal slight fraying. To characterize the compactness of the entire peptide, the time evolution of the radius of gyration (Rg) of the cCP1 is presented (Figure 1c). It is clearly demonstrated that the Rg of the peptide constantly has a large value of 14 Å throughout the simulation. This means that the surface-immobilized cCP1 well maintains a linear structure in the TFE/water mixture, which is different from its compacted bent structure on surfaces in pure water as demonstrated in our previous simulation study.77 Furthermore, we examined the orientational change of the cCP1 by calculating the peptide tilt angle (τ) with respect to the SAM surface normal as a function of time and its distribution (Figure 2a,b). It is obviously shown that the tilt angle τ remains quite stable and only varies inside a small range of 0−10° (Figure 2a). This indicates that cCP1 more or less maintains its initial upright orientation with respect to the SAM surface during the simulation. For immobilized cCP1, although the initial conformation is the same, the MD trajectories of the

N

where ri,t is the current position of atom i at time t, and ri,0 is its initial position at time t = 0 after eliminating the rotation and translation motions; N denotes the total number of atoms of the immobilized peptide backbone. The radius of gyration (Rg) of the peptide backbone was used to evaluate the overall compactness of the peptide, which is defined as the root-mean-square distance from the center of mass of the peptide to every peptide backbone atom: N

∑i = 1 mi(ri − rCOM)2 N

∑i = 1 mi

where N is the number of atoms, mi and ri are the mass and position of atom i, and rCOM is the center-of-mass position of the peptide. The surface orientational behavior of the immobilized peptide can be characterized by the tilt angle (τ). τ is the angle between the principal axis (h) of the peptide and the SAM surface normal (which is along the z-axis), which can be defined as τ = cos−1

V N (r ,dr )mean NANB 4πr 2 dr

where r is the distance between atoms A and B. N(r,dr)mean is the average number of atom A between r and r + dr from atom B, which can be derived from a MD trajectory. NA and NB are the numbers of atom A and atom B, respectively. Here, using the RDF for TFE around peptide permits the structural arrangement of TFE around a peptide to be evaluated.

N

Rg =

× 100

wat where Vtfe mol = 0.07 and Vmol = 0.019 (L/mol), denoting the average occupied volumes of a TFE and water molecule, respectively. The radial distribution function (RDF) for two sets of atoms (e.g., A and B) can be used to characterize the occurrence probability of atom A around atom B at a given distance. The RDF was calculated according to the following equation:

∑i = 1 (ri , t − ri ,0)2

RMSD =

tfe Vmol Ntfe(d) tfe wat VmolNtfe(d) + Vmol Nwat(d)

h·ẑ |h|

where principal axis vector h is calculated as the eigenvector with the smallest eigenvalue of the inertia tensor of the peptide Cα atoms and ẑ is the unit vector along the z-axis direction which is perpendicular to the SAM surface.77,95,96 The local TFE concentration (% v/v) (LTC) around every amino acid residue was calculated from the cumulative number of TFE (Ntfe(d)) and water (Nwat(d)) molecules present within a distance of d from the backbone atom Cα using the following formula:74,76 C

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

Figure 2. (a) Tilt angles of CP1 versus the surface normal in Simulation I as a function of simulation time and (b) its corresponding distribution throughout the 50 ns MD simulation.

sharply increased to 5 Å in the 8−12 ns simulation time range, then it fluctuated around 5 Å with an amplitude of 2 Å during the period of 12−30 ns. After that, it rapidly reached a value of 8 Å at 35 ns and remained steady at 8 Å for the last 15 ns of the simulation time. This result suggests that cCP1 underwent substantial conformational changes on the SAM surface during the 50 ns MD simulation. The average RMSF per amino acid shows that the overall movement of C-terminal residues is much larger than that of the N-terminal residues. Moreover, in the 50 ns period of time, the Rg of cCP1 gradually increased from 10 to 13 Å (Figure 3c), indicating that the immobilized cCP1 tends to extend its initial bent structure to a more linear structure. In addition, the tilt angle of cCP1 gradually decreased from 90° to 50° versus the surface normal (Figure 4), indicating that the peptide tends to orient itself upright with its structure gradually refolding to a linear conformation in the presence of TFE. This trend of conformational transition is well reflected by the representative snapshots from the whole simulation (Figure 5). For cCP1, the population of local secondary structure per amino acid in every 10 ns simulation time except for the Nterminally modified cysteine is presented in Figure 6. In the TFE/water mixture, many amino acid residues adopt a 100% local α-helix conformation during the entire 50 ns simulation time. The C-terminal segment (QGGPR) is less stable, and in the last 10 ns simulation time, it tends to form a turn structure. The local α-helix conformation of the central residue Ser-13 gradually increased its population from 80% to 100% concomitantly with its coil population gradually decreasing to zero, leading to the formation of a linear helical conformation for the entire peptide. This result shows that TFE not only facilitates the stability of the α-helix immobilized on surfaces but also can promote the immobilized peptide to gradually form a linear standing up structure. 3.3. TFE Concentration around the Surface-Immobilized Peptide in Simulation II. To quantify the contacts and

Figure 1. (a) Root mean-square deviations (RMSDs) of the backbone atoms of cCP1 in Simulation I as a function of simulation time. (b) Root-mean-square fluctuations (RMSF, a measure of the average atomic mobility) of the backbone Cα atom of each amino acid residue in cCP1 in Simulation I, averaged over the whole 50 ns MD simulation. (c) Time evolution of the radius of gyration (Rg) of the backbone atoms of cCP1 in Simulation I.

cCP1 generated in the presence of TFE demonstrated completely different structural and orientational behavior of surface immobilized peptide from that generated in pure water.77 This significant difference also indicates that TFE contributes to reinforce the structural and orientational stability of an α-helical peptide immobilized on the SAM surface, which also well agrees with our recent SFG experimental result: The SFG spectra collected from cCP1 showed a much weaker amide I signal in water, whereas a very strong SFG α-helical amide I signal was generated when immersed in a 50% (v/v) TFE/water mixture.62 Such an observed difference in SFG spectra suggests that immobilized cCP1 via the N-terminus tends to maintain a more bent structure lying down on the SAM surface in water but can rapidly regain a more linear αhelical structure standing upright on the SAM surface after the addition of 50% TFE to the contacted water solution. To further reveal conformational dynamics of immobilized cCP1, we performed another 50 ns MD run in Simulation II. 3.2. TFE-Induced Conformational Transition of the Surface-Immobilized Peptide in Simulation II. To further understand the conformational dynamics of surface-immobilized cCP1, the time evolution of RMSD and the RMSF per residue were calculated in Simulation II starting from a different initial state (with peptide in a bent lying down state) to observe the overall time-dependent and local structural variations of the peptide, as shown in Figure 3a,b. For cCP1, the RMSD value D

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

Figure 3. (a) Root mean-square deviations (RMSDs) of the backbone atoms of CP1 in Simulation II as a function of simulation time. (b) Root-mean-square fluctuations (RMSF, a measure of the average atomic mobility) of the backbone Cα atom of each amino acid residue in CP1 in Simulation II, averaged over the whole 50 ns MD simulation. (c) Time evolutions of radius of gyration (Rg) of the backbone atoms of CP1 in Simulation II.

Figure 4. (a) Tilt angles of cCP1 versus the surface normal in Simulation II as a function of simulation time. (b) Tilt angle distribution in the whole 50 ns simulation and (c) the last 10 ns simulation.

specific. Our previous simulation result showed that for the surface immobilized cCP1, its bent lying-down structure in pure water is mainly caused by the hydrophobic interactions of the central and C-terminal residues with the surface.77 Thus, combined with the structural analyses as presented above, we believe that the tight association of TFE around the immobilized peptide should also be correlated with its linear α-helical structure and standing upright surface orientation. To further confirm this, we examined the interactions of TFE molecules with the peptide, especially with the hydrophobic residues such as Leu, Ala, and Ile. 3.4. Interaction of TFE with the Surface-Immobilized Peptide in Simulation II. First, radial density distribution functions (RDFs) were calculated to observe the structural arrangement of the solvent (TFE and water) molecules around the peptide. The left panel of Figure 8a shows the RDF profile of the carbon in the CF3 group of TFE and the water oxygen atom with respect to the Cα atom of the peptide. In the case of water oxygen, the RDF profile shows no obvious peak, indicating a less organized water structure surrounding the peptide. Conversely, a relatively strong hydrophobic interaction was observed between the carbon atom in the TFE CF3 group and the peptide backbone Cα atom within a distance of 4.5 Å. This difference means that the number of TFE molecules near the peptide is significantly higher than that of water indicating

interactions of TFE molecules with the surface immobilized peptide, the LTC per residue of CP1 has also been evaluated in two temporal windows including 0−10 and 40−50 ns, as shown in Figure 7a,b, respectively. The LTC per amino acid in CP1 was evaluated by calculating the relative numbers of TFE and water molecules within a solvation shell with a 6 Å surrounding each residue. In the first 10 ns simulation, the average LTC around each residue is 62% (v/v), whereas in the last 10 ns simulation, it increased to 74% (v/v). This temporal difference in LTC shows that TFE gradually aggregates around the peptide with the increase of simulation time. Additionally, it was found that the average LTC around the nonpolar (hydrophobic) residues (Leu, Ala, Ile, Gly, and Trp) is 83% (v/v), whereas the average LTC around the polar hydrophilic residues is 67% (v/v), indicating that TFE molecules have much stronger solvation preferences for the hydrophobic residues. Moreover, all highest values of LTC (higher than 90% (v/v)) correspond to the most hydrophobic (nonpolar) residues: Ala-7, Leu-10, Ala-14, Ile-18, Ile-22, and Ala-25, further suggesting that the preferred aggregation of TFE is determined by the specific hydrophobic nature of each amino acid. This preferred aggregation behavior has been observed in other studies using MD simulations in bulk solvent74,76 and reveals that aggregation of TFE around the peptide is residueE

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

Figure 5. Representative snapshots of immobilized cCP1 exposed to 50% (v/v) TFE/water mixture (Simulation II). The backbone of cCP1 is represented with green NewCartoon style. All the snapshots are visualized using VMD1.9.99

Figure 8. Radial distribution functions, averaged over the last 10 ns of the simulation time in Simulation-II for (a) (left panel) CF3 carbon of TFE (red line) and water oxygen (blue line) toward backbone Cα atom of CP1 and (right panel) TFE oxygen and water oxygen (red line), terminal (OH) hydrogen of TFE and water oxygen (blue line); (b) (left panel) CF3 carbon of TFE toward hydrophobic (nonpolar) side chain (red line) and hydrophilic (polar) side chain (blue line) of CP1 peptide and (right panel) water oxygen toward hydrophobic (nonpolar) side chain (red line) and hydrophilic (polar) side chain (blue line).

Figure 6. Secondary structure population of each amino acid residue in CP1 in Simulation II calculated with STRIDE (40) (α-helix (red), turns (green), coil (blue), and 3-10 helix (violet)) in the temporal window of (a) 0−10 ns, (b) 10−20 ns, (c) 20−30 ns, (d) 30−40 ns, and (e) 40−50 ns MD simulations in 50% (v/v) TFE/water mixtures.

nearby TFE molecules are much closer to the peptide compared to water, which is also consistent to the above LTC analysis. To further examine the interaction between the TFE molecules and the water molecules surrounding the peptide, the RDF of the distance between the water oxygen and TFE oxygen within a distance of 6 Å from the peptide backbone Cα atom was presented in the right panel of Figure 8a. The first peak probability of the oxygen−oxygen distance occurs at a short distance of 2.75 Å, indicating strong hydrogenbonding interaction between TFE and water molecules, which was accompanied by the occurrence of the first peak at 1.85 Å between the water oxygen and the hydrogen in the TFE hydroxyl group. To examine the overall orientation of TFE with respect to the peptide residues, we calculated the RDF of the CF3 carbon and oxygen of TFE toward the side chains of the hydrophobic and hydrophilic residues, as shown in the left panel and right panel of Figure 8b, respectively. It is clearly shown that the TFE CF3 carbon has a much higher distribution density around the hydrophobic peptide side chains, and a lower distribution density around the hydrophilic peptide side chains (Figure 8b, left panel), indicating that the TFE molecule interacts more favorably with the hydrophobic (nonpolar)

Figure 7. Local TFE concentration (LTC) within a distance of 6 Å from the Cα atom in each residue of CP1 exposed to the 50% (v/v) TFE/water mixture in Simulation II, averaged over (a) the first 10 ns and (b) the last 10 ns MD trajectories, respectively. Each Bar denotes the LTC standard deviations.

F

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

In all RDF profiles, the sharp peaks of the CF3 carbon, CH2 carbon, and oxygen of TFE versus the CH3 carbon are located at 4.5, 5.5, and 6.5 Å, respectively. Accordingly, the CF3 carbon around the Ala, Ile, and Leu also shows a more organized structure compared to the CH2 carbon and oxygen atom of TFE. This result strongly suggests that the TFE is preferentially oriented with the CF3 group pointing to the peptide hydrophobic group. The structure and orientation of TFE in the solvation shell supports the role of its aggregation around the peptide to stabilize the α-helix conformation. This also supports the peptide structural transition from bent to linear in response to the local environmental condition changes. The N-terminally immobilized cCP1 in water has primarily a bent lying down structure. After the addition of TFE to water, the cCP1 molecule tends to transform to a linear α-helical structure for the entire peptide with an upright orientation on the surface, in agreement with our recent SFG experimental result. Our simulation results indicate that the preferred orientation of TFE toward the nonpolar residues can largely promote TFE aggregation around these residues, which further confirms above LTC analysis. Figure 10 shows the final snapshot after 50

amino acids. In the right panel of Figure 8b, the TFE oxygen obviously shows a less organized structure near both hydrophobic and hydrophilic side chains because the RDF profile does not exhibit any obvious peak, indicating that the oxygen of TFE is farther away from the peptide than the terminal carbon of TFE. In addition, the TFE oxygen density around the hydrophobic side chains is on average larger than that around the hydrophilic side chains, suggesting that the TFE molecules adopt a consistent orientation pointing toward the hydrophobic and hydrophilic residues. Thus, we can determine the orientation of the TFE molecule from the above discussion that the hydrophobic CF3 group pointing toward the peptide and the OH group pointing away from the peptide. Because the acidity of the OH group in TFE is increased by the hydrophobic electron-withdrawing trifluoromethyl group, TFE is a good donor but a poor acceptor for hydrogen bonds. Furthermore, in the local solvation shell surrounding the peptide, the hydrogen-bonding association of the TFEhydrogen to the water-oxygen can also bring the water-oxygen closer to the TFE-oxygen. Thus, the water-oxygen in the peptide solvation shell largely points to the peptide. To further explore the nonpolar interactions of TFE with the hydrophobic residues, the RDF profiles for terminal (CF3) carbon, the central carbon and the oxygen of TFE to the CH3 carbon of the Ala, Ile, and Leu residues (corresponding to the highest LTC values) are presented in Figure 9a−c, respectively.

Figure 10. Snapshot after the 50 ns MD run in Simulation II. The backbone of CP1 is represented with green NewCartoon style. Only the solvent molecules (TFE and water) within a distance of 10 Å from the peptide backbone are shown. TFE molecules are represented with red licorice style and water molecules are represented with blue points style. The bottom SAM surface is represented with lines style, colored with its atom nature. The snapshot is visualized using VMD1.9.99

ns of simulation time in Simulation II. It is evident that TFE molecules largely aggregate around the peptide, reducing the interaction of water with the peptide.

4. CONCLUSION In the present study, we performed molecular dynamics simulations to investigate the structural and orientational changes of N-terminally immobilized peptide cCP1 exposed to 50% (v/v) TFE/water mixture. Analyses of the MD trajectories starting from two peptide initial structures, a standing-upright α-helix and a bent structure with a lying down orientation, showed that TFE is able to not only reinforce the structural and orientational stability of an α-helical peptide immobilized on the surface but also promote the peptide to

Figure 9. Radial distribution functions, averaged over the last 10 ns of the simulation time in Simulation II, for CF3 carbon (red line), CH2 carbon (green line), and oxygen (blue line) of TFE toward the CH3 carbon of (top) Ala, (middle) Ile, and (bottom) Leu hydrophobic (nonpolar) residues in the CP1 peptide. G

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

eventually form a linear α-helical structure standing upright on the SAM surface. By analyzing the local TFE concentration around each amino acid residue, we found that almost all the peptide residues were closely surrounded by the TFE molecules and the TFE molecules aggregate more favorably around the nonpolar (hydrophobic) residues of the peptide. Almost each amino acid residue has a preference for contacting TFE (instead of water), particularly for the amino acid residues containing a nonpolar side chain. Furthermore, with the help of the RDF profiles, we further examined the interactions between the solution (water and TFE) and the peptide, confirming that the preferred aggregation of TFE is mainly caused by its hydrophobic interaction with the peptide. Additionally, this result also suggests that a more organized structure and orientation of TFE molecules around the peptide promotes the TFE preferential surrounding around the peptide. In conclusion, our simulation results can be well correlated to the corresponding SFG experimental observations. The study also further revealed that the TFE aggregation around the surface immobilized peptide provides a local hydrophobic microenvironment that contributes to reduce the destabilizing effect on the structure and orientation of the surfaceimmobilized peptide caused by the hydrophobic interactions between the peptide nonpolar amino acids and water. The TFE aggregation also promotes the immobilized peptide to eventually form a linear α-helical structure with an upright orientation on the surface. Such information provides in-depth atomic-level understanding on the effect of chemical environment (e.g., water or TFE/water cosolvent) on the conformational/orientational behavior of biological molecules immobilized on surfaces.



(5) Gregory, K.; Mello, C. M. Immobilization of Escherichia coli Cells by Use of the Antimicrobial Peptide Cecropin P1. Appl. Environ. Microbiol. 2005, 71, 1130−1134. (6) Kulagina, N. V.; Shaffer, K. M.; Anderson, G. P.; Ligler, F. S.; Taitt, C. R. Antimicrobial Peptide-based Array for Escherichia Coli and Salmonella Screening. Anal. Chim. Acta 2006, 575, 9−15. (7) Kulagina, N. V.; Lassman, M. E.; Ligler, F. S.; Taitt, C. R. Antimicrobial Peptides for Detection of Bacteria in Biosensor Assays. Anal. Chem. 2005, 77, 6504−6508. (8) Ulevitch, R. J. Therapeutics Targeting the Innate Immune System. Nat. Rev. Immunol. 2004, 4, 512−520. (9) Mannoor, M. S.; Zhang, S.; Link, A. J.; McAlpine, M. C. Electrical Detection of Pathogenic Bacteria via Immobilized Antimicrobial Peptides. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 19207−19212. (10) Johnson, S.; Evans, D.; Laurenson, S.; Paul, D.; Davies, A. G.; Ko Ferrigno, P.; Wälti, C. Surface Immobilised Peptide Aptamers as Probe Molecules for Protein Detection. Anal. Chem. 2008, 80, 978− 983. (11) North, S. H.; Wojciechowski, J.; Chu, V.; Taitt, C. R. Surface Immobilization Chemistry Influences Peptide-based Detection of Lipopolysaccharide and Lipoteichoic Acid. J. Pept. Sci. 2012, 18, 366−372. (12) Houseman, B. T.; Huh, J. H.; Kron, S. J.; Mrksich, M. Peptide Chips for the Quantitative Evaluation of Protein Kinase Activity. Nat. Biotechnol. 2002, 20, 270−274. (13) Zaytseva, N. V.; Goral, V. N.; Montagna, R. A.; Baeumner, A. J. Development of a Microfluidic Biosensor Module for Pathogen Detection. Lab Chip 2005, 5, 805−811. (14) Kato, R.; Kaga, C.; Kunimatsu, M.; Kobayashi, T.; Honda, H. Peptide Array-based Interaction Assay of Solid-bound Peptides and Anchorage-dependant Cells and Its Effectiveness in Cell-adhesive Peptide Design. J. Biosci. Bioeng. 2006, 101, 485−495. (15) Shriver-Lake, L. C.; North, S. H.; Dean, S. N.; Taitt, C. R. In Designing Receptors for the Next Generation of Biosensors; Piletsky, S. A., Whitcombe, M. J., Eds.; Springer-Verlag: Berlin, Germany, 2013; pp 85−104. (16) Ivanov, I. E.; Morrison, A. E.; Cobb, J. E.; Fahey, C. A.; Camesano, T. A. Creating Antibacterial Surfaces with the Peptide Chrysophsin-1. ACS Appl. Mater. Interfaces 2012, 4, 5891−5897. (17) Soares, J. W.; Kirby, R.; Morin, K. M.; Mello, C. M. Antimicrobial Peptide Preferential Binding of E. coli O157:H7. Protein Pept. Lett. 2008, 15, 1086−1093. (18) Uzarski, J. R.; Mello, C. M. Detection and Classification of Related Lipopolysaccharides via a Small Array of Immobilized Antimicrobial Peptides. Anal. Chem. 2012, 84, 7359−7366. (19) Nakayama, H.; Manaka, T.; Iwamoto, M.; Kimura, S. Vertical Orientation with a Narrow Distribution of Helical Peptides Immobilized on a Quartz Substrate by Stereocomplex Formation. Soft Matter. 2012, 8, 3387−3392. (20) Cha, T.; Guo, A.; Zhu, X. Y. Enzymatic Activity on a Chip: The Critical Role of Protein Orientation. Proteomics 2005, 5, 416−419. (21) Cretich, M.; Damin, F.; Pirri, G.; Chiari, M. Protein and Peptide Arrays: Recent Trends and New Directions. Biomol. Eng. 2006, 23, 77−88. (22) Buck, M. Trifluoroethanol and Colleagues: Cosolvents Come of Age. Recent Studies with Peptides and Proteins. Q. Rev. Biophys. 1998, 31, 297−355. (23) Povey, J. F.; Smales, C. M.; Hassard, S. J.; Howard, M. J. Comparison of the Effects of 2,2,2-trifluoroethanol on Peptide and Protein Structure and Function. J. Struct. Biol. 2007, 157, 329−338. (24) Hong, D. P.; Hoshino, M.; Kuboi, R.; Goto, Y. Clustering of Fluorine-substituted Alcohols as a Factor Responsible for their Marked Effects on Proteins and Peptides. J. Am. Chem. Soc. 1999, 121, 8427− 8433. (25) Kumar, Y.; Muzammil, S.; Tayyab, S. Influence of Fluoro, Chloro and Alkyl Alcohols on the Folding Pathway of Human Serum Albumin. J. Biochem. 2005, 138, 335−341.

AUTHOR INFORMATION

Corresponding Authors

*N. He. E-mail: [email protected]. Phone: 86-25-83790885. *Z. Chen. E-mail: [email protected]. Phone: 1-734-615-4189. *C. L. Brooks. E-mail: [email protected]. Phone: 1-734647-6682. Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This research was supported by the National Key Program for Developing Basic Research of China (2010CB933903 and 2014CB744501), the National Natural Science Foundation of China (61271056 and 21303015), and U.S. Army Research Office (W911NF-11-1-0251). We also acknowledge Shanghai Supercomputer Center for providing computing resources.



REFERENCES

(1) Kulagina, N. V.; Anderson, G. P.; Ligler, F. S.; Shaffer, K. M.; Taitt, C. R. Antimicrobial Peptides: New Recognition Molecules for Detecting Botulinum Toxins. Sensors 2007, 7, 2808−2824. (2) Zhang, L.; Falla, T. J. In Antimicrobial Peptides, Methods in Molecular Biology; Giuliani, A., Rinaldi, A. C., Eds.; Methods in Molecular Biology Vol. 618; Humana: New York, 2010; pp 303−327. (3) Mello, C. M.; Soares, J. In Membrane selectivity of antimicrobial peptides. in Microbial Surfaces; Camesano, T. A., Mello, C. M., Eds.; ACS Symposium Series 984; American Chemical Society: Washington, DC, 2008; pp 52−62. (4) Dathe, M.; Wieprecht, T. Structural Features of Helical Antimicrobial Peptides: Their Potential to Modulate Activity on Model Membranes and Biological Cells. Biochim. Biophys. Acta 1999, 1462, 71−87. H

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

(26) Raigoza, A. F.; Webb, L. J. Molecularly Resolved Images of Peptide-Functionalized Gold Surfaces by Scanning Tunneling Microscopy. J. Am. Chem. Soc. 2012, 134, 19354−19357. (27) Shamsi, F.; Coster, H.; Jolliffe, K. A.; Chilcott, T. Characterization of the Substructure and Properties of Immobilized Peptides on Silicon Surface. Mater. Chem. Phys. 2011, 126, 955−961. (28) Shen, Y. R. The Principles of Nonlinear Optics; WileyInterscience: New York, 1984; pp 475−504. (29) Shen, Y. R. Surface Properties Probed by Second-harmonic and Sum-frequency Generation. Nature 1989, 337, 519−525. (30) Williams, C. T.; Beattie, D. A. Probing Buried Interfaces with Non-linear Optical Spectroscopy. Surf. Sci. 2002, 500, 545−576. (31) Zhuang, X.; Miranda, P. B.; Kim, D.; Shen, Y. R. Mapping Molecular Orientation and Conformation at Interfaces by Surface Nonlinear Optics. Phys. Rev. B 1999, 59, 12632−12640. (32) Lambert, A. G.; Davies, P. B.; Neivandt, D. J. Implementing the Theory of Sum Frequency Generation Vibrational Spectroscopy: A Tutorial Review. Appl. Spectrosc. Rev. 2005, 40, 103−145. (33) Perry, A.; Neipert, C.; Space, B. Theoretical Modeling of Interface Specific Vibrational Spectroscopy: Methods and Applications to Aqueous Interfaces. Chem. Rev. 2006, 106, 1234−1258. (34) Leung, B. O.; Yang, Z.; Wu, S. S. H.; Chou, K. C. Role of Interfacial Water on Protein Adsorption at Cross-linked Polyethylene Oxide Interfaces. Langmuir 2012, 28, 5724−5728. (35) Chen, P.; Kung, K. Y.; Shen, Y. R.; Somorjai, G. A. Sum Frequency Generation Spectroscopic Study of CO/Ethylene Coadsorption on the Pt(111) Surface and CO Poisoning of Catalytic Ethylene Hydrogenation. Surf. Sci. 2001, 494, 289−297. (36) Tong, Y.; Tyrode, E.; Osawa, M.; Yoshida, N.; Watanabe, T.; Nakajima, A.; Ye, S. Preferential Adsorption of Amino-terminated Silane in a Binary Mixed Self-assembled Monolayer. Langmuir 2011, 27, 5420−5426. (37) Ye, H. K.; Gu, Z. Y.; Gracias, D. H. Kinetics of Ultraviolet and Plasma Surface Modification of Poly(dimethylsiloxane) Probed by Sum Frequency Vibrational Spectroscopy. Langmuir 2006, 22, 1863− 1868. (38) Kim, J.; Cremer, P. S. IR-Visible SFG Investigations of Interfacial Water Structure upon Polyelectrolyte Adsorption at the Solid/Liquid Interface. J. Am. Chem. Soc. 2000, 122, 12371−12372. (39) Bordenyuk, A. N.; Jayathilake, H.; Benderskii, A. V. Coherent Vibrational Quantum Beats as a Probe of Langmuir-Blodgett Monolayers. J. Phys. Chem. B 2005, 109, 15941−15949. (40) Baldelli, S. Surface Structure at the Ionic Liquid-electrified Metal Interface. Acc. Chem. Res. 2008, 41, 421−431. (41) Davis, A. P.; Ma, G.; Allen, H. C. Surface Vibrational Sum Frequency and Raman Studies of PAMAM G0, G1 and Acylated PAMAM G0 Dendrimers. Anal. Chim. Acta 2003, 496, 117−131. (42) Voges, A. B.; Al-Abadleh, H. A.; Musorrariti, M. J.; Bertin, P. A.; Nguyen, S. T.; Geiger, F. M. Carboxylic Acid- and Ester-functionalized Siloxane Scaffolds on Glass Studied by Broadband Sum Frequency Generation. J. Phys. Chem. B 2005, 108, 18675−18682. (43) Can, S. Z.; Mago, D. D.; Walker, R. A. Structure and Organization of Hexadecanol Isomers Adsorbed to the Air/Water Interface. Langmuir 2006, 22, 8043−8049. (44) Liu, J.; Conboy, J. C. Direct Measurement of the Transbilayer Movement of Phospholipids by Sum-frequency Vibrational Spectroscopy. J. Am. Chem. Soc. 2004, 126, 8376−8377. (45) Rao, Y.; Comstock, M.; Eisenthal, K. B. Absolute Orientation of Molecules at Interfaces. J. Phys. Chem. B 2006, 110, 1727−1732. (46) Fu, L.; Ma, G.; Yan, E. C. Y. In Situ Misfolding of Human Islet Amyloid Polypeptide at Interfaces Probed by Vibrational Sum Frequency Generation. J. Am. Chem. Soc. 2010, 132, 5405−5412. (47) Fu, L.; Xiao, D.; Wang, Z.; Batista, V. S.; Yan, E. C. Y. Chiral Sum Frequency Generation for In Situ Probing Proton Exchange in Antiparallel β-Sheets at Interfaces. J. Am. Chem. Soc. 2013, 135, 3592− 3598. (48) Weidner, T.; Breen, N. F.; Li, K.; Drobny, G. P.; Castner, D. G. Sum Frequency Generation and Solid-state NMR Study of the

Structure, Orientation, and Dynamics of Polystyrene-adsorbed Peptides. Proc. Natl. Acad. Sci. U. S. A. 2010, 107, 13288−13293. (49) Weidner, T.; Apte, J. S.; Gamble, L. J.; Castner, D. G. Probing the Orientation and Conformation of α-Helix and β-Strand Model Peptides on Self-Assembled Monolayers Using Sum Frequency Generation and NEXAFS Spectroscopy. Langmuir 2010, 26, 3433− 3440. (50) Mermut, O.; Phillips, D. C.; York, R. L.; McCrea, K. R.; Ward, R. S.; Somorjai, G. In situ Adsorption Studies of a 14-amino Acid Leucine-lysine Peptide onto Hydrophobic Polystyrene and Hydrophilic Silica Surfaces Using Quartz Crystal Microbalance, Atomic Force Microscopy, and Sum Frequency Generation Vibrational Spectroscopy. J. Am. Chem. Soc. 2006, 128, 3598−3607. (51) Roeters, S. J.; van Dijk, C. N.; Torres-Knoop, A.; Backus, E. H. G.; Campen, R. k.; Bonn, M.; Woutersen, S. Determining in Situ Protein Conformation and Orientation from the Amide-I Sumfrequency Generation Spectrum: Theory and Experiment. J. Phys. Chem. A 2013, 117, 6311−6322. (52) Li, H.; Ye, S.; Wei, F.; Ma, S.; Luo, Y. In Situ Molecular-level Insights into the Interfacial Structure Changes of Membraneassociated Prion Protein Fragment [118−135] Investigated by Sum Frequency Generation Vibrational Spectroscopy. Langmuir 2012, 28, 16979−16988. (53) Nguyen, K. T.; King, J. T.; Chen, Z. Orientation Determination of Interfacial β-Sheet Structures in Situ. J. Phys. Chem. B 2010, 114, 8291−8300. (54) Ye, S.; Nguyen, K. T.; Boughton, A. P.; Mello, C. M.; Chen, Z. Orientation Difference of Chemically Immobilized and Physically Adsorbed Biological Molecules on Polymers Detected at the Solid/ Liquid Interfaces in Situ. Langmuir 2010, 26, 6471−6477. (55) Nguyen, K. T.; Clair, S. V. L; Ye, S.; Chen, Z. Molecular Interactions between Magainin 2 and Model Membranes in Situ. J. Phys. Chem. B 2009, 113, 12358−12363. (56) Boughton, A. P.; Andricioaei, I.; Chen, Z. Surface Orientation of Magainin 2: Molecular Dynamics Simulation and Sum Frequency Generation Vibrational Spectroscopic Studies. Langmuir 2010, 26, 16031−16036. (57) Yang, P.; Ramamoorthy, A.; Chen, Z. Membrane Orientation of MSI-78 Measured by Sum Frequency Generation Vibrational Spectroscopy. Langmuir 2011, 27, 7760−7767. (58) Han, X.; Soblosky, L.; Slutsky, M.; Mello, C.; Chen, Z. Solvent Effect and Time-Dependent Behavior of C-Terminus-CysteineModified Cecropin P1 Chemically Immobilized on a Polymer Surface. Langmuir 2011, 27, 7042−7051. (59) Boughton, A. P.; Yang, P.; Tesmer, V.; Ding, B.; Tesmer, J.; Chen, Z. Heterotrimeric G protein β1γ2 subunits change orientation upon complex formation with G protein-coupled receptor kinase 2 (GRK2) on a model membrane. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, E667−673. (60) Ye, S.; Li, H.; Wei, F.; Jasensky, J.; Boughton, A. P.; Yang, P.; Chen, Z. Observing a Model Ion Channel Gating Action in Model Cell Membranes in Real Time in Situ: Membrane Potential Change Induced Alamethicin Orientation Change. J. Am. Chem. Soc. 2012, 134, 6237−6243. (61) Nguyen, K. T.; Clair, S. V. L; Ye, S.; Chen, Z. Orientation Determination of Protein Helical Secondary Structure Using Linear and Nonlinear Vibrational Spectroscopy. J. Phys. Chem. B 2009, 113, 12169−12180. (62) Han, X.; Liu, Y.; Wu, F.; Jansensky, J.; Kim, T.; Wang, Z.; Brooks, C. L., III; Wu, J.; Xi, C.; Mello, C. M.; et al. Different Interfacial Behaviors of Peptides chemically Immobilized on Surfaces with Different Linker Lengths and via Different Termini. J. Phys. Chem. B 2014, 118, 2904−2912. (63) Norberg, J.; Nilsson, L. Advances in Biomolecular Simulations: Methodology and Recent Applications. Q. Rev. Biophys. 2003, 36, 257−306. (64) Friedel, M.; Baumketner, A.; Shea, J. E. Effects of Surface Tethering on Protein Folding Mechanisms. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 8396−8401. I

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX

The Journal of Physical Chemistry B

Article

S.; et al. CHARMM: The Biomolecular Simulation Program. J. Comput. Chem. 2009, 30, 1545−1614. (84) Bereźniak, T.; Jäschke, A.; Smith, J. C.; Imhof, P. Stereoselection in the Diels−Alderase Ribozyme: A Molecular Dynamics Study. J. Comput. Chem. 2012, 33, 1603−1614. (85) MacKerell, A. D., Jr.; Feig, M.; Brooks, C. L., III. Extending the Treatment of Backbone Energetics in Protein Force Fields: Limitations of Gas-Phase Quantum Mechanics in Reproducing Protein Conformational Distributions in Molecular Dynamics Simulations. J. Comput. Chem. 2004, 25, 1400−1415. (86) MacKerell, A. D., Jr.; Bashford, D.; Bellott, M.; Dunbrack, R. L., Jr.; 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. (87) Lopes, P. E. M.; Murashov, V.; Tazi, M.; Demchuk, E.; MacKerell, A. D., Jr. Development of an Empirical Force Field for Silica. Application to the Quartz−Water Interface. J. Phys. Chem. B 2006, 110, 2782−2792. (88) 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. (89) Chitra, R.; Smith, P. E. Properties of 2,2,2-Trifluoroethanol and Water Mixtures. J. Chem. Phys. 2001, 114, 426−435. (90) Chen, I.-J.; Yin, D.; MacKerell, A. D., Jr. Combined Abinitio/ Empirical Approach for Optimization of Lennard-Jones Parameters for Polar-Neutral Compounds. J. Comput. Chem. 2002, 23, 199−213. (91) Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kale, L.; Schulten, K. Scalable Molecular Dynamics with NAMD. J. Comput. Chem. 2005, 26, 1781− 1802. (92) Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids; Oxford University Press: New York, 1990; pp 24−27. (93) Andersen, H. C. Rattle: A “Velocity” Version of the Shake Algorithm for Molecular Dynamics Calculations. J. Comput. Phys. 1983, 52, 24−34. (94) Swope, W. C.; Andersen, H. C.; Berens, P. H.; Wilson, K. R. A Computer Simulation Method for the Calculation of Equilibrium Constants for the Formation of Physical Clusters of Molecules: Applications to Small Water Clusters. J. Chem. Phys. 1982, 76, 637− 649. (95) Im, W.; Lee, J.; Kim, T.; Rui, H. Novel Free Energy Calculations to Explore Mechanisms and Energetics of Membrane Protein Structure and Function. J. Comput. Chem. 2009, 30, 1622−1633. (96) Lee, J.; Im, W. Restraint Potential and Free Energy Decomposition Formalism for Helical Tilting. Chem. Phys. Lett. 2007, 441, 132−135. (97) Cramer, C. J. Essentials of Computational Chemistry: Theories and Models, 2nd ed.; John Wiley & Sons: Chicester, U.K., 2004; pp 84−86. (98) Benjamin, G. L.; John, E. S.; Axel, K. Fast Analysis of Molecular Dynamics Trajectories with Graphics Processing UnitsRadial Distribution Function Histogramming. J. Comput. Phys. 2011, 230, 3556−3569. (99) Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graphics 1996, 14, 33−38.

(65) Feng, J.; Wong, K. Y.; Lynch, G. C.; Gao, X. L.; Pettitt, B. M. Peptide Conformations for a Microarray Surface-Tethered Epitope of the Tumor Suppressor p53. J. Phys. Chem. B 2007, 111, 13797−13806. (66) Feng, J.; Wong, K. Y.; Lynch, G. C.; Gao, X.; Pettitt, B. M. Salt Effects on Surface-Tethered Peptides in Solution. J. Phys. Chem. B 2009, 113, 9472−9478. (67) So, C. R.; Kulp, J. L., III; Oren, E. E.; Zareie, H.; Tamerler, C.; Evans, J. S. Molecular Recognition and Supramolecular Self-Assembly of a Genetically Engineered Gold Binding Peptide on Au{111}. ACS Nano 2009, 3, 1525−1531. (68) Soliman, W.; Bhattacharjee, S.; Kaur, K. Adsorption of an Antimicrobial Peptide on Self-Assembled Monolayers by Molecular Dynamics Simulation. J. Phys. Chem. B 2010, 114, 11292−11302. (69) Verde, A. V.; Acres, J. M.; Maranas, J. K. Investigating the Specificity of Peptide Adsorption on Gold Using Molecular Dynamics Simulations. Biomacromolecules 2009, 10, 2118−2128. (70) Wang, Q.; Zhao, C.; Zhao, J.; Wang, J.; Yang, J. C.; Yu, X.; Zheng, J. Comparative Molecular Dynamics Study of Aβ Adsorption on the Self-Assembled Monolayers. Langmuir 2010, 26, 3308−3316. (71) Wang, Q.; Zhao, J.; Yu, X.; Zhao, C.; Li, L.; Zheng, J. Alzheimer Aβ1−42 Monomer Adsorbed on the Self-Assembled Monolayers. Langmuir 2010, 26, 12722−12732. (72) Zheng, J.; Li, L.; Tsao, H. K.; Sheng, Y. J.; Chen, S.; Jiang, S. Strong Repulsive Forces between Protein and Oligo (Ethylene Glycol) Self-Assembled Monolayers: A Molecular Simulation Study. Biophys. J. 2005, 89, 158−166. (73) Brooks, C. L., III; Nilsson, L. Promotion of Helix Formation in Peptides Dissolved in Alcohol and Water-Alcohol Mixtures. J. Am. Chem. Soc. 1993, 115, 11034−11035. (74) Roccatano, D.; Colombo, G.; Fioroni, M.; Mark, A. E. Mechanism by Which 2,2,2-Trifluoroethanolwater Mixtures Stabilize Secondary-Structure Formation in Peptides: a Molecular Dynamics Study. Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 12179−12184. (75) Fioroni, M.; Diaz, M. D.; Burger, K.; Berge, S. Solvation Phenomena of a Tetrapeptide in Water/Trifluoroethanol and Water/ Ethanol Mixtures: A Diffusion NMR, Intermolecular NOE, and Molecular Dynamics Study. J. Am. Chem. Soc. 2002, 124, 7737−7744. (76) Mehrnejad, F.; Naderi-Manesh, H.; Ranjbar, B. The Structural Properties of Magainin in Water, TFE/Water, and Aqueous Urea Solutions: Molecular Dynamics Simulations. Proteins 2007, 67, 931− 940. (77) Wang, Z.; Han, X.; He, N.; Chen, Z.; Brooks, C. L., III. Molecular Structures of C- and N-Terminus Cysteine Modified Cecropin P1 Chemically Immobilized onto Maleimide-Terminated Self-Assembled Monolayers Investigated by Molecular Dynamics Simulation. J. Phys. Chem. B 2014, 118, 5670−5680. (78) Kojio, K.; Ge, S. R.; Takahara, A.; Kajiyama, T. Molecular Aggregation State of n-Octadecyltrichlorosilane Monolayer Prepared at an Air/Water Interface. Langmuir 1998, 14, 971−974. (79) Maboudian, R. Surface Processes in MEMS Technology. Surf. Sci. Rep. 1998, 30, 207−269. (80) Mazyar, O. A.; Jennings, G. K.; Kane, J.; McCabe, C. Frictional Dynamics of Alkylsilane Monolayers on SiO2: Effect of 1-n-Butyl-3methylimidazolium Nitrate as a Lubricant. Langmuir 2009, 25, 5103− 5110. (81) Brooks, B. R.; Bruccoleri, R. E.; Olafson, B. D.; States, D. J.; Swaminathan, S.; Karplus, M. CHARMM: A Program for Macromolecular Energy, Minimization, and Dynamics Calculations. J. Comput. Chem. 1983, 4, 187−217. (82) MacKerell, Jr., A. D.; Brooks, B.; Brooks, C. L., I. I. I.; Nilsson, L.; Roux, B.; Won, Y.; Karplus, M. In CHARMM: The Energy Function and Its Parameterization with an Overview of the Program; Schleyer, P. V. R., Allinger, N. L., Clark, T., Gasteiger, J., Kollman, P. A., Schaefer, H. F., III, Schreiner, P. R., Eds.; The Encyclopedia of Computational Chemistry 1; John Wiley & Sons: Chichester, U.K., 1998; pp 271− 277. (83) Brooks, B. R.; Brooks, C. L., III; Mackerell, A. D., Jr; Nilsson, L.; Petrella, R. J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, J

dx.doi.org/10.1021/jp508550d | J. Phys. Chem. B XXXX, XXX, XXX−XXX