Microscopic Characterization of GRP1 PH Domain Interaction with Anionic Membranes

The pleckstrin homology (PH) domain of general receptor for phosphoionositides 1 (GRP1-PHD) binds specifically to phosphatidylinositol (3,4,5)-triphosphate (PIP₃), and acts as a second messenger. Using an extensive array of molecular dynamics (MD) simulations employing highly mobile membrane mimetic (HMMM) model as well as complementary full membrane simulations, we capture differentiable binding and dynamics of GRP1-PHD in the presence of membranes containing PC, PS, and PIP₃ lipids in varying compositions. While GRP1-PHD forms only transient interactions with pure PC membranes, incorporation of anionic lipids resulted in stable membrane-bound configurations. We report the first observation of two distinct PIP₃ binding modes on GRP1-PHD, involving PIP₃ interactions at a “canonical” and at an “alternate” site, suggesting the possibility of simultaneous binding of multiple anionic lipids. The full membrane simulations confirmed the stability of the membrane bound pose of GRP1-PHD as captured from our HMMM membrane binding simulations. By performing additional steered membrane unbinding simulations and calculating nonequilibrium work associated with the process, as well as metadynamics simulations, on the protein bound to full membranes, allowing for more quantitative examination of the binding strength of the GRP1-PHD to the membrane, we demonstrate that along with the bound PIP₃, surrounding anionic PS lipids increase the energetic cost of unbinding of GRP1-PHD from the canonical mode, causing them to dissociate more slowly than the alternate mode. Our results demonstrate that concurrent binding of multiple anionic lipids by GRP1-PHD contributes to its membrane affinity, which in turn control its signaling activity. © 2019 Wiley Periodicals, Inc.     

Investigating the interactions of the first 17 amino acid residues of Huntingtin with lipid vesicles using mass spectrometry and molecular dynamics

The first 17 amino acid residues of Huntingtin protein (Nt17 of htt) are thought to play an important role in the protein's function; Nt17 is one of two membrane binding domains in htt. In this study the binding ability of Nt17 peptide with vesicles comprised of two subclasses of phospholipids is studied using electrospray ionization ‐ mass spectrometry (ESI‐MS) and molecular dynamics (MD) simulations. Overall, the peptide is shown to have a greater propensity to interact with vesicles of phosphatidylcholine (PC) rather than phosphatidylethanolamine (PE) lipids. Mass spectra show an increase in lipid‐bound peptide adducts where the ordering of the number of such specie is 1,2‐dioleoyl‐sn‐glycero‐3‐phosphocholine (DOPC) > 1‐palmitoyl‐2‐oleoyl‐glycero‐3‐phosphocholine (POPC) > 1‐palmitoyl‐2‐oleoyl‐sn‐glycero‐3 phosphoethanolamine (POPE). MD simulations suggest that the compactness of the bilayer plays a role in governing peptide interactions. The peptide shows greater disruption of the DOPC bilayer order at the surface and interacts with the hydrophobic tails of lipid molecules via hydrophobic residues. Conversely, the POPE vesicle remains ordered and lipids display transient interactions with the peptide through the formation of hydrogen bonds with hydrophilic residues. The POPC system displays intermediate behavior with regard to the degree of peptide‐membrane interaction. Finally, the simulations suggest a helix stabilizing effect resulting from the interactions between hydrophobic residues and the lipid tails of the DOPC bilayer.     

A Mass‐Spectrometry‐Based Approach to Distinguish Annular and Specific Lipid Binding to Membrane Proteins

Membrane proteins engage in a variety of contacts with their surrounding lipids, but distinguishing between specifically bound lipids, and non‐specific, annular interactions is a challenging problem. Applying native mass spectrometry to three membrane protein complexes with different lipid‐binding properties, we explore the ability of detergents to compete with lipids bound in different environments. We show that lipids in annular positions on the presenilin homologue protease are subject to constant exchange with detergent. By contrast, detergent‐resistant lipids bound at the dimer interface in the leucine transporter show decreased k_off rates in molecular dynamics simulations. Turning to the lipid flippase MurJ, we find that addition of the natural substrate lipid‐II results in the formation of a 1:1 protein–lipid complex, where the lipid cannot be displaced by detergent from the highly protected active site. In summary, we distinguish annular from non‐annular lipids based on their exchange rates in solution.     

A Mass-Spectrometry-Based Approach to Distinguish Annular and Specific Lipid Binding to Membrane Proteins

Membrane proteins engage in a variety of contacts with their surrounding lipids, but distinguishing between specifically bound lipids, and non-specific, annular interactions is a challenging problem. Applying native mass spectrometry to three membrane protein complexes with different lipid-binding properties, we explore the ability of detergents to compete with lipids bound in different environments. We show that lipids in annular positions on the presenilin homologue protease are subject to constant exchange with detergent. By contrast, detergent-resistant lipids bound at the dimer interface in the leucine transporter show decreased k_off rates in molecular dynamics simulations. Turning to the lipid flippase MurJ, we find that addition of the natural substrate lipid-II results in the formation of a 1:1 protein–lipid complex, where the lipid cannot be displaced by detergent from the highly protected active site. In summary, we distinguish annular from non-annular lipids based on their exchange rates in solution.     

G-Protein-Coupled Receptor–Membrane Interactions Depend on the Receptor Activation State

G-protein-coupled receptors (GPCRs) are the largest family of human membrane proteins and serve as primary targets of approximately one-third of currently marketed drugs. In particular, adenosine A₁ receptor (A₁AR) is an important therapeutic target for treating cardiac ischemia–reperfusion injuries, neuropathic pain, and renal diseases. As a prototypical GPCR, the A₁AR is located within a phospholipid membrane bilayer and transmits cellular signals by changing between different conformational states. It is important to elucidate the lipid–protein interactions in order to understand the functional mechanism of GPCRs. Here, all-atom simulations using a robust Gaussian accelerated molecular dynamics (GaMD) method were performed on both the inactive (antagonist bound) and active (agonist and G-protein bound) A₁AR, which was embedded in a 1-palmitoyl-2-oleoyl-glycero-3-phosphocholine (POPC) lipid bilayer. In the GaMD simulations, the membrane lipids played a key role in stabilizing different conformational states of the A₁AR. Our simulations further identified important regions of the receptor that interacted distinctly with the lipids in highly correlated manner. Activation of the A₁AR led to differential dynamics in the upper and lower leaflets of the lipid bilayer. In summary, GaMD enhanced simulations have revealed strongly coupled dynamics of the GPCR and lipids that depend on the receptor activation state. © 2019 Wiley Periodicals, Inc.     

Delicate Balance of Non-Covalent Forces Govern the Biocompatibility of Graphitic Carbon Nitride towards Genetic Materials

Despite a plethora of suggested technological and biomedical applications, the nanotoxicity of two-dimensional (2D) graphitic carbon nitride (g-C₃N₄) towards biomolecules remains elusive. To address this issue, we employ all-atom classical molecular dynamics simulations and investigate the interactions between nucleic acids and g-C₃N₄. It is revealed that, toxicity is modulated through a subtle balance between electrostatic and van der Waals interactions. When the exposed nucleobases interact through predominantly short-ranged van der Waals and π–π stacking interactions, they get deviated from their native disposition and adsorb on the surface, leading to loss of self-stacking and intra-quartet H-bonding along with partial disruption of the native structure. In contrast, for the interaction with double-stranded structures of both DNA and RNA, long-range electrostatics govern the adsorption phenomena since the constituent nucleobases are relatively concealed and wrapped, thereby resulting in almost complete preservation of the nucleic acid structures. Construction of free energy landscapes for lateral translation of adsorbed nucleic acids suggests decent targeting specificity owing to their restricted movement on g-C₃N₄. The release times of nucleic acids adsorbed through predominant electrostatics are significantly less than those adsorbed through stacking with the surface. It is therefore proposed that g-C₃N₄ would induce toxicity towards any biomolecule having bare residues available for strong van der Waals and π–π stacking interactions relative to those predominantly interacting through electrostatics.     

Modeling,docking and dynamics simulations of a non-specific lipid transfer protein from Peganum harmala L.

Non-specific lipid transfer proteins (ns-LTPs), ubiquitously found in various types of plants, have been well-known to transfer amphiphilic lipids and promote the lipid exchange between mitochondria and microbody. In this study, an in silico analysis was proposed to study ns-LTP in Peganum harmala L., which may belong to ns-LTP1 family, aiming at constructing its three-dimensional structure. Moreover, we adopted MEGA to analyze ns-LTPs and other species phylogenetically, which brought out an initial sequence alignment of ns-LTPs. In addition, we used molecular docking and molecular dynamics simulations to further investigate the affinities and stabilities of ns-LTP with several ligands complexes. Taken together, our results about ns-LTPs and their ligand-binding activities can provide a better understanding of the lipid–protein interactions, indicating some future applications of ns-LTP-mediated transport.     

Effect of PIP2 on Bilayer Structure and Phase Behavior of DOPC: An X‐ray Scattering Study

Phosphatidylinositol 4,5‐bis‐phosphate (PIP₂) is an important lipid in regulation of several cellular processes, particularly membrane fusion. We use X‐ray diffraction from solid‐supported multilamellar 1,2‐dioleoyl‐sn‐glycero‐3‐phosphocholine (DOPC)/PIP₂ samples to study changes in bilayer structure and the lyotropic phase behavior induced by physiologically relevant concentrations of PIP₂. Electron‐density profiles reconstructed from X‐ray reflectivity measurements indicate that PIP₂ strongly affects structural parameters such as lipid head‐group width, bilayer thickness, and lamellar repeat spacing of DOPC bilayer stacks. In addition, at lower degrees of hydration, a few molar per cent of PIP₂ facilitates stalk‐phase formation and also leads to formation of a hexagonal phase, which is not observed in pure DOPC. These results indicate that the role of PIP₂ in membrane fusion could be, in part, due to its effect on the properties of the lipid bilayer matrix. Furthermore, coexistence of two lamellar phases with different lattice constants is observed in single‐component PIP₂ samples.     

Analysis of lipid surface area in protein-membrane systems combining Voronoi tessellation and Monte Carlo integration methods

All-atom molecular dynamics (MD) simulation has become a powerful research tool to investigate structural and dynamical properties of biological membranes and membrane proteins. The lipid structures of simple membrane systems in recent MD simulations are in good agreement with those obtained by experiments. However, for protein-membrane systems, the complexity of protein-lipid interactions makes investigation of lipid structure difficult. Although the area per lipid is one of the essential structural properties in membrane systems, the area in protein-membrane systems cannot be computed easily by conventional approaches like the Voronoi tessellation method. To overcome this limitation, we propose a new method combining the two-dimensional Voronoi tessellation and Monte Carlo integration methods. This approach computes individual surface areas of lipid molecules not only in bulk lipids but also in proximity to membrane proteins. We apply the method to all-atom MD trajectories of the sarcoplasmic reticulum Ca(2+)-pump and the SecY protein-conducting channel. The calculated lipid surface area is in agreement with experimental values and consistent with other structural parameters of lipid bilayers. We also observe changes in the average area per lipid induced by the conformational transition of the SecY channel. Our method is particularly useful for examining equilibration of lipids around membrane proteins and for analyzing the time course of protein-lipid interactions.     

Coacervation of poly-electrolytes in the presence of lipid bilayers: mutual alteration of structure and morphology

Many intrinsically disordered peptides have been shown to undergo liquid–liquid phase separation and form complex coacervates, which play various regulatory roles in the cell. Recent experimental studies found that such phase separation processes may also occur at the lipid membrane surface and help organize biomolecules during signaling events; in some cases, phase separation of proteins at the membrane surface was also observed to lead to significant remodeling of the membrane morphology. The molecular mechanisms that govern the interactions between complex coacervates and lipid membranes and the impacts of such interactions on their structure and morphology, however, remain unclear. Here we study the coacervation of poly-glutamate (E₃₀) and poly-lysine (K₃₀) in the presence of lipid bilayers of different compositions. We carry out explicit-solvent coarse-grained molecular dynamics simulations by using the MARTINI (v3.0) force-field. We find that more than 20% anionic lipids are required for the coacervate to form stable contact with the bilayer. Upon wetting, the coacervate induces negative curvature to the bilayer and facilitates local lipid demixing, without any peptide insertion. The magnitude of negative curvature, extent of lipid demixing, and asphericity of the coacervate increase with the concentration of anionic lipids. Overall, we observe a decrease in the number of contacts among the polyelectrolytes as the droplet spreads over the bilayer. Therefore, unlike previous suggestions, interactions among polyelectrolytes do not constitute a driving force for the membrane bending upon wetting by the coacervate. Rather, analysis of interaction energy components suggests that bending of the membrane is favored by enhanced interactions between polyelectrolytes with lipids as well as with counterions. Kinetic studies reveal that, at the studied polyelectrolyte concentrations, the coacervate formation precedes bilayer wetting.

Intrinsically disordered polyelectrolytes undergoing liquid–liquid phase separation to form complex coacervates on a membrane, which profoundly alters the membrane morphology.     

Lipid/protein interactions and the membrane/water interfacial region

Despite the importance of lipid/protein interactions in the folding, assembly, stability, and function of membrane proteins, information at an atomic level on how such proteins interact with the lipids that surround them remains sparse. The dynamics and flexible nature of the protein/bilayer interaction make it difficult to study, for example, by crystallographic means. However, based on recent progress in molecular simulations of membranes it is possible to address this problem computationally. This communication reports one of the first attempts to use multiple ns molecular simulations to establish a qualitative picture of the intermolecular interactions between the lipids of a bilayer and two topologically different membrane proteins for which a high resolution (2 A or better) X-ray structure is available.     

CHARMM‐GUI Membrane Builder toward realistic biological membrane simulations

CHARMM‐GUI Membrane Builder, http://www.charmm‐gui.org/input/membrane , is a web‐based user interface designed to interactively build all‐atom protein/membrane or membrane‐only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance‐based algorithm for faster initial ion displacement, (5) CHARMM inputs for P₂₁ image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM‐GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments. © 2014 Wiley Periodicals, Inc.     

Atomistic simulations of an antimicrobial molecule interacting with a model bacterial membrane

The dynamics of an antimicrobial molecule (end-only oligo(p-phenylene ethynylene) or EO-OPE-1 (C3)) interacting with a model bacterial membrane is simulated using all-atom molecular dynamics. It is found that the molecule spontaneously adheres to the membrane at the membrane?Cwater interface, but no insertion into the bilayer was observed within the nanosecond simulation time. However, when the simulations start from an inserted configuration, this molecule aligns with the lipid molecules in the membrane and interacts strongly through electrostatic interactions with the anionic phosphoryl groups of the lipid molecules. Due to the hydrophobic mismatch between the molecule and lipids, the inserted molecule induces the deformation of the membrane in the form of local thinning. When more than one molecule were inserted, self-assembling was observed on a nanosecond scale. However, no transmembrane pore formation was observed, due presumably to the hydrophobic backbone of the molecule. Implications in the biocidal action of this molecule are discussed.     

PHYTOCHROME-MEDIATED RAPID CHANGES IN THE LEVEL OF PHOSPHOINOSITIDES IN ETIOLATED LEAVES OF Zea mays

The role of the active form of phytochrome in Zea mays on the polyphosphoinositide cycle was studied. As little as 15 s of red irradiation of etiolated leaves immediately increased the level of phosphatidylinositol bisphosphate (PIP₂) 3–6-fold compared to unirradiated leaves. The elevated level of PIP₂ decreased with longer red irradiations up to 5 min, but remained higher than in unirradiated leaves. The level of PIP₂ decreased if red irradiation was followed by far-red irradiation. Far-red alone had no effect. Levels of phosphatidylinositol phosphate (PIP) and phosphatidylinositol did not change significantly. Since red irradiation significantly changed PIP, but not PIP, photocontrol appears to be at the PIP kinase and phospholipase level. In related studies of the effect of light on phospholipids, 5 min of red irradiation induced significant decreases in phosphatidylcholine and phosphatidylethanola-mine.     

Molecular dynamics study on the stabilization of dehydrated lipid bilayers with glucose and trehalose

In this study, we use molecular dynamics simulations to investigate and compare the interactions of DPPC bilayers with and without saccharides (glucose or trehalose) under dehydrated conditions. Results from the simulations indicate that unilamellar bilayers lose their structural integrity under dehydrated conditions in the absence of saccharides; however, in the presence of either glucose or trehalose, the bilayers maintain their stability. Hydrogen bond analysis shows that the saccharide molecules displace a significant amount of water surrounding the lipid headgroups. At the same time, the additional hydrogen bonds formed between water and saccharide molecules help to maintain a hydration layer on the lipid bilayer interface. On the basis of the hydrogen bond distributions, trehalose forms more hydrogen bonds with the lipids than glucose, and it is less likely to interact with neighboring saccharide molecules. These results suggest that the interaction between the saccharide and lipid molecules through hydrogen bonds is an essential component of the mechanism for the stabilization of lipid bilayers.     

Probing the Lipid Annular Belt by Gas‐Phase Dissociation of Membrane Proteins in Nanodiscs

Interactions between membrane proteins and lipids are often crucial for structure and function yet difficult to define because of their dynamic and heterogeneous nature. Here, we use mass spectrometry to demonstrate that membrane protein oligomers ejected from nanodiscs in the gas phase retain large numbers of lipid interactions. The complex mass spectra that result from gas‐phase dissociation were assigned using a Bayesian deconvolution algorithm together with mass defect analysis, allowing us to count individual lipid molecules bound to membrane proteins. Comparison of the lipid distributions measured by mass spectrometry with molecular dynamics simulations reveals that the distributions correspond to distinct lipid shells that vary according to the type of protein–lipid interactions. Our results demonstrate that nanodiscs offer the potential for native mass spectrometry to probe interactions between membrane proteins and the wider lipid environment.     

Effect of NaCl and KCl on phosphatidylcholine and phosphatidylethanolamine lipid membranes: insight from atomic-scale simulations for understanding salt-induced effects in the plasma membrane

To gain a better understanding of how monovalent salt under physiological conditions affects plasma membranes, we have performed 200 ns atomic-scale molecular dynamics simulations of phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipid bilayers. These two systems provide representative models for the outer and inner leaflets of the plasma membrane, respectively. The implications of cation-lipid interactions in these lipid systems have been considered in two different aqueous salt solutions, namely NaCl and KCl, and the sensitivity of the results on the details of interactions used for ions is determined by repeating the simulations with two distinctly different force fields. We demonstrate that the main effect of monovalent salt on a phospholipid membrane is determined by cations binding to the carbonyl region of a membrane, while chloride anions mostly stay in the water phase. It turns out that the strength and character of the cation-lipid interactions are quite different for different types of lipids and cations. PC membranes and Na+ ions demonstrate strongest interactions, leading to notable membrane compression. This finding was confirmed by both force fields (Gromacs and Charmm) employed for the ions. The binding of potassium ions to PC membranes (and the overall effect of KCl), in turn, was found to be much weaker mainly due to the larger size of a K+ ion compared to Na+. Furthermore, the effect of KCl on PC membranes was found to be force-field sensitive: The binding of a potassium ion was not observed at all in simulations performed with the Gromacs force-field, which seems to exaggerate the size of a K+ ion. As far as PE lipid bilayers are concerned, they are found to be influenced by monovalent salt to a significantly lesser extent compared to PC bilayers, which is a direct consequence of the ability of PE lipids to form both intra- and intermolecular hydrogen bonds and hence to adopt a more densely packed bilayer structure. Whereas for NaCl we observed weak binding of Na+ cations to the PE lipid-water interface, in the case of KCl we witnessed almost complete lack of cation binding. Overall, our findings indicate that monovalent salt ions affect lipids in the inner and outer leaflets of plasma cell membranes in substantially different ways.