Translocation thermodynamics of linear and cyclic nonaarginine into model DPPC bilayer via coarse-grained molecular dynamics simulation: implications of pore formation and nonadditivity

HP35 Protein in the Mesopore of MIL-101(Cr) MOF: A Model to Study Cotranslocational Unfolding

The immobilization of enzymes in metal-organic framework (MOF) cages is important in biotechnology. In this context, the mechanism of translocation of proteins through the cavities of the MOF and the roles played by confinement and MOF chemistry in giving rise to stable protein intermediates that are otherwise transiently populated in the physiological environment are important questions to be addressed. These unexplored aspects are examined with villin headpiece (HP35) as a model protein confined within a mesopore of MIL-101(Cr) using molecular dynamics simulations. At equilibrium, the protein is located farther from the center of the cavity and closer to the MOF surface. Molecular interactions with the MOF partially unfold helix-1 at its N-terminus. Umbrella sampling simulations inform the range of conformations that HP35 undertakes during translocation from one cavity to another and associated changes in free energy. Relative to its equilibrium state within the cavity, the free energy barrier for the unfolded protein at the cage window is estimated to be 16 kcal/mol. This study of MOF-based protein conformation can also be a general approach to observing intermediates in folding-unfolding pathways.

Oral Delivery of the Vancomycin Derivative FU002 by a Surface-Modified Liposomal Nanocarrier

Oral delivery of peptide therapeutics faces multiple challenges due to their instability in the gastrointestinal tract and low permeation capability. In this study, the aim is to develop a liposomal nanocarrier formulation to enable the oral delivery of the vancomycin-peptide derivative FU002. FU002 is a promising, resistance-breaking, antibiotic which exhibits poor oral bioavailability, limiting its potential therapeutic use. To increase its oral bioavailability, FU002 is incorporated into tetraether lipid-stabilized liposomes modified with cyclic cell-penetrating peptides on the liposomal surface. This liposomal formulation shows strong binding to Caco-2 cells without exerting cytotoxic effects in vitro. Pharmacokinetics studies in vivo in rats reveal increased oral bioavailability of liposomal FU002 when compared to the free drug. In vitro and in vivo antimicrobial activity of FU002 are preserved in the liposomal formulation. As a highlight, oral administration of liposomal FU002 results in significant therapeutic efficacy in a murine systemic infection model. Thus, the presented nanotechnological approach provides a promising strategy for enabling oral delivery of this highly active vancomycin derivative.

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

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.

Cell-Penetrating Peptides

In this introductory chapter, we first define cell-penetrating peptides (CPPs), give short overview of CPP history and discuss several aspects of CPP classification. Next section is devoted to the mechanism of CPP penetration into the cells, where direct and endocytic internalization of CPP is explained. Kinetics of internalization is discussed more extensively, since this topic is not discussed in other chapters of this book. At the end of this section some features of the thermodynamics of CPP interaction with the membrane is also presented. Finally, we present different cargoes that can be transferred into the cells by CPPs and briefly discuss the effect of cargo on the rate and efficiency of penetration into the cells.

Cell-Penetrating Peptides: Correlation between Peptide-Lipid Interaction and Penetration Efficiency

Cell-penetrating peptides are used in the delivery of peptides and biologics, with some cell-penetrating peptides found to be more efficient than others. The exact mechanism of how they interact with the cell membrane and penetrate it, however, remains unclear. This study attempts to investigate the difference in free energy profiles of three cell-penetrating peptides (TAT, CPP1 and CPP9) with a model lipid bilayer (DOPC) using molecular dynamics pulling simulations with umbrella sampling. Potential mean force (PMF) and free energy barrier between the peptides and DOPC are determined using WHAM analysis and MM-PBSA analysis, respectively. CPP9 is found to have the smallest PMF value, followed by CPP1 and TAT, consistent with the experimental data. YDEGE peptide, however, does not give the highest PMF value, although it is a non-cell-permeable peptide. YDEGE is also found to form water pores, alongside with TAT and CPP9, suggesting that it is difficult to distinguish true water pore formation from artefacts arising from pulling simulations. On the contrary, free energy analysis of the peptide-DOPC complex at the lipid-water interface with MM-PBSA provides results consistent with experimental data with CPP9 having the least interaction with DOPC and lowest free energy barrier, followed by CPP1, TAT and YDEGE. These findings suggest that peptide-lipid interaction at the lipid-water interface has a direct correlation with the penetration efficiency of peptides across the lipid bilayer.

Surface charge density and fatty acids enhance the membrane permeation rate of CPP-cargo complexes

The CPP-effect makes reference to the process by which the membrane translocation rate of a cargo is enhanced by chemical functionalization with cell-penetrating peptides (CPPs). In this work we combine a simple kinetic model with free-energy calculations to explore the energetic basis of the CPP-effect. Two polyglicines are selected as model hydrophilic cargoes, and nona-arginine as a prototypical CPP. We assess the cargo carrying efficiency of nona-arginine by comparing the adsorption and insertion energies of the cargoes, the cargo-free CPPs, and the CPP-cargo complexes, into lipid membranes of varying composition. We also analyze the effect of modifying the type and concentration of anionic lipids, and the implication of these factors on the translocation rate of the CPP-cargo complex. Of particular interest is the evaluation of the catalytic role of palmitic acid (palmitate) as a promoter of the CPP-effect. We also analyse the influence of the size of the cargo on the membrane adsorption and insertion energies. Our results show that the efficiency of nona-arginine as a transmembrane carrier of simple hydrophilic molecules is modulated by the size of the cargo, and is strongly enhanced by increasing the concentration of anionic lipids and of ionized fatty acids in the membrane.

Bioelectricity for Drug Delivery: The Promise of Cationic Therapeutics

Biological systems overwhelmingly comprise charged entities generating electrical activity that can have significant impact on biological structure and function. This intrinsic bio-electrical activity can also be harnessed for overcoming the tissue matrix and cell membrane barriers, which have been outstanding challenges for targeted drug delivery, by using rationally designed cationic carriers. The weak and reversible long-range electrostatic interactions with fixed negatively charged groups facilitate electro-diffusive transport of cationic therapeutics through full-tissue thickness to effectively reach intra-tissue, cellular, and intracellular target sites. This article presents a perspective on the promise of using rationally designed cationic biomaterials in targeted drug delivery, the underlying charge-based mechanisms, and bio-transport phenomena while addressing outstanding concerns around toxicity and methods to mitigate them. We also discuss electrically charged drugs that are currently being evaluated in clinical trials and identify areas of further development that have the potential to usher in new treatments.

Thermodynamics of cell penetrating peptides on lipid membranes: sequence and membrane acidity regulate surface binding

Acidic lipids respond to pH in ways that fully promote or deplete the surface accumulation of cell penetrating peptides.

Interaction of aliphatic amino acids with zwitterionic and charged lipid membranes: hydration and dehydration phenomena

Aliphatic amino acids interact differently in order to induce gelation or fluidization in zwitterionic and charged lipid membranes as a result of hydration or dehydration of the membrane surface.

All-Factor Analysis and Correlations on the Transmembrane Process for Arginine-Rich Cell-Penetrating Peptides

Currently, arginine-rich cell-penetrating peptides (CPPs), due to their little cytotoxicity and high transmembrane efficiency, are considered as one of the important intracellular carriers. Although the mechanism of the transmembrane process for arginine-rich CPPs was proposed, the quantitative correlations and the key factors involved in this process still deserve further investigation. In this study, all-atom molecular dynamics and the umbrella sampling technique were employed to study the arginine-rich CPPs transmembrane process. In the adsorption process of CPPs from solution to the surface of the lipid bilayer, the adsorption free energy (ΔG_A) is found to be linearly related to the interaction energy change (ΔE_A): ΔG_A = 0.0426ΔE_A + 36.7, R² = 0.92. In the CPPs transmembrane process, the transmembrane free energy barrier (ΔG_B) is roughly correlated with the corresponding interaction energy change (ΔE_B): ΔG_B = 0.108ΔE_B +135, R² = 0.73. The multiple salt bridges of guanidinium-PO₄ account for 65% of the overall interaction energy, so the increased negative charges of the lipid bilayer or more salt bridges would facilitate CPPs adsorption and transmembrane processes. Also, the increased negative charges of the lipid bilayer would reduce the amount of water to be carried into the pore and further reduce the ΔG_B. The peptide backbone would not have a direct impact on transmembrane efficiency. The ΔG_B is also found to be related to the length of the pore (L): ΔG_B = 46.2L - 31.3, R² = 0.92, which makes the transmembrane efficiency estimable. This work is expected to deliver an in-depth understanding and help the optimization of CPPs.

Membrane potential drives direct translocation of cell-penetrating peptides

Cell-penetrating peptides (CPPs) are frequently employed as drug delivery agents with rapid cellular uptake, however, the uptake mechanism and the detailed translocation pathway are at present not completely understood. Both endocytosis and direct translocation through membrane pores have been observed in experiments and simulations under different conditions. Here we report the molecular dynamics simulations providing evidence for the direct translocation of CPPs across the membrane driven by the membrane electrostatic potential. The local membrane potential can be produced by the ion concentration imbalance across the membrane, which is ubiquitous in biological environments. Moreover, if positively charged CPPs are adsorbed on the membrane, this further enhances the membrane potential, opening membrane pores through which CPPs can be instantly transported in a chain-like configuration. The classical nucleation theory is applied to estimate the translocation time by calculating the changes in the free energy upon transferring CPPs across the membrane at different potentials, showing good agreement with available experimental measurements. The revealed CPP translocation mechanism can be broadly relevant for cellular processes in biology.

The interfacial electrostatic potential modulates the insertion of cell-penetrating peptides into lipid bilayers

Cell-penetrating peptides (CPP) are short sequences of cationic amino-acids that show a surprising ability to traverse lipid bilayers. CPP are considered to be some of the most effective vectors to introduce membrane-impermeable cargos into cells, but the molecular basis of the membrane translocation mechanisms and its dependence on relevant membrane physicochemical properties have yet to be fully determined. In this paper we resort to Molecular Dynamics simulations and experiments to investigate how the electrostatic potential across the lipid/water interface affects the insertion of hydrophilic and amphipathic CPP into two-dimensional lipid structures. Simulations are used to quantify the effect of the transmembrane potential on the free-energy profile associated with the transfer of the CPP across a neutral lipid bilayer. It is found that the electrostatic bias has a relatively small effect on the binding of the peptides to the membrane surface, but that it significantly lowers the permeation barrier. A charge compensation mechanism, arising from the segregation of counter-ions while the peptide traverses the membrane, determines the shape and symmetry of the free-energy curves and underlines relevant mechanistic considerations. Langmuir monolayer experiments performed with a variety of amphiphiles model the incorporation of the CPP into the external membrane leaflet. It is shown that the dipole potential of the monolayer controls the extent of penetration of the CPP into the lipid aggregate, to a greater degree than its surface charge.

Structural effects and translocation of spontaneous membrane-translocating peptides with POPC bilayer

Martini coarse-grained force field simulations have been carried out to estimate the free energy profiles of the spontaneous membrane-translocating peptide TP2 and one negative control peptide ONEG with POPC as the model bilayer. The results show that the free energy minimum of TP2 is [Formula: see text]20[Formula: see text]kJ/mol lower than that of ONEG. In addition, the minimum of TP2 shifts slightly to the bilayer center compared with ONEG. The translocation barrier height for TP2 and ONEG are 119.0[Formula: see text]kJ/mol and 155.7[Formula: see text]kJ/mol, respectively. The lower central energy barrier of TP2 facilitates the transition between two leaflets of POPC. Both translocating peptides induce the formation of funnel-shaped structures at the bilayer center, but TP2 has a more compact structure and brings less perturbation compared with ONEG. Subsequently all atom molecular simulations testify the findings. It is indicated that compared with its negative control ONEG, TP2 binds better with lipid and penetrates deeper into bilayer with less perturbation to the bilayer structure. Our findings may shed light on the design and virtual screening of spontaneous membrane-translocating peptides.

Passive Transmembrane Permeation Mechanisms of Monovalent Ions Explored by Molecular Dynamics Simulations

Passive or unassisted ion permeation through lipid bilayers involves a type of rare events by which cells regulate their salt concentrations and pH. It is important to understand its mechanism in order to develop technologies of, for example, delivering or maintaining small drug-like molecules inside cells. In earlier simulations of passive ion permeations, the commonly used sampling methods usually define the positions of ions relative to the membrane as a measure of permeation, i.e., the collective variable, ignoring the active participations of other particles. Newly defined collective variables involving the movements of ions, lipids, and water molecules allow us to identify the transition paths on the free energy landscape using the 2D umbrella sampling techniques. In this work, this technique was used to study the permeation processes of some well-known ions, sodium, potassium, and chloride. It is found permeations of sodium and potassium are assisted by important lipid bilayer deformations and massive water solvation, while chloride may not. Chloride may have two different possible pathways, in which the energetic favorable one is similar to the solubility-diffusion model. The free energy barriers for the permeation of these ions are in semiquantitative agreement with experiments. Further analyses on the distributions of oxygens and interaction energies suggest the electrostatic interactions between ions and polar headgroups of lipids may greatly influence membrane deformation as well as the water wire and furthermore the free energy barriers of waterwire mediated pathways. For chloride, the nonwaterwire pathway may be energetically favorable.

Thermodynamics of cell-penetrating HIV1 TAT peptide insertion into PC/PS/CHOL model bilayers through transmembrane pores: the roles of cholesterol and anionic lipids

Efficient delivery of pharmaceutically active molecules across cellular membranes using cell penetrating peptides (CPPs), such as the cationic human immunodeficiency virus-1 trans-acting activator of transcription peptide (HIV-1 TAT), continues to attract scientific attention in drug design and disease treatment. Experimental results show that the TAT peptide is not only capable of directly penetrating the biological membrane in a passive manner, but also forming physical, membrane-spanning pores that may facilitate transport. Experiments further show that anionic lipids accelerate peptide permeation within a range of mole percentage composition. In this work, we explored the structures and translocation thermodynamics of the cationic TAT peptide across a series of DPPC/DPPS model membranes with the presence of 0-30 mol% cholesterol. We computed the potentials of the mean force by using umbrella sampling molecular dynamics simulations coupled to the Martini coarse-grained force field. We systematically investigated the roles of cholesterol and anionic lipids (membrane surface charge) in TAT peptide translocation. In qualitative agreement with experimental findings, the barrier heights were significantly reduced in the presence of anionic lipids. A toroidal hydrophilic pore was strongly suggested by membrane structure analysis. Cholesterol stabilizes the liquid-ordered (Lo) phase of membranes and increases the elastic stiffness of bilayers. Consequently, it hinders transmembrane pore formation and thus modulates solute permeability, since the liquid-ordered phase suppresses reorientation of the lipid molecules on simulation time scales. Though cholesterol contributes marginally to the total free energy associated with peptide permeation, the coordination of cholesterol to the peptide weakens more favorable peptide-lipid interactions. The addition of the anionic lipid DPPS to the neutral DPPC bilayer leads to the emergence and further enhancement of an interfacially stable state of the peptide due to the favorable peptide-anionic lipid interactions. Translocation free energy barriers decrease in lockstep with increasing DPPS composition in the model bilayers simulated. Finally, we investigated the size of hydrophilic pores emerging in our simulations, as well as the qualitative mobility of the peptide on the membrane surface.

Prediction of Cell-Penetrating Peptides

The in silico methods for the prediction of the cell-penetrating peptides are reviewed. Those include the multivariate statistical methods, machine-learning methods such as the artificial neural networks and support vector machines, and molecular modeling techniques including molecular docking and molecular dynamics.The applicability of the methods is demonstrated on the basis of the exemplary cases from the literature.

Interaction of peptides with cell membranes: insights from molecular modeling

The investigation of the interaction of peptides with cell membranes is the focus of active research. It can enhance the understanding of basic membrane functions such as membrane transport, fusion, and signaling processes, and it may shed light on potential applications of peptides in biomedicine. In this review, we will present current advances in computational studies on the interaction of different types of peptides with the cell membrane. Depending on the properties of the peptide, membrane, and external environment, the peptide-membrane interaction shows a variety of different forms. Here, on the basis of recent computational progress, we will discuss how different peptides could initiate membrane pores, translocate across the membrane, induce membrane endocytosis, produce membrane curvature, form fibrils on the membrane surface, as well as interact with functional membrane proteins. Finally, we will present a conclusion summarizing recent progress and providing some specific insights into future developments in this field.

Investigating Hydrophilic Pores in Model Lipid Bilayers Using Molecular Simulations: Correlating Bilayer Properties with Pore-Formation Thermodynamics

Cell-penetrating and antimicrobial peptides show a remarkable ability to translocate across physiological membranes. Along with factors such as electric-potential-induced perturbations of membrane structure and surface tension effects, experiments invoke porelike membrane configurations during the solute transfer process into vesicles and cells. The initiation and formation of pores are associated with a nontrivial free-energy cost, thus necessitating a consideration of the factors associated with pore formation and the attendant free energies. Because of experimental and modeling challenges related to the long time scales of the translocation process, we use umbrella sampling molecular dynamics simulations with a lipid-density-based order parameter to investigate membrane-pore-formation free energy employing Martini coarse-grained models. We investigate structure and thermodynamic features of the pore in 18 lipids spanning a range of headgroups, charge states, acyl chain lengths, and saturation. We probe the dependence of pore-formation barriers on the area per lipid, lipid bilayer thickness, and membrane bending rigidities in three different lipid classes. The pore-formation free energy in pure bilayers and peptide translocating scenarios are significantly coupled with bilayer thickness. Thicker bilayers require more reversible work to create pores. The pore-formation free energy is higher in peptide-lipid systems than in peptide-free lipid systems due to penalties to maintain the solvation of charged hydrophilic solutes within the membrane environment.

Reconciling structural and thermodynamic predictions using all-atom and coarse-grain force fields: the case of charged oligo-arginine translocation into DMPC bilayers

Using the translocation of short, charged cationic oligo-arginine peptides (mono-, di-, and triarginine) from bulk aqueous solution into model DMPC bilayers, we explore the question of the similarity of thermodynamic and structural predictions obtained from molecular dynamics simulations using all-atom and Martini coarse-grain force fields. Specifically, we estimate potentials of mean force associated with translocation using standard all-atom (CHARMM36 lipid) and polarizable and nonpolarizable Martini force fields, as well as a series of modified Martini-based parameter sets. We find that we are able to reproduce qualitative features of potentials of mean force of single amino acid side chain analogues into model bilayers. In particular, modifications of peptide-water and peptide-membrane interactions allow prediction of free energy minima at the bilayer-water interface as obtained with all-atom force fields. In the case of oligo-arginine peptides, the modified parameter sets predict interfacial free energy minima as well as free energy barriers in almost quantitative agreement with all-atom force field based simulations. Interfacial free energy minima predicted by a modified coarse-grained parameter set are -2.51, -4.28, and -5.42 for mono-, di-, and triarginine; corresponding values from all-atom simulations are -0.83, -3.33, and -3.29, respectively, all in units of kcal/mol. We found that a stronger interaction between oligo-arginine and the membrane components and a weaker interaction between oligo-arginine and water are crucial for producing such minima in PMFs using the polarizable CG model. The difference between bulk aqueous and bilayer center states predicted by the modified coarse-grain force field are 11.71, 14.14, and 16.53 kcal/mol, and those by the all-atom model are 6.94, 8.64, and 12.80 kcal/mol; those are of almost the same order of magnitude. Our simulations also demonstrate a remarkable similarity in the structural aspects of the ensemble of configurations generated using the all-atom and coarse-grain force fields. Both resolutions show that oligo-arginine peptides adopt preferential orientations as they translocate into the bilayer. The guiding theme centers on charged groups maintaining coordination with polar and charged bilayer components as well as local water. We also observe similar behaviors related with membrane deformations.

Structural and Thermodynamic Insight into Spontaneous Membrane-Translocating Peptides Across Model PC/PG Lipid Bilayers

We present results of Martini coarse-grained force field simulations to estimate the potentials of mean force for a series of recently screened spontaneous membrane-translocating peptides, SMTPs. We consider model bilayer composed of POPC and POPG, the latter providing the anionic component as used in experimental studies. We observe a significant barrier for translocation in the case of the canonical cationic cell-penetrating peptide nona-arginine, ARG9. In the case of the TP1, TP2, and TP3 peptides, potentials of mean force are systematically lower relative to the ARG9 case. Though the barriers predicted by the simulations, on the order of 20 kcal/mol, are still rather large to recapitulate the experimental kinetics of internalization, we emphasize that the qualitative trend of reduction of barrier heights is a significant result. Decomposition of the PMFs indicates that though there is a substantial entropic stability when the peptides reside at bilayer center, barriers as predicted from these force field-based studies are largely determined by enthalpic (potential energy) interactions. We note that the binding of the SMTPs is critically dependent on the mix of hydrophilic and hydrophobic residues that constitute the amino acid motif/sequence of these peptides. For the cationic ARG9 which only contains hydrophilic residues, there is no tight binding observed. The specific motif [Formula: see text] (where [Formula: see text] is a general residue) is a potential sequence in drug/peptide design. The SMTPs with this motif are able to translocate into membrane at a significantly lower free energy cost, compared to the negative control peptides. Finally, we compare the different membrane perturbations induced by the presence of the different peptides in the bilayer center. In some cases, hydrophilic pores are observed to form, thus conferring stability to the internalized state. In other cases, SMTPs are associated only with membrane defects such as induced membrane curvature. These latter observations suggest some influence of membrane rigidity as embodied in the full range of membrane undulatory modes in defining pore-forming propensities in bilayers.