Search

Local structural frustration, the existence of mutually exclusive competing interactions, may explain why some proteins are dynamic while others are rigid. Frustration is thought to underpin biomolecular recognition and the flexibility of protein-binding sites. Here, we show how a small chemical modification, the oxidation of two cysteine thiols to a disulfide bond, during the catalytic cycle of the N-terminal domain of the key bacterial oxidoreductase DsbD (nDsbD), introduces frustration ultimately influencing protein function. In oxidized nDsbD, local frustration disrupts the packing of the protective cap-loop region against the active site allowing loop opening. By contrast, in reduced nDsbD the cap loop is rigid, always protecting the active-site thiols from the oxidizing environment of the periplasm. Our results point toward an intricate coupling between the dynamics of the active-site cysteines and of the cap loop which modulates the association reactions of nDsbD with its partners resulting in optimized protein function., Competing Interests: LS, DM, ES, DG, AB, SF, MS, CR No competing interests declared, (© 2020, Stelzl et al.)

Mechanically activated ion channels underlie touch, hearing, shear-stress sensing, and response to turgor pressure. OSCA/TMEM63s are a newly-identified family of eukaryotic mechanically activated ion channels opened by membrane tension. The structural underpinnings of OSCA/TMEM63 function are not explored. Here, we elucidate high resolution cryo-electron microscopy structures of OSCA1.2, revealing a dimeric architecture containing eleven transmembrane helices per subunit and surprising topological similarities to TMEM16 proteins. We locate the ion permeation pathway within each subunit by demonstrating that a conserved acidic residue is a determinant of channel conductance. Molecular dynamics simulations reveal membrane interactions, suggesting the role of lipids in OSCA1.2 gating. These results lay a foundation to decipher how the structural organization of OSCA/TMEM63 is suited for their roles as MA ion channels., Competing Interests: SJ, KS, SM, CT, MS, AP, AW No competing interests declared, (© 2018, Jojoa-Cruz et al.)

In the present contribution, we report a combined spectroscopic and computational approach aiming to unravel at atomic resolution the effect of the anionic SDS detergent on the structure of two model peptides, the α-helix TrpCage and the β-stranded TrpZip. A detailed characterization of the specific amino acids involved is performed. Monomeric (single molecules) and micellar SDS species differently interact with the α-helix and β-stranded peptides, emphasizing the different mechanisms occurring below and above the critical aggregation concentration (CAC). Below the CAC, the α-helix peptide is fully unfolded, losing its hydrophobic core and its Asp-Arg salt bridge, while the β-stranded peptide keeps its native structure with its four Trp well oriented. Above the CAC, the SDS micelles have the same effect on both peptides, that is, destabilizing the tertiary structure while keeping their secondary structure. Our studies will be helpful to deepen our understanding of the action of the denaturant SDS on peptides and proteins., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Recently developed DNA-based analogues of membrane proteins have advanced synthetic biology. A fundamental question is how hydrophilic nanostructures reside in the hydrophobic environment of the membrane. Here, we use multiscale molecular dynamics (MD) simulations to explore the structure, stability and dynamics of an archetypical DNA nanotube inserted via a ring of membrane anchors into a phospholipid bilayer. Coarse-grained MD reveals that the lipids reorganize locally to interact closely with the membrane-spanning section of the DNA tube. Steered simulations along the bilayer normal establish the metastable nature of the inserted pore, yielding a force profile with barriers for membrane exit due to the membrane anchors. Atomistic, equilibrium simulations at two salt concentrations confirm the close packing of lipid around of the stably inserted DNA pore and its cation selectivity, while revealing localized structural fluctuations. The wide-ranging and detailed insight informs the design of next-generation DNA pores for synthetic biology or biomedicine.

Lipases are enzymes of biotechnological importance that function at the interface formed between hydrophobic and aqueous environments. Hydrophobic interfaces can induce structural transitions in lipases that result in an increase in enzyme activity, although the detailed mechanism of this process is currently not well understood for many lipases. Here, we present a multi-scale molecular dynamics simulation study of how different interfaces affect the conformational dynamics of the psychrophilic lipase M37. Our simulations show that M37 lipase is able to interact both with anionic lipid bilayers and with triglyceride surfaces. Interfacial interactions with triglyceride surfaces promote large-scale motions of the lid region of M37, spanning residues 235-283, revealing an entry pathway to the catalytic site for substrates. Importantly, these results suggest a potential activation mechanism for M37 that deviates from other related enzymes, such as Thermomyces lanuginosus lipase. We also investigated substrate binding in M37 by using steered MD simulations, confirming the open state of this lipase. The exposure of hydrophobic residues within lid and active site flap regions (residues 94-110) during the activation process provides insights into the functional effect of hydrophobic surfaces on lipase activation., (Copyright © 2016. Published by Elsevier B.V.)

Supported phospholipid membrane patches stabilized on graphene surfaces have shown potential in sensor device functionalization, including biosensors and biocatalysis. Lipid dip-pen nanolithography (L-DPN) is a method useful in generating supported membrane structures that maintain lipid functionality, such as exhibiting specific interactions with protein molecules. Here, we have integrated L-DPN, atomic force microscopy, and coarse-grained molecular dynamics simulation methods to characterize the molecular properties of supported lipid membranes (SLMs) on graphene and graphene oxide supports. We observed substantial differences in the topologies of the stabilized lipid structures depending on the nature of the surface (polar graphene oxide vs nonpolar graphene). Furthermore, the addition of water to SLM systems resulted in large-scale reorganization of the lipid structures, with measurable effects on lipid lateral mobility within the supported membranes. We also observed reduced lipid ordering within the supported structures relative to free-standing lipid bilayers, attributed to the strong hydrophobic interactions between the lipids and support. Together, our results provide insight into the molecular effects of graphene and graphene oxide surfaces on lipid bilayer membranes. This will be important in the design of these surfaces for applications such as biosensor devices.

It is desirable that nanopores that are components of biosensors are gated, i.e., capable of controllable switching between closed (impermeable) and open (permeable) states. A central hydrophobic barrier within a nanopore may act as a voltage-dependent gate via electrowetting, i.e., changes in nanopore surface wettability by application of an electric field. We use "computational electrophysiology" simulations to demonstrate and characterize electrowetting of a biomimetic nanopore containing a hydrophobic gate. We show that a hydrophobic gate in a model β-barrel nanopore can be functionally opened by electrowetting at voltages that do not electroporate lipid bilayers. During the process of electrowetting, voltage-induced alignment of water dipoles occurs within the hydrophobic gate region of the nanopore, with water entry preceding permeation of ions through the opened nanopore. When the ionic imbalance that generates a transbilayer potential is dissipated, water is expelled from the hydrophobic gate and the nanopore recloses. The open nanopore formed by electrowetting of a "featureless" β-barrel is anionic selective due to the transmembrane dipole potential resulting from binding of Na + ions to the headgroup regions of the surrounding lipid bilayer. Thus, hydrophobic barriers can provide voltage-dependent gates in designed biomimetic nanopores. This extends our understanding of hydrophobic gating in synthetic and biological nanopores, providing a framework for the design of functional nanopores with tailored gating functionality.

Kindlins co-activate integrins alongside talin. They possess, like talin, a FERM domain (4.1-erythrin-radixin-moiesin domain) comprising F0-F3 subdomains, but with a pleckstrin homology (PH) domain inserted in the F2 subdomain that enables membrane association. We present the crystal structure of murine kindlin-3 PH domain determined at a resolution of 2.23 Å and characterise its lipid binding using biophysical and computational approaches. Molecular dynamics simulations suggest flexibility in the PH domain loops connecting β-strands forming the putative phosphatidylinositol phosphate (PtdInsP)-binding site. Simulations with PtdInsP-containing bilayers reveal that the PH domain associates with PtdInsP molecules mainly via the positively charged surface presented by the β1-β2 loop and that it binds with somewhat higher affinity to PtdIns(3,4,5)P 3 compared with PtdIns(4,5)P 2 Surface plasmon resonance (SPR) with lipid headgroups immobilised and the PH domain as an analyte indicate affinities of 300 µM for PtdIns(3,4,5)P 3 and 1 mM for PtdIns(4,5)P 2 In contrast, SPR studies with an immobilised PH domain and lipid nanodiscs as the analyte show affinities of 0.40 µM for PtdIns(3,4,5)P 3 and no affinity for PtdIns(4,5)P 2 when the inositol phosphate constitutes 5% of the total lipids (∼5 molecules per nanodisc). Reducing the PtdIns(3,4,5)P 3 composition to 1% abolishes nanodisc binding to the PH domain, as does site-directed mutagenesis of two lysines within the β1-β2 loop. Binding of PtdIns(3,4,5)P 3 by a canonical PH domain, Grp1, is not similarly influenced by SPR experimental design. These data suggest a role for PtdIns(3,4,5)P 3 clustering in the binding of some PH domains and not others, highlighting the importance of lipid mobility and clustering for the biophysical assessment of protein-membrane interactions., (© 2017 The Author(s).)

The first transmembrane (TM1) helix in the red cell anion exchanger (AE1, Band 3, or SLC4A1) acts as an internal signal anchor that binds the signal recognition particle and directs the nascent polypeptide chain to the endoplasmic reticulum (ER) membrane where it moves from the translocon laterally into the lipid bilayer. The sequence N-terminal to TM1 forms an amphipathic helix that lies at the membrane interface and is connected to TM1 by a bend at Pro403. Southeast Asian ovalocytosis (SAO) is a red cell abnormality caused by a nine-amino acid deletion (Ala400-Ala408) at the N-terminus of TM1. Here we demonstrate, by extensive (∼4.5 μs) molecular dynamics simulations of TM1 in a model 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine membrane, that the isolated TM1 peptide is highly dynamic and samples the structure of TM1 seen in the crystal structure of the membrane domain of AE1. The SAO deletion not only removes the proline-induced bend but also causes a "pulling in" of the part of the amphipathic helix into the hydrophobic phase of the bilayer, as well as the C-terminal of the peptide. The dynamics of the SAO peptide very infrequently resembles the structure of TM1 in AE1, demonstrating the disruptive effect the SAO deletion has on AE1 folding. These results provide a precise molecular view of the disposition and dynamics of wild-type and SAO TM1 in a lipid bilayer, an important early biosynthetic intermediate in the insertion of AE1 into the ER membrane, and extend earlier results of cell-free translation experiments., Competing Interests: The authors declare no competing financial interest.

Pleckstrin homology (PH) domains are membrane-binding lipid recognition proteins that interact with phosphatidylinositol phosphate (PIP) molecules in eukaryotic cell membranes. Diffusion of PH domains plays a critical role in biological reactions on membrane surfaces. Although diffusivity can be estimated by long-time measurements, it lacks information on the short-time diffusive nature. We reveal two diffusive properties of a PH domain bound to the surface of a PIP-containing membrane using molecular dynamics simulations. One is fractional Brownian motion, attributed to the motion of the lipids with which the PH domain interacts. The other is temporally fluctuating diffusivity; that is, the short-time diffusivity of the bound protein changes substantially with time. Moreover, the diffusivity for short-time measurements is intrinsically different from that for long-time measurements. This fluctuating diffusivity results from dynamic changes in interactions between the PH domain and PIP molecules. Our results provide evidence that the complexity of protein-lipid interactions plays a crucial role in the diffusion of proteins on biological membrane surfaces. Changes in the diffusivity of PH domains and related membrane-bound proteins may in turn contribute to the formation/dissolution of protein complexes in membranes.

Neutron and X-ray reflectivity of model membranes is increasingly used as a tool for the study of membrane structures and dynamics. As the systems under study become more complex, and as long, all-atom molecular-dynamics (MD) simulations of membranes become more available, there is increasing interest in the use of MD simulations in the analysis of reflectometry data from membranes. In order to perform this, it is necessary to produce a model of the complete interface, including not only the MD-derived structure of the membrane, but also the supporting substrate and any other interfacial layers that may be present. Here, it is shown that this is best performed by first producing a model of the occupied volume across the entire interface, and then converting this into a scattering length density (SLD) profile, rather than by splicing together the separate SLD profiles from the substrate layers and the membrane, since the latter approach can lead to discontinuities in the SLD profile and subsequent artefacts in the reflectivity calculation. It is also shown how the MD-derived membrane structure should be corrected to account for lower than optimal coverage and out-of-plane membrane fluctuations. Finally, the method of including the entire membrane structure in the reflectivity calculation is compared with an alternative approach in which the membrane components are approximated by functional forms, with only the component volumes being extracted from the simulation. It is shown that using only the fragment volumes is insufficient for a typical neutron data set of a single deuteration measured at several water contrasts, and that either weighting the model by including more structural information from the fit, or a larger data set involving a range of deuterations, are required to satisfactorily define the problem.

The exchange of ADP and ATP across the inner mitochondrial membrane is a fundamental cellular process. This exchange is facilitated by the adenine nucleotide translocase, the structure and function of which are critically dependent on the signature phospholipid of mitochondria, cardiolipin (CL). Here we employ multiscale molecular dynamics simulations to investigate CL interactions within a membrane environment. Using simulations at both coarse-grained and atomistic resolutions, we identify three CL binding sites on the translocase, in agreement with those seen in crystal structures and inferred from nuclear magnetic resonance measurements. Characterization of the free energy landscape for lateral lipid interaction via potential of mean force calculations demonstrates the strength of interaction compared to those of binding sites on other mitochondrial membrane proteins, as well as their selectivity for CL over other phospholipids. Extending the analysis to other members of the family, yeast Aac2p and mouse uncoupling protein 2, suggests a degree of conservation. Simulation of large patches of a model mitochondrial membrane containing multiple copies of the translocase shows that CL interactions persist in the presence of protein-protein interactions and suggests CL may mediate interactions between translocases. This study provides a key example of how computational microscopy may be used to shed light on regulatory lipid-protein interactions.

Molecular dynamics simulations provide a computational tool to probe membrane proteins and systems at length scales ranging from nanometers to close to a micrometer, and on microsecond timescales. All atom and coarse-grained simulations may be used to explore in detail the interactions of membrane proteins and specific lipids, yielding predictions of lipid binding sites in good agreement with available structural data. Building on the success of protein-lipid interaction simulations, larger scale simulations reveal crowding and clustering of proteins, resulting in slow and anomalous diffusional dynamics, within realistic models of cell membranes. Current methods allow near atomic resolution simulations of small membrane organelles, and of enveloped viruses to be performed, revealing key aspects of their structure and functionally important dynamics., (Copyright © 2016 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

The ease with which a cell membrane can bend and deform is important for a wide range of biological functions. Peripheral proteins that induce curvature in membranes (e.g. BAR domains) have been studied for a number of years. Little is known, however, about the effect of integral membrane proteins on the stiffness of a membrane (characterised by the bending rigidity, K c ). We demonstrate by computer simulation that adding integral membrane proteins at physiological densities alters the stiffness of the membrane. First we establish that the coarse-grained MARTINI forcefield is able to accurately reproduce the bending rigidity of a small patch of 1500 phosphatidyl choline lipids by comparing the calculated value to both experiment and an atomistic simulation of the same system. This enables us to simulate the dynamics of large (ca. 50 000 lipids) patches of membrane using the MARTINI coarse-grained description. We find that altering the lipid composition changes the bending rigidity. Adding integral membrane proteins to lipid bilayers also changes the bending rigidity, whilst adding a simple peripheral membrane protein has no effect. Our results suggest that integral membrane proteins can have different effects, and in the case of the bacterial outer membrane protein, BtuB, the greater the density of protein, the larger the reduction in stiffness.

Characterizing the nanoscale dynamic organization within lipid bilayer membranes is central to our understanding of cell membranes at a molecular level. We investigate phase separation and communication across leaflets in ternary lipid bilayers, including saturated lipids with between 12 and 20 carbons per tail. Coarse-grained molecular dynamics simulations reveal a novel two-step kinetics due to hydrophobic mismatch, in which the initial response of the apposed leaflets upon quenching is to increase local asymmetry (antiregistration), followed by dominance of symmetry (registration) as the bilayer equilibrates. Antiregistration can become thermodynamically preferred if domain size is restricted below ∼20 nm, with implications for the symmetry of rafts and nanoclusters in cell membranes, which have similar reported sizes. We relate our findings to theory derived from a semimicroscopic model in which the leaflets experience a "direct" area-dependent coupling, and an "indirect" coupling that arises from hydrophobic mismatch and is most important at domain boundaries. Registered phases differ in composition from antiregistered phases, consistent with a direct coupling between the leaflets. Increased hydrophobic mismatch purifies the phases, suggesting that it contributes to the molecule-level lipid immiscibility. Our results demonstrate an interplay of competing interleaflet couplings that affect phase compositions and kinetics, and lead to a length scale that can influence lateral and transverse bilayer organization within cells.

The cytoskeleton underlying cell membranes may influence the dynamic organization of proteins and lipids within the bilayer by immobilizing certain transmembrane (TM) proteins and forming corrals within the membrane. Here, we present coarse-grained resolution simulations of a biologically realistic membrane model of asymmetrically organized lipids and TM proteins. We determine the effects of a model of cytoskeletal immobilization of selected membrane proteins using long time scale coarse-grained molecular dynamics simulations. By introducing compartments with varying degrees of restraints within the membrane models, we are able to reveal how compartmentalization caused by cytoskeletal immobilization leads to reduced and anomalous diffusional mobility of both proteins and lipids. This in turn results in a reduced rate of protein dimerization within the membrane and of hopping of membrane proteins between compartments. These simulations provide a molecular realization of hierarchical models often invoked to explain single-molecule imaging studies of membrane proteins.

The Hv1 proton channel is evidently unique among voltage sensor domain proteins in mediating an intrinsic 'aqueous' H + conductance (G AQ ). Mutation of a highly conserved 'gating charge' residue in the S4 helix (R1H) confers a resting-state H + 'shuttle' conductance (G SH ) in VGCs and Ci VSP, and we now report that R1H is sufficient to reconstitute G SH in Hv1 without abrogating G AQ . Second-site mutations in S3 (D185A/H) and S4 (N4R) experimentally separate G SH and G AQ gating, which report thermodynamically distinct initial and final steps, respectively, in the Hv1 activation pathway. The effects of Hv1 mutations on G SH and G AQ are used to constrain the positions of key side chains in resting- and activated-state VS model structures, providing new insights into the structural basis of VS activation and H + transfer mechanisms in Hv1., Competing Interests: The authors declare that no competing interests exist.

Lipid molecules can bind to specific sites on integral membrane proteins, modulating their structure and function. We have undertaken coarse-grained simulations to calculate free energy profiles for glycolipids and phospholipids interacting with modulatory sites on the transmembrane helix dimer of the EGF receptor within a lipid bilayer environment. We identify lipid interaction sites at each end of the transmembrane domain and compute interaction free energy profiles for lipids with these sites. Interaction free energies ranged from ca. -40 to -4 kJ/mol for different lipid species. Those lipids (glycolipid GM3 and phosphoinositide PIP2) known to modulate EGFR function exhibit the strongest binding to interaction sites on the EGFR, and we are able to reproduce the preference for interaction with GM3 over other glycolipids suggested by experiment. Mutation of amino acid residues essential for EGFR function reduce the binding free energy of these key lipid species. The residues interacting with the lipids in the simulations are in agreement with those suggested by experimental (mutational) studies. This approach provides a generalizable tool for characterizing the interactions of lipids that bind to specific sites on integral membrane proteins.

The epidermal growth factor receptor (EGFR) is a dimeric membrane protein that regulates key aspects of cellular function. Activation of the EGFR is linked to changes in the conformation of the transmembrane (TM) domain, brought about by changes in interactions of the TM helices of the membrane lipid bilayer. Using an advanced computational approach that combines Coarse-Grained molecular dynamics and well-tempered MetaDynamics (CG-MetaD), we characterize the large-scale motions of the TM helices, simulating multiple association and dissociation events between the helices in membrane, thus leading to a free energy landscape of the dimerization process. The lowest energy state of the TM domain is a right-handed dimer structure in which the TM helices interact through the N-terminal small-X3-small sequence motif. In addition to this state, which is thought to correspond to the active form of the receptor, we have identified further low-energy states that allow us to integrate with a high level of detail a range of previous experimental observations. These conformations may lead to the active state via two possible activation pathways, which involve pivoting and rotational motions of the helices, respectively. Molecular dynamics also reveals correlation between the conformational changes of the TM domains and of the intracellular juxtamembrane domains, paving the way for a comprehensive understanding of EGFR signaling at the cell membrane.

Dok7 is a peripheral membrane protein that is associated with the MuSK receptor tyrosine kinase. Formation of the Dok7/MuSK/membrane complex is required for the activation of MuSK. This is a key step in the complex exchange of signals between neuron and muscle, which lead to neuromuscular junction formation, dysfunction of which is associated with congenital myasthenic syndromes. The Dok7 structure consists of a Pleckstrin Homology (PH) domain and a Phosphotyrosine Binding (PTB) domain. The mechanism of the Dok7 association with the membrane remains largely unknown. Using multi-scale molecular dynamics simulations we have explored the formation of the Dok7 PH/membrane complex. Our simulations indicate that the PH domain of Dok7 associates with membranes containing phosphatidylinositol phosphates (PIPs) via interactions of the β1/β2, β3/β4, and β5/β6 loops, which together form a positively charged surface on the PH domain and interact with the negatively charged headgroups of PIP molecules. The initial encounter of the Dok7 PH domain is followed by formation of additional interactions with the lipid bilayer, and especially with PIP molecules, which stabilizes the Dok7 PH/membrane complex. We have quantified the binding of the PH domain to the model bilayers by calculating a density landscape for protein/membrane interactions. Detailed analysis of the PH/PIP interactions reveal both a canonical and an atypical site to be occupied by the anionic lipid. PH domain binding leads to local clustering of PIP molecules in the bilayer. Association of the Dok7 PH domain with PIP lipids is therefore seen as a key step in localization of Dok7 to the membrane and formation of a complex with MuSK.

The crystal structure of the dimeric membrane domain of human Band 3(1), the red cell chloride/bicarbonate anion exchanger 1 (AE1, SLC4A1), provides a structural context for over four decades of studies into this historic and important membrane glycoprotein. In this review, we highlight the key structural features responsible for anion binding and translocation and have integrated the following topological markers within the Band 3 structure: blood group antigens, N-glycosylation site, protease cleavage sites, inhibitor and chemical labeling sites, and the results of scanning cysteine and N-glycosylation mutagenesis. Locations of mutations linked to human disease, including those responsible for Southeast Asian ovalocytosis, hereditary stomatocytosis, hereditary spherocytosis, and distal renal tubular acidosis, provide molecular insights into their effect on Band 3 folding. Finally, molecular dynamics simulations of phosphatidylcholine self-assembled around Band 3 provide a view of this membrane protein within a lipid bilayer., (Copyright © 2016 Elsevier B.V. All rights reserved.)

Viruses typically pack their genetic material within a protein capsid. Enveloped viruses also have an outer membrane made up of a lipid bilayer and membrane-spanning glycoproteins. X-ray diffraction and cryoelectron microscopy provide high resolution static views of viral structure. Molecular dynamics (MD) simulations may be used to provide dynamic insights into the structures of viruses and their components. There have been a number of simulations of viral capsids and (in some cases) of the inner core of RNA or DNA packaged within them. These simulations have generally focussed on the structural integrity and stability of the capsid and/or on the influence of the nucleic acid core on capsid stability. More recently there have been a number of simulation studies of enveloped viruses, including HIV-1, influenza A, and dengue virus. These have addressed the dynamic behaviour of the capsid, the matrix, and/or of the outer envelope. Analysis of the dynamics of the lipid bilayer components of the envelopes of influenza A and of dengue virus reveals a degree of biophysical robustness, which may contribute to the stability of virus particles in different environments. Significant computational challenges need to be addressed to aid simulation of complex viruses and their membranes, including the need to integrate structural data from a range of sources to enable us to move towards simulations of intact virions. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov., (Copyright © 2016 The Authors. Published by Elsevier B.V. All rights reserved.)

Latrophilin adhesion-GPCRs (Lphn1-3 or ADGRL1-3) and Unc5 cell guidance receptors (Unc5A-D) interact with FLRT proteins (FLRT1-3), thereby promoting cell adhesion and repulsion, respectively. How the three proteins interact and function simultaneously is poorly understood. We show that Unc5D interacts with FLRT2 in cis, controlling cell adhesion in response to externally presented Lphn3. The ectodomains of the three proteins bind cooperatively. Crystal structures of the ternary complex formed by the extracellular domains reveal that Lphn3 dimerizes when bound to FLRT2:Unc5, resulting in a stoichiometry of 1:1:2 (FLRT2:Unc5D:Lphn3). This 1:1:2 complex further dimerizes to form a larger 'super-complex' (2:2:4), using a previously undescribed binding motif in the Unc5D TSP1 domain. Molecular dynamics simulations, point-directed mutagenesis and mass spectrometry demonstrate the stability and molecular properties of these complexes. Our data exemplify how receptors increase their functional repertoire by forming different context-dependent higher-order complexes.

Understanding the energetics of peripheral protein-membrane interactions is important to many areas of biophysical chemistry and cell biology. Estimating free-energy landscapes by molecular dynamics (MD) simulation is challenging for such systems, especially when membrane recognition involves complex lipids, e.g., phosphatidylinositol phosphates (PIPs). We combined coarse-grained MD simulations with umbrella sampling to quantify the binding of the well-explored GRP1 pleckstrin homology (PH) domain to model membranes containing PIP molecules. The experimentally observed preference of GRP1-PH for PIP3 over PIP2 was reproduced. Mutation of a key residue (K273A) within the canonical PIP-binding site significantly reduced the free energy of PIP binding. The presence of a noncanonical PIP-interaction site, observed experimentally in other PH domains but not previously in GRP1-PH, was also revealed. These studies demonstrate how combining coarse-grained simulations and umbrella sampling can unmask the molecular basis of the energetics of interactions between peripheral membrane proteins and complex cellular membranes.

Porins are water-filled protein channels across the outer membrane of gram-negative bacteria. They facilitate the uptake of nutrients and essential ions. Solutes are filtered by a constriction loop L3 at the mid of a pore. Porins are heat-stable and resistant to toxic agents and detergents. Most porins are trimer, but no clear explanation why trimeric form is preferable. In this work, we thus studied effects of oligomerization on porin structure and function in microscopic detail. A well-studied OmpF (general porin from Escherichia coli) and well-characterised OprP (phosphate-specific pore from Pseudomonas aeruginosa) are used as samples from 2 types of porins found in gram-negative bacteria. MD simulations of trimeric and monomeric pores in pure water and 1M NaCl solution were performed. With a salt solution, the external electric field was applied to mimic a transmembrane potential. Expectedly, OprP is more stable than OmpF. Interestingly, being a monomer turns OmpF into an anion-selective pore. The dislocation of D113's side chain on L3 in OmpF causes the disruption of cation pathway resulting in the reduction of cation influx. In contrast, OprP's structure and function are less dependent on oligomeric states. Both monomeric and trimeric OprP can maintain their anion selectivity. Our findings suggest that trimerization is crucial for both structure and function of general porin OmpF, whereas being trimer in substrate-specific channel OprP supports a pore function., (Copyright © 2016 Elsevier Inc. All rights reserved.)

A macromolecular nanopore inserted into a membrane may perturb the dynamic organization of the surrounding lipid bilayer. To better understand the nature of such perturbations, we have undertaken a systematic molecular dynamics simulation study of lipid bilayer structure and dynamics around three different classes of nanopore: a carbon nanotube, three related cyclic peptide nanotubes differing in the nature of their external surfaces, and a model of a β-barrel nanopore protein. Periodic spatial distributions of several lipid properties as a function of distance from the nanopore were observed. This was especially clear for the carbon nanotube system, for which the density of lipids, the bilayer thickness, the projection of lipid head-to-tail vectors onto the membrane plane, and lipid lateral diffusion coefficients exhibited undulatory behavior as a function of the distance from the surface of the channel. Overall, the differences in lipid behavior as a function of the nanopore structure reveal local adaptation of the bilayer structure and dynamics to different embedded nanopore structures. Both the local structure and dynamic behavior of lipids around membrane-embedded nanopores are sensitive to the geometry and nature of the outer surface of the macromolecule/molecular assembly forming the pore.

The dengue virion is surrounded by an envelope of membrane proteins surrounding a lipid bilayer. We have combined the cryoelectron microscopy structures of the membrane proteins (PDB: 3J27) with a lipid bilayer whose composition is based on lipidomics data for insect cell membranes, to obtain a near-atomic resolution computational model of the envelope of the dengue virion. A coarse-grained molecular dynamics simulation on the microsecond timescale enables analysis of key biophysical properties of the dengue outer envelope. Properties analyzed include area per lipid values (for a spherical virion with a mixed lipid composition), bilayer thickness, and lipid diffusion coefficients. Despite the absence of cholesterol from the lipid bilayer, the virion exhibits biophysical robustness (slow lipid diffusion alongside stable bilayer thickness, virion diameter, and shape) that matches the cholesterol-rich membrane of influenza A, with similarly anomalous diffusion of lipids. Biophysical robustness of the envelope may confer resilience to environmental perturbations., (Copyright © 2016 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Ephs are transmembrane receptors that mediate cell-cell signaling. The N-terminal ectodomain binds ligands and enables receptor clustering, which activates the intracellular kinase. Relatively little is known about the function of the membrane-proximal fibronectin domain 2 (FN2) of the ectodomain. Multiscale molecular dynamics simulations reveal that FN2 interacts with lipid bilayers via a site comprising K441, R443, R465, Q462, S464, S491, W467, F490, and P459-461. FN2 preferentially binds anionic lipids, a preference that is reduced in the mutant K441E + R443E. We confirm these results by measuring the binding of wild-type and mutant FN2 domains to lipid vesicles. In simulations of the complete EphA2 ectodomain plus the transmembrane region, we show that FN2 anchors the otherwise flexible ectodomain at the surface of the bilayer. Altogether, our data suggest that FN2 serves a dual function of interacting with anionic lipids and constraining the structure of the EphA2 ectodomain to adopt membrane-proximal configurations., (Copyright © 2016 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Glycolipids are key components of mammalian cell membranes, influencing a diverse range of cellular functions. For example, a number of receptor tyrosine kinases, including the epidermal growth factor receptor (EGFR), are allosterically regulated by the glycolipid monosialodihexosylganglioside (GM3). Recent advances in molecular dynamics methods, especially the development of coarse-grained models, have enabled simulations of increasingly complex models of cell membranes. We demonstrate these methodological developments via a case study of a coarse-grained model for the ganglioside GM3. This glycolipid is included in simulations of a mixed lipid bilayer model reflecting the compositional complexity of a mammalian cell membrane. The resultant membrane model is used to simulate the interactions of GM3 with the transmembrane domain of the EGFR., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

The resemblance of lipid membrane models to physiological membranes determines how well molecular dynamics (MD) simulations imitate the dynamic behavior of cell membranes and membrane proteins. Physiological lipid membranes are composed of multiple types of phospholipids, and the leaflet compositions are generally asymmetric. Here we describe an approach for self-assembly of a Coarse-Grained (CG) membrane model with physiological composition and leaflet asymmetry using the MARTINI force field. An initial set-up of two boxes with different types of lipids according to the leaflet asymmetry of mammalian cell membranes stacked with 0.5 nm overlap, reliably resulted in the self-assembly of bilayer membranes with leaflet asymmetry resembling that of physiological mammalian cell membranes. Self-assembly in the presence of a fragment of the plasma membrane protein syntaxin 1A led to spontaneous specific positioning of phosphatidylionositol(4,5)bisphosphate at a positively charged stretch of syntaxin consistent with experimental data. An analogous approach choosing an initial set-up with two concentric shells filled with different lipid types results in successful assembly of a spherical vesicle with asymmetric leaflet composition. Self-assembly of the vesicle in the presence of the synaptic vesicle protein synaptobrevin 2 revealed the correct position of the synaptobrevin transmembrane domain. This is the first CG MD method to form a membrane with physiological lipid composition as well as leaflet asymmetry by self-assembly and will enable unbiased studies of the incorporation and dynamics of membrane proteins in more realistic CG membrane models.

Pleckstrin homology (PH) domains are lipid-binding modules present in peripheral membrane proteins which interact with phosphatidyl-inositol phosphates (PIPs) in cell membranes. We use multiscale molecular dynamics simulations to characterize the localization and anomalous dynamics of the DAPP1 PH domain on the surface of a PIP-containing lipid bilayer. Both translational and rotational diffusion of the PH domain on the lipid membrane surface exhibit transient subdiffusion, with an exponent α ≈ 0.5 for times of less than 10 ns. In addition to a PIP3 molecule at the canonical binding site of the PH domain, we observe additional PIP molecules in contact with the protein. Fluctuations in the number of PIPs associated with the PH domain exhibit 1/f noise. We suggest that the anomalous diffusion and long-term correlated interaction of the PH domain with the membrane may contribute to an enhanced probability of encounter with target complexes on cell membrane surfaces.

The ability of proteins to sense membrane tension is pervasive in biology. A higher-resolution structure of the Escherichia coli small-conductance mechanosensitive channel MscS identifies alkyl chains inside pockets formed by the transmembrane helices (TMs). Purified MscS contains E. coli lipids, and fluorescence quenching demonstrates that phospholipid acyl chains exchange between bilayer and TM pockets. Molecular dynamics and biophysical analyses show that the volume of the pockets and thus the number of lipid acyl chains within them decreases upon channel opening. Phospholipids with one acyl chain per head group (lysolipids) displace normal phospholipids (with two acyl chains) from MscS pockets and trigger channel opening. We propose that the extent of acyl-chain interdigitation in these pockets determines the conformation of MscS. When interdigitation is perturbed by increased membrane tension or by lysolipids, the closed state becomes unstable, and the channel gates.

The interactions of membrane proteins are influenced by their lipid environment, with key lipid species able to regulate membrane protein function. Advances in high-resolution microscopy can reveal the organization and dynamics of proteins and lipids within living cells at resolutions <200 nm. Parallel advances in molecular simulations provide near-atomic-resolution models of the dynamics of the organization of membranes of in vivo-like complexity. We explore the dynamics of proteins and lipids in crowded and complex plasma membrane models, thereby closing the gap in length and complexity between computations and experiments. Our simulations provide insights into the mutual interplay between lipids and proteins in determining mesoscale (20-100 nm) fluctuations of the bilayer, and in enabling oligomerization and clustering of membrane proteins.

Recently developed synthetic membrane pores composed of folded DNA enrich the current range of natural and engineered protein pores and of nonbiogenic channels. Here we report all-atom molecular dynamics simulations of a DNA nanotube (DNT) pore scaffold to gain fundamental insight into its atomic structure, dynamics, and interactions with ions and water. Our multiple simulations of models of DNTs that are composed of a six-duplex bundle lead to a coherent description. The central tube lumen adopts a cylindrical shape while the mouth regions at the two DNT openings undergo gating-like motions which provide a possible molecular explanation of a lower conductance state observed in our previous experimental study on a membrane-spanning version of the DNT (ACS Nano 2015, 9, 1117-26). Similarly, the central nanotube lumen is filled with water and ions characterized by bulk diffusion coefficients while the gating regions exhibit temporal fluctuations in their aqueous volume. We furthermore observe that the porous nature of the walls allows lateral leakage of ions and water. This study will benefit rational design of DNA nanopores of enhanced stability of relevance for sensing applications, of nanodevices with tunable gating properties that mimic gated ion channels, or of nanopores featuring defined permeation behavior.

Gram-negative bacteria inhabit a broad range of ecological niches. For Escherichia coli, this includes river water as well as humans and animals, where it can be both a commensal and a pathogen. Intricate regulatory mechanisms ensure that bacteria have the right complement of β-barrel outer membrane proteins (OMPs) to enable adaptation to a particular habitat. Yet no mechanism is known for replacing OMPs in the outer membrane, an issue that is further confounded by the lack of an energy source and the high stability and abundance of OMPs. Here we uncover the process underpinning OMP turnover in E. coli and show it to be passive and binary in nature, in which old OMPs are displaced to the poles of growing cells as new OMPs take their place. Using fluorescent colicins as OMP-specific probes, in combination with ensemble and single-molecule fluorescence microscopy in vivo and in vitro, as well as molecular dynamics simulations, we established the mechanism for binary OMP partitioning. OMPs clustered to form ∼0.5-μm diameter islands, where their diffusion is restricted by promiscuous interactions with other OMPs. OMP islands were distributed throughout the cell and contained the Bam complex, which catalyses the insertion of OMPs in the outer membrane. However, OMP biogenesis occurred as a gradient that was highest at mid-cell but largely absent at cell poles. The cumulative effect is to push old OMP islands towards the poles of growing cells, leading to a binary distribution when cells divide. Hence, the outer membrane of a Gram-negative bacterium is a spatially and temporally organized structure, and this organization lies at the heart of how OMPs are turned over in the membrane.

A gain-of-function mutation (T635A) in the transient receptor potential (TRP) channel TRPC3 results in abnormal channel gating and causes cerebellar ataxia in the dominant Moonwalker (Mwk) mouse mutant. However, the underlying molecular and structural mechanisms are unclear. Here, we used a combined approach of computational modeling and functional characterization of proposed TRPC3 mutants. Our findings support a mechanism by which the hydrogen bonding capability of threonine 635 plays a significant role in maintaining a stable, closed state channel. This capability is lost in the Mwk mutant, suggesting a structural basis for the disease-causing phenotype in the Mwk mouse.

There has been exponential growth in the number of membrane protein structures determined. Nevertheless, these structures are usually resolved in the absence of their lipid environment. Coarse-grained molecular dynamics (CGMD) simulations enable insertion of membrane proteins into explicit models of lipid bilayers. We have automated the CGMD methodology, enabling membrane protein structures to be identified upon their release into the PDB and embedded into a membrane. The simulations are analyzed for protein-lipid interactions, identifying lipid binding sites, and revealing local bilayer deformations plus molecular access pathways within the membrane. The coarse-grained models of membrane protein/bilayer complexes are transformed to atomistic resolution for further analysis and simulation. Using this automated simulation pipeline, we have analyzed a number of recently determined membrane protein structures to predict their locations within a membrane, their lipid/protein interactions, and the functional implications of an enhanced understanding of the local membrane environment of each protein., (Copyright © 2015 The Authors. Published by Elsevier Ltd.. All rights reserved.)

A necessary step prior to starting any membrane protein computer simulation is the creation of a well-packed configuration of protein(s) and lipids. Here, we demonstrate a method, alchembed , that can simultaneously and rapidly embed multiple proteins into arrangements of lipids described using either atomistic or coarse-grained force fields. During a short simulation, the interactions between the protein(s) and lipids are gradually switched on using a soft-core van der Waals potential. We validate the method on a range of membrane proteins and determine the optimal soft-core parameters required to insert membrane proteins. Since all of the major biomolecular codes include soft-core van der Waals potentials, no additional code is required to apply this method. A tutorial is included in the Supporting Information.