CGDB: A database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

CGDB: a database of membrane protein/lipid interactions by coarse-grained molecular dynamics simulations

Membrane protein function and stability has been shown to be dependent on the lipid environment. Recently, we developed a high-throughput computational approach for the prediction of membrane protein/lipid interactions. In the current study, we enhanced this approach with the addition of a new measure of the distortion caused by membrane proteins on a lipid bilayer. This is illustrated by considering the effect of lipid tail length and headgroup charge on the distortion caused by the integral membrane proteins MscS and FLAP, and by the voltage sensing domain from the channel KvAP. Changing the chain length of lipids alters the extent but not the pattern of distortion caused by MscS and FLAP; lipid headgroups distort in order to interact with very similar but not identical regions in these proteins for all bilayer widths investigated. Introducing anionic lipids into a DPPC bilayer containing the KvAP voltage sensor does not affect the extent of bilayer distortion.     

The reconstitution and activity of the small multidrug transporter EmrE is modulated by non-bilayer lipid composition

The ability of multidrug transport proteins within biological membranes to recognise a diverse array of substrates is a fundamental aspect of antibiotic resistance. Detailed information on the mechanisms of recognition and transport can be provided only by in vitro studies in reconstituted bilayer systems. We describe the controlled, efficient reconstitution of the small multidrug transporter EmrE in a simple model membrane and investigate the effect of non-bilayer lipids on this process. Transport activity is impaired, in line with an increase in the lateral pressure within the bilayer. We demonstrate the potential of this lateral pressure modulation method as a general approach to the folding and assembly of membrane proteins in vitro, by recovering functional transporter from a partly denatured state. Our results highlight the importance of optimising reconstitution procedures and bilayer lipid composition in studies of membrane transporters. This is particularly pertinent for multidrug proteins, and we show that the use of a sub-optimal lipid bilayer environment or reconstitution method could lead to incorrect information on protein activity.     

Effect of aminazine and triftazin on the synaptosome membranes of the rat cerebral cortex

Fluorescent probes 1,6-diphenyl-1,3,5-hexatriene (DPH) and pyrene were employed in studying the effect of aminazine and triftazin versus that of imipramine on microviscosity of rat brain cortex synaptosomal membranes. Unlike imipramine, the neuroleptics decrease microviscosity of membrane's lipid bilayer. All drugs decrease fluorescence of endogenous tryptophan, but fail to change fluorescence of L-tryptophan in the solution. It is concluded that neuroleptics induce conformational perturbations in membrane-bound proteins modifying microviscosity of lipid bilayer whereas imipramine changes the surface electric charge of lipid bilayer of synaptosomal membranes.     

Coarse-grained MD simulations of membrane protein-bilayer self-assembly

Complete determination of a membrane protein structure requires knowledge of the protein position within the lipid bilayer. As the number of determined structures of membrane proteins increases so does the need for computational methods which predict their position in the lipid bilayer. Here we present a coarse-grained molecular dynamics approach to lipid bilayer self-assembly around membrane proteins. We demonstrate that this method can be used to predict accurately the protein position in the bilayer for membrane proteins with a range of different sizes and architectures.     

Cholesterol-induced protein sorting: an analysis of energetic feasibility

The mechanism(s) underlying the sorting of integral membrane proteins between the Golgi complex and the plasma membrane remain uncertain because no specific Golgi retention signal has been found. Moreover one can alter a protein's eventual localization simply by altering the length of its transmembrane domain (TMD). M. S. Bretscher and S. Munro (SCIENCE: 261:1280-1281, 1993) therefore proposed a physical sorting mechanism based on the hydrophobic match between the proteins' TMD and the bilayer thickness, in which cholesterol would regulate protein sorting by increasing the lipid bilayer thickness. In this model, Golgi proteins with short TMDs would be excluded from cholesterol-enriched domains (lipid rafts) that are incorporated into transport vesicles destined for the plasma membrane. Although attractive, this model remains unproven. We therefore evaluated the energetic feasibility of a cholesterol-dependent sorting process using the theory of elastic liquid crystal deformations. We show that the distribution of proteins between cholesterol-enriched and cholesterol-poor bilayer domains can be regulated by cholesterol-induced changes in the bilayer physical properties. Changes in bilayer thickness per se, however, have only a modest effect on sorting; the major effect arises because cholesterol changes also the bilayer material properties, which augments the energetic penalty for incorporating short TMDs into cholesterol-enriched domains. We conclude that cholesterol-induced changes in the bilayer physical properties allow for effective and accurate sorting which will be important generally for protein partitioning between different membrane domains.     

Molecular dynamics simulations of discoidal bilayers assembled from truncated human lipoproteins

Human apolipoprotein A-1 (apo A-1) is the major protein component of high-density lipoproteins. The apo A-1 lipid-binding domain was used as a template for the synthesis of amphipathic helical proteins termed membrane scaffold proteins, employed to self-assemble soluble monodisperse discoidal particles called Nanodiscs. In these particles, membrane scaffold proteins surround a lipid bilayer in a belt-like fashion forming bilayer disks of discrete size and composition. Here we investigate the structure of Nanodiscs through molecular dynamics simulations in which Nanodiscs were built from scaffold proteins of various lengths. The simulations showed planar or deformed Nanodiscs depending on optimal length and alignment of the scaffold proteins. Based on mean surface area per lipid calculations, comparison of small-angle x-ray scattering curves, and the relatively planar shape of Nanodiscs made from truncated scaffold proteins, one can conclude that the first 17 to 18 residues of the 200-residue apo A-1 lipid-binding domain are not involved in formation of the protein "belts" surrounding the lipid bilayer. To determine whether the addition of an integral membrane protein has an effect on the overall structure of a Nanodisc, bacteriorhodopsin was embedded into a Nanodisc and simulated using molecular dynamics, revealing a planar disk with a slightly rectangular shape.     

In situ molecular level studies on membrane related peptides and proteins in real time using sum frequency generation vibrational spectroscopy

Sum frequency generation (SFG) vibrational spectroscopy has been demonstrated to be a powerful technique to study the molecular structures of surfaces and interfaces in different chemical environments. This review summarizes recent SFG studies on hybrid bilayer membranes and substrate-supported lipid monolayers and bilayers, the interaction between peptides/proteins and lipid monolayers/bilayers, and bilayer perturbation induced by peptides/proteins. To demonstrate the ability of SFG to determine the orientations of various secondary structures, studies on the interactions between different peptides/proteins (melittin, G proteins, alamethicin, and tachyplesin I) and lipid bilayers are discussed. Molecular level details revealed by SFG in these studies show that SFG can provide a unique understanding on the interactions between a lipid monolayer/bilayer and peptides/proteins in real time, in situ and without any exogenous labeling.     

Lipid-protein interactions

Intrinsic membrane proteins are solvated by a shell of lipid molecules interacting with the membrane-penetrating surface of the protein; these lipid molecules are referred to as annular lipids. Lipid molecules are also found bound between transmembrane α-helices; these are referred to as non-annular lipids. Annular lipid binding constants depend on fatty acyl chain length, but the dependence is less than expected from models based on distortion of the lipid bilayer alone. This suggests that hydrophobic matching between a membrane protein and the surrounding lipid bilayer involves some distortion of the transmembrane α-helical bundle found in most membrane proteins, explaining the importance of bilayer thickness for membrane protein function. Annular lipid binding constants also depend on the structure of the polar headgroup region of the lipid, and hotspots for binding anionic lipids have been detected on some membrane proteins; binding of anionic lipid molecules to these hotspots can be functionally important. Binding of anionic lipids to non-annular sites on membrane proteins such as the potassium channel KcsA can also be important for function. It is argued that the packing preferences of the membrane-spanning α-helices in a membrane protein result in a structure that matches nicely with that of the surrounding lipid bilayer, so that lipid and protein can meet without either having to change very much.     

Effects of chlorpromazine HCl on the structural parameters of bovine brain membranes

Fluorescence probes located in different membrane regions were used to evaluate the effects of chlorpromazine .HCl on structural parameters (transbilayer lateral mobility, annular lipid fluidity, protein distribution, and lipid bilayer thickness) of synaptosomal plasma membrane vesicles (SPMVs) isolated from bovine cerebral cortex. The experimental procedure was based on the selective quenching of 1,3-di(1-pyrenyl)propane (Py-3-Py) by trinitrophenyl groups, radiationless energy transfer from the tryptophan of membrane proteins to Py-3-Py, and energy transfer from Py-3-Py monomers to 1-anilinonaphthalene-8-sulfonic acid (ANS). In this study, chlorpromazine .HCl decreased the lateral mobility of Py-3-Py in a concentration dependent-manner, showed a greater ordering effect on the inner monolayer than on the outer monolayer, decreased annular lipid fluidity in a dose dependent-manner, and contracted the membrane lipid bilayer. Furthermore, the drug was found to have a clustering effect on membrane proteins.     

The effect of propoxycaine.HCl on the physical properties of neuronal membranes

Fluorescent probe techniques were used to evaluate the effect of propoxycaine.HCl on the physical properties (transbilayer asymmetric lateral and rotational mobilities, annular lipid fluidity and protein distribution) of synaptosomal plasma membrane vesicles (SPMVs) isolated from bovine cerebral cortex. An experimental procedure was used based on selective quenching of both 1,3-di(1-pyrenyl)propane (Py-3-Py) and 1,6-diphenyl-1,3,5-hexatriene (DPH) by trinitrophenyl groups, and radiationless energy transfer (RET) from the tryptophans of membrane proteins to Py-3-Py. Propoxycaine.HCl increased the bulk lateral and rotational mobilities, and annular lipid fluidity in SPMVs lipid bilayers, and had a greater fluidizing effect on the inner monolayer than that of the outer monolayer. The magnitude of increasing effect on annular lipid fluidity in SPMVs lipid bilayer induced by propoxycaine.HCl was significantly far greater than magnitude of increasing effect of the drug on the lateral and rotational mobilities of SPMVs lipid bilayer. It also caused membrane proteins to cluster. These effects of propoxycaine.HCl on neuronal membranes may be responsible for some, though not all, of the local anesthetic actions of propoxycaine.HCl.     

The effect of bupivacaine·HCl on the physical properties of neuronal membranes

Fluorescent probe techniques were used to evaluate the effect of bupivacaine.HCl on the physical properties (transbilayer asymmetric lateral and rotational mobilities, annular lipid fluidity and protein distribution) of synaptosomal plasma membrane vesicles (SPMVs) isolated from bovine cerebral cortex. An experimental procedure was used based on selective quenching of both 1,3-di(1-pyrenyl)propane (Py-3-Py) and 1,6-diphenyl-1,3,5-hexatriene (DPH) by trinitrophenyl groups, and radiationless energy transfer (RET) from the tryptophans of membrane proteins to Py-3-Py. Bupivacaine.HCl increased the bulk lateral and rotational mobilities, and annular lipid fluidity in SPMVs lipid bilayers, and had a greater fluidizing effect on the inner monolayer than that of the outer monolayer. The magnitude of increasing effect on annular lipid fluidity in SPMVs lipid bilayer induced by bupivacaine.HCl was significantly far greater than magnitude of increasing effect of the drug on the lateral and rotational mobilities of bulk SPMVs lipid bilayer. It also caused membrane proteins to cluster. These effects of bupivacaine.HCl on neuronal membranes may be responsible for some, though not all, of the local anesthetic actions of bupivacaine.HCl.     

Manipulating the folding of membrane proteins: using the bilayer to our advantage

The folding mechanisms of integral membrane proteins have largely eluded detailed study. This is owing to the inherent difficulties in folding these hydrophobic proteins in vitro, which, in turn, reflects the often apparently insurmountable problem of mimicking the natural membrane bilayer with lipid or detergent mixtures. There is, however, a large body of information on lipid properties and, in particular, on phosphatidylcholine and phosphatidylethanolamine lipids, which are common to many biological membranes. We have exploited this knowledge to develop efficient in vitro lipid-bilayer folding systems for the membrane protein, bacteriorhodopsin. Furthermore, we have shown that a rate-limiting apoprotein folding step and the overall folding efficiency appear to be controlled by particular properties of the lipid bilayer. The properties of interest are the stored curvature elastic energy within the bilayer, and the lateral pressure that the lipid chains exert on the their neighbouring folding proteins. These are generic properties of the bilayer that can be achieved with simple mixtures of biological lipids, and are not specific to the lipids studied here. These bilayer properties also seem to be important in modulating the function of several membrane proteins, as well as the function of membranes in vivo. Thus, it seems likely that careful manipulations of lipid properties will shed light on the forces that drive membrane protein folding, and will aid the development of bilayer folding systems for other membrane proteins.     

High-resolution AFM of membrane proteins directly incorporated at high density in planar lipid bilayer

The heterologous expression and purification of membrane proteins represent major limitations for their functional and structural analysis. Here we describe a new method of incorporation of transmembrane proteins in planar lipid bilayer starting from 1 pmol of solubilized proteins. The principle relies on the direct incorporation of solubilized proteins into a preformed planar lipid bilayer destabilized by dodecyl-beta-maltoside or dodecyl-beta-thiomaltoside, two detergents widely used in membrane biochemistry. Successful incorporations are reported at 20 degrees C and at 4 degrees C with three bacterial photosynthetic multi-subunit membrane proteins. Height measurements by atomic force microscopy (AFM) of the extramembraneous domains protruding from the bilayer demonstrate that proteins are unidirectionally incorporated within the lipid bilayer through their more hydrophobic domains. Proteins are incorporated at high density into the bilayer and on incubation diffuse and segregate into protein close-packing areas. The high protein density allows high-resolution AFM topographs to be recorded and protein subunits organization delineated. This approach provides an alternative experimental platform to the classical methods of two-dimensional crystallization of membrane proteins for the structural analysis by AFM. Furthermore, the versatility and simplicity of the method are important intrinsic properties for the conception of biosensors and nanobiomaterials involving membrane proteins.     

Membrane lipids and transporter function

Transport proteins are essential for cells in allowing the exchange of substances between cells and their environment across the lipid bilayer forming a tight barrier. Membrane lipids modulate the function of transmembrane proteins such as transporters in two ways: Lipids are tightly and specifically bound to transport proteins and in addition they modulate from the bulk of the lipid bilayer the function of transport proteins. This overview summarizes currently available information at the ultrastructural level on lipids tightly bound to transport proteins and the impact of altered bulk membrane lipid composition. Human diseases leading to altered lipid homeostasis will lead to altered membrane lipid composition, which in turn affect the function of transporter proteins.     

Setting up and optimization of membrane protein simulations

In this paper we describe a method for setting up an atomistic simulation of a membrane protein in a hydrated lipid bilayer and report the effect of differing electrostatic parameters on the drift in the protein structure during the subsequent simulation. The method aims to generate a suitable cavity in the interior of a lipid bilayer, using the solvent-accessible surface of the protein as a template, during the course of a short steered molecular dynamics simulation of a solvated lipid membrane. This is achieved by a two-stage process: firstly, lipid molecules whose headgroups are inside a cylindrical volume equivalent to that defined by the protein surface are removed; then the protein-lipid interface is optimized by applying repulsive forces perpendicular to the protein surface, and of gradually increased magnitude, to the remaining lipid atoms inside the volume occupied by the protein surface until it is emptied. The protein itself may then be inserted. Using the bacterial membrane proteins KcsA and FhuA as test cases, we show how the method achieves the formation of a suitable cavity in the interior of a dimyristoylphosphatidylcholine lipid bilayer without perturbing the configuration of the non-interfacial regions of the previously equilibrated lipid bilayer, even in cases of membrane proteins with irregular geometrical shapes. In addition, we compare subsequent simulations in which the long-range electrostatic interactions are treated via either a cut-off or a particle-mesh Ewald method. The results show that the drift from the initial structure is less in the latter case, especially for KcsA and for the non-core secondary structural elements (i.e. surface loops) of both proteins.     

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.     

Insulin-induced changes in mechanical characteristics of lipid bilayers modified by liver plasma membrane fragments

Insulin interaction with BLM with incorporated fragments of rat liver plasma membranes, containing hormone receptors, was studied by determining Young modulus of elasticity of bilayer lipid membranes in direction perpendicular to the surface, E. The presence of membrane proteins in a concentration of 60 micrograms.ml-1 induced a significant decrease in parameter E (to approx. 50%) as compared with values obtained in non-modified membranes during insulin action (concentration interval 10(-11)-10(-9) mol.l-1). The extent of the effect was dependent on the initial phase state of the membrane, on cholesterol content in BLM as well as on membrane proteins concentration in lipid bilayer.     

Effect of diamide on protein oxidation and physico-chemical properties of lipids in erythrocyte membranes

The effect of diamide on the physicochemical state of proteins and lipids of human erythrocyte membrane was studied. It was found that diamide at a concentration of 1 mM decreases the content of the SH-groups of membrane proteins by approximately 50%, resulting in enhanced vesiculation of erythrocytes upon metabolic exhaustion of cells. It was shown using fluorescein isothiocyanate-labeled concanavalin A and 4,4'-diisothiocyano-2,2'-stilbene disulfonate that diamide changes the structural state of the main integral protein of erythrocyte membranes, the band 3 protein. Changes in the microviscosity of the membrane lipid bilayer depending on diamide concentration were determined from the changes in the fluorescence parameters of the lipophilic probes (pyrene and 1,6-diphenyl-3,5-hexatriene). The level of lipid peroxidation products in membranes remained unchanged. It follows from these data that the SH-oxidizing agent diamide does not directly interact with the lipid bilayer of membrane and produces changes in the physicochemical state of lipids presumably by disrupting protein-lipid interactions that take place upon oxidation of the SH-groups and cross-linking of membrane proteins.     

A general mechanism for drug promiscuity: Studies with amiodarone and other antiarrhythmics

Amiodarone is a widely prescribed antiarrhythmic drug used to treat the most prevalent type of arrhythmia, atrial fibrillation (AF). At therapeutic concentrations, amiodarone alters the function of many diverse membrane proteins, which results in complex therapeutic and toxicity profiles. Other antiarrhythmics, such as dronedarone, similarly alter the function of multiple membrane proteins, suggesting that a multipronged mechanism may be beneficial for treating AF, but raising questions about how these antiarrhythmics regulate a diverse range of membrane proteins at similar concentrations. One possible mechanism is that these molecules regulate membrane protein function by altering the common environment provided by the host lipid bilayer. We took advantage of the gramicidin (gA) channels’ sensitivity to changes in bilayer properties to determine whether commonly used antiarrhythmics—amiodarone, dronedarone, propranolol, and pindolol, whose pharmacological modes of action range from multi-target to specific—perturb lipid bilayer properties at therapeutic concentrations. Using a gA-based fluorescence assay, we found that amiodarone and dronedarone are potent bilayer modifiers at therapeutic concentrations; propranolol alters bilayer properties only at supratherapeutic concentration, and pindolol has little effect. Using single-channel electrophysiology, we found that amiodarone and dronedarone, but not propranolol or pindolol, increase bilayer elasticity. The overlap between therapeutic and bilayer-altering concentrations, which is observed also using plasma membrane–like lipid mixtures, underscores the need to explore the role of the bilayer in therapeutic as well as toxic effects of antiarrhythmic agents.