Molecular dynamics studies on structure and dynamics of phospholamban monomer and pentamer in membranes

Phospholamban (PLB) is an integral membrane protein of 52 residues that regulates the activity of the sarcoplasmic reticulum calcium pump in cardiac muscle cells through reversible phosphorylation of Ser16. To explore its possible conformations and dynamics in a monomeric state, we have performed comparative molecular dynamics simulations of unphosphorylated and phosphorylated PLB (pPLB) with various orientations in POPC membranes. The simulations indicate that dynamics of the cytoplasmic domain is highly dependent on its interactions with membranes, that is, large conformational changes in the absence of membrane interactions, but very restricted dynamics in their presence. pPLB shows more structural flexibility in its cytoplasmic domain, which is consistent with experimental observations. We have also performed a simulation of a PLB pentameric structure (the so‐called bellflower model), recently determined in micelles, to investigate its behaviors in a POPC membrane. The cytoplasmic domain in each monomer shows uncorrelated dynamics and undergoes large conformational changes toward the membrane surface during the simulation, which supports the so‐called pinwheel model of the PLB pentamer structure. The hydrophobic nature of the pentameric pore excludes water molecules in the pore region, which illustrates that the pore appears to be an energetic barrier for ion and water translocation. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Role of conformational sampling of Ser16 and Thr17-phosphorylated phospholamban in interactions with SERCA

Phosphorylation of phospholamban (PLB) at Ser16 and/ or Thr17 is believed to release its inhibitory effect on sarcoplasmic reticulum calcium ATPase. Ser16 phosphorylation of PLB has been suggested to cause a conformational change that alters the interaction between the enzyme and protein. Using computer simulations, the conformational sampling of Ser16 phosphorylated PLB in implicit membrane environment is compared here with the unphosphorylated PLB system to investigate these conformational changes. The results suggest that conformational changes in the cytoplasmic domain of PLB upon phosphorylation at Ser16 increase the likelihood of unfavorable interactions with SERCA in the E2 state prompting a conformational switch of SERCA from E2 to E1. Phosphorylation of PLB at Thr17 on the other hand does not appear to affect interactions with SERCA significantly suggesting that the mechanism of releasing the inhibitory effect is different between Thr17 phosphorylated and Ser16 phosphorylated PLB.     

Sampling efficiency in explicit and implicit membrane environments studied by peptide folding simulations

Efficient exploration of the conformational space of peptides embedded in biological membranes is vital to extract converged thermodynamic and kinetic data from simulation studies. Recently developed implicit membrane models promise vast increases in sampling efficiency compared to explicit membrane simulations, allowing for ab initio structure prediction and functional studies. In this study, a previously developed implicit membrane model, based on the generalized Born method, is compared to an explicit di‐palmitoyl‐phosphatidyl‐choline lipid bilayer and an octane slab membrane mimic. The complete folding process of a synthetic 16‐residue peptide is compared using these three setups. Since the comparison requires the entire folding pathway to be captured, individual simulations ranged up to 3 μs of MD. A quantitative sampling comparison using a wide range of performance metrics reveals that the implicit membrane model is at least 2 orders of magnitude more efficient than the simplest explicit setups. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

A Practical Implicit Membrane Potential for NMR Structure Calculations of Membrane Proteins

The highly anisotropic environment of the lipid bilayer membrane imposes significant constraints on the structures and functions of membrane proteins. However, NMR structure calculations typically use a simple repulsive potential that neglects the effects of solvation and electrostatics, because explicit atomic representation of the solvent and lipid molecules is computationally expensive and impractical for routine NMR-restrained calculations that start from completely extended polypeptide templates. Here, we describe the extension of a previously described implicit solvation potential, eefxPot, to include a membrane model for NMR-restrained calculations of membrane protein structures in XPLOR-NIH. The key components of eefxPot are an energy term for solvation free energy that works together with other nonbonded energy functions, a dedicated force field for conformational and nonbonded protein interaction parameters, and a membrane function that modulates the solvation free energy and dielectric screening as a function of the atomic distance from the membrane center, relative to the membrane thickness. Initial results obtained for membrane proteins with structures determined experimentally in lipid bilayer membranes show that eefxPot affords significant improvements in structural quality, accuracy, and precision. Calculations with eefxPot are straightforward to implement and can be used to both fold and refine structures, as well as to run unrestrained molecular-dynamics simulations. The potential is entirely compatible with the full range of experimental restraints measured by various techniques. Overall, it provides a useful and practical way to calculate membrane protein structures in a physically realistic environment.     

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys⁴⁰, residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.     

The role of phosphorylation on the structure and dynamics of phospholamban: a model from molecular simulations

Phospholamban (PLB) is a small membrane protein that regulates the activity of the calcium ATP-ase in the cardiac, slow-twitch, and smooth muscle sarcoplasmic reticulum through the reversible phosphorylation of Ser16. We present here a comparative molecular dynamics study of unmodified and phosphorylated PLB immersed in a phospholipid membrane. The study has been performed under different ionic strength conditions, using the NMR structures of two PLB variants determined in mixed organic solvent and dodecylphosphocholine micelles. The simulations indicate that all PLB forms studied display a highly dynamic behavior of the N-terminal cytoplasmic moiety, with a decrease of its helical content in the phosphorylated forms. The cytoplasmic domain undergoes large collective motions sampling conformations parallel as well as perpendicular to the membrane surface in all the simulations. The transmembrane domain retains a tightly folded helical conformation with a small tilt with respect to the membrane plane probably induced by the presence of Asn30 and Asn34 within the hydrophobic environment. Furthermore, the phosphoric group on Ser16 establishes transient electrostatic interactions with the phospholipid heads. We propose a model in which phosphorylation diminishes the probability of interactions of PLB with residues near Lys400 in the SERCA pump, thus relieving its inhibition.     

Numerical studies of the membrane fluorescent dyes dynamics in ground and excited states

Fluorescence methods are widely used in studies of biological and model membranes. The dynamics of membrane fluorescent markers in their ground and excited electronic states and correlations with their molecular surrounding within the fully hydrated phospholipid bilayer are still not well understood. In the present work, Quantum Mechanical (QM) calculations and Molecular Dynamics (MD) simulations are used to characterize location and interactions of two membrane polarity probes (Prodan; 6-propionyl-2-dimethylaminonaphthalene and its derivative Laurdan; 2-dimethylamino-6-lauroylnaphthalene) with the dioleoylphosphatidylcholine (DOPC) lipid bilayer model. MD simulations with fluorophores in ground and excited states are found to be a useful tool to analyze the fluorescent dye dynamics and their immediate vicinity. The results of QM calculations and MD simulations are in excellent agreement with available experimental data. The calculation shows that the two amphiphilic dyes initially placed in bulk water diffuse within 10 ns towards their final location in the lipid bilayer. Analysis of solvent relaxation process in the aqueous phase occurs on the picoseconds timescale whereas it takes nanoseconds at the lipid/water interface. Four different relaxation time constants, corresponding to different relaxation processes, where observed when the dyes were embedded into the membrane.     

Membrane Curvature Sensing by Amphipathic Helices: Insights from Implicit Membrane Modeling

Sensing and generation of lipid membrane curvature, mediated by the binding of specific proteins onto the membrane surface, play crucial roles in cell biology. A number of mechanisms have been proposed, but the molecular understanding of these processes is incomplete. All-atom molecular dynamics simulations have offered valuable insights but are extremely demanding computationally. Implicit membrane simulations could provide a viable alternative, but current models apply only to planar membranes. In this work, the implicit membrane model 1 is extended to spherical and tubular membranes. The geometric change from planar to curved shapes is straightforward but insufficient for capturing the full curvature effect, which includes changes in lipid packing. Here, these packing effects are taken into account via the lateral pressure profile. The extended implicit membrane model 1 is tested on the wild-types and mutants of the antimicrobial peptide magainin, the ALPS motif of arfgap1, α-synuclein, and an ENTH domain. In these systems, the model is in qualitative agreement with experiments. We confirm that favorable electrostatic interactions tend to weaken curvature sensitivity in the presence of strong hydrophobic interactions but may actually have a positive effect when those are weak. We also find that binding to vesicles is more favorable than binding to tubes of the same diameter and that the long helix of α-synuclein tends to orient along the axis of tubes, whereas shorter helices tend to orient perpendicular to it. Adoption of a specific orientation could provide a mechanism for coupling protein oligomerization to tubule formation.     

Predicting large-scale conformational changes in proteins using energy-weighted normal modes

We report the development of a method to improve the sampling of protein conformational space in molecular simulations. It is shown that a principal component analysis of energy-weighted normal modes in Cartesian coordinates can be used to extract vectors suitable for describing the dynamics of protein substructures. The method can operate with either atomistic or user-defined coarse-grained models of protein structure. An implicit reverse coarse-graining allows the dynamics of all-atoms to be recovered when a coarse-grained model is used. For an external test set of four proteins, it is shown that the new method is more successful than normal mode analysis in describing the large-scale conformational changes observed on ligand binding. The method has potential applications in protein-ligand and protein-protein docking and in biasing molecular dynamics simulations.     

Monte Carlo folding of trans-membrane helical peptides in an implicit generalized Born membrane

An efficient Monte Carlo (MC) algorithm using concerted backbone rotations is combined with a recently developed implicit membrane model to simulate the folding of the hydrophobic transmembrane domain M2TM of the M2 protein from influenza A virus and Sarcolipin at atomic resolution. The implicit membrane environment is based on generalized Born theory and has been calibrated against experimental data. The MC sampling has previously been used to fold several small polypeptides and been shown to be equivalent to molecular dynamics (MD). In combination with a replica exchange algorithm, M2TM is found to form continuous membrane spanning helical conformations for low temperature replicas. Sarcolipin is only partially helical, in agreement with the experimental NMR structures in lipid bilayers and detergent micelles. Higher temperature replicas exhibit a rapidly decreasing helicity, in agreement with expected thermodynamic behavior. To exclude the possibility of an erroneous helical bias in the simulations, the model is tested by sampling a synthetic Alanine-rich polypeptide of known helicity. The results demonstrate there is no overstabilization of helical conformations, indicating that the implicit model captures the essential components of the native membrane environment for M2TM and Sarcolipin.     

Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.     

Effects of cholesterol on the mechanism of fengycin,a biofungicide

Fengycins are a class of antifungal lipopeptides synthesized by the bacteria Bacillus subtilis, commercially available as the primary component of the agricultural fungicide Serenade. They are toxic to fungi but far less to mammalian cells. One key difference between mammalian and fungal cell membranes is the presence of cholesterol only in the former; recent experimental work showed that the presence of cholesterol reduces fengycin-induced membrane leakage. Since our previous all-atom and coarse-grained simulations suggested that aggregation of membrane-bound fengycin is central to its ability to disrupt membranes, we hypothesized that cholesterol might reduce fengycin aggregation. Here, we test this hypothesis using coarse-grained molecular dynamics simulations, with sampling enhanced via the weighted ensemble method. The results indicate that cholesterol subtly alters the size distribution for fengycin aggregates, limits the lateral range of their membrane disordering, and reduces the ability of aggregates to bend the membrane. Taken together, these phenomena may account for cholesterol’s effects on fengycin activity.     

Mode of insertion of chlorophyll a in a lipid layer. A conformational approach

A new procedure for conformational analysis is described to define the orientation of chlorophyll a in model membranes. The method allows to predict the configuration of assembled amphiphilic molecules. The minimal conformational energy is calculated as the sum of the contributions, resulting from the Van der Waals interactions, the torsional potentials, the electrostatic interactions and the transfer energy. In the calculated conformation, the porphyrin ring is orientated at an angle of 45° ± 5° to the membrane interface, with the phytol chain inserted into the lipid layer. This orientation of the porphyrin ring is in excellent agreement with the experimental value obtained with model membranes. The method could allow to define the orientation of other components in the thylakoid membrane.     

Molecular dynamics simulations of the interaction of phospholipid bilayers with polycaprolactone

The molecular interaction between common polymer chains and the cell membrane is unknown. Molecular dynamics simulations offer an emerging tool to characterise the nature of the interaction between common degradable polymer chains used in biomedical applications, such as polycaprolactone, and model cell membranes. Herein we characterise with all-atomistic and coarse-grained molecular dynamics simulations the interaction between single polycaprolactone chains of varying chain lengths with a phospholipid membrane. We find that the length of the polymer chain greatly affects the nature of interaction with the membrane, as well as the membrane properties. Furthermore, we next utilise advanced sampling techniques in molecular dynamics to characterise the two-dimensional free energy surface for the interaction of varying polymer chain lengths (short, intermediate, and long) with model cell membranes. We find that the free energy minimum shifts from the membrane-water interface to the hydrophobic core of the phospholipid membrane as a function of chain length. Finally, we perform coarse-grained molecular dynamics simulations of slightly larger membranes with polymers of the same length and characterise the results as compared with all-atomistic molecular dynamics simulations. These results can be used to design polymer chain lengths and chemistries to optimise their interaction with cell membranes at the molecular level.     

Probing the Periplasmic-Open State of Lactose Permease in Response to Sugar Binding and Proton Translocation

Based on the crystal structure of lactose permease (LacY) open to the cytoplasm, a hybrid molecular simulation approach with self-guided Langevin dynamics is used to describe conformational changes that lead to a periplasmic-open state. This hybrid approach consists of implicit (IM) and explicit (EX) membrane simulations and requires self-guided Langevin dynamics to enhance protein motions during the IM simulations. The pore radius of the lumen increases by 3.5 Å on the periplasmic side and decreases by 2.5 Å on the cytoplasmic side (relative to the crystal structure), suggesting a lumen that is fully open to the periplasm to allow for extracellular sugar transport and closed to the cytoplasm. Based on our simulations, the mechanism that triggers this conformational change to the periplasmic-open state is the protonation of Glu269 and binding of the disaccharide. Then, helix packing is destabilized by breaking of several side chains involved in hydrogen bonding (Asn245, Ser41, Glu374, Lys42, and Gln242). For the periplasmic-open conformations obtained from our simulations, helix-helix distances agree well with experimental measurements using double electron-electron resonance, fluorescence resonance energy transfer, and varying sized cross-linkers. The periplasmic-open conformations are also in compliance with various substrate accessibility/reactivity measurements that indicate an opening of the protein lumen on the periplasmic side on sugar binding. The comparison with these measurements suggests a possible incomplete closure of the cytoplasmic half in our simulations. However, the closure is sufficient to prevent the disaccharide from transporting to the cytoplasm, which is in accordance with the well-established alternating access model. Ser53, Gln60, and Phe354 are determined to be important in sugar transport during the periplasmic-open stage of the sugar transport cycle and the sugar is found to undergo an orientational change in order to escape the protein lumen.     

De novo folding of membrane proteins: an exploration of the structure and NMR properties of the fd coat protein

De novo folding simulations of the major pVIII coat protein from filamentous fd bacteriophage, using a newly developed implicit membrane generalized Born model and replica-exchange molecular dynamics, are presented and discussed. The quality of the predicted structures, judged by comparison of the root-mean-square deviations of a room temperature ensemble of conformations from the replica-exchange simulations and experimental structures from both solid-state NMR in lipid bilayers and solution-phase NMR on the protein in micelles, was quite good, reinforcing the general quality of the folding simulations. The transmembrane helical segment of the protein was well defined in comparison with experiment and the amphipathic helical fragment remained at the membrane/aqueous phase boundary while undergoing significant conformational flexibility due to the loop connecting the two helical segments of the protein. Additional comparisons of computed solid-state NMR properties, the 15N chemical shift and 15N-1H dipolar coupling constants, showed semi-quantitative agreement with the corresponding measurements. These findings suggest an emerging potential for the de novo investigation of integral membrane peptides and proteins and a mechanism to assist experimental approaches to the characterization and structure determination of these important systems.     

Simulations of a membrane-anchored peptide: structure, dynamics, and influence on bilayer properties

A three-dimensional structure of a model decapeptide is obtained by performing molecular dynamics simulations of the peptide in explicit water. Interactions between an N-myristoylated form of the folded peptide anchored to dipalmitoylphosphatidylcholine fluid phase lipid membranes are studied at different applied surface tensions by molecular dynamics simulations. The lipid membrane environment influences the conformational space explored by the peptide. The overall secondary structure of the anchored peptide is found to deviate at times from its structure in aqueous solution through reversible conformational transitions. The peptide is, despite the anchor, highly mobile at the membrane surface with the peptide motion along the bilayer normal being integrated into the collective modes of the membrane. Peptide anchoring moderately alters the lateral compressibility of the bilayer by changing the equilibrium area of the membrane. Although membrane anchoring moderately affects the elastic properties of the bilayer, the model peptide studied here exhibits conformational flexibility and our results therefore suggest that peptide acylation is a feasible way to reinforce peptide-membrane interactions whereby, e.g., the lifetime of receptor-ligand interactions can be prolonged.     

Flexibility of the Bacterial Chaperone Trigger Factor in Microsecond-Timescale Molecular Dynamics Simulations

The bacterial chaperone trigger factor (TF) is the first chaperone to be encountered by a nascent protein chain as it emerges from the ribosome exit tunnel. Experimental results suggest that TF possesses considerable conformational flexibility, and in an attempt to provide an atomic-level view of this flexibility, we have performed independent 1.5-μs molecular dynamics simulations of TF in explicit solvent using two different simulation force fields (OPLS-AA/L and AMBER ff99SB-ILDN). Both simulations indicate that TF possesses tremendous flexibility, with huge excursions from the crystallographic conformation caused by reorientations of the protein’s constituent domains; both simulations also predict the formation of extensive contacts between TF’s PPIase domain and the Arm 1 domain that is involved in nascent-chain binding. In the OPLS simulation, however, TF rapidly settles into a very compact conformation that persists for at least 1 μs, whereas in the AMBER simulation, it remains highly dynamic; additional simulations in which the two force fields were swapped suggest that these differences are at least partly attributable to sampling issues. The simulation results provide potential rationalizations of a number of experimental observations regarding TF’s conformational behavior and have implications for using simulations to model TF’s function on translating ribosomes.     

A multidomain outer membrane protein from Pasteurella multocida: Modelling and simulation studies of PmOmpA

PmOmpA is a two-domain outer membrane protein from Pasteurella multocida. The N-terminal domain of PmOmpA is a homologue of the transmembrane β-barrel domain of OmpA from Escherichia coli, whilst the C-terminal domain of PmOmpA is a homologue of the extra-membrane Neisseria meningitidis RmpM C-terminal domain. This enables a model of a complete two domain PmOmpA to be constructed and its conformational dynamics explored via MD simulations of the protein embedded within two different phospholipid bilayers (DMPC and DMPE). The conformational stability of the transmembrane β-barrel is similar to that of a homology model of OprF from Pseudomonas aeruginosa in bilayer simulations. There is a degree of water penetration into the interior of the β-barrel, suggestive of a possible transmembrane pore. Although the PmOmpA model is stable over 20 ns simulations, retaining its secondary structure and fold integrity throughout, substantial flexibility is observed in a short linker region between the N- and the C-terminal domains. At low ionic strength, the C-terminal domain moves to interact electrostatically with the lipid bilayer headgroups. This study demonstrates that computational approaches may be applied to more complex, multi-domain outer membrane proteins, rather than just to transmembrane β-barrels, opening the possibility of in silico proteomics approaches to such proteins.