Link to publicationCitation for published version (APA): Sun, D., Forsman, J., Lund, M., & Woodward, C. E. (2014). Effect of arginine-rich cell penetrating peptides on membrane pore formation and life-times: a molecular simulation study. Physical Chemistry Chemical Physics, 16(38), 20785-20795. https://doi.org/10.1039/c4cp02211d General rights Copyright and moral rights for the publications made accessible in the public portal are retained by the authors and/or other copyright owners and it is a condition of accessing publications that users recognise and abide by the legal requirements associated with these rights.• Users may download and print one copy of any publication from the public portal for the purpose of private study or research.• You may not further distribute the material or use it for any profit-making activity or commercial gain • You may freely distribute the URL identifying the publication in the public portal Take down policy If you believe that this document breaches copyright please contact us providing details, and we will remove access to the work immediately and investigate your claim. 2 AbstractThe molecular basis for the effectiveness of arginine-rich cell penetrating peptides (ARCPPs) traversing cell membrane barrier is not well established. The fact that a threshold concentration of ARCPPs is required for efficient translocation in model membranes suggests cooperative action by ARCPPs. We used umbrella sampling simulations to calculate the free energies for membrane pore formation. The membrane-bound octaarginine (ARG8) peptides showed little cooperativity in lowering the free energy barrier to generate membrane pores by direct peptide translocation or by lipid flip-flop. Instead, high concentrations of ARG8 peptides were found to expand the surface area of the lipid bilayer due to the deep partitioning of guanidinium ions into the lipid glycerol regions. Surface-bound ARG8 peptides can also insert the arginine side chain into one existing transient membrane pore, and the lifetime of the transient membrane pore is significantly extended by arginine. This suggests a cooperative kinetic mechanism may act above a threshold adsorption concentration to facilitate the rapid uptake of these peptides.3

Voltage-gated sodium channels are essential for carrying electrical signals throughout the body, and mutations in these proteins are responsible for a variety of disorders, including epilepsy and pain syndromes. As such, they are the target of a number of drugs used for reducing pain or combatting arrhythmias and seizures. However, these drugs affect all sodium channel subtypes found in the body. Designing compounds to target select sodium channel subtypes will provide a new therapeutic pathway and would maximize treatment efficacy while minimizing side effects. Here, we examine the binding preferences of nine compounds known to be sodium channel pore blockers in molecular dynamics simulations. We use the approach of replica exchange solute tempering (REST) to gain a more complete understanding of the inhibitors' behavior inside the pore of NavMs, a bacterial sodium channel, and NavPas, a eukaryotic sodium channel. Using these simulations, we are able to show that both charged and neutral compounds partition into the bilayer, but neutral forms more readily cross it. We show that there are two possible binding sites for the compounds: () a site on helix 6, which has been previously determined by many experimental and computational studies, and () an additional site, occupied by protonated compounds in which the positively charged part of the drug is attracted into the selectivity filter. Distinguishing distinct binding poses for neutral and charged compounds is essential for understanding the nature of pore block and will aid the design of subtype-selective sodium channel inhibitors.

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.