Search

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.

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.

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.

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.

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.

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.

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 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.

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.

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.

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.

Nanopores in membranes have a range of potential applications. Biomimetic design of nanopores aims to mimic key functions of biological pores within a stable template structure. Molecular dynamics simulations have been used to test whether a simple β-barrel protein nanopore can be modified to incorporate a hydrophobic barrier to permeation. Simulations have been used to evaluate functional properties of such nanopores, using water flux as a proxy for ionic conductance. The behavior of these model pores has been characterized as a function of pore size and of the hydrophobicity of the amino acid side chains lining the narrow central constriction of the pore. Potential of mean force calculations have been used to calculate free energy landscapes for water and for ion permeation in selected models. These studies demonstrate that a hydrophobic barrier can indeed be designed into a β-barrel protein nanopore, and that the height of the barrier can be adjusted by modifying the number of consecutive rings of hydrophobic side chains. A hydrophobic barrier prevents both water and ion permeation even though the pore is sterically unoccluded. These results both provide insights into the nature of hydrophobic gating in biological pores and channels, and furthermore demonstrate that simple design features may be computationally transplanted into β-barrel membrane proteins to generate functionally complex nanopores.

EphA1 is a receptor tyrosine kinase (RTK) that plays a key role in developmental processes, including guidance of the migration of axons and cells in the nervous system. EphA1, in common with other RTKs, contains an N-terminal extracellular domain, a single transmembrane (TM) α-helix, and a C-terminal intracellular kinase domain. The TM helix forms a dimer, as seen in recent NMR studies. We have modeled the EphA1 TM dimer using a multiscale approach combining coarse-grain (CG) and atomistic molecular dynamics (MD) simulations. The one-dimensional potential of mean force (PMF) for this system, based on interhelix separation, has been calculated using CG MD simulations. This provides a view of the free energy landscape for helix-helix interactions of the TM dimer in a lipid bilayer. The resulting PMF profiles suggest two states, consistent with a rotation-coupled activation mechanism. The more stable state corresponds to a right-handed helix dimer interacting via an N-terminal glycine zipper motif, consistent with a recent NMR structure (2K1K). A second metastable state corresponds to a structure in which the glycine zipper motif is not involved. Analysis of unrestrained CG MD simulations based on representative models from the PMF calculations or on the NMR structure reveals possible pathways of interconversion between these two states, involving helix rotations about their long axes. This suggests that the interaction of TM helices in EphA1 dimers may be intrinsically dynamic. This provides a potential mechanism for signaling whereby extracellular events drive a shift in the repopulation of the underlying TM helix dimer energy landscape.

The interactions of transmembrane (TM) α-helices with the phospholipid membrane and with one another are central to understanding the structure and stability of integral membrane proteins. These interactions may be analyzed via coarse grained molecular dynamics (CGMD) simulations. To obtain statistically meaningful analysis of TM helix interactions, large (N ca. 100) ensembles of CGMD simulations are needed. To facilitate the running and analysis of such ensembles of simulations, we have developed Sidekick, an automated pipeline software for performing high throughput CGMD simulations of α-helical peptides in lipid bilayer membranes. Through an end-to-end approach, which takes as input a helix sequence and outputs analytical metrics derived from CGMD simulations, we are able to predict the orientation and likelihood of insertion into a lipid bilayer of a given helix of a family of helix sequences. We illustrate this software via analyses of insertion into a membrane of short hydrophobic TM helices containing a single cationic arginine residue positioned at different positions along the length of the helix. From analyses of these ensembles of simulations, we estimate apparent energy barriers to insertion which are comparable to experimentally determined values. In a second application, we use CGMD simulations to examine the self-assembly of dimers of TM helices from the ErbB1 receptor tyrosine kinase and analyze the numbers of simulation repeats necessary to obtain convergence of simple descriptors of the mode of packing of the two helices within a dimer. Our approach offers a proof-of-principle platform for the further employment of automation in large ensemble CGMD simulations of membrane proteins.

The phosphatase and tensin homologue (PTEN) and the Ciona intestinalis voltage sensitive phosphatase (Ci-VSP) are both phosphatidylinositol phosphate (PIP) phosphatases that contain a C2 domain. PTEN is a tumor suppressor protein that acts as a phosphatase on PIP3 in mammalian cell membranes. It contains two principal domains: a phosphatase domain (PD) and a C2 domain. Despite detailed structural and functional characterization, less is known about its mechanism of interaction with PIP-containing lipid bilayers. Ci-VSP consists of an N-terminal transmembrane voltage sensor domain and a C-terminal PTEN domain, which in turn contains a PD and a C2 domain. The nature of the interaction of the PTEN domain of Ci-VSP with membranes has not been well established. We have used multiscale molecular dynamics simulations to define the interaction mechanisms of PTEN and of the Ci-VSP PTEN domains with PIP-containing lipid bilayers. Our results suggest a novel mechanism of association of the PTEN with such bilayers, in which an initial electrostatics-driven encounter of the protein and bilayer is followed by reorientation of the protein to optimize its interactions with PIP molecules in the membrane. Although a PIP3 molecule binds close to the active site of PTEN, our simulations suggest a further conformational change of the protein may be required for catalytically productive binding to occur. Ci-VSP interacted with membranes in an orientation comparable to that of PTEN but bound directly to PIP-containing membranes without a subsequent reorientation step. Again, PIP3 bound close to the active site of the Ci-VSP PD, but not in a catalytically productive manner. Interactions of Ci-VSP with the bilayer induced clustering of PIP molecules around the protein.

Receptor tyrosine kinases are single-pass membrane proteins that form dimers within the membrane. The interactions of their transmembrane domains (TMDs) play a key role in dimerization and signaling. Fibroblast growth factor receptor 3 (FGFR3) is of interest as a G380R mutation in its TMD is the underlying cause of ~99% of the cases of achondroplasia, the most common form of human dwarfism. The structural consequences of this mutation remain uncertain: the mutation shifts the position of the TMD relative to the lipid bilayer but does not alter the association free energy. We have combined coarse-grained and all-atom molecular dynamics simulations to study the dimerization of wild-type, heterodimer, and mutant FGFR3 TMDs. The simulations reveal that the helices pack together in the dimer to form a flexible interface. The primary packing mode is mediated by a Gx3G motif. There is also a secondary dimer interface that is more highly populated in heterodimer and mutant configurations that may feature in the molecular mechanism of pathology. Both coarse-grained and atomistic simulations reveal a significant shift of the G380R mutant dimer TMD relative to the bilayer to allow interactions of the arginine side chain with lipid headgroup phosphates.

Potassium ion channels form pores in cell membranes, allowing potassium ions through while preventing the passage of sodium ions. Despite numerous high-resolution structures, it is not yet possible to relate their structure to their single molecule function other than at a qualitative level. Over the past decade, there has been a concerted effort using molecular dynamics to capture the thermodynamics and kinetics of conduction by calculating potentials of mean force (PMF). These can be used, in conjunction with the electro-diffusion theory, to predict the conductance of a specific ion channel. Here, we calculate seven independent PMFs, thereby studying the differences between two potassium ion channels, the effect of the CHARMM CMAP forcefield correction, and the sensitivity and reproducibility of the method. Thermodynamically stable ion-water configurations of the selectivity filter can be identified from all the free energy landscapes, but the heights of the kinetic barriers for potassium ions to move through the selectivity filter are, in nearly all cases, too high to predict conductances in line with experiment. This implies it is not currently feasible to predict the conductance of potassium ion channels, but other simpler channels may be more tractable.

Although extensively studied, it has proved difficult to describe in detail how potassium ion channels conduct cations and water. We present a computational study that, by using stratified umbrella sampling, examines nearly an entire conduction event of the Kv1.2/2.1 paddle chimera and thereby identifies the expected stable configurations of ions and waters in the selectivity filter of the channel. We describe in detail the motions of the ions and waters during a conduction event, focusing on how waters and ions enter the filter, the rotation of water molecules inside the filter, and how potassium ions are coordinated as they move from a water to a protein environment. Finally, we analyze the small conformational changes undergone by the protein, showing that the stable configurations are most similar to the experimental crystal structure.

Protein-lipid interactions regulate many membrane protein functions. Using a multiscale approach that combines coarse-grained and atomistic molecular dynamics simulations, we have predicted the binding site for the anionic phospholipid phosphatidylinositol 4,5-bisphosphate (PIP(2)) on the Kir2.2 inwardly rectifying (Kir) potassium channel. Comparison of the predicted binding site to that observed in the recent PIP(2)-bound crystal structure of Kir2.2 reveals good agreement between simulation and experiment. In addition to providing insight into the mechanism by which PIP(2) binds to Kir2.2, these results help to establish the validity of this multiscale simulation approach and its future application in the examination of novel membrane protein-lipid interactions in the increasing number of high-resolution membrane protein structures that are now available.

Membrane proteins interact with their lipid bilayer environment via both a transmembrane helix and juxtamembrane regions. The effect of juxtamembrane regions and membrane lipid composition on these interactions has been explored by multiscale molecular dynamics simulations. The consequences of anionic lipids within the inner leaflet of a membrane were studied in combination with membrane spanning protein models differing in their juxtamembrane domains. The simulations reveal sensitivity of the protein-lipid interactions to membrane lipid composition and charged amino acid side chains. Basic residues on the intracellular side of the protein facilitated interactions with anionic lipids. Protein systems without basic residues do not show selectivity for anionic compared with zwitterionic lipids. This reveals the sensitivity to the composition of both the membrane and the protein system when studying membrane-embedded proteins. The results presented here illustrate how even a simple transmembrane domain is able to induce lipid reorganization in a mixed asymmetric bilayer.

Fusion peptides of type I fusion glycoproteins are structural elements of several enveloped viruses which enable the fusion between host and virus membranes. It is generally suggested that these peptides can promote the early fusion steps by inducing membrane curvature and that they adopt a tilted helical conformation in membranes. Although this property has been the subject of several experimental and in silico studies, an extensive sampling of the membrane peptide interaction has not yet been done. In this study, we performed coarse-grained molecular dynamic simulations in which the lipid bilayer self-assembles around the peptide. The simulations indicate that the SIV fusion peptide can adopt two different orientations in a DPPC bilayer, a major population which adopts a tilted interfacial orientation and a minor population which is perpendicular to the bilayer. The simulations also indicate that for the SIV mutant that does not induce fusion in vitro the tilt is abolished.

Both Trp and Arg in transmembrane protein domains make important interactions with lipids at the membrane/water interface, but at different depths. Derivatives of the designed peptide GWALP23, acetyl-GGALW(5)LALALALALALALW(19)LAGA-amide, with single Trp anchors, have proven useful for characterizing such interactions. Indeed, previous work revealed quite different effects emanating from Arg substitutions at positions 12 and 14 within GWALP23, with the R12 peptide exhibiting multiple positions and orientations with respect to DOPC bilayer membranes (Vostrikov et al. J. Am. Chem. Soc. 2010, 132, 5803-5811). To gain further understanding of the multistate behavior, we moved the Trp "anchor" residues to more outer positions 3 and 21 in GWALP23 itself, and in the R12 and R14 derivatives. The locations and orientations of the peptides with respect to lipid bilayer membranes of differing thickness were investigated by means of solid-state (2)H NMR spectroscopy, using labeled alanines, and coarse-grained molecular dynamics simulations. Interestingly, relatively intense and narrow (2)H resonances from selected backbone C(α) deuterons were observed over quite narrow ranges of frequency and sample orientation. The backbone resonances reflect dynamic complexities and at the same time provide important contributions for the analysis of peptide transmembrane orientation. With the Trp(3,21) anchors relatively far from the peptide and bilayer center, the results indicate significantly large apparent tilt angles, for example, close to 30° for the new R12 and R14 peptides with respect to the bilayer normal of DLPC membranes. The R12 side chain indeed is "rescued" to a stable position, where it is accommodated within the transmembrane helix, when the Trp anchors are moved outward and to another face of the helix. At the same time, the R14 side chain of transmembrane GW(3,21)ALP23 also retains a stable favored position.

The antimicrobial peptide maculatin 1.1 (M1.1) is an amphipathic α-helix that permeabilizes lipid bilayers. In coarse-grained molecular dynamics (CG MD) simulations, M1.1 has previously been shown to form membrane-spanning aggregates in DPPC bilayers. In this study, a simple multiscale methodology has been applied to allow sampling of important regions of the free energy surface at higher resolution. Thus, by back-converting the CG configurations to atomistic representations, it is shown that water is able to permeate through the M1.1 aggregates. Investigation of aggregate stoichiometry shows that at least six peptides are required for water permeation. The aggregates are dynamically disordered structures, and water flux occurs through irregular, fluctuating channels. The results are discussed in relation to experimental data and other simulations of antimicrobial peptides.

We developed molecular models for the cystic fibrosis transmembrane conductance regulator chloride channel based on the prokaryotic ABC transporter, Sav1866. Here we analyze predicted pore geometry and side-chain orientations for TM3, TM6, TM9, and TM12, with particular attention being paid to the location of the rate-limiting barrier for anion conduction. Side-chain orientations assayed by cysteine scanning were found to be from 77 to 90% in accord with model predictions. The predicted geometry of the anion conduction path was defined by a space-filling model of the pore and confirmed by visualizing the distribution of water molecules from a molecular dynamics simulation. The pore shape is that of an asymmetric hourglass, comprising a shallow outward-facing vestibule that tapers rapidly toward a narrow "bottleneck" linking the outer vestibule to a large inner cavity extending toward the cytoplasmic extent of the lipid bilayer. The junction between the outer vestibule and the bottleneck features an outward-facing rim marked by T338 in TM6 and I1131 in TM12, consistent with the observation that cysteines at both of these locations reacted with both channel-permeant and channel-impermeant, thiol-directed reagents. Conversely, cysteines substituted for S341 in TM6 or T1134 in TM12, predicted by the model to lie below the rim of the bottleneck, were found to react exclusively with channel-permeant reagents applied from the extracellular side. The predicted dimensions of the bottleneck are consistent with the demonstrated permeation of Cl(-), pseudohalide anions, water, and urea.

Combining a high degree of selectivity and nanoscale dimensions, biological pores are attractive potential components for nanotechnology devices and applications. Biomimetic design will facilitate production of stable synthetic nanopores with defined functionality. Bacterial porins offer a good source of possible templates for such nanopores as they form stable, selective pores in lipid bilayers. Molecular dynamics simulations have been used to design simple model nanopores with permeation free energy profiles that can be made to mimic a template protein, the OprP porin, which forms pores selective for anions. In particular, we explored the effects of varying the nature of pore-lining groups on free energy profiles for phosphate and chloride ions along the pore axis and the total charge of the permeation pathway of the pore. Cationic side chains lining the model nanopore are required to model the local detail of the OprP permeation landscape, whereas the total charge contributes to its magnitude. These studies indicate that a locally accurate biomimetic design has captured the essentials of the structure/function relationship of the parent protein.

Coarse-grained molecular dynamics provides a means for simulating the assembly and the interactions of membrane protein/lipid complexes at a reduced level of representation, allowing longer and larger simulations. We describe a fragment-based protocol for converting membrane simulation systems, comprising a membrane protein embedded in a phospholipid bilayer, from coarse-grained to atomistic resolution, for further refinement and analysis via atomistic simulations. Overall, this provides a method for generating an accurate and well equilibrated membrane protein/lipid complex. We exemplify the protocol using the acid-sensing/amiloride-sensitive ion channel protein (ASIC) channel protein, a trimeric integral membrane protein. The method is further evaluated using a test set of 10 different membrane proteins of differing size and complexity. Simulations are assessed in terms of protein conformational drift, lipid/protein interactions, and lipid dynamics.

Independent experimental and computational approaches show agreement concerning arginine/membrane interactions when a single arginine is introduced at selected positions within the membrane-spanning region of acetyl-GGALW(5)LALALAL(12)AL(14)ALALW(19)LAGA-ethanolamide, designated GWALP23. Peptide sequence isomers having Arg in position 12 or position 14 display markedly different behaviors, as deduced by both solid-state NMR experiments and coarse-grained molecular dynamics (CG-MD) simulations. With respect to the membrane normal of DOPC or DPPC lipid bilayer membranes, GWALP23-R14 shows one major state whose apparent average tilt is approximately 10 degrees greater than that of GWALP23. The presence of R14 furthermore induces bilayer thinning and peptide displacement to "lift" the charged guanidinium toward the bilayer surface. By contrast, GWALP23-R12 exhibits multiple states that are in slow exchange on the NMR time scale, with CG-MD simulations indicating two distinct positions with different screw rotation angles in the membrane, along with an increased tendency to exit the lipid bilayer.

Membrane proteins account for nearly a quarter of all genes, but their structure and function remain incompletely understood. Most membrane proteins have transmembrane (TM) domains made up of bundles of hydrophobic alpha-helices. The lateral association of TM helices within the lipid bilayer is a key stage in the folding of membrane proteins. It may also play a role in signaling across cell membranes. Dimerization of TM helices is a simple example of such lateral association. Molecular dynamics (MD) simulations have been used for over a decade to study membrane proteins in a lipid bilayer environment. However, direct atomistic (AT) MD simulation of self-assembly of a TM helix bundle remains challenging. AT-MD may be complemented by coarse-grained (CG) simulations, in which small numbers of atoms are grouped together into particles. In this Account, we demonstrate how CG-MD may be used to simulate formation of dimers of TM helices. We also show how a serial combination of CG and AT simulation provides a multiscale approach for generating and refining models of TM helix dimers. The glycophorin A (GpA) TM helix dimer represents a paradigm for helix-helix packing, mediated by a GxxxG sequence motif. It is well characterized experimentally and so is a good test case for evaluating computational methods. CG-MD simulations in which two separate TM helices are inserted in a lipid bilayer result in spontaneous formation of a right-handed GpA dimer, in agreement with NMR structures. CG-MD models were evaluated via comparison with data on destabilizing mutants of GpA. Such mutants increased the conformational flexibility and the dissociation constants of helix dimers. GpA dimers have been used to evaluate a multiscale approach: A CG model is converted to an AT model, which is used as the basis of an AT-MD simulation. Comparison of three AT-MD simulations of GpA, one starting from a CG model and two starting from NMR structures, leads to convergence to a common refined structure for the dimer. CG-MD self-assembly has also been used to model dimerization of the TM domain of the syndecan-2 receptor protein. This TM helix contains a GxxxG motif, which mediates right-handed helix packing comparable to that of the GxxxG motif in GpA. The multiscale approach has been applied to a more complex system, the heterodimeric alphaIIb/beta3 integrin TM helix dimer. In CG-MD, both right-handed and left-handed structures were formed. Subsequent AT-MD simulations showed that the right-handed structure was more stable, yielding a dimer in which the GxxxG motif of the alphaIIb TM helix packed against a hydrophobic surface of the beta3 helix in a manner comparable to that observed in two recent NMR studies. This work demonstrates that the multiscale simulation approach can be used to model simple membrane proteins. The method may be applied to more complex proteins, such as the influenza M2 channel protein. Future refinements, such as extending the multiscale approach to a wider range of scales (from CG through QM/MM simulations, for example), will expand the range of applications and the accuracy of the resultant models.

Computing free energies of complex biomolecular systems via atomistic (AT) molecular dynamics (MD) simulations remains a challenge due to the need for adequate sampling and convergence. Recent coarse-grained (CG) methodology allows simulations of significantly larger systems (∼10(6) to 10(8) atoms) over longer (μs/ms) time scales. Such CG models appear to be capable of making semiquantitative predictions. However, their ability to reproduce accurate thermodynamic quantities remains uncertain. We have recently used CG MD simulations to compute the potential of mean force (PMF) or free energy profile of a small peptide toxin interacting with a lipid bilayer along a 1D reaction coordinate. The toxin studied was VSTx1 (Voltage Sensor Toxin 1) from spider venom which inhibits the archeabacterial voltage-gated potassium (Kv) channel KvAP by binding to the voltage-sensor (VS) domains. Here, we re-estimate this PMF profile using (i) AT MD simulations with explicit membrane and solvent and (ii) an implicit membrane and solvent (generalized Born; GBIM) model where only the peptide was explicit. We used the CG MD free energy simulations to guide the setup of the corresponding AT MD simulations. The aim was to avoid local minima in the AT simulations which would be difficult over shorter AT time scales. A cross-comparison of the PMF profiles revealed a conserved topology, although there were differences in the magnitude of the free energies. The CG and AT simulations predicted a membrane/water interface free energy well of -27 and -23 kcal/mol, respectively (with respect to water). The GBIM model, however, gave a reduced interfacial free energy well (-12 kcal/mol). In addition, the CG and GBIM models predicted a free energy barrier of +61 and +96 kcal/mol, respectively, for positioning the toxin at the center of the bilayer, which was considerably smaller in the AT simulations (+26 kcal/mol). Thus, we present a framework for serially combining CG and AT simulations to estimate the free energy of peptide/membrane interactions. Such approaches for combining simulations at different levels of granularity will become increasingly important in future studies of complex membrane/protein systems.

Phosphatidylinositol bisphosphate (PIP(2)) is an activator of mammalian inwardly rectifying potassium (Kir) channels. Multiscale simulations, via a sequential combination of coarse-grained and atomistic molecular dynamics, enabled exploration of the interactions of PIP(2) molecules within the inner leaflet of a lipid bilayer membrane with possible binding sites on Kir channels. Three Kir channel structures were investigated: X-ray structures of KirBac1.1 and of a Kir3.1-KirBac1.3 chimera and a homology model of Kir6.2. Coarse-grained simulations of the Kir channels in PIP(2)-containing lipid bilayers identified the PIP(2)-binding site on each channel. These models of the PIP(2)-channel complexes were refined by conversion to an atomistic representation followed by molecular dynamics simulation in a lipid bilayer. All three channels were revealed to contain a conserved binding site at the N-terminal end of the slide (M0) helix, at the interface between adjacent subunits of the channel. This binding site agrees with mutagenesis data and is in the proximity of the site occupied by a detergent molecule in the Kir chimera channel crystal. Polar contacts in the coarse-grained simulations corresponded to long-lived electrostatic and H-bonding interactions between the channel and PIP(2) in the atomistic simulations, enabling identification of key side chains.

BM2 is a small integral membrane protein from influenza B virus which forms proton-permeable channels. Coarse-grained (CG) molecular dynamics simulations have been used to produce a model of the BM2 channel by self-assembly of a tetrameric bundle of BM2 transmembrane helices in a lipid bilayer. The BM2 channel model is conformationally stable on a 5 mus time scale. This CG model was converted to atomistic resolution to refine interhelix and channel-water interactions. Atomistic molecular dynamics simulations indicate that the BM2 channel is closed when no more than two of the four His19 residues are protonated. Protonating a third His19 side chain initiates a conformational change that opens the channel. In summary, our simulations suggest a common mechanism for BM2 and A/M2, whereby changes in helix packing play a functional role in channel gating.

The Kir2.1 potassium channel owes its inward-rectifying behavior to blocking by multivalent ions, e.g., magnesium and spermine, which access the channel from the cytoplasm and are thought to bind within the pore. To investigate the pathway followed by these ions from the cytoplasm through the pore, we have used multiscale modeling (via continuum electrostatics calculations, docking, and molecular dynamics simulations) to identify possible binding sites en route. On its way to eventually binding in the cavity, magnesium interacts extensively with Glu299, which lines the pore in the center of the intracellular domain. Interaction sites for spermine are formed by Asp255, Glu299, and Glu224. Entropic factors seem to favor interactions of spermine within the center of the cytoplasmic domain.

Coarse-grained molecular dynamics (CG-MD) may be applied as part of a multiscale modeling approach to protein-protein interactions. The cohesin-dockerin interaction provides a valuable test system for evaluation of the use of CG-MD, as structural (X-ray) data indicate a dual binding mode for the cohesin-dockerin pair. CG-MD simulations (of 5 μs duration) of the association of cohesin and dockerin identify two distinct binding modes, which resemble those observed in X-ray structures. For each binding mode, ca. 80% of interfacial residues are predicted correctly. Furthermore, each of the binding modes identified by CG-MD is conformationally stable when converted to an atomistic model and used as the basis of a conventional atomistic MD simulation of duration 20 ns.

9(10),16-Dihydroxypalmitic acid (diHPA) is a particularly interesting polyhydroxylated fatty acid (1) because it is the main monomer of cutin, the most abundant biopolyester in nature, and (2) because the presence of a terminal and a secondary hydroxyl group in midchain positions provides an excellent model to study their intermolecular interactions in a confined phase such as self-assembled layers. In this study we have combined atomic force microscopy (AFM), X-ray photoelectron spectroscopy (XPS), attenuated total reflection Fourier transform infrared (ATR-FT-IR) spectroscopy, as well as molecular dynamics (MD) simulations to conclude that the self-assembling of diHPA molecules on mica is a layer by layer process following a Brunauer-Emmett-Teller (BET) type isotherm and with the first layer growing much faster than the rest. Interactions between secondary hydroxyls reinforce the cohesive energy of the monolayer, while the presence of the terminal hydroxyl group is necessary to trigger the multilayered growth. Besides, XPS and ATR-FT-IR spectroscopies clearly indicate that spontaneous self-esterification occurs upon self-assembling. The esterification reaction is a prerequisite to propose a self-assembly route for the biosynthesis of cutin in nature. Molecular dynamics simulations have shown that internal molecular reorganization within the self-assembled layers provides the appropriate intermolecular orientation to facilitate the nucleophilic attack and the release of a water molecule required by the esterification reaction.

Molecular dynamics simulations have been used to explore the interactions of two PI(3)P-binding domains with their PI ligands and with a phospholipid bilayer. Three simulations each of the EEA1-FYVE domain and the p40(phox)-PX domain have been compared: with the protein in an apo state, with a bound Ins(1,3)P(2) molecule, and bound to a PI(3)P molecule embedded in a lipid bilayer. Two main questions were addressed in analysis of the simulations: (i) the location of these domains relative to the lipid bilayer and (ii) their interactions with the lipids, both specific interactions via bound PI(3)P and nonspecific interactions with bilayer phospholipids. Both domains underwent a decrease in dynamic flexibility on binding to the ligand and to the membrane, this being more pronounced for the FYVE domain. Compared to their starting locations [docked to a membrane-inserted PI(3)P molecule], each of the domains penetrated more deeply into the lipid bilayer. For FYVE, nonspecific protein-lipid interactions were formed mainly by the N-terminal hydrophobic region of the protein. For PX, both the alpha1-alpha2 and the beta1-beta2 regions penetrated the bilayer. There appeared to be more marked dynamic fluctuations in hydrogen bonds between basic side chains and PI(3)P for FYVE than for PX, but for both domains, such interactions were maintained throughout the simulations. The simulations agree well with available biophysical data, suggesting this computational method may be used to predict protein-bilayer interactions for other PI-binding proteins.

Monotopic membrane proteins bind tightly to cell membranes but do not generally span the lipid bilayer. Their interactions with lipid bilayers may be studied via coarse-grained molecular dynamics (CG-MD) simulations. Understanding such interactions is important as monotopic enzymes frequently act on hydrophobic substrates, while X-ray structures rarely provide direct information about their interactions with membranes. CG-MD self-assembly simulations enable prediction of the orientation and depth of insertion into a lipid bilayer of a monotopic protein, and also of the interactions of individual protein residues with lipid molecules. The CG-MD method has been evaluated via comparison with extended (>30 ns) atomistic simulations of monoamine oxidase, revealing good agreement between the results of coarse-grained and atomistic simulations. CG-MD simulations have been applied to a set of 11 monotopic proteins for which three-dimensional structures are available. These proteins may be divided into two groups on the basis of the results of the simulations. One group consists of those proteins which are inserted into the lipid bilayer to a limited extent, interacting mainly at the phospholipid-water interface. The second group consists of those which are inserted more deeply into the bilayer. Those monotopic proteins which are inserted more deeply cause significant local perturbation of bilayer properties such as bilayer thickness. Deeper insertion seems to correlate with a greater number of basic residues in the "foot" whereby a monotopic protein interacts with the membrane.

Experimental and computational studies have indicated that hydrophobicity plays a key role in driving the insertion of transmembrane alpha-helices into lipid bilayers. Molecular dynamics simulations allow exploration of the nature of the interactions of transmembrane alpha-helices with their lipid bilayer environment. In particular, coarse-grained simulations have considerable potential for studying many aspects of membrane proteins, ranging from their self-assembly to the relation between their structure and function. However, there is a need to evaluate the accuracy of coarse-grained estimates of the energetics of transmembrane helix insertion. Here, three levels of complexity of model system have been explored to enable such an evaluation. First, calculated free energies of partitioning of amino acid side chains between water and alkane yielded an excellent correlation with experiment. Second, free energy profiles for transfer of amino acid side chains along the normal to a phosphatidylcholine bilayer were in good agreement with experimental and atomistic simulation studies. Third, estimation of the free energy profile for transfer of an arginine residue, embedded within a hydrophobic alpha-helix, to the center of a lipid bilayer gave a barrier of approximately 15 kT. Hence, there is a substantial barrier to membrane insertion for charged amino acids, but the coarse-grained model still underestimates the corresponding free energy estimate (approximately 29 kT) from atomistic simulations (Dorairaj, S., and Allen, T. W. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 4943-4948). Coarse-grained simulations were then used to predict the free energy profile for transfer of a simple model transmembrane alpha-helix (WALP23) across a lipid bilayer. The results indicated that a transmembrane orientation was favored by about -70 kT.

Oligomerization of transmembrane (TM) helices is a key stage in the folding of membrane proteins. Glycophorin A (GpA) is a well-documented test system for this process. Coarse-grained molecular dynamics (CG-MD) allows us to simulate the self-assembly of TM helices into dimers, for both wild-type (WT) and mutant GpA sequences. For the WT sequences, dimers formed rapidly and remained stable in all simulations. The resultant dimers exhibited right-handed crossing and the same interhelix contacts as in NMR structures. Simulations of disruptive mutants revealed the dimers were less stable, with values of DeltaDelta G dimerization consistent with experimental data. The dimers of disruptive mutants were distorted relative to the WT and showed left-handed crossing of their helices. CG-MD can therefore be used to explore the interactions of TM helices, an important stage in the folding of membrane proteins. In particular, CG-MD has been shown to be sensitive enough to detect disruptions introduced by mutation. Future refinement of such models via atomistic simulations will enable a multiscale approach to predict the folding of membrane proteins.

Carbon nanotubes (CNTs) are possible nanoinjectors for the introduction of therapeutic agents into cells. To explore their interactions with a lipid bilayer membrane and to model the nanoinjection process, we used coarse-grained molecular dynamics to simulate the penetration of dipalmitoylphosphatidylcholine (DPPC) bilayers by single-walled CNTs. Lipids are extracted from a bilayer during CNT penetration and reside on both the inner and the outer tube surfaces. Lipids that interact with the CNT interior wall spread out and hence can "block" the tube. However, the degree of lipid lining of the inner surface is strongly dependent upon the tube penetration velocity, with fewer lipids extracted from the bilayer at higher rates. There is no apparent effect on bilayer integrity after CNT penetration, with the bilayer able to self-seal. Our findings reveal some of the complexities of the interactions of lipids with CNT nanoinjectors and suggest a need to further characterize the influence of, for example, CNT functionalization and cargo on lipid blocking of CNTs.

Potassium (K (+)) channels can regulate ionic conduction through their pore by a mechanism, involving the selectivity filter, known as C-type inactivation. This process is rapid in the hERG K (+) channel and is fundamental to its physiological role. Although mutations within hERG are known to remove this process, a structural basis for the inactivation mechanism has yet to be characterized. Using MD simulations based on homology modeling, we observe that the carbonyl of the filter aromatic, Phe627, forming the S 0 K (+) binding site, swiftly rotates away from the conduction axis in the wild-type channel. In contrast, in well-characterized non-inactivating mutant channels, this conformational change occurs less frequently. In the non-inactivating channels, interactions with a water molecule located behind the selectivity filter are critical to the enhanced stability of the conducting state. We observe comparable conformational changes in the acid sensitive TASK-1 channel and propose a common mechanism in these channels for regulating efflux of K (+) ions through the selectivity filter.

The computational costs associated with performing molecular dynamics (MD) simulations are still somewhat prohibitive and therefore limit the time and length scales that can be currently achieved. One approach to overcoming the limited size and duration of a simulation is to reduce the amount of detail when representing a system of interest, generally termed "coarse-graining". An alternative approach is via more efficient sampling methods that offer an enhanced search of a complex multidimensional energy landscape. One could also combine enhanced sampling methods with a coarse-grained (CG) force field. Here, we apply generalized shadow hybrid Monte Carlo (GSHMC), a recently proposed simulation protocol, to a biomolecular system of moderate size and show that GSHMC offers improved sampling compared to standard MD simulation. Our test system is a CG representation of a small peptide toxin interacting with a phospholipid bilayer. Specifically, we show that GSHMC allows for a quicker localization of the toxin to its equilibrium location of interaction at the headgroup/water interface of the bilayer. GSHMC therefore potentially allows for future exploration of larger and more complex systems over longer periods, which would otherwise be impractical to perform using conventional simulation methodology.

The mechanism of interaction of pleckstrin homology (PH) domains with phosphatidylinositol 4,5-bisphosphate (PIP 2)-containing lipid bilayers remains uncertain. While crystallographic studies have emphasized PH-inositol 1,4,5-trisphosphate (IP 3) interactions, biophysical studies indicate a degree of less specific protein-bilayer interactions. We have used molecular dynamics simulations to characterize the interactions of the PH domain from phospholipase C-delta1 with IP 3 and with PIP 2, the latter in lipid bilayers and in detergent micelles. Simulations of the PH domain in water reveal a reduction in protein flexibility when IP 3 is bound. Simulations of the PH domain bound to PIP 2 in lipid bilayers indicate a tightening of ligand-protein interactions relative to the PH-IP 3 complex, alongside formation of H-bonds between PH side chains and lipid (PC) headgroups, and a degree of penetration of hydrophobic side chains into the core of the bilayer. Comparison with simulations of the PH-bound domain to a PC bilayer in the absence of PIP 2 suggests that the presence of PIP 2 increases the extent of PH-membrane interactions. Thus, comparative molecular dynamics simulations reveal how a PI-binding domain undergoes changes in conformational dynamics on binding to a PIP 2-containing membrane and how interactions additional to those with the PI headgroup are formed.

Saposin B is a water soluble alpha-helical protein which can bind to membranes and extract selected lipids, especially cerebroside sulfates. The X-ray structure of saposin B is homodimeric. There are two conformations of the dimer in the crystal-one with a closed central cavity (the AB dimer) and one (the CD dimer) with a more open cavity. We have conducted a series of short (5 ns) molecular dynamics simulations of saposin B, starting from both the AB and CD conformations and with/without bound lipid and/or water molecules within the central hydrophobic cavity. The more open (CD) dimer showed greater conformational drift than the AB dimer. The conformational drift was also somewhat higher in the absence of bound lipid. Two more extended (30 ns) simulations of AB and CD dimers were performed and analyzed in terms of changes in intersubunit packing within the dimers. The AB dimer remained largely unchanged in conformation over the duration of the extended simulation. In contrast, the CD dimer underwent a substantial conformational change corresponding to a 'scissor' motion of the two monomers so as to compress the central cavity to a more closed conformation than that seen in the AB dimer structure. A H-bond between the Q53 and Y54 side chains of the alpha3 helices of the two opposing monomers seems to hold the dimer in this 'scissor-closed' conformation. We suggest that a cycle of conformational changes, expanding and compressing the central cavity of the saposin B dimer, may play a key role in facilitating lipid extraction from bilayers.

Detergent interactions with carbon nanotubes are of potential importance in a number of bionanotechnology applications. We investigate the interaction of lysophospholipids with single-walled carbon nanotubes via coarse-grained molecular dynamics. We present compelling evidence that the mechanism of adsorption of these detergents onto a carbon nanotube is dependent upon detergent concentration. Furthermore, the chirality of the carbon nanotube influences the detergent wrapping angle for low detergent concentration. These findings advance our understanding of the mechanism of carbon nanotube solubilization via detergent molecules.

Inward rectifier potassium (Kir) channels regulate cell excitability and transport K+ ions across membranes. Homotetrameric models of three mammalian Kir channels (Kir1.1, Kir3.1, and Kir6.2) have been generated, using the KirBac3.1 transmembrane and rat Kir3.1 intracellular domain structures as templates. All three models have been explored by 10 ns molecular dynamics simulations in phospholipid bilayers. Analysis of the initial structures revealed conservation of potential lipid interaction residues (Trp/Tyr and Arg/Lys side chains near the lipid headgroup-water interfaces). Examination of the intracellular domains revealed key structural differences between Kir1.1 and Kir6.2 which may explain the difference in channel inhibition by ATP. The behavior of all three models in the MD simulations revealed that they have conformational stability similar to that seen for comparable simulations of, for example, structures derived from cryoelectron microscopy data. Local distortions of the selectivity filter were seen during the simulations, as observed in previous simulations of KirBac and in simulations and structures of KcsA. These may be related to filter gating of the channel. The intracellular hydrophobic gate does not undergo any substantial changes during the simulations and thus remains functionally closed. Analysis of lipid-protein interactions of the Kir models emphasizes the key role of the M0 (or "slide") helix which lies approximately parallel to the bilayer-water interface and forms a link between the transmembrane and intracellular domains of the channel.

Monotopic proteins make up a class of membrane proteins that bind tightly to, but do not span, cell membranes. We examine and compare how two monotopic proteins, monoamine oxidase B (MAO-B) and cyclooxygenase-2 (COX-2), interact with a phospholipid bilayer using molecular dynamics simulations. Both enzymes form between three and seven hydrogen bonds with the bilayer in our simulations with basic side chains accounting for the majority of these interactions. By analyzing lipid order parameters, we show that, to a first approximation, COX-2 disrupts only the upper leaflet of the bilayer. In contrast, the top and bottom halves of the lipid tails surrounding MAO-B are more and less ordered, respectively, than in the absence of the protein. Finally, we identify which residues are important in binding individual phospholipids by counting the number and type of lipid atoms that come close to each amino acid residue. The existing models that explain how these proteins bind to bilayers were proposed following inspection of the X-ray crystallographic structures. Our results support these models and suggest that basic residues contribute significantly to the binding of these monotopic proteins to bilayers through the formation of hydrogen bonds with phospholipids.