Computer simulations of protein-membrane systems.
Ontology highlight
ABSTRACT:
SUBMITTER: Loschwitz J
PROVIDER: S-EPMC7109768 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The conformational sampling of monomeric, membrane-bound phospholamban is described from computer simulations. Phospholamban (PLB) plays a key role as a regulator of sarcoplasmic reticulum calcium ATPase. An implicit membrane model is used in conjunction with replica exchange molecular dynamics simulations to reach mus-ms timescales. The implicit membrane model was also used to study the effect of different membrane thicknesses by scaling the low-dielectric region. The conformational sampling with the membrane model mimicking dipalmitoylphosphatidylcholine bilayers is in good agreement overall with experimental measurements, but consists of a wide variety of different conformations including structures not described previously. The conformational ensemble shifts significantly in the presence of thinner or thicker membranes. This has implications for the structure and dynamics of PLB in physiological membranes and offers what we believe to be a new interpretation of previous experimental measurements of PLB in detergents and microsomal membrane.
| S-EPMC2830431 | biostudies-literature
Project description:Transient protein-protein interactions occur frequently under the crowded conditions encountered in biological environments. Tryptophan-cysteine quenching is introduced as an experimental approach with minimal labeling for characterizing such interactions between proteins due to its sensitivity to nano- to microsecond dynamics on subnanometer length scales. The experiments are paired with computational modeling at different resolutions including fully atomistic molecular dynamics simulations for interpretation of the experimental observables and to gain molecular-level insights. This approach is applied to model systems, villin variants and the drkN SH3 domain, in the presence of protein G crowders. It is demonstrated that Trp-Cys quenching experiments can differentiate between overall attractive and repulsive interactions between different proteins, and they can discern variations in interaction preferences at different protein surface locations. The close integration between experiment and simulations also provides an opportunity to evaluate different molecular force fields for the simulation of concentrated protein solutions.
| S-EPMC9870652 | biostudies-literature
Project description:A theoretical framework is elaborated to account for the effect of a transmembrane potential in computer simulations. It is shown that a simulation with a constant external electric field applied in the direction normal to the membrane is equivalent to the influence of surrounding infinite baths maintained to a voltage difference via ion-exchanging electrodes connected to an electromotive force. It is also shown that the linearly-weighted displacement charge within the simulation system tracks the net flow of charge through the external circuit comprising the electromotive force and the electrodes. Using a statistical mechanical reduction of the degrees of freedom of the external system, three distinct theoretical routes are formulated and examined for the purpose of characterizing the free energy of a protein embedded in a membrane that is submitted to a voltage difference. The W-route is constructed from the variations in the voltage-dependent potential of mean force along a reaction path connecting two conformations of the protein. The Q-route is based on the average displacement charge as a function of the conformation of the protein. Finally, the G-route considers the relative charging free energy of specific residues, with and without applied membrane potentials. The theoretical formulation is illustrated with a simple model of an ion crossing a vacuum slab surrounded by two aqueous bulk phases and with a fragment of the voltage-sensor of the KvAP potassium channel.
| S-EPMC2567939 | biostudies-literature
Project description:A general methodology for calculating the equilibrium binding constant of a flexible ligand to a protein receptor is formulated on the basis of potentials of mean force. The overall process is decomposed into several stages that can be computed separately: the free ligand in the bulk is first restrained into the conformation it adopts in the bound state, position, and orientation by applying biasing potentials, then it is translated into the binding site, where it is released completely. The conformational restraining potential is based on the root-mean-square deviation of the peptide coordinates relative to its average conformation in the bound complex. Free energy contributions from each stage are calculated by means of free energy perturbation potential of mean force techniques by using appropriate order parameters. The present approach avoids the need to decouple the ligand from its surrounding (bulk solvent and receptor protein) as is traditionally performed in the double-decoupling scheme. It is believed that the present formulation will be particularly useful when the solvation free energy of the ligand is very large. As an application, the equilibrium binding constant of the phosphotyrosine peptide pYEEI to the Src homology 2 domain of human Lck has been calculated. The results are in good agreement with experimental values.
| S-EPMC1100764 | biostudies-literature
Project description: Not available
| S-EPMC5504919 | biostudies-literature
Project description:Using the endosomal sorting complex required for transport (ESCRT)-III membrane remodeling complex as an example, we analyze three popular coarse-grained models (the regular MARTINI, polarizable MARTINI (POL-MARTINI), and big multipole water MARTINI (BMW-MARTINI)) for the description of membrane curvature sensing and generation activities of peripheral proteins. Although the three variants of the MARTINI model provide consistent descriptions for the protein-protein interface in a linear filament model of ESCRT-III, they differ considerably in terms of protein-membrane interface and therefore membrane curvature sensing and generation behaviors. In particular, BMW-MARTINI provides the most consistent description of the protein-membrane interface as compared to all-atom simulations, whereas the regular MARTINI is most consistent with atomistic simulations in terms of the qualitative sign of membrane curvature sensing and generation. With POL-MARTINI, the ESCRT-III model interacts weakly with the membrane and therefore does not exhibit any curvature-sensitive activities. Analysis suggests that the incorrect membrane curvature activities predicted by BMW-MARTINI are due to overestimated insertion depth of an amphipathic helix and incorrect sign for the spontaneous curvature of anionic lipids. These results not only point to ways that coarse-grained models can be improved but also explicitly highlight local lipid composition and insertion depth of protein motifs as essential regulatory factors for membrane curvature sensing and generation.
| S-EPMC8390925 | biostudies-literature
Project description:The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.
| S-EPMC4183268 | biostudies-literature
Project description:rhodamines are dyes widely used as fluorescent tags in cell imaging, probing of mitochondrial membrane potential, and as P-glycoprotein model substrates. In all these applications, detailed understanding of the interaction between rhodamines and biomembranes is fundamental. we combined atomistic molecular dynamics (MD) simulations and fluorescence spectroscopy to characterize the interaction between rhodamines 123 and B (Rh123 and RhB, respectively) and POPC bilayers. while the xanthene moiety orients roughly parallel to the membrane plane in unrestrained MD simulations, variations on the relative position of the benzoic ring (below the xanthene for Rh123, above it for RhB) were observed, and related to the structure of the two dyes and their interactions with water and lipids. Subtle distinctions were found among different ionization forms of the probes. Experimentally, RhB displayed a lipid/water partition coefficient more than two orders of magnitude higher than Rh123, in agreement with free energy profiles obtained from umbrella sampling MD. this work provided detailed insights on the similarities and differences in the behavior of bilayer-inserted Rh123 and RhB, related to the structure of the probes. The much higher affinity of RhB for the membranes increases the local concentration and explains its higher apparent affinity for P-glycoprotein reconstituted in model membranes.
| S-EPMC8876248 | biostudies-literature
Project description:Hemagglutinin (HA) mediated fusion of influenza virus envelope with host lipid membrane is a critical step warrantying virus entry to the cell. Despite tremendous advances in structural biology methods, the knowledge concerning the details of HA2 subunit insertion into the target membrane and its subsequent bilayer perturbing effect is still rather limited. Herein, based on a set of molecular dynamics simulations, we investigate the structure and interaction with lipid membrane of the N-terminal HA2 region comprising a trimer of fusion peptides (HAfps) tethered by flexible linkers to a fragment of coiled-coil stem structure. We find that, prior to insertion into the membrane, HAfps within the trimers do not sample space individually but rather associate into a compact hydrophobic aggregate. Once within the membrane, they fold into tight helical hairpins, which remain at the lipid-water interface. However, they can also assume stable, membrane-spanning configurations of significantly increased membrane-perturbing potential. In this latter case, HAfps trimers centre around the well-hydrated transmembrane channel-forming distinct, symmetric assemblies, whose wedge-like shape may play a role in promoting membrane curvature. We also demonstrate that, following HAfps insertion, the coiled-coil stem spontaneously tilts to almost membrane-parallel orientation, reflecting experimentally observed configuration adopted in the course of membrane fusion by complete HA2 units at the rim of membrane contact zones.
| S-EPMC8830744 | biostudies-literature
Project description:Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at Lo/Ld interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.
| S-EPMC5937199 | biostudies-literature