Voltage-dependent insertion of alamethicin at phospholipid/water and octane/water interfaces

Equivalence of Charge Imbalance and External Electric Fields during Free Energy Calculations of Membrane Electroporation

Electric fields across lipid membranes play important roles in physiology, medicine, and biotechnology, rationalizing the wide interest in modeling transmembrane potentials in molecular dynamics simulations. Transmembrane potentials have been implemented with external electric fields or by imposing charge imbalance between the two water compartments of a stacked double-membrane system. We compare the two methods in the context of membrane electroporation, which involves a large change of membrane structure and capacitance. We show that, given that Ewald electrostatics are defined with tinfoil boundary conditions, the two methods lead to (i) identical potentials of mean force (PMFs) of pore formation and expansion at various potentials, demonstrating that the two methods impose equivalent driving forces for large-scale transitions at membranes, and (ii) to identical polarization of water within thin water wires or open pores, suggesting that the two methods furthermore impose equivalent local electric fields. Without tinfoil boundary conditions, effects from external fields on pore formation are spuriously suppressed or even removed. Together, our study shows that both methods, external fields and charge imbalance, are well suitable for studying large-scale transitions of lipid membranes that involve changes of membrane capacitance. However, using charge imbalance is technically more challenging for maintaining a constant transmembrane potential since it requires updating of the charge imbalance as the membrane capacitance changes.

Structure Changes of a Membrane Polypeptide under an Applied Voltage Observed with Surface-Enhanced 2D IR Spectroscopy

The structures of many membrane-bound proteins and polypeptides depend on the membrane potential. However, spectroscopically studying their structures under an applied field is challenging, because a potential is difficult to generate across more than a few bilayers. We study the voltage-dependent structures of the membrane-bound polypeptide, alamethicin, using a spectroelectrochemical cell coated with a rough, gold film to create surface plasmons. The plasmons sufficiently enhance the 2D IR signal to measure a single bilayer. The film is also thick enough to conduct current and thereby apply a potential. The 2D IR spectra resolve features from both 3₁₀- and α-helical structures and cross-peaks connecting the two. We observe changes in the peak intensity, not their frequencies, upon applying a voltage. A similar change occurs with pH, which is known to alter the angle of alamethicin relative to the surface normal. The spectra are modeled using a vibrational exciton Hamiltonian, and the voltage-dependent spectra are consistent with a change in angle of the 3₁₀- and α-helices in the membrane from 55 to 44°and from 31 to 60°, respectively. The 3₁₀- and α-helices are coupled by approximately 10 cm^-1. These experiments provide new structural information about alamethicin under a potential difference and demonstrate a technique that might be applied to voltage-gated membrane proteins and compared to molecular dynamics structures.

Modification of lipid membrane compressibility induced by an electric field

Changes in membrane deformation and compressibility, induced by an external electric field, are investigated using coarse-grained martini force field simulations in a salt-free environment. We observe changes in the area of the membrane above a critical electric field. Below this value, the membrane compressibility modulus is found to decrease monotonically. For higher electric fields, the membrane projected area remains constant while the net interfacial area increases, with the corresponding compressibility moduli, show the opposite behavior. We find that the mechanical parameters, surface tension and bending modulus, of a freely floating membrane in the absence of explicit ions, are unaffected by the presence of the electric field. We believe these results have a bearing on our understanding of the electroformation of uncharged lipids in a salt-free environment.

Monitoring the Orientational Changes of Alamethicin during Incorporation into Bilayer Lipid Membranes

Antimicrobial peptides (AMPs) are the first line of defense after contact of an infectious invader, for example, bacterium or virus, with a host and an integral part of the innate immune system of humans. Their broad spectrum of biological functions ranges from cell membrane disruption over facilitation of chemotaxis to interaction with membrane-bound or intracellular receptors, thus providing novel strategies to overcome bacterial resistances. Especially, the clarification of the mechanisms and dynamics of AMP incorporation into bacterial membranes is of high interest, and different mechanistic models are still under discussion. In this work, we studied the incorporation of the peptaibol alamethicin (ALM) into tethered bilayer lipid membranes on electrodes in combination with surface-enhanced infrared absorption (SEIRA) spectroscopy. This approach allows monitoring the spontaneous and potential-induced ion channel formation of ALM in situ. The complex incorporation kinetics revealed a multistep mechanism that points to peptide-peptide interactions prior to penetrating the membrane and adopting the transmembrane configuration. On the basis of the anisotropy of the backbone amide I and II infrared absorptions determined by density functional theory calculations, we employed a mathematical model to evaluate ALM reorientations monitored by SEIRA spectroscopy. Accordingly, ALM was found to adopt inclination angles of ca. 69°-78° and 21° in its interfacially adsorbed and transmembrane incorporated states, respectively. These orientations can be stabilized efficiently by the dipolar interaction with lipid head groups or by the application of a potential gradient. The presented potential-controlled mechanistic study suggests an N-terminal integration of ALM into membranes as monomers or parallel oligomers to form ion channels composed of parallel-oriented helices, whereas antiparallel oligomers are barred from intrusion.

Structural Behavior of the Peptaibol Harzianin HK VI in a DMPC Bilayer: Insights from MD Simulations

Microsecond molecular dynamics simulations of harzianin HK VI (HZ) interacting with a dimyristoylphosphatidylcholine bilayer were performed at the condition of low peptide-to-lipid ratio. Two orientations of HZ molecule in the bilayer were found and characterized. In the orientation perpendicular to the bilayer surface, HZ induces a local thinning of the bilayer. When inserted into the bilayer parallel to its surface, HZ is located nearly completely within the hydrophobic region of the bilayer. A combination of solid-state NMR and circular dichroism experiments found the latter orientation to be dominant. An extended sampling simulation provided qualitative results and showed the same orientation to be a global minimum of free energy. The secondary structure of HZ was characterized, and it was found to be located in the 3₁₀-helical family. The specific challenges of computer simulation of nonpolar peptides are discussed briefly.

A molecular insight into the electro-transfer of small molecules through electropores driven by electric fields

The transport of chemical compounds across the plasma membrane into the cell is relevant for several biological and medical applications. One of the most efficient techniques to enhance this uptake is reversible electroporation. Nevertheless, the detailed molecular mechanism of transport of chemical species (dyes, drugs, genetic materials, …) following the application of electric pulses is not yet fully elucidated. In the past decade, molecular dynamics (MD) simulations have been conducted to model the effect of pulsed electric fields on membranes, describing several aspects of this phenomenon. Here, we first present a comprehensive review of the results obtained so far modeling the electroporation of lipid membranes, then we extend these findings to study the electrotransfer across lipid bilayers subject to microsecond pulsed electric fields of Tat11, a small hydrophilic charged peptide, and of siRNA. We use in particular a MD simulation protocol that allows to characterize the transport of charged species through stable pores. Unexpectedly, our results show that for an electroporated bilayer subject to transmembrane voltages in the order of 500mV, i.e. consistent with experimental conditions, both Tat11 and siRNA can translocate through nanoelectropores within tens of ns. We discuss these results in comparison to experiments in order to rationalize the mechanism of drug uptake by cells. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Insights into the function of ion channels by computational electrophysiology simulations

Ion channels are of universal importance for all cell types and play key roles in cellular physiology and pathology. Increased insight into their functional mechanisms is crucial to enable drug design on this important class of membrane proteins, and to enhance our understanding of some of the fundamental features of cells. This review presents the concepts behind the recently developed simulation protocol Computational Electrophysiology (CompEL), which facilitates the atomistic simulation of ion channels in action. In addition, the review provides guidelines for its application in conjunction with the molecular dynamics software package GROMACS. We first lay out the rationale for designing CompEL as a method that models the driving force for ion permeation through channels the way it is established in cells, i.e., by electrochemical ion gradients across the membrane. This is followed by an outline of its implementation and a description of key settings and parameters helpful to users wishing to set up and conduct such simulations. In recent years, key mechanistic and biophysical insights have been obtained by employing the CompEL protocol to address a wide range of questions on ion channels and permeation. We summarize these recent findings on membrane proteins, which span a spectrum from highly ion-selective, narrow channels to wide diffusion pores. Finally we discuss the future potential of CompEL in light of its limitations and strengths. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Simulations of outer membrane channels and their permeability

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Molecular dynamics simulations on PGLa using NMR orientational constraints

NMR data obtained by solid state NMR from anisotropic samples are used as orientational constraints in molecular dynamics simulations for determining the structure and dynamics of the PGLa peptide within a membrane environment. For the simulation the recently developed molecular dynamics with orientational constraints technique (MDOC) is used. This method introduces orientation dependent pseudo-forces into the COSMOS-NMR force field. Acting during a molecular dynamics simulation these forces drive molecular rotations, re-orientations and folding in such a way that the motional time-averages of the tensorial NMR properties are consistent with the experimentally measured NMR parameters. This MDOC strategy does not depend on the initial choice of atomic coordinates, and is in principle suitable for any flexible and mobile kind of molecule; and it is of course possible to account for flexible parts of peptides or their side-chains. MDOC has been applied to the antimicrobial peptide PGLa and a related dimer model. With these simulations it was possible to reproduce most NMR parameters within the experimental error bounds. The alignment, conformation and order parameters of the membrane-bound molecule and its dimer were directly derived with MDOC from the NMR data. Furthermore, this new approach yielded for the first time the distribution of segmental orientations with respect to the membrane and the order parameter tensors of the dimer systems. It was demonstrated the deuterium splittings measured at the peptide to lipid ratio of 1/50 are consistent with a membrane spanning orientation of the peptide.

Multiscale modeling of droplet interface bilayer membrane networks

Droplet interface bilayer (DIB) networks are considered for the development of stimuli-responsive membrane-based materials inspired by cellular mechanics. These DIB networks are often modeled as combinations of electrical circuit analogues, creating complex networks of capacitors and resistors that mimic the biomolecular structures. These empirical models are capable of replicating data from electrophysiology experiments, but these models do not accurately capture the underlying physical phenomena and consequently do not allow for simulations of material functionalities beyond the voltage-clamp or current-clamp conditions. The work presented here provides a more robust description of DIB network behavior through the development of a hierarchical multiscale model, recognizing that the macroscopic network properties are functions of their underlying molecular structure. The result of this research is a modeling methodology based on controlled exchanges across the interfaces of neighboring droplets. This methodology is validated against experimental data, and an extension case is provided to demonstrate possible future applications of droplet interface bilayer networks.

Structure and orientation of antibiotic peptide alamethicin in phospholipid bilayers as revealed by chemical shift oscillation analysis of solid state nuclear magnetic resonance and molecular dynamics simulation

The structure, topology and orientation of membrane-bound antibiotic alamethicin were studied using solid state nuclear magnetic resonance (NMR) spectroscopy. (13)C chemical shift interaction was observed in [1-(13)C]-labeled alamethicin. The isotropic chemical shift values indicated that alamethicin forms a helical structure in the entire region. The chemical shift anisotropy of the carbonyl carbon of isotopically labeled alamethicin was also analyzed with the assumption that alamethicin molecules rotate rapidly about the bilayer normal of the phospholipid bilayers. It is considered that the adjacent peptide planes form an angle of 100° or 120° when it forms α-helix or 310-helix, respectively. These properties lead to an oscillation of the chemical shift anisotropy with respect to the phase angle of the peptide plane. Anisotropic data were acquired for the 4 and 7 sites of the N- and C-termini, respectively. The results indicated that the helical axes for the N- and C-termini were tilted 17° and 32° to the bilayer normal, respectively. The chemical shift oscillation curves indicate that the N- and C-termini form the α-helix and 310-helix, respectively. The C-terminal 310-helix of alamethicin in the bilayer was experimentally observed and the unique bending structure of alamethicin was further confirmed by measuring the internuclear distances of [1-(13)C] and [(15)N] doubly-labeled alamethicin. Molecular dynamics simulation of alamethicin embedded into dimyristoyl phophatidylcholine (DMPC) bilayers indicates that the helical axes for α-helical N- and 310-helical C-termini are tilted 12° and 32° to the bilayer normal, respectively, which is in good agreement with the solid state NMR results.

Electric Field Effect on Phospholipid Monolayers at an Aqueous-Organic Liquid-Liquid Interface

The electric potential difference across cell membranes, known as the membrane potential, plays an important role in the activation of many biological processes. To investigate the effect of the membrane potential on the molecular ordering of lipids within a biomimetic membrane, a self-assembled monolayer of 1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine (SOPC) lipids at an electrified 1,2-dichloroethane/water interface is studied with X-ray reflectivity and interfacial tension. Measurements over a range of electric potential differences, -150 to +130 mV, that encompass the range of typical biomembrane potentials demonstrate a nearly constant and stable structure whose lipid interfacial density is comparable to that found in other biomimetic membrane systems. Measurements at higher positive potentials, up to 330 mV, illustrate a monotonic decrease in the lipid interfacial density and accompanying variations in the interfacial configuration of the lipid. Molecular dynamics simulations, designed to mimic the experimental conditions, show that the measured changes in lipid configuration are due primarily to the variation in area per lipid with increasing applied electric field. Rotation of the SOPC dipole moment by the torque from the applied electric field appears to be negligible, except at the highest measured potentials. The simulations confirm in atomistic detail the measured potential-dependent characteristics of SOPC monolayers. Our hybrid study sheds light on phospholipid monolayer stability under different membrane potentials, which is important for understanding membrane processes. This study also illustrates the use of X-ray surface scattering to probe the ordering of surfactant monolayers at an electrified aqueous-organic liquid-liquid interface.

Water channel formation and ion transport in linear and branched lipid bilayers

The lipid bilayer stability and water channel morphologies are affected by the presence of methyl branches on lipid tails.

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8Å pore and the octamer in an 11Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted "barrel-stave" model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.

Constant electric field simulations of the membrane potential illustrated with simple systems

Advances in modern computational methods and technology make it possible to carry out extensive molecular dynamics simulations of complex membrane proteins based on detailed atomic models. The ultimate goal of such detailed simulations is to produce trajectories in which the behavior of the system is as realistic as possible. A critical aspect that requires consideration in the case of biological membrane systems is the existence of a net electric potential difference across the membrane. For meaningful computations, it is important to have well validated methodologies for incorporating the latter in molecular dynamics simulations. A widely used treatment of the membrane potential in molecular dynamics consists of applying an external uniform electric field E perpendicular to the membrane. The field acts on all charged particles throughout the simulated system, and the resulting applied membrane potential V is equal to the applied electric field times the length of the periodic cell in the direction perpendicular to the membrane. A series of test simulations based on simple membrane-slab models are carried out to clarify the consequences of the applied field. These illustrative tests demonstrate that the constant-field method is a simple and valid approach for accounting for the membrane potential in molecular dynamics studies of biomolecular systems. This article is part of a Special Issue entitled: Membrane protein structure and function.

Coarse grained model for exploring voltage dependent ion channels

The relationship between the membrane voltage and the gating of voltage activated ion channels and other systems have been a problem of great current interest. Unfortunately, reliable molecular simulations of external voltage effects present a major challenge, since meaningful converging microscopic simulations are not yet available and macroscopic treatments involve major uncertainties in terms of the dielectric used and other key features. This work extends our coarse grained (CG) model to simulations of membrane/protein systems under external potential. Special attention is devoted to a consistent modeling of the effect of external potential due to the electrodes, emphasizing semimacroscopic description of the electrolytes in the solution regions between the membranes and the electrodes, as well as the coupling between the combined potential from the electrodes plus the electrolytes and the protein ionized groups. We also provide a clear connection to microscopic treatment of the electrolytes and thus can explore possible conceptual problems that are hard to resolve by other current approaches. For example, we obtain a clear description of the charge distribution in the entire electrolyte system, including near the electrodes in membrane/electrodes systems (where continuum models do not seem to provide the relevant results). Furthermore, the present treatment provides an insight on the distribution of the electrolyte charges before and after equilibration across the membrane, and thus on the nature of the gating charge. The different aspects of the model have been carefully validated by considering problems ranging for the simple Debye-Huckel, and the Gouy-Chapman models to the evaluation of the electrolyte distribution between two electrodes, as well as the effect of extending the simulation system by periodic replicas. Overall the clear connection to microscopic descriptions combined with the power of the CG modeling seems to offer a powerful tool for exploring the balance between the protein conformational energy and the interaction with the external potential in voltage activated channels. To illustrate these features we present a preliminary study of the gating charge in the voltage activated Kv1.2 channel, using the actual change in the electrolyte charge distribution rather than the conventional macroscopic estimate. We also discuss other special features of the model, which include the ability to capture the effect of changes in the protonation states of the protein residues during the close to open voltage induced transition. This article is part of a Special Issue entitled: Membrane protein structure and function.

Toward understanding protocell mechanosensation

Mechanosensitive (MS) channels can prevent bacterial bursting during hypo-osmotic shocks by responding to increases in lateral tension at the membrane level through an integrated and coordinated opening mechanism. Mechanical regulation in protocells could have been one of the first mechanisms to evolve in order to preserve their integrity against changing environmental conditions. How has the rich functional diversity found in present cells been created throughout evolution, and what did the primordial MS channels look like? This review has been written with the aim of identifying which factors may have been important for the appearance of the first osmotic valve in a prebiotic context, and what this valve may have been like. It highlights the mechanical properties of lipid bilayers, the association of peptides as aggregates in membranes, and the conservation of sequence motifs as central aspects to understand the evolution of proteins that gate below the tension required for spontaneous pore formation and membrane rupture. The arguments developed here apply to both MscL and MscS homologs, but could be valid to mechano-susceptible proteins in general.

The structural and topological analysis of membrane-associated polypeptides by oriented solid-state NMR spectroscopy: established concepts and novel developments

Solid-state NMR spectroscopy is a powerful technique for the investigation of membrane-associated peptides and proteins as well as their interactions with lipids, and a variety of conceptually different approaches have been developed for their study. The technique is unique in allowing for the high-resolution investigation of liquid disordered lipid bilayers representing well the characteristics of natural membranes. Whereas magic angle solid-state NMR spectroscopy follows approaches that are related to those developed for solution NMR spectroscopy the use of static uniaxially oriented samples results in angular constraints which also provide information for the detailed analysis of polypeptide structures. This review introduces this latter concept theoretically and provides a number of examples. Furthermore, ongoing developments combining solid-state NMR spectroscopy with information from solution NMR spectroscopy and molecular modelling as well as exploratory studies using dynamic nuclear polarization solid-state NMR will be presented.

Antimicrobial peptides in toroidal and cylindrical pores

Antimicrobial peptides (AMPs) are small, usually cationic peptides, which permeabilize biological membranes. Their mechanism of action is still not well understood. Here we investigate the preference of alamethicin and melittin for pores of different shapes, using molecular dynamics (MD) simulations of the peptides in pre-formed toroidal and cylindrical pores. When an alamethicin hexamer is initially embedded in a cylindrical pore, at the end of the simulation the pore remains cylindrical or closes if glutamines in the N-termini are not located within the pore. On the other hand, when a melittin tetramer is embedded in toroidal pore or in a cylindrical pore, at the end of the simulation the pore is lined both with peptides and lipid headgroups, and, thus, can be classified as a toroidal pore. These observations agree with the prevailing views that alamethicin forms barrel-stave pores whereas melittin forms toroidal pores. Both alamethicin and melittin form amphiphilic helices in the presence of membranes, but their net charge differs; at pH approximately 7, the net charge of alamethicin is -1 whereas that of melittin is +5. This gives rise to stronger electrostatic interactions of melittin with membranes than those of alamethicin. The melittin tetramer interacts more strongly with lipids in the toroidal pore than in the cylindrical one, due to more favorable electrostatic interactions.

Residue-specific information about the dynamics of antimicrobial peptides from (1)H-(15)N and (2)H solid-state NMR spectroscopy

We present a new method to obtain information about the conformational dynamics of membrane-proteins using solid-state NMR experiments of oriented samples. By measuring the orientation-dependent (1)H-(15)N dipole-dipole coupling, (15)N anisotropic chemical shift, and (2)H quadrupole coupling parameters for a single residue, it is possible to obtain information about the local dynamics of each residue in the protein. This may be interpreted on an individual basis or through models extended to study conformational motion of membrane-protein segments. The method is demonstrated for the antimicrobial peptaibol alamethicin for which combined analysis of anisotropic interactions for the Aib(8) residue provides detailed information about helix-tilt angle, wobbling, and oscillatory rotation around the helix axis in the membrane bound state. This information is in very good agreement with coarse-grained MD simulations of the peptide in lipid bilayers.

Antimicrobial peptides bind more strongly to membrane pores

Antimicrobial peptides (AMPs) are small, usually cationic peptides, which permeabilize bacterial membranes. Understanding their mechanism of action might help design better antibiotics. Using an implicit membrane model, modified to include pores of different shapes, we show that four AMPs (alamethicin, melittin, a magainin analogue, MG-H2, and piscidin 1) bind more strongly to membrane pores, consistent with the idea that they stabilize them. The effective energy of alamethicin in cylindrical pores is similar to that in toroidal pores, whereas the effective energy of the other three peptides is lower in toroidal pores. Only alamethicin intercalates into the membrane core; MG-H2, melittin and piscidin are located exclusively at the hydrophobic/hydrophilic interface. In toroidal pores, the latter three peptides often bind at the edge of the pore, and are in an oblique orientation. The calculated binding energies of the peptides are correlated with their hemolytic activities. We hypothesize that one distinguishing feature of AMPs may be the fact that they are imperfectly amphipathic which allows them to bind more strongly to toroidal pores. An initial test on a melittin-based mutant seems to support this hypothesis.

Lipids on the move: simulations of membrane pores, domains, stalks and curves

In this review we describe the state-of-the-art of computer simulation studies of lipid membranes. We focus on collective lipid-lipid and lipid-protein interactions that trigger deformations of the natural lamellar membrane state, showing that many important biological processes including self-aggregation of membrane components into domains, the formation of non-lamellar phases, and membrane poration and curving, are now amenable to detailed simulation studies.

Pore formation induced by an antimicrobial peptide: electrostatic effects

We investigate the mode of action of Cateslytin, an antimicrobial peptide, on zwitterionic biomembranes by performing numerical simulations and electrophysiological measurements on membrane vesicles. Using this natural beta-sheet antimicrobial peptide secreted during stress as a model we show that a single peptide is able to form a stable membrane pore of 1 nm diameter of 0.25 nS conductance found both from calculation and electrical measurements. The resulting structure does not resemble the barrel-stave or carpet models earlier predicted, but is very close to that found in the simulation of alpha-helical peptides. Based on the simulation of a mutated peptide and the effects of small external electric fields, we conclude that electrostatic forces play a crucial role in the process of pore formation.

Simulations of lipid adsorption on TiO2 surfaces in solution

Molecular dynamics simulations are carried out to study the adsorption of three lipids, namely, DOPC, DOPS, and DMTAP, on TiO2(110) rutile surfaces and the influence of the interface on their conformational properties. Three types of rutile (110) surfaces, characterized by a different degree of hydroxylation (the neutral nonhydroxylated and hydroxylated surfaces and a partially hydroxylated surface with charge density corresponding to physiological pH) are investigated using force fields derived from ab initio calculations and experimental data. It is found that the stability of the adsorbate and the strength of the attachment are strictly connected with the nature of both the lipid and the surface. Direct coordination of the phosphate or carbonyl oxygens of the lipids with available titanium sites, observed in the case of partially or nonhydroxylated layers, determines stronger adsorption and, as a consequence, reduced dynamics. For a given hydration state of the surface, the adsorption strengths are in the order DOPS > DOPC >> DMTAP, in agreement with experimental data according to which the presence of DOPS units inside lipid bilayers favors stronger adsorption and lower mobility. The adsorption geometry, the hydration state of the lipid headgroups, and the dynamical processes (detachment, diffusion, etc.) occurring at the lipid/oxide interface are analyzed in detail, putting on a roughly quantitative basis time scales and energy barriers of the latter processes.

Molecular dynamics simulations of Na+/Cl(-)-dependent neurotransmitter transporters in a membrane-aqueous system

We have performed molecular dynamics simulations of a homology model of the human serotonin transporter (hSERT) in a membrane environment and in complex with either the natural substrate 5-HT or the selective serotonin reuptake inhibitor escitalopram. We have also included a transporter homologue, the Aquifex aeolicus leucine transporter (LeuT), in our study to evaluate the applicability of a simple and computationally attractive membrane system. Fluctuations in LeuT extracted from simulations are in good agreement with crystallographic B factors. Furthermore, key interactions identified in the X-ray structure of LeuT are maintained throughout the simulations indicating that our simple membrane system is suitable for studying the transmembrane protein hSERT in complex with 5-HT or escitalopram. For these transporter complexes, only relatively small fluctuations are observed in the ligand-binding cleft. Specific interactions responsible for ligand recognition, are identified in the hSERT-5HT and hSERT-escitalopram complexes. Our findings are in good agreement with predictions from mutagenesis studies.

Peptide aggregation and pore formation in a lipid bilayer: a combined coarse-grained and all atom molecular dynamics study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-micros CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 micros is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the alpha-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.

The membrane potential and its representation by a constant electric field in computer simulations

A theoretical framework is elaborated to account for the effect of a transmembrane potential in computer simulations. It is shown that a simulation with a constant external electric field applied in the direction normal to the membrane is equivalent to the influence of surrounding infinite baths maintained to a voltage difference via ion-exchanging electrodes connected to an electromotive force. It is also shown that the linearly-weighted displacement charge within the simulation system tracks the net flow of charge through the external circuit comprising the electromotive force and the electrodes. Using a statistical mechanical reduction of the degrees of freedom of the external system, three distinct theoretical routes are formulated and examined for the purpose of characterizing the free energy of a protein embedded in a membrane that is submitted to a voltage difference. The W-route is constructed from the variations in the voltage-dependent potential of mean force along a reaction path connecting two conformations of the protein. The Q-route is based on the average displacement charge as a function of the conformation of the protein. Finally, the G-route considers the relative charging free energy of specific residues, with and without applied membrane potentials. The theoretical formulation is illustrated with a simple model of an ion crossing a vacuum slab surrounded by two aqueous bulk phases and with a fragment of the voltage-sensor of the KvAP potassium channel.

Probing the conformation and dynamics of allatostatin neuropeptides: a structural model for functional differences

Allatostatins are a family of related neuropeptides that play an important role in development, reproduction, and digestion in insects. The cockroach Diploptera punctata has 13 allatostatin neuropeptides, with pleiotropic functions, two of which are: inhibition of juvenile hormone (JH) production and inhibition of gut muscle contraction. In this study, the conformation and dynamics of D. punctata allatostatin 5 (Dippu-AST 5) and allatostatin 8 (Dippu-AST 8) are investigated by CD, NMR, and molecular dynamics simulations. These peptides contain eight and nine residues, respectively, and the identical six-residue C-terminal motif. Yet Dippu-AST 5 and Dippu-AST 8 affect juvenile hormone production and hindgut contraction with different potencies. Dippu-AST 5 is one of the most potent inhibitors of juvenile hormone production and one of the least potent inhibitors of gut contraction, whereas Dippu-AST 8 has the opposite potencies with respect to these tissues. From the NMR structure, it is clear that Dippu-AST 5 has a 3(10) helix involving three of its residues and a "gamma" turn at the end of its C-terminal motif. In contrast Dippu-AST 8 has an open "pi" turn among five of its central residues. In addition, the orientation preferences within the membrane of the two peptides were simulated. Our simulation results show that the C-terminal segment of Dippu-AST 5 orients in the membrane surface with an average angle of 17.5 degrees, whereas Dippu-AST 8 orients with an average angle of 5.1 degrees. Taken together, from the structures and orientation preferences of these peptides within the membrane, it appears that these peptides may interact with the receptor very differently.

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

The recombinant HIV-1 Tat protein contains a small region corresponding to residues (47)YGRKKRRQRR(57)R, which is capable of translocating cargoes of different molecular sizes, such as proteins, DNA, RNA, or drugs, across the cell membrane in an apparently energy-independent manner. The pathway that these peptides follow for entry into the cell has been the subject of strong controversy for the last decade. This peptide is highly basic and hydrophilic. Therefore, a central question that any candidate mechanism has to answer is how this highly hydrophilic peptide is able to cross the hydrophobic barrier imposed by the cell membrane. We propose a mechanism for the spontaneous translocation of the Tat peptides across a lipid membrane. This mechanism involves strong interactions between the Tat peptides and the phosphate groups on both sides of the lipid bilayer, the insertion of charged side chains that nucleate the formation of a transient pore, followed by the translocation of the Tat peptides by diffusing on the pore surface. This mechanism explains how key ingredients, such as the cooperativity among the peptides, the large positive charge, and specifically the arginine amino acids, contribute to the uptake. The proposed mechanism also illustrates the importance of membrane fluctuations. Indeed, mechanisms that involve large fluctuations of the membrane structure, such as transient pores and the insertion of charged amino acid side chains, may be common and perhaps central to the functions of many membrane protein functions.

Nanosecond field alignment of head group and water dipoles in electroporating phospholipid bilayers

To investigate the mechanism of biological cell membrane electroporation at the nanosecond and nanometer scale, we tracked pore-forming lipids and water in molecular dynamics simulations of a palmitoyloleoylphosphatidylcholine bilayer in a minimum porating electric field. Although the field-generated torque tilts the mean head group dipole a few degrees away from its equilibrium, zero-field position relative to the bilayer plane, this change in conformation does not appear to contribute directly to the development of the pore-initiating aggregation of lipid head groups and water that leads to the formation of a membrane-spanning hydrophilic pore. Field-directed rotation of the head group dipoles in the plane of the incipient pore wall, in combination with water dipole and solvation interactions at the aqueous-lipid interface, is one component in the coordinated ensemble of electroporation events.

Computational modeling of poly(alkylthiophene) conductive polymer insertion into phospholipid bilayers

We have previously demonstrated that some poly(alkylthiophenes) (PATs) are able to increase the electrical conductance of unsupported phospholipid bilayers and have hypothesized that this effect is due to the ability of some PAT side chains to permit stable insertion into the bilayer. We have further proposed the development of long-term intracellular electrodes based on that phenomenon. In this article, we apply molecular dynamics techniques to study the insertion of two model PATs into a patch of a lipid bilayer. Steered molecular dynamics is used to obtain potential trajectories of insertion, followed by umbrella sampling to determine the free-energy change upon insertion. Our results indicate that both branched-side-chain poly(3-(2-ethylhexyl)thiophene) (EHPT) and straight-side-chain poly(3-hexylthiophene) (HPT) are able to enter the bilayer but only EHPT can cross the center of the membrane and establish an electrical bridge. HPT penetrates the head groups but is not able to enter the alkyl tail phase. These findings support the feasibility of our electrode concept and raise questions regarding the mechanisms by which branched side chains grant PATs greater solubility in a lipid bilayer environment. The parameters and methods used in this study establish a novel framework for studying these and similar systems, and the results hold promise for the use of EHPT in biosensing and neural interfacing.

All-atom molecular dynamics simulations using orientational constraints from anisotropic NMR samples

Orientational constraints obtained from solid state NMR experiments on anisotropic samples are used here in molecular dynamics (MD) simulations for determining the structure and dynamics of several different membrane-bound molecules. The new MD technique is based on the inclusion of orientation dependent pseudo-forces in the COSMOS-NMR force field. These forces drive molecular rotations and re-orientations in the simulation, such that the motional time-averages of the tensorial NMR properties approach the experimentally measured parameters. The orientational-constraint-driven MD simulations are universally applicable to all NMR interaction tensors, such as chemical shifts, dipolar couplings and quadrupolar interactions. The strategy does not depend on the initial choice of coordinates, and is in principle suitable for any flexible molecule. To test the method on three systems of increasing complexity, we used as constraints some deuterium quadrupolar couplings from the literature on pyrene, cholesterol and an antimicrobial peptide embedded in oriented lipid bilayers. The MD simulations were able to reproduce the NMR parameters within experimental error. The alignment of the three membrane-bound molecules and some aspects of their conformation were thus derived from the NMR data, in good agreement with previous analyses. Furthermore, the new approach yielded for the first time the distribution of segmental orientations with respect to the membrane and the order parameter tensors of all three systems.