Alamethicin in lipid bilayers: Combined use of X-ray scattering and MD simulations

We study fully hydrated bilayers of two di-monounsaturated phospholipids diC18:1PC (DOPC) and diC22:1PC with varying amounts of alamethicin (Alm). We combine the use of X-ray diffuse scattering and molecular dynamics simulations to determine the orientation of alamethicin in model lipids. Comparison of the experimental and simulated form factors shows that Alm helices are inserted transmembrane at high humidity and high concentrations, in agreement with earlier results. The X-ray scattering data and the MD simulations agree that membrane thickness changes very little up to 1/10 Alm/DOPC. In contrast, the X-ray data indicate that the thicker diC22:1PC membrane thins with added Alm, a total decrease in thickness of 4 Å at 1/10 Alm/diC22:1PC. The different effect of Alm on the thickness changes of the two bilayers is consistent with Alm having a hydrophobic thickness close to the hydrophobic thickness of 27 Å for DOPC; Alm is then mismatched with the 7 Å thicker diC22:1PC bilayer. The X-ray data indicate that Alm decreases the bending modulus (K_C) by a factor of ∼ 2 in DOPC and a factor of ∼ 10 in diC22:1PC membranes (P/L ∼ 1/10). The van der Waals and fluctuational interactions between bilayers are also evaluated through determination of the anisotropic B compressibility modulus.     

Outer membrane protein A of Escherichia coli forms temperature-sensitive channels in planar lipid bilayers

The temperature dependence of single-channel conductance and open probability for outer membrane protein A (OmpA) of Escherichia coli were examined in planar lipid bilayers. OmpA formed two interconvertible conductance states, small channels, 36-140 pS, between 15 and 37 degrees C, and large channels, 115-373 pS, between 21 and 39 degrees C. Increasing temperatures had strong effects on open probabilities and on the ratio of large to small channels, particularly between 22 and 34 degrees C, which effected sharp increases in average conductance. The data infer that OmpA is a flexible temperature-sensitive protein that exists as a small pore structure at lower temperatures, but refolds into a large pore at higher temperatures.     

Interaction between dihydropyridines and phospholipid bilayers: a molecular dynamics simulation

Interaction of the calcium-channel antagonist dihydropyridines (DHPs), lacidipine and nifedipine, with a phospholipid bilayer was studied using 600 ps molecular dynamic simulations. We have constructed a double layer membrane model composed of 42 dimirystoyl-phosphatidylcholine molecules. The DHP molecules locate at about 7 Å from the centre of the membrane, inducing an asymmetry in the bilayer. While lacidipine did not induce significant local perturbations as judged by the gauche-trans isomerisation rate, nifedipine significantly decreased this rate, probably by producing a local rigidity of the membrane in the vicinity of the DHP.     

Lipid enrichment and selectivity of integral membrane proteins in two-component lipid bilayers

A model recently used to study lipid-protein interactions in one-component lipid bilayers (Sperotto and Mouritsen, 1991 a, b) has been extended in order to include two different lipid species characterized by different acyl-chain lengths. The model, which is a statistical mechanical lattice model, assumes that hydrophobic matching between lipid-bilayer hydrophobic thickness and hydrophobic length of the integral protein is an important aspect of the interactions. By means of Monte Carlo simulation techniques, the lateral distribution of the two lipid species near the hydrophobic protein-lipid interface in the fluid phase of the bilayer has been derived. The results indicate that there is a very structured and heterogeneous distribution of the two lipid species near the protein and that the protein-lipid interface is enriched in one of the lipid species. Out of equilibrium, the concentration profiles of the two lipid species away from the protein interface are found to develop a long-range oscillatory behavior. Such dynamic membrane heterogeneity may be of relevance for determining the physical factors involved in lipid specificity of protein function.     

Co-translational association of cell-free expressed membrane proteins with supplied lipid bilayers

Abstract

Routine strategies for the cell-free production of membrane proteins in the presence of detergent micelles and for their efficient co-translational solubilization have been developed. Alternatively, the expression in the presence of rationally designed lipid bilayers becomes interesting in particular for biochemical studies. The synthesized membrane proteins would be directed into a more native-like environment and cell-free expression of transporters, channels or other membrane proteins in the presence of supplied artificial membranes could allow their subsequent functional analysis without any exposure to detergents. In addition, lipid-dependent effects on activity and stability of membrane proteins could systematically be studied. However, in contrast to the generally efficient detergent solubilization, the successful stabilization of membrane proteins with artificial membranes appears to be more difficult. A number of strategies have therefore been explored in order to optimize the co-translational association of membrane proteins with different forms of supplied lipid bilayers including liposomes, bicelles, microsomes or nanodiscs. In this review, we have compiled the current state-of-the-art of this technology and we summarize parameters which have been indicated as important for the co-translational association of cell-free synthesized membrane proteins with supplied membranes.     

Outer membrane complex proteins of Chlamydia pneumoniae

Abstract The protein composition of the outer membrane complex (OMC) of Chlamydia pneumoniae strain AR-39 was analyzed by metabolic labeling with [³⁵S]methionine and [³⁵S]cysteine. Cysteine-rich proteins with molecular masses of 98, 60 doublet, 39.5 (MOMP) and 15.5 kDa were found in the OMC of C. pneumoniae . The cysteine-rich proteins of the OMCs of the threee Chlamydia species showed specific reaction patterns by immunoassay and autoradiography to rabbit or turkey immune sera. Recognition of the MOMP and 60-kDa proteins of the three species was cross-reactive. However, the C. pneumoniae 98-kDa protein was recognized by anti- C. pneumoniae (AR-39) and anti- C. psittaci (TT3) immune sera. None of the immunee sera recognized the 12-kDa cysteine-rich complex.     

Coarse-grained molecular dynamics simulations of membrane proteins and peptides

Molecular dynamics (MD) simulations provide a valuable approach to the dynamics, structure, and stability of membrane-protein systems. Coarse-grained (CG) models, in which small groups of atoms are treated as single particles, enable extended (>100 ns) timescales to be addressed. In this study, we explore how CG-MD methods that have been developed for detergents and lipids may be extended to membrane proteins. In particular, CG-MD simulations of a number of membrane peptides and proteins are used to characterize their interactions with lipid bilayers. CG-MD is used to simulate the insertion of synthetic model membrane peptides (WALPs and LS3) into a lipid (PC) bilayer. WALP peptides insert in a transmembrane orientation, whilst the LS3 peptide adopts an interfacial location, both in agreement with experimental biophysical data. This approach is extended to a transmembrane fragment of the Vpu protein from HIV-1, and to the coat protein from fd phage. Again, simulated protein/membrane interactions are in good agreement with solid state NMR data for these proteins. CG-MD has also been applied to an M3-M4 fragment from the CFTR protein. Simulations of CFTR M3-M4 in a detergent micelle reveal formation of an alpha-helical hairpin, consistent with a variety of biophysical data. In an I231D mutant, the M3-M4 hairpin is additionally stabilized via an inter-helix Q207/D231 interaction. Finally, CG-MD simulations are extended to a more complex membrane protein, the bacterial sugar transporter LacY. Comparison of a 200 ns CG-MD simulation of LacY in a DPPC bilayer with a 50 ns atomistic simulation of the same protein in a DMPC bilayer shows that the two methods yield comparable predictions of lipid-protein interactions. Taken together, these results demonstrate the utility of CG-MD simulations for studies of membrane/protein interactions.     

Characterization of the solubilization of lipid bilayers by surfactants

This communication addresses the state of aggregation of lipid-detergent mixed dispersions. Analysis of recently published data suggest that for any given detergent-lipid mixture the most important factor in determining the type of aggregates (mixed vesicles or mixed micelles) and the size of the aggregate is the detergent to lipid molar ratio in these aggregates, herein denoted the effective ratio, R_e. For mixed bilayers this effective ratio has been previously shown to be a function of the lipid and detergent concentrations and of an equilibrium partition coefficient, K, which describes the distribution of the detergent between the bilayers and the aqueous phase. We show that, similar to mixed bilayers, the size of mixed micelles is also a function of the effective ratio, but for these dispersions the distribution of detergent between the mixed micelles and the aqueous medium obeys a much higher partition coefficient. In practical terms, the detergent concentration in the mixed micelles is equal to the difference between the total detergent concentration and the critical micelle concentration (cmc). Thus, the effective ratio is equal to this difference divided by the lipid concentration. Transformation of mixed bilayers to mixed micelles, commonly denoted solubilization, occurs when the surfactant to lipid effective ratio reaches a critical value. Experimental evaluation of this critical ratio can be based on the linear dependence of detergent concentration, required for solubilization, on the lipid concentration. According to the ‘equilibrium partition model’, the dependence of the ‘solubilizing detergent concentration’ on the lipid concentration intersects with the lipid axis at −1/K, while the slope of this dependence is the critical effective ratio. On the other hand, assuming that when solubilization occurs the detergent concentration in the aqueous phase is approximately equal to the critical micelle concentration, implies that the above dependence intersects with the detergent axis at the critical micelle concentration, while its slope, again, is equal to the critical effective ratio. Analysis of existing data suggests that within experimental error both these distinctively different approaches are valid, indicating that the critical effective ratio at which solubilization occurs is approximately equal to the product of the critical micelle concentration and the distribution coefficient K. Since the nature of detergent affects K and the critical micelle concentration in opposite directions, the critical (‘solubilizing’) effective ratio depends upon the nature of detergent less than any of these two factors.     

Outer membrane proteins as major antigens of Fusobacterium nucleatum

Abstract The immunochemical reactions of rabbit polyclonal antibodies directed to different preparations of Fusobacterium nucleatum i.e, whole cells, peptidoglycan associated proteins, a peptidoglycan-protein complex and a purified 40 kiloDalton (kDa) protein, were investigated on outer membrane preparations of Fusobacterium species and a restricted number of Leptotrichia buccalis after their separation on sodium dodecyl sulphate polyacrylamide gels and electrotransfer to nitrocellulose. All F. nucleatum strains had identical reaction patterns with the immune sera tested. Surface exposed parts of a restricted number of proteins with apparent molecular weights at 70 kDa (a doublet band), 60 kDa, 55 kDa and 40 kDa seemed to be major immunogens. Antigenic related proteins either of identical or slightly deviating electrophoretic mobilities to the 40-kDa protein were observed with the other members of Bacteroidaceae tested. The characteristic 70-kDa protein doublet seemed to be restricted to F. nucleatum although single protein bands of near identical molecular weights belonging to the other species tested also reacted. The data also indicate that the 60-kDa and 55-kDa polypeptides might be present in other species of Fusobacterium .     

Effect of magnesium sulfate in oxidized lipid bilayers properties by using molecular dynamics

Magnesium sulfate (MgSO₄) has been used as a protector agent for many diseases related to oxidative stress. The effect of MgSO₄ on the oxidized lipid bilayer has not yet been studied using molecular dynamics calculations. In this work, the effects of oxidation were evaluated by using a POPC membrane model at different concentrations of its aldehyde (-CHO) and hydroperoxide (-OOH) derivatives with and without MgSO₄. Several quantitative and qualitative properties were evaluated, such as membrane thickness, area per lipid, area compressibility modulus, snapshots after simulation finish, density distributions, time evolutions of oxidized group positions, and radial distributions of oxidized group concerning Mg. Results indicate that in the absence of MgSO₄ the mobility of oxidized groups, particularly –CHO, toward the surface interface is high. At a low oxidation level of the bilayer there is an increase in the compressibility modulus as compared to the unoxidized bilayer. MgSO₄, at a low oxidation level, tends to lessen the oxidation effects by lowering the dispersion in the distribution of oxidized species toward the membrane surface and the water region. However, MgSO₄ does not change the trends of decreasing membrane thickness and area compressibility modulus and increasing area per lipid upon oxidation. In this regard, MgSO₄ diminishes the electrostatic long-distance attractive interactions between the oxidized groups and the charged headgroups of the interface, owing to the Mg⁺² and SO₄^-2 screening effects and an electrostatic stabilization of the headgroups, preventing the pore formation, which is well-known to occur in oxidized membranes.     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.     

Fluorescent membrane probes’ behavior in lipid bilayers: insights from molecular dynamics simulations

Fluorescence spectroscopy and microscopy have been used as tools to study membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Two major points of concern arise regarding this matter, namely the incomplete understanding of the probe behavior inside the bilayer and the perturbation of the latter resulting from probe incorporation. To this effect, molecular dynamics (MD) simulations, by providing detailed atomic-scale information, represent a valuable way to characterize the location and dynamics of bilayer-inserted membrane probes, as well as the magnitude of perturbation they induce on the host lipid structure, and several important classes of reporter molecules have been studied in recent years. This article reviews the state of the art of MD simulations of bilayer-inserted fluorescent probes, focusing on the information that has been obtained from previous studies and hinting at future perspectives in this rapidly emerging field.     

Expression of outer membrane proteins by Salmonella enteritidis relating to pH

Abstract The pH of the environment influenced the expression of outer membrane protein by S. enteritidis PT4 growing in broth. Growth in broth at pH 5 to 7 resulted in variation in expression of outer membrane proteins of 18 to 22 kDa. Bacteria became acid-fixed and non-viable following prolonged incubation in broth with a pH below 5, and expression of flagella was repressed.     

OmpA: Gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of ∼ 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.     

NMR of membrane proteins in micelles and bilayers: the FXYD family proteins

Determining the atomic resolution structures of membrane proteins is of particular interest in contemporary structural biology. Helical membrane proteins constitute one-third of the expressed proteins encoded in a genome, many drugs have membrane-bound proteins as their receptors, and mutations in membrane proteins result in human diseases. Although integral membrane proteins provide daunting technical challenges for all methods of protein structure determination, nuclear magnetic resonance (NMR) spectroscopy can be an extremely versatile and powerful method for determining their structures and characterizing their dynamics, in lipid environments that closely mimic the cell membranes. Once milligram amounts of isotopically labeled protein are expressed and purified, micelle samples can be prepared for solution NMR analysis, and lipid bilayer samples can be prepared for solid-state NMR analysis. The two approaches are complementary and can provide detailed structural and dynamic information. This paper describes the steps for membrane protein structure determination using solution and solid-state NMR. The methods for protein expression and purification, sample preparation and NMR experiments are described and illustrated with examples from the FXYD proteins, a family of regulatory subunits of the Na,K-ATPase.     

Characterization of domain instabilities in lipid bilayers by Karhunen-Loeve analysis

We study the stability of lipid bilayers with artificial domains. In investigating different domain structures, we identify scenarios of stable and unstable arrangements of patches of mixed phospholipids. These are then characterized using Karhunen-Loeve Expansion (KLE), a special form of Principal Components Analysis (PCA). The simulation data are interrogated using KLE to reveal spatiotemporal patterns that explain relevant motions in the bilayer system. By projecting the high-dimensional dataset onto a small number of key modes, KLE reveals specific dynamic signatures that can help distinguish and characterize various domain instability mechanisms. We find that typically very few modes are responsible for describing a mechanism of instability to a reasonable extent and can clearly distinguish between stable and unstable arrangements. Different instability modes are characterized as they exhibit unique features like global deformation or local mixing modes.     

Rhodopsin and other proteins in artificial lipid membranes

Some basic aspects of incorporation of hydrophobic peptides and proteins in artificial lipid membranes are discussed. As examples valinomycin as a carrier model and gramicidin A as a channel former in lipid vesicles and in planar lipid membranes are presented.In the second part of the lecture some examples of incorporation of membrane proteins into lipid vesicles and planar lipid membranes are reported. The interaction with artificial lipid membranes of the Ca⁺ ATPase from the sarcoplasmic reticulum, of Rhodopsin, and of Bacteriorhodopsin is presented.Presented at the EMBO-Workshop on Transduction Mechanism of Photoreceptors, Jülich, Germany, October 4–8, 1976     

Setting up and running molecular dynamics simulations of membrane proteins

Molecular dynamics simulations have become a popular and powerful technique to study lipids and membrane proteins. We present some general questions and issues that should be considered prior to embarking on molecular dynamics simulation studies of membrane proteins and review common simulation methods. We suggest a practical approach to setting up and running simulations of membrane proteins, and introduce two new (related) methods to embed a protein in a lipid bilayer. Both methods rely on placing lipids and the protein(s) on a widely spaced grid and then 'shrinking' the grid until the bilayer with the protein has the desired density, with lipids neatly packed around the protein. When starting from a grid based on a single lipid structure, or several potentially different lipid structures (method 1), the bilayer will start well-packed but requires more equilibration. When starting from a pre-equilibrated bilayer, either pure or mixed, most of the structure of the bilayer stays intact, reducing equilibration time (method 2). The main advantages of these methods are that they minimize equilibration time and can be almost completely automated, nearly eliminating one time consuming step in MD simulations of membrane proteins.     

Purification of squid rhodopsin and reassembly into lipid bilayers

Rhodopsin from squid photoreceptor membranes was solubilized in octyl glucoside and purified to a single band on SDS-polyacrylamide gels of $\text{M}{}_{\mathrm{<mtext>r</mtext>}}$ 46 000. Purified rhodopsin was recombined with phospholipids to form vesicles by detergent dialysis. Spectroscopic analysis of the rhodopsin-lipid vesicles showed that the interconversion between acid and basic metarhodopsin had a $\text{p}\text{K}$ of 8. Furthermore, rhodopsin in the vesicles could be photoregenerated from metarhodopsin in solutions of either neutral or alkaline pH. These two spectroscopic properties are comparable to those for rhodopsin in photoreceptor membranes. The results indicate that the native conformation of rhodopsin is preserved during purification and after recombination with phospholipids into vesicles. This preparation is, therefore, an active starting point for functional reconstitution studies.