Lipid-protein interactions of integral membrane proteins: a comparative simulation study

Probing the Hydrophobic Region of a Lipid Bilayer at Specific Depths Using Vibrational Spectroscopy

A novel spectroscopic approach for studying the flexibility and mobility in the hydrophobic interior of lipid bilayers at specific depths is proposed. A set of test compounds featuring an azido moiety and a cyano or carboxylic acid moiety, connected by an alkyl chain of different lengths, was synthesized. FTIR data and molecular dynamics calculations indicated that the test compounds in a bilayer are oriented so that the cyano or carboxylic acid moiety is located in the lipid head-group region, while the azido group stays inside the bilayer at the depth determined by its alkyl chain length. We found that the asymmetric stretching mode of the azido group (νN3) can serve as a reporter of the membrane interior dynamics. FTIR and two-dimensional infrared (2DIR) studies were performed at different temperatures, ranging from 22 to 45 °C, covering the Lβ-Lα phase transition temperature of dipalmitoylphosphatidylcholine (∼41 °C). The width of the νN3 peak was found to be very sensitive to the phase transition and to the temperature in general. We introduced an order parameter, SN3, which characterizes restrictions to motion inside the bilayer. 2DIR spectra of νN3 showed different extents of inhomogeneity at different depths in the bilayer, with the smallest inhomogeneity in the middle of the leaflet. The spectral diffusion dynamics of the N3 peak was found to be dependent on the depth of the N3 group location in the bilayer. The obtained results enhance our understanding of the bilayer dynamics and can be extended to investigate membranes with more complex compositions.

Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

Volleying plasma membrane proteins from birth to death: Role of J-domain proteins

The function, stability, and turnover of plasma membrane (PM) proteins are crucial for cellular homeostasis. Compared to soluble proteins, quality control of plasma membrane proteins is extremely challenging. Failure to meet the high quality control standards is detrimental to cellular and organismal health. J-domain proteins (JDPs) are among the most diverse group of chaperones that collaborate with other chaperones and protein degradation machinery to oversee cellular protein quality control (PQC). Although fragmented, the available literature from different models, including yeast, mammals, and plants, suggests that JDPs assist PM proteins with their synthesis, folding, and trafficking to their destination as well as their degradation, either through endocytic or proteasomal degradation pathways. Moreover, some JDPs interact directly with the membrane to regulate the stability and/or functionality of proteins at the PM. The deconvoluted picture emerging is that PM proteins are relayed from one JDP to another throughout their life cycle, further underscoring the versatility of the Hsp70:JDP machinery in the cell.

Lipid-mediated prestin organization in outer hair cell membranes and its implications in sound amplification

Prestin is a high-density motor protein in the outer hair cells (OHCs), whose conformational response to acoustic signals alters the shape of the cell, thereby playing a major role in sound amplification by the cochlea. Despite recent structures, prestin's intimate interactions with the membrane, which are central to its function remained unresolved. Here, employing a large set (collectively, more than 0.5 ms) of coarse-grained molecular dynamics simulations, we demonstrate the impact of prestin's lipid-protein interactions on its organization at densities relevant to the OHCs and its effectiveness in reshaping OHCs. Prestin causes anisotropic membrane deformation, which mediates a preferential membrane organization of prestin where deformation patterns by neighboring copies are aligned constructively. The resulting reduced membrane rigidity is hypothesized to maximize the impact of prestin on OHC reshaping. These results demonstrate a clear case of protein-protein cooperative communication in membrane, purely mediated by interactions with lipids.

The effects of molecular weight and orientation on the membrane permeation and partitioning of polycyclic aromatic hydrocarbons: a computational study

Membrane permeation and the partitioning of polycyclic aromatic hydrocarbons (PAHs) are crucial aspects affecting their carcinogenicity and mutagenicity. However, a clear understanding of these processes is still rare due to the difficulty of determining the details experimentally. Here, the interactions between PAHs and lipid bilayers were studied by molecular simulations, mainly to check the influence of molecular weight and orientation. The liposome-water partition coefficient (K_LW), transmembrane time (τ), and permeability coefficient (P) of the PAHs were calculated by integrating free energy profiles from umbrella sampling. For selected PAHs, the membrane adsorption is a spontaneous process. The preferred location is near the CC bond and the orientation is related to the molecular structure. The P values of all the PAHs are basically the same order of magnitude, which means that the molecular weight contributes little to the process. As for K_LW and τ, they show obvious increases with different molecular weights. Unconstrained simulations showed that a flat orientation on the membrane surface would prevent PAHs from being transported through the membrane. Highly hydrophobic driving forces are not always good for the absorption of PAHs, especially the formation of aggregates. In addition, the orientations and energetic barriers of PAHs near the midplane of the lipid bilayer explain the different transitions of high- and low-weight PAHs. This work provides molecular level details relating to the interactions of PAHs with lipid membranes, with significance for understanding the health effects of PAHs.

The role of the envelope protein in the stability of a coronavirus model membrane against an ethanolic disinfectant

Ethanol is highly effective against various enveloped viruses and can disable the virus by disintegrating the protective envelope surrounding it. The interactions between the coronavirus envelope (E) protein and its membrane environment play key roles in the stability and function of the viral envelope. By using molecular dynamics simulation, we explore the underlying mechanism of ethanol-induced disruption of a model coronavirus membrane and, in detail, interactions of the E-protein and lipids. We model the membrane bilayer as N-palmitoyl-sphingomyelin and 1-palmitoyl-2-oleoylphosphatidylcholine lipids and the coronavirus E-protein. The study reveals that ethanol causes an increase in the lateral area of the bilayer along with thinning of the bilayer membrane and orientational disordering of lipid tails. Ethanol resides at the head-tail region of the membrane and enhances bilayer permeability. We found an envelope-protein-mediated increase in the ordering of lipid tails. Our simulations also provide important insights into the orientation of the envelope protein in a model membrane environment. At ∼25 mol. % of ethanol in the surrounding ethanol-water phase, we observe disintegration of the lipid bilayer and dislocation of the E-protein from the membrane environment.

A Single Mutation in the Outer Lipid-Facing Helix of a Pentameric Ligand-Gated Ion Channel Affects Channel Function Through a Radially-Propagating Mechanism

Pentameric ligand-gated ion channels (pLGICs) mediate fast synaptic transmission and are crucial drug targets. Their gating mechanism is triggered by ligand binding in the extracellular domain that culminates in the opening of a hydrophobic gate in the transmembrane domain. This domain is made of four α-helices (M1 to M4). Recently the outer lipid-facing helix (M4) has been shown to be key to receptor function, however its role in channel opening is still poorly understood. It could act through its neighboring helices (M1/M3), or via the M4 tip interacting with the pivotal Cys-loop in the extracellular domain. Mutation of a single M4 tyrosine (Y441) to alanine renders one pLGIC-the 5-HT3A receptor-unable to function despite robust ligand binding. Using Y441A as a proxy for M4 function, we here predict likely paths of Y441 action using molecular dynamics, and test these predictions with functional assays of mutant receptors in HEK cells and Xenopus oocytes using fluorescent membrane potential sensitive dye and two-electrode voltage clamp respectively. We show that Y441 does not act via the M4 tip or Cys-loop, but instead connects radially through M1 to a residue near the ion channel hydrophobic gate on the pore-lining helix M2. This demonstrates the active role of the M4 helix in channel opening.

Cholesterol content in the membrane promotes key lipid-protein interactions in a pentameric serotonin-gated ion channel

Pentameric ligand-gated ion channels (pLGICs), embedded in the lipid membranes of nerve cells, mediate fast synaptic transmission and are major pharmaceutical targets. Because of their complexity and the limited knowledge of their structure, their working mechanisms have still to be fully unraveled at the molecular level. Over the past few years, evidence that the lipid membrane may modulate the function of membrane proteins, including pLGICs, has emerged. Here, we investigate, by means of molecular dynamics simulations, the behavior of the lipid membrane at the interface with the 5-HT_3A receptor (5-HT_3AR), a representative pLGIC which is the target of nausea-suppressant drugs, in a nonconductive state. Three lipid compositions are studied, spanning different concentrations of the phospholipids, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine, and of cholesterol, hence a range of viscosities. A variety of lipid interactions and persistent binding events to different parts of the receptor are revealed in the investigated models, providing snapshots of the dynamical environment at the membrane-receptor interface. Some of these events result in lipid intercalation within the transmembrane domain, and others reach out to protein key sections for signal transmission and receptor activation, such as the Cys-loop and the M2-M3 loop. In particular, phospholipids, with their long hydrophobic tails, play an important role in these interactions, potentially providing a bridge between these two structures. A higher cholesterol content appears to promote lipid persistent binding to the receptor.

Dimerization of Aβ40 inside dipalmitoylphosphatidylcholine bilayer and its effect on bilayer integrity: Atomistic simulation at three temperatures

Amyloid-beta (Aβ) protein is related to Alzheimer disease (AD), and various experiments have shown that oligomers as small as dimers are cytotoxic. Recent studies have concluded that interactions of Aβ with neuronal cell membranes lead to disruption of membrane integrity and toxicity and they play a key role in the development of AD. Molecular dynamics (MD) simulations have been used to investigate Aβ in aqueous solution and membranes. We have previously studied monomeric Aβ40 embedded in dipalmitoylphosphatidylcholine (DPPC) membrane using MD simulations. Here, we explore interactions of two Aβ40 peptides in DPPC bilayer and its consequences on dimer distribution in a lipid bilayer and on the secondary structure of the peptides. We explored that N-terminals played an important role in dimeric Aβ peptide aggregations and Aβ-bilayer interactions, while C-terminals bound peptides to bilayer like anchors. We did not observe exiting of peptides in our simulations although we observed insertion of peptides into the core of bilayer in some of our simulations. So it seems that the presence of Aβ on membrane surface increases its aggregation rate, and as diffusion occurs in two dimensions, it can increase the probability of interpeptide interactions. We found that dimeric Aβ, like monomeric one, had the ability to cause structural destabilization of DPPC membrane, which in turn might ultimately lead to cell death in an in vivo system. This information could have important implications for understanding the affinity of Aβ oligomers (here dimer) for membranes and the mechanism of Aβ oligomer toxicity in AD.

Assessing the extent of the structural and dynamic modulation of membrane lipids due to pore forming toxins: insights from molecular dynamics simulations

Infections caused by many virulent bacterial strains are triggered by the release of pore forming toxins (PFTs), which form oligomeric transmembrane pore complexes on the target plasma membrane. The spatial extent of the perturbation to the surrounding lipids during pore formation is relatively unexplored. Using all-atom molecular dynamics simulations, we investigate the changes in the structure and dynamics of lipids in a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer in the presence of contrasting PFTs. Cytolysin A (ClyA), an α toxin with its inserted wedge shaped bundle of inserted α helices, induces significant asymmetry across the membrane leaflets in comparison with α hemolysin (AHL), a β toxin. Despite the differences in hydrophobic mismatch and uniquely different topologies of the two oligomers, perturbations to lipid order as reflected in the tilt angle and order parameters and membrane thinning are short ranged, lying within ∼2.5 nm from the periphery of either pore complex, and commensurate with distances typically associated with van der Waals forces. In contrast, the spatial extent of perturbations to the lipid dynamics extends outward to at least 4 nm for both proteins, and the continuous survival probabilities reveal the presence of a tightly bound shell of lipids in this region. Displacement probability distributions show long tails and the distinctly non-Gaussian features reflect the induced dynamic heterogeneity. A detailed profiling of the protein-lipid contacts with tyrosine, tryptophan, lysine and arginine residues shows increased non-polar contacts in the cytoplasmic leaflet for both PFTs, with a higher number of atomic contacts in the case of AHL in the extracellular leaflet due to the mushroom-like topology of the pore complex. The short ranged nature of the perturbations observed in this simple one component membrane suggests inherent plasticity of membrane lipids enabling the recovery of the structure and membrane fluidity even in the presence of these large oligomeric transmembrane protein assemblies. This observation has implications in membrane repair processes such as budding or vesicle fusion events used to mitigate PFT virulence, where the underlying lipid dynamics and fluidity in the vicinity of the pore complex are expected to play an important role.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Influence of membrane lipid composition on the structure and activity of γ-secretase

γ-Secretase (GS) is a multi-subunit membrane-embedded aspartyl protease that cleaves more than 80 integral membrane proteins, including the amyloid precursor protein (APP) to produce the amyloid-β (Aβ) peptide. Oligomerization and aggregation of the 42-amino acid length Aβ isoform in the brain has been associated with the development and progression of Alzheimer's disease (AD). Based on recent experimental structural studies and using multiscale computational modeling approaches, the conformational states and protein-membrane interactions of the GS complex embedded in six homogeneous and six heterogeneous lipid bilayers were characterized. In order to identify potential lipid and cholesterol binding sites, GS regions with high lipid/cholesterol occupancy values were analyzed using atomistic and coarse-grained simulations. Long lipid residence times were observed to be correlated with a large number of hydrogen bonds between the charged headgroups and key GS amino acids. This observation provides a plausible explanation for the inhibition of GS by charged lipids observed in previous experimental studies. Computed lateral pressure profiles suggest that higher transmembrane pressures favor active state conformations of the catalytic subunit. A probable mechanism for the regulation of the local stress response in cholesterol-rich multicomponent lipid bilayers is identified. Finally, it is demonstrated that interactions between the nicastrin extracellular domain and lipid headgroups leads to a compact structural conformation of the GS complex. Overall, this study provides valuable insight into the effect of bilayer lipid composition on the GS structural ensemble and its function.

The Role of Lipid Interactions in Simulations of the α-Hemolysin Ion-Channel-Forming Toxin

Molecular dynamics simulations were performed to describe the function of the ion-channel-forming toxin α-hemolysin (αHL) in lipid membranes that were composed of either 1,2-diphytanoyl-sn-glycero-3-phospho-choline or 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline. The simulations highlight the importance of lipid type in maintaining αHL structure and function, enabling direct comparison to experiments for biosensing applications. We determined that although the two lipids studied are similar in structure, 1,2-diphytanoyl-sn-glycero-3-phospho-choline membranes better match the hydrophobic thickness of αHL compared to 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-choline membranes. This hydrophobic match is essential to maintaining proper alignment of β-sheet loops at the trans entrance of αHL, which, when disrupted, creates an additional constriction to ion flow that decreases the channel current below experimental values and creates greater variability in channel conductance. Agreement with experiments was further improved with sufficient lipid membrane equilibration and allowed the discrimination of subtle αHL conduction states with lipid type. Finally, we explore the effects of truncating the extramembrane cap of αHL and its role in maintaining proper alignment of αHL in the membrane and channel conductance. Our results demonstrate the essential role of lipid type and lipid-protein interactions in simulations of αHL and will considerably improve the interpretation of experimental data.

Structural Insights into the Yersinia pestis Outer Membrane Protein Ail in Lipid Bilayers

Yersinia pestis the causative agent of plague, is highly pathogenic and poses very high risk to public health. The outer membrane protein Ail (Adhesion invasion locus) is one of the most highly expressed proteins on the cell surface of Y. pestis, and a major target for the development of medical countermeasures. Ail is essential for microbial virulence and is critical for promoting the survival of Y. pestis in serum. Structures of Ail have been determined by X-ray diffraction and solution NMR spectroscopy, but the protein's activity is influenced by the detergents in these samples, underscoring the importance of the surrounding environment for structure-activity studies. Here we describe the backbone structure of Ail, determined in lipid bilayer nanodiscs, using solution NMR spectroscopy. We also present solid-state NMR data obtained for Ail in membranes containing lipopolysaccharide (LPS), a major component of the bacterial outer membranes. The protein in lipid bilayers, adopts the same eight-stranded β-barrel fold observed in the crystalline and micellar states. The membrane composition, however, appears to have a marked effect on protein dynamics, with LPS enhancing conformational order and slowing down the ¹⁵N transverse relaxation rate. The results provide information about the way in which an outer membrane protein inserts and functions in the bacterial membrane.

Structural features and lipid binding domain of tubulin on biomimetic mitochondrial membranes

Dimeric tubulin, an abundant water-soluble cytosolic protein known primarily for its role in the cytoskeleton, is routinely found to be associated with mitochondrial outer membranes, although the structure and physiological role of mitochondria-bound tubulin are still unknown. There is also no consensus on whether tubulin is a peripheral membrane protein or is integrated into the outer mitochondrial membrane. Here the results of five independent techniques-surface plasmon resonance, electrochemical impedance spectroscopy, bilayer overtone analysis, neutron reflectometry, and molecular dynamics simulations-suggest that α-tubulin's amphipathic helix H10 is responsible for peripheral binding of dimeric tubulin to biomimetic "mitochondrial" membranes in a manner that differentiates between the two primary lipid headgroups found in mitochondrial membranes, phosphatidylethanolamine and phosphatidylcholine. The identification of the tubulin dimer orientation and membrane-binding domain represents an essential step toward our understanding of the complex mechanisms by which tubulin interacts with integral proteins of the mitochondrial outer membrane and is important for the structure-inspired design of tubulin-targeting agents.

The molecular behavior of a single β-amyloid inside a dipalmitoylphosphatidylcholine bilayer at three different temperatures: An atomistic simulation study: Aβ interaction with DPPC: Atomistic simulation

The behavior of a single Aβ40 molecule within a dipalmitoylphosphatidylcholine (DPPC) bilayer was studied by all-atom molecular dynamics simulations. The effect of membrane structure was investigated on Aβ40 behavior, secondary structure, and insertion depth. Simulations were performed at three temperatures (323, 310, and 300 K) to probe three different bilayer fluidities. Results show that at all above temperatures, the peptide contains two short helices, coil, bend, and turn structures. At 300 K, the peptide contains a region with β structure in C-terminal region. Our results also show that Aβ decreases the bilayer thickness and the order of lipids in its vicinity which leads to water insertion into the bilayer and concomitant increase in the local fluidity. The peptide remains embedded in the bilayer at all temperatures, and become inserted into the bilayer up to several residues at 323 and 310 K. At 310 and 300 K, the dominant interaction energy between Aβ and bilayer changes from electrostatic to van der Waals. It can be proposed that at higher temperatures (e.g., 323 K), Lys28 and the C-terminal region of the peptide play the role of two anchors that keep Aβ inside the top leaflet. This study demonstrates that Aβ molecule can perturb the integrity of cellular membranes. Proteins 2017; 85:1298-1310. © 2017 Wiley Periodicals, Inc.

Direct Interaction between the Voltage Sensors Produces Cooperative Sustained Deactivation in Voltage-gated H+ Channel Dimers

The voltage-gated H(+) channel (Hv) is a voltage sensor domain-like protein consisting of four transmembrane segments (S1-S4). The native Hv structure is a homodimer, with the two channel subunits functioning cooperatively. Here we show that the two voltage sensor S4 helices within the dimer directly cooperate via a π-stacking interaction between Trp residues at the middle of each segment. Scanning mutagenesis showed that Trp situated around the original position provides the slow gating kinetics characteristic of the dimer's cooperativity. Analyses of the Trp mutation on the dimeric and monomeric channel backgrounds and analyses with tandem channel constructs suggested that the two Trp residues within the dimer are functionally coupled during Hv deactivation but are less so during activation. Molecular dynamics simulation also showed direct π-stacking of the two Trp residues. These results provide new insight into the cooperative function of voltage-gated channels, where adjacent voltage sensor helices make direct physical contact and work as a single unit according to the gating process.

Glu106 targeted inhibitors of ORAI1 as potential Ca2+ release-activated Ca2+ (CRAC) channel blockers - molecular modeling and docking studies

Calcium release-activated calcium modulator 1(ORAI1) is an integral component of the calcium release-activated calcium channel (CRAC) channel complex and plays a central role in regulating Ca^2 +concentrations in T-lymphocytes. It is critical for many physiological processes, including cell-proliferation, cytokine production and activation of the immune system. Loss of ORAI1 function is linked with rheumatoid arthritis (RA) and hence pharmacological blockers of ORAI1 could be potential therapeutic agents for the treatment of RA. In this study, we have used a high-throughput screening approach to inhibit the binding of Ca²⁺toward ORAI1 and the interactions are verified through induced fit docking. The results hint that these compounds act by possibly binding with, and thereby blocking Ca²⁺-binding with ORAI1 (E106). The molecular dynamics (MD) simulations shows strong support toward the hit compounds by showing the ligand potency throughout the simulation timescale of 30 ns. We have thus identified a novel class of highly stable, potential lead compounds that directly bind with the selectivity filter region E106 and block Ca²⁺ binding on ORAI1. This resulting alteration in the pore geometry of ORAI1 due to the strong blocking mechanism of lead compounds will greatly diminish its function and the downstream activities that result from the same including decreased production of cytokines in autoimmune disorders. This study may lay the foundation for finding novel lead compounds for clinical trials that could positively modulate the course of autoimmune disorders with ORAI1 as its specific target.

Protein-fluctuation-induced water-pore formation in ion channel voltage-sensor translocation across a lipid bilayer membrane

We have applied a combined fluorescence microscopy and single-ion-channel electric current recording approach, correlating with molecular dynamics (MD) simulations, to study the mechanism of voltage-sensor domain translocation across a lipid bilayer. We use the colicin Ia ion channel as a model system, and our experimental and simulation results show the following: (1) The open-close activity of an activated colicin Ia is not necessarily sensitive to the amplitude of the applied cross-membrane voltage when the cross-membrane voltage is around the resting potential of excitable membranes; and (2) there is a significant probability that the activation of colicin Ia occurs by forming a transient and fluctuating water pore of ∼15 Å diameter in the lipid bilayer membrane. The location of the water-pore formation is nonrandom and highly specific, right at the insertion site of colicin Ia charged residues in the lipid bilayer membrane, and the formation is intrinsically associated with the polypeptide conformational fluctuations and solvation dynamics. Our results suggest an interesting mechanistic pathway for voltage-sensitive ion channel activation, and specifically for translocation of charged polypeptide chains across the lipid membrane under a transmembrane electric field: the charged polypeptide domain facilitates the formation of hydrophilic water pore in the membrane and diffuses through the hydrophilic pathway across the membrane; i.e., the charged polypeptide chain can cross a lipid membrane without entering into the hydrophobic core of the lipid membrane but entirely through the aqueous and hydrophilic environment to achieve a cross-membrane translocation. This mechanism sheds light on the intensive and fundamental debate on how a hydrophilic and charged peptide domain diffuses across the biologically inaccessible high-energy barrier of the hydrophobic core of a lipid bilayer: The peptide domain does not need to cross the hydrophobic core to move across a lipid bilayer.

Membrane-induced structural rearrangement and identification of a novel membrane anchor in talin F2F3

Experimental challenges associated with characterization of the membrane-bound form of talin have prevented us from understanding the molecular mechanism of its membrane-dependent integrin activation. Here, utilizing what we believe to be a novel membrane mimetic model, we present a reproducible model of membrane-bound talin observed across multiple independent simulations. We characterize both local and global membrane-induced structural transitions that successfully reconcile discrepancies between biochemical and structural studies and provide insight into how talin might modulate integrin function. Membrane binding of talin, captured in unbiased simulations, proceeds through three distinct steps: initial electrostatic recruitment of the F2 subdomain to anionic lipids via several basic residues; insertion of an initially buried, conserved hydrophobic anchor into the membrane; and association of the F3 subdomain with the membrane surface through a large, interdomain conformational change. These latter two steps, to our knowledge, have not been observed or described previously. Electrostatic analysis shows talin F2F3 to be highly polarized, with a highly positive underside, which we attribute to the initial electrostatic recruitment, and a negative top face, which can help orient the protein optimally with respect to the membrane, thereby reducing the number of unproductive membrane collision events.

All-atom molecular dynamics simulation of a photosystem i/detergent complex

All-atom molecular dynamics (MD) simulation was used to investigate the solution structure and dynamics of the photosynthetic pigment-protein complex photosystem I (PSI) from Thermosynechococcus elongatus embedded in a toroidal belt of n-dodecyl-β-d-maltoside (DDM) detergent. Evaluation of root-mean-square deviations (RMSDs) relative to the known crystal structure show that the protein complex surrounded by DDM molecules is stable during the 200 ns simulation time, and root-mean-square fluctuation (RMSF) analysis indicates that regions of high local mobility correspond to solvent-exposed regions such as turns in the transmembrane α-helices and flexible loops on the stromal and lumenal faces. Comparing the protein-detergent complex to a pure detergent micelle, the detergent surrounding the PSI trimer is found to be less densely packed but with more ordered detergent tails, contrary to what is seen in most lipid bilayer models. We also investigated any functional implications for the observed conformational dynamics and protein-detergent interactions, discovering interesting structural changes in the psaL subunits associated with maintaining the trimeric structure of the protein. Importantly, we find that the docking of soluble electron mediators such as cytochrome c6 and ferredoxin to PSI is not significantly impacted by the solubilization of PSI in detergent.

Assembly and stability of Salmonella enterica ser. Typhi TolC protein in POPE and DMPE

In this work we assessed the suitability of two different lipid membranes for the simulation of a TolC protein from Salmonella enterica serovar Typhi. The TolC protein family is found in many pathogenic Gram-negative bacteria including Vibrio cholera and Pseudomonas aeruginosa and acts as an outer membrane channel for expulsion of drug and toxin from the cell. In S. typhi, the causative agent for typhoid fever, the TolC outer membrane protein is an antigen for the pathogen. The lipid environment is an important modulator of membrane protein structure and function. We evaluated the conformation of the TolC protein in the presence of DMPE and POPE bilayers using molecular dynamics simulation. The S. typhi TolC protein exhibited similar conformational dynamics to TolC and its homologues. Conformational flexibility of the protein is seen in the C-terminal, extracellular loops, and α-helical region. Despite differences in the two lipids, significant similarities in the motion of the protein in POPE and DMPE were observed, including the rotational motion of the C-terminal residues and the partially open extracellular loops. However, analysis of the trajectories demonstrated effects of hydrophobic matching of the TolC protein in the membrane, particularly in the lengthening of the lipids and subtle movements of the protein's β-barrel towards the lower leaflet in DMPE. The study exhibited the use of molecular dynamics simulation in revealing the differential effect of membrane proteins and lipids on each other. In this study, POPE is potentially a more suitable model for future simulation of the S. typhi TolC protein.

Free energy change for insertion of charged, monolayer-protected nanoparticles into lipid bilayers

Charged, monolayer-protected gold nanoparticles (AuNPs) with core diameters smaller than 10 nm have recently emerged as a prominent class of nanomaterial for use in targeted drug delivery and biosensing. In particular, recent experimental studies showed that AuNPs protected by a binary mixture of purely hydrophobic and anionic, end-functionalized alkanethiol ligands were able to spontaneously penetrate through cell membranes via a non-endocytic, non-disruptive mechanism. The critical step in the penetration process is a fusion step during which the AuNPs insert into the hydrophobic core of the bilayer. This fusion step is driven by hydrophobic forces as inserted AuNPs minimize their exposed hydrophobic surface area and thereby lower their free energy compared to particles in the bulk. Here, we explore the effect of the large parameter space of composition, size, ligand length, morphology, and hydrophobicity strength on the change in the free energy upon insertion. Using a newly developed implicit bilayer, implicit solvent simulation model, our work shows that there is a size cutoff for insertion that has a strong dependence on surface composition and ligand chemistry. Our results agree well with previous experimental findings for a particular value of the hydrophobicity strength. This work provides physical insight that may be used to both understand the insertion of AuNPs into bilayers and guide the design of monolayers to either encourage or inhibit insertion.

Molecular dynamics simulations reveal the HIV-1 Vpu transmembrane protein to form stable pentamers

The human immunodeficiency virus type I (HIV-1) Vpu protein is 81 residues long and has two cytoplasmic and one transmembrane (TM) helical domains. The TM domain oligomerizes to form a monovalent cation selective ion channel and facilitates viral release from host cells. Exactly how many TM domains oligomerize to form the pore is still not understood, with experimental studies indicating the existence of a variety of oligomerization states. In this study, molecular dynamics (MD) simulations were performed to investigate the propensity of the Vpu TM domain to exist in tetrameric, pentameric, and hexameric forms. Starting with an idealized α-helical representation of the TM domain, a thorough search for the possible orientations of the monomer units within each oligomeric form was carried out using replica-exchange MD simulations in an implicit membrane environment. Extensive simulations in a fully hydrated lipid bilayer environment on representative structures obtained from the above approach showed the pentamer to be the most stable oligomeric state, with interhelical van der Waals interactions being critical for stability of the pentamer. Atomic details of the factors responsible for stable pentamer structures are presented. The structural features of the pentamer models are consistent with existing experimental information on the ion channel activity, existence of a kink around the Ile17, and the location of tetherin binding residues. Ser23 is proposed to play an important role in ion channel activity of Vpu and possibly in virus propagation.

Identification of novel cholesterol-binding regions in Kir2 channels

Inwardly rectifying potassium (Kir) channels play an important role in setting the resting membrane potential and modulating membrane excitability. We have recently shown that cholesterol regulates representative members of the Kir family and that in the majority of the cases, cholesterol suppresses channel function. Furthermore, recent data indicate that cholesterol regulates Kir channels by specific sterol-protein interactions, yet the location of the cholesterol binding site in Kir channels is unknown. Using a combined computational-experimental approach, we show that cholesterol may bind to two nonanular hydrophobic regions in the transmembrane domain of Kir2.1 located between adjacent subunits of the channel. The location of the binding regions suggests that cholesterol modulates channel function by affecting the hinging motion at the center of the pore-lining transmembrane helix that underlies channel gating either directly or through the interface between the N and C termini of the channel.

Structural properties of model phosphatidylcholine flippases

Lipid translocation from one lipid bilayer leaflet to the other, termed flip-flop, is required for the distribution of newly synthesized phospholipids during membrane biogenesis. However, a dedicated biogenic lipid flippase has not yet been identified. Here, we show that the efficiency by which model transmembrane peptides facilitate flip of reporter lipids with different headgroups critically depends on their content of helix-destabilizing residues, the charge state of polar flanking residues, and the composition of the host membrane. In particular, increased backbone dynamics of the transmembrane helix relates to its increased ability to flip lipids with phosphatidylcholine and phosphatidylserine headgroups, whereas a more rigid helix favors phosphatidylethanolamine flip. Further, the transmembrane domains of many SNARE protein subtypes share essential features with the dynamic model peptides. Indeed, recombinant SNAREs possess significant lipid flippase activity.

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

Cholesterol and POPC segmental order parameters in lipid membranes: solid state 1H-13C NMR and MD simulation studies

The concentration of cholesterol in cell membranes affects membrane fluidity and thickness, and might regulate different processes such as the formation of lipid rafts. Since interpreting experimental data from biological membranes is rather intricate, investigations on simple models with biological relevance are necessary to understand the natural systems. We study the effect of cholesterol on the molecular structure of multi-lamellar vesicles (MLVs) composed of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a phospholipid ubiquitous in cell membranes, with compositions in the range 0-60 mol% cholesterol. Order parameters, |S(CH)|, are experimentally determined by using (1)H-(13)C solid-state nuclear magnetic resonance (NMR) spectroscopy with segmental detail for all parts of both the cholesterol and POPC molecules, namely the ring system and alkyl chain of the sterol, as well as the glycerol backbone, choline headgroup and the sn-1 and sn-2 acyl chains of POPC. With increasing cholesterol concentration the acyl chains gradually adopt a more extended conformation while the orientation and dynamics of the polar groups are rather unaffected. Additionally, we perform classical molecular dynamics simulations on virtual bilayers mimicking the POPC-cholesterol MLVs investigated by NMR. Good agreement between experiments and simulations is found for the cholesterol alignment in the bilayer and for the |S(CH)| profiles of acyl chains below 15 mol% cholesterol. Deviations occur for the choline headgroup and glycerol backbone parts of POPC, as well as for the phospholipid and cholesterol alkyl chains at higher cholesterol concentrations. The unprecedented detail of the NMR data enables a more complete comparison between simulations and experiments on POPC-cholesterol bilayers and may aid in developing more realistic model descriptions of biological membranes.

Lipid recognition propensities of amino acids in membrane proteins from atomic resolution data

Background

Protein-lipid interactions play essential roles in the conformational stability and biological functions of membrane proteins. However, few of the previous computational studies have taken into account the atomic details of protein-lipid interactions explicitly.

Results

To gain an insight into the molecular mechanisms of the recognition of lipid molecules by membrane proteins, we investigated amino acid propensities in membrane proteins for interacting with the head and tail groups of lipid molecules. We observed a common pattern of lipid tail-amino acid interactions in two different data sources, crystal structures and molecular dynamics simulations. These interactions are largely explained by general lipophilicity, whereas the preferences for lipid head groups vary among individual proteins. We also found that membrane and water-soluble proteins utilize essentially an identical set of amino acids for interacting with lipid head and tail groups.

Conclusions

We showed that the lipophilicity of amino acid residues determines the amino acid preferences for lipid tail groups in both membrane and water-soluble proteins, suggesting that tightly-bound lipid molecules and lipids in the annular shell interact with membrane proteins in a similar manner. In contrast, interactions between lipid head groups and amino acids showed a more variable pattern, apparently constrained by each protein's specific molecular function.

Understanding functional residues of the cannabinoid CB1

The brain cannabinoid (CB(1)) receptor that mediates numerous physiological processes in response to marijuana and other psychoactive compounds is a G protein coupled receptor (GPCR) and shares common structural features with many rhodopsin class GPCRs. For the rational development of therapeutic agents targeting the CB(1) receptor, understanding of the ligand-specific CB(1) receptor interactions responsible for unique G protein signals is crucial. For a more than a decade, a combination of mutagenesis and computational modeling approaches has been successfully employed to study the ligand-specific CB(1) receptor interactions. In this review, after a brief discussion about recent advances in understanding of some structural and functional features of GPCRs commonly applicable to the CB(1) receptor, the CB(1) receptor functional residues reported from mutational studies are divided into three different types, ligand binding (B), receptor stabilization (S) and receptor activation (A) residues, to delineate the nature of the binding pockets of anandamide, CP55940, WIN55212-2 and SR141716A and to describe the molecular events of the ligand-specific CB(1) receptor activation from ligand binding to G protein signaling. Taken these CB(1) receptor functional residues, some of which are unique to the CB(1) receptor, together with the biophysical knowledge accumulated for the GPCR active state, it is possible to propose the early stages of the CB(1) receptor activation process that not only provide some insights into understanding molecular mechanisms of receptor activation but also are applicable for identifying new therapeutic agents by applying the validated structure-based approaches, such as virtual high throughput screening (HTS) and fragment-based approach (FBA).

Molecular dynamics simulations of the interactions of kinin peptides with an anionic POPG bilayer

We have performed molecular dynamics simulations of peptide hormone bradykinin (BK) and its fragment des-Arg9-BK in the presence of an anionic lipid bilayer, with an aim toward delineating the mechanism of action related to their bioactivity. Starting from the initial aqueous environment, both of the peptides are quickly adsorbed and stabilized on the cell surface. Whereas BK exhibits a stronger interaction with the membrane and prefers to stay on the interface, des-Arg9-BK, with the loss of C-terminal Arg, penetrates further. The heterogeneous lipid-water interface induces β-turn-like structure in the otherwise inherently flexible peptides. In the membrane-bound state, we observed C-terminal β-turn formation in BK, whereas for des-Arg9-BK, with the deletion of Arg9, turn formation occurred in the middle of the peptide. The basic Arg residues anchor the peptide to the bilayer by strong electrostatic interactions with charged lipid headgroups. Simulations with different starting orientations of the peptides with respect to the bilayer surface lead to the same observations, namely, the relative positioning of the peptides on the membrane surface, deeper penetration of the des-Arg9-BK, and the formation of turn structures. The lipid headgroups adjacent to the bound peptides become substantially tilted, causing bilayer thinning near the peptide contact region and increase the degree of disorder in nearby lipids. Again, because of hydrogen bonding with the peptide, the neighboring lipid's polar heads exhibit considerably reduced flexibility. Corroborating findings from earlier experiments, our results provide important information about how the lipid environment promotes peptide orientation/conformation and how the peptide adapts to the environment.

The potassium channel KcsA: a model protein in studying membrane protein oligomerization and stability of oligomeric assembly?

Many membrane proteins are functional as stable oligomers. An understanding of the conditions that elicit and enhance oligomerization is important in many therapeutics. In this regard, protein-protein and protein-lipid interactions play crucial roles in the assembly and stability of oligomeric complexes. Recent years have seen a rapid increase in the mechanistic information on the importance of cytoplasmic termini in determining subunit assembly and stability of oligomeric complexes. In addition, the role of specific protein-lipid interaction between anionic phospholipids and "hot spots" on the protein surface has also become evident in stabilizing oligomeric assemblies. This review focuses on several contemporary developments of membrane proteins that stabilize oligomers by taking the potassium channel KcsA as an exemplary ion channel.

Membrane protein dynamics in different environments: simulation study of the outer membrane protein X in a lipid bilayer and in a micelle

The bacterial outer membrane protein OmpX from Escherichia coli has been investigated by molecular dynamics simulations when embedded in a phospholipid bilayer and as a protein-micelle aggregate. The resulting simulation trajectories were analysed in terms of structural and dynamic properties of the membrane protein. In agreement with experimental observations, highest relative stability was found for the β-barrel region that is embedded in the lipophilic phase, whereas an extracellular protruding β-sheet, which is a unique structural feature of OmpX that supposedly plays an important role in cell adhesion and invasion, shows larger structure fluctuations. Additionally, we investigated water permeation into the core of the β-barrel protein, which contains a tight salt-bridge and hydrogen-bond network, so that extensive water flux is unlikely. Differences between the bilayer and the micellar system were observed in the length of the barrel and its position inside the lipid environment, and in the protein interactions with the hydrophilic part of the lipids near the lipid/water interface. Those variations suggest that micelles and other detergent environments might not offer a wholly membrane-like milieu to promote adoption of the physiological conformational state by OmpX.

Characterization of membrane-protein interactions for the leucine transporter from Aquifex aeolicus by molecular dynamics calculations

Multinanosecond molecular dynamics (MD) simulations have been employed to characterize the interaction of an integral membrane protein (IMP), the leucine transmitter from Aquifex aeolicus (Yamashita et al., Nature 2005, 437, 215-223), with hydrated lipid bilayer membranes in their physiologically relevant liquid crystalline phases. Analysis of the MD trajectories for dimyristoyl phosphatidylcholine (DMPC), 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), and 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE) focused on the contacts between aromatic and basic side chains of the IMP with the lipid head groups and water. Structural fluctuations of the IMP were investigated as well as the contact dynamics of neighboring lipids. In characterizing the IMP-membrane systems, the behaviors of the protein's cytoplasmic and periplasmic parts are considered separately. All three lipid membranes show a rather similar overall level of association with the IMP. However, for DMPC there is a better matching of the membrane core to the hydrophobic transmembrane portion of the IMP. The closed cytoplasmic end of the IMP exhibits a higher degree of association with lipids than the more open periplasmic end, an observation which correlates with the more compact structure and a slower dynamics of surrounding lipids.

Identification of specific lipid-binding sites in integral membrane proteins

Protein-lipid interactions are increasingly recognized as central to the structure and function of membrane proteins. However, with the exception of simplified models, specific protein-lipid interactions are particularly difficult to highlight experimentally. Here, we used molecular dynamics simulations to identify a specific protein-lipid interaction in lactose permease, a prototypical model for transmembrane proteins. The interactions can be correlated with the functional dependence of the protein to specific lipid species. The technique is simple and widely applicable to other membrane proteins, and a variety of lipid matrices can be used.

Cause and effect of melittin-induced pore formation: a computational approach

Melittin embedded in a palmitoyl oleyl phosphatidylcholine bilayer at a high peptide/lipid ratio (1:30) was simulated in the presence of explicit water and ions. The simulation results indicate the incipience of an ion-permeable water pore through collective membrane perturbation by bound peptides. The positively charged residues of melittin not only act as "anchors" but also disrupt the membrane, leading to cell lysis. A detailed analysis of the lipid tail order parameter profile depicts localized membrane perturbation. The lipids in the vicinity of the aqueous cavity adopt a tilted conformation, which allows local bilayer thinning. The prepore thus formed can be considered as the melittin-induced structural defects in the bilayer membrane. Because of the strong cationic nature, the melittin-induced prepore exhibits selectivity toward anions over cations. As Cl(-) ions entered into the prepore, they are electrostatically entrapped by positively charged residues located at its wall. The confined motion of the Cl(-) ions in the membrane interior is obvious from calculated diffusion coefficients. Moreover, reorientation of the local lipids occurs in such a way that few lipid heads along with peptide helices can line the surface of the penetrating aqueous phase. The flipping of lipids argued in favor of melittin-induced toroidal pore over a barrel-stave mechanism. Thus, our result provides atomistic level details of the mechanism of membrane disruption by antimicrobial peptide melittin.

Perturbation of membranes by the amyloid beta-peptide--a molecular dynamics study

The etiology of Alzheimer's disease is considered to be linked to interactions between amyloid beta-peptide (Abeta) and neural cell membranes. Membrane disruption and increased ion conductance have been observed in vitro in the presence of Abeta, and it is assumed that these same phenomena occur in the brain of an individual afflicted with Alzheimer's. The effects of Abeta on lipid behavior have been characterized experimentally, but details are lacking regarding how Abeta induces these effects. Simulations of Abeta in a bilayer environment can provide the resolution necessary to explain how the peptide interacts with the surrounding lipids. In the present study, we present an extensive analysis of lipid parameters for a model dipalmitoylphosphatidylcholine bilayer in the presence of the 40-residue Abeta peptide (Abeta40). The simulated systems examine the effects of the insertion depth of the peptide, temperature, the protonation state of the peptide, and ionic strength on the features of the lipid bilayer. The results show that Abeta40 is capable of disordering nearby lipids, as well as decreasing bilayer thickness and area per lipid headgroup. These phenomena arise as a result of the unfolding process of the peptide, which leads to a disordered, extended conformation that is capable of extensive electrostatic and hydrogen-bonding interactions between the peptide and the lipid headgroups. Comparisons are made using melittin-dipalmitoylphosphatidylcholine systems as positive controls of a membrane-disrupting peptide because these systems have previously been characterized experimentally as well as in molecular dynamics simulations.

One membrane protein, two structures and six environments: a comparative molecular dynamics simulation study of the bacterial outer membrane protein PagP

PagP is a bacterial outer membrane protein consisting of an 8 stranded transmembrane beta-barrel and an N-terminal alpha-helix. It is an enzyme which catalyses transfer of a palmitoyl chain from a phospholipid to lipid A. Molecular dynamics simulations have been used to compare the dynamic behaviour in simulations starting from two different structures (X-ray vs. NMR) and in six different environments (detergent micelles formed by dodecyl phosphocholine and by octyl glucoside, vs. four species of phospholipid bilayer). Analysis of interactions between the protein and its environment reveals the role played by the N-terminal alpha-helix, which interacts with the lipid headgroups to lock the PagP molecule into the bilayer. The PagP beta-barrel adopts a tilted orientation in lipid bilayers, facilitating access of lipid tails into the mouth of the central binding pocket. In simulations starting from the X-ray structure in lipid bilayer, the L1 and L2 loops move towards one another, leading to the formation of a putative active site by residues H33, D76 and S77 coming closer together.

Ion channels in microbes

Studies of ion channels have for long been dominated by the animalcentric, if not anthropocentric, view of physiology. The structures and activities of ion channels had, however, evolved long before the appearance of complex multicellular organisms on earth. The diversity of ion channels existing in cellular membranes of prokaryotes is a good example. Although at first it may appear as a paradox that most of what we know about the structure of eukaryotic ion channels is based on the structure of bacterial channels, this should not be surprising given the evolutionary relatedness of all living organisms and suitability of microbial cells for structural studies of biological macromolecules in a laboratory environment. Genome sequences of the human as well as various microbial, plant, and animal organisms unambiguously established the evolutionary links, whereas crystallographic studies of the structures of major types of ion channels published over the last decade clearly demonstrated the advantage of using microbes as experimental organisms. The purpose of this review is not only to provide an account of acquired knowledge on microbial ion channels but also to show that the study of microbes and their ion channels may also hold a key to solving unresolved molecular mysteries in the future.

How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers

We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.

United-atom acyl chains for CHARMM phospholipids

In all-atom simulations of lipid membranes, explicit hydrogen atoms contained in the hydrocarbon region are described by a large number of degrees of freedom, although they convey only limited physical information. We propose an implicit-hydrogen model for saturated and monounsaturated acyl chains, aimed at complementing the all-atom CHARMM27 model for phospholipid headgroups. Torsional potentials and nonbonded parameters were fitted to reproduce experimental data and free energy surfaces of all-atom model systems. Comparative simulations of fluid-phase POPC bilayers were performed using the all-hydrogen force field and the present model. The hybrid model accelerates a typical bilayer simulation by about 50% while sacrificing a minimal amount of detail with respect to the fully atomistic description. In addition, the united-atom description is energetically compatible with all-atom CHARMM models, making it suitable for simulations of complex membrane systems.