Search

Free energy profiles for insertion of a hydrophobic transmembrane protein α-helix (M2 from CFTR) into a lipid bilayer have been calculated using coarse-grained molecular dynamics simulations and umbrella sampling to yield potentials of mean force along a reaction path corresponding to translation of a helix across a lipid bilayer. The calculated free energy of insertion is smaller when a bilayer with a thinner hydrophobic region is used. The free energies of insertion from the potentials of mean force are compared with those derived from a number of hydrophobicity scales and with those derived from translocon-mediated insertion. This comparison supports recent models of translocon-mediated insertion and in particular suggests that: 1), helices in an about-to-be-inserted state may be located in a hydrophobic region somewhat thinner than the core of a lipid bilayer; and/or 2), helices in a not-to-be-inserted state may experience an environment more akin (e.g., in polarity/hydrophobicity) to the bilayer/water interface than to bulk water.

Computing free energies of complex biomolecular systems via atomistic (AT) molecular dynamics (MD) simulations remains a challenge due to the need for adequate sampling and convergence. Recent coarse-grained (CG) methodology allows simulations of significantly larger systems (~106 to 108 atoms) over longer (s/ms) time scales. Such CG models appear to be capable of making semiquantitative predictions. However, their ability to reproduce accurate thermodynamic quantities remains uncertain. We have recently used CG MD simulations to compute the potential of mean force (PMF) or free energy profile of a small peptide toxin interacting with a lipid bilayer along a 1D reaction coordinate. The toxin studied was VSTx1 (Voltage Sensor Toxin 1) from spider venom which inhibits the archeabacterial voltage-gated potassium (Kv) channel KvAP by binding to the voltage-sensor (VS) domains. Here, we reestimate this PMF profile using (i) AT MD simulations with explicit membrane and solvent and (ii) an implicit membrane and solvent (generalized Born; GBIM) model where only the peptide was explicit. We used the CG MD free energy simulations to guide the setup of the corresponding AT MD simulations. The aim was to avoid local minima in the AT simulations which would be difficult over shorter AT time scales. A cross-comparison of the PMF profiles revealed a conserved topology, although there were differences in the magnitude of the free energies. The CG and AT simulations predicted a membrane/water interface free energy well of -27 and -23 kcal/mol, respectively (with respect to water). The GBIM model, however, gave a reduced interfacial free energy well (-12 kcal/mol). In addition, the CG and GBIM models predicted a free energy barrier of +61 and +96 kcal/mol, respectively, for positioning the toxin at the center of the bilayer, which was considerably smaller in the AT simulations (+26 kcal/mol). Thus, we present a framework for serially combining CG and AT simulations to estimate the free energy of peptide/membrane interactions. Such approaches for combining simulations at different levels of granularity will become increasingly important in future studies of complex membrane/protein systems. Copyright © 2010 American Chemical Society.

The outer membrane protein OprP from Pseudomonas aeruginosa forms a phosphate selective pore. To understand the mechanism of phosphate permeation and selectivity, we used three simulation techniques [equilibrium molecular dynamics simulations, steered molecular dynamics, and calculation of a potential of mean force (PMF)]. The PMF for phosphate reveals a deep free energy well midway along the OprP channel. Two adjacent phosphate-binding sites (W1 and W2), each with a well depth of ≈8 kT , are identified close to the L3 loop in the most constricted region of the pore. A dissociation constant for phosphate of 6 μM is computed from the PMF, within the range of reported experimental values. The transfer of phosphate between sites W1 and W2 is correlated with changes in conformation of the sidechain of K121, which serves as a “charged brush” to facilitate phosphate passage between the two subsites. OprP also binds chloride, but less strongly than phosphate, as calculated from a Cl − PMF. The difference in affinity and hence selectivity is due to the “Lys-cluster” motif, the positive charges of which interact strongly with a partially dehydrated phosphate ion but are shielded from a Cl − by the hydration shell of the smaller ion. Our simulations suggest that OprP does not conform to the conventional picture of a channel with relatively flat energy landscape for permeant ions, but rather resembles a membrane-inserted binding protein with a high specificity that allows access to a centrally located binding site from both the extracellular and the periplasmic spaces.

A number of membrane proteins act via binding at the water/lipid bilayer interface. An important example of such proteins is provided by the gating-modifier toxins that act on voltage-gated potassium (Kv) channels. They are thought to partition to the headgroup region of lipid bilayers, and so provide a good system for probing the nature of interactions of a protein with the water/bilayer interface. We used coarse-grained molecular dynamics simulations to compute the one-dimensional potential of mean force (i.e., free energy) profile that governs the interaction between a Kv channel gating-modifier toxin (VSTx1) and model phospholipid bilayers. The reaction coordinate sampled corresponds to the position of the toxin along the bilayer normal. The course-grained representation of the protein and lipids enabled us to explore extended time periods, revealing aspects of toxin/bilayer dynamics and energetics that would be difficult to observe on the timescales currently afforded by atomistic molecular dynamics simulations. In particular, we show for this model system that the bilayer deforms as it interacts with the toxin, and that such deformations perturb the free energy profile. Bilayer deformation therefore adds an additional layer of complexity to be addressed in investigations of membrane/protein systems. In particular, one should allow for local deformations that may arise due to the spatial array of charged and hydrophobic elements of an interfacially located membrane protein.

A number of membrane-active enzymes act in a complex environment formed by the interface between a lipid bilayer and bulk water. Although x-ray diffraction studies yield structures of isolated enzyme molecules, a detailed characterization of their interactions with the interface requires a measure of how deeply such a membrane-associated protein penetrates into a lipid bilayer. Here, we apply coarse-grained (CG) molecular dynamics (MD) simulations to probe the interaction of porcine pancreatic phospholipase A2 (PLA2) with a lipid bilayer containing palmitoyl-oleoyl-phosphatidyl choline and palmitoyl-oleoyl-phosphatidyl glycerol molecules. We also used a configuration from a CG-MD trajectory to initiate two atomistic (AT) MD simulations. The results of the CG and AT simulations are evaluated by comparison with available experimental data. The membrane-binding surface of PLA2 consists of a patch of hydrophobic residues surrounded by polar and basic residues. We show this proposed footprint interacts preferentially with the anionic headgroups of the palmitoyl-oleoyl-phosphatidyl glycerol molecules. Thus, both electrostatic and hydrophobic interactions determine the location of PLA2 relative to the bilayer. From a general perspective, this study demonstrates that CG-MD simulations may be used to reveal the orientation and location of a membrane-surface-bound protein relative to a lipid bilayer, which may subsequently be refined by AT-MD simulations to probe more detailed interactions.

Gating-modifier toxins inhibit voltage-gated ion channels by binding the voltage sensors (VS) and altering the energetics of voltage-dependent gating. These toxins are thought to gain access to the VS via the membrane (i.e., by partitioning from water into the membrane before binding the VS). We used serial multiscale molecular-dynamics (MD) simulations, via a combination of coarse-grained (CG) and atomistic (AT) simulations, to study how the toxin VSTx1, which inhibits the archeabacterial voltage-gated potassium channel KvAP, interacts with an isolated membrane-embedded VS domain. In the CG simulations, VSTx1, which was initially located in water, partitioned into the headgroup/water interface of the lipid bilayer before binding the VS. The CG configurations were used to generate AT representations of the system, which were subjected to AT-MD to further evaluate the stability of the complex and refine the predicted VS/toxin interface. VSTx1 interacted with a binding site on the VS formed by the C-terminus of S1, the S1-S2 linker, and the N-terminus of S4. The predicted VS/toxin interactions are suggestive of toxin-mediated perturbations of the interaction between the VS and the pore domain of Kv channels, and of the membrane. Our simulations support a membrane-access mechanism of inhibition of Kv channels by VS toxins. Overall, the results show that serial multiscale MD simulations may be used to model a two-stage process of protein-bilayer and protein-protein interactions within a membrane.

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.

Experimental and computational studies have indicated that hydrophobicity plays a key role in driving the insertion of transmembrane alpha-helices into lipid bilayers. Molecular dynamics simulations allow exploration of the nature of the interactions of transmembrane alpha-helices with their lipid bilayer environment. In particular, coarse-grained simulations have considerable potential for studying many aspects of membrane proteins, ranging from their self-assembly to the relation between their structure and function. However, there is a need to evaluate the accuracy of coarse-grained estimates of the energetics of transmembrane helix insertion. Here, three levels of complexity of model system have been explored to enable such an evaluation. First, calculated free energies of partitioning of amino acid side chains between water and alkane yielded an excellent correlation with experiment. Second, free energy profiles for transfer of amino acid side chains along the normal to a phosphatidylcholine bilayer were in good agreement with experimental and atomistic simulation studies. Third, estimation of the free energy profile for transfer of an arginine residue, embedded within a hydrophobic alpha-helix, to the center of a lipid bilayer gave a barrier of approximately 15 kT. Hence, there is a substantial barrier to membrane insertion for charged amino acids, but the coarse-grained model still underestimates the corresponding free energy estimate (approximately 29 kT) from atomistic simulations (Dorairaj, S., and Allen, T. W. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 4943-4948). Coarse-grained simulations were then used to predict the free energy profile for transfer of a simple model transmembrane alpha-helix (WALP23) across a lipid bilayer. The results indicated that a transmembrane orientation was favored by about -70 kT.

SGTx1 is a gating-modifier toxin that has been shown to inhibit the voltage-gated potassium channel Kv2.1. SGTx1 is thought to bind to the S3b-S4a region of the voltage-sensor, and is believed to alter the energetics of gating. Gating-modifier toxins such as SGTx1 are of interest as they can be used to probe the structure and dynamics of their target channels. Although there are experimental data for SGTx1, its interaction with lipid bilayer membranes remains to be characterized. We performed atomistic and coarse-grained molecular dynamics simulations to study the interaction of SGTx1 with a POPC and a 3:1 POPE/POPG lipid bilayer membrane. We reveal the preferential partitioning of SGTx1 into the water/membrane interface of the bilayer. We also show that electrostatic interactions between the charged residues of SGTx1 and the lipid headgroups play an important role in stabilizing SGTx1 in a bilayer environment. © 2007 by the Biophysical Society.

VSTx1 is a tarantula venom toxin which binds to the archaebacterial voltage-gated potassium channel KvAP. VSTx1 is thought to access the voltage sensor domain of the channel via the lipid bilayer phase. In order to understand its mode of action and implications for the mechanism of channel activation, it is important to characterize the interactions of VSTx1 with lipid bilayers. Molecular dynamics (MD) simulations (for a total simulation time in excess of 0.2 micros) have been used to explore VSTx1 localization and interactions with zwitterionic (POPC) and with anionic (POPE/POPG) lipid bilayers. In particular, three series of MD simulations have been used to explore the net drift of VSTx1 relative to the center of a bilayer, starting from different locations of the toxin. The preferred location of the toxin is at the membrane/water interface. Although there are differences between POPC and POPE/POPG bilayers, in both cases the toxin forms favorable interactions at the interface, maximizing H-bonding to lipid headgroups and to water molecules while retaining interactions with the hydrophobic core of the bilayer. A 30 ns unrestrained simulation reveals dynamic partitioning of VSTx1 into the interface of a POPC bilayer. The preferential location of VSTx1 at the interface is discussed in the context of Kv channel gating models and provides support for a mode of action in which the toxin interacts with the Kv voltage sensor "paddle" formed by the S3 and S4 helices.

The outer membrane protein OprP from Pseudomonas aeruginosa forms an anion-selective pore, especially selective for phosphate ions. The protein is homo-trimeric, with each pore lined by three positively charged loops (L3, L5, and T7) folded into its lumen. OprP plays a key role in high-affinity phosphate uptake under the condition of phosphate starvation. To better understand the mechanism of phosphate-selective permeation, we employed three simulation techniques: (i) equilibrium molecular dynamics simulations (MD); (ii) steered MD (SMD); (iii) umbrella sampling to calculate a potential of mean force (PMF) for phosphate and chloride ions. The PMFs reveal a deep energy well midway along the OprP channel. Two adjacent phosphate-binding sites (W1 and W2), each with a well depth of ∼8kT, are identified close to the L3 loop in the most constricted region of the pore. The transfer of phosphate between sites W1 and W2 is correlated with changes in conformation of the sidechain of K121, which serves as a ‘charged brush’ to facilitate phosphate passage between the two subsites. The PMF for chloride has also been computed and can be compared with that of phosphate. Our simulations suggest that OprP does not conform to the conventional picture of a channel with a relatively flat energy landscape for permeant ions, but rather resembles a membrane-inserted binding protein with a high specificity that allows access to a centrally located binding site from both the extracellular and the periplasmic spaces.

Coarse-grained molecular dynamics simulations have been used to explore the interactions of C(60) and its derivatives with lipid bilayers. Pristine C(60) partitions into the bilayer core, whilst C(60)(OH)(20) experiences a central energetic barrier to permeation across the bilayer. For intermediate levels of derivatization, e.g. C(60)(OH)(10), this central barrier is smaller and there is an energetic well at the bilayer/water interface, thus promoting entry into cells via bilayer permeation whilst maintaining solubility in water.

The S4 helix is a major element of the voltage-sensor of voltage-sensitive ion channels. This helix contains an array of positively charged sidechains and yet adopts a transmembrane orientation within the voltage sensor of a voltage-gated channel. Thus, from both mechanistic and a biosynthetic perspectives, the question of how the S4 helix may be stabilized in a membrane environment is of some importance. We have performed coarse-grained (CG) molecular dynamics (MD) simulations to calculate: (1) the free energy of insertion of a S4 helix; and (2) the free energy cost of driving a S4 helix through an angular motion in model membranes. Our results suggest that it is possibly to meta-stably insert a S4 helix in a TM orientation in a lipid bilayer. In this orientation, the helix is stabilized local bilayer deformation and by snorkelling of the sidechains of the positively-charged residues of S4 to interact with lipid phosphates and waters.