Soluble Regions of GlpG Influence Protein-Lipid Interactions

10 hours ago - The GlpG rhomboid protease from E. coli is a well-characterized intramembrane protease that cleaves transmembrane substrates inside the...
0 downloads 0 Views 849KB Size
Subscriber access provided by University of Glasgow Library

B: Biophysics; Physical Chemistry of Biological Systems and Biomolecules

Soluble Regions of GlpG Influence ProteinLipid Interactions and Lipid Distribution Yasser Almeida-Hernandez, and Henning Tidow J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.9b06943 • Publication Date (Web): 03 Sep 2019 Downloaded from pubs.acs.org on September 3, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Soluble Regions of GlpG Influence Protein-Lipid Interactions and Lipid Distribution Yasser Almeida-Hernandez1,2* and Henning Tidow1,2*

1

Department of Chemistry, Institute for Biochemistry and Molecular Biology, University of

Hamburg, Martin-Luther-King-Platz 6, D-20146 Hamburg, Germany 2

The Hamburg Centre for Ultrafast Imaging, University of Hamburg, Luruper Chaussee 149,

D-22761 Hamburg, Germany

* Corresponding authors Henning Tidow Department of Chemistry, Institute for Biochemistry and Molecular Biology, University of Hamburg, Martin-Luther-King-Platz 6, D-20146 Hamburg, Germany Tel: +49 40428388984 e-mail: [email protected] Yasser Almeida-Hernandez Department of Chemistry, Institute for Biochemistry and Molecular Biology, University of Hamburg, Martin-Luther-King-Platz 6, D-20146 Hamburg, Germany e-mail: [email protected]

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 21

Abstract The GlpG rhomboid protease from E. coli is a well-characterized intramembrane protease that cleaves transmembrane substrates inside the lipid bilayer. Most studies have focused on the GlpG transmembrane domain containing the catalytic site, while the full-length protein, also containing a soluble cytoplasmic domain, a linker region, and a small positively charged Cterminal fragment remains poorly understood. In this work, we used coarse-grained molecular dynamics (CGMD) simulations to investigate full-length GlpG embedded in a native-like model of the E. coli membrane. We identified differences in the distribution and clustering of phosphoglycerol(PG)-based lipids around GlpG in both leaflets depending on whether the soluble regions are present or absent. These data suggest a possible role of the cytoplasmic extensions of GlpG in the regulation of the lipid environment around GlpG, which may influence the activity of GlpG in vivo.

2 ACS Paragon Plus Environment

Page 3 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Introduction Intramembrane proteases are present across all species and are involved in key physiologic functions such as cell signaling, pathogen invasion, cell growth, and epithelial homeostasis 1– 6.

These proteins perform their proteolytic function by cleaving peptide bonds inside the

hydrophobic environment of the lipid bilayer 7,8. How these enzymes perform their hydrolytic reaction in a lipidic environment has been investigated for 20 years but still remains poorly understood at a molecular level. The GlpG rhomboid protease from E. coli is one of the most extensively studied intramembrane proteases. GlpG is a serine protease with Ser201 and His254 forming the catalytic site. The protein contains a transmembrane domain (TMD) consisting of six αhelices, a linker region (Ln) and a soluble cytosolic N-terminal domain (CytoD) 9. Several crystal structures of wild-type and mutants of the TMD have been determined, either in apo form and in complex with inhibitors

10–17,

providing detailed structural information about the

catalytic site and its arrangement within the TMD. As the TMD shows activity in the absence of the rest of the protein 18, the influence of the Ln and the CytoD domain have been widely ignored. The atomic structure of the CytoD has been determined by NMR and X-ray crystallography, showing a compact α-β fold that can also form domain-swapped dimers in the presence of detergents

18,19.

The function of CytoD is still unknown. It has been shown

that cytosolic extensions of rhomboid protease Rho4 from Drosophila melanogaster regulate the activity by modulating the substrate gating by calcium binding

20,

while in yeast, the

cytosolic domain of Rbd2 rhomboid protease controls the lipids organization and actin assembly during clathrin-mediated endocytosis 21. In the course of the standard purification of TMD or CytoD by protease cleavage of full-length GlpG 10 the Ln region is eliminated; hence there is little structural information about this region, although it is expected to be flexible and

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 21

unstructured. In addition, GlpG contains a small positively charged, C-terminal fragment, that it is eliminated during the purification of TMD. Using atomistic MD simulations of GlpG TMD inserted in POPE or POPC model membranes, differences in the hydrogen bonding between lipid head groups and the protein surface were found, which influence the orientation of GlpG in the bilayer 22. More recently, full-length GlpG was reconstituted in SMALP nanodiscs to study the influence of the native membrane lipids composition in the dynamics and accessibility of the protein using hydrogendeuterium exchange mass spectrometry (HDX-MS) 23. In this study, we used the previously determined lipid ratios in E. coli membranes

23

to build

systems of full-length GlpG (GlpG) and GlpG-TMD embedded in a model membrane composed of PE-, PG-derived lipids, and cardiolipin. Using coarse-grained molecular dynamics (CGMD) simulations we identified GlpG/lipid interaction sites and distinct lipid clustering in both leaflets, with differences between GlpG and GlpG-TMD. Our simulations revealed specific accumulation features of various lipids (differing with respect to saturation level and head group) to GlpG. While some lipids accumulate around the protein (annular lipids) and show lower presence at higher distances from the protein (e.g. DOPG), other lipids are less likely to interact with the protein and are more homogeneously distributed in the bilayer (bulk, DPPG).

Methods Systems setup A structural model for full-length GlpG was generated using the Phyre server 24. This model and a GlpG-TMD structure (PDB: 2XOV) were converted to CG models using the program martinize

25.

The proteins were inserted with the program insane

26

in 12 x 12 nm solvated

model membranes, containing PE- (DPPE: di-C16:0 PE, DOPE: di-C18:1(9c) PE, POPE: 4 ACS Paragon Plus Environment

Page 5 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

C16:0/18:1 PE), PG- (DPPG: di-C16:0 PE, di-C18:1(9c) PG, POPG: C16:0/18:1 PG) based lipids, and cardiolipin, at a ratio of 70.5 % / 27 % / 2.5 %, respectively (Table 1). Lipids topologies and coordinates were obtained from the MARTINI force field web site 27. Na+ and Cl- were added at concentration of 0.15 M for neutralization. 5 independent systems were built for each protein, where the bilayers were constructed placing the lipids randomly.

Simulations The simulations were performed with the GROMACS suite 5.1.2 28, using the MARTINI 2.2 force field

29.

Elastic network RubberBand, similar to ElNeDyn

30,

and implemented in

martinize, was used to restrain the entire protein during the simulations. We used all the default parameters in Martini force field defined in the martinize script. We built the ELNEDYN elastic network CG system with the default cut-off distance RC = 0.9 nm and the elastic bond constant KSPRING = 500 kJ mol-1 nm-2 30. We used the default reaction-field with the parametrized short-range shifted electrostatic interaction, with Coulomb cut-off distance rCoulomb = 1.2 nm. The systems were minimized with the steepest descent method up to a maximum of 10000 steps. Next, equilibration runs were performed using a Berendsen barostat 31

with a coupling time of 10 ps. Independent production runs for each system were executed

during 10 s using semi-isotropic pressure coupling to a reference pressure of 1.0 bar with the Parrinello-Rahman barostat

32

and compressibility of 3.4 x 10-4 bar-1. The temperature was

controlled at 300 K using the velocity rescaling thermostat

33

with a time step of 10 fs. Five

independently built trajectories were simulated for each system. Analysis In order to allow sufficient equilibration of the systems, all analyses were performed on the last 1 s of the simulations. Before the analysis, all the trajectories were structurally aligned

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 21

to the last frame of the GlpG in the first run. XY radial distribution functions (RDF) were computed with the xy_rdf tool implemented in the LOOS suite (v2.3.2)

34

using the center of

mass of the transmembrane domain in both systems. 2D density maps were obtained with the densmap tool from GROMACS, with a grid size of 1 Å. Proteins were rotationally aligned and used as reference to obtain the maps in the same coordinate system. Protein/PO4 bead contacts were calculated with the tool select (implemented in GROMACS) and mapped onto the CG structure with in house scripts. Contacts were defined as protein backbone BB beads and lipid PO4 beads at a cutoff distance of a smaller than 6 Å. Structures and trajectories were visualized and analyzed with UCSF Chimera 35 and VMD 36.

Results and Discussion As the transmembrane domain (TMD) of GlpG retains proteolytic activity most previous studies on GlpG focused on the TMD, and only a few studies of the full-length protein are available. We performed coarse-grained molecular dynamics simulations of full-length GlpG as well as GlpG-TMD embedded in a model membrane composed of PE-, PG-derived lipids, and cardiolipin (Figure 1, Figure S1). The investigated lipids differ with respect to their head group and to the number of unsaturations (double-bonds) in their acyl chains. We computed XY radial distribution functions (RDF) for the phosphate group (coarse-grained PO4 bead) for each phospholipid, in the plane of both leaflets, for GlpG and GlpG-TMD (Figure 2). RDFs show clear differences between lipids with respect to the interaction with GlpG, depending on the presence of CytoD. DOPG dominates the interaction with the protein in both leaflets and systems, although this effect is more prominent in the bottom leaflet. POPG also shows a strong interaction with GlpG in the bottom leaflet, which is lost when CytoD is eliminated. The RDF plot for GlpG-TMD in the bottom leaflet shows that almost all

6 ACS Paragon Plus Environment

Page 7 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

lipids (in particular DOPG) are located closer to the center of mass of the transmembrane domain, filling the space that was occupied by the CytoD in GlpG. In order to study the spatial distribution of the lipids in the XY plane of the membrane surrounding GlpG and GlpG-TMD, we computed 2D density maps, averaged along the last 1 s of the simulations for PE- and PG-based lipids in both leaflets (Figure 3 and S2). Each pixel of the maps corresponds to a spatial bin of 1 Å2 and contains the number of events of the presence of a certain PO4 bead in a XY coordinate bin, averaged among all simulations for each system. The inspection of the maps revealed interesting features about the interaction of the lipids with the protein in each leaflet, in particular for POPG and DOPG (Figure 2). In the bottom leaflet, the absence of CytoD affects the protein/lipid interactions, in particular for POPG. We observed a clustering of POPG to two regions of full-length GlpG (Figure 3). The first region corresponds to the residues G170, S171, G172 and A272 (TMD), and the second to residues R81, R82 (Ln) and R227 (TMD). These two clusters were affected when the soluble region was eliminated. In particular, the second cluster is formed by the interaction of Arg residues from the TMD and the Ln region, which provides a positively charged, responsible for the clustering of the lipid in this area. However, when the soluble region is removed, the clustering and number of contacts between POPG and R227 were significantly reduced (Figure 4). For DOPG, there was also a strong clustering in the bottom leaflet in GlpG, in a region corresponding to the Ln segment, as well as an accumulation around the protein forming an annulus. This clustering was diminished in the area corresponding to the Ln segment. This indicated that, although residues located in the TMD seem to be sufficient for the lipid-interaction, the presence of certain residues in the cytoplasmic domain and the linker region provides further stabilization to retain the protein/lipid interaction. In this respect, the

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 21

amount of positively charges residues located in Ln and the region of TMD facing the bottom leaflet may play a role. The differences between POPG and DOPG in the bottom leaflet of GlpG-TMD can be attributed to unsaturations and the hydrophobic mismatch effect. It has been reported that unsaturated lipids show stronger interactions with several membrane proteins through this effect, which can subsequently affect their activity

37–41.

Compared to DPPG, which is

completely saturated, POPG and DOPG both show increased interactions with the protein (Figures 3 and S2). In addition, the two double bonds of DOPG make the length of the acyl chains shorter compared to POPG (contains one double bond), leading to membrane thinning. This behaviour of DOPG could make it more prone to interact with GlpG by matching the hydrophobic surface area of GlpG forming the before mentioned annulus, which is thinner than the average thickness of the membrane (hydrophobic mismatch)

42.

POPG, however,

with only one unsaturation can be more extended and interacts with the protein in defined residues. The hydrophobic mismatch effect on GlpG and its impact on the protein diffusion in the membrane were recently also experimentally demonstrated by Kreutzberger and coworkers 43. Interestingly, we also detected some changes in POPG and DOPG clustering in the top leaflet (Figure 3), in a region corresponding to an interaction with the L5 loop, which is responsible for the protein gating to allow substrates to enter in the catalytic site. This effect was also evident in the number of protein/PO4 bead contacts (Figures 4,5 and S3). Previous studies have shown the importance of loop L5 in vitro and in vivo for the gating during substrate binding, and mutations in this region completely abolished the activity 44–46. These variations in the lipid clustering and contacts around the L5 loop may suggest that these could regulate the access of substrates to the active site and hence the activity of the protein. Similar behavior has recently been described in atomistic simulations, where lipid head groups could

8 ACS Paragon Plus Environment

Page 9 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

compete with the substrate in the membrane by transient contacts in the active site of the protein

47.

Furthermore, crystal structures of GlpG-TMD show a detergent molecule in the

active site as well as lipid molecules interacting with residues in the cytoplasmatic side of the molecule 48. We also found that the positively charged C-terminal fragment interacts with all lipids (Figure 5). This fragment together with CytoD and Ln is usually lost during the standard proteolytic purification of TMD

10

and thus not present in the available crystal

structures. In summary, our simulations reveal that GlpG preferentially interacts with unsaturated and PG-based lipids. These lipids preferentially accumulate in the annular belt of GlpG, which may lead to membrane thinning due to hydrophobic mismatch 42. In addition, residues in the linker and bottom region of GlpG are involved in the regulation of the lipid distribution in the bottom leaflet of the membrane. Although our simulations did not show large conformational changes for the linker region, it may still be flexible and adopt extended conformations, in which the cytoplasmic domain and the linker region are more distant from the membrane, resulting in an arrangement resembling the GlpG-TMD simulated in this study. In that case, the interactions with the lipids in the bottom leaflet would change, affecting their distribution in both leaflets. We propose this as a putative a mechanism to regulate the activity of GlpG, where the soluble regions of GlpG act as a relay by regulating the lipid distribution around the protein. As stated before, this mechanism has been demonstrated before in yeast, where the cytosolic domain of Rbd2 rhomboid protease influences the lipids organization, which in turn regulates actin assembly during clathrin-mediated endocytosis 21. Our results indicate that the cytoplasmatic extensions of GlpG affect the lipidic environment around its transmembrane domain.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 21

Conclusion We performed coarse-grained molecular dynamics simulations to study the interaction of GlpG with lipids of the cell membrane. Although most structural and functional studies of GlpG to date have been performed using the transmembrane domain only, our data indicate that cytoplasmic regions of GlpG (N-terminal cytoplasmic domain, linker region and Cterminal extension) play a role in the regulation of the activity of the protein by modifying the surrounding lipid environment. In addition, we show that unsaturated PG-based lipids interact preferentially with GlpG compared to PE-based and saturated lipids. Although further experimental work is required to prove this hypothesis, our simulations contribute to a better understanding of the possible regulatory role of the soluble fragments of GlpG, widely ignored to date.

Supporting Information Available 

Figure S1: Top and bottom views of GlpG and GlpG-TMD



Figure S2: Density maps of PE-based lipids and DPPG for GlpG and GlpG-TMD



Figure S3: GlpG/GlpG-TMD contacts with PO4 beads along the simulation

Author contributions Conceptualization, YAH and HT; Investigation, YAH; Analysis and Interpretation, YAH and HT; Manuscript Writing, YAH and HT; Funding Acquisition, HT; Supervision, HT.

Acknowledgments We thank Regionales Rechenzentrum (RRZ) at the University of Hamburg for access to the Hummel computer cluster. This research was funded by the excellence cluster 'The Hamburg

10 ACS Paragon Plus Environment

Page 11 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Centre for Ultrafast Imaging - Structure, Dynamics and Control of Matter at the Atomic Scale' of the Deutsche Forschungsgemeinschaft (DFG EXC 1074).

Competing Financial Interest statement We declare no competing financial interests.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 21

References (1)

Urban, S.; Lee, J. R.; Freeman, M. Drosophila Rhomboid-1 Defines a Family of Putative Intramembrane Serine Proteases. Cell 2001, 107 (2), 173–182.

(2)

Urban, S. Making the Cut: Central Roles of Intramembrane Proteolysis in Pathogenic Microorganisms. Nat. Rev. Microbiol. 2009, 7 (6), 411–423.

(3)

Baker, R. P.; Wijetilaka, R.; Urban, S. Two Plasmodium Rhomboid Proteases Preferentially Cleave Different Adhesins Implicated in All Invasive Stages of Malaria. PLoS Pathog. 2006, 2 (10), 0922–0932.

(4)

Brossier, F.; Jewett, T. J.; Sibley, L. D.; Urban, S. A Spatially Localized Rhomboid Protease Cleaves Cell Surface Adhesins Essential for Invasion by Toxoplasma. Proc. Natl. Acad. Sci. 2005, 102 (11), 4146–4151.

(5)

Johnson, N.; Březinová, J.; Stephens, E.; Burbridge, E.; Freeman, M.; Adrain, C.; Strisovsky, K. Quantitative Proteomics Screen Identifies a Substrate Repertoire of Rhomboid Protease RHBDL2 in Human Cells and Implicates It in Epithelial Homeostasis. Sci. Rep. 2017, 7 (1), 7283.

(6)

Urban, S.; Dickey, S. W. The Rhomboid Protease Family: A Decade of Progress on Function and Mechanism. Genome Biol. 2011, 12 (10).

(7)

Brown, M. S.; Ye, J.; Rawson, R. B.; Goldstein, J. L. Regulated Intramembrane Proteolysis: A Control Mechanism Conserved from Bacteria to Humans. Cell. 2000, 100 (4)391–398.

(8)

Urban, S. Taking the Plunge: Integrating Structural, Enzymatic and Computational Insights into a Unified Model for Membrane-Immersed Rhomboid Proteolysis. Biochem. J. 2010, 425, 501–512.

(9)

Cournia, Z.; Allen, T. W.; Andricioaei, I.; Antonny, B.; Baum, D.; Brannigan, G.; Buchete, N.-V.; Deckman, J. T.; Delemotte, L.; del Val, C.; et al. Membrane Protein Structure, Function, and Dynamics: A Perspective from Experiments and Theory. J. Membr. Biol. 2015, 248 (4), 611–640.

(10)

Wang, Y.; Zhang, Y.; Ha, Y. Crystal Structure of a Rhomboid Family Intramembrane Protease. Nature 2006, 444 (7116), 179–183.

(11)

Wu, Z.; Yan, N.; Feng, L.; Oberstein, A.; Yan, H.; Baker, R. P.; Gu, L.; Jeffrey, P. D.; Urban, S.; Shi, Y. Structural Analysis of a Rhomboid Family Intramembrane Protease Reveals a Gating Mechanism for Substrate Entry. Nat. Struct. Mol. Biol. 2006, 13 (12), 1084–1091.

(12)

Ben-Shem, A.; Fass, D.; Bibi, E. Structural Basis for Intramembrane Proteolysis by Rhomboid Serine Proteases. Proc. Natl. Acad. Sci. 2007, 104 (2), 462–466.

12 ACS Paragon Plus Environment

Page 13 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(13)

Wang, Y.; Ha, Y. Open-Cap Conformation of Intramembrane Protease GlpG. Proc. Natl. Acad. Sci. 2007, 104 (7), 2098–2102.

(14)

Xue, Y.; Chowdhury, S.; Liu, X.; Akiyama, Y.; Ellman, J.; Ha, Y. Conformational Change in Rhomboid Protease GlpG Induced by Inhibitor Binding to Its S′ Subsites. Biochemistry 2012, 51 (18), 3723–3731.

(15)

Zoll, S.; Stanchev, S.; Began, J.; Škerle, J.; Lepšík, M.; Peclinovská, L.; Majer, P.; Strisovsky, K. Substrate Binding and Specificity of Rhomboid Intramembrane Protease Revealed by Substrate–Peptide Complex Structures. EMBO J. 2014, 33 (20), 2408– 2421.

(16)

Cho, S.; Dickey, S. W.; Urban, S. Crystal Structures and Inhibition Kinetics Reveal a Two-Stage Catalytic Mechanism with Drug Design Implications for Rhomboid Proteolysis. Mol. Cell 2016, 61 (3), 329–340.

(17)

Tichá, A.; Stanchev, S.; Vinothkumar, K. R.; Mikles, D. C.; Pachl, P.; Began, J.; Škerle, J.; Švehlová, K.; Nguyen, M. T. N.; Verhelst, S. H. L.; et al. General and Modular Strategy for Designing Potent, Selective, and Pharmacologically Compliant Inhibitors of Rhomboid Proteases. Cell Chem. Biol. 2017, 24 (12), 1523-1536.e4.

(18)

Lazareno-Saez, C.; Arutyunova, E.; Coquelle, N.; Lemieux, M. J. Domain Swapping in the Cytoplasmic Domain of the Escherichia Coli Rhomboid Protease. J. Mol. Biol. 2013, 425 (7), 1127–1142.

(19)

Sherratt, A. R.; Blais, D. R.; Ghasriani, H.; Pezacki, J. P.; Goto, N. K. Activity-Based Protein Profiling of the Escherichia Coli GlpG Rhomboid Protein Delineates the Catalytic Core. Biochemistry 2012, 51 (39), 7794–7803.

(20)

Baker, R. P.; Urban, S. Cytosolic Extensions Directly Regulate a Rhomboid Protease by Modulating Substrate Gating. Nature 2015, 523 (7558), 101–105.

(21)

Cortesio, C. L.; Lewellyn, E. B.; Drubin, D. G. Control of Lipid Organization and Actin Assembly during Clathrin-Mediated Endocytosis by the Cytoplasmic Tail of the Rhomboid Protein Rbd2. Mol. Biol. Cell 2015, 26 (8), 1509–1522.

(22)

Bondar, A.-N.; del Val, C.; White, S. H. Rhomboid Protease Dynamics and Lipid Interactions. Structure 2009, 17 (3), 395–405.

(23)

Reading, E.; Hall, Z.; Martens, C.; Haghighi, T.; Findlay, H.; Ahdash, Z.; Politis, A.; Booth, P. J. Interrogating Membrane Protein Conformational Dynamics within Native Lipid Compositions. Angew. Chemie - Int. Ed. 2017, 56 (49), 15654–15657.

(24)

Kelley, L. A; Mezulis, S.; Yates, C. M.; Wass, M. N.; Sternberg, M. J. E. The Phyre2 Web Portal for Protein Modeling, Prediction and Analysis. Nat. Protoc. 2015, 10 (6), 845–858.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 21

(25)

de Jong, D. H.; Singh, G.; Bennett, W. F. D.; Arnarez, C.; Wassenaar, T. a.; Schäfer, L. V.; Periole, X.; Tieleman, D. P.; Marrink, S. J. Improved Parameters for the Martini Coarse-Grained Protein Force Field. J. Chem. Theory Comput. 2013, 9 (1), 687–697.

(26)

Wassenaar, T. a.; Ingólfsson, H. I.; Böckmann, R. a; Tieleman, D. P.; Marrink, S. J. Computational Lipidomics with Insane: A Versatile Tool for Generating Custom Membranes for Molecular Simulations. J. Chem. Theory Comput. 2015, 11, 2144−2155.

(27)

Marrink, S. J. Coarse Grain Forcefield for Biomolecules http://cgmartini.nl/index.php/force-field-parameters/lipids/ (accessed August 8, 2017)

(28)

Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindah, E. Gromacs: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1–2, 19–25.

(29)

Ravikumar, K. M.; Huang, W.; Yang, S. Coarse-Grained Simulations of ProteinProtein Association: An Energy Landscape Perspective. Biophys. J. 2012, 103 (4), 837–845.

(30)

Periole, X.; Cavalli, M.; Marrink, S. J.; Ceruso, M. a. Combining an Elastic Network with a Coarse-Grained Molecular Force Field: Structure, Dynamics, and Intermolecular Recognition. J. Chem. Theory Comput. 2009, 5 (9), 2531–2543.

(31)

Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, a; Haak, J. R. Molecular Dynamics with Coupling to an External Bath. J. Chem. Phys. 1984, 81 (8), 3684–3690.

(32)

Parrinello, M. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52 (12), 7182.

(33)

Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126 (1), 1–7.

(34)

Romo, T. D.; Leioatts, N.; Grossfield, A. Lightweight Object Oriented Structure Analysis: Tools for Building Tools to Analyze Molecular Dynamics Simulations. J. Comput. Chem. 2014, 35 (32), 2305–2318.

(35)

Pettersen, E. F.; Goddard, T. D.; Huang, C. C.; Couch, G. S.; Greenblatt, D. M.; Meng, E. C.; Ferrin, T. E. UCSF Chimera - A Visualization System for Exploratory Research and Analysis. J. Comput. Chem. 2004, 25 (13), 1605–1612.

(36)

Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 1996, 14 (18), 33–38.

(37)

Yazdi, S.; Stein, M.; Elinder, F.; Andersson, M.; Lindahl, E. The Molecular Basis of Polyunsaturated Fatty Acid Interactions with the Shaker Voltage-Gated Potassium Channel. PLOS Comput. Biol. 2016, 12 (1), e1004704.

14 ACS Paragon Plus Environment

Page 15 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(38)

Grossfield, A.; Feller, S. E.; Pitman, M. C. A Role for Direct Interactions in the Modulation of Rhodopsin by Omega-3 Polyunsaturated Lipids. Proc. Natl. Acad. Sci. U. S. A. 2006, 103 (13), 4888–4893.

(39)

Niu, S.-L.; Mitchell, D. C.; Litman, B. J. Optimization of Receptor-G Protein Coupling by Bilayer Lipid Composition II. J. Biol. Chem. 2001, 276 (46), 42807–42811.

(40)

Oliver, D.; Lien, C.-C.; Soom, M.; Baukrowitz, T.; Jonas, P.; Fakler, B. Functional Conversion between A-Type and Delayed Rectifier K+ Channels by Membrane Lipids. Science 2004, 304 (5668), 265–270.

(41)

Pitman, M. C.; Grossfield, A.; Suits, F.; Feller, S. E. Role of Cholesterol and Polyunsaturated Chains in Lipid−Protein Interactions: Molecular Dynamics Simulation of Rhodopsin in a Realistic Membrane Environment. J. Am. Chem. Soc. 2005, 127 (13), 4576–4577.

(42)

Strisovsky, K. Structural and Mechanistic Principles of Intramembrane Proteolysis-Lessons from Rhomboids. FEBS J. 2013, 280 (7), 1579–1603.

(43)

Kreutzberger, A. J. B.; Ji, M.; Aaron, J.; Mihaljević, L.; Urban, S. Rhomboid Distorts Lipids to Break the Viscosity-Imposed Speed Limit of Membrane Diffusion. Science. 2019, 363 (6426), eaao0076.

(44)

Baker, R. P.; Young, K.; Feng, L.; Shi, Y.; Urban, S. Enzymatic Analysis of a Rhomboid Intramembrane Protease Implicates Transmembrane Helix 5 as the Lateral Substrate Gate. Proc. Natl. Acad. Sci. 2007, 104 (20), 8257–8262.

(45)

Urban, S.; Baker, R. P. In Vivo Analysis Reveals Substrate-Gating Mutants of a Rhomboid Intramembrane Protease Display Increased Activity in Living Cells. Biol. Chem. 2008, 389 (8), 1107–1115.

(46)

Xue, Y.; Ha, Y. Large Lateral Movement of Transmembrane Helix S5 Is Not Required for Substrate Access to the Active Site of Rhomboid Intramembrane Protease. J. Biol. Chem. 2013, 288 (23), 16645–16654.

(47)

Bondar, A.-N. Mechanisms by Which Lipids Influence Conformational Dynamics of the GlpG Intramembrane Protease. J. Phys. Chem. B 2019, 114 (3), 4159-4172.

(48) Vinothkumar, K. R. Structure of Rhomboid Protease in a Lipid Environment. J. Mol. Biol. 2011, 407 (2), 232–247.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 21

Figure legends Figure 1

Full-length model of GlpG embedded in a phospholipid bilayer. A) Coarse-grained (CG) representation of GlpG embedded in the bilayer. Regions of the protein are shown: CytoD (magenta), Ln (cyan), TMD (green), L5 loop (blue) and C-terminal fragment (red). The model was generated with the Phyre server using the crystal structure of the transmembrane domain (pdb:2XOV) and the NMR-structure of the cytoplasmic domain (pdb:2LEP), and transformed with martinize

25.

B) CG models of lipids used with the beads for every atoms group:

ethanolamine (blue), PO4 (orange), glycerol (magenta), acyl chain (gray). The white bead represents the unsaturation of the acyl chain.

16 ACS Paragon Plus Environment

Page 17 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 2

Lateral radial distribution function (XY-RDF) analysis. Displayed is the XY-RDF of the lipid phosphate beads in the top and bottom leaflets of the model bilayer, with embedded GlpG and GlpG-TMD. Curves for each lipid represent the average between five independent runs.

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 21

Figure 3

2D density maps of the phosphate beads of POPG and DOPG, with embedded GlpG and GlpG-TMD. Each pixel of the maps represents a bin of 1 Å2 and contains the count of PO4 beads in each grid point along the last s of the simulations averaged between five independent runs. The red arrowheads (top leaflet) indicate the lipid accumulation close to the L5 loop and the substrate-gating site. Dark grey arrowhead (bottom leaflet, POPG) indicates an interaction cluster with G170, S171, G172 and A272 (TMD). Light arrow indicates an interaction cluster with R81, R82 (Ln) and R227 (TMD).

18 ACS Paragon Plus Environment

Page 19 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 4

Residue-based protein-lipid contacts. Number of contacts between backbone BB beads of the protein and PO4 beads of each lipid. Contacts were defined as the PO4 at a cutoff distance of smaller than 6 Å. The colored bars represent the regions in the protein. Magenta: CytoD, Cyan: Ln region, Green: TMD, Blue: L5 loop.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 21

Figure 5

3D mapping of the protein-lipid (backbone BB bead/lipid PO4 bead) contacts on the structure. The color scale and the size of the BB beads represent the contact frequency of each lipid PO4 bead with the protein residues, normalized against the total number of contacts of each lipid. The color of the bonds represents each region of the protein. Magenta: CytoD, Cyan: Ln region, Green: TMD, Blue: L5 loop.

20 ACS Paragon Plus Environment

Page 21 of 21 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Table 1: Simulations setups GlpG full length

GlpG-TMD

Lipid

TOP

BOTTOM

TOP

BOTTOM

DPPE

47

48

49

48

POPE

47

48

49

48

DOPE

47

48

49

48

DPPG

18

18

18

18

DOPG

18

18

18

18

POPG

18

18

18

18

CDL

3

3

3

3

Solvent Water

8966-9036

9224-9239

Na+

133

138

Cl-

66

68

TOC Graphic

ACS Paragon Plus Environment

21