Characterization of the outer membrane protein OprF of Pseudomonas aeruginosa in a lipopolysaccharide membrane by computer simulation

The N-terminal domain of outer membrane protein OprF of Pseudomonas aeruginosa forms a membrane spanning eight-stranded antiparallel beta-barrel domain that folds into a membrane channel with low conductance. The structure of this protein has been modeled after the crystal structure of the homologous protein OmpA of Escherichia coli. A number of molecular dynamics simulations have been carried out for the homology modeled structure of OprF in an explicit molecular model for the rough lipopolysaccharide (LPS) outer membrane of P. aeruginosa. The structural stability of the outer membrane model as a result of the strong electrostatic interactions compared with simple lipid bilayers is restricting both the conformational flexibility and the lateral diffusion of the porin in the membrane. Constricting side-chain interactions within the pore are similar to those found in reported simulations of the protein in a solvated lipid bilayer membrane. Because of the strong interactions between the loop regions of OprF and functional groups in the saccharide core of the LPS, the entrance to the channel from the extracellular space is widened compared with the lipid bilayer simulations in which the loops are extruding in the solvent. The specific electrostatic signature of the LPS membrane, which results in a net intrinsic dipole across the membrane, is found to be altered by the presence of OprF, resulting in a small electrically positive patch at the position of the channel.     

Regulation of channel function due to physical energetic coupling with a lipid bilayer

Regulation of membrane protein functions due to hydrophobic coupling with a lipid bilayer has been investigated. An energy formula describing interactions between lipid bilayer and integral ion channels with different structures, which is based on the screened Coulomb interaction approximation, has been developed. Here the interaction energy is represented as being due to charge-based interactions between channel and lipid bilayer. The hydrophobic bilayer thickness channel length mismatch is found to induce channel destabilization exponentially while negative lipid curvature linearly. Experimental parameters related to channel dynamics are consistent with theoretical predictions. To measure comparable energy parameters directly in the system and to elucidate the mechanism at an atomistic level we performed molecular dynamics (MD) simulations of the ion channel forming peptide–lipid complexes. MD simulations indicate that peptides and lipids experience electrostatic and van der Waals interactions for short period of time when found within each other’s proximity. The energies from these two interactions are found to be similar to the energies derived theoretically using the screened Coulomb and the van der Waals interactions between peptides (in ion channel) and lipids (in lipid bilayer) due to mainly their charge properties. The results of in silico MD studies taken together with experimental observable parameters and theoretical energetic predictions suggest that the peptides induce ion channels inside lipid membranes due to peptide–lipid physical interactions. This study provides a new insight helping better understand of the underlying mechanisms of membrane protein functions in cell membrane leading to important biological implications.     

Transmembrane helix uniformity examined by spectral mapping of torsion angles

The environment and unique balance of molecular forces within lipid bilayers has a profound impact upon the structure, dynamics, and function of membrane proteins. We describe the biophysical foundations for the remarkable uniformity of many transmembrane helices that result from the molecular interactions within lipid bilayers. In fact, the characteristic uniformity of transmembrane helices leads to unique spectroscopic opportunities allowing for phi,psi torsion angles to be mapped directly onto solid state nuclear magnetic resonance (NMR) PISEMA spectra. Results from spectral simulations, the solid state NMR-derived structure of the influenza A M2 proton channel transmembrane domain, and high-resolution crystal structures of 27 integral membrane proteins demonstrate that transmembrane helices tend to be more uniform than previously thought. The results are discussed through the definition of a preferred range of backbone varphi,psi torsion angles for transmembrane alpha helices and are presented with respect to improving biophysical characterizations of integral membrane proteins.     

Interactions of the pleckstrin homology domain with phosphatidylinositol phosphate and membranes: characterization via molecular dynamics simulations

The mechanism of interaction of pleckstrin homology (PH) domains with phosphatidylinositol 4,5-bisphosphate (PIP 2)-containing lipid bilayers remains uncertain. While crystallographic studies have emphasized PH-inositol 1,4,5-trisphosphate (IP 3) interactions, biophysical studies indicate a degree of less specific protein-bilayer interactions. We have used molecular dynamics simulations to characterize the interactions of the PH domain from phospholipase C-delta1 with IP 3 and with PIP 2, the latter in lipid bilayers and in detergent micelles. Simulations of the PH domain in water reveal a reduction in protein flexibility when IP 3 is bound. Simulations of the PH domain bound to PIP 2 in lipid bilayers indicate a tightening of ligand-protein interactions relative to the PH-IP 3 complex, alongside formation of H-bonds between PH side chains and lipid (PC) headgroups, and a degree of penetration of hydrophobic side chains into the core of the bilayer. Comparison with simulations of the PH-bound domain to a PC bilayer in the absence of PIP 2 suggests that the presence of PIP 2 increases the extent of PH-membrane interactions. Thus, comparative molecular dynamics simulations reveal how a PI-binding domain undergoes changes in conformational dynamics on binding to a PIP 2-containing membrane and how interactions additional to those with the PI headgroup are formed.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

跨膜蛋白CD3ε与磷脂酰肌醇二磷酸在脂筏域体系中蛋白质-脂质相互作用的粗粒化模拟研究

细胞膜局部区域可形成富含饱和脂质、胆固醇、鞘脂的脂筏域作为其信号转导调控平台。传统实验手段在研究脂筏及其功能时受到系统复杂度高及区域结构瞬时性强等制约。近年来,分子动力学模拟技术为细胞膜的组织原则提供了重要的理论支撑,从简单的单一组分模型到多组分系统转变,最终形成了越来越多的细胞膜仿真模型。其中,粗粒化模拟由于其简化模型,可大副拓展模拟体系的复杂程度与模拟时间,在细胞膜以及蛋白质-脂质相互作用相关研究中得到了广泛应用。本文采用MARTINI粗粒化力场模拟,构建了一种含有阴离子脂质磷脂酰肌醇二磷酸(phosphatidylinositol diphosphate, PIP2)的混合膜体系。模拟结果表明,该体系在适当温度及饱和度条件下,能自发分层形成脂筏域;膜厚度、膜组分分布、膜组分流动性等多种参数均表明,脂筏结构形成且符合其结构特征;少量PIP2添加不影响分层特性且PIP2对脂筏具有显著亲和性。此外,利用该模型以跨膜信号蛋白CD3ε为例研究了脂筏域体系中蛋白质-脂质相互作用。结果表明,PIP2-CD3ε胞内区相互作用可能是脂筏招募CD3ε的驱动力,且该过程可受钙离子调控。本工作体现了粗粒化模拟在仿真膜相关研究中的巨大优势及良好应用前景。     

How Does a Voltage Sensor Interact with a Lipid Bilayer? Simulations of a Potassium Channel Domain

The nature of voltage sensing by voltage-activated ion channels is a key problem in membrane protein structural biology. The way in which the voltage-sensor (VS) domain interacts with its membrane environment remains unclear. In particular, the known structures of Kv channels do not readily explain how a positively charged S4 helix is able to stably span a lipid bilayer. Extended (2 x 50 ns) molecular dynamics simulations of the high-resolution structure of the isolated VS domain from the archaebacterial potassium channel KvAP, embedded in zwitterionic and in anionic lipid bilayers, have been used to explore VS/lipid interactions at atomic resolution. The simulations reveal penetration of water into the center of the VS and bilayer. Furthermore, there is significant local deformation of the lipid bilayer by interactions between lipid phosphate groups and arginine side chains of S4. As a consequence of this, the electrostatic field is "focused" across the center of the bilayer.     

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.     

Investigating lipid composition effects on the mechanosensitive channel of large conductance (MscL) using molecular dynamics simulations

Previous experimental work has shown that the functional properties of the mechanosensitive channel of large conductance (MscL) are affected by variations in lipid composition. Here, we utilize molecular dynamics simulations of Mycobacterium tuberculosis MscL to investigate such lipid composition effects on a molecular level. In particular, two sets of simulations were performed. In the first, trajectories using lipids with different headgroups (phosphatidylcholine and phosphatidylethanolamine) were compared. Protein-lipid interactions were clearly altered by the headgroup changes, leading to conformational differences in the C-terminal region of M. tuberculosis MscL. In the second set of simulations, lipid tails were gradually shortened, thinning the membrane over a molecular dynamics trajectory. These simulations showed evidence of hydrophobic matching between MscL and the lipid membrane, as previously proposed. For all simulations, protein-lipid interaction energies in the second transmembrane region were correlated to mutagenic data, emphasizing the importance of lipid interactions for proper MscL function.     

Effects of PI(4,5)P2 concentration on the F-BAR domain membrane binding as revealed by coarse-grained simulations

Bin/Amphyphysin/Rvs (BAR) domain proteins form a key link between membrane remodeling and cytoskeleton dynamics. They are dimers that bind to membranes via electrostatic interactions with different preferences toward negatively charged lipids. In the present article, we examine the interactions of the F-BAR domain of nervous wreck (Nwk) with phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂)-containing membranes using coarse-grained molecular dynamics. We demonstrated PI(4,5)P₂ concentration effects, identified the sequence of events that underlies the protein binding and identified amino acids involved in protein–lipid interactions. Our simulations point out the primary role of the basic stretch at the tips of the dimer, which anchors the protein to the membrane and initiates the binding process. When the PI(4,5)P₂ concentration is high, the protein stably associates with the membrane by its concave surface or by the opposite side. At low PI(4,5)P₂ concentration, the former orientation becomes more favorable; also a state with only one tip bound is observed, due to the weaker attachment and more pronounced association/dissociation events. Our results provide a theoretical model that describes the lipid-binding behavior of Nwk observed in vitro.     

Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK.     

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.     

Pulling single bacteriorhodopsin out of a membrane: Comparison of simulation and experiment

Mechanical unfolding of single bacteriorhodopsins from a membrane bilayer is studied using molecular dynamics simulations. The initial conformation of the lipid membrane is determined through all-atom simulations and then its coarse-grained representation is used in the studies of stretching. A Go-like model with a realistic contact map and with Lennard-Jones contact interactions is applied to model the protein-membrane system. The model qualitatively reproduces the experimentally observed differences between force-extension patterns obtained on bacteriorhodopsin at different temperatures and predicts a lack of symmetry in the choice of the terminus to pull by. It also illustrates the decisive role of the interactions of the protein with the membrane in determining the force pattern and thus the stability of transmembrane proteins.     

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.     

Structure and dynamics of K channel pore-lining helices: a comparative simulation study

Isolated pore-lining helices derived from three types of K-channel have been analyzed in terms of their structural and dynamic features in nanosecond molecular dynamics (MD) simulations while spanning a lipid bilayer. The helices were 1) M1 and M2 from the bacterial channel KcsA (Streptomyces lividans), 2) S5 and S6 from the voltage-gated (Kv) channel Shaker (Drosophila melanogaster), and 3) M1 and M2 from the inward rectifier channel Kir6.2 (human). In the case of the Kv and Kir channels, for which x-ray structures are not known, both short and long models of each helix were considered. Each helix was incorporated into a lipid bilayer containing 127 palmitoyloleoylphosphatidylcholine molecules, which was solvated with approximately 4000 water molecules, yielding approximately 20, 000 atoms in each system. Nanosecond MD simulations were used to aid the definition of optimal lengths for the helix models from Kv and Kir. Thus the study corresponds to a total simulation time of 10 ns. The inner pore-lining helices (M2 in KcsA and Kir, S6 in Shaker) appear to be slightly more flexible than the outer pore-lining helices. In particular, the Pro-Val-Pro motif of S6 results in flexibility about a molecular hinge, as was suggested by previous in vacuo simulations (, Biopolymers. 39:503-515). Such flexibility may be related to gating in the corresponding intact channel protein molecules. Analysis of H-bonds revealed interactions with both water and lipid molecules in the water/bilayer interfacial region. Such H-bonding interactions may lock the helices in place in the bilayer during the folding of the channel protein (as is implicit in the two-stage model of membrane protein folding). Aromatic residues at the extremities of the helices underwent complex motions on both short (<10 ps) and long (>100 ps) time scales.     

Multiscale modeling of lipids and lipid bilayers

A multiscale modeling approach is applied for simulations of lipids and lipid assemblies on mesoscale. First, molecular dynamics simulation of initially disordered system of lipid molecules in water within all-atomic model was carried out. On the next stage, structural data obtained from the molecular dynamics (MD) simulation were used to build a coarse-grained (ten sites) lipid model, with effective interaction potentials computed by the inverse Monte Carlo method. Finally, several simulations of the coarse-grained model on longer length- and time-scale were performed, both within Monte Carlo and molecular dynamics simulations: a periodical sample of lipid molecules ordered in bilayer, a free sheet of such bilayer without periodic boundary conditions, formation of vesicle from a plain membrane, process of self-assembly of lipids randomly dispersed in volume. It was shown that the coarse-grained model, developed exclusively from all-atomic simulation data, reproduces well all the basic features of lipids in water solution.     

Molecular dynamics study of aquaporin-1 water channel in a lipid bilayer

The aquaporin-1 water channel was modeled in a palmitoyl-oleoyl-phosphatidyl-choline lipid bilayer, by means of molecular dynamics simulations. Interaction of the protein with the membrane and inter-monomer interactions were analyzed. Structural features of the channel important for its biological function, including the Asn-Pro-Ala (NPA) motifs, and the diffusion of water molecules into the channels, were investigated. Simulations revealed the formation of single file water inside the channels for certain relative positions of the NPA motifs.     

Reversible Effects of Peptide Concentration and Lipid Composition on H-Ras Lipid Anchor Clustering

Dynamic clusters of lipid-anchored Ras proteins are important for high-fidelity signal transduction in cells. The average size of Ras nanoclusters was reported to be independent of protein expression levels, and cholesterol depletion is commonly used to test the raft-preference of nanoclusters. However, whether protein concentration and membrane domain stability affect Ras clustering in a reversible manner is not well understood. We used coarse-grained molecular dynamics simulations to examine the reversibility of the effects of peptide and cholesterol concentrations as well as a lipid domain-perturbing nanoparticle (C₆₀) on the dynamics and stability of H-Ras lipid-anchor nanoclusters. By comparing results from these simulations with previous observations from the literature, we show that effects of peptide/cholesterol concentrations on the dynamics and stability of H-Ras peptide nanoclusters are reversible. Our results also suggest a correlation between the stabilities of lipid domains and Ras nanoclusters, which is supported by our finding that C₆₀ penetrates into the liquid-disordered domain of the bilayer, destabilizing lipid domains and thereby the stability of the nanoclusters.     

Lipid Fingerprints and Cofactor Dynamics of Light-Harvesting Complex II in Different Membranes

Plant light-harvesting complex II (LHCII) is the key antenna complex for plant photosynthesis. We present coarse-grained molecular dynamics simulations of monomeric and trimeric LHCII in a realistic thylakoid membrane environment based on the Martini force field. The coarse-grained protein model has been optimized with respect to atomistic reference simulations. Our simulations provide detailed insights in the thylakoid lipid fingerprint of LHCII which compares well with experimental data from membrane protein purification. Comparing the monomer and trimeric LHCII reveals a stabilizing effect of trimerization on the chromophores as well as the protein. Moreover, the average chromophore distance shortens in the trimer leading to stronger excitonic couplings. When changing the native thylakoid environment to a model membrane the protein flexibility remains constant, whereas the chromophore flexibility is reduced. Overall, the presented LHCII model lays the foundation to investigate the μs dynamics of this key antenna protein of plants.