Functional dynamics and allosteric modulation of TRPA1.

The cation channel TRPA1 is a potentially important drug target, and characterization of TRPA1 functional dynamics might help guide structure-based drug design. Here, we present results from long-timescale molecular dynamics simulations of TRPA1 with an allosteric activator, allyl isothiocyanate (AITC), in which we observed spontaneous transitions from a closed, non-conducting channel conformation into an open, conducting conformation. Based on these transitions, we propose a gating mechanism in which movement of a regulatory TRP-like domain allosterically translates into pore opening in a manner reminiscent of pore opening in voltage-gated ion channels. In subsequent experiments, we found that mutations that disrupt packing of the S4-S5 linker-TRP-like domain and the S5 and S6 helices also affected channel activity. In simulations, we also observed A-967079, a known allosteric inhibitor, binding between helices S5 and S6, suggesting that A-967079 may suppress activity by stabilizing a non-conducting pore conformation-a finding consistent with our proposed gating mechanism.

Structure and Dynamics of Cinnamycin-Lipid Complexes: Mechanisms of Selectivity for Phosphatidylethanolamine Lipids.

Cinnamycin is a lantibiotic peptide, which selectively binds to and permeabilizes membranes containing phosphatidylethanolamine (PE) lipids. As PE is a major component of many bacterial cell membranes, cinnamycin has considerable potential for destroying these. In this study, molecular dynamics simulations are used to elucidate the structure of a lipid-cinnamycin complex and the origin of selective lipid binding. The simulations reveal that cinnamycin selectively binds to PE by forming an extensive hydrogen-bonding network involving all three hydrogen atoms of the primary ammonium group of the PE head group. The substitution of a single hydrogen atom with a methyl group on the ammonium nitrogen destabilizes this hydrogen-bonding network. In addition to binding the primary ammonium group, cinnamycin also interacts with the phosphate group of the lipid through a previously uncharacterized phosphate-binding site formed by the backbone Phe10-Abu11-Phe12-Val13 moieties (Abu = 1-α-aminobutyric acid). In addition, hydroxylation of Asp15 at Cß plays a role in selective binding of PE due to its tight interaction with the charged amine of the lipid head group. The simulations reveal that the position and orientation of the peptide in the membrane depend on the type of lipid to which it binds, suggesting a reason for why cinnamycin selectively permeabilizes PE-containing membranes.

Cholesterol Interaction Sites on the Transmembrane Domain of the Hedgehog Signal Transducer and Class F G Protein-Coupled Receptor Smoothened.

Transduction of Hedgehog signals across the plasma membrane is facilitated by the class F G-protein-coupled-receptor (GPCR) Smoothened (SMO). Recent studies suggest that SMO is modulated via interactions of its transmembrane (TM) domain with cholesterol. We apply molecular dynamics simulations of SMO embedded in cholesterol containing lipid bilayers, revealing a direct interaction of cholesterol with the TM domain at regions distinct from those observed in class A GPCRs. In particular the extracellular tips of helices TM2 and TM3 form a well-defined cholesterol interaction site. Potential of mean force calculations yield a free energy landscape for cholesterol binding. Alongside analysis of equilibrium cholesterol occupancy, this reveals the existence of a dynamic "greasy patch" interaction with the TM domain of SMO, which may be compared with previously identified lipid interaction sites on other membrane proteins. These predictions provide molecular-level insights into cholesterol interactions with a class F GPCR, suggesting potential druggable sites.

Interactions of a Bacterial Cu(I)-ATPase with a Complex Lipid Environment.

Phospholipids and sterols play multiple roles in cells. In addition to establishing barriers between compartments, they also provide the matrix for assembly and function of a large variety of catalytic processes. Lipid composition is a highly regulated feature of biological membranes, yet its implications for membrane proteins are difficult problems to approach. One obstacle is the inherent complexity of observing and describing these interactions and their dynamics at a molecular and atomic level. However, lipid interactions are pivotal for membrane protein function and should be acknowledged. The enzymatic activity of several different P-type ATPases, one of the major families of ion pumping primary active transporters, has previously been shown to exhibit a strong dependence on phospholipids; however, distinguishing the effects of annular and specific lipid interactions is challenging. Here we show that the hydrolytic activity of a bacterial Cu(I)-transporting P-type ATPase (LpCopA) is stimulated by the bacterial, anionic phospholipid cardiolipin and to some extent by phosphatidylglycerol. Furthermore, multiscale molecular dynamics simulations pinpoint lipid hot spots on the membrane-spanning domain of LpCopA. Thus, using two independent methods, our study shows converging evidence that the lipid membrane composition plays an important role for LpCopA.

Multiscale molecular dynamics simulations of lipid interactions with P-glycoprotein in a complex membrane.

P-glycoprotein (P-gp) can transport a wide range of very different hydrophobic organic molecules across the membrane. Its ability to extrude molecules from the cell creates delivery problems for drugs that target proteins in the central nervous system (CNS) and also causes drug-resistance in many forms of cancer. Whether a drug will be susceptible to export by P-gp is difficult to predict and currently this is usually assessed with empirical and/or animal models. Thus, there is a need to better understand how P-gp works at the molecular level in order to fulfil the 3Rs: Refinement, reduction and replacement of animals in research. As structural information increasingly becomes available, our understanding at the molecular level improves. Proteins like P-gp are however very dynamic entities and thus one of the most appropriate ways to study them is with molecular dynamics simulations, especially as this can capture the influence of the surrounding environment. Recent parameterization developments have meant that it is now possible to simulate lipid bilayers that more closely resemble in vivo membranes in terms of their composition. In this report we construct a complex lipid bilayer that mimics the composition of brain epithelial cells and examine the interactions of it with P-gp. We find that the negatively charged phosphatidylserine lipids in the inner leaflet of the membrane tend to form an annulus around P-gp. We also observed the interaction of cholesterol with three distinct areas of the P-gp. Potential of mean force (PMF) calculations suggest that a crevice between transmembrane helices 10 and 12 has particularly favourable interaction energy for cholesterol.

Protein crowding and lipid complexity influence the nanoscale dynamic organization of ion channels in cell membranes.

Cell membranes are crowded and complex environments. To investigate the effect of protein-lipid interactions on dynamic organization in mammalian cell membranes, we have performed coarse-grained molecular dynamics simulations containing >100 copies of an inwardly rectifying potassium (Kir) channel which forms specific interactions with the regulatory lipid phosphatidylinositol 4,5-bisphosphate (PIP2). The tendency of protein molecules to cluster has the effect of organizing the membrane into dynamic compartments. At the same time, the diversity of lipids present has a marked effect on the clustering behavior of ion channels. Sub-diffusion of proteins and lipids is observed. Protein crowding alters the sub-diffusive behavior of proteins and lipids such as PIP2 which interact tightly with Kir channels. Protein crowding also affects bilayer properties, such as membrane undulations and bending rigidity, in a PIP2-dependent manner. This interplay between the diffusion and the dynamic organization of Kir channels may have important implications for channel function.

Multiscale molecular dynamics simulations of lipid interactions with P-glycoprotein in a complex membrane.

P-glycoprotein (P-gp) can transport a wide range of very different hydrophobic organic molecules across the membrane. Its ability to extrude molecules from the cell creates delivery problems for drugs that target proteins in the central nervous system (CNS) and also causes drug-resistance in many forms of cancer. Whether a drug will be susceptible to export by P-gp is difficult to predict and currently this is usually assessed with empirical and/or animal models. Thus, there is a need to better understand how P-gp works at the molecular level in order to fulfil the 3Rs: Refinement, reduction and replacement of animals in research. As structural information increasingly becomes available, our understanding at the molecular level improves. Proteins like P-gp are however very dynamic entities and thus one of the most appropriate ways to study them is with molecular dynamics simulations, especially as this can capture the influence of the surrounding environment. Recent parameterization developments have meant that it is now possible to simulate lipid bilayers that more closely resemble in vivo membranes in terms of their composition. In this report we construct a complex lipid bilayer that mimics the composition of brain epithelial cells and examine the interactions of it with P-gp. We find that the negatively charged phosphatidylserine lipids in the inner leaflet of the membrane tend to form an annulus around P-gp. We also observed the interaction of cholesterol with three distinct areas of the P-gp. Potential of mean force (PMF) calculations suggest that a crevice between transmembrane helices 10 and 12 has particularly favourable interaction energy for cholesterol.

Simulations of membrane-bound diglycosylated human prion protein reveal potential protective mechanisms against misfolding.

Prion diseases are associated with the misfolding of the prion protein (PrP) from its normal cellular form (PrPC ) to its infectious scrapie form (PrPSc ). Post-translational modifications in PrP in vivo can play an important role in modulating the process of misfolding. To gain more insight into the effects of post-translational modifications in PrP structure and dynamics and to test the hypothesis that such modifications can interact with the protein, we have performed molecular dynamics simulations of diglycosylated human PrPC bound to a lipid bilayer via a glycophosphatidylinositol anchor. Multiple simulations were performed at three different pH ranges to explore pH effects on structure and dynamics. In contrast to simulations of protein-only PrPC , no large effects were observed upon lowering the pH of the system. The protein tilted toward the membrane surface in all of the simulations and the putative PrPSc oligomerization sites became inaccessible, thereby offering a possible protective mechanism against PrPSc -induced misfolding of PrPC .

Interfacial activation of M37 lipase: A multi-scale simulation study.

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.

Ligand Binding in the Extracellular Vestibule of the Neurotransmitter Transporter Homologue LeuT.

The human monoamine transporters (MATs) facilitate the reuptake of monoamine neurotransmitters from the synaptic cleft. MATs are linked to a number of neurological diseases and are the targets of both therapeutic and illicit drugs. Until recently, no high-resolution structures of the human MATs existed, and therefore, studies of this transporter family have relied on investigations of the homologues bacterial transporters such as the leucine transporter LeuT, which has been crystallized in several conformational states. A two-substrate transport mechanism has been suggested for this transporter family, which entails that high-affinity binding of a second substrate in an extracellular site is necessary for the substrate in the central binding site to be transported. Compelling evidence for this mechanism has been presented, however, a number of equally compelling accounts suggest that the transporters function through a mechanism involving only a single substrate and a single high-affinity site. To shed light on this apparent contradiction, we have performed extensive molecular dynamics simulations of LeuT in the outward-occluded conformation with either one or two substrates bound to the transporter. We have also calculated the substrate binding affinity in each of the two proposed binding sites through rigorous free energy simulations. Results show that substrate binding is unstable in the extracellular vestibule and the substrate binding affinity within the suggested extracellular site is very low (0.2 and 3.3 M for the two dominant binding modes) compared to the central substrate binding site (14 nM). This suggests that for LeuT in the outward-occluded conformation only a single high-affinity substrate binding site exists.

Membrane stiffness is modified by integral membrane proteins.

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

Lipid-Loving ANTs: Molecular Simulations of Cardiolipin Interactions and the Organization of the Adenine Nucleotide Translocase in Model Mitochondrial Membranes.

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.

Membrane Compartmentalization Reducing the Mobility of Lipids and Proteins within a Model Plasma Membrane.

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.

Structural Basis for Simvastatin Competitive Antagonism of Complement Receptor 3.

The complement system is an important part of the innate immune response to infection but may also cause severe complications during inflammation. Small molecule antagonists to complement receptor 3 (CR3) have been widely sought, but a structural basis for their mode of action is not available. We report here on the structure of the human CR3 ligand-binding I domain in complex with simvastatin. Simvastatin targets the metal ion-dependent adhesion site of the open, ligand-binding conformation of the CR3 I domain by direct contact with the chelated Mg(2+) ion. Simvastatin antagonizes I domain binding to the complement fragments iC3b and C3d but not to intercellular adhesion molecule-1. By virtue of the I domain's wide distribution in binding kinetics to ligands, it was possible to identify ligand binding kinetics as discriminator for simvastatin antagonism. In static cellular experiments, 15-25 µm simvastatin reduced adhesion by K562 cells expressing recombinant CR3 and by primary human monocytes, with an endogenous expression of this receptor. Application of force to adhering monocytes potentiated the effects of simvastatin where only a 50-100 nm concentration of the drug reduced the adhesion by 20-40% compared with untreated cells. The ability of simvastatin to target CR3 in its ligand binding-activated conformation is a novel mechanism to explain the known anti-inflammatory effects of this compound, in particular because this CR3 conformation is found in pro-inflammatory environments. Our report points to new designs of CR3 antagonists and opens new perspectives and identifies druggable receptors from characterization of the ligand binding kinetics in the presence of antagonists.

Synthesis and inhibitory evaluation of 3-linked imipramines for the exploration of the S2 site of the human serotonin transporter.

The human serotonin transporter is the primary target of several antidepressant drugs, and the importance of a primary, high affinity binding site (S1) for antidepressant binding is well documented. The existence of a lower affinity, secondary binding site (S2) has, however, been debated. Herein we report the synthesis of 3-position coupled imipramine ligands from clomipramine using a copper free Sonogashira reaction. Ligand design was inspired by results from docking and steered molecular dynamics simulations, and the ligands were utilized in a structure-activity relationship study of the positional relationship between the S1 and S2 sites. The computer simulations suggested that the S2 site does indeed exist although with lower affinity for imipramine than observed within the S1 site. Additionally, it was possible to dock the 3-linked imipramine analogs into positions which occupy the S1 and the S2 site simultaneously. The structure activity relationship study showed that the shortest ligands were the most potent, and mutations enlarging the proposed S2 site were found to affect the larger ligands positively, while the smaller ligands were mostly unaffected.

Free Energy Landscape of Lipid Interactions with Regulatory Binding Sites on the Transmembrane Domain of the EGF Receptor.

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.

Molecular simulations of glycolipids: Towards mammalian cell membrane models.

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

Monoamine transporters: insights from molecular dynamics simulations.

The human monoamine transporters (MATs) facilitate the reuptake of the neurotransmitters serotonin, dopamine, and norepinephrine from the synaptic cleft. Imbalance in monoaminergic neurotransmission is linked to various diseases including major depression, attention deficit hyperactivity disorder, schizophrenia, and Parkinson's disease. Inhibition of the MATs is thus an important strategy for treatment of such diseases. The MATs are sodium-coupled transport proteins belonging to the neurotransmitter/Na(+) symporter (NSS) family, and the publication of the first high-resolution structure of a NSS family member, the bacterial leucine transporter LeuT, in 2005, proved to be a major stepping stone for understanding this family of transporters. Structural data allows for the use of computational methods to study the MATs, which in turn has led to a number of important discoveries. The process of substrate translocation across the membrane is an intrinsically dynamic process. Molecular dynamics simulations, which can provide atomistic details of molecular motion on ns to ms timescales, are therefore well-suited for studying transport processes. In this review, we outline how molecular dynamics simulations have provided insight into the large scale motions associated with transport of the neurotransmitters, as well as the presence of external and internal gates, the coupling between ion and substrate transport, and differences in the conformational changes induced by substrates and inhibitors.

Organization and Dynamics of Receptor Proteins in a Plasma Membrane.

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.

Insights to ligand binding to the monoamine transporters-from homology modeling to LeuBAT and dDAT.

Understanding of drug binding to the human biogenic amine transporters (BATs) is essential to explain the mechanism of action of these pharmaceuticals but more importantly to be able to develop new and improved compounds to be used in the treatment of depression or drug addiction. Until recently no high resolution structure was available of the BATs and homology modeling was a necessity. Various studies have revealed experimentally validated binding modes of numerous ligands to the BATs using homology modeling. Here we examine and discuss the similarities between the binding models of substrates, antidepressants, psychostimulants, and mazindol in homology models of the human BATs and the recently published crystal structures of the Drosophila dopamine transporter and the engineered protein, LeuBAT. The comparison reveals that careful computational modeling combined with experimental data can be utilized to predict binding of molecules to proteins that agree very well with crystal structures.