Double bilayers and transmembrane gradients: a molecular dynamics study of a highly charged peptide

Translocation Processes of Pt(II)-Based Drugs through Human Breast Cancer Cell Membrane: In Silico Experiments

Breast cancer is one of the most frequent modalities of cancer worldwide, with notable mortality. The medication based on platinum drugs (cisplatin (cddp), carboplatin (cpx), and oxaliplatin (oxa)) is a conventional chemotherapy despite severe side effects and the development of drug resistance. In order to provide a deeper molecular description of the influx and efflux processes of platinum drugs through breast cancer tissues, this study focuses on molecular dynamics (MD) simulations of the passive translocation process through a realistic plasma membrane prototype of human breast cancer cell (c_memb). The results showed that the permeation events were mainly mediated by neutral lipids (DOPC, DOPE, and cholesterol), producing a low and temporary membrane deformation. The drug insertion in the region of polar heads was the most favorable stage of the translocation mechanism, especially for cddp and oxa with potential wells of -8.6 and -9.8 kcal mol-1, respectively. However, the potentials of mean force (PMF) revealed unfavorable kinetics for the permeation of these drugs through lipid tails, with energy barriers of 28.3 (cddp), 32.2 (cpx), and 30.4 kcal mol-1 (oxa). The low permeability coefficients (P) of cpx and oxa, which were 3 and 1 orders of magnitude inferior than for cddp, resulted from the high energy barriers for their traslocation processes through the membrane. The obtained results provide a more accurate picture of the permeation of Pt(II)-based drugs through breast cancer cells, which may be relevant for the design and evaluation of new platinum complexes.

Efficient Simulations of Solvent Asymmetry Across Lipid Membranes Using Flat-Bottom Restraints

The routinely employed periodic boundary conditions complicate molecular simulations of physiologically relevant asymmetric lipid membranes together with their distinct solvent environments. Therefore, separating the extracellular fluid from its cytosolic counterpart has often been performed using a costly double-bilayer setup. Here, we demonstrate that the lipid membrane and solvent asymmetry can be efficiently modeled with a single lipid bilayer by applying an inverted flat-bottom potential to ions and other solute molecules, thereby restraining them to only interact with the relevant leaflet. We carefully optimized the parameters of the suggested method so that the results obtained using the flat-bottom and double-bilayer approaches become mutually indistinguishable. Then, we apply the flat-bottom approach to lipid bilayers with various compositions and solvent environments, covering ions and cationic peptides to validate the approach in a realistic use case. We also discuss the possible limitations of the method as well as its computational efficiency and provide a step-by-step guide on how to set up such simulations in a straightforward manner.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

Functional Annotation of Ion Channel Structures by Molecular Simulation

Ion channels play key roles in cell membranes, and recent advances are yielding an increasing number of structures. However, their functional relevance is often unclear and better tools are required for their functional annotation. In sub-nanometer pores such as ion channels, hydrophobic gating has been shown to promote dewetting to produce a functionally closed (i.e., non-conductive) state. Using the serotonin receptor (5-HT₃R) structure as an example, we demonstrate the use of molecular dynamics to aid the functional annotation of channel structures via simulation of the behavior of water within the pore. Three increasingly complex simulation analyses are described: water equilibrium densities; single-ion free-energy profiles; and computational electrophysiology. All three approaches correctly predict the 5-HT₃R crystal structure to represent a functionally closed (i.e., non-conductive) state. We also illustrate the application of water equilibrium density simulations to annotate different conformational states of a glycine receptor.

Energetics of Multi-Ion Conduction Pathways in Potassium Ion Channels

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

Membrane insertion of a voltage sensor helix

Most membrane proteins contain a transmembrane (TM) domain made up of a bundle of lipid-bilayer-spanning α-helices. TM α-helices are generally composed of a core of largely hydrophobic amino acids, with basic and aromatic amino acids at each end of the helix forming interactions with the lipid headgroups and water. In contrast, the S4 helix of ion channel voltage sensor (VS) domains contains four or five basic (largely arginine) side chains along its length and yet adopts a TM orientation as part of an independently stable VS domain. Multiscale molecular dynamics simulations are used to explore how a charged TM S4 α-helix may be stabilized in a lipid bilayer, which is of relevance in the context of mechanisms of translocon-mediated insertion of S4. Free-energy profiles for insertion of the S4 helix into a phospholipid bilayer suggest that it is thermodynamically favorable for S4 to insert from water to the center of the membrane, where the helix adopts a TM orientation. This is consistent with crystal structures of Kv channels, biophysical studies of isolated VS domains in lipid bilayers, and studies of translocon-mediated S4 helix insertion. Decomposition of the free-energy profiles reveals the underlying physical basis for TM stability, whereby the preference of the hydrophobic residues of S4 to enter the bilayer dominates over the free-energy penalty for inserting charged residues, accompanied by local distortion of the bilayer and penetration of waters. We show that the unique combination of charged and hydrophobic residues in S4 allows it to insert stably into the membrane.

APBSmem: a graphical interface for electrostatic calculations at the membrane

Electrostatic forces are one of the primary determinants of molecular interactions. They help guide the folding of proteins, increase the binding of one protein to another and facilitate protein-DNA and protein-ligand binding. A popular method for computing the electrostatic properties of biological systems is to numerically solve the Poisson-Boltzmann (PB) equation, and there are several easy-to-use software packages available that solve the PB equation for soluble proteins. Here we present a freely available program, called APBSmem, for carrying out these calculations in the presence of a membrane. The Adaptive Poisson-Boltzmann Solver (APBS) is used as a back-end for solving the PB equation, and a Java-based graphical user interface (GUI) coordinates a set of routines that introduce the influence of the membrane, determine its placement relative to the protein, and set the membrane potential. The software Jmol is embedded in the GUI to visualize the protein inserted in the membrane before the calculation and the electrostatic potential after completing the computation. We expect that the ease with which the GUI allows one to carry out these calculations will make this software a useful resource for experimenters and computational researchers alike. Three examples of membrane protein electrostatic calculations are carried out to illustrate how to use APBSmem and to highlight the different quantities of interest that can be calculated.

From the gating charge response to pore domain movement: initial motions of Kv1.2 dynamics under physiological voltage changes

Recent structures of the potassium channel provide an essential beginning point for explaining how the pore is gated between open and closed conformations by changes in membrane voltage. Yet, the molecular details of this process and the connections to transmembrane gradients are not understood. To begin addressing how changes within a membrane environment lead to the channel's ability to sense shifts in membrane voltage and to gate, we performed double-bilayer simulations of the Kv1.2 channel. These double-bilayer simulations enable us to simulate realistic voltage drops from resting potential conditions to depolarized conditions by changes in the bath conditions on each side of the bilayer. Our results show how the voltage sensor domain movement responds to differences in transmembrane potential. The initial voltage sensor domain movement, S4 in particular, is modulated by the gating charge response to changes in voltage and is initially stabilized by the lipid headgroups. We show this response is directly coupled to the initial stages of pore domain motion. Results presented here provide a molecular model for how the pre-gating process occurs in sequential steps: Gating charge response, movement and stabilization of the S4 voltage sensor domain, and movement near the base of the S5 region to close the pore domain.

Cause and effect of melittin-induced pore formation: a computational approach

Melittin embedded in a palmitoyl oleyl phosphatidylcholine bilayer at a high peptide/lipid ratio (1:30) was simulated in the presence of explicit water and ions. The simulation results indicate the incipience of an ion-permeable water pore through collective membrane perturbation by bound peptides. The positively charged residues of melittin not only act as "anchors" but also disrupt the membrane, leading to cell lysis. A detailed analysis of the lipid tail order parameter profile depicts localized membrane perturbation. The lipids in the vicinity of the aqueous cavity adopt a tilted conformation, which allows local bilayer thinning. The prepore thus formed can be considered as the melittin-induced structural defects in the bilayer membrane. Because of the strong cationic nature, the melittin-induced prepore exhibits selectivity toward anions over cations. As Cl(-) ions entered into the prepore, they are electrostatically entrapped by positively charged residues located at its wall. The confined motion of the Cl(-) ions in the membrane interior is obvious from calculated diffusion coefficients. Moreover, reorientation of the local lipids occurs in such a way that few lipid heads along with peptide helices can line the surface of the penetrating aqueous phase. The flipping of lipids argued in favor of melittin-induced toroidal pore over a barrel-stave mechanism. Thus, our result provides atomistic level details of the mechanism of membrane disruption by antimicrobial peptide melittin.