This Article Abstract Full Text (PDF) All Versions of this Article: biophysj.104.052787v1 88/2/939    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Lensink, M. F. Articles by Rosseneu, M. Search for Related Content PubMed PubMed Citation Articles by Lensink, M. F. Articles by Rosseneu, M.

Penetratin-Membrane Association: W48/R52/W56 Shield the Peptide from the Aqueous Phase

M. F. Lensink ^*, B. Christiaens ^*, J. Vandekerckhove ^* , A. Prochiantz  and M. Rosseneu ^*

^* Department of Lipoprotein Chemistry, Faculty of Medicine and Health Sciences, Ghent, Belgium; Department of Medical Protein Research, Flanders Interuniversity Institute for Biotechnology, Ghent, Belgium; and Development and Neuropharmacology, Ecole Normale Supérieure, Paris, France

Correspondence: Address reprint requests to M. F. Lensink, Service de Conformation de Macromolécules Biologiques et de Bioinformatique, Université Libre de Bruxelles, Boulevard du Triomphe – CP 263, B-1050, Brussels, Belgium. Tel.: 32-2-650-2013; Fax: 32-2-650-5425; E-mail: lensink{at}scmbb.ulb.ac.be.

	ABSTRACT

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Using molecular dynamics simulations, we studied the mode of association of the cell-penetrating peptide penetratin with both a neutral and a charged bilayer. The results show that the initial peptide-lipid association is a fast process driven by electrostatic interactions. The homogeneous distribution of positively charged residues along the axis of the helical peptide, and especially residues K46, R53, and K57, contribute to the association of the peptide with lipids. The bilayer enhances the stability of the penetratin helix. Oriented parallel to the lipid-water interface, the subsequent insertion of the peptide through the bilayer headgroups is significantly slower. The presence of negatively charged lipids considerably enhances peptide binding. Lateral side-chain motion creates an opening for the helix into the hydrophobic core of the membrane. The peptide aromatic residues form a -stacking cluster through W48/R52/W56 and F49/R53, protecting the peptide from the water phase. Interaction with the penetratin peptide has only limited effect on the overall membrane structure, as it affects mainly the conformation of the lipids which interact directly with the peptide. Charge matching locally increases the concentration of negatively charged lipids, lateral lipid diffusion locally decreases. Lipid disorder increases, through decreased order parameters of the lipids interacting with the penetratin side chains. Penetratin molecules at the membrane surface do not seem to aggregate.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Cell-penetrating peptides (CPPs) have been proposed as carriers for cellular uptake of various molecules, including proteins, oligonucleotides, and drugs. The homeodomain of Antennapedia is able to translocate through a cellular membrane (Joliot et al., 1991; Le Roux et al., 1993), mainly thanks to its third helix (Derossi et al., 1994, 1996), consisting of residues 43–58. The so-called pAntp peptide, or penetratin, was proposed as a vehicle for intracellular delivery of hydrophilic cargo molecules (Prochiantz, 1996, 1998; Pooga et al., 1998; Braun et al., 2002). Penetratin was the first member of an increasingly large family of natural or synthetic CPPs (Lindgren et al., 2000; Langel, 2002). Its uptake mechanism is to date not clear. Most membrane-associating peptides are amphipathic and translocate across membranes through either pore formation (Cornut et al., 1993; Bechinger, 1997; Matsuzaki, 1998) or in a potential-dependent manner (Maduke and Roise, 1993; Leenhouts et al., 1996). Charged residues in these peptides are mostly lysines. Molecular dynamics simulations of the N-terminal region of human surfactant protein-B (SP-B_1–25) showed partial diffusion of the N-terminal part of the peptide into the hydrophobic phase of a bilayer on a nanosecond timescale (Kaznessis et al., 2002). The membrane-interacting residues in SP-B_1–25 are hydrophobic, whereas the C-terminal charged residues remain within the lipid interface. The N-terminal end of penetratin is mostly hydrophobic, although not to the same extent. Its fluctuating helical hydrophobic moment (calculated over a seven-residue window), related to its DNA-binding ability, differentiates penetratin from the amphipathic helical peptide family (Thorén et al., 2000; Drin et al., 2001a; Binder and Lindblom, 2003a). Since cellular penetratin uptake occurs both at 37°C and 4°C, internalization through endocytosis was excluded (Derossi et al., 1994), whereas cellular import of reverse forms of various CPPs suggested that it is a receptor-independent process (Brugidou et al., 1995; Derossi et al., 1996; Mitchell et al., 2000; Wender et al., 2000). The nontoxic penetratin translocates across cellular membranes without membrane damage (Terrone et al., 2003), even when coupled to drugs (Hosotani et al., 2002), suggesting direct interaction with membrane lipids (Berlose et al., 1996; Derossi et al., 1996; Thorén et al., 2000; Terrone et al., 2003; Vivès et al., 2003). However, the mode of interaction and of lipid reorganization that allow for internalization without leakage is not known. Formation of inverted lipid micelles was suggested (Derossi et al., 1994, 1998; Prochiantz, 1996), implying that the peptide remains in an aqueous environment. Such a mechanism would facilitate the transport of hydrophilic compounds linked to the peptide. According to recent reports, penetratin internalization might at least partially occur through an energy-dependent process, involving endocytosis (Richard et al., 2003; Thorén et al., 2003; Terrone et al., 2003; Console et al., 2003; Vivès et al., 2003). On the other hand, the uptake of other arginine-rich peptides does not seem to be affected by endocytosis inhibitors (Suzuki et al., 2002; Console et al., 2003), suggesting the existence of different competing pathways for CPPs (Hällbrink et al., 2001; Vivès et al., 2003).

The penetratin sequence (RQIKIWFQNRRMKWKK) consists of 16 residues out of which four lysines and three arginines. Electrostatic interactions are therefore likely to play a key role in the association process (Christiaens et al., 2002). However, no single residue critical for membrane interaction could be identified (Drin et al., 2001a). Negatively charged lipids promote the transfer of penetratin from a hydrophilic to a hydrophobic environment (Dom et al., 2003), and electrostatic effects probably contribute not only to binding, but further to peptide translocation (Binder and Lindblom, 2003a). Penetratin is not sufficiently hydrophobic to insert deeply into phospholipid model membranes (Drin et al., 2001a; Brattwall et al., 2003). This suggests that charge neutralization is required for a deeper insertion of the peptide into the hydrophobic core of the membrane (Dom et al., 2003; Console et al., 2003). The contribution of the hydrophobic residues, and especially of the aromatic tryptophan and phenylalanine residues is probably crucial for internalization. Mutation of either tryptophan decreases internalization, whereas double substitution completely inhibits peptide internalization (Derossi et al., 1994; Fischer et al., 2000; Drin et al., 2001b; Christiaens et al., 2002; Lindberg et al., 2003).

We selected molecular dynamics (MD) simulations as a method to investigate the atomic interactions that lie at the basis of the penetratin-membrane association and internalization. MD is a well-established methodology for the study of the dynamics of proteins (Berendsen, 1996; Karplus, 2002), membranes (Tieleman et al., 1997), and peptide-membrane interactions (La Rocca et al., 1999; Shepherd et al., 2003) that enables a direct observation of dynamical events (Berendsen, 2001). The penetratin-membrane association was studied by careful analysis of the obtained trajectories. We examined the residues responsible for the peptide-lipid binding, the local environment around the peptide, the behavior of peptide aromatic residues, possible peptide aggregation, the effect of K/RA mutations, and the penetratin structure in solution. The simulation results are summarized under the Discussion. The computational details are described in the following section.

	COMPUTATIONAL DETAILS

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Molecular dynamics
MD simulations were performed with the Gromacs 3.1.4 package (Berendsen et al., 1995; Lindahl et al., 2001), using the Gromos96 43a2 force field (Van Gunsteren et al., 1996), with an extension for the lipid phosphate and acyl chain parameters (Berger et al., 1997). The peptide structure was taken by cutting residues 43–58 (the 3rd helix) from the NMR structure of the Antennapedia homeodomain (pdb-code 1ahd (Billeter et al., 1993), model 1). The C-terminus was capped with an amine group. Peptide-membrane simulations were prepared by placing the peptide horizontally on top of an equilibrated bilayer box, at a distance of 2 nm from the bilayer surface; the axis normal to the bilayer plane was extended to include the peptide and the formed vacuum was filled with water molecules, resulting in 4500 water molecules in a box 6 x 6 x 8.5 nm³. Single point charge water was used (Berendsen et al., 1981). An equilibrated bilayer of 128 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) molecules was taken as the starting point (Tieleman et al., 1999), with at random 8 POPC molecules in each layer exchanged for the negatively charged alternatives 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphate (POPA), 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-rac-(1-glycerol)] (POPG), or 1-palmitoyl-2-oleoyl-sn-glycero-3-[phospho-L-serine] (POPS). Coordinates of the phosphate groups and acyl chains were kept in this process, adding water molecules in a formed vacuum (POPA) or deleting them upon overlap with the new headgroup function. All systems, including the bilayer systems without the peptide yet present, were made electrostatically neutral by adding the required amount of counterions Na⁺ or Cl^– and then submitted to an energy minimization, 10 ps position restraint MD and 2 ns free MD before the production run was started. Weak coupling to a temperature (300 K) and pressure (isotropic, 1.0 bar) bath was employed (Berendsen et al., 1984), using coupling constants of 0.1 and 1.0 ps, respectively. Coulomb interactions were treated with fast particle mesh Ewald (Essman et al., 1995), using a grid spacing of 0.12 nm and fourth-order interpolation. A spherical cutoff of 1.0 nm was applied to the Van der Waals interactions. Bond lengths were constrained using the LINCS algorithm (Hess et al., 1997). Equations of motion for the water atoms were solved analytically with the SETTLE algorithm (Miyamoto and Kollman, 1992). Dummy atoms were used to limit high-frequency vibrations involving hydrogen atoms (Feenstra et al., 1999). A time step of 4 fs was employed, removing center of mass motion every step and updating the neighbor list every 10 steps. Table 1 shows a list of the performed simulations. The simulations were performed on our own linux cluster consisting of five dual-CPU machines (eight Athlon 2000+ (AMD, Sunnyvale, CA) and two Xeon processors (Intel, Santa Clara, CA)).

Analyses
The axis of the penetratin helix is determined using a rotational least-squares fitting method mapping the C's of the helix onto itself, but one residue out of phase, i.e., residue i is mapped onto residue i + 1 (a screw-transform superimposes the two helices) (Christopher et al., 1996). A quaternion-based method is used to find the helical axis as the axis of rotation necessary to superimpose these atoms (Mackay, 1984).

The orientation of the aromatic residues is described in terms of two order parameters (Tieleman et al., 1998). S_N is defined as with the angle between the bilayer normal and the normal to the plane of the aromatic ring. S_L is defined in the same way but with the angle between the bilayer normal and the vector from C to C or C₂ for Phe and Trp, respectively. For a more in-depth description, we refer to Fig. 2 of Tieleman et al. (1998).

View larger version (42K): [in this window] [in a new window]   FIGURE 2  Minimal distance between the hydrogen bond donor atoms of the penetratin side chains and any of the phosphorus atoms of the lipid bilayer.

View larger version (22K): [in this window] [in a new window]   FIGURE 5  Atom density profile, calculated over a time frame of 100–230 ns, as a function of distance in the z-coordinate to the bilayer center. Density profiles for both POPG (A) and POPC (B) lipids are plotted. Solid lines indicate the peptide-containing layer, dashed lines are for the layer without peptide, mirrored in the center of the bilayer. Color coding is black for headgroup functions (choline or glycerol), red for phosphate, green for carbonyl, and blue for double-bond groups. Penetratin density is indicated with the cyan solid line, scaled (divided by 7) in A.

View larger version (50K): [in this window] [in a new window]   FIGURE 6  (A) Number of lipids interacting with penetratin, as defined in Computational Details. (B) Average distance to the bilayer center (in the z-coordinate), for penetratin and selected lipid subgroups. Color coding corresponds to Fig. 5, with interacting phosphate, carbonyl, and double-bond groups in brown, dark green, and magenta, respectively.

The programs to calculate 1), the axis of a helix, 2), the orientation of aromatic residues, and 3), lipids that are under a peptide have been written by the author (M.F.L.) and are available upon request.

	RESULTS

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
A number of simulations were performed, as listed in Table 1. The long simulations P-PC and P-PG, with neutral and 12.5% negatively charged bilayer were used for all major analyses. The P-PC and P-PC^ß simulations were used to test the peptide-membrane association for a neutral bilayer, using a different initial orientation of the peptide (rotated 60 and 120 degrees about the axis of the helix). The P-PA and P-PS simulations were used to investigate the influence a different negatively charged lipid has on the overall peptide-lipid interactions. As no significant difference in orientation was found during membrane binding with respect to the reference P-PC and P-PG simulations, these simulations show the same characteristics as their P-PC and P-PG counterparts, unless stated otherwise. Moreover, as the results from both these simulations were also similar, all figures refer to the P-PG simulations, unless stated otherwise in the figure caption. Once the peptide was well docked in the bilayer, we took a key frame of the P-PG simulation and performed additional simulations: P_{52A, 55A}-PG and P_{53A, 57A}-PG are simulations with R52A/K55A and R53A/K57A mutations, respectively. Aggregation of penetratin was investigated by adding additional penetratin peptides, yielding simulations P₂-PG and P₃-PG.

Penetratin-membrane association
Fig. 1 shows the approach of penetratin toward the membrane surface, as the z-component of the distance between the average coordinates of the peptide and the phosphorus atoms of the half-side lipid bilayer, that is closest to the peptide. The P-PG simulation (Table 1) in Fig. 1 A shows a rapid penetratin-membrane association that takes place during the equilibration phase of the simulation, as the mean distance decreases from 2 to <1 nm within 2 ns. It takes about a third of the remaining 230 ns of the simulation to reach an equilibrium between peptide and membrane (Fig. 1, B and C), when penetratin has docked between the hydrophilic headgroups of the lipids. The concomitant secondary structure evolution of the peptide (data not shown) shows that penetratin forms a stable -helix when bound to the lipid bilayer.

View larger version (36K): [in this window] [in a new window]   FIGURE 1  Penetratin-membrane distance and interaction energy for the P-PG (A–C) and P-PC (D–F) simulations. (A, B, D, and E) Distance (z-coordinate) between the mean coordinates of penetratin and the phosphorus atoms of the lipid bilayer half interacting with the peptide. Data are given for both the equilibration (A and D) and production runs (B and E), separated at t = 0 ns. (C and F) Peptide-membrane interaction energy, split into the Coulomb and Lennard-Jones contributions.

Fig. 1, D and E, show the penetratin-membrane distance for the simulation of the peptide association with a neutral POPC membrane (P-PC). Peptide-membrane association is significantly slower in this case as binding is 40 ns longer than the equilibration phase, compared to P-PG. Interaction energies are similar to those in Fig. 1 C, but equilibrium is reached only after 160 ns. A slightly more favorable electrostatic energy component, and a deeper positioning of the peptide in the bilayer, is detected. The structure of the peptide bound to a neutral bilayer remains -helical.

To characterize the peptide orientation during its approach toward the bilayer, we plotted the minimal distance between the nitrogen atoms of the seven positively charged penetratin residues and any of the lipid phosphorus atoms (Fig. 2). This figure identifies residues K46, R53, K57, and K58 as those accounting for the mutual electrostatic binding. The first three residues are located on the same side of the helical peptide (see also Fig. 3 B); K58 is located at the C-terminal end of the peptide and has therefore higher mobility. Fig. 3 A shows binding of penetratin to a phospholipid bilayer, from left to right in three concatenated snapshots. Whereas the left-hand image shows the initial structure before equilibration, the center image—taken at the end of the equilibration phase—shows K46, R53, and K57 directed toward the bilayer surface, making hydrogen bonds and salt bridges with the lipids. The right-hand image shows the position of penetratin in the lipid bilayer, as recorded toward the end of the simulation. Throughout this process, penetratin remains horizontal, as shown by the angle between the axis of the helix and the normal to the bilayer plane (Fig. 4 A), with the positively charged residues K46, R53, and K57 directed toward the membrane surface. Fig. 4 A shows that during practically the entire P-PG simulation, penetratin remains parallel, or slightly tilted, to the membrane surface. The interaction of penetratin with a neutral lipid membrane occurs through the C-terminal end of the peptide, and during the initial binding penetratin is perpendicular to the bilayer surface (Fig. 4 C). Subsequent rotations perpendicular to and about the helix axis enable interactions of the same residues K46, R53, and K57 with lipids. Once the peptide is settled in the membrane, the longitudinal vector (the helix axis) assumes a value of 80° (90° for P-PC). A steady value of 75° (Fig. 4 B) is found (60° for P-PC, Fig. 4 D) for the angle between the z-axis and the lateral vector, excluding further rotation about the helix axis. Fig. 3 B shows the position of the helix in the membrane, with the z-axis as in Fig. 3 A.

View larger version (55K): [in this window] [in a new window]   FIGURE 3  (A) Binding of penetratin to a lipid bilayer. Snapshots are taken from the P-PG simulation at t = –1.85 ns, t = –50 ps, and t = 200 ns (concatenated left to right); only peptide and membrane are plotted. The arrow indicates the z axis and the normal to the bilayer plane. Membrane lipids are drawn in wireframe, the peptide, charged and aromatic residues only, with ball-and-sticks. Color coding for lipids: P, violet; O, red; N, blue; other, gray; color coding for the peptide: Lys and Arg, blue; Trp and Phe, orange. The figure was prepared with MOLSCRIPT (Kraulis, 1991). (B) Edmundson wheel representation of penetratin. The figure is plotted such that the orientation in the z-coordinate corresponds to A, i.e., the C of Q50 is located deepest in the membrane.

View larger version (28K): [in this window] [in a new window]   FIGURE 4  Alignment of penetratin with the bilayer surface during the P-PG (A and B) and P-PC (C and D) simulations. (A and C) Longitudinal alignment: angle between the axis of the helix (NC) and the normal of the bilayer plane. (B and D) Lateral alignment: angle between the vector from F49C to the geometric center of I47C and N51C, and the normal of the bilayer plane.

Influence on the bilayer structure
Fig. 5 shows the density profile for the lipid bilayer halves in the presence and absence of penetratin, averaged over a time frame of 100–230 ns. The density peaks are broader and deeper for the peptide-containing layer (solid lines), for both lipid headgroup function, phosphate, carbonyl, and double-bond groups (for definitions see Computational Details), due to the presence of penetratin, docked between the lipid headgroups. The overall bilayer structure shows little change, as the peak maxima are not significantly shifted, but peaks are broadened as some lipids are pushed toward the core of the membrane by the bound penetratin. This applies in particular to the lower end of the phosphate and carbonyl group densities of POPG. The peak broadening is more pronounced for the POPG glycerol group, which is also clearly shifted, than for the POPC choline group. The POPG headgroups (glycerol and phosphate), including those which are not in the peptide-containing layer, are set deeper in the membrane core than the neutral POPC lipids.

Fig. 6 B shows a plot of the average mass-weighted distance (z-coordinate) to the membrane center of selected groups of the bilayer half containing the peptide, and of penetratin. Subtraction of the red phosphate line from the turquoise penetratin line yields Fig. 1 B. Fig. 6 A represents the number of lipids interacting with the peptide as defined by our criteria in Computational Details. This number is fluctuating around 15, corresponding to 25% of the bilayer half, including lipids whose acyl chains are partly under the peptide, but whose headgroups are beside it. During the first 10 ns of the simulation, when penetratin is still at an average distance of >5 Å from the bilayer surface, the interacting lipids are at roughly the same distance from the membrane center as the entire layer. Penetratin binding to the bilayer surface immediately affects the lipids under the peptide, as they are significantly pushed down toward the lipid core. This applies not only to the groups that directly interact with the peptide, i.e., the phosphate or carbonyl groups, but also to the double bonds further away.

The conformation of lipids is usually described by deuterium order parameters S_CD, measured by ²H-NMR experiments, or calculated from the lipid acyl chain C–C dihedral angles (Merz and Roux, 1996). Fig. 7 shows the order parameters for the sn2 and sn1 acyl chains (Fig. 7, A and B, respectively). The shape of the order parameter curves resembles that obtained from ²H-NMR experiments, with a plateau near the carbonyl groups and increasing random structure toward the bilayer core, including the zig-zag appearance for the unsaturated oleoyl chain and the dip due to the sp₂ instead of sp₃ hybridization around the double bond (Seelig and Seelig, 1980). The order parameters calculated for the peptide-interacting lipids (dashed lines and open diamonds) are shifted downward over the entire range of the acyl chain, showing that the presence of penetratin introduces disorder over the entire length of the lipids.

View larger version (13K): [in this window] [in a new window]   FIGURE 7  Calculated deuterium order parameters for the lipid acyl chains. (A) sn2 Oleoyl chain. (B) sn1 Palmitoyl chain. Solid lines and circles are order parameters calculated over all lipids, dashed lines and open diamonds denote the order parameters for the lipids that are interacting with penetratin. Order parameters were calculated over the entire simulation length.

View larger version (59K): [in this window] [in a new window]   FIGURE 8  (A) Scatter plot showing the diffusion of the lipid phosphate groups in the peptide-binding membrane layer. The (x, y) coordinates are plotted, with black and blue dots indicating a POPC and POPG phosphorus atom, respectively, and red dots a penetratin C atom. Coordinates are collected over the time frame 100–230 ns, with an increment of 100 ps. All coordinates are given relative to penetratin; during the simulation the lateral diffusion constant D of penetratin is (1.8 ± 8.3) x 10–11 m2/s. This value, although taken from the simulation of a single penetratin molecule, agrees with experimental values (Andersson et al., 2004). Penetratin orientation in this figure, is with R43 at the lower left and the associating residues K46, R53, and K57 to the right, corresponding to a top view of Fig. 3 B. (B) Root mean-square fluctuation of penetratin side chains, in solution (thin-dashed bars), during initial binding to the bilayer (thick-dashed bars), and in complex with phospholipids (solid bars). The RMSF values are calculated from the P-WAT simulation, from the third added penetratin in the P3-PG simulation, and from the P-PG simulation, respectively.

View larger version (51K): [in this window] [in a new window]   FIGURE 9  (A, B, and C) Orientation and position of the aromatic residues W48 (top), F49 (center), and W56 (bottom). (A and B) Orientational parameters SN and SL, respectively. (C) Average distance (z-coordinate) between the aromatic ring systems and the phosphate groups of peptide-associated lipids. (D) Distance of the cation- interactions between aromatic residues W48, F49, W56 and either R52 or R53. (Top) W48-R52; (center) F49-R53; (bottom) W56-R52. The distance is measured from the geometric center of the two arginine N atoms to the center of the aromatic 6-ring group.

Association between penetratin peptides at the bilayer surface
A second penetratin molecule was added to the P-PG simulation at t = 100 ns, yielding simulation P₂-PG, to investigate any association of penetratin molecules docked to the lipid membrane. After 25 ns of simulation, when the second peptide was docked to the bilayer, a third penetratin was added, yielding simulation P₃-PG.

Binding of the second penetratin molecule occurred onto the bilayer side where no peptide was yet present (we apply periodic boundary conditions), with the same residues K46, R53, and K57 (and K58) accounting for initial peptide-lipid binding. However, binding of a third penetratin was significantly slower, resembling the P-PC simulation, due to prior neutralization of POPG negative charges by the bound penetratin peptides. The third peptide associates first with the lipids through its C-terminal end, perpendicular to the lipid bilayer, as observed for the neutral bilayer P-PC simulation. This binding occurs at a distance from the first and second peptides, and no peptide-peptide interaction is observed.

Binding of penetratin variants to the lipid bilayers
Fig. 3 B shows that the pairs of charged residues R52/K55 and R53/K57 are located at different sides of the helical peptide. Double mutants were synthesized to study the contributions of these residues to penetratin interactions with model membranes and cells (Christiaens et al., 2004). We looked at the effect of these mutations upon penetratin interactions with the lipid bilayer, by applying the R52A/K55A and R53A/K57A mutations to the P-PG simulation at t = 100 ns, and pursuing the simulations for 20 ns (simulations P_{52A, 55A}-PG and P_{53A, 57A}-PG).

The structure, position, and orientation of the peptides were not affected by the mutations. The -helical structure was retained, and the lipid structure did not change significantly, in agreement with the experimental data (Christiaens et al., 2004).

Table 2 shows the interaction energies of the mutant peptides both with the membrane, expressed as separate electrostatic (Coulomb) and hydrophobic (Lennard-Jones) energies, and with the aqueous phase. The R/KA mutations have little effect on the LJ energy, whereas the Coulomb energy was strongly affected. Peptide-membrane interactions increase for the R52A/K55A mutant, whereas they decrease for the R53A/K57A mutant. The peptide-solvent interaction decreases slightly for the R53A/K57A mutant, and more strongly for the R52A/K55A mutant.

Penetratin structure in water
The original penetratin peptide is -helical, as it originates from the third helix of the Antennapedia homeodomain. Circular dichroism measurements show that penetratin in buffer is mostly random, and that its -helical content increases upon binding to phospholipid vesicles (Persson et al., 2001, 2003; Magzoub et al., 2003; Christiaens et al., 2004). We performed a 120-ns simulation of the helical penetratin to investigate its stability in aqueous solution. The peptide assumes a random conformation in water, in agreement with the experimental observations. The -helical structure remains stable enough during the beginning of the simulation (10 ns) to prevent unfolding during the bilayer simulations, before lipid association. Moreover, our P-PC simulation shows that the presence of a neutral bilayer has a stabilizing effect on the helix, preventing unfolding of the peptide during the first 40 ns, before binding to the bilayer (Fig. 1 E). The simulation of the penetratin structure in an aqueous solution, using a cluster analysis based on backbone root mean-square deviation, shows a largest-populated cluster of structures with an unfolded C-terminus. W48, W56, and R53 form a -cluster through R53/W56 cation- and R53/W48 NH- interactions. The other residues are all solvent-directed.

The association with phospholipids significantly lowers the freedom of motion of the peptide residues. Fig. 8 B shows the root mean-square fluctuation (RMSF) of the individual residues for the peptide in solution (thin dashed bars) and for the P-PG simulation (solid bars). For the peptide in water, the lower RMSF on the N-terminal side suggests that this end of the peptide remains -helical, as is confirmed by the secondary structure evolution (data not shown). In contrast, when bound to a membrane, the peptide retains little flexibility, even for the terminal residues R43 and K58. The flexibility of a peptide approaching the membrane, such as the third peptide in the P₃-PG simulation, is intermediate (thick dashed bars).

	SUMMARY AND DISCUSSION

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
In this computational approach we showed that negatively charged lipids significantly increase the speed of association between penetratin and a lipid membrane, as peptide-lipid association occurs within a few nanoseconds for a partially negatively charged membrane, whereas it is slower by an order of magnitude for a neutral bilayer. Penetratin approaches a charged lipid bilayer horizontally, compared to a perpendicular approach toward a neutral bilayer. Peptide-lipid association occurs through formation of salt bridges between the positively charged residues K46, R53, and K57 and the lipid phosphate groups. Tryptophan fluorescence studies previously showed the importance of peptide positively charged residues for the initial binding to negatively charged vesicles, since double R/KA mutations involving the residues K46/R52/R53/K55/K57 significantly decreased the binding affinity (Christiaens et al., 2002, 2004). Our results, which identify three of these residues as the ones responsible for initial binding, account for both the lower binding affinity of mutants including double mutations of binding residues, and for the relatively higher binding affinity for the R52A/K55A variant. Attraction between penetratin and bilayer lipids is bidirectional, i.e., penetratin is not only attracted toward the membrane, but a few lipids are further partially pulled out of the bilayer. However, as the displacement of a lipid out of a bilayer is an energetically unfavorable process (Cevc and Marsh, 1987; Marrink et al., 1998), and established hydrogen bonds are maintained, penetratin comes close to the bilayer. No preference for negatively charged phospholipids POPG, POPA, and POPS was found during initial binding.

After initial binding, penetratin is oriented parallel to the bilayer surface. Subsequent docking of the peptide between the lipid headgroups involves only side chain movement. This process takes 100 ns for a (negatively) charged bilayer and up to 200 ns for a neutral bilayer. In peptide uptake experiments, using a mixed POPC/POPG bilayer, an increase in neutral/charged vesicles ratio from 1:1 to 3:1 is found to decrease internalization by a factor of 2 (Terrone et al., 2003). This is in agreement with our simulations, where the absence of a negatively charged phospholipid headgroup slows down peptide insertion by increased repulsion with penetratin. In addition, this repulsion is also responsible for the deeper insertion of penetratin in a neutral bilayer, an effect also observed in the binding of magainin to lipid membranes (Wieprecht et al., 1999). Kinetics of binding and insertion decrease as POPA > POPG > POPS > POPC.

The peptide remains fully -helical during the entire process, with the Trp residues oriented toward the solvent. No rotation around the helical axis occurs during docking. The abundance of aromatic residues in membrane proteins has long been observed (Von Heijne, 1997), together with their propensity to reside at membrane interfaces (Yau et al., 1998). A role of Trp as translocation determinant of peptides has been proposed (Schiffer et al., 1992), and mutation of both tryptophans in penetratin was found to abolish internalization (Derossi et al., 1994). Trp fluorescence quenching and oriented circular dichroism studies showed that Trp residues in penetratin remain close to the water-lipid interface (Magzoub et al., 2003; Christiaens et al., 2004). Although the preference of tryptophan for a membrane interface is generally attributed to an amphipathic and dipolar interaction, its flat and rigid shape limits its access to the hydrophobic core of the membrane, whereas its -electronic structure favors its position at the interface (Yau et al., 1998). Previous studies have suggested a solvent-exposed role for both tryptophans in penetratin (Fragneto et al., 2000; Magzoub et al., 2002). Our results show that both tryptophans lie preferentially at the interface, in a slightly tilted orientation, acting as a shield between the peptide and the aqueous phase above. R52 forms mutually exclusive -stacking interactions with W48 (cation-) and W56 (NH- hydrogen bond), with a strong preference for W48. A NH- hydrogen bond interaction exists also between F49 and R53. When R53 forms hydrogen bonds with the phosphate groups of two lipids, F49 becomes part of the W48-R52 cluster through CH- hydrogen bonds to W48. The neighboring residues W48 and F49 were previously proposed to belong to a structural motif contributing to peptide internalization (Le Roux et al., 1993). This hypothesis is supported by our simulations showing that these residues can assume a T-stacking conformation. However, the W48/R52 and F49/R53 pairs seem more relevant, and thus might account for decreased internalization of mutants lacking either R52 or R53 (Christiaens et al., 2004).

Cellular uptake experiments previously demonstrated the crucial role of arginines in cell-penetrating peptides (Mitchell et al., 2000; Thorén et al., 2003). Our results further stress the importance of arginines R52 and R53 for penetratin-membrane interactions, with the latter being mainly responsible for the initial electrostatic binding, where it can form bidentate hydrogen bonds with lipid phosphate groups, whereas the former contributes to the subsequent insertion of the peptide into the lipid bilayer and possibly to its translocation through its cation- interaction with both tryptophans. Although quantum-mechanical in nature, cation- interactions can be described using classical mechanical force fields (Donini and Weaver, 1998). Cation- interactions, especially between arginine and tryptophan, significantly contribute to the structure and function of biomolecules (Ma and Dougherty, 1997; Gallivan and Dougherty, 1999). In most proteins, aromatic residues forming cation- interactions are found at the protein surface (Flocco and Mowbray, 1994). This is also observed for the third helix of the Antennapedia homeodomain (Billeter et al., 1993), i.e., penetratin, where W48 and R52 are in a stacking conformation between the protein and DNA. A parallel configuration of these residues, with the NH₂ groups of arginine still available for hydrogen bonding, is energetically favorable (Mitchell et al., 1994; Gallivan and Dougherty, 1999), in agreement with the results of our calculations for the R52-W48 combination. The lower probability for the R52-W56 stacking might be due to destabilizing effects of the three lysines K55, K57, and K58 surrounding W56. In conclusion, we find that the aromatic residues do not contribute to the initial binding, but rather to the subsequent insertion of the peptide between the bilayer headgroups, when they shield the peptide from the aqueous phase. As W48, F49, and W56 behave as aromatic -acceptors, whereas R52 and R53 act as donors, this shield is stretched over the entire length of the peptide through cation- interactions and NH- hydrogen bonding. The unique combination of hydrophobicity with (weak) hydrogen-bond donor and acceptor capacities might be the key parameter for the transition from a hydrophilic to a hydrophobic phase. A WF mutation would still be functional, though less efficient. The W48F mutation would shift the cation- interaction from W48/R52 to W56/R52, i.e., from the center of the peptide toward its C-terminal end. In contrast, the W56F mutation strengthens the W48/R52 interaction, as F56 would be a less favorable alternative, in agreement with the crucial role of W48 for the internalization process (Derossi et al., 1994).

The preferred orientation of penetratin is almost parallel to the bilayer, at a slight angle (from C to N) of 80–90° with the normal to the bilayer surface, as observed before (Magzoub et al., 2003), and as could be expected for cationic peptides (Zhang et al., 2001). Under this angle, W48 inserts deeper into the bilayer than W56, in agreement with previous reports (Lindberg et al., 2003; Salamon et al., 2003). Trp fluorescence-quenching measurements supported a deeper insertion of the W56F compared to the W48F variant, whereas in the R52A/K55A and R53A/K57A mutants Trp residues inserted deeper into the lipids compared to wild-type penetratin (Christiaens et al., 2004). The R52A/K55A double mutation has a twofold effect: removal of the W48/R52 cation- interaction favors insertion of both tryptophan residues into the membrane hydrophobic phase, whereas the repulsion of penetratin and water is increased by the mutation of the solvent-directed residues R52 and K55 to an alanine. In the R53A/K57A double mutant, where two of the positively charged residues interacting with the lipid phosphate headgroups are mutated to an alanine, binding might occur through residues R52 and K55. However, the deeper insertion that was found for the R52A/K55A compared to the R53A/K57A mutant (Christiaens et al., 2004) suggests similar initial peptide-lipid binding, where both tryptophans remain solvent-accessible. Such a mode of binding might involve residues R43/K46/Q50. In both double mutants, decreased interaction of W56 with the surrounding lysines promotes its insertion into the hydrophobic membrane phase. In a neutral membrane (simulations P-PC and P₃-PG) we find that W56 inserts deeper than W48, whereas the peptide remains parallel to the membrane surface. This is possible due to an increased freedom of motion for W56 as it is surrounded by positively charged residues that might otherwise interact with negatively charged lipids, preventing passage of the solvent-accessible W56 into the lipid bilayer core.

When the peptide docks between the lipid headgroups, the positively charged residues K46, R53, and K57, which were responsible for initial binding, do not penetrate further into the bilayer, but rather move sideways, together with Q50 and N51. The lipid-peptide hydrogen bonds remain intact, lipid headgroups move along, and the hydrophobic membrane core is exposed to the helical peptide surface. Upon deeper insertion into the membrane core, the charge of R52 might be compensated through a bidentate hydrogen bond with an additional lipid phosphate group, thereby enhancing the hydrophobicity of the surrounding tryptophan residues, and of the peptide as a whole. The lipid phase is condensed, and lipids are pushed down, with shifted and broader density peaks compared to the inner membrane. This effect is more significant for the negatively charged lipids, than for POPC. With respect to the inner membrane, positively charged POPC choline groups are directed outward, whereas their neutral POPG glycerol counterparts are directed inward. The presence of penetratin decreases the order of the lipid acyl chains in contact with the peptide, as deuterium order parameter curves are shifted down. The differences are short-ranged, and lipids that are not in the immediate vicinity of penetratin are not affected by its presence. Association with penetratin, moreover, stabilizes the membrane bilayer, through hydrogen bonding with nonsolvent-directed penetratin side chains. A study of penetratin diffusion showed that a negatively charged lipid surface, but not a neutral one, restricted motional flexibility (Andersson et al., 2004). The presence of negatively charged lipids around penetratin supports a mechanism involving charge compensation, whereas lipid association around the peptide might account for the lack of calcein leakage upon peptide translocation (Christiaens et al., 2004).

The complete internalization process occurs on a millisecond timescale (Hällbrink et al., 2001; Terrone et al., 2003), and is therefore out of range of conventional MD simulations. Moreover, although internalization does not need a transmembrane potential (Thorén et al., 2000), such a potential is required to obtain a physiologically relevant process (Drin et al., 2001a; Kramer and Wunderli-Allenspach, 2003; Terrone et al., 2003). We are currently investigating the effect a transbilayer potential has on the structure and positioning of penetratin in a bilayer.

Although these results clearly identify the charged residues that lie at the basis of the peptide-membrane association, and hence confirm the importance of electrostatics for initial binding, there are additional effects that play a role as well. The classical hydrophobic effect can be understood as the release of the hydration shell around the peptide upon membrane incorporation. This is an entropy-driven process (Tanford, 1973) that is strongly dependent on temperature (Gill and Wadsö, 1976) and as such was found not to be the dominating contribution to penetratin binding (Binder and Lindblom, 2003b). The classical hydrophobic effect should not be confused with the so-called nonclassical hydrophobic effect (Seelig and Ganz, 1991), where an exothermic binding heat is found due to favorable peptide-lipid and lipid-lipid interactions (Wieprecht et al., 1999). Our simulations show an alteration in the conformation of the lipid headgroups and an increase in the lipid packing density. Both effects have also been observed using surface plasmon resonance and impedance measurements (Binder and Lindblom, 2003b). A third, enthalpically favorable, effect is the transition to a more ordered (-helical or ß-sheet) conformation upon membrane binding (White and Wimley, 1998). According to our simulations, penetratin is mainly unfolded in water, whereas it becomes 60% -helical when associated with lipids. This increased stability of the -helical structure is due to both backbone hydrogen bond shielding at the lipid/water interface by R52, W48, and W56, and to the decreased dielectric constant in the hydrophobic membrane core. These effects lower the hydration of backbone NH and CO groups, resulting in a more stable and less hydrophilic helical structure (Avbelj et al., 2000; García and Sanbonmatsu, 2002; Roccatano et al., 2002).

We found no aggregation of penetratin peptides at the membrane surface, as reported in solution (Magzoub et al., 2002; Christiaens et al., 2004), suggesting that the negatively charged bilayer does not induce cooperative binding. A theoretical study of protein adsorption on a mixed membrane shows that the negatively charged lipids in the lipid bilayer adjust their local concentration to achieve optimal charge matching with the protein or peptide (May et al., 2000). Subsequent binding of additional peptides in the vicinity of the first one is significantly less favorable due to lateral repulsive interactions between the adsorbed peptides and the effective charge reversal as solvent-directed positive charges of the peptide are not matched by negatively charged lipids. A penetratin peptide, when bound to a bilayer, interacts with 15–20 lipids, suggesting that above a molar peptide:lipid ratio of 0.05 free peptides remain in the aqueous phase, as observed for the circular dichroism spectra of penetratin associated with 70:30 PC/PS vesicles (Christiaens et al., 2004). The lack of cooperation and aggregation reported here does not exclude a cooperative effect in the final phase of the internalization, as only the initial steps of this process were studied. Further internalization might involve cooperation between several penetratin peptides, to reduce surface tension and enable cell-penetrating peptides, and its cargo, to enter the membrane. Membrane thinning was shown for alamethicin (He et al., 1996) and many other peptide-lipid systems (Ludtke et al., 1995; Tieleman et al., 1998; Heller et al., 2000; Chen et al., 2003); at low concentrations, the peptides adsorb onto the bilayer surface, whereas above a critical lipid-dependent concentration, a fraction of the peptides insert into the membrane. Penetratin might induce membrane thinning when entering the hydrophobic inner bilayer and dragging phospholipids along. Motion of the phospholipid headgroups together with penetratin secures the current phase separation consisting of lipid headgroups that could tightly encompass any hydrophilic cargo, resulting in minimal leakage. This mechanism would include flip-flop of certain phospholipids (Matsuzaki et al., 1996), but excludes the formation of pores as with some antimicrobial peptides (Zasloff, 2002). Our simulations showed favorable interactions between penetratin and the hydrophobic core of the bilayer, which together with hydrogen bonding to negatively charged lipids, might destabilize the phospholipid bilayer, finally resulting in membrane translocation of the penetratin peptide.

	ACKNOWLEDGEMENTS

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
M.F.L. acknowledges support from European Economic Community contract NPRN-CT-2001-0242 "Messenger proteins: mechanisms of action and biological significance". B.C. was supported by a grant from the Institute for the Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen).

Submitted on September 14, 2004; accepted for publication November 2, 2004.

	REFERENCES

TOP ABSTRACT INTRODUCTION COMPUTATIONAL DETAILS RESULTS SUMMARY AND DISCUSSION ACKNOWLEDGEMENTS REFERENCES

 
Andersson, A., J. Almqvist, F. Hagn, and L. Mäler. 2004. Diffusion and dynamics of penetratin in different membrane mimicking media. Biochim. Biophys. Acta. 1661:18–24.[Medline]

Avbelj, F., P. Luo, and R. L. Baldwin. 2000. Energetics of the interaction between water and the helical peptide group and its role in determining helix propensities. Proc. Natl. Acad. Sci. USA. 97:10786–10791.[Abstract/Free Full Text]

Bechinger, B. 1997. Structure and functions of channel-forming peptides: magainins, cecropins, melittin, and alamethicin. J. Membr. Biol. 156:197–211.[CrossRef][Medline]

Berendsen, H. J. C. 1996. Bio-molecular dynamics comes of age. Science. 271:954–955.[CrossRef][Medline]

Berendsen, H. J. C. 2001. Reality simulation: Observe while it happens. Science. 294:2304–2305.[Free Full Text]

Berendsen, H. J. C., J. P. M. Postma, A. DiNola, and J. R. Haak. 1984. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 181:3684–3690.

Berendsen, H. J. C., J. P. M. Postma, W. F. Van Gunsteren, and J. Hermans. 1981. Interaction models for water in relation to protein hydration. In Intermolecular Forces. B. Pullman, editor. Reidel, Dordrecht. 331–342.

Berendsen, H. J. C., D. Van der Spoel, and R. Van Drunen. 1995. GROMACS: a message-passing parallel molecular dynamics implementation. Comput. Phys. Commun. 91:43–56.[CrossRef]

Berger, O., O. Edholm, and F. Jähnig. 1997. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys. J. 72:2002–2013.[Abstract/Free Full Text]

Berlose, J. P., O. Convert, D. Derossi, A. Brunissen, and G. Chassaing. 1996. Conformational and associative behaviours of the third helix of antennapedia homeodomain in membrane-mimetic environment. Eur. J. Biochem. 242:372–386.[Medline]

Billeter, M., Y. Q. Qian, G. Otting, M. Müller, W. J. Gehring, and K. Wütrich. 1993. Determination of the NMR solution structure of an antennapedia homeodomain-DNA complex. J. Mol. Biol. 234:1084–1097.[CrossRef][Medline]

Binder, H., and G. Lindblom. 2003a. Charge-dependent translocation of the trojan peptide penetratin across lipid membranes. Biophys. J. 85:982–995.[Abstract/Free Full Text]

Binder, H., and G. Lindblom. 2003b. Interaction of the trojan peptide penetratin with anionic lipid membranes: a calorimetric study. Phys. Chem. Chem. Phys. 5:5108–5117.[CrossRef]

Brattwall, C. E., P. Lincoln, and B. Nordén. 2003. Orientation and conformation of cell-penetrating peptide penetratin in phospholipid vesicle membranes determined by polarized-light spectroscopy. J. Am. Chem. Soc. 125:14214–14215.[CrossRef][Medline]

Braun, K., P. Peschke, R. Pipkorn, S. Lampel, M. Wachsmuth, W. Waldeck, E. Friedrich, and J. Debus. 2002. A biological transporter for the delivery of peptide nucleic acids (PNAs) to the nuclear compartment of living cells. J. Mol. Biol. 318:237–243.[CrossRef][Medline]

Brugidou, J., C. Legrand, J. Méry, and A. Rabié. 1995. The retro-inverso form of a homeobox-derived short peptide is rapidly internalised by cultured neurones: A new basis for an efficient intracellular delivery system. Biochem. Biophys. Res. Commun. 214:685–693.[CrossRef][Medline]

Cevc, G., and D. Marsh. 1987. Phospholipid Bilayers: Physical Principles and Models. John Wiley & Sons, New York.

Chen, F.-Y., M.-T. Lee, and H. W. Huang. 2003. Evidence for membrane thinning effect as the mechanism for peptide-induced pore formation. Biophys. J. 84:3751–3758.[Abstract/Free Full Text]

Christiaens, B., J. Grooten, M. Reusens, A. Joliot, M. Goethals, J. Vandekerckhove, A. Prochiantz, and M. Rosseneu. 2004. Membrane interaction and cellular internalization of penetratin peptides. Eur. J. Biochem. 271:1187–1197.[Medline]

Christiaens, B., S. Symoens, S. Vanderheyden, Y. Engelborghs, A. Joliot, A. Prochiantz, J. Vandekerckhove, M. Rosseneu, and B. Vanloo. 2002. Tryptophan fluorescence study of the interaction of penetratin peptides with model membranes. Eur. J. Biochem. 269:2918–2926.[Medline]

Christopher, J. A., R. Swanson, and T. O. Baldwin. 1996. Algorithms for finding the axis of a helix: Fast rotational and parametric least-squares methods. Comput. Chem. 20:339–345.[CrossRef][Medline]

Console, S., C. Marty, C. García-Echeverría, R. Schwendener, and K. Ballmer-Hofer. 2003. Antennapedia and HIV transactivator of transcription (TAT) "protein transduction domains" promote endocytosis of high molecular weight cargo upon binding to cell surface glycosaminoglycans. J. Biol. Chem