Transient receptor potential vanilloid subtype 1 (TRPV1) is a heat-sensitive ion channel also involved in pain sensation, and is the receptor for capsaicin, the active ingredient of hot chili peppers. The recent structures of TRPV1 revealed putative ligand density within the S1 to S4 voltage-sensor-like domain of the protein. However, questions remain regarding the dynamic role of the lipid bilayer in ligand binding to TRPV1. Molecular dynamics simulations were used to explore behavior of capsaicin in a 1-palmitoyl-2-oleoyl phosphatidylcholine bilayer and with the target S1-S4 transmembrane helices of TRPV1. Equilibrium simulations reveal a preferred interfacial localization for capsaicin. We also observed a capsaicin molecule flipping from the extracellular to the intracellular leaflet, and subsequently able to access the intracellular TRPV1 binding site. Calculation of the potential of mean force (i.e., free energy profile) of capsaicin along the bilayer normal confirms that it prefers an interfacial localization. The free energy profile indicates that there is a nontrivial but surmountable barrier to the flipping of capsaicin between opposing leaflets of the bilayer. Molecular dynamics of the S1-S4 transmembrane helices of the TRPV1 in a lipid bilayer confirm that Y511, known to be crucial to capsaicin binding, has a distribution along the bilayer normal similar to that of the aromatic group of capsaicin. Simulations were conducted of the TRPV1 S1-S4 transmembrane helices in the presence of capsaicin placed in the aqueous phase, in the lipid, or docked to the protein. No stable interaction between ligand and protein was seen for simulations initiated with capsaicin in the bilayer. However, interactions were seen between TRPV1 and capsaicin starting from the cytosolic aqueous phase, and capsaicin remained stable in the majority of simulations from the docked pose. We discuss the significance of capsaicin flipping from the extracellular to the intracellular leaflet and mechanisms of binding site access by capsaicin.

Structural studies of membrane proteins have highlighted the likely influence of membrane mimetic environments (i.e., lipid bilayers versus detergent micelles) on the conformation and dynamics of small α-helical membrane proteins. We have used molecular dynamics simulations to compare the conformational dynamics of BM2 (a small α-helical protein from the membrane of influenza B) in a model phospholipid bilayer environment with its behavior in protein-detergent complexes with either the zwitterionic detergent dihexanoylphosphatidylcholine (DHPC) or the nonionic detergent dodecylmaltoside (DDM). We find that DDM more closely resembles the lipid bilayer in terms of its interaction with the protein, while the short-tailed DHPC molecule forms "nonphysiological" interactions with the protein termini. We find that the intrinsic micelle properties of each detergent are conserved upon formation of the protein-detergent complex. This implies that simulations of detergent micelles may be used to help select optimal conditions for experimental studies of membrane proteins.

Multiple nanosecond duration molecular dynamics simulations were performed on the transmembrane region of the Torpedo nicotinic acetylcholine receptor embedded within a bilayer mimetic octane slab. The M2 helices and M2-M3 loop regions were free to move, whereas the outer (M1, M3, M4) helix bundle was backbone restrained. The M2 helices largely retain their hydrogen-bonding pattern throughout the simulation, with some distortions in the helical end and loop regions. All of the M2 helices exhibit bending motions, with the hinge point in the vicinity of the central hydrophobic gate region (corresponding to residues alphaL251 and alphaV255). The bending motions of the M2 helices lead to a degree of dynamic narrowing of the pore in the region of the proposed hydrophobic gate. Calculations of Born energy profiles for various structures along the simulation trajectory suggest that the conformations of the M2 bundle sampled correspond to a closed conformation of the channel. Principal components analyses of each of the M2 helices, and of the five-helix M2 bundle, reveal concerted motions that may be relevant to channel function. Normal mode analyses using the anisotropic network model reveal collective motions similar to those identified by principal components analyses.

Molecular simulations are an invaluable tool for understanding membrane proteins. Improvements to both hardware and simulation methods have allowed access to physiologically relevant timescales and have permitted the simulation of large multimeric complexes. This, coupled to the recent expansion in membrane protein structures, provides a means to elucidate the relationship between protein structure and function. In this review, we discuss the progress in using simulations to understand the complex processes that occur at the boundary of a cell, ranging from the transport of solutes and the interactions of ligands with ion channels to the conformational rearrangements required for gating of channels and the signaling by membrane-associated complexes.

Ion channels play key roles in cell membranes, and recent advances are yielding an increasing number of structures. However, their functional relevance is often unclear and better tools are required for their functional annotation. In sub-nanometer pores such as ion channels, hydrophobic gating has been shown to promote dewetting to produce a functionally closed (i.e., non-conductive) state. Using the serotonin receptor (5-HT3R) structure as an example, we demonstrate the use of molecular dynamics to aid the functional annotation of channel structures via simulation of the behavior of water within the pore. Three increasingly complex simulation analyses are described: water equilibrium densities; single-ion free-energy profiles; and computational electrophysiology. All three approaches correctly predict the 5-HT3R crystal structure to represent a functionally closed (i.e., non-conductive) state. We also illustrate the application of water equilibrium density simulations to annotate different conformational states of a glycine receptor.

Specific interactions of lipids with membrane proteins contribute to protein stability and function. Multiple lipid interactions surrounding a membrane protein are often identified in molecular dynamics (MD) simulations and are, increasingly, resolved in cryo-electron microscopy (cryo-EM) densities. Determining the relative importance of specific interaction sites is aided by determination of lipid binding affinities using experimental or simulation methods. Here, we develop a method for determining protein-lipid binding affinities from equilibrium coarse-grained MD simulations using binding saturation curves, designed to mimic experimental protocols. We apply this method to directly obtain affinities for cholesterol binding to multiple sites on a range of membrane proteins and compare our results with free energies obtained from density-based equilibrium methods and with potential of mean force calculations, getting good agreement with respect to the ranking of affinities for different sites. Thus, our binding saturation method provides a robust, high-throughput alternative for determining the relative consequence of individual sites seen in, e.g., cryo-EM derived membrane protein structures surrounded by an array of ancillary lipid densities.

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (alpha-helical and beta-barrel), and the seven different bilayer systems range in thickness from approximately 28 to approximately 43 A. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.

Drug resistance in gram-negative bacteria may be conferred via efflux through a tripartite complex of an inner membrane pump, an outer membrane pore, and a periplasmic adaptor protein. These are AcrB, TolC, and AcrA, respectively, in Escherichia coli. In Pseudomonas aerugonisa, their homologs are MexB, OprM, and MexA. Defining the interdomain dynamics of the adaptor protein is essential to understanding the mechanism of complex formation. Extended (25 ns) molecular dynamics simulations of MexA have been performed to determine such interdomain dynamics. Analysis of conformational drift demonstrates substantial motions of the three domains of MexA relative to one another. Principal components analysis reveals a hinge-bending motion and rotation of the alpha-helical hairpin relative to the other domains to be the two dominant motions. These two motions provide an element of considerable flexibility which is likely to be exploited in the adaptor function of MexA.

Integrins are heterodimeric (αβ) cell surface receptors that are potential therapeutic targets for a number of diseases. Despite the existence of structural data for all parts of integrins, the structure of the complete integrin receptor is still not available. We have used available structural data to construct a model of the complete integrin receptor in complex with talin F2-F3 domain. It has been shown that the interactions of integrins with their lipid environment are crucial for their function but details of the integrin/lipid interactions remain elusive. In this study an integrin/talin complex was inserted in biologically relevant bilayers that resemble the cell plasma membrane containing zwitterionic and charged phospholipids, cholesterol and sphingolipids to study the dynamics of the integrin receptor and its effect on bilayer structure and dynamics. The results of this study demonstrate the dynamic nature of the integrin receptor and suggest that the presence of the integrin receptor alters the lipid organization between the two leaflets of the bilayer. In particular, our results suggest elevated density of cholesterol and of phosphatidylserine lipids around the integrin/talin complex and a slowing down of lipids in an annulus of ~30 Å around the protein due to interactions between the lipids and the integrin/talin F2-F3 complex. This may in part regulate the interactions of integrins with other related proteins or integrin clustering thus facilitating signal transduction across cell membranes.

ATP-binding cassette (ABC) transporters mediate the movement of molecules across cell membranes in both prokaryotes and eukaryotes. In ABC transporters, solute translocation occurs after ATP is either bound or hydrolyzed at the intracellular nucleotide-binding domains (NBDs). Molecular dynamics (MD) simulations have been employed to study the interactions of nucleotide with NBD. The results of extended (approximately 20 ns) MD simulations of HisP (total simulation time approximately 80 ns), the NBD of the histidine transporter HisQMP2J from Salmonella typhimurium, are presented. Analysis of the MD trajectories reveals conformational changes within HisP that are dependent on the presence of ATP in the binding pocket of the protein, and are sensitive to the presence/absence of Mg ions bound to the ATP. These changes are predominantly confined to the alpha-helical subdomain of HisP. Specifically there is a rotation of three alpha-helices within the subdomain, and a movement of the signature sequence toward the bound nucleotide. In addition, there is considerable conformational flexibility in a conserved glutamine-containing loop, which is situated at the interface between the alpha-helical subdomain and the F1-like subdomain. These results support the mechanism for ATP-induced conformational transitions derived from the crystal structures of other NBDs.

C2 domains facilitate protein interactions with lipid bilayers in either a Ca2+-dependent or -independent manner. We used molecular dynamics (MD) simulations to explore six Ca2+-independent C2 domains, from KIBRA, PI3KC2α, RIM2, PTEN, SHIP2, and Smurf2. In coarse-grained MD simulations these C2 domains formed transient interactions with zwitterionic bilayers, compared with longer-lived interactions with anionic bilayers containing phosphatidylinositol bisphosphate (PIP2). Type I C2 domains bound non-canonically via the front, back, or side of the β sandwich, whereas type II C2 domains bound canonically, via the top loops. C2 domains interacted strongly with membranes containing PIP2, causing bound anionic lipids to cluster around the protein. Binding modes were refined via atomistic simulations. For PTEN and SHIP2, CG simulations of their phosphatase plus C2 domains with PIP2-containing bilayers were also performed, and the roles of the two domains in membrane localization compared. These studies establish a simulation protocol for membrane-recognition proteins.

A molecular simulation pipeline for determining the mode of interaction of pleckstrin homology (PH) domains with phosphatidylinositol phosphate (PIP)-containing lipid bilayers is presented. We evaluate our methodology for the GRP1 PH domain via comparison with structural and biophysical data. Coarse-grained simulations yield a 2D density landscape for PH/membrane interactions alongside residue contact profiles. Predictions of the membrane localization and interactions of 13 PH domains reveal canonical, non-canonical, and dual PIP-binding sites on the proteins. Thus, the PH domains associate with the PIP molecules in the membrane via a highly positively charged loop. Some PH domains exhibit modes of interaction with PIP-containing membranes additional to this canonical binding mode. All 13 PH domains cause a degree of local clustering of PIP molecules upon binding to the membrane. This provides a global picture of PH domain interactions with membranes. The high-throughput approach could be extended to other families of peripheral membrane proteins.

α-Synuclein (αS) is a presynaptic protein that binds to cell membranes and is linked to Parkinson's disease (PD). Binding of αS to membranes is a likely first step in the molecular pathophysiology of PD. The αS molecule can adopt multiple conformations, being largely disordered in water, adopting a β-sheet conformation when present in amyloid fibrils, and forming a dynamic multiplicity of α-helical conformations when bound to lipid bilayers and related membrane-mimetic surfaces. Multiscale molecular dynamics simulations in conjunction with nuclear magnetic resonance (NMR) and cross-linking mass spectrometry (XLMS) measurements are used to explore the interactions of αS with an anionic lipid bilayer. The simulations and NMR measurements together reveal a break in the helical structure of the central non-amyloid-β component (NAC) region of αS in the vicinity of residues 65-70, which may facilitate subsequent oligomer formation. Coarse-grained simulations of αS starting from the structure of αS when bound to a detergent micelle reveal the overall pattern of protein contacts to anionic lipid bilayers, while subsequent all-atom simulations provide details of conformational changes upon membrane binding. In particular, simulations and NMR data for liposome-bound αS indicate incipient β-strand formation in the NAC region, which is supported by intramolecular contacts seen via XLMS and simulations. Markov state models based on the all-atom simulations suggest a mechanism of conformational change of membrane-bound αS via a dynamic helix break in the region of residue 65 in the NAC region. The emergent dynamic model of membrane-interacting αS advances our understanding of the mechanism of PD, potentially aiding the design of novel therapeutic approaches.

Water molecules within biological ion channels are in a nanoconfined environment and therefore exhibit behaviors which differ from that of bulk water. Here, we investigate the phenomenon of hydrophobic gating, the process by which a nanopore may spontaneously dewet to form a "vapor lock" if the pore is sufficiently hydrophobic and/or narrow. This occurs without steric occlusion of the pore. Using molecular dynamics simulations with both rigid fixed-charge and polarizable (AMOEBA) force fields, we investigate this wetting/dewetting behavior in the transmembrane protein 175 ion channel. We examine how a range of rigid fixed-charge and polarizable water models affect wetting/dewetting in both the wild-type structure and in mutants chosen to cover a range of nanopore radii and pore-lining hydrophobicities. Crucially, we find that the rigid fixed-charge water models lead to similar wetting/dewetting behaviors, but that the polarizable water model resulted in an increased wettability of the hydrophobic gating region of the pore. This has significant implications for molecular simulations of nanoconfined water, as it implies that polarizability may need to be included if we are to gain detailed mechanistic insights into wetting/dewetting processes. These findings are of importance for the design of functionalized biomimetic nanopores (e.g., sensing or desalination) as well as for furthering our understanding of the mechanistic processes underlying biological ion channel function.

FhuA is one of the more complex members of the superfamily of bacterial outer membrane proteins. Its primary function is to provide a binding site on the outer membrane surface for siderophores, such as ferrichrome, and subsequently to facilitate their energy-dependent transport across the membrane, presumably powered by the TonB-ExbBD protein complex that resides in the cytoplasmic membrane. Crystal structures of FhuA with and without a bound ferrichrome molecule have provided some clues as to the initial stages of the siderophore transport mechanism. In the current study, we have employed 10-ns duration molecular dynamics simulations of FhuA and of the FhuA-ferrichrome complex, both embedded in a phospholipid bilayer, to probe the short timescale dynamics of this integral membrane protein, and to explore possible mechanistic implications of this dynamic behavior. Analysis of the dynamics of the protein suggests that the extracellular loops move as relatively rigid entities relative to the transmembrane beta-barrel. Comparison of the two simulations (with and without bound ferrichrome) revealed some ligand-induced changes in loop mobility. Specifically, loop L8 appears to be involved in a mechanism whereby the binding site is gated closed upon ligand binding. Analysis of the dynamics of water molecules within the core of the FhuA protein provided no evidence for a water-permeable protopore through which the ferrichrome might pass without a major perturbation of the FhuA protein. Overall, these simulations support the proposal that binding of ferrichrome initiates a signaling mechanism that ultimately leads to the TonB-mediated partial or total removal of the core domain from the beta-barrel, thus opening up a permeable pore. These simulations are among the longest that have been performed on a complex membrane protein. However, a simple analysis of sampling reveals that the description of protein motions is far from complete.

ABC transporters are a large family of membrane proteins involved in a variety of cellular processes, including multidrug and tumor resistance and ion channel regulation. Advances in the structural and functional understanding of ABC transporters have revealed that hydrolysis at the two canonical nucleotide-binding sites (NBSs) is co-operative and non-simultaneous. A conserved core architecture of bacterial and eukaryotic ABC exporters has been established, as exemplified by the crystal structure of the homodimeric multidrug exporter Sav1866. Currently, it is unclear how sequential ATP hydrolysis arises in a symmetric homodimeric transporter, since it implies at least transient asymmetry at the NBSs. We show by molecular dynamics simulation that the initially symmetric structure of Sav1866 readily undergoes asymmetric transitions at its NBSs in a pre-hydrolytic nucleotide configuration. MgATP-binding residues and a network of charged residues at the dimer interface are shown to form a sequence of putative molecular switches that allow ATP hydrolysis only at one NBS. We extend our findings to eukaryotic ABC exporters which often consist of two non-identical half-transporters, frequently with degeneracy substitutions at one of their two NBSs. Interestingly, many residues involved in asymmetric conformational switching in Sav1866 are substituted in degenerate eukaryotic NBS. This finding strengthens recent suggestions that the interplay of a consensus and a degenerate NBS in eukaroytic ABC proteins pre-determines the sequence of hydrolysis at the two NBSs.

Ion channel-forming peptides enable us to study the conformational dynamics of a transmembrane helix as a function of sequence and environment. Molecular dynamics simulations are used to study the conformation and dynamics of three 22-residue peptides derived from the second transmembrane domain of the glycine receptor (NK4-M2GlyR-p22). Simulations are performed on the peptide in four different environments: trifluoroethanol/water; SDS micelles; DPC micelles; and a DMPC bilayer. A hierarchy of alpha-helix stabilization between the different environments is observed such that TFE/water < micelles < bilayers. Local clustering of trifluoroethanol molecules around the peptide appears to help stabilize an alpha-helical conformation. Single (S22W) and double (S22W,T19R) substitutions at the C-terminus of NK4-M2GlyR-p22 help to stabilize a helical conformation in the micelle and bilayer environments. This correlates with the ability of the W22 and R19 side chains to form H-bonds with the headgroups of lipid or detergent molecules. This study provides a first atomic resolution comparison of the structure and dynamics of NK4-M2GlyR-p22 peptides in membrane and membrane-mimetic environments, paralleling NMR and functional studies of these peptides.

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG >> PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.

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.

The bacterial outer membrane protein OmpA is composed of an N-terminal 171-residue beta-barrel domain (OmpA(171)) that spans the bilayer and a periplasmic, C-terminal domain of unknown structure. OmpA has been suggested to primarily serve a structural role, as no continuous pore through the center of the barrel can be discerned in the crystal structure of OmpA(171). However, several groups have recorded ionic conductances for bilayer-reconstituted OmpA(171). To resolve this apparent paradox we have used molecular dynamics (MD) simulations on OmpA(171) to explore the conformational dynamics of the protein, in particular the possibility of transient formation of a central pore. A total of 19 ns of MD simulations of OmpA(171) have been run, and the results were analyzed in terms of 1) comparative behavior of OmpA(171) in different bilayer and bilayer-mimetic environments, 2) solvation states of OmpA(171), and 3) pore characteristics in different MD simulations. Significant mobility was observed for residues and water molecules within the beta-barrel. A simulation in which putative gate region side chains of the barrel interior were held in a non-native conformation led to an open pore, with a predicted conductance similar to experimental measurements. The OmpA(171) pore has been shown to be somewhat more dynamic than suggested by the crystal structure. A gating mechanism is proposed to explain its documented channel properties, involving a flickering isomerization of Arg138, forming alternate salt bridges with Glu52 (closed state) and Glu128 (open state).