The Influence of 1-Alkanols and External Pressure on the Lateral Pressure Profiles of Lipid Bilayers

AbstractThe suggestion by Robert Cantor, that drug-induced pressure changes in lipid bilayers can change the conformational equilibrium between open and closed states of membrane proteins and thereby cause anesthesia, attracted much attention lately. Here, we studied the effect of both large external pressure and of 1-alkanols of different chain lengths—some of them anesthetics, others not—on the lateral pressure profiles across dimyristoylphosphatidylcholine (DMPC) bilayers by molecular dynamics simulations. For a pure DMPC bilayer, high pressure both reduced and broadened the tension at the interface hydrophobic/hydrophilic and diminished the repulsion between the phospholipid headgroups. Whereas the effect of ethanol on the lateral pressure profile was similar to the effect of a large external pressure on a DMPC bilayer, long-chain 1-alkanols significantly amplified local maxima and minima in the lateral pressure profile. For most 1-alkanols, external pressure had moderate effects and did not reverse the changes 1-alkanols exerted on the pressure profile. Nevertheless, assuming the bent helix model as a simple geometric model for the transmembrane region of a membrane protein, protein conformational equilibria were shifted in opposite directions by addition of 1-alkanols and additional application of external pressure

Characteristics of Sucrose Transport through the Sucrose-Specific Porin ScrY Studied by Molecular Dynamics Simulations

Sucrose-specific porin (ScrY) is a transmembrane protein that allows for the uptake of sucrose under growth-limiting conditions. The crystal structure of ScrY was resolved before by X-ray crystallography, both in its uncomplexed form and with bound sucrose. However, little is known about the molecular characteristics of the transport mechanism of ScrY. To date, there has not yet been any clear demonstration for sucrose transport through the ScrY.Here, the dynamics of the ScrY trimer embedded in a phospholipid bilayer as well as the characteristics of sucrose translocation were investigated by means of atomistic molecular dynamics (MD) simulations. The potential of mean force (PMF) for sucrose translocation through the pore showed two main energy barriers within the constriction region of ScrY. Energy decomposition allowed to pinpoint three aspartic acids as key residues opposing the passage of sucrose, all located within the L3 loop. Mutation of two aspartic acids to uncharged residues resulted in an accordingly modified electrostatics and decreased PMF barrier. The chosen methodology and results will aid in the design of porins with modified transport specificities

Structural Model of the mIgM B-Cell Receptor Transmembrane Domain From Self-Association Molecular Dynamics Simulations

Antigen binding to B-cell antigen receptors (BCRs) followed by signaling initiates the humoral immune response. The signaling is intimately coupled to nanoclustering of BCRs and their sorting to specific membrane domains, a process that is ruled by interactions between the BCR transmembrane domain and lipids. While the structure of the extracellular domains of BCRs has been resolved, little is known about the configuration of the constituting four immunoglobulin domains spanning the membrane. Here, we modeled the structure of the transmembrane (TM) domain of the IgM B-cell receptor using self-assembly coarse-grained molecular dynamics simulations. The obtained quaternary structure was validated against available experimental data and atomistic simulations. The IgM-BCR-TM domain configuration shows a 1:1 stoichiometry between the homodimeric membrane-bound domain of IgM (mIgM) and a Ig-α/Ig-β heterodimer. The mIgM homodimer is based on an asymmetric association of two mIgM domains. We show that a specific site of the Ig-α/Ig-β heterodimer is responsible for the association of IgM-BCRs with lipid rafts. Our results further suggest that this site is blocked in small-sized IgM-BCR clusters. The BCR TM structure provides a molecular basis for the previously suggested dissociation activation model of B-cell receptors. Self-assembly molecular dynamics simulations at the coarse-grained scale here proved as a versatile tool in the study of receptor complexes

The Molecular Switching Mechanism at the Conserved D(E)RY Motif in Class-A GPCRs

The disruption of ionic and H-bond interactions between the cytosolic ends of transmembrane helices TM3 and TM6 of class-A (rhodopsin-like) G protein-coupled receptors (GPCRs) is a hallmark for their activation by chemical or physical stimuli. In the bovine photoreceptor rhodopsin, this is accompanied by proton uptake at Glu134 in the class-conserved D(E)RY motif. Studies on TM3 model peptides proposed a crucial role of the lipid bilayer in linking protonation to stabilization of an active state-like conformation. However, the molecular details of this linkage could not be resolved and have been addressed in this study by molecular dynamics (MD) simulations on TM3 model peptides in a bilayer of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). We show that protonation of the conserved glutamic acid alters the peptide insertion depth in the membrane, its side-chain rotamer preferences, and stabilizes the C-terminal helical structure. These factors contribute to the rise of the side-chain pKa (> 6) and to reduced polarity around the TM3 C terminus as confirmed by fluorescence spectroscopy. Helix stabilization requires the protonated carboxyl group; unexpectedly, this stabilization could not be evoked with an amide in MD simulations. Additionally, time-resolved Fourier transform infrared (FTIR) spectroscopy of TM3 model peptides revealed a different kinetics for lipid ester carbonyl hydration, suggesting that the carboxyl is linked to more extended H-bond clusters than an amide. Remarkably, this was seen as well in DOPC-reconstituted Glu134- and Gln134-containing bovine opsin mutants and demonstrates that the D(E)RY motif is a hydrated microdomain. The function of the D(E)RY motif as a proton switch is suggested to be based on the reorganization of the H-bond network at the membrane interface

How Thioredoxin Dissociates Its Mixed Disulfide

The dissociation mechanism of the thioredoxin (Trx) mixed disulfide complexes is unknown and has been debated for more than twenty years. Specifically, opposing arguments for the activation of the nucleophilic cysteine as a thiolate during the dissociation of the complex have been put forward. As a key model, the complex between Trx and its endogenous substrate, arsenate reductase (ArsC), was used. In this structure, a Cys29Trx-Cys89ArsC intermediate disulfide is formed by the nucleophilic attack of Cys29Trx on the exposed Cys82ArsC-Cys89ArsC in oxidized ArsC. With theoretical reactivity analysis, molecular dynamics simulations, and biochemical complex formation experiments with Cys-mutants, Trx mixed disulfide dissociation was studied. We observed that the conformational changes around the intermediate disulfide bring Cys32Trx in contact with Cys29Trx. Cys32Trx is activated for its nucleophilic attack by hydrogen bonds, and Cys32Trx is found to be more reactive than Cys82ArsC. Additionally, Cys32Trx directs its nucleophilic attack on the more susceptible Cys29Trx and not on Cys89ArsC. This multidisciplinary approach provides fresh insights into a universal thiol/disulfide exchange reaction mechanism that results in reduced substrate and oxidized Trx

Interaction Pattern of Arg 62 in the A-Pocket of Differentially Disease-Associated HLA-B27 Subtypes Suggests Distinct TCR Binding Modes

The single amino acid replacement Asp116His distinguishes the two subtypes HLA-B*2705 and HLA-B*2709 which are, respectively, associated and non-associated with Ankylosing Spondylitis, an autoimmune chronic inflammatory disease. The reason for this differential association is so far poorly understood and might be related to subtype-specific HLA:peptide conformations as well as to subtype/peptide-dependent dynamical properties on the nanoscale. Here, we combine functional experiments with extensive molecular dynamics simulations to investigate the molecular dynamics and function of the conserved Arg62 of the α1-helix for both B27 subtypes in complex with the self-peptides pVIPR (RRKWRRWHL) and TIS (RRLPIFSRL), and the viral peptides pLMP2 (RRRWRRLTV) and NPflu (SRYWAIRTR). Simulations of HLA:peptide systems suggest that peptide-stabilizing interactions of the Arg62 residue observed in crystal structures are metastable for both B27 subtypes under physiological conditions, rendering this arginine solvent-exposed and, probably, a key residue for TCR interaction more than peptide-binding. This view is supported by functional experiments with conservative (R62K) and non-conservative (R62A) B*2705 and B*2709 mutants that showed an overall reduction in their capability to present peptides to CD8+ T cells. Moreover, major subtype-dependent differences in the peptide recognition suggest distinct TCR binding modes for the B*2705 versus the B*2709 subtype

Spontaneous Formation of Detergent Micelles around the Outer Membrane Protein OmpX

The structure and flexibility of the outer membrane protein X (OmpX) in a water-detergent solution and in pure water are investigated by molecular dynamics simulations on the 100-ns timescale and compared with NMR data. The simulations allow for an unbiased determination of the structure of detergent micelles and the protein-detergent mixed micelle. The short-chain lipid dihexanoylphosphatidylcholine, as a detergent, aggregates into pure micelles of ∼18 molecules, or alternatively, it binds to the protein surface. The detergent binds in the form of a monolayer ring around the hydrophobic β-barrel of OmpX rather than in a micellar-like oblate; ∼40 dihexanoylphosphatidylcholine lipids are sufficient for an effective suppression of water from the surface of the β-barrel region. The phospholipids bind also on the extracellular, protruding β-sheet. Here, polar interactions between charged amino acids and phosphatidylcholine headgroups act as condensation seed for detergent micelle formation. The polar protein surface remains accessible to water molecules. In total, ∼90–100 detergent molecules associate within the protein-detergent mixed micelle, in agreement with experimental estimates. The simulation results indicate that OmpX is not a water pore and support the proposed role of the protruding β-sheet as a “fishing rod”

Membrane-mediated oligomerization of G protein coupled receptors and its implications for GPCR function

The dimerization or even oligomerization of G protein coupled receptors (GPCRs) causes ongoing, controversial debates about its functional role and the coupled biophysical, biochemical or biomedical implications. A continously growing number of studies hints to a relation between oligomerization and function of GPCRs and strengthens the assumption that receptor assembly plays a key role in the regulation of protein function. Additionally, progress in the structural analysis of GPCR-G protein and GPCR-ligand interactions allows to distinguish between actively functional and non-signalling complexes. Recent findings further suggest that the surrounding membrane, i.e. its lipid composition may modulate the preferred dimerization interface and as a result the abundance of distinct dimeric conformations. In this review, the association of GPCRs and the role of the membrane in oligomerization will be discussed. An overview of the different reported oligomeric interfaces is provided and their capability for signaling discussed. The currently available data is summarized with regard to the formation of GPCR oligomers, their structures and dependency on the membrane microenvironment as well as the coupling of oligomerization to receptor function