Molecular dynamics simulations suggest a mechanism for translocation of the HIV-1 TAT peptide across lipid membranes

Ion Channels Made from a Single Membrane-Spanning DNA Duplex

Because of their hollow interior, transmembrane channels are capable of opening up pathways for ions across lipid membranes of living cells. Here, we demonstrate ion conduction induced by a single DNA duplex that lacks a hollow central channel. Decorated with six porpyrin-tags, our duplex is designed to span lipid membranes. Combining electrophysiology measurements with all-atom molecular dynamics simulations, we elucidate the microscopic conductance pathway. Ions flow at the DNA-lipid interface as the lipid head groups tilt toward the amphiphilic duplex forming a toroidal pore filled with water and ions. Ionic current traces produced by the DNA-lipid channel show well-defined insertion steps, closures, and gating similar to those observed for traditional protein channels or synthetic pores. Ionic conductances obtained through simulations and experiments are in excellent quantitative agreement. The conductance mechanism realized here with the smallest possible DNA-based ion channel offers a route to design a new class of synthetic ion channels with maximum simplicity.

Molecular Mechanism of Ionic-Liquid-Induced Membrane Disruption: Morphological Changes to Bilayers, Multilayers, and Vesicles

The application of ionic liquids (ILs) in many industrially relevant processes provides an urgent need to better understand their molecular interactions with biological systems. A detailed understanding of the cytotoxicity mechanism of ILs can be helpful in facilitating the molecular design of nontoxic ILs. Using coarse-grained molecular dynamics (MD) simulations, we investigate the effects of imidazolium-based ILs on several lipid bilayer morphologies. Our results demonstrate that the asymmetric insertion of IL cations into one side of a lipid bilayer leaflet enhances the leaflet strain, which upon reaching a critical value triggers a morphological disruption in the bilayer. Consistently, the bending modulus of the bilayer is reduced by 1 to 2 orders of magnitude relative to that of an IL-free planar bilayer prior to the disruption event. Our results suggest that ILs that can easily insert into the lipid bilayer without diffusing across or inducing lipid flip-flop can be more disruptive to a lipid biomembrane.

Cell-penetrating peptides: Possible transduction mechanisms and therapeutic applications

Cell-penetrating peptides (CPPs), also known as protein transduction domains, are a class of diverse peptides with 5-30 amino acids. CPPs are divided into cationic, amphipathic and hydrophobic CPPs. They are able to carry small molecules, plasmid DNA, small interfering RNA, proteins, viruses, imaging agents and other various nanoparticles across the cellular membrane, resulting in internalization of the intact cargos. However, the mechanisms of CPP internalization remain to be elucidated. Recently, CPPs have received considerable attention due to their high transduction efficiency and low cytotoxicity. These peptides have a significant potential for diagnostic and therapeutic applications, such as delivery of fluorescent or radioactive compounds for imaging, delivery of peptides and proteins for therapeutic application, and delivery of molecules into induced pluripotent stem cells for directing differentiation. The present study reviews the classifications and transduction mechanisms of CPPs, as well as their potential applications.

Ultrafast Killing and Self-Gelling Antimicrobial Imidazolium Oligomers

Infectious diseases and the increasing threat of worldwide pandemics have underscored the importance of antibiotics and hygiene. Intensive efforts have been devoted to developing new antibiotics to meet the rapidly growing demand. In particular, advancing the knowledge of the structure-property-activity relationship is critical to expedite the design and development of novel antimicrobial with the needed potential and efficacy. Herein, a series of new antimicrobial imidazolium oligomers are developed with the rational manipulation of terminal group's hydrophobicity. These materials exhibit superior activity, excellent selectivity, ultrafast killing (>99.7% killing within 30 s), and desirable self-gelling properties. Molecular dynamic simulations reveal the delicate effect of structural changes on the translocation motion across the microbial cell membrane. The energy barrier of the translocation process analyzed by free energy calculations provides clear kinetic information to suggest that the spontaneous penetration requires a very short timescale of seconds to minutes for the new imidazolium oligomers.

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.

Maximally asymmetric transbilayer distribution of anionic lipids alters the structure and interaction with lipids of an amyloidogenic protein dimer bound to the membrane surface

We used molecular dynamics simulations to explore the effects of asymmetric transbilayer distribution of anionic phosphatidylserine (PS) lipids on the structure of a protein on the membrane surface and subsequent protein-lipid interactions. Our simulation systems consisted of an amyloidogenic, beta-sheet rich dimeric protein (D42) absorbed to the phosphatidylcholine (PC) leaflet, or protein-contact PC leaflet, of two membrane systems: a single-component PC bilayer and double PC/PS bilayers. The latter comprised of a stable but asymmetric transbilayer distribution of PS in the presence of counterions, with a 1-component PC leaflet coupled to a 1-component PS leaflet in each bilayer. The maximally asymmetric PC/PS bilayer had a non-zero transmembrane potential (TMP) difference and higher lipid order packing, whereas the symmetric PC bilayer had a zero TMP difference and lower lipid order packing under physiologically relevant conditions. Analysis of the adsorbed protein structures revealed weaker protein binding, more folding in the N-terminal domain, more aggregation of the N- and C-terminal domains and larger tilt angle of D42 on the PC leaflet surface of the PC/PS bilayer versus the PC bilayer. Also, analysis of protein-induced membrane structural disruption revealed more localized bilayer thinning in the PC/PS versus PC bilayer. Although the electric field profile in the non-protein-contact PS leaflet of the PC/PS bilayer differed significantly from that in the non-protein-contact PC leaflet of the PC bilayer, no significant difference in the electric field profile in the protein-contact PC leaflet of either bilayer was evident. We speculate that lipid packing has a larger effect on the surface adsorbed protein structure than the electric field for a maximally asymmetric PC/PS bilayer. Our results support the mechanism that the higher lipid packing in a lipid leaflet promotes stronger protein-protein but weaker protein-lipid interactions for a dimeric protein on membrane surfaces.

How transmembrane peptides insert and orientate in biomembranes: a combined experimental and simulation study

After the synthesis of transmembrane peptides/proteins (TMPs), their insertion into a lipid bilayer is a fundamental biophysical process.

Effect of lipid head group interactions on membrane properties and membrane-induced cationic β-hairpin folding

Stages in membrane induced SVS-1 folding.

Penetration of HIV-1 Tat47-57 into PC/PE Bilayers Assessed by MD Simulation and X-ray Scattering

The interactions of the basic, cell-penetrating region (Y47GRKKRRQRRR57) of the HIV-1 Tat protein with dioleoylphosphatidylcholine (DOPC) bilayers were previously assessed by comparing experimental X-ray diffuse scattering with atomistic molecular dynamics simulations. Here, we extend this investigation by evaluating the influence of phosphatidylethanolamine (PE) lipids. Using experimental bilayer form factors derivedfrom X-ray diffuse scattering data as a guide, our simulations indicate that Tat peptides localize close to the carbonyl-glycerol group in the headgroup region of bilayers composed of either DOPC or DOPC:DOPE (1:1) lipid. Our results also suggest that Tat peptides may more frequently insert into the hydrophobic core of bilayers composed of PC:PE (1:1) lipids than into bilayers composed entirely of PC lipids. PE lipids may facilitate peptide translocation across a lipid bilayer by stabilizing intermediate states in which hydrated peptides span the bilayer.

Ligand-Modified Human Serum Albumin Nanoparticles for Enhanced Gene Delivery

The development of nonviral gene delivery systems is a great challenge to enable safe gene therapy. In this study, ligand-modified nanoparticles based on human serum albumin (HSA) were developed and optimized for an efficient gene therapy. Different glutaraldehyde cross-linking degrees were investigated to optimize the HSA nanoparticles for gene delivery. The peptide sequence arginine-glycine-aspartate (RGD) and the HIV-1 transactivator of transduction sequence (Tat) are well-known as promising targeting ligands. Plasmid DNA loaded HSA nanoparticles were covalently modified on their surface with these different ligands. The transfection potential of the obtained plasmid DNA loaded RGD- and Tat-modified nanoparticles was investigated in vitro, and optimal incubation conditions for these preparations were studied. It turned out that Tat-modified HSA nanoparticles with the lowest cross-linking degree of 20% showed the highest transfection potential. Taken together, ligand-functionalized HSA nanoparticles represent promising tools for efficient and safe gene therapy.

Cell-penetrating peptides transport therapeutics into cells

Nearly 30years ago, certain small, relatively nontoxic peptides were discovered to be capable of traversing the cell membrane. These cell-penetrating peptides, as they are now called, have been shown to not only be capable of crossing the cell membrane themselves but can also carry many different therapeutic agents into cells, including small molecules, plasmid DNA, siRNA, therapeutic proteins, viruses, imaging agents, and other various nanoparticles. Many cell-penetrating peptides have been derived from natural proteins, but several other cell-penetrating peptides have been developed that are either chimeric or completely synthetic. How cell-penetrating peptides are internalized into cells has been a topic of debate, with some peptides seemingly entering cells through an endocytic mechanism and others by directly penetrating the cell membrane. Although the entry mechanism is still not entirely understood, it seems to be dependent on the peptide type, the peptide concentration, the cargo the peptide transports, and the cell type tested. With new intracellular disease targets being discovered, cell-penetrating peptides offer an exciting approach for delivering drugs to these intracellular targets. There are hundreds of cell-penetrating peptides being studied for drug delivery, and ongoing studies are demonstrating their success both in vitro and in vivo.

Membrane permeation of a peptide: it is better to be positive

A joint experimental and computational study investigates the translocation of a tryptophan molecule through a phospholipid membrane. Time dependent spectroscopy of the tryptophan side chain determines the rate of permeation into 150 nm phospholipid vesicles. Atomically detailed simulations are conducted to calculate the free energy profiles and the permeation coefficient. Different charging conditions of the peptide (positive, negative, or zwitterion) are considered. Both experiment and simulation reproduce the qualitative trend and suggest that the fastest permeation is when the tryptophan is positively charged. The permeation mechanism, which is revealed by molecular dynamics simulations, is of a translocation assisted by a local defect. The influence of long-range electrostatic interactions, such as the membrane dipole potential on the permeation process, is not significant.

Modeling kinetics and equilibrium of membranes with fields: milestoning analysis and implication to permeation

Coarse graining of membrane simulations by translating atomistic dynamics to densities and fields with Milestoning is discussed. The space of the membrane system is divided into cells and the different cells are characterized by order parameters presenting the number densities. The dynamics of the order parameters are probed with Milestoning. The methodology is illustrated here for a phospholipid membrane system (a hydrated bilayer of DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine) lipid molecules). Significant inhomogeneity in membrane internal number density leads to complex free energy landscape and local maps of transition times. Dynamics and distributions of cavities within the membrane assist the permeation of nonpolar solutes such as xenon atoms. It is illustrated that quantitative and detailed dynamics of water transport through DOPC membrane can be analyzed using Milestoning with fields. The reaction space for water transport includes at least two slow variables: the normal to the membrane plane, and the water density.

NKCS, a Mutant of the NK-2 Peptide, Causes Severe Distortions and Perforations in Bacterial, But Not Human Model Lipid Membranes

NKCS is an improved mutant of the bioactive peptide NK-2, which shows strong activity against Escherichia coli and low toxicity towards human cells. The different activity demonstrates the relevance of the physico-chemical nature of the target membrane for the biological effect of this peptide. We studied the effect of this potent antimicrobial peptide on model membranes by activity studies, differential scanning calorimetry, single molecule tracking and tracer efflux experiments. We found that NKCS severely distorted, penetrated and perforated model lipid membranes that resembled bacterial membranes, but not those that were similar to human cell membranes. The interactions of NKCS with phosphatidylethanolamine, which is abundant in bacterial membranes, were especially strong and are probably responsible for its antimicrobial activity.

Physical chemistry of nanomedicine: understanding the complex behaviors of nanoparticles in vivo

Nanomedicine is an interdisciplinary field of research at the interface of science, engineering, and medicine, with broad clinical applications ranging from molecular imaging to medical diagnostics, targeted therapy, and image-guided surgery. Despite major advances during the past 20 years, there are still major fundamental and technical barriers that need to be understood and overcome. In particular, the complex behaviors of nanoparticles under physiological conditions are poorly understood, and detailed kinetic and thermodynamic principles are still not available to guide the rational design and development of nanoparticle agents. Here we discuss the interactions of nanoparticles with proteins, cells, tissues, and organs from a quantitative physical chemistry point of view. We also discuss insights and strategies on how to minimize nonspecific protein binding, how to design multistage and activatable nanostructures for improved drug delivery, and how to use the enhanced permeability and retention effect to deliver imaging agents for image-guided cancer surgery.

Evaluating Force Fields for the Computational Prediction of Ionized Arginine and Lysine Side-Chains Partitioning into Lipid Bilayers and Octanol

Abundant peptides and proteins containing arginine (Arg) and lysine (Lys) amino acids can apparently permeate cell membranes with ease. However, the mechanisms by which these peptides and proteins succeed in traversing the free energy barrier imposed by cell membranes remain largely unestablished. Precise thermodynamic studies (both theoretical and experimental) on the interactions of Arg and Lys residues with model lipid bilayers can provide valuable clues to the efficacy of these cationic peptides and proteins. We have carried out molecular dynamics simulations to calculate the interactions of ionized Arg and Lys side-chains with the zwitterionic 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayer for 10 widely used lipid/protein force fields: CHARMM36/CHARMM36, SLIPID/AMBER99SB-ILDN, OPLS-AA/OPLS-AA, Berger/OPLS-AA, Berger/GROMOS87, Berger/GROMOS53A6, GROMOS53A6/GROMOS53A6, nonpolarizable MARTINI, polarizable MARTINI, and BMW MARTINI. We performed umbrella sampling simulations to obtain the potential of mean force for Arg and Lys side-chains partitioning from water to the bilayer interior. We found significant differences between the force fields, both for the interactions between side-chains and bilayer surface, as well as the free energy cost for placing the side-chain at the center of the bilayer. These simulation results were compared with the Wimley-White interfacial scale. We also calculated the free energy cost for transferring ionized Arg and Lys side-chains from water to both dry and wet octanol. Our simulations reveal rapid diffusion of water molecules into octanol whereby the equilibrium mole fraction of water in the wet octanol phase was ∼25%. Surprisingly, our free energy calculations found that the high water content in wet octanol lowered the water-to-octanol partitioning free energies for cationic residues by only 0.6 to 0.7 kcal/mol.

Molecular modeling of peptides

This article presents a review of the field of molecular modeling of peptides. The main focus is on atomistic modeling with molecular mechanics potentials. The description of peptide conformations and solvation through potentials is discussed. Several important computer simulation methods are briefly introduced, including molecular dynamics, accelerated sampling approaches such as replica-exchange and metadynamics, free energy simulations and kinetic network models like Milestoning. Examples of recent applications for predictions of structure, kinetics, and interactions of peptides with complex environments are described. The reliability of current simulation methods is analyzed by comparison of computational predictions obtained using different models with each other and with experimental data. A brief discussion of coarse-grained modeling and future directions is also presented.

Cell penetrating peptide TAT can kill cancer cells via membrane disruption after attachment of camptothecin

Attachment of traditional anticancer drugs to cell penetrating peptides is an effective strategy to improve their application in cancer treatment. In this study, we designed and synthesized the conjugates TAT-CPT and TAT-2CPT by attaching camptothecin (CPT) to the N-terminus of the cell penetrating peptide TAT. Interestingly, we found that TAT-CPT and especially TAT-2CPT could kill cancer cells via membrane disruption, which is similar to antimicrobial peptides. This might be because that CPT could perform as a hydrophobic residue to increase the extent of membrane insertion of TAT and the stability of the pores. In addition, TAT-CPT and TAT-2CPT could also kill cancer cells by the released CPT after they entered cells. Taken together, attachment of CPT could turn cell penetrating peptide TAT into an antimicrobial peptide with a dual mechanism of anticancer action, which presents a new strategy to develop anticancer peptides based on cell penetrating peptides.

Chirality switching within an anionic cell-penetrating peptide inhibits translocation without affecting preferential entry

Multiple substitution of d- for l-amino acids decreases the intracellular uptake of cationic cell penetrating peptides (CPP) in a cell line-dependent manner. We show here that a single d-amino acid substitution can decrease the overall uptake of the anionic, amphipathic CPP, p28, into cancer and histologically matched normal cell lines, while not altering the preferential uptake of p28 into cancer cells. The decrease appears dependent on the position of the d-substitution within the peptide and the ability of the substituted d-amino acid to alter chirality. We also suggest that when d-substitution alters the ratio of α-helix to β-sheet content of an anionic CPP, its translocation across the cell membrane is altered, reducing overall entry. These observations may have a significant effect on the design of future d-substituted analogues of cell penetrating peptides.

Fundamental molecular mechanism for the cellular uptake of guanidinium-rich molecules

Guanidinium-rich molecules, such as cell-penetrating peptides, efficiently enter living cells in a non-endocytic energy-independent manner and transport a wide range of cargos, including drugs and biomarkers. The mechanism by which these highly cationic molecules efficiently cross the hydrophobic barrier imposed by the plasma membrane remains a fundamental open question. Here, a combination of computational results and in vitro and live-cell experimental evidence reveals an efficient energy-independent translocation mechanism for arginine-rich molecules. This mechanism unveils the essential role of guanidinium groups and two universal cell components: fatty acids and the cell membrane pH gradient. Deprotonated fatty acids in contact with the cell exterior interact with guanidinium groups, leading to a transient membrane channel that facilitates the transport of arginine-rich peptides toward the cell interior. On the cytosolic side, the fatty acids become protonated, releasing the peptides and resealing the channel. This fundamental mechanism appears to be universal across cells from different species and kingdoms.

Microsecond molecular dynamics simulations of lipid mixing

Molecular dynamics (MD) simulations of membranes are often hindered by the slow lateral diffusion of lipids and the limited time scale of MD. In order to study the dynamics of mixing and characterize the lateral distribution of lipids in converged mixtures, we report microsecond-long all-atom MD simulations performed on the special-purpose machine Anton. Two types of mixed bilayers, POPE:POPG (3:1) and POPC:cholesterol (2:1), as well as a pure POPC bilayer, were each simulated for up to 2 μs. These simulations show that POPE:POPG and POPC:cholesterol are each fully miscible at the simulated conditions, with the final states of the mixed bilayers similar to a random mixture. By simulating three POPE:POPG bilayers at different NaCl concentrations (0, 0.15, and 1 M), we also examined the effect of salt concentration on lipid mixing. While an increase in NaCl concentration is shown to affect the area per lipid, tail order, and lipid lateral diffusion, the final states of mixing remain unaltered, which is explained by the largely uniform increase in Na(+) ions around POPE and POPG. Direct measurement of water permeation reveals that the POPE:POPG bilayer with 1 M NaCl has reduced water permeability compared with those at zero or low salt concentration. Our calculations provide a benchmark to estimate the convergence time scale of all-atom MD simulations of lipid mixing. Additionally, equilibrated structures of POPE:POPG and POPC:cholesterol, which are frequently used to mimic bacterial and mammalian membranes, respectively, can be used as starting points of simulations involving these membranes.

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.

A C-terminal membrane anchor affects the interactions of prion proteins with lipid membranes

Membrane attachment via a C-terminal glycosylphosphatidylinositol anchor is critical for conversion of PrP(C) into pathogenic PrP(Sc). Therefore the effects of the anchor on PrP structure and function need to be deciphered. Three PrP variants, including full-length PrP (residues 23-231, FL_PrP), N-terminally truncated PrP (residues 90-231, T_PrP), and PrP missing its central hydrophobic region (Δ105-125, ΔCR_PrP), were equipped with a C-terminal membrane anchor via a semisynthesis strategy. Analyses of the interactions of lipidated PrPs with phospholipid membranes demonstrated that C-terminal membrane attachment induces a different binding mode of PrP to membranes, distinct from that of non-lipidated PrPs, and influences the biochemical and conformational properties of PrPs. Additionally, fluorescence-based assays indicated pore formation by lipidated ΔCR_PrP, a variant that is known to be highly neurotoxic in transgenic mice. This finding was supported by using patch clamp electrophysiological measurements of cultured cells. These results provide new evidence for the role of the membrane anchor in PrP-lipid interactions, highlighting the importance of the N-terminal and the central hydrophobic domain in these interactions.

Atomistic simulations of pore formation and closure in lipid bilayers

Cellular membranes separate distinct aqueous compartments, but can be breached by transient hydrophilic pores. A large energetic cost prevents pore formation, which is largely dependent on the composition and structure of the lipid bilayer. The softness of bilayers and the disordered structure of pores make their characterization difficult. We use molecular-dynamics simulations with atomistic detail to study the thermodynamics, kinetics, and mechanism of pore formation and closure in DLPC, DMPC, and DPPC bilayers, with pore formation free energies of 17, 45, and 78 kJ/mol, respectively. By using atomistic computer simulations, we are able to determine not only the free energy for pore formation, but also the enthalpy and entropy, which yields what is believed to be significant new insights in the molecular driving forces behind membrane defects. The free energy cost for pore formation is due to a large unfavorable entropic contribution and a favorable change in enthalpy. Changes in hydrogen bonding patterns occur, with increased lipid-water interactions, and fewer water-water hydrogen bonds, but the total number of overall hydrogen bonds is constant. Equilibrium pore formation is directly observed in the thin DLPC lipid bilayer. Multiple long timescale simulations of pore closure are used to predict pore lifetimes. Our results are important for biological applications, including the activity of antimicrobial peptides and a better understanding of membrane protein folding, and improve our understanding of the fundamental physicochemical nature of membranes.

HIV-1 Tat membrane interactions probed using X-ray and neutron scattering, CD spectroscopy and MD simulations

We report the effect on lipid bilayers of the Tat peptide Y47GRKKRRQRRR57 from the HIV-1 virus transactivator of translation (Tat) protein. Synergistic use of low-angle X-ray scattering (LAXS) and atomistic molecular dynamic simulations (MD) indicate Tat peptide binding to neutral dioleoylphosphocholine (DOPC) lipid headgroups. This binding induced the local lipid phosphate groups to move 3Å closer to the center of the bilayer. Many of the positively charged guanidinium components of the arginines were as close to the center of the bilayer as the locally thinned lipid phosphate groups. LAXS data for DOPC, DOPC/dioleoylphosphoethanolamine (DOPE), DOPC/dioleoylphosphoserine (DOPS), and a mimic of the nuclear membrane gave similar results. Generally, the Tat peptide decreased the bilayer bending modulus KC and increased the area/lipid. Further indications that Tat softens a membrane, thereby facilitating translocation, were provided by wide-angle X-ray scattering (WAXS) and neutron scattering. CD spectroscopy was also applied to further characterize Tat/membrane interactions. Although a mechanism for translation remains obscure, this study suggests that the peptide/lipid interaction makes the Tat peptide poised to translocate from the headgroup region.

The importance of membrane defects-lessons from simulations

The defects and pores within lipid membranes are scientifically interesting and have a number of biological applications. Although lipid bilayers are extremely thin hydrophobic barriers, just ∼3 nm thick, they include diverse chemistry and have complex structures. Bilayers are soft and dynamic, and as a result, they can bend and deform in response to different stimuli by means of structural changes in their component lipids. Though defects occur within these structures, their transience and small size have made it difficult to characterize them. However, with recent advances in computer power and computational modeling techniques, researchers can now use simulations as a powerful tool to probe the mechanism and energies of defect and pore formation in a number of situations. In this Account, we present results from our detailed molecular dynamics computer simulations of hydrophilic pores and related defects in lipid bilayers at an atomistic level. Electroporation can be used to increase the permeability of cellular membranes, with potential therapeutic applications. Atomistic simulations of electroporation have illustrated the molecular details of this process, including the importance of water dipole interactions at the water-membrane interface. Characterization of the lipid-protein interactions provides an important tool for understanding transmembrane protein structure and thermodynamic stability. Atomistic simulations give a detailed picture of the free energies of model peptides and side chains in lipid membranes; the energetic cost of defect formation strongly influences the energies of interactions between lipids and polar and charged residues. Many antimicrobial peptides form hydrophilic pores in lipid membranes, killing bacteria or cancer cells. On the basis of simulation data, at least some of these peptides form defects and pores near the center of the bilayer, with a common disordered structure where hydrated headgroups form an approximately toroidal shape. The localization and trafficking of lipids supports general membrane structure and a number of important signaling cascades, such as those involving ceramide, diacylglycerol, and cholesterol. Atomistic simulations have determined the rates and free energies of lipid flip-flop. During the flip-flop of most phosphatidylcholine lipids, a hydrophilic pore forms when the headgroup moves near the center of the bilayer. Simulations have provided novel insight into many features of defects and pores in lipid membranes. Simulation data from very different systems and models show how water penetration and defect formation can determine the free energies of many membrane processes. Bilayers can deform and allow transient defects and pores when exposed to a diverse range of stimuli. Future work will explore many aspects of membrane defects with increased resolution and scope, including the study of more complex lipid mixtures, membrane domains, and large-scale membrane remodeling. Such studies will examine processes including vesicle budding and fusion, non-bilayer lipid phases, and interactions between lipid bilayers and other biomolecules. Simulations provide information that complements experimental studies, allowing microscopic insight into experimental observations and suggesting novel hypotheses and experiments. These studies should enable a deeper understanding of the role of lipid bilayers in cellular biology and support the development of future lipid-based biotechnology.

Theoretical insight into the relationship between the structures of antimicrobial peptides and their actions on bacterial membranes

Antimicrobial peptides with diverse cationic charges, amphiphathicities, and secondary structures possess a variety of antimicrobial activities against bacteria, fungi, and other generalized targets. To illustrate the relationship between the structures of these peptide and their actions at microscopic level, we present systematic coarse-grained dissipative particle dynamics simulations of eight types of antimicrobial peptides with different secondary structures interacting with a lipid bilayer membrane. We find that the peptides use multiple mechanisms to exert their membrane-disruptive activities: A cationic charge is essential for the peptides to selectively target negatively charged bacterial membranes. This cationic charge is also responsible for promoting electroporation. A significant hydrophobic portion is necessary to disrupt the membrane through formation of a permeable pore or translocation. Alternatively, the secondary structure and the corresponding rigidity of the peptides determine the pore structure and the translocation pathway.

Peptide-ligand binding modeling of siRNA with cell-penetrating peptides

The binding affinity of a series of cell-penetrating peptides (CPP) was modeled through docking and making use of the number of intermolecular hydrogen bonds, lipophilic contacts, and the number of sp3 molecular orbital hybridization carbons. The new ranking of the peptides is consistent with the experimentally determined efficiency in the downregulation of luciferase activity, which includes the peptides' ability to bind and deliver the siRNA into the cell. The predicted structures of the complexes of peptides to siRNA were stable throughout 10 ns long, explicit water molecular dynamics simulations. The stability and binding affinity of peptide-siRNA complexes was related to the sidechains and modifications of the CPPs, with the stearyl and quinoline groups improving affinity and stability. The reranking of the peptides docked to siRNA, together with explicit water molecular dynamics simulations, appears to be well suited to describe and predict the interaction of CPPs with siRNA.

Early endosomal escape of a cyclic cell-penetrating peptide allows effective cytosolic cargo delivery

Cyclic heptapeptide cyclo(FΦRRRRQ) (cFΦR4, where Φ is l-2-naphthylalanine) was recently found to be efficiently internalized by mammalian cells. In this study, its mechanism of internalization was investigated by perturbing various endocytic events through the introduction of pharmacologic agents and genetic mutations. The results show that cFΦR4 binds directly to membrane phospholipids, is internalized into human cancer cells through endocytosis, and escapes from early endosomes into the cytoplasm. Its cargo capacity was examined with a wide variety of molecules, including small-molecule dyes, linear and cyclic peptides of various charged states, and proteins. Depending on the nature of the cargos, they may be delivered by endocyclic (insertion of cargo into the cFΦR4 ring), exocyclic (attachment of cargo to the Gln side chain), or bicyclic approaches (fusion of cFΦR4 and cyclic cargo rings). The overall delivery efficiency (i.e., delivery of cargo into the cytoplasm and nucleus) of cFΦR4 was 4-12-fold higher than those of nonaarginine, HIV Tat-derived peptide, or penetratin. The higher delivery efficiency, coupled with superior serum stability, minimal toxicity, and synthetic accessibility, renders cFΦR4 a useful transporter for intracellular cargo delivery and a suitable system for investigating the mechanism of endosomal escape.

Antimicrobial potency and selectivity of simplified symmetric-end peptides

Because antimicrobial peptides (AMPs) are potentially useful for the treatment of multidrug-resistant infections, more attention is being paid to the structural modification and structure-function relationship of both naturally occurring and synthetic AMPs. Previous studies indicated that Protegrin-1 (PG-1), isolated from porcine leukocytes, exhibited considerable antimicrobial activity and cytotoxicity. The β-turn of PG-1 floated on the surface of bacterial membrane, while its β-strand inserted into the bacterial membrane and formed pores that were dedicated to producing cytotoxicity. For reducing cytotoxicity and improving cells selectivity, we designed a series of simplified symmetric-end peptides by combining the β-turn of PG-1 with simple amino acid repeat sequences. The sequence of designed symmetric-end peptides is (XR)nH(RX)n, (n = 1,2; X represents I, F, W and P; H represents CRRRFC). The symmetric-end peptides displayed antimicrobial activity against both gram-positive and gram-negative bacteria. In particular, (XR)2H(RX)2 (X here is I, F and W) showed greater antimicrobial potency than PG-1. Hemolysis activity and cytotoxicity, detected by using human red blood cells (RBCs) and human embryonic lung fibroblasts MRC-5 cells, were observably lower than the native peptide PG-1. (IR)2H(RI)2 (IR2), folded into β-sheet structures, displayed the highest therapeutic index, suggesting its great cell selectivity. The fluorescence spectroscopy, flow cytometry, and electron microscopy observation indicated that IR2 exhibited great membrane penetration potential by inducing membrane blebbing, disruption and lysis. Collectively, generating symmetric-end β-sheet peptides is a promising strategy for designing effective AMPs with great antimicrobial activities and cell selectivity.

Water-mediated ion-ion interactions are enhanced at the water vapor-liquid interface

There is overwhelming evidence that ions are present near the vapor-liquid interface of aqueous salt solutions. Charged groups can also be driven to interfaces by attaching them to hydrophobic moieties. Despite their importance in many self-assembly phenomena, how ion-ion interactions are affected by interfaces is not understood. We use molecular simulations to show that the effective forces between small ions change character dramatically near the water vapor-liquid interface. Specifically, the water-mediated attraction between oppositely charged ions is enhanced relative to that in bulk water. Further, the repulsion between like-charged ions is weaker than that expected from a continuum dielectric description and can even become attractive as the ions are drawn to the vapor side. We show that thermodynamics of ion association are governed by a delicate balance of ion hydration, interfacial tension, and restriction of capillary fluctuations at the interface, leading to nonintuitive phenomena, such as water-mediated like charge attraction. "Sticky" electrostatic interactions may have important consequences on biomolecular structure, assembly, and aggregation at soft liquid interfaces. We demonstrate this by studying an interfacially active model peptide that changes its structure from α-helical to a hairpin-turn-like one in response to charging of its ends.

Molecular Simulation of the Transport of Drugs across Model Membranes

This Perspective describes recent progress in the area of the molecular simulation of the interactions of hydrophobic and hydrophilic solutes with membranes. The ability to predict drug solubility prior to synthesis is an extremely desirable goal for pharmaceutical research. A major advantage of molecular dynamics is the ability to computationally probe both the drug solubility as well as the pathway for the transport of drugs across membranes. Computational methods to predict the interaction free energy of solutes with membranes have advanced significantly in recent years and can characterize the intra- and intermolecular state of the drug as well as perturbations of the drug to the membrane environment itself. In addition to a brief review on computational methods to characterize the transport of drugs across membranes, we will highlight recent discoveries and discuss the challenges and opportunities in the field.