Electrostatic binding of proteins to membranes. Theoretical predictions and experimental results with charybdotoxin and phospholipid vesicles

Optimal Dielectric Boundary for Binding Free Energy Estimates in the Implicit Solvent

Accuracy of binding free energy calculations utilizing implicit solvent models is critically affected by parameters of the underlying dielectric boundary, specifically, the atomic and water probe radii. Here, a multidimensional optimization pipeline is used to find optimal atomic radii, specifically for binding calculations in the implicit solvent. To reduce overfitting, the optimization target includes separate, weighted contributions from both binding and hydration free energies. The resulting five-parameter radii set, OPT_BIND5D, is evaluated against experiment for binding free energies of 20 host-guest (H-G) systems, unrelated to the types of structures used in the training. The resulting accuracy for this H-G test set (root mean square error of 2.03 kcal/mol, mean signed error of -0.13 kcal/mol, mean absolute error of 1.68 kcal/mol, and Pearson's correlation of r = 0.79 with the experimental values) is on par with what can be expected from the fixed charge explicit solvent models. Best agreement with the experiment is achieved when the implicit salt concentration is set equal or close to the experimental conditions.

Coarse-grained molecular dynamics simulation of cation distribution profiles on negatively charged lipid membranes during phase separation

Revealing the ion distributions on a charged lipid membrane in aqueous solution under the influence of long-range interactions is essential for understanding the origin of the stability of the bilayer structure and the interaction between biomembranes and various electrolytes. However, the ion distributions and their dynamics associated with the phase separation process of the lipid bilayer membrane are still unclear. We perform coarse-grained molecular dynamics simulations to reveal the Na⁺ and Cl^- distributions on charged phospholipid bilayer membranes during phase separation. During the phase separation, cations closely follow the position of negatively charged lipids on a microsecond timescale and are rapidly redistributed parallel to the lipid bilayer. In the homogenous mixture of zwitterionic and negatively charged lipids, cations weakly follow negatively charged lipids, indicating the strong interaction between cations and negatively charged lipids. We also compare cation concentrations as a function of surface charge density obtained by our simulation with those obtained by a modified Poisson-Boltzmann theory. Including the ion finite size makes the statistical results consistent, suggesting the importance of the ion-ion interactions in aqueous solution. Our simulation results advance our understanding of ion distribution during phase separation.

Lipid–Protein Interactions in Plasma Membrane Organization and Function

Lipid–protein interactions in cells are involved in various biological processes, including metabolism, trafficking, signaling, host–pathogen interactions, and transmembrane transport. At the plasma membrane, lipid–protein interactions play major roles in membrane organization and function. Several membrane proteins have motifs for specific lipid binding, which modulate protein conformation and consequent function. In addition to such specific lipid–protein interactions, protein function can be regulated by the dynamic, collective behavior of lipids in membranes. Emerging analytical, biochemical, and computational technologies allow us to study the influence of specific lipid–protein interactions, as well as the collective behavior of membranes on protein function. In this article, we review the recent literature on lipid–protein interactions with a specific focus on the current state-of-the-art technologies that enable novel insights into these interactions.

Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Phosphatidylinositol synthesis, its selective salvage, and inter-regulation of anionic phospholipids in Toxoplasma gondii

Phosphatidylinositol (PtdIns) serves as an integral component of eukaryotic membranes; however, its biosynthesis in apicomplexan parasites remains poorly understood. Here we show that Toxoplasma gondii-a common intracellular pathogen of humans and animals-can import and co-utilize myo-inositol with the endogenous CDP-diacylglycerol to synthesize PtdIns. Equally, the parasite harbors a functional PtdIns synthase (PIS) containing a catalytically-vital CDP-diacylglycerol phosphotransferase motif in the Golgi apparatus. Auxin-induced depletion of PIS abrogated the lytic cycle of T. gondii in human cells due to defects in cell division, gliding motility, invasion, and egress. Isotope labeling of the PIS mutant in conjunction with lipidomics demonstrated de novo synthesis of specific PtdIns species, while revealing the salvage of other lipid species from the host cell. Not least, the mutant showed decline in phosphatidylthreonine, and elevation of selected phosphatidylserine and phosphatidylglycerol species, indicating a rerouting of CDP-diacylglycerol and homeostatic inter-regulation of anionic phospholipids upon knockdown of PIS. In conclusion, strategic allocation of own and host-derived PtdIns species to gratify its metabolic demand features as a notable adaptive trait of T. gondii. Conceivably, the dependence of T. gondii on de novo lipid synthesis and scavenging can be exploited to develop new anti-infectives.

Multivalent assembly of KRAS with the RAS-binding and cysteine-rich domains of CRAF on the membrane

Membrane anchoring of farnesylated KRAS is critical for activation of RAF kinases, yet our understanding of how these proteins interact on the membrane is limited to isolated domains. The RAS-binding domain (RBD) and cysteine-rich domain (CRD) of RAF engage KRAS and the plasma membrane, unleashing the kinase domain from autoinhibition. Due to experimental challenges, structural insight into this tripartite KRAS:RBD-CRD:membrane complex has relied on molecular dynamics simulations. Here, we report NMR studies of the KRAS:CRAF RBD-CRD complex. We found that the nucleotide-dependent KRAS-RBD interaction results in transient electrostatic interactions between KRAS and CRD, and we mapped the membrane interfaces of the CRD, RBD-CRD, and the KRAS:RBD-CRD complex. RBD-CRD exhibits dynamic interactions with the membrane through the canonical CRD lipid-binding site (CRD β7-8), as well as an alternative interface comprising β6 and the C terminus of CRD and β2 of RBD. Upon complex formation with KRAS, two distinct states were observed by NMR: State A was stabilized by membrane association of CRD β7-8 and KRAS α4-α5 while state B involved the C terminus of CRD, β3-5 of RBD, and part of KRAS α5. Notably, α4-α5, which has been proposed to mediate KRAS dimerization, is accessible only in state B. A cancer-associated mutation on the state B membrane interface of CRAF RBD (E125K) stabilized state B and enhanced kinase activity and cellular MAPK signaling. These studies revealed a dynamic picture of the assembly of the KRAS-CRAF complex via multivalent and dynamic interactions between KRAS, CRAF RBD-CRD, and the membrane.

Influence of phenothiazine molecules on the interactions between positively charged poly-l-lysine and negatively charged DPPC/DPPG membranes

Phenothiazines are very effective antipsychotic drugs, which also have anticancer and antimicrobial activities. Despite being used in human treatment, the molecular mechanism of the biological actions of these molecules is not yet understood in detail. The role of the interactions between phenothiazines and proteins or lipid membranes has been much discussed. Herein, fourier-transform infrared (FTIR) spectroscopic studies were used to investigate the effect of three phenothiazines: fluphenazine (FPh); chlorpromazine (ChP); and propionylpromazine (PP) on the structures of a positively charged poly-l-lysine (PLL) peptide, a negatively charged dipalmitoylphosphatidylcholine/dipalmitoylphosphatidylglycerol (DPPC/DPPG) membrane, and on the mutual interactions between electrostatically associated PLL molecules and DPPC/DPPG membranes. Phenothiazine-induced alterations in the secondary structure of PLL, the conformational state (trans/gauche) of the hydrocarbon lipid chains, and the hydration of the DPPC/DPPG membrane interface were studied on the basis of amide I' vibrations, antisymmetric and symmetric stretching vibrations of the CH₂ groups of the lipid hydrocarbon chains (ν_sCH₂), and stretching vibrations of the lipid C=O groups (νC = O), respectively. It was shown that in the presence of negatively charged DPPC/DPPG membranes, the phenothiazines were able to modify the secondary structure of charged PLL molecules. Additionally, the effect of PLL on the structure of DPPC/DPPG membranes was also altered by the presence of the phenothiazine molecules.

Coarse-grained molecular dynamics simulation for uptake of nanoparticles into a charged lipid vesicle dominated by electrostatic interactions

We use a coarse-grained molecular dynamics simulation to investigate the interaction between neutral or charged nanoparticles (NPs) and a vesicle consisting of neutral and negatively charged lipids. We focus on the interaction strengths of hydrophilic and hydrophobic attraction and electrostatic interactions between a lipid molecule and an NP. A neutral NP passes through the lipid membrane when the hydrophobic interaction is sufficiently strong. As the valence of the positively charged NP increases, the membrane permeation speed of the NP is increased compared with the neutral NP and charged lipids are accumulated around the charged NP. A charged NP with a high valence passes through the lipid membrane via a transient channel formed by charged lipids or transportlike endocytosis. These permeation processes can be classified based on analyses of the density correlation function. When the nonelectrostatic interaction parameters are large enough, a negatively charged NP can be adsorbed on the membrane and a neutral lipid-rich region is formed directly below the NP. The NP is spontaneously incorporated into the vesicle under various conditions and the incorporation is mediated by the membrane curvature. We reveal how the NP's behavior depends on the NP valence, size, and the nonelectrostatic interaction parameters.

Influence of resveratrol on interactions between negatively charged DPPC/DPPG membranes and positively charged poly-l-lysine

Resveratrol (Res), a natural polyphenol present in different plants and vegetables, exhibits potential therapeutic activity with cardioprotective, antineurodegenerative, antioxidant, and antitumor action. In this study, the effect of Res on the mutual interactions between positively charged poly-l-lysine (PLL) and negatively charged dipalmitoylphosphatidylcholine/dipalmitoylphosphatidylglycerol (DPPC/DPPG) membranes was studied using Fourier-transform infrared (FTIR) spectroscopy supported by principal component analysis (PCA). The interactions between PLL and DPPC/DPPG membranes were strongly affected by the presence of Res molecules. Depending on the Res concentration and method of its supply (through the water or lipid phase) to the studied peptide-membrane systems, the membrane-induced transition of PLL from an α-helix to an extended left-handed polyproline II helix (PPII) occurred at different temperatures, with different cooperativity, or was even completely suppressed. The influence of PLL on the conformational (trans/gauche) state of the hydrocarbon chain region of the lipid membranes and the hydration state of the polar/apolar membrane interface was also modulated by Res, depending on the membrane phase state.

The cytosolic domain of T-cell receptor ζ associates with membranes in a dynamic equilibrium and deeply penetrates the bilayer

Interactions between lipid bilayers and the membrane-proximal regions of membrane-associated proteins play important roles in regulating membrane protein structure and function. The T-cell antigen receptor is an assembly of eight single-pass membrane-spanning subunits on the surface of T lymphocytes that initiates cytosolic signaling cascades upon binding antigens presented by MHC-family proteins on antigen-presenting cells. Its ζ-subunit contains multiple cytosolic immunoreceptor tyrosine-based activation motifs involved in signal transduction, and this subunit by itself is sufficient to couple extracellular stimuli to intracellular signaling events. Interactions of the cytosolic domain of ζ (ζcyt) with acidic lipids have been implicated in the initiation and regulation of transmembrane signaling. ζcyt is unstructured in solution. Interaction with acidic phospholipids induces structure, but its disposition when bound to lipid bilayers is controversial. Here, using surface plasmon resonance and neutron reflection, we characterized the interaction of ζcyt with planar lipid bilayers containing mixtures of acidic and neutral lipids. We observed two binding modes of ζcyt to the bilayers in dynamic equilibrium: one in which ζcyt is peripherally associated with lipid headgroups and one in which it penetrates deeply into the bilayer. Such an equilibrium between the peripherally bound and embedded forms of ζcyt apparently controls accessibility of the immunoreceptor tyrosine-based activation signal transduction pathway. Our results reconcile conflicting findings of the ζ structure reported in previous studies and provide a framework for understanding how lipid interactions regulate motifs to tyrosine kinases and may regulate the T-cell antigen receptor biological activities for this cell-surface receptor system.

Tyrosine Hydroxylase Binding to Phospholipid Membranes Prompts Its Amyloid Aggregation and Compromises Bilayer Integrity

Tyrosine hydroxylase (TH), a rate-limiting enzyme in the synthesis of catecholamine neurotransmitters and hormones, binds to negatively charged phospholipid membranes. Binding to both large and giant unilamellar vesicles causes membrane permeabilization, as observed by efflux and influx of fluorescence dyes. Whereas the initial protein-membrane interaction involves the N-terminal tail that constitutes an extension of the regulatory ACT-domain, prolonged membrane binding induces misfolding and self-oligomerization of TH over time as shown by circular dichroism and Thioflavin T fluorescence. The gradual amyloid-like aggregation likely occurs through cross-β interactions involving aggregation-prone motives in the catalytic domains, consistent with the formation of chain and ring-like protofilaments observed by atomic force microscopy in monolayer-bound TH. PC12 cells treated with the neurotoxin 6-hydroxydopamine displayed increased TH levels in the mitochondrial fraction, while incubation of isolated mitochondria with TH led to a decrease in the mitochondrial membrane potential. Furthermore, cell-substrate impedance and viability assays showed that supplementing the culture media with TH compromises cell viability over time. Our results revealed that the disruptive effect of TH on cell membranes may be a cytotoxic and pathogenic factor if the regulation and intracellular stability of TH is compromised.

Nutrient-Dependent Efficacy of the Antifungal Protein YvgO Correlates to Cellular Proliferation Rate in Candida albicans 3153A and Byssochlamys fulva H25

YvgO is a recently characterized antifungal protein isolated from Bacillus thuringiensis SF361 that exhibits a broad spectrum of activity and pH stability. Customized colorimetric metabolic assays based on standard broth microdilution techniques were used to determine the variable tolerance of Byssochlamys fulva H25 and Candida albicans 3153A to YvgO exposure under select matrix conditions impacting cellular proliferation. Normalization of the solution pH after antifungal challenge expanded the available pH range under consideration allowing for a comprehensive in vitro assessment of YvgO efficacy. Indicator susceptibility was examined across an array of elementary growth-modifying conditions, including media pH, incubation temperature, ionic strength, and carbohydrate supplementation. Under suboptimal temperature and pH conditions, the indicator growth rate reduced, and YvgO-mediated susceptibility was attenuated. While YvgO association but not efficacy was somewhat influenced by solution ionic strength, carbohydrate supplementation was shown to be the most influential susceptibility factor, particularly for C. albicans. Although the specific choice of carbohydrate/nutrient supplement dictated the extent of enhanced YvgO efficacy, D-glucose additionally improved the association between antifungal and target. Indeed, when exposed to YvgO under conditions that lead to increased cellular proliferation, both indicators displayed a stronger association and susceptibility to YvgO when compared to carbohydrate-deprived media or suboptimal incubation environments. With further study, YvgO may have the capacity to function as a prophylaxis for food safety and preservation, as well as a pharmaceutical agent against opportunistic fungal pathogens either independently or in combination with other established treatments applied to both livestock and human health concerns.

Method for measuring the unbinding energy of strongly-bound membrane-associated proteins

We describe a new method to measure the activation energy for unbinding (enthalpy ΔH*_u and free energy ΔG*_u) of a strongly-bound membrane-associated protein from a lipid membrane. It is based on measuring the rate of release of a liposome-bound protein during centrifugation on a sucrose gradient as a function of time and temperature. The method is used to determine ΔH*_u and ΔG*_u for the soluble dengue virus envelope protein (sE) strongly bound to 80:20 POPC:POPG liposomes at pH5.5. ΔH*_u is determined from the Arrhenius equation whereas ΔG*_u is determined by fitting the data to a model based on mean first passage time for escape from a potential well. The binding free energy ΔG_b of sE was also measured at the same pH for the initial, predominantly reversible, phase of binding to a 70:30 PC:PG lipid bilayer. The unbinding free energy (20±3kcal/mol, 20% PG) was found to be roughly three times the binding energy per monomer, (7.8±0.3kcal/mol for 30% PG, or est. 7.0kcal/mol for 20% PG). This is consistent with data showing that free sE is a monomer at pH5.5, but assembles into trimers after associating with membranes. This new method to determine unbinding energies should be useful to understand better the complex interactions of integral monotopic proteins and strongly-bound peripheral membrane proteins with lipid membranes.

Structural Thermodynamics of myr-Src(2-19) Binding to Phospholipid Membranes

Many proteins are anchored to lipid bilayer membranes through a combination of hydrophobic and electrostatic interactions. In the case of the membrane-bound nonreceptor tyrosine kinase Src from Rous sarcoma virus, these interactions are mediated by an N-terminal myristoyl chain and an adjacent cluster of six basic amino-acid residues, respectively. In contrast with the acyl modifications of other lipid-anchored proteins, the myristoyl chain of Src does not match the host lipid bilayer in terms of chain conformation and dynamics, which is attributed to a tradeoff between hydrophobic burial of the myristoyl chain and repulsion of the peptidic moiety from the phospholipid headgroup region. Here, we combine thermodynamic information obtained from isothermal titration calorimetry with structural data derived from (2)H, (13)C, and (31)P solid-state nuclear magnetic resonance spectroscopy to decipher the hydrophobic and electrostatic contributions governing the interactions of a myristoylated Src peptide with zwitterionic and anionic membranes made from lauroyl (C12:0) or myristoyl (C14:0) lipids. Although the latter are expected to enable better hydrophobic matching, the Src peptide partitions more avidly into the shorter-chain lipid analog because this does not require the myristoyl chain to stretch extensively to avoid unfavorable peptide/headgroup interactions. Moreover, we find that Coulombic and intrinsic contributions to membrane binding are not additive, because the presence of anionic lipids enhances membrane binding more strongly than would be expected on the basis of simple Coulombic attraction.

Regulation of membrane-shape transitions induced by I-BAR domains

I-BAR proteins are well-known actin-cytoskeleton adaptors and have been observed to be involved in the formation of plasma membrane protrusions (filopodia). I-BAR proteins contain an all-helical, crescent-shaped IRSp53-MIM domain (IMD) dimer that is believed to be able to couple with a membrane shape. This coupling could involve the sensing and even the generation of negative plasma membrane curvature. Indeed, the in vitro studies have shown that IMDs can induce inward tubulation of liposomes. While N-BAR domains, which generate positive membrane curvature, have received a considerable amount of attention from both theory and experiments, the mechanisms of curvature coupling through IMDs are comparatively less studied and understood. Here we used a membrane-shape stability assay developed recently in our lab to quantitatively characterize IMD-induced membrane-shape transitions. We determined a membrane-shape stability diagram for IMDs that reveals how membrane tension and protein density can comodulate the generation of IMD-induced membrane protrusions. From comparison to analytical theory, we determine three key parameters that characterize the curvature coupling of IMD. We find that the curvature generation capacity of IMDs is significantly stronger compared to that of endophilin, an N-BAR protein known to be involved in plasma membrane shape transitions. Contrary to N-BAR domains, where amphipathic helix insertion is known to promote its membrane curvature generation, for IMDs we find that amphipathic helices inhibit membrane shape transitions, consistent with the inverse curvature that IMDs generate. Importantly, in both of these types of BAR domains, electrostatic interactions affect membrane-binding capacity, but do not appear to affect the curvature generation capacity of the protein. These two types of BAR domain proteins show qualitatively similar membrane shape stability diagrams, suggesting an underlying ubiquitous mechanism by which peripheral proteins regulate membrane curvature.

Characterization of the binding between a 70-kDa heat shock protein, HspA1A, and phosphoinositides

HspA1A, a seventy-kilodalton heat shock protein, binds to specific anionic lipids and this interaction regulates important physiological phenomena like apoptosis, tumor growth, and lysosomal rescue. However, whether HspA1A binds to phosphoinositides has yet to be established and quantified. Therefore, in this study, we determined the binding affinity of HspA1A to several phosphoinositides and characterized five aspects of their molecular interaction. First, we established that HspA1A binds phosphatidylinositol monophosphates with higher affinity than di- and triphosphorylated inositides. Second, using high concentrations of potassium we found that HSPA1A embeds within the lipid bilayer of all phosphoinositides tested. However, the effects of the high salt concentrations were significantly different between the different phosphoinositides. Third, using calcium and reaction buffers equilibrated at different pH values we found that these differentially affected HspA1A-phosphoinositide binding, revealing a lipid-specific pattern of binding. Fourth, by assessing the binding properties of the two HspA1A domains, the nucleotide-binding domain and the substrate-binding domain, we determined that in most cases the full-length protein is necessary for binding to phosphoinositides. Fifth, by including in the reactions nucleotides and protein substrates we determined that they minimally and differentially affected phosphoinositide-binding. Collectively, these findings strongly suggest that the HspA1A-phosphoinositide binding is complex yet specific, is mediated by both electrostatic and hydrophobic interactions, is not related to the lipid-head charge, and depends on the physicochemical properties of the lipid.

Biochemical characterization of the interaction between HspA1A and phospholipids

Seventy-kilodalton heat shock proteins (Hsp70s) are molecular chaperones essential for maintaining cellular homeostasis. Apart from their indispensable roles in protein homeostasis, specific Hsp70s localize at the plasma membrane and bind to specific lipids. The interaction of Hsp70s with lipids has direct physiological outcomes including lysosomal rescue, microautophagy, and promotion of cell apoptosis. Despite these essential functions, the Hsp70-lipid interactions remain largely uncharacterized. In this study, we characterized the interaction of HspA1A, an inducible Hsp70, with five phospholipids. We first used high concentrations of potassium and established that HspA1A embeds in membranes when bound to all anionic lipids tested. Furthermore, we found that protein insertion is enhanced by increasing the saturation level of the lipids. Next, we determined that the nucleotide-binding domain (NBD) of the protein binds to lipids quantitatively more than the substrate-binding domain (SBD). However, for all lipids tested, the full-length protein is necessary for embedding. We also used calcium and reaction buffers equilibrated at different pH values and determined that electrostatic interactions alone may not fully explain the association of HspA1A with lipids. We then determined that lipid binding is inhibited by nucleotide-binding, but it is unaffected by protein-substrate binding. These results suggest that the HspA1A lipid-association is specific, depends on the physicochemical properties of the lipid, and is mediated by multiple molecular forces. These mechanistic details of the Hsp70-lipid interactions establish a framework of possible physiological functions as they relate to chaperone regulation and localization.

Free Energy Calculations for the Peripheral Binding of Proteins/Peptides to an Anionic Membrane. 1. Implicit Membrane Models

The binding of peptides and proteins to the surface of complex lipid membranes is important in many biological processes such as cell signaling and membrane remodeling. Computational studies can aid experiments by identifying physical interactions and structural motifs that determine the binding affinity and specificity. However, previous studies focused on either qualitative behaviors of protein/membrane interactions or the binding affinity of small peptides. Motivated by this observation, we set out to develop computational protocols for bimolecular binding to charged membranes that are applicable to both peptides and large proteins. In this work, we explore a method based on an implicit membrane/solvent model (generalized Born with a simple switching in combination with the Gouy-Chapman-Stern model for a charged interface), which we expect to lead to useful results when the binding does not implicate significant membrane deformation and local demixing of lipids. We show that the binding free energy can be efficiently computed following a thermodynamic cycle similar to protein-ligand binding calculations, especially when a Bennett acceptance ratio based protocol is used to consider both the membrane bound and solution conformational ensembles. Test calculations on a series of peptides show that our computational approach leads to binding affinities in encouraging agreement with experimental data, including for the challenging example of the bringing of flexible MARCKS-ED peptides to membranes. The calculations highlight that for a membrane with a significant fraction of anionic lipids, it is essential to include the effect of ion adsorption using the Stern model, which significantly modifies the effective surface charge. This implicit membrane model based computational protocol helps lay the groundwork for more systematic analysis of protein/peptide binding to membranes of complex shape and composition.

Octreotide inhibits the bilirubin carriers organic anion transporting polypeptides 1B1 and 1B3 and the multidrug resistance-associated protein 2

The somatostatin analog octreotide can lead to hyperbilirubinemia without evidence of liver injury. Here we investigate whether octreotide inhibits the main sinusoidal/canalicular bilirubin carriers and whether it is a transport substrate. Octreotide showed the most potent inhibitory effect toward OATP1B1-mediated transport and weaker inhibition for OATP1B3- and MRP2-mediated transport. Octreotide had no effect on OATP2B1-mediated transport. Octreotide inhibited [(3)H]estradiol-17-β-glucuronide (E17βG) influx mediated by OATP1B1, 1B3, and multidrug resistance-associated protein 2 (MRP2) in a concentration-dependent manner, and the IC50 values were computed to be 23 μM (95% confidence interval [CI] 18-29), 68 μM (95% CI 50-91), and 116.6 μM (95% CI 74.5-182.4), respectively. The interaction between octreotide and OATP1B1 was further studied. Inhibition of [(3)H]E17βG OATP1B1-mediated transport was purely competitive with no changes in maximum transport capacity (Vmax) and a twofold Km increase when the influx kinetics of [(3)H]E17βG were measured in the presence of octreotide (8.8 ± 3.1 versus 4.4 ± 1.2 μM, P = 0.03). The inhibition constant (Ki) of octreotide for the transport of [(3)H]E17βG was calculated at 33.5 ± 5.5 μM. Uptake of radiolabeled octreotide by OATP1B1-CHO cells was higher than in wild-type CHO cells and nonlabeled octreotide at the extracellular compartment was able to trans-stimulate the OATP1B1-mediated efflux of intracellular [(3)H]E17βG, suggesting that octreotide is a substrate of OATP1B1. In summary, this study shows interaction of octreotide on the human hepatocellular bilirubin transporters OATP1B1, OATP1B3, and MRP2, notably OATP1B1. These findings are in line with the clinical observation that a fraction of patients under treatment with octreotide exhibit hyperbilirubinemia.

Tarantula toxins use common surfaces for interacting with Kv and ASIC ion channels

Tarantula toxins that bind to voltage-sensing domains of voltage-activated ion channels are thought to partition into the membrane and bind to the channel within the bilayer. While no structures of a voltage-sensor toxin bound to a channel have been solved, a structural homolog, psalmotoxin (PcTx1), was recently crystalized in complex with the extracellular domain of an acid sensing ion channel (ASIC). In the present study we use spectroscopic, biophysical and computational approaches to compare membrane interaction properties and channel binding surfaces of PcTx1 with the voltage-sensor toxin guangxitoxin (GxTx-1E). Our results show that both types of tarantula toxins interact with membranes, but that voltage-sensor toxins partition deeper into the bilayer. In addition, our results suggest that tarantula toxins have evolved a similar concave surface for clamping onto α-helices that is effective in aqueous or lipidic physical environments.

The HIV-1 nucleocapsid protein recruits negatively charged lipids to ensure its optimal binding to lipid membranes

The HIV-1 Gag polyprotein precursor composed of the matrix (MA), capsid (CA), nucleocapsid (NC), and p6 domains orchestrates virus assembly via interactions between MA and the cell plasma membrane (PM) on one hand and NC and the genomic RNA on the other hand. As the Gag precursor can adopt a bent conformation, a potential interaction of the NC domain with the PM cannot be excluded during Gag assembly at the PM. To investigate the possible interaction of NC with lipid membranes in the absence of any interference from the other domains of Gag, we quantitatively characterized by fluorescence spectroscopy the binding of the mature NC protein to large unilamellar vesicles (LUVs) used as membrane models. We found that NC, either in its free form or bound to an oligonucleotide, was binding with high affinity (∼ 10(7) M(-1)) to negatively charged LUVs. The number of NC binding sites, but not the binding constant, was observed to decrease with the percentage of negatively charged lipids in the LUV composition, suggesting that NC and NC/oligonucleotide complexes were able to recruit negatively charged lipids to ensure optimal binding. However, in contrast to MA, NC did not exhibit a preference for phosphatidylinositol-(4,5)-bisphosphate. These results lead us to propose a modified Gag assembly model where the NC domain contributes to the initial binding of the bent form of Gag to the PM.

IMPORTANCE

The NC protein is a highly conserved nucleic acid binding protein that plays numerous key roles in HIV-1 replication. While accumulating evidence shows that NC either as a mature protein or as a domain of the Gag precursor also interacts with host proteins, only a few data are available on the possible interaction of NC with lipid membranes. Interestingly, during HIV-1 assembly, the Gag precursor is thought to adopt a bent conformation where the NC domain may interact with the plasma membrane. In this context, we quantitatively characterized the binding of NC, as a free protein or as a complex with nucleic acids, to lipid membranes and showed that the latter constitute a binding platform for NC. Taken together, our data suggest that the NC domain may play a role in the initial binding events of Gag to the plasma membrane during HIV-1 assembly.

Salt as a catalyst in the mitochondria: returning cytochrome c to its native state after it misfolds on the surface of cardiolipin containing membranes

Cytochrome c binding to cardiolipin receptors on the surface of TOCL/DOPC(20% : 80%) liposomes induces a conformational change, which is not reversible after the protein's dissociation at low ionic strength. Addition of 100 mM NaCl switches the protein back to its native state.

Lipid domains in HIV-1 assembly

In CD⁺₄ T cells, HIV-1 buds from the host cell plasma membrane. The viral Gag polyprotein is mainly responsible for this process. However, the intimate interaction of Gag and lipids at the plasma membrane as well as its consequences, in terms of lipids lateral organization and virus assembly, is still under debate. In this review we propose to revisit the role of plasma membrane lipids in HIV-1 Gag targeting and assembly, at the light of lipid membranes biophysics and literature dealing with Gag-lipid interactions.

Implications of salivary protein binding to commensal and pathogenic bacteria

An important function of salivary proteins is to interact with microorganisms that enter the oral cavity. For some microbes, these interactions promote microbial colonization. For others, these interactions are deleterious and result in the elimination of the microbe from the mouth, This paper reviews recent studies of the interaction of salivary proteins with two model bacteria; the commensal species Streptococcus gordonii, and the facultative pathogen Staphylococcus aureus. These organisms selectively interact with a variety of salivary proteins to influence important functions such as bacterial adhesion to surfaces, evasion of host defense, bacterial nutrition and metabolism and gene expression.

Binding of cationic pentapeptides with modified side chain lengths to negatively charged lipid membranes: Complex interplay of electrostatic and hydrophobic interactions

Basic amino acids play a key role in the binding of membrane associated proteins to negatively charged membranes. However, side chains of basic amino acids like lysine do not only provide a positive charge, but also a flexible hydrocarbon spacer that enables hydrophobic interactions. We studied the influence of hydrophobic contributions to the binding by varying the side chain length of pentapeptides with ammonium groups starting with lysine to lysine analogs with shorter side chains, namely omithine (Orn), alpha, gamma-diaminobutyric acid (Dab) and alpha, beta-diaminopropionic acid (Dap). The binding to negatively charged phosphatidylglycerol (PG) membranes was investigated by calorimetry, FT-infrared spectroscopy (FT-IR) and monolayer techniques. The binding was influenced by counteracting and sometimes compensating contributions. The influence of the bound peptides on the lipid phase behavior depends on the length of the peptide side chains. Isothermal titration calorimetry (ITC) experiments showed exothermic and endothermic effects compensating to a different extent as a function of side chain length. The increase in lipid phase transition temperature was more significant for peptides with shorter side chains. FTIR-spectroscopy revealed changes in hydration of the lipid bilayer interface after peptide binding. Using monolayer techniques, the contributions of electrostatic and hydrophobic effects could clearly be observed. Peptides with short side chains induced a pronounced decrease in surface pressure of PG monolayers whereas peptides with additional hydrophobic interactions decreased the surface pressure much less or even lead to an increase, indicating insertion of the hydrophobic part of the side chain into the lipid monolayer.

Lateral dynamics of charged lipids and peripheral proteins in spatially heterogeneous membranes: comparison of continuous and Monte Carlo approaches

Biological membranes are complex environments whose physico-chemical properties are of utmost importance for the understanding of many crucial biological processes. Much attention has been given in the literature to the description of membranes along the z-axis perpendicular to the membrane. Here, we instead consider the lateral dynamics of lipids and peripheral proteins due to their electrostatic interaction. Previously, we constructed a Monte Carlo automaton capable of simulating mutual diffusive dynamics of charged lipids and associated positively charged peptides. Here, we derive and numerically analyze a system of Poisson-Boltzmann-Nernst-Planck (PBNP) equations that provide a mean-field approximation compatible with our Monte Carlo model. The thorough comparison between the mean-field PBNP equations and Monte Carlo simulations demonstrates that both the approaches are in a good qualitative agreement in all tested scenarios. We find that the two methods quantitatively deviate when the local charge density is high, presumably because the Poisson-Boltzmann formalism is applicable in the so-called weak coupling limit, whose validity is restricted to low charge densities. Nevertheless, we conclude that the mean-field PBNP approach provides a good approximation for the considerably more detailed Monte Carlo model at only a fraction of the associated computational cost and allows simulation of the membrane lateral dynamics on the space and time scales relevant for the realistic biological problems.

OPM database and PPM web server: resources for positioning of proteins in membranes

The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topology and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approximately 350 organisms that represent approximately 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielectric constant and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calculating spatial positions in membranes of newly determined proteins structures or theoretical models.

Thermodynamics and mechanics of membrane curvature generation and sensing by proteins and lipids

Research investigating lipid membrane curvature generation and sensing is a rapidly developing frontier in membrane physical chemistry and biophysics. The fast recent progress is based on the discovery of a plethora of proteins involved in coupling membrane shape to cellular membrane function, the design of new quantitative experimental techniques to study aspects of membrane curvature, and the development of analytical theories and simulation techniques that allow a mechanistic interpretation of quantitative measurements. The present review first provides an overview of important classes of membrane proteins for which function is coupled to membrane curvature. We then survey several mechanisms that are assumed to underlie membrane curvature sensing and generation. Finally, we discuss relatively simple thermodynamic/mechanical models that allow quantitative interpretation of experimental observations.

Sum-frequency spectroscopic study of Langmuir monolayers of lipids having oppositely charged headgroups

Sum-frequency vibrational spectroscopy, with the help of surface pressure-area (π-A) isotherm, was used to study lipid Langmuir monolayers composed of molecules with positively and negatively charged headgroups as well as a 1:1 neutral mixture of the two. The spectral profiles of the CH(x) stretch vibrations are similar for all monolayers in the liquid-condensed (LC) phase. They suggest a monolayer structure of closely packed alkyl chains that are nearly all-trans and well oriented along the surface normal. In the liquid-expanded (LE) phase, the spectra of all monolayers appear characteristic of loosely packed chains with significant gauche defects. The OH stretch spectra of interfacial water for both positively and negatively charged monolayers are significantly enhanced in comparison with a neutral water interface, but the phase measurement of SFVS indicates that OH in the two cases points toward the bulk and the interface, respectively. The enhancement results mainly from surface-field-induced polar ordering of interfacial water molecules. For a charge-neutral monolayer composed of an equal number of positively and negatively charged lipid molecules, no such enhancement is observed. This mixed monolayer exhibits a wide range of LC/LE coexistence region extended to very low surface pressure and its CH(x) spectral profile in the coexistence region resembles that of the LC phase. This result suggests that in the LC/LE coexistence region, the mixed monolayer consists of coexisting LC and LE patches in which oppositely charged lipid molecules are homogeneously mixed and dispersed.

On the role of anionic lipids in charged protein interactions with membranes

We investigate the role of anionic lipids in the binding to, and subsequent movement of charged protein groups in lipid membranes, to help understand the role of membrane composition in all membrane-active protein sequences. We demonstrate a small effect of phosphatidylglycerol (PG) lipids on the ability of an arginine (Arg) side chain to bind to, and cross a lipid membrane, despite possessing a neutralizing charge. We observe similar membrane deformations in lipid bilayers composed of phosphatidylcholine (PC) and PC/PG mixtures, with comparable numbers of water and lipid head groups pulled into the bilayer hydrocarbon core, and prohibitively large ~20 kcal/mol barriers for Arg transfer across each bilayer, dropping by just 2-3 kcal/mol due to the binding of PG lipids. We explore the causes of this small effect of introducing PG lipids and offer an explanation in terms of the limited membrane interaction for the choline groups of PC lipids bound to the translocating ion. Our calculations reveal a surprising lack of preference for Arg binding to PG lipids themselves, but a small increase in interfacial binding affinity for lipid bilayers containing PG lipids. These results help to explain the nature of competitive lipid binding to charged protein sequences, with implications for a wide range of membrane binding domains and cell perturbing peptides.

Protein folding in membranes

Separation of cells and organelles by bilayer membranes is a fundamental principle of life. Cellular membranes contain a baffling variety of proteins, which fulfil vital functions as receptors and signal transducers, channels and transporters, motors and anchors. The vast majority of membrane-bound proteins contain bundles of alpha-helical transmembrane domains. Understanding how these proteins adopt their native, biologically active structures in the complex milieu of a membrane is therefore a major challenge in today's life sciences. Here, we review recent progress in the folding, unfolding and refolding of alpha-helical membrane proteins and compare the molecular interactions that stabilise proteins in lipid bilayers. We also provide a critical discussion of a detergent denaturation assay that is increasingly used to determine membrane-protein stability but is not devoid of conceptual difficulties.

Molecular dynamics simulations of PIP2 and PIP3 in lipid bilayers: determination of ring orientation, and the effects of surface roughness on a Poisson-Boltzmann description

Molecular dynamics (MD) simulations of phosphatidylinositol (4,5)-bisphosphate (PIP2) and phosphatidylinositol (3,4,5)-trisphosphate (PIP3) in 1-palmitoyl 2-oleoyl phosphatidylcholine (POPC) bilayers indicate that the inositol rings are tilted approximately 40 degrees with respect to the bilayer surface, as compared with 17 degrees for the P-N vector of POPC. Multiple minima were obtained for the ring twist (analogous to roll for an airplane). The phosphates at position 1 of PIP2 and PIP3 are within an Angström of the plane formed by the phosphates of POPC; lipids in the surrounding shell are depressed by 0.5-0.8 A, but otherwise the phosphoinositides do not substantially perturb the bilayer. Finite size artifacts for ion distributions are apparent for systems of approximately 26 waters/lipid, but, based on simulations with a fourfold increase of the aqueous phase, the phosphoinositide positions and orientations do not show significant size effects. Electrostatic potentials evaluated from Poisson-Boltzmann (PB) calculations show a strong dependence of potential height and ring orientation, with the maxima on the -25 mV surfaces (17.1 +/- 0.1 A for PIP2 and 19.4 +/- 0.3 A for PIP3) occurring near the most populated orientations from MD. These surfaces are well above the background height of 10 A estimated for negatively charged cell membranes, as would be expected for lipids involved in cellular signaling. PB calculations on microscopically flat bilayers yield similar maxima as the MD-based (microscopically rough) systems, but show less fine structure and do not clearly indicate the most probable regions. Electrostatic free energies of interaction with pentalysine are also similar for the rough and flat systems. These results support the utility of a rigid/flat bilayer model for PB-based studies of PIP2 and PIP3 as long as the orientations are judiciously chosen.

EGFR juxtamembrane domain, membranes, and calmodulin: kinetics of their interaction

Calcium/calmodulin (Ca/CaM) binds to the intracellular juxtamembrane domain (JMD) of the epidermal growth factor receptor (EGFR). The basic JMD also binds to acidic lipids in the inner leaflet of the plasma membrane, and this interaction may contribute an extra level of autoinhibition to the receptor. Binding of a ligand to the EGFR produces a rapid increase in intracellular calcium, [Ca2+]i, and thus Ca/CaM. How does Ca/CaM compete with the plasma membrane for the JMD? Does Ca/CaM directly pull the JMD off the membrane or does Ca/CaM only bind to the JMD after it has dissociated spontaneously from the bilayer? To answer this question, we studied the effect of Ca/CaM on the rate of dissociation of fluorescent JMD peptides from phospholipid vesicles by making kinetic stop-flow measurements. Ca/CaM increases the rate of dissociation: an analysis of the differential equations that describe the dissociation shows that Ca/CaM must directly pull the basic JMD peptide off the membrane surface. These measurements lead to a detailed atomic-level mechanism for EGFR activation that reconciles the existence of preformed EGFR dimers/oligomers with the Kuriyan allosteric model for activation of the EGFR kinase domains.

A comparison of the membrane binding properties of C1B domains of PKCgamma, PKCdelta, and PKCepsilon

The C1 domains of classical and novel PKCs mediate their diacylglycerol-dependent translocation. Using fluorescence resonance energy transfer, we studied the contribution of different negatively charged phospholipids and diacylglycerols to membrane binding. Three different C1B domains of PKCs were studied (the classical gamma, and the novel delta and epsilon), together with different lipid mixtures containing three types of acidic phospholipids and three types of activating diacylglycerols. The results show that C1Bgamma and C1Bepsilon exhibit a higher affinity to bind to vesicles containing 1-palmitoyl-2-oleoyl-sn-phosphatidic acid, 1-palmitoyl-2-oleoyl-sn-phoshatidylserine, or 1-palmitoyl-2-oleoyl-sn-phosphatidylglycerol, with C1Bepsilon being the most relevant case because its affinity for POPA-containing vesicles increased by almost two orders of magnitude. When the effect of the diacylglycerol fatty acid composition on membrane binding was studied, the C1Bepsilon domain showed the highest binding affinity to membranes containing 1-stearoyl-oleoyl-sn-glycerol or 1,2-sn-dioleoylglycerol with POPA as the acidic phospholipid. Of the three diacylglycerols used in this study, 1,2-sn-dioleoylglycerol and 1-stearoyl-oleoyl-sn-glycerol showed the highest affinities for each isoenzyme, whereas 1,2-sn-dipalmitoylglycerol; showed the lowest affinity. DSC experiments showed this to be a consequence of the nonfluid conditions of 1,2-sn-dipalmitoylglycerol;-containing systems.

Electrostatic and lipid anchor contributions to the interaction of transducin with membranes: mechanistic implications for activation and translocation

The heterotrimeric G protein transducin is a key component of the vertebrate phototransduction cascade. Transducin is peripherally attached to membranes of the rod outer segment, where it interacts with other proteins at the membrane-cytosol interface. However, upon sustained activation by light, the dissociated G_tα and Gβ₁γ₁ subunits of transducin translocate from the outer segment to other parts of the rod cell. Here we used a computational approach to analyze the interaction strength of transducin and its subunits with acidic lipid bilayers, as well as the range of orientations that they are allowed to occupy on the membrane surface. Our results suggest that the combined constraints of electrostatics and lipid anchors substantially limit the rotational degrees of freedom of the membrane-bound transducin heterotrimer. This may contribute to a faster transducin activation rate by accelerating transducin-rhodopsin complex formation. Notably, the membrane interactions of the dissociated transducin subunits are very different from those of the heterotrimer. As shown previously, Gβ₁γ₁ experiences significant attractive interactions with negatively charged membranes, whereas our new results suggest that G_tα is electrostatically repelled by such membranes. We suggest that this repulsion could facilitate the membrane dissociation and intracellular translocation of G_tα. Moreover, based on similarities in sequence and electrostatic properties, we propose that the properties described for transducin are common to its homologs within the G_i subfamily. In a broader view, this work exemplifies how the activity-dependent association and dissociation of a G protein can change both the affinity for membranes and the range of allowed orientations, thereby modulating G protein function.