Chze Ling Wee

Coarse-grained molecular dynamics simulations of the energetics of helix insertion into a lipid bilayer.

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.

The interaction of C60 and its derivatives with a lipid bilayer via molecular dynamics simulations

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.

Simulations of anion transport through OprP reveal the molecular basis for high affinity and selectivity for phosphate

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.

The interaction of phospholipase A2 with a phospholipid bilayer: coarse-grained molecular dynamics simulations.

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.

Interactions Between a Voltage Sensor and a Toxin via Multiscale Simulations

Abstract 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.

Lipid Bilayer Deformation and the Free Energy of Interaction of a Kv Channel Gating-Modifier Toxin ☆

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.

The Energetics of Transmembrane Helix Insertion into a Lipid Bilayer

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.

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.

Membrane/Toxin Interaction Energetics via Serial Multiscale Molecular Dynamics Simulations.

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 (∼10(6) to 10(8) 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 re-estimate 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.