SignificanceThe passive translocation mechanism of arginine-rich cell-penetrating peptides has puzzled the scientific community for more than 20 y. In this study we propose a hitherto unrecognized mechanism of passive cell entry involving fusion of multilamellar structures generated by the cell-penetrating peptides. The geometry of entry for this mechanism is completely different from previously suggested direct translocation mechanisms, leading to another paradigm for designing molecular carriers for drug delivery to the cell.

It is a textbook knowledge that charges of the same polarity repel each other. For two monovalent ions in the gas phase at a close contact this repulsive interaction amounts to hundreds of kilojoules per mole. In aqueous solutions, however, this Coulomb repulsion is strongly attenuated by a factor equal to the dielectric constant of the medium. The residual repulsion, which now amounts only to units of kilojoules per mole, may be in principle offset by attractive interactions. Probably the smallest cationic pair, where a combination of dispersion and cavitation forces overwhelms the Coulomb repulsion, consists of two guanidinium ions in water. Indeed, by a combination of molecular dynamics with electronic structure calculations and electrophoretic, as well as spectroscopic, experiments, we have demonstrated that aqueous guanidinium cations form (weakly) thermodynamically stable like-charge ion pairs. The importance of pairing of guanidinium cations in aqueous solutions goes beyond a mere physical curiosity, since it has significant biochemical implications. Guanidinium chloride is known to be an efficient and flexible protein denaturant. This is due to the ability of the orientationally amphiphilic guanidinium cations to disrupt various secondary structural motifs of proteins by pairing promiscuously with both hydrophobic and hydrophilic groups, including guanidinium-containing side chains of arginines. The fact that the cationic guanidinium moiety forms the dominant part of the arginine side chain implies that the like-charge ion pairing may also play a role for interactions between peptides and proteins. Indeed, arginine-arginine pairing has been frequently found in structural protein databases. In particular, when strengthened by a presence of negatively charged glutamate, aspartate, or C-terminal carboxylic groups, this binding motif helps to stabilize peptide or protein dimers and is also found in or near active sites of several enzymes. The like-charge pairing of the guanidinium side-chain groups may also hold the key to the understanding of the arginine "magic", that is, the extraordinary ability of arginine-rich polypeptides to passively penetrate across cellular membranes. Unlike polylysines, which are also highly cationic but lack the ease in crossing membranes, polyarginines do not exhibit mutual repulsion. Instead, they accumulate at the membrane, weaken it, and might eventually cross in a concerted, "train-like" manner. This behavior of arginine-rich cell penetrating peptides can be exploited when devising smart strategies how to deliver in a targeted way molecular cargos into the cell.

The surface tension of the air-water interface increases upon addition of inorganic salts, implying a negative surface excess of ionic species. Most acids, however, induce a decrease in surface tension, indicating a positive surface excess of hydrated protons. In combination with the apparent negative charge at pure air-water interfaces derived from electrokinetic experiments, this experimental observation has been a source of intense debate since the mid-19th century. Herein, we calculate surface tensions and ionic surface propensities at air-water interfaces from classical, thermodynamically consistent molecular dynamics simulations. The surface tensions of NaOH, HCl, and NaCl solutions show outstanding quantitative agreement with experiment. Of the studied ions, only H O adsorbs to the air-water interface. The adsorption is explained by the deep potential well caused by the orientation of the H O dipole in the interfacial electric field, which is confirmed by ab initio simulations.

Lipid nanodiscs are small synthetic lipid bilayer structures that are stabilized in solution by special circumscribing (or scaffolding) proteins or polymers. Because they create native-like environments for transmembrane proteins, lipid nanodiscs have become a powerful tool for structural determination of this class of systems when combined with cryo-electron microscopy or nuclear magnetic resonance. The elastic properties of lipid bilayers determine how the lipid environment responds to membrane protein perturbations, and how the lipid in turn modifies the conformational state of the embedded protein. However, despite the abundant use of nanodiscs in determining membrane protein structure, the elastic material properties of even pure lipid nanodiscs (i.e., without embedded proteins) have not yet been quantitatively investigated. A major hurdle is due to the inherently nonlocal treatment of the elastic properties of lipid systems implemented by most existing methods, both experimental and computational. In addition, these methods are best suited for very large “infinite” size lipidic assemblies, or ones that contain periodicity, in the case of simulations. We have previously described a computational analysis of molecular dynamics simulations designed to overcome these limitations, so it allows quantification of the bending rigidity ( K C ) and tilt modulus (κ t ) on a local scale even for finite, nonperiodic systems, such as lipid nanodiscs. Here we use this computational approach to extract values of K C and κ t for a set of lipid nanodisc systems that vary in size and lipid composition. We find that the material properties of lipid nanodiscs are different from those of infinite bilayers of corresponding lipid composition, highlighting the effect of nanodisc confinement. Nanodiscs tend to show higher stiffness than their corresponding macroscopic bilayers, and moreover, their material properties vary spatially within them. For small-size MSP1 nanodiscs, the stiffness decreases radially, from a value that is larger in their center than the moduli of the corresponding bilayers by a factor of ∼2–3. The larger nanodiscs (MSP1E3D1 and MSP2N2) show milder spatial changes of moduli that are composition dependent and can be maximal in the center or at some distance from it. These trends in moduli correlate with spatially varying structural properties, including the area per lipid and the nanodisc thickness. Finally, as has previously been reported, nanodiscs tend to show deformations from perfectly flat circular geometries to varying degrees, depending on size and lipid composition. The modulations of lipid elastic properties that we find should be carefully considered when making structural and functional inferences concerning embedded proteins.

We study a series of intermolecular hydrogen-bonded 1 : 1 complexes formed by chloroacetic acid with 19 substituted pyridines and one aliphatic amine dissolved in CDCl at low temperature by H andC NMR and FTIR spectroscopy. The hydrogen bond geometries in these complexes vary from molecular (O-HN) to zwitterionic (OH-N) ones, while NMR spectra show the formation of short strong hydrogen bonds in intermediate cases. Analysis of C[double bond, length as m-dash]O stretching and asymmetric CO stretching bands in FTIR spectra reveal the presence of proton tautomerism. On the basis of these data, we construct the overall proton transfer pathway. In addition to that, we also study by use of ab initio molecular dynamics the complex formed by chloroacetic acid with 2-methylpyridine, surrounded by 71 CDCl molecules, revealing a dual-maximum distribution of hydrogen bond geometries in solution. The analysis of the calculated trajectory shows that the proton jumps between molecular and zwitterionic forms are indeed driven by dipole-dipole solvent-solute interactions, but the primary cause of the jumps is the formation/breaking of weak CHO bonds from solvent molecules to oxygen atoms of the carboxylate group.

No abstract

Ab initio free energy calculations of guanidinium pairing in aqueous solution confirm the counterintuitive conjecture that the like-charge ion pair is thermodynamically stable. Transferring the guanidinium pair to the inside of a POPC lipid bilayer, like-charge ion pairing is found to occur also inside the membrane defect. It is found to contribute to the nonadditivity of ion transfer, thereby facilitating the presence of ions inside the bilayer. The effect is quantified by free energy decomposition and comparison with ammonium ions, which do not form a stable pair. The presence of two charges inside the center of the bilayer leads to the formation of a pore. Potential consequences for cell penetrating peptides and ion conduction are drawn.

A protocol for the ab initio construction of a realistic cylindrical pore in amorphous silica, serving as a geometric nanoscale confinement for liquids and solutions, is presented. Upon filling the pore with liquid water at different densities, the structure and dynamics of the liquid inside the confinement can be characterized. At high density, the pore introduces long-range oscillations into the water density profile, which makes the water structure unlike that of the bulk across the entire pore. The tetrahedral structure of water is also affected up to the second solvation shell of the pore wall. Furthermore, the effects of the confinement on hydrogen bonding and diffusion, resulting in a weakening and distortion of the water structure at the pore walls and a slowdown in diffusion, are characterized.