Next Article in Journal
Accelerated Evolution of Cytochrome c in Higher Primates, and Regulation of the Reaction between Cytochrome c and Cytochrome Oxidase by Phosphorylation
Next Article in Special Issue
Ocular Delivery of Therapeutic Agents by Cell-Penetrating Peptides
Previous Article in Journal
Multiomics Analysis of a DNAH5-Mutated PCD Organoid Model Revealed the Key Role of the TGF-β/BMP and Notch Pathways in Epithelial Differentiation and the Immune Response in DNAH5-Mutated Patients
Previous Article in Special Issue
Functional Peptides from SARS-CoV-2 Binding with Cell Membrane: From Molecular Dynamics Simulations to Cell Demonstration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations

1
School of Public Courses, Bengbu Medical College, Bengbu 233030, China
2
Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructure, Department of Physics, Nanjing University, Nanjing 210093, China
3
High Performance Computing Center, Nanjing University, Nanjing 210093, China
*
Authors to whom correspondence should be addressed.
Cells 2022, 11(24), 4016; https://doi.org/10.3390/cells11244016
Submission received: 16 November 2022 / Revised: 9 December 2022 / Accepted: 10 December 2022 / Published: 12 December 2022

Abstract

:
With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure–activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide–membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP–membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.

1. Introduction

Cell-penetrating peptides (CPPs), also named protein-transduction domains (PTDs), are usually amphipathic and/or cationic molecules consisting of 5~30 amino acid residues, and rich in basic amino acids such as arginine and lysine. They exhibit their biological activity by interacting with cellular membranes, and belong to the category of membrane active peptides (MAPs) as well as antimicrobial peptides (AMPs). Such peptides are typically unstructured in aqueous solutions, but fold into a well-defined secondary structure near the membrane [1]. Compared to AMPs, CPPs can penetrate cellular membranes without disrupting their integrity. Furthermore, CPPs can not only cross cell membranes, but also transport covalently or non-covalently linked bioactive substances, such as proteins [2,3], peptides, DNA, RNA, fluorescent labels, nanoparticles [4], viruses, liposomes [5,6], oligonucleotides [7], and drugs, into the cells. Regarding the capacity of delivering bioactive active substances into the eukaryotic cells, CPPs exhibit brilliant prospects in the field of therapeutics delivery [8]. Although CPPs can promote the delivery of cargo molecules, the nonspecific characteristics of CPPs themselves to cells may lead to a decrease in the specificity of the CPP-cargo complex and an increase in cytotoxicity. How to target CPPs and cargo molecules to specific cells effectively while maintaining their biological stability and biocompatibility in complex environments has always been an active area of research in the pharmaceutical realm. In recent years, a great number of studies have reported the potential of CPPs as carriers for the treatment of various diseases, such as cancers [9,10,11,12,13,14], diabetes [15], and Alzheimer’s disease [16]. It has been found that some CPPs have exhibited greater selectivity for cancer cells in preference to noncancer cells [17]. For example, poly-L-arginine could facilitate the cellular uptake of doxorubicin and increase the cytotoxic effect of doxorubicin in human prostate cancer DU145 cells [18], and polyarginine (R11) conjugated with fluorescein isothiocyanate FITC-R11 has also been proved to be a more specific molecular probe for bladder cancer and has shown potential applications as imaging agents in cancer treatment [19].
Over the past 30 years, the outstanding properties of CPPs have attracted considerable attention and significant progress has been achieved. A variety of experimental methods, such as nuclear magnetic resonance (NMR) experiments, circular dichroism (CD) spectroscopy, isothermal titration calorimetry (ITC), single-molecule fluorescence microscopy, fluorescence spectroscopy, dynamic light scattering, and lamellar neutron diffraction, have been used to explore the internalization process, penetration mechanism and potential applications of CPPs. Even so, differences in experimental strategies, experimental conditions, and analytical methods may result in contradictory results. In addition, current experimental techniques cannot provide atomic information with enough temporal and spatial resolution about the internalization process of CPPs. As a consequence, the uptake mechanism of CPPs is still unclear despite having a host of alternative models.
With the development of computational resources and methods, molecular dynamics (MD) simulation has become a powerful tool to explore the internalization mechanisms, which can provide kinetic clues to peptide–membrane interactions and cellular uptake processes at the atomistic level [20,21,22,23,24]. In addition, we can monitor various effects during the process of penetration by changing simulation conditions such as temperature, pH condition, initial conformation or orientation of the peptides, and the lipid composition of the membrane. For example, by decreasing the transmembrane potential, Trofimenko et al. coined a megapolarization situation, which could trigger the formation of water pores, and enable the direct translocation of CPP into cells [25] (Figure 1). As an effective complement to the experimental results, MD simulations have been used to study the CPPs with great success.
In this review, we begin with a brief introduction to the general properties of CPPs, including the CPP families, categories and cellular uptake mechanisms. Then, we describe the progress of MD simulations on CPPs from four aspects, including the computational methods and technical improvements in MD simulations, the research achievements in the uptake mechanism, the decoration and coupling of CPPs, and the influences and reactions of membrane models during the penetration process.

2. Categories and Cellular Uptake Mechanisms of CPPs

Since the discovery of the first CPP in 1988, thousands of CPPs have been discovered or designed artificially. CPPsite 2.0 [2], the database of cell-penetrating peptides, contains a total of 1855 entries of CPPs along with their secondary and tertiary structures now. Some of these CPPs have been well studied. Taking the HIV-1 Tat peptide (GRKKRRQRRRPPQ) as an example [26,27], the peptide contains eight positively charged residues and has no regular structure both in solution and on the membrane. However, it has been found that this kind of peptide can facilitate the uptake of large, biologically active molecules into the mammalian cells, and can mediate the delivery of different cargoes into tissues throughout living organisms. Nowadays, the Tat peptides have been successfully used for intracellular delivery of a broad variety of cargoes, including various nanoparticulate pharmaceutical carriers (liposomes, micelles, nanoparticles) [28,29]. Another typical example is penetratin (RQIKIWFQNRRMKWKK) [30,31,32], which is derived from the third helix of the Drosophila Antennapedia homeodomain. Penetratin has seven positively charged residues, and can penetrate directly into the giant monolithic vesicles [33]. It is worthy to note that both Tat and penetratin are arginine-rich peptides, suggesting that the arginine plays an important role in the uptake mechanism due to the bidentate hydrogen-bonding interaction between the guanidinium groups of arginine residues and the phosphate groups in the membrane, e.g., the polyarginine sequence (Arg8) has been demonstrated to be sufficient to drive molecules into cells [34]. Besides arginine, lysine also plays an important role in the CPP penetration process, such as Transportan 10 (TP10) (AGYLLGKINLKALAALAKKIL) [35,36], which is a lysine-rich chimeric CPP and adopts an α-helical secondary structure on the surface of membranes. Compared to Tat and penetratin, TP10 has more helical structure and shows higher efficiency in cellular uptake as the increase in CPPs’ helicity can promote the ability of cell-penetration. Since the previous findings, more and more CPPs have been designed to facilitate the penetration of cells across membranes into the cytoplasm.
As we know, CPPs have a common ability to penetrate cell membranes, although they are different in length, electric charge, hydrophobicity/hydrophilicity, flexibility, and so on. There are a variety of methods for classification of the CPPs [37]. The first one was based on the origin of the CPPs. In this manner, CPPs are classified into natural CPPs, designed CPPs and CPPs derived from peptide libraries. The natural CPPs are derived from natural proteins or peptides, such as the HIV-1 Tat peptide which is derived from the transcription transactivating protein of HIV-1 and the penetratin derived from the Antennapedia homeotic transcription factor. The designed CPPs are designed according to specific requirements or to mimic the structures of the known CPPs. A typical example for designed CPPs is transportan, which was designed by fusing amphipathic peptide mastoparan from bee venom to the fragment of human neuropeptide galanin. In recent years, a number of design methods were proposed to further enhance CPPs’ properties through integrating several different protein fragments into chimeric sequences, or incorporating of unnatural amino acids into peptides, and create numerous efficient CPPs, such as CADY (GLWRALWRLLRSLWRLLWRA) [38] and MAP (KLALKLALKALKAALKLA) [39]. Due to the lack of quantitative understanding of the sequence-structure-function relationship of CPPs and their complex mechanisms of action, it is difficult to achieve rational design and optimization of them. One possible approach is to obtain CPPs from the peptide libraries (large, diverse collections of related molecules) through high-throughput screening for desirable properties, such as targeting and delivery. The second classification method is based on the physicochemical properties (Figure 2), and CPPs can be divided into the cationic, amphipathic and hydrophobic ones. The cationic CPPs, such as HIV-1 Tat, penetratin and the Arg9 (RRRRRRRRR) [40] peptide, contain usually both arginine and lysine residues, and carry net positive charges which are important for the first contact with the membrane. The amphipathic CPPs that have both hydrophilic and hydrophobic regions, such as transportan, MAP and pVEC (LLIILRRRIRKQAHAHSK) [41], are the most numerous category, accounting for more than 40% of the known CPPs to date. They are strongly associated with neutral and anionic membranes, and their transmembrane transport patterns are closely related to their interactions with polar head and lipid chain regions of the membranes. The amphipathic CPPs can be further divided into two subclasses: primary amphipathic CPPs and secondary amphipathic ones [42]. The primary amphipathic CPPs normally contain more than 20 amino acids, have sequentially hydrophobic and hydrophilic domains along their amino acid sequence, and bind to both neutral and anionic lipid membranes with strong affinity, such as Pep-1 (KETWWETWWTEWSQPKKKRKV) [43] and TP10. They can penetrate into membranes with both high and low anionic lipid ratios and induce membrane leakage at low micromolar concentrations. The secondary amphipathic CPPs, such as MAP and penetratin, show a stronger affinity to anionic membranes than electrically neutral membranes. They show their amphipathic property through the α-helical or β-sheet folding, i.e., the distribution of their hydrophobic and hydrophilic domains occurs through the secondary structure formation. For example, penetratin has a random coil conformation in aqueous solution, while upon interacting with the anionic lipids it can form the α-helical or the β-sheet amphipathic secondary structures at lower or higher peptide-to-lipid ratios, respectively. As for the hydrophobic CPPs, such as K-FGF (AAVLLPVLLAAP) [44], they contain only apolar residues, or few charged amino acid residues in the sequence. The internalization mechanism of the hydrophobic CPPs are less well studied, but there is an opinion that their hydrophobic sequences are often critical for CPPs to penetrate cellular membranes. The third classification method is based on the CPPs’ cytotoxicity, which divides CPPs into the non-toxic and the antimicrobial ones. As with other chemical substances, the cytotoxicity of the CPPs may depend on the concentration, the amino-acid sequence, the composition of the membrane and so on. In addition, CPPs can also be classified as the nonspecific or the targeted ones depending on the capacity of cellular uptake. How to effectively target CPPs and cargo molecules to specific cells has always been a topic of interest in the area of drug development.
Although CPPs can penetrate biological membranes in a variety of ways, the uptake and internalization mechanisms have not been well understood [46]. In general, the cellular uptake mechanisms can be divided into two categories: energy-dependent endocytosis and energy-independent direct translocation (Figure 3). For the first category, there are three main pathways of endocytosis, namely macropinocytosis, caveolae-mediated endocytosis, and clathrin-mediated endocytosis [47]. Interestingly, some studies have found that many CPPs cross cellular membranes by endocytosis at low concentrations, while at high concentrations they penetrate membranes through the direct translocation mechanism [48]. For example, as the peptide concentration increases, the way pTat and R8 peptides enter HeLa cells switches from endocytosis to direct translocation [49]. Although the hydrophobic core of biological membranes is the major barrier to charged CPPs, when the peptides aggregate on the membrane surface, the interactions between the positive charged CPPs and negative charged membrane components result in the lipid flipping or in pore formation, and CPPs subsequently penetrate the hydrophobic lipid core of the membrane with limited toxicity, suggesting that direct penetration is an ideal approach for designing intracellular delivery vectors. The transition of uptake mechanisms from endocytosis to direct penetration can also be evoked by lowering the temperature. For example, it did not abrogate HIV-1 Tat uptake by knock-down of clathrin-mediated and caveolin-mediated endocytosis, or even by incubation of the cells at 4 ℃ [50], suggesting that the Tat peptides penetrate the membrane directly in this case, and the direct translocation mechanism of CPPs is related to the instability of membrane to some extent. There are a host of models to explain the direct translocation, such as the pore formation model [51,52], inverted micelle model [53,54], carpet model [55,56], and membrane thinning model [57]. It is worth mentioning that there is no unified conclusion about the uptake mechanisms of CPPs, and it is agreed that the cellular uptake pathway depends on many factors, including the physicochemical properties such as charge distribution and peptide length [58], the secondary structures of peptides [1], CPP concentration [59], cell type [60], lipid composition of the membrane [61], the ability of CPP to interact with the membrane, the property of the cargo molecule and its coupling type with CPP [62], experimental conditions [63], environmental factors like pH condition and ambient temperature [59].

3. Techniques for CPP Simulations

The earliest MD simulations were used in 1957 to study the motion of a simplified hard-sphere system [64]. In 1964, Rahman simulated the systems of 864 Ar atoms with the Lennard-Jones potential [65], obtained the radial distribution function, diffusion constant, velocity autocorrelation function, and mean azimuth shift of the system, and compared them with experimental results. Over the decades, molecular simulation methods have made great progress with the development of simulation algorithms, force fields and computational resources. MD simulations allow us to rigorously control the composition and condition of the system, monitor the evolution of the system, and analyze the interactions at the atomic or molecular level. Thus, through the detailed simulations, we can obtain lots of microscopic details for CPP internalization by monitoring various indicators [66] such as evolution of the peptide conformation and orientation, time dependence of the distance and the contact area between peptide and membrane, depth of insertion of the peptide into the membrane, radial distribution functions, density profile of the system, charge distribution, deformation and reorganization of the lipid bilayer induced by the peptide insertion. Moreover, we can also obtain accurate peptide–membrane interactions and free energy landscape by analyzing hydrophobic/hydrophilic properties, electrostatic interactions, hydrogen bonds, interaction energies, the potential of mean force (PMF), etc. MD simulations can complement our understanding with information that is hard to gain from experiments, and excite new ideas on experiments to further examine the findings from simulations. Nevertheless, we should also keep in mind that the results of MD simulations depend to some extent on initial conditions, degree of equilibrium, system size, and the conditions under which the simulations are performed.

3.1. All-Atom MD Simulations, Coarse-Grained Simulations and Implicit Membrane Models

In general, MD is a method to simulate molecular motion based on Newtonian mechanics. The forces acting on each atom can be obtained according to the molecular potential function, and by given the initial velocity and position, we can obtain the trajectory of particles as a function of time. MD simulation can provide dynamic and thermodynamic information of the system in atomistic detail, and has become an effective tool for understanding the structure and function of biosystems. Classical all-atom MD simulations have higher resolution but are very time-consuming and computational cost expensive. In order to decrease the computational costs, improve the sampling efficiency and speed up the computation, some researchers adopt the coarse-grained models to perform the simulations. For example, Hu et al. explored the dependence of translocation free energetics on cyclic and linear Arg9 peptide conformation using coarse-grained MD simulations with the MARTINI force field (Figure 4) [67]. The results were consistent with the all-atom force field simulations, showing that the current coarse-grained force field provides sufficiently robust and reliable free energy differences. For coarse-grained simulations [68], groups of atoms are mapped onto different interaction beads, this allows the simulations to scale up to larger systems and longer timescale, but at the expense of the model accuracy, some behaviors and properties, such as the conformational transition of peptides, the protein folding events, the directionality of hydrogen-bonding patterns that underlies protein conformational stability, will be difficult to be investigated [69]. Given the balance between computational cost and simulation accuracy, another option is to consider implicit membrane models. Implicit membrane models attempt to account for the influence of lipids and water through the solvation free energy term in the energy function, which have many advantages such as fast speed, rapid equilibrium, and extensive exploration of configurational space. One of the problems with the implicit membrane models is that these models are only suitable for plane membranes because the membrane cannot be deformed, which introduces flaws to peptide–membrane interaction studies. To overcome this shortcoming, Lazaridis et al. incorporated the data from both all-atom simulations and experiments, and estimated the elasticity of lipid membrane by improving the implicit membrane model IMM1 and adding the free energy of membrane deformation [70]. For a helical secondary structure, the minimum energy pathways and transition states were determined by mapping the effective energy of the peptide as a function of vertical position and tilt angle. They used this model to study the transport of peptides, such as S4 helix, TP10W and DL1, across the membrane bilayers, and obtained results consistent with the experiment. In addition, the free energy predicted by IMM1 for the transfer of amino acid sidechains from bulk water to the membrane center was quite close to the calculated results based on all-atom simulations. However, the shape of the free energy profile of the charged side chains needed to be modified to reflect the all-atom simulation results. Later, the same group extended the implicit membrane model 1 (IMM1) to spherical and tubular membranes [71]. The expanded IMM1 was tested in a series of systems, including the wild-types and mutants of the antimicrobial peptide magainin, and the results were qualitatively consistent with experimental results.

3.2. Enhanced Sampling Approaches

Despite MD simulations providing an effective method as a special complement to the experimental approach, there are still many limitations, such as the difficulty of adequate sampling of energy landscapes, over-simplification of the membrane model, and limited accuracy of the force fields. Not only that, but the time it takes for CPPs to spontaneously translocate through the cell membrane is on the order of minutes, which is almost impossible to achieve in MD simulations unless the internalization process is accelerated artificially. Therefore, in the studies of the peptide–membrane system, some enhanced sampling approaches are often employed, such as steered MD (SMD), umbrella sampling, metadynamics, and replica exchange methods. SMD usually accelerates the process by providing an external steering potential or force to move the system along a specified path [72], and in peptide–membrane system, SMD is commonly used to pull the peptides across the membrane to speed up the penetration process. For example, Yesylevskyy et al. used a harmonic potential with a force constant of 1000 kJ/mol/nm2 to pull a single penetratin molecule through a DPPC bilayer to mimic the translocation of penetratin [73] (Figure 5). The umbrella sampling method can calculate the PMF along the reaction coordinate by the weighted histogram analysis method (WHAM) based on the data from a series of windows (each of which performs an MD simulation) [74,75]. Metadynamics can reconstruct the free energy surface as a function of few collective variables (CV) by driving the system to keep away from previously visited areas and to sample the unfavorable regions of the free energy landscape by a history-dependent bias potential along the CVs [76,77]. Replica exchange molecular dynamics (REMD) is a method which combines MD simulations with a Monte Carlo algorithm [78]. In REMD simulations, multiple replicas of the same system are simulated in parallel, with occasional exchange of temperatures or interaction energies between adjacent copies. In this way, REMD can overcome high-energy barriers and sample the conformational space sufficiently [79,80,81]. Another method to calculate the free energy surface of the reaction coordinate system using non-equilibrium dynamics is adaptively biased molecular dynamics (ABMD) method [82], which is characterized by the need for fewer control parameters, the flooding timescale and the kernel width, as well as the favorable O(t) scaling of molecular dynamics time. Therefore, it is suitable for long-term biomolecular MD simulation. The ABMD method also allows extensions based on replica exchange, which greatly improves the speed and accuracy of the method. Gimenez-Dejoz et al. used ABMD to determine the free-energy landscapes for the CPPs containing unnatural amino acids crossing membranes [83]. These simulation techniques mentioned above can provide information on the energetics and dynamics of various stages of CPP penetration, such as peptide binding to the water-lipid interface, peptide insertion into the membrane, and pore formation on the membrane. Another noteworthy approach is the Markov state models (MSMs) [84]. They are commonly used to determine events such as membrane penetration or translocation from long unbiased simulations, especially when good order parameters are not known a priori, and are therefore suitable for molecular models of CPP-membrane interactions. The MSM is mathematically expressed as an n × n transition probability matrix, where the entry in row i and column j represents the probability of transitioning from state i to state j within the lag time. If the system is in one of n discrete substates, after a fixed time, the system will be found in any of these n discrete substates. Thus, the long-term statistical dynamics of molecules can be modeled as a Markov chain on a discrete state space, and by running many simulations of short trajectories from each state, the transition rates between states can be precisely determined, resulting in long-term dynamical information [85,86]. MSM has been widely used to study a variety of biophysical problems, such as protein dynamics [87], protein folding [88], protein–ligand association [89], molecular self-assembly and aggregation [90]. For example, Xu et al. demonstrated that MSM analyses could be used to delineate the variation of free energy during the self-assembly process of a typical amphiphilic DPPC lipid [91].
Since umbrella sampling restraints the peptide through successively pulling, the calculated free energy should be considered as free energy along a particular path. Although it is difficult to obtain a full free-energy profile, comparisons between free energies along different paths can still shed insight on the translocation mechanism. Huang et al. calculated the free energies of translocating Arg9 into a neutral DOPC lipid bilayer along two different paths using the umbrella sampling MD simulations [92], and the results showed that the free-energy barrier along the pore path is 80 kJ/mol lower than the one along a pore-free path.
For the enhanced sampling techniques, the key factor in effectively obtaining the free energy profile is to choose a good CV. Kabelka et al. developed a new CV that can describe peptide insertion, local membrane deformation, and the internal degrees of freedom of the peptide associated with its charged groups [93]. By comparing different CVs, they demonstrated that all of these components are necessary to accurately describe peptide translocation, and the developed CV can reliably and effectively calculate the free energy of peptide translocation.

4. Simulations on CPP Internalization Mechanisms

In order to better develop and optimize CPP application strategies, it is important to understand the cellular uptake mechanisms of CPPs. Although all CPPs share some common properties of cell penetration, different CPPs have different interaction partners with membranes, different membrane responses induced by peptide insertion, and different translocation mechanisms [94] (Figure 6). The interactions between CPPs and the membrane may come from various contributions such as electrostatic interactions, hydrophobic effects, and hydrogen bonding. In addition, depending on experimental conditions, most CPPs can enter cells through multiple pathways. Although currently there are a variety of experimental methods to explore the proposed CPP internalization mechanism and manifold theoretical models, it is still necessary to verify these findings and investigate the mechanistic details from the perspective of simulation. MD simulations provide insights into these peptide–membrane interactions, help us to analyze the dynamics of internalization processes in depth, provide information on charge, hydrogen bonding, orientation, and helical properties of peptides which can be accurate to specific residues and demonstrate the importance of certain specific residues or events during the penetration process.

4.1. Role of Arginine Residues

In general, CPPs have no sequence homologies, and different residues and sequence arrangements may play a role in the process of penetration. This structural diversity leads to different uptake patterns. In most CPPs uptake mechanisms, the usual initial step is the electrostatic interaction of positively charged CPPs with negatively charged membrane components such as phosphate groups, which affects the organization of lipid molecules and leads to membrane perturbation. It has been reported that the effective penetration of CPPs requires at least eight positive charges [58]. It is well known that positively charged arginine residues play an important role in the uptake process, due to the strong interaction between the guanidinium group of Arg and the phosphate groups of lipids, that is, the formation of bidentate hydrogen bonds and electrostatic interactions between the guanidine and the sulfate, the phosphate or the carboxylate moieties of membrane, polyarginine CPPs can cross membrane more efficiently. This strong interaction will destabilize the packing of the lipid and distort the membrane structure, and the charged amino acid residues route water molecules into the hydrophobic core of the lipid bilayer, thereby inducing a hydrophilic pore in the membrane, and facilitating the transfer of CPPs. As we know, antimicrobial peptides are able to translocate across the membrane and form stable pores at the millisecond level during the translocation [24,95]. The arginine-rich CPPs can also form a short-living hydrophilic channel in the plasma membrane of cells upon translocation [96,97,98]. In Ref. [52], Herce and Garcia successfully observed the pore formation caused by polyarginine peptides, and proposed a mechanism for the spontaneous translocation of the arginine-rich HIV-1 Tat peptides across a DOPC lipid membrane. MD simulations in the work showed that on the timescale of 100–200 ns, strong interactions between the Tat peptides and the phosphate groups of the lipid bilayer induced the insertion of charged side chains and the formation of a transient pore, followed by the translocation of the Tat peptides. But the results are controversial because the pore formation in the simulations was observed at a high ratio of peptide to lipid of 1:18 and at elevated temperature of 343 K and 363 K, and meanwhile without counter ions. Particularly suspicious is that the time scale on which peptide penetration occurs is too short compared to experimental observations (on the order of minutes). Yesylevskyy et al. studied the interactions of penetratin and the Tat peptides with DPPC and DOPC phospholipid bilayers by MD simulations at a lower ratio of peptide to lipid of 1:64 [73]. In contrast to the previous simulation work, no transmembrane pores were observed, suggesting that the formation of pores can only be induced above a certain concentration threshold. They further demonstrated that the free energy barrier obtained by the umbrella sampling simulations of inserting a single penetratin peptide into a DPPC bilayer was estimated to be ~75 kJ/mol, suggesting that the spontaneous penetration of a single peptide takes at least a few seconds to a few minutes on a timescale, which is consistent with experimental observations. Huang et al. also demonstrated that water pores can form in the membrane from an energetical point of view [92]. They calculated the free energies of translocating a cycle cationic peptide Arg9 into a neutral DOPC lipid bilayer along two different paths by umbrella sampling MD simulations. The results showed that the free-energy barrier along a pore path was 80 kJ/mol lower than that along a pore-free path. Similar conclusion for the same system was also confirmed in a coarse-grained simulation [67]. Another simulation of penetratin with DPPC lipid bilayer also confirmed that when the penetratin penetrated into the bilayer, a rather significant thinning of the bilayer at the contact site was observed, which may be the first stage of pore formation [99]. Yao et al. studied the transmembrane process of arginine-rich CPP Arg9 and tried to understand the quantitative correlation between the hydrophobic barrier of the plasma membrane and the free energy barrier of the CPP transmembrane using MD simulations and the umbrella sampling methods [100]. The energy analysis showed that multiple salt bridges of guanidinium-phosphate accounted for 65% of the overall interaction energy, whereas the interaction energy from the backbone of the peptide could be negligible, so the increased negative charges of the lipid bilayer or more salt bridges would reduce the transmembrane free energy barrier and further facilitate CPP adsorption and transmembrane processes.

4.2. Role of the Secondary Structures of CPPs

It is generally believed that the role of positive charges is critical in the uptake mechanism. The importance of positive charges, especially those of arginine residues, in the uptake mechanism has been demonstrated in [40,101], but the high transport efficiency of TP10, which is lysine-rich and arginine-deficient, suggests that this is not always the case. The work by Magzoub et al. showed that penetratin preferentially interacts with negatively charged membrane components, whereas transportan interacts with membranes independently of charge [102]. In addition, the secondary structures of CPPs, especially the structures adopted when interacting with the membrane, are also important for penetration. The α-helical peptides have proved to be better penetrators than the peptides with random coil structures, especially for pore formation models. Taking the arginine-deficient but lysine-rich Pep-1 for example, CD analysis revealed a transition from unstructured to helical conformation upon increasing the concentration, which is conducive to the insertion of Pep-1 into the membrane by forming a transient transmembrane pore structure [103]. Dunkin et al. used MD simulations to study the interaction between the α-helical amphipathic TP10 and a zwitterionic POPC bilayer [104], and assessed the plausibility of the “sinking raft” model, in which the peptide oriented parallel to the membrane surface, resulting in a mass imbalance across the bilayer, and this perturbation caused the peptide to sink deeper into the bilayer. The results showed that TP10 preferentially binds to the membrane in a parallel orientation, which was energetically favorable because the Lys residues are distributed along the sequence on the same side of the peptide, and this orientation kept the Lys residues in a polar environment during insertion. The simulation results also demonstrated the importance of the amphiphilic structure of TP10, which differs from some peptides (such as Tat) with unstructured conformations at the interface, and the analysis further revealed that the formation of salt bridges between the Lys residues of TP10 and the phosphate groups of the lipids is a key factor in determining the orientation of the peptide as well as stabilizing the peptide–membrane interactions, and the salt bridges can be compared with the Arg-phosphate ones of the HIV-1 Tat peptides determined by NMR [105]. Song et al. further evaluated the effects of structure and charge on the translocation ability of TP10 and its analogs [106]. The results showed that disrupting the helical structure or performing Arg substitution could remarkably reduce the cellular uptake ability of TP10, while increasing the positive charge and the amphipathicity was an effective strategy for enhancing the translocation ability of TP10.

4.3. Role of Hydrophobic Aromatic Residues

As for penetratin, which consists of seven basic residues (3 Arg and 4 Lys) and three aromatic residues (1 Phe and 2 Trp), the conformational changes are more complex. Polyansky et al. explored the intramolecular and intermolecular interaction characteristics of penetratin at the water–lipid interface and the conformational plasticity of the peptide by simulating the binding of penetratin to the zwitterionic DOPC and anionic DOPS lipid bilayers [107]. The simulation results showed that penetratin exhibited complex conformational behavior, that is, it maintained helical conformation in an aqueous solution and zwitterionic bilayer, while exhibiting high conformational flexibility in negatively charged membranes. It was also revealed that the conformation of penetratin depended on its initial orientation relative to the membrane, especially when it was in anionic DOPS bilayers. The major contribution of the peptide–membrane interaction was the salt bridges and hydrogen bonds formed by the basic Lys and Arg residues with the lipid polar heads. To evaluate the role of specific intramolecular interactions in penetratin, free energy perturbation (FEP) calculations were performed on aromatic–cationic residual–residue contacts of several model pairs. The results showed that the basic residues Arg and Lys could form energetically favorable pairs with aromatic residues Phe and Trp in the aqueous and apolar environments, which facilitated the insertion of the peptide and enhanced the stability of the membrane-bound state [108]. The formation of the above pairs was important because it was found that the energy gain extrapolated from the FEP calculations might be comparable to that observed upon the creation of backbone hydrogen bonds, which agreed well with the NMR experiment results of penetratin in the membrane mimic media-TFE-water mixture [109] and SDS micelles [110]. Pourmousa et al. also demonstrated the importance of the basic and aromatic residues of penetratin on the binding mode of penetratin to a zwitterionic DPPC lipid bilayer by MD simulations [99]. They found that penetratin exhibited a high degree of conformational flexibility during the simulations, it formed α-helical or β-like conformations in independent different simulation runs, which indicated that penetratin might fine-tune its structure in different membrane surroundings. The authors interpreted that all the different conformations could exist in the penetration process, and the peptide might switch between them when they penetrated into the bilayer. The simulation results also showed that for penetratin, the hydrogen bonding and charge–pair interactions between the bilayer and the Arg and Lys residues could help to localize the charged residues in the hydrophobic region of the bilayer, which was critical for the penetratin binding mode. As for the aromatic residues, Trp could form hydrogen bonds with the phosphate groups of lipids, while Phe formed no hydrogen bonds with lipids.
Other simulations have also observed the contact of basic residues with hydrophobic aromatic residues such as tryptophan during the internalization process of amphiphilic CPPs. Jobin et al. used experimental approaches as well as MD simulations to investigate the membrane insertion of the amphipathic peptide RW16 (RRWRRWWRRWWRRWRR) into the zwitterionic DOPC bilayer and the orientation properties of the peptide during the process [111]. The trend of RW16 insertion was determined by the hydrophobic interaction with the fatty acid chains via Trp residues and the electrostatic interaction with the lipid phosphates via Arg residues. Such interactions are important in the cellular uptake mechanisms because they create a small curvature of the membrane, which can induce the invagination phenomena [112,113] and contribute to the tubulation and internal vesicle formation as induced by RW16 on giant unilamellar vesicles (GUVs) [114]. Moreover, data obtained by NMR spectroscopy and MD simulations both revealed that side chain contacts between Arg and Trp residues contributed greatly to the conformational stability and the biomolecule functions [115,116]. Walrant et al. studied the binding, insertion and orientation of RW9 (RRWWRRWRR) and RL9 (RRLLRRLRR) into membranes [117], where RW9 was internalized into cells effectively, whereas RL9 was not internalized but bound to the cell membrane. During the simulation, the orientation of the peptides changed, and many hydrogen bonds and salt bridges formed between Arg side chains and the lipid phosphate groups, providing anchors for the peptides. It is worth noting that RL9 actually inserted deeper than RW9 and the result was consistent with what was observed experimentally in SDS and DPC micelles. The reason for the difference is that the Leu residues are more hydrophobic and may promote better penetration. However, deep insertion may impair more superficial interactions which could trigger a direct translocation [97], that is, the deeper insertion of RL9 may result in peptide trapping in the membrane. The difference in penetration between RW9 and RL9 might be due to their different effects on lipid packing. RW9 can decrease the lipid chain ordering, allowing lipids to adopt different supramolecular organizations to promote the passage of the peptides, where the Trp residues play a decisive role.

4.4. Other Factors Affecting the Uptake Mechanism

There are also several other factors that affect the uptake mechanism, such as initial configuration, changes in CPP conformation, concentration, and length of CPP sequences. On the effect of peptide length, He et al. performed coarse-grained MD simulations to investigate the interaction mechanism between polyarginine peptides with 4, 8, 12, 16 Arg, respectively, and asymmetric membranes including DPPC, DPPE and DPPS [118]. The results indicated that the peptides could penetrate through the lipid bilayer by inducing hydrophilic hole formation. Furthermore, long peptide chain length and high membrane asymmetry level have a positive effect on peptide penetration. In many cases, CPP concentration has shown a significant effect on the pathways and efficiency of penetration. In addition, some studies have shown that polymerized CPPs have better penetration ability than non-polymerized CPPs due to their stronger interactions with the membrane components. Some peptides exhibit their activities depending on the environmental conditions. A typical example is the pH low-insertion peptide (pHLIP), which forms a transmembrane helix and insert into the lipid bilayer under acidic pH conditions, and has shown promising applications in tumor-targeted delivery of cargoes [119,120]. Wei et al. investigated the interaction of a 35-amino acid pHLIP peptide (AEQNPIYWARYADWLFTTPLLLLDLALLVDADEGT) [121] with the POPC bilayers under acidic and basic pH conditions, using all-atom MD simulations with different initial configurations [122]. The results showed that when pHLIP was placed parallel to the bilayer surface, it could insert into the POPC bilayer spontaneously at acidic pH, while at basic pH it bound to the membrane surface without insertion. When pHLIP partially pre-inserted in the POPC bilayer, it could insert deep into the bilayer at acidic pH, while at basic pH it moved toward the bilayer surface.

4.5. Comparison and Evaluation of the Characteristics of Different CPPs

In order to compare and evaluate the characteristics of different CPPs during the penetration, Ramaker et al. tested the penetration ability of 474 CPPs from the CPP database for HeLa cells, and the results showed that there were significant differences in penetration efficiency, and none of the CPPs performed best in all conditions [123]. However, CPPs that form stable α-helical structures and have positive charges are good candidates for efficient cargo transport. Kumara et al. investigated the free energy profiles of a single peptide inserting into a neutral DOPC lipid bilayer membrane using coarse-grained MD simulations [124]. The selected 17 CPPs span cationic, hydrophobic and amphipathic categories. The results indicated that the adsorption of hydrophobic peptides was the strongest, and the free energy minima of these peptides were close to the bilayer interface. Among the selected CPPs, K-FGF (AAVALLPAVLLALLAP) acted as the best hydrophobic CPP, consistent with the fluorescence and CD spectroscopy studies [125]. It should be noted that, given the need to release into the cytosol, the strong adsorption of CPPs might be a disadvantage. However, cationic CPPs were found to exhibit only moderate adsorption capacity, that contradicted experimental results in which cationic peptides had been shown to be effective in translocating the membrane, possibly due to the fact that the charge on the bilayer and the concentration of the peptide were not taken into account. The best performers confirmed by experimental results among the cationic peptides were penetratin and R10 (RRRRRRRRRR) [123]. In addition, simulation results showed that Transportan [126] and C6 were the best amphipathic CPPs for passive translocation. In MD simulation studies, the penetration behavior of CPPs is mainly achieved through peptide induced membrane perturbations, resulting in local membrane thinning or the formation of hydrophilic pores. For example, MD simulations of TP10 and its analogs in the DPPC bilayer suggested that the degree of local membrane thinning could reflect the extent of bilayer disturbance, and higher membrane disturbances led to higher translocation activity and higher cytotoxicity [106].

4.6. Larger Systems

If we adopt larger systems with more lipids, it is hopeful to get richer responses, such as changes in membrane curvature. Yesylevskyy et al. used MD simulations to study the interactions of penetratin and the Tat peptide with the zwitterionic DPPC and DOPC phospholipid bilayers. The extended system consisted of 8 penetratin or the Tat peptide molecules and large bilayer patches including 512 lipids [73]. The simulation results showed that the transmembrane pores did not form spontaneously, but instead the bilayer patches underwent large-scale deformation. The peptides could enter the cell by micropinocytosis, then multiple peptides induced deformations of the lipid bilayer, which subsequently led to the formation of small vesicles that encapsulate the peptides (Figure 7). This mechanism might be mainly attributed to the high charge density on the membrane surface, which was induced by the aggregation of the peptides. Even if this potential micropinocytosis mechanism could be observed, computer simulations of endocytosis are still lacking to date due to the limitations of the temporal and spatial scales of simulations, as well as the fact that endocytosis often requires the participation of membrane proteins. Compared to the pore formation model, the inverted micelle model and carpet model require a larger system, thus at this stage coarse-grained methods might be a better choice for these two models. Kawamoto et.al used coarse-grained simulations to investigate the dynamics of Tat and the bilayer membrane with 512 lipid molecules [127]. They found that the membrane curvatures or invaginations could lead to the formation of inverted micelle when the peptide and the lipid heads were strongly attracted, and the simulation results were in good agreement with the NMR and TEM studies [128]. They therefore proposed a mechanism in which the transition from bilayer to inverted micelle stems from the competition between the attractive potential energy of the peptide and lipid heads and the bending energy of the bilayer.

4.7. Design Strategies of CPPs

Although the mechanisms of cellular uptake of CPPs depend on a large number of factors, and the sequences of different CPPs vary greatly (Table 1), there are some similar structural properties that can be referred to in the design of CPPs [129]. Regardless of the type of cargo and the composition of the membrane, the first factor that needs to be considered is the net charge of the peptide. Many CPPs are rich in positively charged arginine residues, and their guanidine groups can form electrostatic interactions, bidentate hydrogen-bonding interactions and salt bridges with the phosphate groups of lipids. Compared with lysine, the guanidine groups on arginine help CPPs cross the membrane more effectively [130]. The second factor is the secondary structure of the peptide [1], especially the structure adopted after the peptide interacts with the cell membrane. For amphiphilic CPPs, the translocation is the result of amphiphilicity rather than positive charge, although most of these CPPs are positively charged. For example, if the MAPs are amphiphically conserved, the cellular penetration ability of neutral MAP17 (QLALQLALQALQAALQLA) and anionic MAP12 (LKTLTETLKELTKTLTEL) can be obtained by replacing lysine with other polar residues [45,131,132], and the results suggest that α-helical structure is more favorable for penetration. The third factor is the peptide concentration, which can trigger different uptake pathways. It is generally believed that amphiphilic peptides penetrate the membrane directly at low concentrations, while non-amphiphilic CPPs are taken up by endocytosis [94]. Besides the above factors, there are other factors that need to be considered when designing CPPs, such as hydrophobicity and the length of the peptide. In addition to the cell penetration properties, we also need to consider the penetration efficiency, the aggregation, and, especially, the cytotoxicity. Therefore, the design of CPPs is a complex task. In contrast, chemical modification of the structure of a known CPP or fusion of effective sequences may be a more economical way. MD simulations of CPP decoration will be discussed in the next section.

5. Simulations of CPP Decoration and Coupling

Over the past 30 years, the family of CPPs has expanded rapidly, and the search for new sequences continues. In recent years, numerous studies have reported the potential of CPPs as carriers in drug delivery, and the research field has shifted to the design of peptides and CPP-cargo complexes [129] with higher transduction efficiency, biostability, specificity, and lower cytotoxicity. A comparison between different experiments is relatively difficult due to many factors, such as cell type, cargo type and size, coupling modes of CPP and cargo, and experimental conditions. MD simulations offer a promising approach to facilitate the design of CPPs by decorating or substituting amino acid sequences, fusing valid sequences, or conjugating with small bioactive molecules. For example, improved analogues have been obtained by shuffling the amino acid sequence of penetratin [134]. Although these modified or decorated CPPs have not been accepted for clinical applications, their potential application value has received considerable attention and discussion [135,136].

5.1. Chemical Modifications of CPPs

Functionalization or chemical modification of CPPs, such as to include unnatural and modified residues [137] or to substitute the amino acid sequence of CPPs, has been proposed to improve the delivery of cargoes that successfully target specific cells or tissues. MD simulations also reported some studies on the modification or substitution of certain amino acids in CPPs. SAP(E) (VELPPPVELPPPVELPPP) [133] is known as the first synthetic anionic CPP by replacing the positive arginines with glutamic acids from a well-known CPP SAP (VRLPPPVRLPPPVRLPPP). SAP(E) showed similar uptake efficiency and action mechanism to SAP, suggesting that a positive net charge is not a necessary factor for the activity of CPPs. Based on the above reason, Antunes et al. presented a synthetic anionic peptide LE10 (LELELELELELELELELELE), which strongly interacts with model membranes and has demonstrated the properties of CPPs and AMPs both in MD simulations with the POPC and cholesterol bilayers and in in vitro studies with liposomes and mammalian cells [138]. The MD simulations showed that LE10 rapidly interacted with the bilayer surface and led to membrane bending, and at the end of simulations, the membrane involved the peptide, such as one half of a micelle. Simulations of larger systems with a lateral size of 26 nm showed that the curved region of the membrane eventually coalesced into a vesicular bud or a micelle, similar to the early stage of the endocytosis mechanism. Tsai et al. designed the SAP10 (RRWKFFPWRR) peptide, an almost mirror-symmetric peptide derived from the cationic antimicrobial peptide indolicidin (IL) (ILPWKWPWWPWRR), to improve gene delivery safety and efficiency [139]. All-atom MD simulations were performed to understand the association between SAP10 and POPC lipid bilayers. The results showed that the inserted IL and SAP10 peptides exhibited different configurations although they both disturbed the lipid bilayers. Because the cytotoxicity of IL had been reported to associate with Trp residues, the reduction of Trp in the SAP10 peptide resulted in reduced membrane perturbation and improved biocompatibility. In vitro experiments have demonstrated that the designed SAP10 is a safe and effective peptide that promotes PEI-mediated gene delivery. Another example comes from Gimenez-Dejoz’s team [83]. It was reported that the incorporation of nonproteinogenic amino acids, such as α-aminoisobutyric acid (Aib), could increase the helicity, biostability and penetration efficiency of CPPs [140]. They used MD simulations with enhanced sampling technique to study two Aib-containing CPPs, poly ( LysAibAla ) 3 (KaibA) and poly ( LysAibGly ) 3 (KaibG), and they both exhibited improved biostability and internalization abilities [141]. Compared to some widely known CPPs such as BP100 and R9, kAibA and kAibG showed the lowest internalization energies due to their compact amphipathic structures, and the Aib residues facilitated internalization by forming hydrophobic contacts with the lipid tails. There are a variety of modification methods, but it is important to note that each modification must be carefully designed to avoid possible problems, such as low solubility, aggregation or toxicity [129,142,143,144].
Another promising approach to improve the efficiency of CPP internalization is to fuse effective sequences, such as the design of the well-known TP10. Selective targeting of CPPS can be achieved through chimerism [66,145,146]. Snini et al. investigated the conformational preference of the fusion peptides Vim-TBS (58-81)-p10 (GGAYVTRSSAVRLRSSVPGVRLLQ-RQTSMTDFYHSKRRLIFS) and Tat (48-60)-p10 (GRKKRRQRRRPPQ-RQTSMTDFYHSKRRLIFS) [147]. These fusion peptides are derived from the Vimetin Tublin binding site (TBS) peptide, the Tat peptide, and a pro-apoptogenic peptide of the p21/WAFI protein, respectively, and the peptide p10 has frequently been used as a model cargo to evaluate cell-penetrating properties. It was found that the above fusion peptides will form helical structures during MD simulations. Compared to the Tat-p10 peptide, Vim-TBS-p10 was predicted to show better cell penetration properties because of its relatively better amphipathic nature and the hydrophilic face formed by the positively charged residues, that facilitated a better interaction with the membrane. The above finding has been well supported by the experimental observations showing that this fusion peptide accumulates and distributes in glioblastoma cells [148]. In the aspect of supramolecular assembly of CPP nanostructures, Ashwanikumar et al. created cell-penetrating self-assembling peptide nanomaterials (CSPNs) with the capability to translocate inside cells through sequential ligation of peptide building blocks [149]. They devised a Tat peptide-based triblock array, in which the Tat peptides were conjugated to an amphiphilic linker ( RADA ) 2 , followed by the sequential addition of Phe residues. This design enhances the hydrophobicity of the CSPN, and makes it a versatile host that can accommodate guest molecules. The designed CSPNs containing two, three or four Phe residues would self-assemble into “nanodrill-like structures”, which had high capacity to encapsulate hydrophobic guest molecules and enhanced the internalization. MD simulations of the Fmoc FF ( RADA ) 2 Tat (2F-RT) peptide aggregate were also carried out in the work to study the evolution of secondary structure and intermolecular interactions. During the simulation, the Tat segment lost its α-helical nature and formed random coils, while the amphiphilic ( RADA ) 2 retained its native helical structure for the most part. In addition, the peptides in the core of the aggregate fluctuated significantly less than their surface counterparts, which allowed them to form intermolecular salt bridges. From this work, we can see that MD simulations provide microscopic insights into the structural properties and mechanisms of action of fusion peptides, which is of great significance for guiding and facilitating the design of peptide drugs with cell penetration functions.

5.2. CPP Coupling with Cargoes

Designing a perfect analytical method to characterize the transport dynamics of CPPs and cargoes is a difficult task because the internalization process of CPPs is very complex and strongly depends on the type, size and binding method of the cargo molecules. In different experiments, CPP might function as a monomer or form complex with its payload molecule. For the latter, CPPs can form complexes with cargo molecules through covalent or non-covalent bonds. In the case of covalent conjugation, CPPs are coupled to neutral cargoes via disulfide bonds, amide bonds, etc., while in the case of non-covalent conjugation, CPPs are coupled to negatively charged cargoes by electrostatic and/or hydrophobic interactions. The conjugation of peptides offers a possible pathway for the design of new drug delivery systems and thus holds great potential in future treatments. For the CPP-dominated cargo delivery system, the internalization efficiency can be changed by modifying the properties of the cargoes and their binding methods. Alfonso T. Garcia-Sosa et al. studied the stability and binding affinity of the complexes of CPPs NF51, PF3, PF6, and TP10 to siRNA [150]. In the work, they found that the bound structures of all complexes were stable throughout the MD simulations, and the stability and binding affinity were related to the sidechains and modifications of the CPPs. The guanidino groups on arginine residues might bind to nucleic acids through electrostatic, hydrogen bonding, cation-π, and π-π interactions. The stearyl groups presented in PF3, PF6, and NF51 and the quinoline groups presented in PF6 increased the binding affinity and stability of their complexes with siRNA compared to TP10, which lacks these groups. CPPs with conjugated fatty acid moieties have shown excellent potential as novel delivery excipients as the increased hydrophobicity of the CPP carriers promotes stronger association with the lipophilic cell membrane [151]. Furthermore, this modification of CPPs might increase the self-assembling tendency of the peptides and then diminish the interactions with the membrane. To explore the impact of site and degree of lipidation on the membrane interactions of a cationic CPP, Hedegaard et al. designed four peptides lipidated on different lysine residues using capric acid (C10) based on the penetratin analogue PenetraMax (KWFKIQMQIRRWKNKR) [152]. Experimental results and coarse-grained MD simulations demonstrated that all CPP conjugates interacted with the membrane by inserting lipid chains into the core of the bilayer, and that the site and degree of lipidation affected the behavior of the carriers in their interactions with the membrane and thus altered the penetration efficiency of CPPs. In addition, membrane thinning effects and CPP induced membrane bending were also observed. The work revealed that the CPP performance could indeed be tuned by lipidation to facilitate membrane insertion. In most experimental studies, CPPs are tethered to fluorescent dyes in order to allow tracking of peptide transduction events under the microscope. Studies have shown that different fluorophore labeled peptides may also affect uptake mechanisms, intracellular distribution and cytotoxicity. To explore the possible role of dyes in assisting translocation, Akhunzada et al. used MD simulations to analyze the Tat peptides in monomeric and dimeric forms across a purposely-created membrane pore, whether or not covalently-linked to the TARMA dye [153]. The results showed that in comparison to the monomers, Tat dimers translocated much more slowly, thus their capability to stabilize membrane pores was enhanced due to the strong coulomb interactions between arginine residues and surrounding lipid head groups. Combined with other studies, a picture for the role that the TAMRA dye played in the process of Tat11 membrane transduction was drawn: it promotes peptide dimerization and, by doing this, it remarkably slows down the translocation kinetics, thus enhancing the lifetime of the membrane pores.
The coupling of CPPs with nanoparticles has also been studied by many MD simulations. Using coarse-grained MD simulations, Hu et al. showed that the transmembrane efficiency of polyarginine (R8) peptides was improved by conjugating with small hydrophilic nanoparticles (NPs) with appropriate linkers [154]. At high concentrations, the R8-NP conjugates could directly translocate across the asymmetry membrane through water pores. The translocation efficiency was closely related to the length of the linkers. When the linker length was about half of the membrane thickness, the translocation efficiency reached a maximum. However, the R8-NP conjugates with overlong linkers not only reducing the transmembrane efficiency due to the blockage of NPs in the water pores, but also causing cytotoxicity because the water pores were not closed for too long. Magnetic iron oxide nanoparticles have also been actively investigated for their possible applications in medicine, as they represent the next generation of targeted drug delivery technologies, and silica is the most common compound for the preparation of coated iron oxide nanoparticles. Grasso et al. investigated the adsorption of penetratin on silica and magnetite (MAG) surfaces, and estimated the penetratin-surface adsorption free energy quantitatively by umbrella sampling MD simulations [155]. The results indicated that Arg and Lys played a key role in surface binding, despite their physicochemical properties and surface charge densities differ, and this result was consistent with similar computational studies focusing on a spontaneous membrane penetrating peptide (PLIYLRLLRGQFC) adsorption to the silica surface [156]. The simulation study from Grasso’s team suggested a competitive mechanism between magnetic nanoparticles (MNPs) and cell membranes, both of which interact with the same sites of penetratin and may partially inhibit penetratin from carrying out its membrane penetrating function. In order to overcome the difficulty of low cellular internalization efficiency of gold nanoparticles (AuNPs), decorating AuNPs with functional moieties such as CPPs has been considered as a promising strategy. Quan et al. explored the interactions of the Tat peptides and their conjugated AuNPs with lipid membranes by coarse-grained MD simulations [157]. The results showed that the translocation of AuNPs across the membrane was significantly affected by the Tat concentration on the particle surface (Figure 8). When a small amount of peptides were decorated on the AuNP surface, the transmembrane efficiency of AuNPs might even be reduced, because the peptides adhering to the outer membrane surface might restrict the movement of AuNPs into the membrane. However, when the number of grafted peptides increased to a threshold, a transient hydrophilic pore formed in the lipid membrane due to the transmembrane electrostatic potential difference, so that the Tat-AuNP complex could translocate across the membrane in a pore-mediated way.

5.3. The Design of Self-Complementary Peptides

It is also noted that ionic-complementary peptides can be modified to exhibit similar properties to CPPs by adjusting the amino acid sequence and chain length. As a new type of nanomaterial with wide application value, the side chains of ionic self-complementary peptides are complementary to each other, that is, positively charged and negatively charged amino acid residues can form ionic pairs. This special structure determines that they can spontaneously assemble to form stable nanostructures [158], and the self-assembly is mainly driven by intra- and intermolecular hydrogen bonding, hydrophobic, and electrostatic interactions. In previous work, we performed MD simulations to investigate the adsorption and self-assembly behaviors of the amphiphilic ionic complementary peptide EAK16-II (AEAEAKAKAEAEAKAK) on the hydrophobic highly ordered pyrolytic graphite (HOPG) surface [159,160]. Simulations showed that the peptides tended to form hydrophobic contacts with the HOPG surface rather than to form a hydrophobic core, and the solution pH could affect the adsorption rate of EAK16-II. Under neutral pH condition, the interchain electrostatic attraction was favorable for adsorption, while under acidic and basic conditions, the electrostatic repulsion force slowed down the adsorption because of the protonation and deprotonation of glutamic acid and lysine residues. Hydrogen bonds did not play an obvious role during the adsorption, but several long-lasting interchain hydrogen bonds were observed after the peptides depositing on the surface, indicating that hydrogen bonds played an important role in the subsequent peptide self-assembly process.
Molecular self-assembly is a widespread phenomenon in nature. Self-assembling peptides are a kind of biomaterials that are relatively easy to design and synthesize, and have many potential applications such as templates for nanowire fabrication and delivery of peptide medicines. Studies have shown that abnormal folding and aggregation of proteins can lead to related diseases, such as Alzheimer’s, Parkinson’s, Down’s syndrome, and the nanostructure formed by peptides depends on the influence of the amino acid sequence and solution pH [161]. By regulating and modifying the amino acid sequences of peptides, the nanostructures generated by self-assembly of peptides can be controlled, which can be used in the design and application of new nanomaterials. Chen et al. proposed a de novo amino acid pairing (AAP) design principle, that is, using the hydrogen bond pairing, ionic and hydrophobic pairing strategies to construct new self-assembled peptides (SAPs), which could self-assemble into β-sheet rich nanofibers [162]. This kind of peptide showed the ability to stabilize and deliver hydrophobic anticancer agent ellipticine in an aqueous solution. Moreover, the peptide–drug complexes exhibited good antitumor activity against human lung carcinoma cells A549 and breast cancer cells MCF-7. For example, they put forward a systematic design of a histidine-rich lipidated peptide sequence named NP1 (Figure 9), which contains a hydrophobic stearic acid, histidines and positively charged arginines, where histidines promoted the proton sponge effect and arginines had a high affinity to bind to the cell surface for intracellular delivery [163]. The NP1 peptide successfully exfoliated graphite flakes, exhibited long-term dispersion stability in aqueous solution, and functionalized the graphene nanosheets, which not only exhibited good biocompatibility with cells, but also enhanced cell uptake of cancer drugs. The details of peptide–peptide and peptide–graphite interactions were investigated by experiment and all-atom MD simulations. MD results showed that the strength of the repulsive interactions might cause a single-layer peptide pattern with ∼0.2 nm height, and a multi-layered (or compact) structure with ∼3.1 nm height under the acidic and neutral conditions, respectively, as observed in AFM images. In another study, they tried to modulate the morphologies and internal cohesion of peptide amphiphile (PA) self-assemblies and their resultant functions by tuning the molecular interactions [164]. The coarse-grained MD simulations were applied in the work and revealed that Palm-CR3 (C16-CSNSNSRRR) self-assembled into cylindrical supramolecular fibers in water with a diameter around 8.0 nm, whereas the fiber diameter decreased to around 7.0 nm in a high ionic strength environment. Around the positively charged surface of the fiber, abundant chloride ions were observed while sodium ions were rare, resulting in the screened electrostatic repulsion. When electrostatic repulsion was screened, the peptide amphiphiles were supposed to interact with one another in a more compact manner than their “repulsion on” counterparts. The widths of Palm-CR3 fibers measured from the TEM images confirmed the decrease of fiber diameter when electrostatic repulsion is switched off. The simulations also found that the n-dodecyl-β-D-maltoside (DDM) surfactant molecules interacted with Palm-CR3 molecules and took part in the formation of variform clusters.

6. Simulations on Membrane Modification and Simulations of Multi-Component Membranes

6.1. The Effect of the Membrane Compositions on CPP Uptake Process

Besides the physicochemical properties and concentrations of the CPPs themselves, other factors that affect the uptake pathways of CPPs are the composition and properties of the membrane. The cell membranes constitute the first barrier for molecules and ions to enter the cells. Lipid membranes have complex and heterogeneous structures composed of diverse components such as different lipids, head groups, membrane proteins, and play an important role in the normal function of cells. Normally, the lipid bilayer can be divided into the choline group (head), the phosphate group (PO4), the glycerol ester group (GE), and the hydrocarbon chain (HC). Such a structure leads to the complexity of interactions between the membrane and CPPs and also affects the way CPPs are inserted. For example, negatively charged components of the membrane, such as phosphate groups and glycosaminoglycans, contribute to the accumulation of positively charged CPPs on the membrane surface, leading to the membrane destabilization, allowing the peptides to traverse the membrane. Mukherjee et al. studied eight different model membranes containing zwitterionic POPC, anionic POPG and POPS, and neutral POPE, as well as various heterogeneous lipid bilayers, to understand the effect of lipid composition on the antimicrobial peptide MSI-594 (GIGKFLKKAKKGIGAVLKVLTTG) using accelerated MD simulation, and confirmed that the lipid composition indeed affected the peptide–membrane interaction [165]. During internalization, the formation of pores is associated with a non-trivial free energy cost, Hu et al. used umbrella-sampling MD simulations with a lipid-density-based order parameter to investigate the membrane-pore-formation free energies of cyclic nona-arginine peptide Arg9 in 18 different lipid bilayer systems [166]. The lipid bilayer systems in the work are classified into three broad classes based on the length, the saturation of hydrocarbon acyl chain and the charge state of the lipid headgroup. The results showed that nona-arginine could induce transmembrane pores, and the free energy of peptide translocation along the pore pathways was relatively lower than that of the pore-free pathways. It was also revealed in the work that the pore-formation barriers were dependent on area per lipid, lipid bilayer thickness, and membrane bending rigidities. Niemel a ¨ et al. also studied the properties of three-component lipid bilayers, so-called lipid rafts, rich in cholesterol (CHOL) and sphingolipids using MD simulations [167]. Their work showed that the membrane elasticity and dynamic properties were strongly dependent on lipid composition due to the local interactions between lipid species. The lateral pressure profile of the membrane could be altered by changes in lipid content, which in turn regulated the function of some membrane proteins. Another work by Crosio et al. showed that the peptide accumulation on the membrane surface was maximal for the DOPC/CHOL/DOPG composition [61]. However, the peptides did not insert into the monolayer containing CHOL, indicating that although the anionic membrane concentrated CPPs on its surface, the reduced permeability prevented CPP translocation. The PMF calculated from umbrella sampling simulations using the coarse-grained force field revealed that the translocation free energy barriers decreased synchronously as DPPS composition increased due to the favorable peptide-anionic lipid interactions. Moreover, the addition of CHOL into the membrane increased the barriers of peptide translocation. CHOL not only coordinated with the peptide leading to a decrease in the peptide-lipid interactions, but also stabilized the liquid-ordered phase of membranes, which would suppress the reorientation of the lipid molecules and lead to the increase of the elastic stiffness of bilayers. Consequently, cholesterol hindered the formation of transmembrane pores, which was in consistent with experimental observations [168]. Similar conclusions were also obtained in MD simulations of the HIV-1 Tat peptide translocation in DPPC/DPPS membranes with the presence of 0–30 mol% cholesterol [169], that is, the cholesterol stabilized the liquid-ordered phase of the membrane and increased the elastic stiffness of the bilayer, thereby hindered the formation of transmembrane pores. Recent work by Tang et al. provides another perspective on the role of cholesterol [170]. The plasma membrane of cancer cells is rich in cholesterol, which impairs T-cell mediated cytotoxicity because cancer cells soften cortical structures through cholesterol enrichment in the plasma membrane. Cancer cells can be stiffened through cholesterol depletion, thus enhancing the cytotoxicity of T-cell and the efficacy of adoptive T-cell transfer (ACT) therapies. The results show the importance of the composition and mechanical properties of the plasma membrane, which are worthy of further study from the perspective of simulation. Polyansky et al. used MD simulations to investigate the spatial heterogeneous distribution of hydrophobic/hydrophilic “mosaic” character of the bilayer surface and their influence on the binding mode of penetratin [171]. They calculated and mapped the molecular hydrophobicity potential (MHP) which has been used to characterize the polarity of peptide–membrane interactions and found that the polar properties of the local environment of peptide residues at the interface could vary significantly, determining the detailed characteristics of the insertion. In anionic membranes, the initial adsorption of penetratin was strongly depended on the “complementarity” between the polar properties of the peptide and its local interfacial environment. If a high degree of complementarity was established, penetratin would penetrate deeply into the membrane without significantly disrupting its initial secondary structure. Such effects explained the complicated behaviors of CPPs, especially if the target membrane surface was of distinctly mosaic nature, and depending on the microscopic properties of the water–lipid interface, the peptides were able to adopt different pathways to exert their biological activities.
Since lipids in the outer membrane of eukaryotic cells are mainly zwitterionic, and the membranes have anion properties owing to the presence of glycosaminoglycans, for the choice of membrane models, most MD simulation studies on the interactions of CPPs with membranes employed simplified anionic or zwitterionic lipid bilayers in the past. For example, Crosio et al. investigated the effect of mechanical properties of anionic membranes with different fluidity and rigidity on the adsorption and penetration of poly-arginine [61]. With the improvement of computing methods and the development of computing resources, computational models of cell membranes have undergone a transition from simple models to multi-component models, such as mixed bilayers with different compositions [172,173]. Efforts are being made to realize more and more realistic models in which different groups in lipids can be distinguished [69], as well as models that can be used for larger systems. For example, Wang et.al proposed an implicit solvent coarse-grained model for quantitative simulations of POPC bilayers [174], by matching the structural and mechanical properties from experiments and all-atom bilayer simulations, the interactions were systematically adjusted to make the model more useful for studies of large-scale phenomena in membranes.

6.2. The Effect of Membrane Tension and Transmembrane Potential

Additional constraints or forces can also be applied to membranes to observe the effect on CPP penetration. For example, in many cases, osmotic pressure and external forces can stretch the membrane and position it under tension. He et al. used coarse-grained MD simulations to study the translocation of the polyarginine peptides Arg8 across the asymmetric human erythrocyte membranes, including DPPC, DPPE and DPPS lipids under tension [175]. The results showed that the membrane tension could inhibit the penetration of peptides, and with the increase of membrane tension, the inhibition effect strengthened as the membrane became thinner, the lipid tails were more disordered, and the potential of intertwining between neighboring lipid tails increased. Higher concentrations of CPPs or the application of a transmembrane potential could destabilize the membrane, thereby promoting the internalization of CPPs. Gao et al.’s MD simulation work provided evidence for the direct translocation of CPPs across membranes driven by membrane electrostatic potentials, which are ubiquitous in biological systems [176]. The local membrane potential could be generated by an imbalance in the concentration of transmembrane ions, and if positively charged CPPs adsorbed on the membrane, the membrane potential would further increase, resulting in the opening of membrane pores through which CPPs could be transported. The authors applied the classical nucleation theory to estimate the translocation time by calculating the changes in free energy when CPPs transferred across the membrane at different potentials, and the results were in good agreement with experimental measurements [177].

6.3. Peptide-Induced Membrane Response

Another issue of concern is the possible membrane response to CPP binding, including peptide-induced changes in membrane structure, such as membrane thickness, curvature, area per lipid, the lipid rearrangements, lateral diffusion constants of lipid molecules, and membrane mechanical performance. These parameters are important for understanding the peptide-lipid interactions. For example, the order parameter represents the mobility of the C-H bond of the lipid tails, and is calculated by averaging over time and all lipid tails. During interaction with CPPs, peptide-induced membrane responses could reflect some important aspects of the CPP penetration mechanism, such as the insertion of polar side chains into the hydrophobic core of the membrane to induce local membrane deformations [178,179], which had been thought to be a possible translocation mechanism for Arg-rich peptides. To understand how CPPs affect the mechanical properties of membranes, Grasso et al. used MD simulations to investigate the interactions between different types of CPPs embedded in the DOPC lipid bilayers [180]. The findings highlighted that the presence of CPPs increased local lipid disorder caused by the amount of water molecules conveyed by the peptides within the lipid bilayers, and the lipid disorder was directly related to the membrane stiffness. As a consequence, water penetration promoted by CPPs led to a local decrease in lipid order, which emerged macroscopically as a decrease in the modulus of membrane bending. Meanwhile, other simulations also found that the phenomenon of membrane bending is concentration dependent [138]. However, there’s still a lot that’s unclear for us to understand the membrane response to CPP binding, more research efforts need to be made in this field.

7. Conclusions and Future Perspectives

Cell-penetrating peptides have become promising candidates for next-generation intracellular drug delivery carriers due to their ability to facilitate the transmembrane transport of cargo molecules. However, how to improve the transmembrane efficiency and selective permeation of CPPs while limiting their cytotoxicity is a long-standing challenge for scientists. Considering the complicated physiological environment, the factors affecting the penetration process and internalization efficiency are not isolated. On one hand, the transmembrane efficiency depends on the choice of membrane model; on the other hand, the insertion of CPPs leads to changes in membrane thickness, curvature, etc., or the formation of hydrophilic pores, which in turn affects the internalization process. In addition, the effect of cargo molecules on the internalization is critical, and their physicochemical properties may greatly affect the performance of CPPs. In this review, the summarized computational methods have been shown to provide feasible pathways for elucidating the details of peptide–membrane interactions. All-atom MD simulations provide atomistic information on the intra- and inter-molecular interactions that drive peptides binding to the water-lipid interface, and characterize the cooperative mechanisms of membrane destabilization and the thermodynamics of these processes. Such details are especially complementary to experiments and are helpful in guiding the design of new experiments. We also introduced the research progress of MD simulation methods on CPP-membrane interactions in this review, including MD simulation methods and techniques, CPP penetration mechanisms, decoration or coupling of CPPs, the effect of different membrane models on the penetration process, and peptide-induced membrane reactions, as well as the comparison of simulation and experimental results. Investigating peptide–membrane interactions through MD simulations can let us focus on the influence of one or more specific factors, help us to understand the translocation mechanism at atomic scale, and provide new insights into designing more efficient translocation sequences.
For the mechanisms of CPP internalization, most of the simulation results can be explained by the pore formation model. Many CPPs are rich in positively charged arginine residues, and their guanidine groups can form electrostatic interactions, bidentate hydrogen-bonding interactions, and salt bridges with the phosphate groups of lipids. Such strong interactions destabilize the packing of the lipids and distort the membrane structure, and the charged amino acid residues will route water molecules into the hydrophobic core of the lipid bilayer, thereby inducing a hydrophilic pore in the membrane, and facilitating the transfer of CPPs. For some amphiphilic CPPs lacking arginine, peptides with the α-helical conformation have been shown to have better penetration when interacting with the cell membranes, especially in the pore formation model. In some simulations, the peptides are oriented parallel to the membrane surface, resulting in a mass imbalance across the bilayer, and this perturbation causes the peptides to sink deeper into the bilayer, which can be explained by the membrane thinning model. If we adopt larger systems with more lipids, it is hopefully to obtain richer responses, such as changes in membrane curvature. When the peptides and the lipid heads are strongly attracted, the membrane curvatures or invaginations may lead to the formation of inverted micelles. Some related simulations have shown that the aggregation of peptides induces the deformation of the lipid bilayer to form small vesicles that encapsulate the peptides, and the bilayer patches undergo large-scale deformation, which might be a potential microendocytosis mechanism. Computer simulations of endocytosis are still lacking to date due to the limitations of the temporal and spatial scales of simulations, as well as the fact that the endocytosis often requires the participation of membrane proteins.
More efforts are being made to improve the application availability of CPP simulations. Studies on the translocation mechanism of CPPs have shown that the secondary structures of CPPs and specific residues, such as positively charged arginine and lysine, as well as aromatic amino acids are of significant importance during uptake. Identifying intrinsic properties and conformational preferences of the peptides by MD simulations provides valuable clues for more flexible engineering design of novel CPPs. Minor modifications to the CPP structure may significantly alter its cellular uptake capacity, as this bottom-up approach enables regulation of the secondary structure of the peptide, making it promising for clinical applications such as cancer imaging and therapy. The optimization of the structure can be achieved by sequence recombination or amino acids substitution of the widely studied CPP templates. The affinity and stability of ligand docking and the transport efficiency of the complexes can also be improved by modifying the nature of the cargo and the CPP-cargo coupling mode, which is another effective way to design novel peptide-based drug delivery systems with high translocation efficiency and low biotoxicity. Studies have shown that the presence of cargoes affects the internalization pathway of CPPs, and the coupling with cargoes also alters the chemical characteristics of the CPPs, which in turn alters their interactions with the membrane.
In the end, some limitations should be addressed for the applications of MD simulations. Due to the limited accuracy of force fields, the difficulty in adequate sampling of the full energy landscape, the simplification of the membrane model, and the convergence problem, the interpretation of simulation results must be very cautious. The outcome of the simulations depends largely on the initial conformations and orientations of the peptides, the extent of equilibration, and the conditions under which the simulations are performed. It is also possible to miss the lowest free energy translocation path due to inappropriate CV selection, resulting in sampling of a local minimum energy region. Therefore, it is necessary to perform the simulations starting from the significant conformations obtained in experiments such as NMR, X-ray crystallography and cryo-electron microscopy, and to run multiple MD trajectories under independent start-up conditions to understand the average behavior, and the convergence of the simulations must be evaluated carefully. One possible approach is to incorporate experimental data into model development to make the simulation models more quantitative consistent with experiments. With the development of new sampling methods, more accurate force fields and membrane models, it is now promising, with the help of MD simulations, to extract the structure-activity relationship of CPPs, to adjust the nature of existing peptides and the coupling mode with cargo molecules, and to predict and generate bioactive CPPs with high transmembrane efficiency and cell selectivity.

Author Contributions

Conceptulizazion, Y.S. and W.W.; investigation, J.O. and Y.S.; writing—original draft preparation, J.O.; writing—review and editing, Y.S.; supervision, Y.S. and W.W.; funding acquisition, Y.S. and W.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (NSFC), grant numbers 11874197 and 11934008.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eiriksdottir, E.; Konate, K.; Langel, U.; Divita, G.; Deshayes, S. Secondary structure of cell-penetrating peptides controls membrane interaction and insertion. Biochim. Biophys. Acta-Biomembr. 2010, 1798, 1119–1128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Snyder, E.L.; Dowdy, S.F. Cell penetrating peptides in drug delivery. Pharm. Res. 2004, 21, 389–393. [Google Scholar] [CrossRef] [PubMed]
  3. Niesner, U.; Halin, C.; Lozzi, L.; Gunthert, M.; Neri, P.; Wunderli-Allenspach, H.; Zardi, L.; Neri, D. Quantitation of the tumor-targeting properties of antibody fragments conjugated to cell-permeating HIV-1 TAT peptides. Bioconjug. Chem. 2002, 13, 729–736. [Google Scholar] [CrossRef]
  4. Silva, S.; Almeida, A.J.; Vale, N. Combination of Cell-Penetrating Peptides with Nanoparticles for Therapeutic Application: A Review. Biomolecules 2019, 9, 22. [Google Scholar] [CrossRef] [Green Version]
  5. Tseng, Y.L.; Liu, J.J.; Hong, R.L. Translocation of liposomes into cancer cells by cell-penetrating peptides penetratin and TAT: A kinetic and efficacy study. Mol. Pharmacol. 2002, 62, 864–872. [Google Scholar] [CrossRef] [Green Version]
  6. Torchilin, V.P.; Levchenko, T.S.; Rammohan, R.; Volodina, N.; Papahadjopoulos-Sternberg, B.; D’Souza, G.G.M. Cell transfection in vitro and in vivo with nontoxic TAT peptide-liposome-DNA complexes. Proc. Natl. Acad. Sci. USA 2003, 100, 1972–1977. [Google Scholar] [CrossRef] [Green Version]
  7. Astriab-Fisher, A.; Sergueev, D.; Fisher, M.; Shaw, B.R.; Juliano, R.L. Conjugates of antisense oligonucleotides with the Tat and antennapedia cell-penetrating peptides: Effects on cellular uptake, binding to target sequences, and biologic actions. Pharm. Res. 2002, 19, 744–754. [Google Scholar] [CrossRef]
  8. Heitz, F.; Morris, M.C.; Divita, G. Twenty years of cell-penetrating peptides: From molecular mechanisms to therapeutics. Br. J. Pharmacol. 2009, 157, 195–206. [Google Scholar] [CrossRef] [Green Version]
  9. Borrelli, A.; Tornesello, A.L.; Tornesello, M.L.; Buonaguro, F.M. Cell Penetrating Peptides as Molecular Carriers for Anti-Cancer Agents. Molecules 2018, 23, 295. [Google Scholar] [CrossRef] [Green Version]
  10. Gessner, I.; Neundorf, I. Nanoparticles Modified with Cell-Penetrating Peptides: Conjugation Mechanisms, Physicochemical Properties, and Application in Cancer Diagnosis and Therapy. Int. J. Mol. Sci. 2020, 21, 2536. [Google Scholar] [CrossRef]
  11. Grau, M.; Walker, P.R.; Derouazi, M. Mechanistic insights into the efficacy of cell penetrating peptide-based cancer vaccines. Cell. Mol. Life Sci. 2018, 75, 2887–2896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Habault, J.; Poyet, J.L. Recent Advances in Cell Penetrating Peptide-Based Anticancer Therapies. Molecules 2019, 24, 927. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Nam, S.H.; Jang, J.; Cheon, D.H.; Chong, S.E.; Ahn, J.H.; Hyun, S.; Yu, J.; Lee, Y. pH-Activatable cell penetrating peptide dimers for potent delivery of anticancer drug to triple-negative breast cancer. J. Control. Release 2021, 330, 898–906. [Google Scholar] [CrossRef]
  14. Wu, J.; Li, J.S.; Wang, H.; Liu, C.B. Mitochondrial-targeted penetrating peptide delivery for cancer therapy. Expert Opin. Drug Deliv. 2018, 15, 951–964. [Google Scholar] [CrossRef]
  15. Fukuoka, Y.; Khafagy, E.; Goto, T.; Kamei, N.; Takayama, K.; Peppas, N.A.; Takeda-Morishita, M. Combination Strategy with Complexation Hydrogels and Cell-Penetrating Peptides for Oral Delivery of Insulin. Biol. Pharm. Bull. 2018, 41, 811–814. [Google Scholar] [CrossRef] [Green Version]
  16. Cheng, Y.S.; Chen, Z.T.; Liao, T.Y.; Lin, C.; Shen, H.C.H.; Wang, Y.H.; Chang, C.W.; Liu, R.S.; Chen, R.P.Y.; Tu, P.H. An intranasally delivered peptide drug ameliorates cognitive decline in Alzheimer transgenic mice. EMBO Mol. Med. 2017, 9, 703–715. [Google Scholar] [CrossRef]
  17. Tansi, F.L.; Filatova, M.P.; Koroev, D.O.; Volpina, O.M.; Lange, S.; Schumann, C.; Teichgraber, U.K.; Reissmann, S.; Hilger, I. New generation CPPs show distinct selectivity for cancer and noncancer cells. J. Cell. Biochem. 2019, 120, 6528–6541. [Google Scholar] [CrossRef] [PubMed]
  18. Movafegh, B.; Jalal, R.; Mohanunadi, Z.; Aldaghi, S.A. Poly-L-arginine: Enhancing Cytotoxicity and Cellular Uptake of Doxorubicin and Necrotic Cell Death. Anti-Cancer Agents Med. Chem. 2018, 18, 1448–1456. [Google Scholar] [CrossRef]
  19. Du, Y.Q.; Wang, L.; Wang, W.Y.; Guo, T.; Zhang, M.Z.; Zhang, P.; Zhang, Y.N.; Wu, K.J.; Li, A.X.; Wang, X.Y.; et al. Novel Application of Cell Penetrating R11 Peptide for Diagnosis of Bladder Cancer. J. Biomed. Nanotechnol. 2018, 14, 161–167. [Google Scholar] [CrossRef]
  20. Berneche, S.; Nina, M.; Roux, B. Molecular dynamics simulation of melittin in a dimyristoylphosphatidylcholine bilayer membrane. Biophys. J. 1998, 75, 1603–1618. [Google Scholar] [CrossRef]
  21. Biggin, P.C.; Sansom, M.S.P. Interactions of alpha-helices with lipid bilayers: A review of simulation studies. Biophys. Chem. 1999, 76, 161–183. [Google Scholar] [CrossRef] [PubMed]
  22. La Rocca, P.; Biggin, P.C.; Tieleman, D.P.; Sansom, M.S.P. Simulation studies of the interaction of antimicrobial peptides and lipid bilayers. Biochim. Biophys. Acta-Biomembr. 1999, 1462, 185–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Karplus, M.; McCammon, J.A. Molecular dynamics simulations of biomolecules. Nat. Struct. Biol. 2002, 9, 646–652. [Google Scholar] [CrossRef] [PubMed]
  24. Leontiadou, H.; Mark, A.E.; Marrink, S.J. Antimicrobial Peptides in Action. J. Am. Chem. Soc. 2006, 128, 12156–12161. [Google Scholar] [CrossRef] [Green Version]
  25. Trofimenko, E.; Grasso, G.; Heulot, M.; Chevalier, N.; Deriu, M.A.; Dubuis, G.; Arribat, Y.; Serulla, M.; Michel, S.; Vantomme, G.; et al. Genetic, cellular, and structural characterization of the membrane potential-dependent cell-penetrating peptide translocation pore. Elife 2021, 10, e69832. [Google Scholar] [CrossRef]
  26. Frankel, A.D.; Pabo, C.O. Cellular uptake of the tat protein from human immunodeficiency virus. Cell 1988, 55, 1189–1193. [Google Scholar] [CrossRef] [PubMed]
  27. Fawell, S.; Seery, J.; Daikh, Y.; Moore, C.; Chen, L.L.; Barsoum, P.J. Tat-mediated delivery of heterologous proteins into cells. Proc. Natl. Acad. Sci. USA 1994, 91, 664–668. [Google Scholar] [CrossRef] [Green Version]
  28. Torchilin, V.P. Tat peptide-mediated intracellular delivery of pharmaceutical nanocarriers. Adv. Drug Deliv. Rev. 2008, 60, 548–558. [Google Scholar] [CrossRef]
  29. Tao, F.; Johns, R.A. Tat-Mediated Peptide Intervention in Analgesia and Anesthesia. Drug Dev. Res. 2010, 71, 99–105. [Google Scholar] [CrossRef] [Green Version]
  30. Derossi, D.; Joliot, A.H.; Chassaing, G.; Prochiantz, A. The 3rd helix of the Antennapedia homeodomain translocates through biological-membranes. J. Biol. Chem. 1994, 269, 10444–10450. [Google Scholar] [CrossRef]
  31. Derossi, D.; Chassaing, G.; Prochiantz, A. Trojan peptides: The penetratin system for intracellular delivery. Trends Cell Biol. 1998, 8, 84–87. [Google Scholar] [CrossRef] [PubMed]
  32. Lindberg, M.; Biverstahl, H.; Graslund, A.; Maler, L. Structure and positioning comparison of two variants of penetratin in two different membrane mimicking systems by NMR. Eur. J. Biochem. 2003, 270, 3055–3063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Thoren, P.E.G.; Persson, D.; Karlsson, M.; Norden, B. The Antennapedia peptide penetratin translocates across lipid bilayers—The first direct observation. FEBS Lett. 2000, 482, 265–268. [Google Scholar] [CrossRef] [Green Version]
  34. Futaki, S.; Ohashi, W.; Suzuki, T.; Niwa, M.; Tanaka, S.; Ueda, K.; Harashima, H.; Sugiura, Y. Stearylated arginine-rich peptides: A new class of transfection systems. Bioconj. Chem. 2001, 12, 1005–1011. [Google Scholar] [CrossRef]
  35. Pooga, M.; Hallbrink, M.; Zorko, M.; Langel, U. Cell penetration by transportan. FASEB J. 1998, 12, 67–77. [Google Scholar] [CrossRef] [PubMed]
  36. Soomets, U.; Lindgren, M.; Gallet, X.; Hallbrink, M.; Elmquist, A.; Balaspiri, L.; Zorko, M.; Pooga, M.; Brasseur, R.; Langel, U. Deletion analogues of transportan. Biochim. Biophys. Acta-Biomembr. 2000, 1467, 165–176. [Google Scholar] [CrossRef] [Green Version]
  37. Pooga, M.; Langel, U. Classes of Cell-Penetrating Peptides. In Cell-Penetrating Peptides: Methods and Protocols, 2nd ed.; Langel, U., Ed.; Humana: Atlanta, GA, USA, 2015; Volume 1324, pp. 3–28. [Google Scholar]
  38. Crombez, L.; Aldrian-Herrada, G.; Konate, K.; Nguyen, Q.N.; McMaster, G.K.; Brasseur, R.; Heitz, F.; Divita, G. A New Potent Secondary Amphipathic Cell-penetrating Peptide for siRNA Delivery Into Mammalian Cells. Mol. Ther. 2009, 17, 95–103. [Google Scholar] [CrossRef]
  39. Oehlke, J.; Scheller, A.; Wiesner, B.; Krause, E.; Beyermann, M.; Klauschenz, E.; Melzig, M.; Bienert, M. Cellular uptake of an alpha-helical amphipathic model peptide with the potential to deliver polar compounds into the cell interior non-endocytically. Biochim. Biophys. Acta-Biomembr. 1998, 1414, 127–139. [Google Scholar] [CrossRef] [Green Version]
  40. Futaki, S.; Suzuki, T.; Ohashi, W.; Yagami, T.; Tanaka, S.; Ueda, K.; Sugiura, Y. Arginine-rich peptides. An abundant source of membrane-permeable peptides having potential as carriers for intracellular protein delivery. J. Biol. Chem. 2001, 276, 5836–5840. [Google Scholar] [CrossRef] [Green Version]
  41. Elmquist, A.; Lindgren, M.; Bartfai, T.; Langel, U. VE-cadherin-derived cell-penetrating peptide, pVEC, with carrier functions. Exp. Cell Res. 2001, 269, 237–244. [Google Scholar] [CrossRef]
  42. Ziegler, A. Thermodynamic studies and binding mechanisms of cell-penetrating peptides with lipids and glycosaminoglycans. Adv. Drug Deliv. Rev. 2008, 60, 580–597. [Google Scholar] [CrossRef] [PubMed]
  43. Morris, M.C.; Depollier, J.; Mery, J.; Heitz, F.; Divita, G. A peptide carrier for the delivery of biologically active proteins into mammalian cells. Nat. Biotechnol. 2001, 19, 1173–1176. [Google Scholar] [CrossRef]
  44. Lin, Y.Z.; Yao, S.Y.; Veach, R.A.; Torgerson, T.R.; Hawiger, J. Inhibition of nuclear translocation of transcription factor NF-kappa B by a synthetic peptide containing a cell membrane-permeable motif and nuclear localization sequence. J. Biol. Chem. 1995, 270, 14255–14258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Milletti, F. Cell-penetrating peptides: Classes, origin, and current landscape. Drug Discov. Today 2012, 17, 850–860. [Google Scholar] [CrossRef] [PubMed]
  46. Di Pisa, M.; Chassaing, G.; Swiecicki, J.M. Translocation Mechanism(s) of Cell-Penetrating Peptides: Biophysical Studies Using Artificial Membrane Bilayers. Biochemistry 2015, 54, 194–207. [Google Scholar] [CrossRef] [PubMed]
  47. Doherty, G.J.; McMahon, H.T. Mechanisms of Endocytosis. Annu. Rev. Biochem. 2009, 78, 857–902. [Google Scholar] [CrossRef] [Green Version]
  48. Kauffman, W.B.; Fuselier, T.; He, J.; Wimley, W.C. Mechanism Maters: A Taxonomy of Cell Penetrating Peptices. Trends Biochem. Sci. 2015, 40, 749–764. [Google Scholar] [CrossRef] [Green Version]
  49. Duchardt, F.; Fotin-Mleczek, M.; Schwarz, H.; Fischer, R.; Brock, R. A comprehensive model for the cellular uptake of cationic cell-penetrating peptides. Traffic 2007, 8, 848–866. [Google Scholar] [CrossRef]
  50. Ter-Avetisyan, G.; Tuennemann, G.; Nowak, D.; Nitschke, M.; Herrmann, A.; Drab, M.; Cardoso, M.C. Cell Entry of Arginine-rich Peptides Is Independent of Endocytosis. J. Biol. Chem. 2009, 284, 3370–3378. [Google Scholar] [CrossRef] [Green Version]
  51. Gazit, E.; Lee, W.J.; Brey, P.T.; Shai, Y.C. Mode of Action of the Antibacterial Cecropin B2: A Spectrofluorometric Study. Biochemistry 1994, 33, 10681–10692. [Google Scholar] [CrossRef]
  52. Herce, H.D.; Garcia, A.E. Molecular dynamics simulations suggest a mechanism for translocation of the HIV-1 TAT peptide across lipid membranes. Proc. Natl. Acad. Sci. USA 2007, 104, 20805–20810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Derossi, D.; Calvet, S.; Trembleau, A.; Brunissen, A.; Chassaing, G.; Prochiantz, A. Cell internalization of the third helix of the antennapedia homeodomain is receptor-independent. J. Biol. Chem. 1996, 271, 18188–18193. [Google Scholar] [CrossRef] [Green Version]
  54. Prochiantz, A. Getting hydrophilic compounds into cells: Lessons from homeopeptides—Commentary. Curr. Opin. Neurobiol. 1996, 6, 629–634. [Google Scholar] [CrossRef] [PubMed]
  55. Pouny, Y.; Rapaport, D.; Mor, A.; Nicolas, P.; Shai, Y. Interaction of Antimicrobial Dermaseptin and Its Fluorescently Labeled Analogues with Phospholipid Membranes. Biochemistry 1992, 31, 12416–12423. [Google Scholar] [CrossRef] [PubMed]
  56. Taylor, C.T.; Furuta, G.T.; Synnestvedt, K.; Colgan, S.P. Phosphorylation-dependent targeting of cAMP response element binding protein to the ubiquitin/proteasome pathway in hypoxia. Proc. Natl. Acad. Sci. USA 2000, 97, 12091–12096. [Google Scholar] [CrossRef] [Green Version]
  57. Lee, M.T.; Hung, W.C.; Chen, F.Y.; Huang, H.W. Many-body effect of antimicrobial peptides: On the correlation between lipid’s spontaneous curvature and pore formation. Biophys. J. 2005, 89, 4006–4016. [Google Scholar] [CrossRef] [Green Version]
  58. Wang, F.; Wang, Y.; Zhang, X.; Zhang, W.; Guo, S.; Jin, F. Recent progress of cell-penetrating peptides as new carriers for intracellular cargo delivery. J. Control. Release 2014, 174, 126–136. [Google Scholar] [CrossRef]
  59. Fretz, M.M.; Penning, N.A.; Al-Taei, S.; Futaki, S.; Takeuchi, T.; Nakase, I.; Storm, G.; Jones, A.T. Temperature-, concentration- and cholesterol-dependent translocation of L- and D-octa-arginine across the plasma and nuclear membrane of CD34(+) leukaemia cells. Biochem. J. 2007, 403, 335–342. [Google Scholar] [CrossRef] [Green Version]
  60. Mueller, J.; Kretzschmar, I.; Volkmer, R.; Boisguerin, P. Comparison of Cellular Uptake Using 22 CPPs in 4 Different Cell Lines. Bioconj. Chem. 2008, 19, 2363–2374. [Google Scholar] [CrossRef]
  61. Crosio, M.A.; Via, M.A.; Camara, C.I.; Mangiarotti, A.; Del Popolo, M.G.; Wilke, N. Interaction of a Polyarginine Peptide with Membranes of Different Mechanical Properties. Biomolecules 2019, 9, 625. [Google Scholar] [CrossRef]
  62. Tunnemann, G.; Martin, R.M.; Haupt, S.; Patsch, C.; Edenhofer, F.; Cardoso, M.C. Cargo-dependent mode of uptake and bioavailability of TAT-containing proteins and peptides in living cells. FASEB J. 2006, 20, 1775–1784. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Jones, A.T.; Sayers, E.J. Cell entry of cell penetrating peptides: Tales of tails wagging dogs. J. Control. Release 2012, 161, 582–591. [Google Scholar] [CrossRef] [PubMed]
  64. Alder, B.J.; Wainwright, T.E. Phase Transition for a Hard Sphere System. J. Chem. Phys. 1957, 27, 1208–1209. [Google Scholar] [CrossRef] [Green Version]
  65. Rahman, A. Correlations in the Motion of Atoms in Liquid Argon. Phys. Rev. 1964, 136, A405–A411. [Google Scholar] [CrossRef]
  66. Lensink, M.F. Membrane-Associated Proteins and Peptides. In Molecular Modeling of Proteins, 2nd ed.; Kukol, A., Ed.; Humana: New York, NY, USA, 2015; pp. 109–124. [Google Scholar]
  67. Hu, Y.; Liu, X.R.; Sinha, S.K.; Patel, S. Trans location Thermodynamics of Linear and Cyclic Nonaarginine into Model DPPC Bilayer via Coarse-Grained Molecular Dynamics Simulation: Implications of Pore Formation and Nonadditivity. J. Phys. Chem. B 2014, 118, 2670–2682. [Google Scholar] [CrossRef]
  68. Saunders, M.G.; Voth, G.A. Coarse-Graining Methods for Computational Biology. In Annual Review of Biophysics; Dill, K.A., Ed.; Annual Reviews: San Mateo, CA, USA, 2013; Volume 42, pp. 73–93. [Google Scholar]
  69. Marrink, S.J.; Corradi, V.; Souza, P.C.T.; Ingolfsson, H.I.; Tieleman, D.P.; Sansom, M.S.P. Computational Modeling of Realistic Cell Membranes. Chem. Rev. 2019, 119, 6184–6226. [Google Scholar] [CrossRef] [Green Version]
  70. Lazaridis, T.; Leveritt, J.M.; PeBenito, L. Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers. Biochim. Biophys. Acta-Biomembr. 2014, 1838, 2149–2159. [Google Scholar] [CrossRef] [Green Version]
  71. Nepal, B.; Leveritt, J.; Lazaridis, T. Membrane Curvature Sensing by Amphipathic Helices: Insights from Implicit Membrane Modeling. Biophys. J. 2018, 114, 2128–2141. [Google Scholar] [CrossRef] [Green Version]
  72. Park, S.; Schulten, K. Calculating potentials of mean force from steered molecular dynamics simulations. J. Chem. Phys. 2004, 120, 5946–5961. [Google Scholar] [CrossRef]
  73. Yesylevskyy, S.; Marrink, S.J.; Mark, A.E. Alternative Mechanisms for the Interaction of the Cell-Penetrating Peptides Penetratin and the TAT Peptide with Lipid Bilayers. Biophys. J. 2009, 97, 40–49. [Google Scholar] [CrossRef]
  74. Torrie, G.M.; Valleau, J.P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys. 1977, 23, 187–199. [Google Scholar] [CrossRef]
  75. Kstner, J. Umbrella sampling. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 932–942. [Google Scholar] [CrossRef]
  76. Laio, A.; Parrinello, M. Escaping free-energy minima. Proc. Natl. Acad. Sci. USA 2002, 99, 12562–12566. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Barducci, A.; Bonomi, M.; Parrinello, M. Metadynamics. Wiley Interdiscip. Rev.-Comput. Mol. Sci. 2011, 1, 826–843. [Google Scholar] [CrossRef]
  78. Sugita, Y.; Okamoto, Y. Replica-exchange molecular dynamics method for protein folding. Chem. Phys. Lett. 1999, 314, 141–151. [Google Scholar] [CrossRef]
  79. Reid, L.M.; Verma, C.S.; Essex, J.W. The role of molecular simulations in understanding the mechanisms of cell-penetrating peptides. Drug Discov. Today 2019, 24, 1821–1835. [Google Scholar] [CrossRef]
  80. Nymeyer, H.; Woolf, T.B.; Garcia, A.E. Folding is not required for bilayer insertion: Replica exchange simulations of an alpha-helical peptide with an explicit lipid bilayer. Proteins-Struct. Funct. Bioinform. 2005, 59, 783–790. [Google Scholar] [CrossRef]
  81. Alaybeyoglu, B.; Akbulut, B.S.; Ozkirimli, E. pVEC hydrophobic N-terminus is critical for antibacterial activity. J. Pept. Sci. 2018, 24, e3083. [Google Scholar] [CrossRef]
  82. Babin, V.; Roland, C.; Sagui, C. Adaptively biased molecular dynamics for free energy calculations. J. Chem. Phys. 2008, 128, 134101. [Google Scholar] [CrossRef] [Green Version]
  83. Gimenez-Dejoz, J.; Numata, K. Molecular dynamics study of the internalization of cell-penetrating peptides containing unnatural amino acids across membranes. Nanoscale Adv. 2022, 4, 397–407. [Google Scholar] [CrossRef]
  84. Noe, F.; Fischer, S. Transition networks for modeling the kinetics of conformational change in macromolecules. Curr. Opin. Struct. Biol. 2008, 18, 154–162. [Google Scholar] [CrossRef] [PubMed]
  85. Chodera, J.D.; Swope, W.C.; Pitera, J.W.; Dill, K.A. Long-time protein folding dynamics from short-time molecular dynamics simulations. Multiscale Model. Simul. 2006, 5, 1214–1226. [Google Scholar] [CrossRef] [Green Version]
  86. Prinz, J.-H.; Wu, H.; Sarich, M.; Keller, B.; Senne, M.; Held, M.; Chodera, J.D.; Schuette, C.; Noe, F. Markov models of molecular kinetics: Generation and validation. J. Chem. Phys. 2011, 134, 174105. [Google Scholar] [CrossRef] [PubMed]
  87. Qiao, Q.; Bowman, G.R.; Huang, X. Dynamics of an Intrinsically Disordered Protein Reveal Metastable Conformations That Potentially Seed Aggregation. J. Am. Chem. Soc. 2013, 135, 16092–16101. [Google Scholar] [CrossRef] [PubMed]
  88. Noe, F.; Schuette, C.; Vanden-Eijnden, E.; Reich, L.; Weikl, T.R. Constructing the equilibrium ensemble of folding pathways from short off-equilibrium simulations. Proc. Natl. Acad. Sci. USA 2009, 106, 19011–19016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  89. Wang, D.; Weng, J.; Wang, W. An unconventional ligand-binding mechanism of substrate-binding proteins: MD simulation and Markov state model analysis of BtuF. J. Comput. Chem. 2019, 40, 1440–1448. [Google Scholar] [CrossRef]
  90. Sengupta, U.; Carballo-Pacheco, M.; Strodel, B. Automated Markov state models for molecular dynamics simulations of aggregation and self-assembly. J. Chem. Phys. 2019, 150, 115101. [Google Scholar] [CrossRef]
  91. Weng, J.; Yang, M.; Wang, W.; Xu, X.; Tian, Z. Revealing Thermodynamics and Kinetics of Lipid Self-Assembly by Markov State Model Analysis. J. Am. Chem. Soc. 2020, 142, 21344–21352. [Google Scholar] [CrossRef]
  92. Huang, K.; Garcia, A.E. Free Energy of Translocating an Arginine-Rich Cell-Penetrating Peptide across a Lipid Bilayer Suggests Pore Formation. Biophys. J. 2013, 104, 412–420. [Google Scholar] [CrossRef] [Green Version]
  93. Kabelka, I.; Brozek, R.; Vacha, R. Selecting Collective Variables and Free-Energy Methods for Peptide Translocation across Membranes. J. Chem. Inf. Model. 2021, 61, 819–830. [Google Scholar] [CrossRef]
  94. Madani, F.; Lindberg, S.; Langel, Ü.; Futaki, S.; Gräslund, A. Mechanisms of Cellular Uptake of Cell-Penetrating Peptides. J. Biophys. 2011, 2011, 414729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Sengupta, D.; Leontiadou, H.; Mark, A.E.; Marrink, S.-J. Toroidal pores formed by antimicrobial peptides show significant disorder. Biochim. Biophys. Acta-Biomembr. 2008, 1778, 2308–2317. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Herce, H.D.; Garcia, A.E.; Cardoso, M.C. Fundamental Molecular Mechanism for the Cellular Uptake of Guanidinium-Rich Molecules. J. Am. Chem. Soc. 2014, 136, 17459–17467. [Google Scholar] [CrossRef]
  97. Mishra, A.; Lai, G.H.; Schmidt, N.W.; Sun, V.Z.; Rodriguez, A.R.; Tong, R.; Tang, L.; Cheng, J.J.; Deming, T.J.; Kamei, D.T.; et al. Translocation of HIV TAT peptide and analogues induced by multiplexed membrane and cytoskeletal interactions. Proc. Natl. Acad. Sci. USA 2011, 108, 16883–16888. [Google Scholar] [CrossRef] [Green Version]
  98. Chen, X.C.; Sa’adedin, F.; Deme, B.; Rao, P.F.; Bradshaw, J. Insertion of TAT peptide and perturbation of negatively charged model phospholipid bilayer revealed by neutron diffraction. Biochim. Biophys. Acta-Biomembr. 2013, 1828, 1982–1988. [Google Scholar] [CrossRef] [Green Version]
  99. Pourmousa, M.; Karttunen, M. Early stages of interactions of cell-penetrating peptide penetratin with a DPPC bilayer. Chem. Phys. Lipids 2013, 169, 85–94. [Google Scholar] [CrossRef]
  100. Yao, C.; Kang, Z.Z.; Yu, B.; Chen, Q.; Liu, Y.C.; Wang, Q. All-Factor Analysis and Correlations on the Transmembrane Process for Arginine-Rich Cell-Penetrating Peptides. Langmuir 2019, 35, 9286–9296. [Google Scholar] [CrossRef]
  101. Mitchell, D.J.; Kim, D.T.; Steinman, L.; Fathman, C.G.; Rothbard, J.B. Polyarginine enters cells more efficiently than other polycationic homopolymers. J. Pept. Res. 2000, 56, 318–325. [Google Scholar] [CrossRef]
  102. Magzoub, M.; Kilk, K.; Eriksson, L.E.; Langel, U.; Graslund, A. Interaction and structure induction of cell-penetrating peptides in the presence of phospholipid vesicles. Biochim. Biophys. Acta 2001, 1512, 77–89. [Google Scholar] [CrossRef] [Green Version]
  103. Deshayes, S.; Heitz, A.; Morris, M.C.; Charnet, P.; Divita, G.; Heitz, F. Insight into the mechanism of internalization of the cell-penetrating carrier peptide Pep-1 through conformational analysis. Biochemistry 2004, 43, 1449–1457. [Google Scholar] [CrossRef]
  104. Dunkin, C.M.; Pokorny, A.; Almeida, P.F.; Lee, H.S. Molecular Dynamics Studies of Transportan 10 (Tp10) Interacting with a POPC Lipid Bilayer. J. Phys. Chem. B 2011, 115, 1188–1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  105. Su, Y.C.; Waring, A.J.; Ruchala, P.; Hong, M. Membrane-Bound Dynamic Structure of an Arginine-Rich Cell-Penetrating Peptide, the Protein Transduction Domain of HIV TAT, from Solid-State NMR. Biochemistry 2010, 49, 6009–6020. [Google Scholar] [CrossRef] [Green Version]
  106. Song, J.J.; Kai, M.; Zhang, W.; Zhang, J.D.; Liu, L.W.; Zhang, B.Z.; Liu, X.; Wang, R. Cellular uptake of transportan 10 and its analogs in live cells: Selectivity and structure-activity relationship studies. Peptides 2011, 32, 1934–1941. [Google Scholar] [CrossRef]
  107. Polyansky, A.A.; Volynsky, P.E.; Arseniev, A.S.; Efremov, R.G. Adaptation of a Membrane-active Peptide to Heterogeneous Environment. I. Structural Plasticity of the Peptide. J. Phys. Chem. B 2009, 113, 1107–1119. [Google Scholar] [CrossRef]
  108. Tsou, L.K.; Tatko, C.D.; Waters, M.L. Simple cation-pi interaction between a phenyl ring and a protonated amine stabilizes an alpha-helix in water. J. Am. Chem. Soc. 2002, 124, 14917–14921. [Google Scholar] [CrossRef]
  109. Czajlik, A.; Mesko, E.; Penke, B.; Perczel, A. Investigation of penetratin peptides Part 1. The environment dependent conformational properties of penetratin and two of its derivatives. J. Pept. Sci. 2002, 8, 151–171. [Google Scholar] [CrossRef]
  110. Balayssac, S.; Burlina, F.; Convert, O.; Bolbach, G.; Chassaing, G.; Lequin, O. Comparison of penetratin and other homeodomain-derived cell-penetrating peptides: Interaction in a membrane-mimicking environment and cellular uptake efficiency. Biochemistry 2006, 45, 1408–1420. [Google Scholar] [CrossRef]
  111. Jobin, M.L.; Vamparys, L.; Deniau, R.; Grelard, A.; Mackereth, C.D.; Fuchs, P.F.J.; Alves, I.D. Biophysical Insight on the Membrane Insertion of an Arginine-Rich Cell-Penetrating Peptide. Int. J. Mol. Sci. 2019, 20, 4441. [Google Scholar] [CrossRef] [Green Version]
  112. Rothbard, J.B.; Jessop, T.C.; Wender, P.A. Adaptive translocation: The role of hydrogen bonding and membrane potential in the uptake of guanidinium-rich transporters into cells. Adv. Drug Deliv. Rev. 2005, 57, 495–504. [Google Scholar] [CrossRef]
  113. Tang, M.; Waring, A.J.; Hong, M. Phosphate-mediated arginine insertion into lipid membranes and pore formation by a cationic membrane peptide from solid-state NMR. J. Am. Chem. Soc. 2007, 129, 11438–11446. [Google Scholar] [CrossRef]
  114. Lamaziere, A.; Burlina, F.; Wolf, C.; Chassaing, G.; Trugnan, G.; Ayala-Sanmartin, J. Non-Metabolic Membrane Tubulation and Permeability Induced by Bioactive Peptides. PLoS ONE 2007, 2, e201. [Google Scholar] [CrossRef] [PubMed]
  115. Smith, J.S.; Scholtz, J.M. Energetics of polar side-chain interactions in helical peptides: Salt effects on ion pairs and hydrogen bonds. Biochemistry 1998, 37, 33–40. [Google Scholar] [CrossRef] [PubMed]
  116. Andrew, C.D.; Penel, S.; Jones, G.R.; Doig, A.J. Stabilizing nonpolar/polar side-chain interactions in the alpha-helix. Proteins-Struct. Funct. Genet. 2001, 45, 449–455. [Google Scholar] [CrossRef]
  117. Walrant, A.; Vogel, A.; Correia, I.; Lequin, O.; Olausson, B.E.S.; Desbat, B.; Sagan, S.; Alves, I.D. Membrane interactions of two arginine-rich peptides with different cell internalization capacities. Biochim. Biophys. Acta-Biomembr. 2012, 1818, 1755–1763. [Google Scholar] [CrossRef]
  118. He, X.C.; Qu, Z.G.; Xu, F. Simulation study of interaction mechanism between peptide and asymmetric membrane. Mol. Simul. 2017, 43, 34–41. [Google Scholar] [CrossRef]
  119. Burns, K.E.; Robinson, M.K.; Thevenin, D. Inhibition of Cancer Cell Proliferation and Breast Tumor Targeting of pHLIP-Monomethyl Auristatin E Conjugates. Mol. Pharm. 2015, 12, 1250–1258. [Google Scholar] [CrossRef] [Green Version]
  120. Viola-Villegas, N.T.; Carlin, S.D.; Ackerstaff, E.; Sevak, K.K.; Divilov, V.; Serganova, I.; Kruchevsky, N.; Anderson, M.; Blasberg, R.G.; Andreev, O.A.; et al. Understanding the pharmacological properties of a metabolic PET tracer in prostate cancer. Proc. Natl. Acad. Sci. USA 2014, 111, 7254–7259. [Google Scholar] [CrossRef] [Green Version]
  121. Andreev, O.A.; Karabadzhak, A.G.; Weerakkody, D.; Andreev, G.O.; Engelman, D.M.; Reshetnyak, Y.K. pH (low) insertion peptide (pHLIP) inserts across a lipid bilayer as a helix and exits by a different path. Proc. Natl. Acad. Sci. USA 2010, 107, 4081–4086. [Google Scholar] [CrossRef] [Green Version]
  122. Deng, Y.H.; Qian, Z.Y.; Luo, Y.; Zhang, Y.; Mu, Y.G.; Wei, G.H. Membrane Binding and Insertion of a pHLIP Peptide Studied by All-Atom Molecular Dynamics Simulations. Int. J. Mol. Sci. 2013, 14, 14532–14549. [Google Scholar] [CrossRef] [Green Version]
  123. Ramaker, K.; Henkel, M.; Krause, T.; Roeckendorf, N.; Frey, A. Cell penetrating peptides: A comparative transport analysis for 474 sequence motifs. Drug Deliv. 2018, 25, 928–937. [Google Scholar] [CrossRef]
  124. Kumara, B.T.; Wijesiri, N.K.; Rathnayake, P.; Ranatunga, R. A Re-evaluation of the Free Energy Profiles for Cell-Penetrating Peptides across DOPC Membranes. Int. J. Pept. Res. Ther. 2021, 27, 2931–2943. [Google Scholar] [CrossRef]
  125. Ruzza, P.; Biondi, B.; Marchiani, A.; Antolini, N.; Calderan, A. Cell-Penetrating Peptides: A Comparative Study on Lipid Affinity and Cargo Delivery Properties. Pharmaceuticals 2010, 3, 1045–1062. [Google Scholar] [CrossRef] [Green Version]
  126. Langel, U. Cell-Penetrating Peptides and Transportan. Pharmaceutics 2021, 13, 987. [Google Scholar] [CrossRef]
  127. Kawamoto, S.; Takasu, M.; Miyakawa, T.; Morikawa, R.; Oda, T.; Futaki, S.; Nagao, H. Inverted micelle formation of cell-penetrating peptide studied by coarse-grained simulation: Importance of attractive force between cell-penetrating peptides and lipid head group. J. Chem. Phys. 2011, 134, 095103. [Google Scholar] [CrossRef]
  128. Afonin, S.; Frey, A.; Bayerl, S.; Fischer, D.; Wadhwani, P.; Weinkauf, S.; Ulrich, A.S. The cell-penetrating peptide TAT(48–60) induces a non-lamellar phase in DMPC membranes. Chemphyschem 2006, 7, 2134–2142. [Google Scholar] [CrossRef]
  129. Copolovici, D.M.; Langel, K.; Eriste, E.; Langel, U. Cell-Penetrating Peptides: Design, Synthesis, and Applications. ACS Nano 2014, 8, 1972–1994. [Google Scholar] [CrossRef]
  130. Li, L.B.; Vorobyov, I.; Allen, T.W. The Different Interactions of Lysine and Arginine Side Chains with Lipid Membranes. J. Phys. Chem. B 2013, 117, 11906–11920. [Google Scholar] [CrossRef]
  131. Shimonishi, Y.; Oehlke, J.; Scheller, A.; Janek, K.; Bienert, M. Rapid Translocation of Amphipathic α-Helical and β-Sheet-Forming Peptides through Plasma Membranes of Endothelial Cells; Peptide Science: Present and Future; Springer: Berlin/Heidelberg, Germany, 2002. [Google Scholar]
  132. Scheller, A.; Oehlke, J.; Wiesner, B.; Dathe, M.; Krause, E.; Beyermann, M.; Melzig, M.; Bienert, M. Structural requirements for cellular uptake of alpha-helical amphipathic peptides. J. Pept. Sci. 1999, 5, 185–194. [Google Scholar] [CrossRef]
  133. Martin, I.; Teixido, M.; Giralt, E. Design, Synthesis and Characterization of a New Anionic Cell-Penetrating Peptide: SAP(E). Chembiochem 2011, 12, 896–903. [Google Scholar] [CrossRef]
  134. Khafagy, E.S.; Morishita, M.; Ida, N.; Nishio, R.; Isowa, K.; Takayama, K. Structural requirements of penetratin absorption enhancement efficiency for insulin delivery. J. Control. Release 2010, 143, 302–310. [Google Scholar] [CrossRef]
  135. Kurrikoff, K.; Langel, U. Recent CPP-based applications in medicine. Expert Opin. Drug Deliv. 2019, 16, 1183–1191. [Google Scholar] [CrossRef] [PubMed]
  136. Xie, J.; Bi, Y.; Zhang, H.; Dong, S.; Teng, L.; Lee, R.J.; Yang, Z. Cell-Penetrating Peptides in Diagnosis and Treatment of Human Diseases: From Preclinical Research to Clinical Application. Front. Pharmacol. 2020, 11, 697. [Google Scholar] [CrossRef] [PubMed]
  137. Farrera-Sinfreu, J.; Giralt, E.; Royo, M.; Albericio, F. Cell-Penetrating Proline-Rich Peptidomimetics. In Peptide Characterization and Application Protocols; Fields, G.B., Ed.; Humana: Atlanta, GA, USA, 2007; Volume 386, pp. 241–267. [Google Scholar]
  138. Antunes, E.; Azoia, N.G.; Matama, T.; Gomes, A.C.; Cavaco-Paulo, A. The activity of LE10 peptide on biological membranes using molecular dynamics, in vitro and in vivo studies. Colloids Surf. B-Biointerfaces 2013, 106, 240–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  139. Tsai, C.W.; Lin, Z.W.; Chang, W.F.; Chen, Y.F.; Hu, W.W. Development of an indolicidin-derived peptide by reducing membrane perturbation to decrease cytotoxicity and maintain gene delivery ability. Colloids Surf. B-Biointerfaces 2018, 165, 18–27. [Google Scholar] [CrossRef]
  140. Karle, I.L. Controls exerted by the Aib residue: Helix formation and helix reversal. Biopolymers 2001, 60, 351–365. [Google Scholar] [CrossRef]
  141. Fujita, S.; Motoda, Y.; Kigawa, T.; Tsuchiya, K.; Numata, K. Peptide-Based Polyion Complex Vesicles That Deliver Enzymes into Intact Plants to Provide Antibiotic Resistance without Genetic Modification. Biomacromolecules 2021, 22, 1080–1090. [Google Scholar] [CrossRef]
  142. Tuennemann, G.; Ter-Avetisyan, G.; Martin, R.M.; Stoeckl, M.; Herrmann, A.; Cardoso, C. Live-cell analysis of cell penetration ability and toxicity of oligo-arginines. J. Pept. Sci. 2008, 14, 469–476. [Google Scholar] [CrossRef]
  143. Bode, S.A.; Wallbrecher, R.; Brock, R.; van Hest, J.C.M.; Lowik, D.W.P.M. Activation of cell-penetrating peptides by disulfide bridge formation of truncated precursors. Chem. Commun. 2014, 50, 415–417. [Google Scholar] [CrossRef]
  144. Oskolkov, N.; Arukuusk, P.; Copolovici, D.-M.; Lindberg, S.; Margus, H.; Padari, K.; Pooga, M.; Langel, U. NickFects, Phosphorylated Derivatives of Transportan 10 for Cellular Delivery of Oligonucleotides. Int. J. Pept. Res. Ther. 2011, 17, 147–157. [Google Scholar] [CrossRef]
  145. Zhang, L.L.; Hoffman, J.A.; Ruoslahti, E. Molecular profiling of heart endothelial cells. Circulation 2005, 112, 1601–1611. [Google Scholar] [CrossRef]
  146. Stewart, K.M.; Horton, K.L.; Kelley, S.O. Cell-penetrating peptides as delivery vehicles for biology and medicine. Org. Biomol. Chem. 2008, 6, 2242–2255. [Google Scholar] [CrossRef] [PubMed]
  147. Saini, A.; Jaswal, R.R.; Negi, R.; Nandel, F.S. Insights on the structural characteristics of Vim-TBS (58-81) peptide for future applications as a cell penetrating peptide. Biosci. Trends 2013, 7, 209–220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  148. Balzeau, J.; Peterson, A.; Eyer, J. The vimentin-tubulin binding site peptide (Vim-TBS.58-81) crosses the plasma membrane and enters the nuclei of human glioma cells. Int. J. Pharm. 2012, 423, 77–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  149. Ashwanikumar, N.; Plaut, J.S.; Mostofian, B.; Patel, S.; Kwak, P.; Sun, C.; McPhail, K.; Zuckerman, D.M.; Esener, S.C.; Sahay, G. Supramolecular self assembly of nanodrill-like structures for intracellular delivery. J. Control. Release 2018, 282, 76–89. [Google Scholar] [CrossRef]
  150. Garcia-Sosa, A.T.; Tulp, I.; Langel, K.; Langel, U. Peptide-Ligand Binding Modeling of siRNA with Cell-Penetrating Peptides. Biomed Res. Int. 2014, 2014, 257040. [Google Scholar] [CrossRef] [Green Version]
  151. Oh, D.; Shirazi, A.N.; Northup, K.; Sullivan, B.; Tiwari, R.K.; Bisoffi, M.; Parang, K. Enhanced Cellular Uptake of Short Polyarginine Peptides through Fatty Acylation and Cyclization. Mol. Pharm. 2014, 11, 2845–2854. [Google Scholar] [CrossRef]
  152. Hedegaard, S.F.; Bruhn, D.S.; Khandelia, H.; Cardenas, M.; Nielsen, H.M. Shuffled lipidation pattern and degree of lipidation determines the membrane interaction behavior of a linear cationic membrane-active peptide. J. Colloid Interface Sci. 2020, 578, 584–597. [Google Scholar] [CrossRef]
  153. Akhunzada, M.J.; Chandramouli, B.; Bhattacharjee, N.; Macchi, S.; Cardarelli, F.; Brancato, G. The role of Tat peptide self-aggregation in membrane pore stabilization: Insights from a computational study. Phys. Chem. Chem. Phys. 2017, 19, 27603–27610. [Google Scholar] [CrossRef]
  154. Hu, J.M.; Lou, Y.M.; Wu, F.M. Improved Intracellular Delivery of Polyarginine Peptides with Cargoes. J. Phys. Chem. B 2019, 123, 2636–2644. [Google Scholar] [CrossRef]
  155. Grasso, G.; Deriu, M.A.; Prat, M.; Rimondini, L.; Verne, E.; Follenzi, A.; Danani, A. Cell Penetrating Peptide Adsorption on Magnetite and Silica Surfaces: A Computational Investigation. J. Phys. Chem. B 2015, 119, 8239–8246. [Google Scholar] [CrossRef]
  156. Kubiak-Ossowska, K.; Burley, G.; Patwardhan, S.V.; Mulheran, P.A. Spontaneous Membrane-Translocating Peptide Adsorption at Silica Surfaces: A Molecular Dynamics Study. J. Phys. Chem. B 2013, 117, 14666–14675. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  157. Quan, X.B.; Sun, D.L.; Zhou, J. Molecular mechanism of HIV-1 TAT peptide and its conjugated gold nanoparticles translocating across lipid membranes. Phys. Chem. Chem. Phys. 2019, 21, 10300–10310. [Google Scholar] [CrossRef]
  158. Zhang, S.G.; Holmes, T.C.; Dipersio, C.M.; Hynes, R.O.; Su, X.; Rich, A. Self-complementary oligopeptide matrices support mammalian cell attachment. Biomaterials 1995, 16, 1385–1393. [Google Scholar] [CrossRef] [PubMed]
  159. Sheng, Y.B.; Wang, W.; Chen, P. Interaction of an ionic complementary peptide with a hydrophobic graphite surface. Protein Sci. 2010, 19, 1639–1648. [Google Scholar] [CrossRef] [Green Version]
  160. Sheng, Y.B.; Wang, W.; Chen, P. Adsorption of an Ionic Complementary Peptide on the Hydrophobic Graphite Surface. J. Phys. Chem. C 2010, 114, 454–459. [Google Scholar] [CrossRef]
  161. Hong, Y.S.; Legge, R.L.; Zhang, S.; Chen, P. Effect of amino acid sequence and pH on nanofiber formation of self-assembling peptides EAK16-II and EAK16-IV. Biomacromolecules 2003, 4, 1433–1442. [Google Scholar] [CrossRef]
  162. Fung, S.Y.; Yang, H.; Sadatmousavi, P.; Sheng, Y.; Mamo, T.; Nazarian, R.; Chen, P. Amino Acid Pairing for De Novo Design of Self-Assembling Peptides and Their Drug Delivery Potential. Adv. Funct. Mater. 2011, 21, 2456–2464. [Google Scholar] [CrossRef]
  163. Zhang, L.; Sheng, Y.B.; Yazdi, A.Z.; Sarikhani, K.; Wang, F.; Jiang, Y.S.; Liu, J.W.; Zheng, T.; Wang, W.; Ouyang, P.K.; et al. Surface-assisted assembly of a histidine-rich lipidated peptide for simultaneous exfoliation of graphite and functionalization of graphene nanosheets. Nanoscale 2019, 11, 2999–3012. [Google Scholar] [CrossRef]
  164. Lu, S.; Cui, W.J.; Li, J.; Sheng, Y.B.; Chen, P. Functional Control of Peptide Amphiphile Assemblies via Modulation of Internal Cohesion and Surface Chemistry Switch. Chem. Eur. J. 2018, 24, 13931–13937. [Google Scholar] [CrossRef]
  165. Mukherjee, S.; Kar, R.K.; Nanga, R.P.R.; Mroue, K.H.; Ramamoorthy, A.; Bhunia, A. Accelerated molecular dynamics simulation analysis of MSI-594 in a lipid bilayer. Phys. Chem. Chem. Phys. 2017, 19, 19289–19299. [Google Scholar] [CrossRef]
  166. Hu, Y.; Sinha, S.K.; Patel, S. Investigating Hydrophilic Pores in Model Lipid Bilayers Using Molecular Simulations: Correlating Bilayer Properties with Pore-Formation Thermodynamics. Langmuir 2015, 31, 6615–6631. [Google Scholar] [CrossRef] [Green Version]
  167. Niemela, P.S.; Ollila, S.; Hyvonen, M.T.; Karttunen, M.; Vattulainen, I. Assessing the nature of lipid raft membranes. PLoS Comput. Biol. 2007, 3, e34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  168. Sharmin, S.; Islam, M.Z.; Karal, M.A.; Shibly, S.U.; Dohra, H.; Yamazaki, M. Effects of Lipid Composition on the Entry of Cell-Penetrating Peptide Oligoarginine into Single Vesicles. Biochemistry 2016, 55, 4154–4165. [Google Scholar] [CrossRef] [PubMed]
  169. Hu, Y.; Patel, S. Thermodynamics of cell-penetrating HIV1 TAT peptide insertion into PC/PS/CHOL model bilayers through transmembrane pores: The roles of cholesterol and anionic lipids. Soft Matter 2016, 12, 6716–6727. [Google Scholar] [CrossRef] [PubMed]
  170. Lei, K.W.; Kurum, A.; Kaynak, M.; Bonati, L.; Han, Y.L.; Cencen, V.; Gao, M.; Xie, Y.Q.; Guo, Y.G.; Hannebelle, M.T.M.; et al. Cancer-cell stiffening via cholesterol depletion enhances adoptive T-cell immunotherapy. Nat. Biomed. Eng. 2021, 5, 1411–1425. [Google Scholar] [CrossRef]
  171. Polyansky, A.A.; Volynsky, P.E.; Arseniev, A.S.; Efremov, R.G. Adaptation of a Membrane-active Peptide to Heterogeneous Environment. II. The Role of Mosaic Nature of the Membrane Surface. J. Phys. Chem. B 2009, 113, 1120–1126. [Google Scholar] [CrossRef]
  172. Orsi, M.; Essex, J.W. Physical properties of mixed bilayers containing lamellar and nonlamellar lipids: Insights from coarse-grain molecular dynamics simulations. Faraday Discuss. 2013, 161, 249–272. [Google Scholar] [CrossRef]
  173. Ding, W.; Palaiokostas, M.; Wang, W.; Orsi, M. Effects of Lipid Composition on Bilayer Membranes Quantified by All-Atom Molecular Dynamics. J. Phys. Chem. B 2015, 119, 15263–15274. [Google Scholar] [CrossRef] [Green Version]
  174. Wang, Z.J.; Deserno, M. A Systematically Coarse-Grained Solvent-Free Model for Quantitative Phospholipid Bilayer Simulations. J. Phys. Chem. B 2012, 116, 3907. [Google Scholar] [CrossRef]
  175. He, X.C.; Lin, M.; Sha, B.Y.; Feng, S.S.; Shi, X.H.; Qu, Z.G.; Xu, F. Coarse-grained molecular dynamics studies of the translocation mechanism of polyarginines across asymmetric membrane under tension. Sci. Rep. 2015, 5, 12808. [Google Scholar] [CrossRef]
  176. Gao, X.L.; Hong, S.; Liu, Z.P.; Yue, T.T.; Dobnikar, J.; Zhang, X.R. Membrane potential drives direct translocation of cell-penetrating peptides. Nanoscale 2019, 11, 1949–1958. [Google Scholar] [CrossRef]
  177. Rao, Y.; Kwok, S.J.J.; Lombardi, J.; Turro, N.J.; Eisenthal, K.B. Label-free probe of HIV-1 TAT peptide binding to mimetic membranes. Proc. Natl. Acad. Sci. USA 2014, 111, 12684–12688. [Google Scholar] [CrossRef] [Green Version]
  178. Dorairaj, S.; Allen, T.W. On the thermodynamic stability of a charged arginine side chain in a transmembrane helix. Proc. Natl. Acad. Sci. USA 2007, 104, 4943–4948. [Google Scholar] [CrossRef] [Green Version]
  179. MacCallum, J.L.; Bennett, W.F.D.; Tieleman, D.P. Partitioning of amino acid side chains into lipid bilayers: Results from computer simulations and comparison to experiment. J. Gen. Physiol. 2007, 129, 371–377. [Google Scholar] [CrossRef]
  180. Grasso, G.; Muscat, S.; Rebella, M.; Morbiducci, U.; Audenino, A.; Danani, A.; Deriu, M.A. Cell penetrating peptide modulation of membrane biomechanics by Molecular dynamics. J. Biomech. 2018, 73, 137–144. [Google Scholar] [CrossRef]
Figure 1. Coarse-grained simulations show that the hyperpolarization favors the formation of ~2-nm-wide water pores, which would be used by CPPs to directly translocate across membranes within a few tens of nanoseconds. Membrane hyperpolarization was achieved by setting an ion imbalance through a net charge difference of 30 positive ions (corresponding to a V m of ~2 V) between the intracellular and extracellular spaces. The CPPs used in the simulations were TAT-RasGAP317–326, TAT, R9 and Penetratin, respectively, and the simulation system also contained a natural cell membrane-like composition (for both inner and outer leaflets). Image reprinted from Ref. [25].
Figure 1. Coarse-grained simulations show that the hyperpolarization favors the formation of ~2-nm-wide water pores, which would be used by CPPs to directly translocate across membranes within a few tens of nanoseconds. Membrane hyperpolarization was achieved by setting an ion imbalance through a net charge difference of 30 positive ions (corresponding to a V m of ~2 V) between the intracellular and extracellular spaces. The CPPs used in the simulations were TAT-RasGAP317–326, TAT, R9 and Penetratin, respectively, and the simulation system also contained a natural cell membrane-like composition (for both inner and outer leaflets). Image reprinted from Ref. [25].
Cells 11 04016 g001
Figure 2. Distribution of CPPs by net average charge and category. Most of the CPPs have net positive charges. Many cationic CPPs are highly charged peptides, without any amphipathic arrangement or hydrophobic character. Anionic CPPs do not form a separate class but can be classified as hydrophobic or amphipathic CPPs. Image reprinted from Ref. [45].
Figure 2. Distribution of CPPs by net average charge and category. Most of the CPPs have net positive charges. Many cationic CPPs are highly charged peptides, without any amphipathic arrangement or hydrophobic character. Anionic CPPs do not form a separate class but can be classified as hydrophobic or amphipathic CPPs. Image reprinted from Ref. [45].
Cells 11 04016 g002
Figure 3. Mechanisms of CPP uptake across the cellular membranes. The cellular uptake mechanisms are mainly divided into two categories: energy-dependent endocytosis and energy-independent direct translocation. Image reprinted from Ref. [12].
Figure 3. Mechanisms of CPP uptake across the cellular membranes. The cellular uptake mechanisms are mainly divided into two categories: energy-dependent endocytosis and energy-independent direct translocation. Image reprinted from Ref. [12].
Cells 11 04016 g003
Figure 4. Bead representation of coarse-grained (a) polarizable water; (b) DPPC lipid; (c) cyclic Arg9; and (d) linear Arg9 molecules. All of the water beads are red. W, WP, and WM are neutral, positively, and negatively charged beads of polarizable coarse-grained water. The tails and head groups of coarse-grained lipid molecules are yellow and green. The NC3, PO4, GL1, and GL2 on lipid molecule are marked for choline, phosphate, and two carbonyl beads. The backbone and non-backbone beads of coarse-grained Arg9 peptide are blue and purple. The marked strings such as BB, SC1, SC2, and SCP correspond to the backbone, two non-charged beads, and the positively charged beads, respectively. Image reprinted from Ref. [67].
Figure 4. Bead representation of coarse-grained (a) polarizable water; (b) DPPC lipid; (c) cyclic Arg9; and (d) linear Arg9 molecules. All of the water beads are red. W, WP, and WM are neutral, positively, and negatively charged beads of polarizable coarse-grained water. The tails and head groups of coarse-grained lipid molecules are yellow and green. The NC3, PO4, GL1, and GL2 on lipid molecule are marked for choline, phosphate, and two carbonyl beads. The backbone and non-backbone beads of coarse-grained Arg9 peptide are blue and purple. The marked strings such as BB, SC1, SC2, and SCP correspond to the backbone, two non-charged beads, and the positively charged beads, respectively. Image reprinted from Ref. [67].
Cells 11 04016 g004
Figure 5. Key stages during the pulling of a single penetratin molecule through a DPPC bilayer: (a) the initial system; (b) the peptide binds to the upper monolayer; (c) the water-filled pore is formed; (d) the peptide reaches the center of the bilayer; (e) the peptide reaches the surface of the lower monolayer; and (f) the peptide is pulled out of the membrane. The arrow indicates the direction of the pulling force. The semitransparent yellow surface corresponds to the lipid-water interface. Image reprinted from Ref. [73].
Figure 5. Key stages during the pulling of a single penetratin molecule through a DPPC bilayer: (a) the initial system; (b) the peptide binds to the upper monolayer; (c) the water-filled pore is formed; (d) the peptide reaches the center of the bilayer; (e) the peptide reaches the surface of the lower monolayer; and (f) the peptide is pulled out of the membrane. The arrow indicates the direction of the pulling force. The semitransparent yellow surface corresponds to the lipid-water interface. Image reprinted from Ref. [73].
Cells 11 04016 g005
Figure 6. Different mechanisms of CPPs interacting with the membrane. Image reprinted from Ref. [46].
Figure 6. Different mechanisms of CPPs interacting with the membrane. Image reprinted from Ref. [46].
Cells 11 04016 g006
Figure 7. Stages of the deformation of a DPPC bilayer with 512 lipids and 8 bound penetratin peptides. Image reprinted from Ref. [73].
Figure 7. Stages of the deformation of a DPPC bilayer with 512 lipids and 8 bound penetratin peptides. Image reprinted from Ref. [73].
Cells 11 04016 g007
Figure 8. Representative snapshots of (a) the unmodified AuNP (N = 0) and (bd) the Tat–AuNP complexes (N = 1, 4 and 8) interacting with the lipid membrane (N denotes the peptide number). Image reprinted from Ref. [157].
Figure 8. Representative snapshots of (a) the unmodified AuNP (N = 0) and (bd) the Tat–AuNP complexes (N = 1, 4 and 8) interacting with the lipid membrane (N denotes the peptide number). Image reprinted from Ref. [157].
Cells 11 04016 g008
Figure 9. Molecular structure of NP1 which contains a hydrophobic stearic acid, histidines and positively charged arginines. Image reprinted from Ref. [163].
Figure 9. Molecular structure of NP1 which contains a hydrophobic stearic acid, histidines and positively charged arginines. Image reprinted from Ref. [163].
Cells 11 04016 g009
Table 1. Overview of some representative CPPs and their properties.
Table 1. Overview of some representative CPPs and their properties.
PeptideSequenceTypeNumber of Arginine ResiduesChargeConformation *Concentration
Dependence **
Ref.
HIV-1 TatGRKKRRQRRRPPQCationic6+8random coil+[26]
Arg9RRRRRRRRRCationic9+9random coil+[40]
PenetratinRQIKIWFQNRRMKWKKCationic and
Amphipathic
3+7α-helix or β-strand+[30]
TP10AGYLLGKINLKALAALAKKILAmphipathic0+4α-helix+[35]
pVECLLIILRRRIRKQAHAHSKAmphipathic4+8β-strand or random coil+[41]
Pep-1KETWWETWWTEWSQPKKKRKVAmphipathic0+3α-helix+[43]
CADYGLWRALWRLLRSLWRLLWRAAmphipathic5+5α-helix+[38]
MAPKLALKLALKALKAALKLAAmphipathic0+5α-helix+[39]
SAP(E)VELPPPVELPPPVELPPPAmphipathic0-3α-helix-[133]
K-FGFAAVLLPVLLAAPHydrophobic00random coil+[44]
* Conformation: conformations adopted when CPPs interact with cell membranes. ** Concentration dependence: whether the uptake mode of CPPs is concentration-dependent, “+” indicates concentration-dependent and “-” indicates concentration-independent.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ouyang, J.; Sheng, Y.; Wang, W. Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations. Cells 2022, 11, 4016. https://doi.org/10.3390/cells11244016

AMA Style

Ouyang J, Sheng Y, Wang W. Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations. Cells. 2022; 11(24):4016. https://doi.org/10.3390/cells11244016

Chicago/Turabian Style

Ouyang, Jun, Yuebiao Sheng, and Wei Wang. 2022. "Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations" Cells 11, no. 24: 4016. https://doi.org/10.3390/cells11244016

APA Style

Ouyang, J., Sheng, Y., & Wang, W. (2022). Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations. Cells, 11(24), 4016. https://doi.org/10.3390/cells11244016

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop