Subscriber access provided by UNIV OF DURHAM
Biological and Environmental Phenomena at the Interface
Charged Surface Regulates the Molecular Interactions of Electrostatically Repulsive Peptides by Inducing Oriented Alignment Lin Zhang, and Yan Sun Langmuir, Just Accepted Manuscript • DOI: 10.1021/acs.langmuir.7b04308 • Publication Date (Web): 22 Mar 2018 Downloaded from http://pubs.acs.org on March 23, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
Charged Surface Regulates the Molecular
2
Interactions of Electrostatically Repulsive Peptides
3
by Inducing Oriented Alignment
4
Lin Zhang, Yan Sun*
5
Department of Biochemical Engineering and Key Laboratory of Systems Bioengineering of the
6
Ministry of Education, School of Chemical Engineering and Technology, Tianjin University,
7
Tianjin 300072, China
8 9 10
ACS Paragon Plus Environment
1
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
ABSTRACT: Regulation of molecular orientation of charged dipeptides and involved
2
interactions by electrostatic repulsion from like-charged surfaces were studied using all-atom
3
molecular dynamics simulations. It was found that a charged surface can induce oriented
4
alignment of like-charged peptides, and the oriented alignment leads to enhanced electrostatic
5
repulsion between the peptide molecules. The findings are consistent with previous experimental
6
results about the inhibition of charged protein aggregation using like-charged ion-exchange resin.
7
Furthermore, the simulations provided molecular insights into this process, and demonstrated the
8
distinct regulation effect of like-charged surfaces on the molecular interactions between peptides
9
that possess an electric dipole structure. Both the charged surface and the electric dipole structure
10
of peptides were confirmed crucial for the regulation. The research are expected to facilitate the
11
rational design of surfaces or devices to regulate the behavior of amphoteric molecules such as
12
proteins for both in vivo and in vitro applications, which would contribute to the regulation of
13
protein-protein interactions and its application in life science and biotechnology.
Page 2 of 31
14 15
KEYWORDS: charged surface; like-charged peptide, oriented alignment; electrostatic
16
repulsion; molecular dynamics simulation
17
ACS Paragon Plus Environment
2
Page 3 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
Langmuir
1.
INTRODUCTION
2
Regulation of protein-protein interactions, protein adsorption and protein desorption
3
(repulsion) at various surfaces is of fundamental importance, because proteins are the most
4
abundant organic molecules in the cell and their interactions at surfaces are crucial in life science
5
and biotechnology,1-3 including but not limited to the biotransformation, transportation of matter,
6
energy transduction, cell-to-cell interactions and metabolic control.4,5 Molecular behavior of
7
proteins adsorbed at surfaces has been extensively studied because of its wide applications
8
ranging from material science to chemical biology,6-9 focusing on the orientation and
9
conformation of adsorbed proteins and consequent functionality.10-16 For instance, protein
10
adsorption at charged surfaces has been extensively investigated using molecular simulation,17-21
11
including the adsorption of charged protein to a like-charged surface.
12
Besides these adsorbed proteins, proteins repulsed from various surfaces, are also important in
13
both life sciences and bioengineering,22-26 but are not well understood. For instance, crowding
14
effect22-24 has long been considered effective in regulating protein interactions in vivo to assist
15
protein folding. The repulsive interactions between a protein and like-charged crowding agents
16
was speculated as a reason for the increased thermodynamic stability of the protein.23,27
17
However, possible enhancement on protein-protein interactions under crowded conditions is still
18
in dispute and thus an obstacle for the utilization of macromolecular crowding. Investigating the
19
influence of various surfaces on protein-protein interactions would be helpful but is usually
20
overlooked. Electrostatic forces in solutions of like-charged colloidal particles were investigated,
21
and electrostatic attraction was usually more focused.28-31 For electrostatic repulsion, a successful
22
regulation on protein-protein interactions using like-charged surfaces in vitro was experimentally
23
found32 in our previous work, which facilitated the on-pathway folding of proteins expressed as
ACS Paragon Plus Environment
3
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
inclusion bodies.33-35 The working mechanism was supposed as that the charged surfaces
2
oriented alignment of like-charged protein molecules inhibited protein-protein interactions.
3
However, experimental verification of such a mechanism, especially the electrostatic repulsion
4
oriented alignment36 was a challenge. Molecular dynamics (MD) simulation, as a powerful tool
5
that can offer clear microscopic information in a direct manner,6,37 even the orientational
6
alignment of proteins,38 was then used to explore the molecular insights into this process in our
7
previous work.39 Standing orientation of lysozyme molecules near a like-charged surface was
8
confirmed. However, complicated structure of protein molecules7,14,40,41 limited the investigation
9
on the oriented alignment at surfaces and consequent regulation on the molecular interactions.
10
Therefore, dipeptides were designed and used instead of protein molecules in this work.
11
Chemically, protein molecule is composed of one or more polypeptide chains with complicated
12
structures. The dipeptide has much simpler structure than that of a protein, but it can mimic the
13
amphoteric feature of proteins, which is expected important for the formation of orientation
14
controlled by electrostatic interactions. So the dipeptide was used as a dipole model of protein to
15
examine the protein orientation and molecular interactions at like-charged surface. Herein, four
16
dipeptides were designed, including KL, LK, LL, and LE, as shown in Figure 1a. Among these
17
peptides, both KL and LK are positively charged, but the former has larger electric dipole than
18
the latter; LL has several charged sites but is overall neutral; LE is negatively charged. Three
19
surfaces were designed for the simulations (Figure 1b), including one neutral, one positively
20
charged and one negatively charged. All-atom MD simulations were then performed to examine
21
the behavior of the peptides near the surfaces, focusing on the molecular orientations and
22
molecular interactions between peptides.
ACS Paragon Plus Environment
Page 4 of 31
4
Page 5 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1 2
Figure 1. All-atom models of dipeptides (a) and surfaces (b). Four dipeptides were considered.
3
KL has a net charge of +1 at pH 7 because it has 2 positive charges and 1 negative charge. LK
4
has also a net charge of +1 contributed by 1 positive charge at the N-terminal. LL is neutral. LE
5
has a net charge of -1. Because the dipeptide has heterogeneous charge distribution, the atoms in
6
each dipeptide are divided into two groups, one includes all atoms with positive charge and the
7
other includes all atoms with negative charge. The center of mass of each group was calculated.
8
The equilibrated distance between these two groups was calculated to evaluate the electric
9
dipole, as shown below each peptide. All atoms were colored as CPK model and prepared using
10
Rasmol.42
11
2. MODELS AND METHODS
12
2.1. Model Construction. The AA models of dipeptides were constructed and the three
13
dimensional structures were determined using CharMM program (http://www.charmm-gui.org/),
14
as shown in Figure 1a. The net charge and electric dipole of each dipeptide are shown in Figure
15
1a.
16
To construct the simulation system in aqueous solution, two dipeptides were randomly placed
17
into a simulation box with a size of 6.0 × 6.0 × 6.0 nm3. The solvent molecules were then added.
ACS Paragon Plus Environment
5
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 6 of 31
1
The water molecule was treated using SPC model. The system was neutralized by adding Na+ or
2
Cl- as counter ions.
3
Charged surface was modeled as a neutral planar surface with immobilized positively or
4
negatively charged ligands, following the structure of ion-exchange resin.32 The surface is
5
composed of repetitive hexagonal unit43,44 with a side length of 0.152 nm (Figure 1b). For the
6
ligand, its topology was generated by the Dundee PRODRG 2.5 server (beta)45
7
(http://davapc1.bioch.dundee.ac.uk/cig-bin/prodrg_beta).
8
distribution were defined based on the parameters of similar structure in literature46 with minor
9
modification. Thereafter, ligands were immobilized to the matrix through a C-O-C bridge
10
(Figure 1b) with a ligand density of 0.21 mmol/mL according to the experimental data.32 In
11
positively/negatively charged surface, each ligand has a protonation state with a net charge of
12
+1/-1. Meanwhile, a surface with neutral ligands (with a net charge of zero) was constructed.
Its
charge
group
and
charge
13
To construct the simulation system containing charged/neutral surface, the surface was placed
14
at the bottom of a simulation box with a size of 7.3 × 6.3 × 6.0 nm3, and an inert surface was
15
placed at the top of the box to keep all molecules in the box. Two dipeptides equilibrated in
16
aqueous solution were used and put into the simulation box with the equilibrated conformations
17
in aqueous solution. That is, the dipeptides were firstly equilibrated in aqueous solution, and then
18
the equilibrated structures were used as the initial conformations at surface. Herein, rotation was
19
performed to ensure that each dipeptide had a same initial distance of 1.5 nm from the surface.
20
Solvent molecules around dipeptides in the equilibration stage were kept except those out of the
21
new simulation box. Supplement of solvent molecules was then performed to make correct
22
density in the simulation box. The system was neutralized by adding Na+ or Cl- as counter ions.
ACS Paragon Plus Environment
6
Page 7 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
These were then placed in the center of a cuboid box of size 7.3 × 6.3 × 100 nm3 for MD
2
simulation.
3
2.2. Simulation Method. MD simulations in the NVT ensemble were performed using
4
Gromacs 4.5.3 package (http://www.gromacs.org/)47,48 and Gromos96 43a1 force field.49 A 2 fs
5
time step was used to integrate the equations of motion. The cut-off of the Lennard-Jones (LJ)
6
potential was set to 1.2 nm. Standard PME function was used to deal with the Coulomb potential.
7
Periodic boundary was used in the x, y and z directions. Neighbor list was updated every 10
8
steps. The temperature was controlled at 298.15 K by the Berendsen method with a time constant
9
of 0.5 ps. The initial velocities of beads were generated according to the Maxwell distribution at
10
system temperature. Simulations were run for 20 ns in aqueous solution firstly, to examine the
11
behavior of dipeptides in bulk solution, and to obtain the equilibrated structures to provide the
12
initial conformations for the following 5 ns simulation at surfaces, as described in Section 2.1.
13
Sixty-four independent simulations were performed for each set of conditions. The snapshots
14
were prepared using Rasmol42 (http://www.umass.edu/microbio/rasmol/).
15
2.3. Dipeptide Orientation Analysis. The angle (θ) between the direction of dipeptide
16
dipole and z axis was calculated to evaluate the orientation of dipeptide, as reported in previous
17
work.39,50 When θ=0, the dipole is perpendicular to the surface with the negative pole facing the
18
surface. When θ=-90 or 90, the dipole is parallel to the surface (the dipole with positive x
19
coordinate is defined as plus, the other minus). θ=180° or -180° indicates a dipole also
20
perpendicular to the surface but with the positive pole facing the surface. A parameter P⊥ was
21
defined to represent the probability of dipeptide molecule with a θ between -30° to 30° to
22
evaluate the dipeptides perpendicular to the positively charged surface. It should be noted that P⊥
23
has been normalized by the volume fraction at certain θ in three dimensional space (see
ACS Paragon Plus Environment
7
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 8 of 31
1
Appendix S1, Supporting Information). Over negatively charged surface, P⊥ was defined to θ150° instead.
3
2.4. Dipeptide Movement Analysis. The z coordinate of the center of mass (com) of the
4
dipeptide (zcom) was monitored to depict the location of dipeptide. The distribution probability of
5
dipeptide, P(zcom,θ), was then calculated and shown by the filling colors. Herein, P(zcom,θ) was
6
calculated based on a statistics in the whole simulation to take all possible angle values. That is,
7
the dipeptide molecules with a certain zcom and θ were counted, and the obtained number was
8
divided by the total number of dipeptides sampled in the whole simulation to calculate P(zcom,θ).
9
Therefore, in the distribution map, at a certain zcom, the range of all possible θ was marked. The
10
conformational change of the peptides at the exclusion was evaluated using root mean square
11
deviation (RMSD) and radius of gyration (Rg). RMSD representing the structural change of a
12
molecule at time t with respect to a reference structure (herein the initial structure is used), was
13
calculated using g_rms provided by Gromacs. Rg representing the compactness of a molecule
14
was calculated using g_gyrate provided by Gromacs.
15
2.5. Molecular Interaction Analysis. The potential energies, including LJ potential and
16
Coulomb potential, were calculated using auxiliary program g_energy provided by Gromacs to
17
evaluate the molecular interactions between the dipeptide molecules or between the dipeptide
18
and the surface. Herein, hydrophobic interaction was considered and included in the LJ potential
19
energy. Hydrogen bond was analyzed using auxiliary program g_hbond provided by Gromacs.
20
The LJ potential between atoms i and j is calculated by Eq. (1).
VLJ (rij ) =
Cij12 rij12
−
Cij6
(1)
rij6
ACS Paragon Plus Environment
8
Page 9 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
rij
1
where
2
atom types.
3
is the distance between atoms i and j;
and
Cij6
are parameters depend on pairs of
The Coulomb potential between atoms i and j is calculated by Eq. (2). VC (rij ) = f
4
Cij12
qi q j
(2)
ε r rij
where f is electric conversion factor; ε r is dielectric constant; qi is the charge carried in atom i.
5 6
3. RESULTS AND DISCUSSION
7
3.1. Oriented Alignment of Dipeptides over Like-charged Surface. Two dipeptides,
8
equilibrated in aqueous solution, were put onto different surfaces with an initial distance of 1.5
9
nm to examine the regulation of protein behaviors by surfaces.
10
Quantitative description based on a statistics of 64 independent simulations was obtained
11
through P ⊥ and zcom values (Section 2.3). Low P ⊥ was observed in the aqueous solution,
12
indicating no preferred orientation but only random distribution (Figure S1a in Supporting
13
Information). At the neutral surface (Figure 2a), little change of zcom was observed, indicating
14
that the peptides randomly moved around their initial positions. The P⊥ values were small with
15
random changes, confirming that no oriented alignment was formed. This was confirmed by the
16
orientation and location of the peptides in typical snapshots (Figure S2) in a direct manner. No
17
obvious movement of the peptides was observed at a neutral surface (Figure S2a); KL moved
18
around the initial position and rotated randomly, and no preferred orientation was observed in the
19
snapshots. In contrary, at the like-charged surface, the P⊥ of KL changed from initial values to
20
approximately 1.0 immediately (within 0.02 ns), confirming the formation of oriented alignment
21
where the peptides were perpendicular to the surface with the negative pole facing the surface,
ACS Paragon Plus Environment
9
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
consistent with the snapshots in Figure S2b (snapshot at 0.02 ns). Slight conformational change
2
of the dipeptide was observed, as indicated by the increase of RMSD and Rg values (Figure S3).
3
It was attributed to the simultaneous electrostatic attraction at the negative pole of the dipeptide
4
and the electrostatic repulsion at the positive pole in the opposite directions, as a result of the
5
electric dipole structure. Meanwhile, the peptides were repulsed and then went far away, which
6
was consistent with the theoretical calculation of the equilibrated position of electric dipole near
7
a charged dot (see Appendix S2). In this process (rapid increase of zcom value), the P⊥ decreased
8
and reached an average value of 0.26 at 5 ns with fluctuations, indicating the diminishing of
9
oriented alignment following the fast exclusion of the peptides by the like-charged surface.
10
Decrease of Rg was also observed (Figure S3), indicating the diminishing of conformational
11
change, which is consistent with the protein behavior at desorption stage in chromatography in
12
our previous work.13,51-54
13
The distribution probability of the peptides was drawn for further evaluation of the orientation
14
based on a statistics of 64 independent simulations. Random distribution in aqueous solution was
15
observed (Figure S1b), indicated by the wide distribution with low probabilities. Similarly, wide
16
distribution of θ with low probabilities was observed at the neutral surface (Figure 2b),
17
indicating no preferred θ for the peptides. Meanwhile, the distribution region of the peptides at
18
the neutral surface is located at zcom = [0.5,1.5], which is near the initial value of zcom. That is, the
19
peptides maintained its initial state without significant change of orientation, consistent with the
20
snapshots in Figure S2. However, at the like-charged surface (Figure 2c), a dominant state with
21
zcom = [1.5,2.5] and θ = [-20°,20°] was observed, indicating the formation of oriented alignment
22
of the peptides. The preferred angle distribution near the like-charged surface located at the
23
range corresponding to the oriented perpendicular direction. As compared to the distribution
ACS Paragon Plus Environment
Page 10 of 31
10
Page 11 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
probability at neutral surface, an obvious change of peptides to approximately perpendicular
2
orientation can be concluded. However, this oriented alignment diminished with increased
3
distance from the surface. That is, the convergent angle distribution at small zcom became wider at
4
larger zcom (see the distribution in red, Figure 2c). According to the Coulomb potential [Eq. (2)],
5
the electrostatic interaction decreases with the increase of distance. Then the diminishing of
6
oriented alignment further highlighted the important role of electrostatic interaction from the
7
charged surface in this process.
8 9
Figure 2. Molecular behaviors of KL (a-c) and LE (d-f) at a neutral surface and a like-charged
10
surface. Time courses of P⊥ and zcom are shown in (a) for KL and (d) for LE, the distribution
11
probability near neutral surface represented as a zcom vs θ plot is shown in (b) for KL and (e) for
ACS Paragon Plus Environment
11
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
LE, and the distribution probability at the like-charged surface is shown in (c) for KL and (f) for
2
LE. The distribution probability is shown by the filling colors. Herein, because of the negatively
3
charged surface, the P⊥ for LE was defined as the probability of the dipeptide with a θ ≤ -150° or
4
≥ 150°, corresponding to an oriented alignment where the dipole is perpendicular to the surface
5
but the positive pole facing the surface. Schematic diagrams composed of representative
6
snapshots are provided in the figure top, where the dipoles of the dipeptides, at a direction from
7
red (negative pole) to blue (positive pole), is shown for better examination of their orientations.
8
The dipole moment is enlarged for clear view.
9
For further verification, the molecular behavior of the negatively charged LE at a negatively
10
charged surface was examined, as shown in Figures 2d to 2f. The P⊥ of peptides changed to 1.0
11
immediately (Figure 2d), indicating the formation of oriented alignment where the peptides were
12
perpendicular to the surface with the positive pole facing the surface, similar with the oriented
13
alignment of KL over the positively charged surface (Figure 2a). Thereafter, P⊥ decreased with
14
the rapid increase of zcom values, and reached an average value of 0.18 at 5 ns with fluctuations,
15
indicating the diminishing of oriented alignment following the fast exclusion of peptides by the
16
like-charged surface. In the distribution probability of LE at negatively charged surface (Figure
17
2f), a dominant state with zcom = [1.7,2.5] and θ ≤ -170° or ≥ 160° was observed, proving the
18
formation of like-charged surface oriented alignment.
19
3.2. Regulation of Intermolecular Interactions by Charged Surface. Molecular
20
interactions between dipeptides were monitored, including the electrostatic interactions,
21
hydrophobic interactions and hydrogen bonds. However, the value of LJ potential energy was
22
much smaller than that of Coulomb potential energy (see Figure S4), because the LJ potential
23
energy decreases much faster with the increase of distance, as shown in Eqs. (1) and (2).
ACS Paragon Plus Environment
Page 12 of 31
12
Page 13 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
Hydrogen bond is also only effective within a small distance of about 0.35 nm. Therefore,
2
hydrophobic interactions and hydrogen bonds were not discussed herein. The Coulomb potential
3
energy and the distance between two dipeptides were used to evaluate the effect of charged
4
surface on the molecular interactions between the two peptides. Herein, △ d=d-d0 was used
5
instead of d to avoid the disturbance of initial distance between dipeptides (d0). EC around zero
6
was observed with significant fluctuations in aqueous solution (Figure S1c), indicating random
7
collision of the peptides in bulk solution. At the neutral surface (Figure 3), EC had significant
8
fluctuations near zero but without a definite trend, indicating that random movement did not
9
influence the molecular interactions between the peptides. In contrary, at the like-charged surface,
10
for both positively charged KL (Figure 3a) and negatively charged LE (Figure 3b), immediate
11
increases of EC were observed, indicating that the oriented alignments (Figure 2) indeed
12
enhanced the repulsion between the peptides, which was also confirmed by the larger △d than
13
that at the neutral surface. However, as the discussion refer to Figure 1, the oriented alignment
14
diminished when the dipeptides went away from the surface, because the electrostatic interaction
15
decreased with increased distance according to the Coulomb potential [Eq. (2)]. As a result, the
16
enhanced repulsion between the peptides diminished. △ d of KL (LE) dipeptides near like-
17
charged surface decreased and became similar as that above neutral surface at the last part of
18
simulation. This further confirmed that the like-charged surface oriented alignment was the
19
reason for increased molecular repulsion between the peptides, and this can explain the
20
experimental observation of reduced protein aggregation in a suspension with like-charged
21
particles.32,35,55,56 From the experimental results, enhanced refolding of charged protein at like-
22
charged surface was confirmed. However, the exploration of molecular mechanism was still a
23
challenge using experimental approaches. From the simulation results, the like-charged surface
ACS Paragon Plus Environment
13
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
oriented alignment and consequent enhancement of electrostatic repulsion between the
2
dipeptides were confirmed. That is, the simulation results can explore the molecular mechanism
3
about the function of charged surface directly. Based on these findings, charged surface can be
4
designed and used to induce microscopic oriented alignment of peptides, which leads to the
5
enhanced electrostatic repulsion between peptides, and accounted for the macroscopic
6
experimental results of reduced protein aggregation by like-charged particles.
Page 14 of 31
7 8
Figure 3. Molecular interactions between KL (a) or LE (b) at a neutral or a like-charged surface,
9
as shown by the time courses of EC and △d. △d is the increase of d as compared to the initial
10
value. Schematic diagrams composed of representative snapshots are provided in the figure top.
11
3.3. The Dependence of Oriented Alignment on Electric Dipole. As mentioned above, a
12
like-charged surface was necessary to realize the regulation on protein orientation. Besides this,
13
the electric dipole of peptides may also be crucial. Simulations using other dipeptides with
14
different dipoles and net charges were therefore performed to evaluate the contribution of the
ACS Paragon Plus Environment
14
Page 15 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
electric dipole and net charge. LK with the same net charge of +1 as KL but lower charge density
2
was used, together with a neutral dipeptide LL (see Figure 1).
3 4
Figure 4. Molecular behaviors of LK and LL over a positively charged surface, where the time
5
courses of P⊥ and zcom are shown in (a), the distribution probability represented as zcom vs θ plot
6
is shown in (b) for LK and (c) for LL, and the time courses of EC and △d values are shown in
7
(d).
8
The P⊥ and zcom were calculated, as shown in Figure 4a. Similar results were observed for LK
9
as those for KL (Figure 2). Oriented alignment was formed immediately and then diminished in
10
the exclusion process of LK. For LL, however, a significant difference was observed. The
ACS Paragon Plus Environment
15
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
peptides were adsorbed onto the surface very quickly rather than being repulsed, as indicated by
2
the zcom. This was consistent with the theoretical calculation in Appendix S2 (Supporting
3
Information). Obvious oriented alignment of the adsorbed LL molecules was observed at the
4
surface. That was attributed to the heterogeneous charge distribution of LL. That is, LL should
5
also be considered as an electric dipole with distinct negatively charged part (negative pole) and
6
positively charged part (positive pole). The attraction of negative part at surface induced the
7
formation of oriented alignment. Moreover, LL has an overall neutral charge, so the attraction
8
became stronger due to smaller distance as compared with the repulsion of positively charged
9
part (positive pole), leading to the adsorption rather than repulsion of the whole LL molecule.
10
The distribution probability was calculated to evaluate the orientation of the peptides at the
11
charged surface based on a statistics of 64 independent simulations, as shown by the filling
12
colors in Figures 4b and 4c. For LK, a dominant distribution located in the region of θ = [-30°,
13
30°] was observed, indicating the formation of a perfect oriented alignment where the peptides
14
were perpendicular to the surface with the negative pole facing the surface. More converged state
15
around θ = 0° was observed for LL with a small zcom, indicating the formation of oriented
16
alignment on the surface as well as the dominant attraction. The more converged state further
17
confirmed the enhanced electrostatic attraction due to the small distance from the surface.
18
The Coulomb potential energy and the distance between two dipeptides were monitored to
19
examine the contribution of charged surface on inhibiting the aggregation of the peptides, as
20
shown in Figure 4d. For LK, rapid increase of △d was observed, followed by an immediate
21
increase of EC, indicating the oriented alignment indeed enhanced the repulsion between the
22
peptides. As compared to KL (Figure 3a), however, the increase of EC was smaller. This clearly
23
indicates that the smaller dipole resulted in a weaker enhancement of repulsion between the
ACS Paragon Plus Environment
Page 16 of 31
16
Page 17 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
peptides at like-charged surface. For LL, the increase of EC was also observed, but without
2
significant increase of △d, because the peptides were adsorbed onto the surface and could not
3
move freely (see zcom in Figure 4a). Similar results were obtained for the positively charged KL
4
at a negatively charged surface (see Figure S5).
5
Therefore, the results indicated that the net charge of protein determines the adsorption or
6
repulsion at a charged surface. Charged protein is repulsed from a like-charged surface, such as
7
the repulsion of positively charged KL by a positively charged surface, or negatively charged LE
8
by a negatively charged surface (Figure 2). In contrary, protein without effective repulsion by the
9
surface will be adsorbed, such as the adsorption of neutral LL on a positively charged surface
10
(Figure 4), or positively charged KL at a negatively charged surface (Figure S5). The
11
dependence of adsorption or repulsion on the net charge of protein is consistent with the
12
theoretical calculation in Appendix S2 (Supporting Information). Moreover, the regulation of
13
intermolecular interactions of adsorbed peptides was different from that for repulsed peptides.
14
Because the electrostatic interaction decreased with increased distance according to the Coulomb
15
potential [Eq. (2)], the diminishing of regulation was observed during the excluding process of
16
repulsed peptides, as the discussion refer to Figure 1. For the adsorbed peptides, the maintained
17
small distance from the surface caused large electrostatic interaction and strong regulation. For
18
example, LL adsorbed at the surface make it obtain a large P ⊥ above 0.6 and a large △ d,
19
indicating a good spatial isolation of peptides. The spatial isolation of proteins on ion-exchange
20
chromatography (IEC) was widely described in experimental results, and was considered as a
21
driving force for the protein refolding in IEC.57-59 Protein adsorption driven by electrostatic
22
attraction from the IEC resin is the key to separate the bound protein molecules and helpful for
23
the protein refolding.
ACS Paragon Plus Environment
17
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
Moreover, the results indicated that the electric dipole of a protein is crucial for the formation
2
of charged surface oriented alignment. A protein with heterogeneous charge distribution can
3
possess preferred orientation at a charged surface. As evidenced by the results in the present
4
study, all the four dipeptides, KL, LK, LL and EK, can form oriented alignment at a charged
5
surface due to their electric dipoles, regardless of being adsorbed or repulsed. The attraction from
6
the surface to the part of a dipeptide with opposite charge is the determinant for the preferred
7
adsorption sites. For example, the negative pole of LL determined the preferred orientation of LL
8
and consequent oriented alignment. The adsorption of protein on a charged surface due to its
9
heterogeneous surface has been extensively reported in the simulation results.17-21 The preferred
10
adsorption sites and orientation have been discussed. The simulation results of LL on a charged
11
surface, featured by the preferred adsorption sites and oriented alignment, are consistent with the
12
results reported in literature.57-59 Moreover, the regulation of molecular interactions between
13
peptides by the oriented alignment was highlighted herein, which accounted for the experimental
14
observation about the inhibition of protein aggregation by like-charged resins.32 This may be the
15
reason for protein folding rather than aggregation in crowded conditions, which can interpret the
16
crowding effect. In crowded conditions, proteins are always near various surfaces (including
17
various charged surfaces) provided by the crowding agents. The induced oriented alignment of
18
protein by such charged surfaces can be expected due to the heterogeneous charge distribution
19
(electric dipole) of protein molecules, which would then enhance the repulsion between proteins
20
and thus help enhance the inhibition of their aggregation and improve the on-pathway protein
21
folding.
Page 18 of 31
22 23
ACS Paragon Plus Environment
18
Page 19 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
Langmuir
4. CONCLUSIONS
2
In summary, charged surfaces can induce the formation of oriented alignment of peptides,
3
leading to enhanced electrostatic repulsion between peptides or proteins, which is considered as
4
the reason for the inhibition of aggregation and the enhancement of on-pathway protein folding.
5
The oriented alignment of peptides diminished with increasing distance from the like-charged
6
surface, and can be enhanced by increasing electric dipole of the peptides. The simulation results
7
explained the experimental observation of reduced protein aggregation in a suspension with like-
8
charged particles, and confirmed the like-charged surface oriented alignment and consequent
9
enhanced molecular repulsion between the peptides at molecular level. These results have
10
explored the working mechanism of charged surfaces on the regulation of molecular interactions
11
between peptides, which would be helpful for interpreting the crowding effect in vivo and the
12
enhanced on-pathway folding of high concentration proteins by like-charged additives in vitro.
13
The research will thus facilitate the rational design of surfaces or devices for regulating the
14
interfacial behaviors of protein molecules both in vitro and in vivo.
15 16
ASSOCIATED CONTENT
17
Supporting Information.
18
Appendix about the derivation of volume fraction in three dimensional space, appendix about the
19
calculation of the equilibrated position of electric dipole near a charged dot, molecular behaviors
20
of KL and LE in aqueous solution, snapshots of two KL dipeptides over a neutral surface and a
21
charged surface, LJ potential energies of KL and LE dipeptides at a like-charged surface, the
22
time courses of RMSD and Rg of KL and LE dipeptides at a like-charged surface, and molecular
ACS Paragon Plus Environment
19
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1
behaviors of KL at a negatively charged surface. This material is available free of charge via the
2
Internet at http://pubs.acs.org.
3
AUTHOR IMFORMATION
4
Corresponding author
5
*Tel & Fax: +86 22 27403389; E-mail addresses:
[email protected] (Y. Sun).
6
Author Contributions
7
The manuscript was written through contributions of all authors. All authors have given approval
8
to the final version of the manuscript.
9
Notes
10
The authors declare no competing financial interest.
11
ACKNOWLEDGMENTS
12
This work was supported by the Natural Science Foundation of China (Nos. 91534119,
13
21236005, 21376173 and 21621004) and the Innovation Foundation of Tianjin University.
14 15
REFERENCES
16
(1) Belfort, G.; Lee, C. S. Attractive and Repulsive Interactions Between and within Adsorbed Ribonuclease a
17 18 19 20 21 22 23
Page 20 of 31
Layers. P. Natl. Acad. Sci. U.S.A. 1991, 88, 9146–9150. (2) Scott, D. E.; Bayly, A. R.; Abell, C.; Skidmore, J. Small Molecules, Big Targets: Drug Discovery Faces the Protein-Protein Interaction Challenge. Nat. Rev. Drug Discov. 2016, 15, 533–550. (3) Keskin, O.; Gursoy, A.; Ma, B.; Nussinov, R. Principles of Protein-Protein Interactions: What are the Preferred Ways for Proteins to Interact? Chem. Rev. 2008, 108, 1225–1244. (4) Yang, H.; Yuan, B.; Zhang, X.; Scherman, O. A. Supramolecular Chemistry at Interfaces: Host-Guest Interactions for Fabricating Multifunctional Biointerfaces. Accounts Chem. Res. 2014, 47, 2106–2115.
ACS Paragon Plus Environment
20
Page 21 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1 2 3 4 5 6 7 8
Langmuir
(5) Rao, V. S.; Srinivas, K.; Sujini, G. N.; Kumar, G. N. S. Protein-Protein Interaction Detection: Methods and Analysis. International Journal of Proteomics 2014, 2014, 147648–147648. (6) Ozboyaci, M.; Kokh, D. B.; Corni, S.; Wade, R. C. Modeling and Simulation of Protein-Surface Interactions: Achievements and Challenges. Q. Rev. Biophys. 2016, 49, 1–45. (7) Zhang, L.; Sun, Y. Molecular Simulation of Adsorption and its Implications to Protein Chromatography: A Review. Biochem. Eng. J. 2010, 48, 408–415. (8) Duval, F.; van Beek, T. A.; Zuilhof, H. Key Steps Towards the Oriented Immobilization of Antibodies Using Boronic Acids. Analyst 2015, 140, 6467–6472.
9
(9) Murugan, P.; Krishnamurthy, M.; Jaisankar, S. N.; Samanta, D.; Mandal, A. B. Controlled Decoration of the
10
Surface with Macromolecules: Polymerization On a Self-Assembled Monolayer (Sam). Chem. Soc. Rev. 2015,
11
44, 3212–3243.
12
(10) Penna, M. J.; Mijajlovic, M.; Biggs, M. J. Molecular-Level Understanding of Protein Adsorption at the
13
Interface Between Water and a Strongly Interacting Uncharged Solid Surface. J. Am. Chem. Soc. 2014, 136,
14
5323–5331.
15 16
(11) Puddu, V.; Perry, C. C. Peptide Adsorption On Silica Nanoparticles: Evidence of Hydrophobic Interactions. ACS Nano 2012, 6, 6356–6363.
17
(12) Lin, D. Q.; Tong, H. F.; Wang, H. Y.; Shao, S.; Yao, S. J. Molecular Mechanism of Hydrophobic Charge-
18
Induction Chromatography: Interactions Between the Immobilized 4-Mercaptoethyl-Pyridine Ligand and Igg.
19
J. Chromatogr. A 2012, 1260, 143–153.
20
(13) Zhang, L.; Zhao, G. F.; Sun, Y. Molecular Insight Into Protein Conformational Transition in Hydrophobic
21
Charge Induction Chromatography: A Molecular Dynamics Simulation. J. Phys. Chem. B 2009, 113, 6873–
22
6880.
23 24 25 26 27 28
(14) Watkins, A. M.; Bonneau, R.; Arora, P. S. Side-Chain Conformational Preferences Govern Protein-Protein Interactions. J. Am. Chem. Soc. 2016, 138, 10386–10389. (15) Fears, K. P.; Clark, T. D.; Petrovykh, D. Y. Residue-Dependent Adsorption of Model Oligopeptides On Gold. J. Am. Chem. Soc. 2013, 135, 15040–15052. (16) McUmber, A. C.; Randolph, T. W.; Schwartz, D. K. Electrostatic Interactions Influence Protein Adsorption (but Not Desorption) at the Silica-Aqueous Interface. J. Phys. Chem. Lett. 2015, 6, 2583–2587.
ACS Paragon Plus Environment
21
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1 2
(17) Liu, J.; Peng, C.; Yu, G.; Zhou, J. Molecular Simulation Study of Feruloyl Esterase Adsorption On Charged Surfaces: Effects of Surface Charge Density and Ionic Strength. Langmuir 2015, 31, 10751–10763.
3
(18) Peng, C.; Liu, J.; Xie, Y.; Zhou, J. Molecular Simulations of Cytochrome C Adsorption On Positively Charged
4
Surfaces: The Influence of Anion Type and Concentration. Phys. Chem. Chem. Phys. 2016, 18, 9979–9989.
5
(19) Romanowska, J.; Kokh, D. B.; Wade, R. C. When the Label Matters: Adsorption of Labeled and Unlabeled
6 7 8
Proteins On Charged Surfaces. Nano Lett. 2015, 15, 7508–7513. (20) Kubiak-Ossowska, K.; Mulheran, P. A. Mechanism of Hen Egg White Lysozyme Adsorption On a Charged Solid Surface. Langmuir 2010, 26, 15954–15965.
9
(21) Kubiak-Ossowska, K.; Jachimska, B.; Mulheran, P. A. How Negatively Charged Proteins Adsorb to Negatively
10
Charged Surfaces: A Molecular Dynamics Study of Bsa Adsorption On Silica. J. Phys. Chem. B 2016, 120,
11
10463–10468.
12 13 14 15 16 17 18 19
(22) Guigas, G.; Weiss, M. Effects of Protein Crowding On Membrane Systems. Bba.-Biomembranes 2016, 1858, 2441–2450. (23) Mittal, S.; Chowhan, R. K.; Singh, L. R. Macromolecular Crowding: Macromolecules Friend Or Foe. Bba.Gen. Subjects 2015, 1850, 1822-1831. (24) Sapir, L.; Harries, D. Is the Depletion Force Entropic? Molecular Crowding Beyond Steric Interactions. Curr. Opin. Colloid in 2015, 20, 3–10. (25) Shao, Q.; Jiang, S. Y. Molecular Understanding and Design of Zwitterionic Materials. Adv. Mater. 2015, 27, 15–26.
20
(26) Zheng, J.; Li, L. Y.; Tsao, H. K.; Sheng, Y. J.; Chen, S. F.; Jiang, S. Y. Strong Repulsive Forces Between
21
Protein and Oligo (Ethylene Glycol) Self-Assembled Monolayers: A Molecular Simulation Study. Biophys. J.
22
2005, 89, 158–166.
23 24 25 26 27 28
Page 22 of 31
(27) Sarkar, M.; Lu, J.; Pielak, G. J. Protein Crowder Charge and Protein Stability. Biochemistry-US 2014, 53, 1601–1606. (28) Linse, P.; Lobaskin, V. Electrostatic Attraction and Phase Separation in Solutions of Like-Charged Colloidal Particles. Phys. Rev. Lett. 1999, 83, 4208–4211. (29) Angelescu, D. G.; Linse, P. Monte Carlo Simulation of the Mean Force Between Two Like-Charged Macroions with Simple 1 : 3 Salt Added. Langmuir 2003, 19, 9661–9668.
ACS Paragon Plus Environment
22
Page 23 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
1
(30) Tiraferri, A.; Maroni, P.; Borkovec, M. Adsorption of Polyelectrolytes to Like-Charged Substrates Induced by
2
Multivalent Counterions as Exemplified by Poly(Styrene Sulfonate) and Silica. Phys. Chem. Chem. Phys.
3
2015, 17, 10348–10352.
4
(31) Luque-Caballero, G.; Martin-Molina, A.; Quesada-Perez, M. Polyelectrolyte Adsorption Onto Like-Charged
5
Surfaces Mediated by Trivalent Counterions: A Monte Carlo Simulation Study. J. Chem. Phys. 2014, 140,
6
174701.
7 8 9 10 11 12
(32) Wang, G.; Dong, X.; Sun, Y. Ion–Exchange Resins Greatly Facilitate Refolding of Like–Charged Proteins at High Concentrations. Biotechnol. Bioeng. 2011, 108, 1068–1077. (33) Yu, L. L.; Dong, X. Y.; Sun, Y. Ion-Exchange Resins Facilitate Like-Charged Protein Refolding: Effects of Porous Solid Phase Properties. J. Chromatogr. A 2012, 1225, 168–173. (34) Yang, C. Y.; Li, M.; Dong, X. Y.; Sun, Y. A Double-Modification Strategy for Enhancing Charge Density of Mono-Sized Beads for Facilitated Refolding of Like-Charged Protein. J. Chromatogr. A 2013, 1299, 85–93.
13
(35) Liu, H.; Dong, X.; Sun, Y. Grafting Iminodiacetic Acid On Silica Nanoparticles for Facilitated Refolding of
14
Like-Charged Protein and its Metal-Chelate Affinity Purification. J. Chromatogr. A 2016, 1429, 277–283.
15
(36) Liu, B.; Besseling, T. H.; Blaaderen, A. V.; Imhof, A. Confinement Induced Plastic Crystal-to-Crystal
16 17 18 19 20 21 22 23 24 25 26
Transitions in Rodlike Particles with Long-Ranged Repulsion. Phys. Rev. Lett. 2015, 115, 078301. (37) Akimov, A. V.; Prezhdo, O. V. Large-Scale Computations in Chemistry: A Bird's Eye View of a Vibrant Field. Chem. Rev. 2015, 115, 5797–5890. (38) Schröder, C.; Steinhauser, O.; Sasisanker, P.; Weingärtner, H. Orientational Alignment of Amyloidogenic Proteins in Pre-Aggregated Solutions. Phys. Rev. Lett. 2015, 114, 128101. (39) Bai, S.; Li, H.; Zhang, L. Standing Orientation of Lysozymes Induced by Electrostatically Repulsive Surface. Acta Phys.-Chim. Sin. 2013, 29, 849–857. (40) Zhang, L.; Lua, L.; Middelberg, A.; Sun, Y.; Connors, N. K. Biomolecular Engineering of Virus-Like Particles Aided by Computational Chemistry Methods. Chem. Soc. Rev. 2015, 44, 8608–8618. (41) Dijk, E. V.; Varilly, P.; Knowles, T. P.; Frenkel, D.; Abeln, S. Consistent Treatment of Hydrophobicity in Protein Lattice Models Accounts for Cold Denaturation. Phys. Rev. Lett. 2016, 116, 078101.
27
(42) Sayle, R.; Milnerwhite, E. Rasmol - Biomolecular Graphics for All. Trends Biochem. Sci. 1995, 20, 374–376.
28
(43) Dresselhaus, M. S.; Avouris, P. Introduction to Carbon Materials Research. Top. Appl. Phys. 2001, 80, 1–9.
ACS Paragon Plus Environment
23
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1 2 3 4
(44) Hummer, G.; Rasaiah, J. C.; Noworyta, J. P. Water Conduction through the Hydrophobic Channel of a Carbon Nanotube. Nature 2001, 414, 188–190. (45) Schuttelkopf, A. W.; van, A. D. Prodrg: A Tool for High-Throughput Crystallography of Protein-Ligand Complexes. Acta Crystallographica Section D-Biological Crystallography 2004, 60, 1355–1363.
5
(46) Chandrasekhar, I.; Kastenholz, M.; Lins, R. D.; Oostenbrink, C.; Schuler, L. D.; Tieleman, D. P.; van
6
Gunsteren, W. F. A Consistent Potential Energy Parameter Set for Lipids: Dipalmitoylphosphatidylcholine as a
7
Benchmark of the Gromos96 45a3 Force Field. Eur. Biophys. J. Biophy. 2003, 32, 67–77.
8 9 10 11 12 13 14 15 16 17
(47) Berendsen, H. J.; Vanderspoel, D.; Vandrunen, R. Gromacs - A Message-Passing Parallel Molecular-Dynamics Implementation. Comput. Phys. Commun. 1995, 91, 43–56. (48) Lindahl, E.; Hess, B.; van, S. D. Gromacs 3.0: A Package for Molecular Simulation and Trajectory Analysis. J. Mol. Model. 2001, 7, 306–317. (49) Mackerell, A. D. Empirical Force Fields for Biological Macromolecules: Overview and Issues. J. Comput. Chem. 2004, 25, 1584–1604. (50) Yu, G.; Liu, J.; Zhou, J. Mesoscopic Coarse-Grained Simulations of Hydrophobic Charge Induction Chromatography (Hcic) for Protein Purification. AIChE J. 2015, 61, 2035–2047. (51) Zhang, L.; Zhao, G. F.; Sun, Y. Effects of Ligand Density On Hydrophobic Charge Induction Chromatography: Molecular Dynamics Simulation. J. Phys. Chem. B 2010, 114, 2203–2211.
18
(52) Zhang, L.; Bai, S.; Sun, Y. Molecular Dynamics Simulation of the Effect of Ligand Homogeneity On Protein
19
Behavior in Hydrophobic Charge Induction Chromatography. Journal of Molecular Graphics and Modelling
20
2010, 28, 863–869.
21
(53) Zhang, L.; Zhao, G. F.; Sun, Y. Molecular Dynamics Simulation and Experimental Validation of the Effect of
22
Ph On Protein Desorption in Hydrophobic Charge Induction Chromatography. Mol. Simulat. 2010, 36, 1096–
23
1103.
24
(54) Zhang, L.; Sun, Y. Effect of Ligand Chain Length On Hydrophobic Charge Induction Chromatography
25
Revealed by Molecular Dynamics Simulations. Frontiers of Chemical Science and Engineering 2013, 7, 456–
26
463.
27 28
Page 24 of 31
(55) Yang, C. Y.; Yu, L. L.; Dong, X. Y.; Sun, Y. Mono-Sized Microspheres Modified with Poly(Ethylenimine) Facilitate the Refolding of Like-Charged Lysozyme. React. Funct. Polym. 2012, 72, 889–896.
ACS Paragon Plus Environment
24
Page 25 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
1 2 3 4 5 6 7 8
Langmuir
(56) Yang, C.; Dong, X.; Sun, Y. Mechanistic Studies of Protein Refolding Facilitated by Like-Charged Polymers. React. Funct. Polym. 2013, 73, 1405–1411. (57) Langenhof, M.; Leong, S.; Pattenden, L. K.; Middelberg, A. Controlled Oxidative Protein Refolding Using an Ion-Exchange Column. J. Chromatogr. A 2005, 1069, 195–201. (58) Machold, C.; Schlegl, R.; Buchinger, W.; Jungbauer, A. Matrix Assisted Refolding of Proteins by Ion Exchange Chromatography. J. Biotechnol. 2005, 117, 83–97. (59) Freydell, E. J.; van der Wielen, L.; Eppink, M.; Ottens, M. Ion-Exchange Chromatographic Protein Refolding. J. Chromatogr. A 2010, 1217, 7265–7274.
9
ACS Paragon Plus Environment
25
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Page 26 of 31
Table of contents only
ACS Paragon Plus Environment
26
Page 27 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
Figure 1 64x23mm (300 x 300 DPI)
ACS Paragon Plus Environment
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 2 134x102mm (300 x 300 DPI)
ACS Paragon Plus Environment
Page 28 of 31
Page 29 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
Figure 3 94x49mm (300 x 300 DPI)
ACS Paragon Plus Environment
Langmuir 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Figure 4 133x99mm (300 x 300 DPI)
ACS Paragon Plus Environment
Page 30 of 31
Page 31 of 31 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
Langmuir
TOC 56x37mm (300 x 300 DPI)
ACS Paragon Plus Environment