Low Free Energy Barrier for Ion Permeation Through Double-Helical Gramicidin

Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.

A Polarizable Atomic Multipole-Based Force Field for Molecular Dynamics Simulations of Anionic Lipids

In all of the classical force fields, electrostatic interaction is simply treated and explicit electronic polarizability is neglected. The condensed-phase polarization, relative to the gas-phase charge distributions, is commonly accounted for in an average way by increasing the atomic charges, which remain fixed throughout simulations. Based on the lipid polarizable force field DMPC and following the same framework as Atomic Multipole Optimized Energetics for BiomoleculAr (AMOEBA) simulation, the present effort expands the force field to new anionic lipid models, in which the new lipids contain DMPG and POPS. The parameters are compatible with the AMOEBA force field, which includes water, ions, proteins, etc. The charge distribution of each atom is represented by the permanent atomic monopole, dipole and quadrupole moments, which are derived from the ab initio gas phase calculations. Many-body polarization including the inter- and intramolecular polarization is modeled in a consistent manner with distributed atomic polarizabilities. Molecular dynamics simulations of the two aqueous DMPG and POPS membrane bilayer systems, consisting of 72 lipids with water molecules, were then carried out to validate the force field parameters. Membrane width, area per lipid, volume per lipid, deuterium order parameters, electron density profile, electrostatic potential difference between the center of the bilayer and water are all calculated, and compared with limited experimental data.

Effective lipid-detergent system for study of membrane active peptides in fluid liposomes

The structure of peptide antibiotic gramicidin A (gA) was studied in phosphatidylcholin liposomes modified by nonionic detergent Triton X-100. First, the detergent : lipid ratio at which the saturation of lipid membrane by Triton X-100 occurs (Re (sat)), was determined by light scattering. Measurements of steady-state fluorescence anisotropy of 1,6-diphenyl-1,3,5-hexatriene at sublytic concentrations of detergent showed that after saturation of the membrane by Triton X-100 microviscosity of lipid bilayer is reduced by 20%. The equilibrium conformational state of gA in phosphatidylcholine liposomes at Re (sat) was studied by CD spectroscopy. It was found that the conformational state of this channel-forming peptide changed crucially when Triton X-100 induced transition to more fluid membranes. The gA single-channel measurements were made with Triton X-100 containing bilayers. Tentative assignment of the channel type and gA structures was made by correlation of CD data with conductance histograms. Lipid-detergent system with variable viscosity developed in this work can be used to study the structure and folding of other membrane-active peptides.

Small-Anion Selective Transmembrane "Holes" Induced by an Antimicrobial Peptide Too Short to Span Membranes

Whereas many membrane-destabilization modes have been suggested for membrane-spanning antimicrobial peptides (AMPs), few are available for those too short to span membrane thickness. Here we show that ORB-1, a 15-residue disulfide-bridged AMP that is only ∼20 Å long even when fully stretched like a hairpin, may act by inducing small anion-selective transmembrane "holes" of negative mean curvature. In model membranes of Gram-negative bacteria, ORB-1 induces chloride transmembrane transport and formation of transmembrane channels of negative mean curvature, whereas the inactive analogue, ORB-N, does not, suggesting a correlation between antibacterial activity and ability to induce transmembrane channels. Given that ORB-N is the C-terminus amidated form of ORB-1, our results further suggest that formation of membrane-spanning dimers may be required to initiate the observed channel induction. Moreover, ORB-1 renders model bacterial membranes permeable to anions with effective hydration diameters of <1 nm (e.g., Cl(-) and NO3(-)), but not cations of similar sizes (e.g., H3O(+)), indicative of anion-selective transmembrane channels with an effective inner diameter of ≤1 nm. In addition, negative-intrinsic-curvature (NIC) lipids such as phosphoethanolamine (PE) may facilitate the membrane-destabilization process of ORB-1. Our findings may expand current understandings on how AMPs destabilize membranes and facilitate the pharmaceutical development of ORB-1.

High-Throughput Simulations of Dimer and Trimer Assembly of Membrane Proteins. The DAFT Approach

Interactions between membrane proteins are of great biological significance and are consequently an important target for pharmacological intervention. Unfortunately, it is still difficult to obtain detailed views on such interactions, both experimentally, where the environment hampers atomic resolution investigation, and computationally, where the time and length scales are problematic. Coarse grain simulations have alleviated the later issue, but the slow movement through the bilayer, coupled to the long life times of nonoptimal dimers, still stands in the way of characterizing binding distributions. In this work, we present DAFT, a Docking Assay For Transmembrane components, developed to identify preferred binding orientations. The method builds on a program developed recently for generating custom membranes, called insane (INSert membrANE). The key feature of DAFT is the setup of starting structures, for which optimal periodic boundary conditions are devised. The purpose of DAFT is to perform a large number of simulations with different components, starting from unbiased noninteracting initial states, such that the simulations evolve collectively, in a manner reflecting the underlying energy landscape of interaction. The implementation and characteristic features of DAFT are explained, and the efficacy and relaxation properties of the method are explored for oligomerization of glycophorin A dimers, polyleucine dimers and trimers, MS1 trimers, and rhodopsin dimers. The results suggest that, for simple helices, such as GpA and polyleucine, in POPC/DOPC membranes series of 500 simulations of 500 ns each allow characterization of the helix dimer orientations and allow comparing associating and nonassociating components. However, the results also demonstrate that short simulations may suffer significantly from nonconvergence of the ensemble and that using too few simulations may obscure or distort features of the interaction distribution. For trimers, simulation times exceeding several microseconds appear needed, due to the increased complexity. Similarly, characterization of larger proteins, such as rhodopsin, takes longer time scales due to the slower diffusion and the increased complexity of binding interfaces. DAFT and its auxiliary programs have been made available from http://cgmartini.nl/ , together with a working example.