Molecular Details of the PH Domain of ACAP1BAR-PH Protein

Subscriber access provided by UNIV OF CALIFORNIA SAN DIEGO LIBRARIES

Article BAR-

Molecular Details of the PH Domain of ACAP1 Protein Binding to PIP-Containing Membrane PH

Chun Chan, Lanyuan Lu, Fei Sun, and Jun Fan J. Phys. Chem. B, Just Accepted Manuscript • DOI: 10.1021/acs.jpcb.6b09563 • Publication Date (Web): 16 Jan 2017 Downloaded from http://pubs.acs.org on January 17, 2017

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry B 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 48

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

Molecular Details of the PH domain of ACAP1BAR-PH Protein Binding to PIP-containing Membrane Kevin Chun Chan1, Lanyuan Lu2, Fei Sun3,4, Jun Fan*,1,5,6 1

Department of Physics and Materials Science, City University of Hong Kong, Hong Kong, China;

2

School of Biological Sciences, Nanyang Technological University, Singapore;

3

National Laboratory of Biomacromolecules, Institute of Biophysics, Chinese Academy of Sciences, Beijing 100101, China;

4

Center for Biological Imaging, Institute of Biophysics, Chinese Academy of Sciences, Beijing 100101, China;

5

City University of Hong Kong Shenzhen Research Institute, Shenzhen 518057, China

6

Center for Advanced Nuclear Safety and Sustainable Development, City University of Hong Kong, Hong Kong, China;

1 ACS Paragon Plus Environment


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

ABSTRACT ACAP1 proteins were previously reported to specifically bind PIP2-containing cell membranes and form well-structured protein lattices in order to conduct membrane tubulation. We carried out molecular dynamics simulations to characterize orientation of the PH domains with respect to the BAR domains inside the protein dimer. Followed by molecular dynamics simulations, we present a comprehensive orientation analysis of PH domain under different states including unbound and bound with lipids. We have examined two binding pockets on the PH domain and present PMF profiles of the two pockets to account for their preference to PIP2 lipids. Combining orientation analysis and studies of binding pockets, our simulations results reveal valuable molecular basis for protein-lipid interactions of ACAP1 proteins during membrane remodeling process.


Page 2 of 48

Page 3 of 48

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


INTRODUCTION

Understanding the association of peripheral proteins to cell membranes is critical for many aspects of cellular function1,2. One of the fundamental questions would be how proteins recognize and bind to specific lipid molecules3-6. Addressing this question requires knowledge of structure, energetics, and dynamics of the complex protein-lipid systems. Molecular dynamics (MD) simulations have been widely employed as one of the best tools for studying such complex biological systems on both temporal and spatial scales not accessible to current experimental methods7-12. MD simulations provide not only molecular details of the structures and dynamics of the protein-lipid complexes, but also valuable energetics information when combining with sophisticated free-energy calculation methodologies13-18. Arfgap with Coiled coil, Ankyrin repeat, and PH domain protein 1 (ACAP1) was initially characterized as a GTPase-activating protein (GAP) for ADP-ribosylation factor 6 (ARF6)19. Subsequently, it was found to also involve in endocytic recycling as an ARF6 effector20,21. Being a coat component, four domains of ACAP1 were found to be responsible for different functions during endocytosis: ArfGAP and ankryin repeat (ANK) domains have been well studied for its cargo binding abilities22; Bin-Amphiphysin-Rvs (BAR) and Pleckstrin Homology (PH) domains were found to oligomerize and promote membrane curvatures in order to generate transport carriers. BAR domain is one of the best studied protein domain motif related to membrane curvature recognition and generation23-31. In previous studies, however, we found that the PH domain of ACAP1 protein dimer take over the role of the BAR domain when binding to cell membranes32. The PH domain was first identified in 1993 and has been found in numerous proteins involved in cellular signaling33,34. The view that the PH domain directs membrane targeting of their host proteins by binding to phosphoinositides (PIs) has been established since then35-37. PIs play crucial roles in both cell signaling and membrane traffic38-41. Phosphatidylinositol-4,5-bisphosphate (PI(4,5)P2 or more briefly PIP2), being the best studied species among eukaryotic PIs, were reported to be bound by a significant number of known PH domains with high affinity and specificity42,43. In



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

addition to PIs, PH domains were also reported to bind to other proteins in signal transduction pathways44,45. Multiple functions of the PH domain revealed by experiments initiate a number of computational studies9-12,46. While most of the previous studies were restricted to the PH domain of GRP1 protein, we recently performed experimental and computational studies on BAR-PH domain of ACAP1 protein (ACAP1BAR-PH) to address novel features of the PH domain over-taking the role of the BAR domain during membrane curvature generation32. In current studies, we further investigate the molecular details of the protein-lipid interactions of ACAP1BAR-PH protein with the assistance of both restrained and unbiased atomistic molecular dynamics (MD) simulations and sophisticated freeenergy calculation methodologies.


Page 4 of 48

Page 5 of 48

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


METHODS Molecular dynamics simulations. A single ACAP1BAR-PH dimer with and without lipid bilayers underneath were constructed in a similar manner. Initial structures of the protein were based on the crystal structure of ACAP1BAR-PH protein (PDB: 4NSW). The lipid

bilayer

was

constructed

using

CHARMM-GUI47

with

a

composition

of

DOPC/DOPE/POPS/PI(4,5)P2 (4:3:2:1) and contained 600 lipid molecules in total. Details of the membrane construction can be referred to our previous work32. The protein was placed on the top of the membrane for the protein-membrane system. Orientations of the protein were chosen so that the concave surface of the BAR domain faces the underneath membrane. The distance between the protein and the membrane was chosen to be 10 A, within the cutoff of the non-bonded interactions. The systems were then solvated in explicit TIP3P water molecules48, with potassium and chloride ions included at a final concentration of 0.18 M using the solvate and auto-ionize tools in VMD49. Periodic boundary conditions were introduced to both systems. Overall, three protein-solvent simulations (210x100x110 Å3, ~220k atoms, simulation 1-3) and two protein-lipid simulations (270x105x155 Å3, ~300k atoms, simulation 4-5) have been performed. The MD simulations were performed using NAMD 2.11 package50. CHARMM 22/27 force field51 with CMAP correction52,53 was used for proteins and CHARMM 36 force field54 was used for lipids. Electrostatic interactions were calculated using the particle mesh Ewald sum method55 with a cutoff of 12 Å. All hydrogen-containing covalent bonds were constrained by the SHAKE algorithm56, therefore allowing an integration time step of 2 fs. Initially, the system was minimized using conjugate gradient algorithm followed by an equilibrium in canonical ensemble at 310K. Protein backbone and oxygen atoms of water were harmonically restrained by a spring constant (1 kcal mol-1 Å-2). Simulations were then continued in the constant NPT ensemble (310K and 1 atm). Langevin thermostats with a damping coefficient of 0.5 ps-1 were used to control the system temperature. A Langevin-piston barostat57 with a piston period of 2 ps and a damping time of 2 ps was used to control the pressure. Anisotropic fluctuations were allowed for the protein-lipid systems. Constraints were next released step-wise 5 ACS Paragon Plus Environment


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

(with spring constant gradually decrease from 1 kcal mol-1 Å2 to 0 by steps of 0.1 kcal mol-1 Å2) before starting the production runs. A total of 420 ns and 2 x 380 ns of trajectory were generated for the proteinsolvent systems respectively. A total of 420 ns and 250 ns of trajectories were generated for the two protein-lipid systems respectively. In both simulations, we have observed that only one of the PH domains bound to the membrane underneath, i.e. the asymmetrical binding event. Equilibrium of the systems were judged by examining a number of thermodynamic quantities, such as simulation box size, energy, temperature, pressure, and Root-mean-square deviation (RMSD) of the protein backbone. All the quantities reached a plateau after equilibrium, while the RMSD profile being the latest one (Fig. S1). Only data after the system were equilibrated were taken for further analysis. For the protein-solvent systems, the first 180 ns of trajectories were excluded for analysis. For the 420 ns and 250 ns protein-lipid system, the first 220 and 150 ns of trajectories were excluded for analysis. All quantities presented in this article are averaged values over the equilibrated simulation windows. Protein-lipid contacts were defined as a cutoff distance of 12 Å between the center of mass (COM) of a residue and any lipid backbone phosphate. Such distance was chosen according to the cutoff distance of the non-bonded interactions.

Umbrella sampling simulations. The PH domain of the ACAP1BAR-PH structure from the last frame of the proteinlipid simulation 1 served as the initial PIP2-bound structure for umbrella sampling. As we have PIP2 lipid found at both binding pockets on the PH domain, we generated two structures, in which PIP2 lipid bound to only one of the pockets, by replacing one of the PIP2 with a DOPC molecule. We then placed the new PIP2-bound structure inside DOPC lipid bilayers, so that the bound PIP2 is the only PIP2 lipid in the membrane. In this way, we may remove the effects of surrounding charged lipids. POPS-bound structures were prepared in a similar manner by replacing the PIP2 lipid inside the desired pocket into a POPS lipid. Using a steered molecular dynamics (SMD) simulation, the PH domain was pulled away from the bound PIP2 lipid by fixing the COM to a reference point (with restraint constant k=1000 kJ mol-1 nm-2) which was moved along 6 ACS Paragon Plus Environment

Page 6 of 48

Page 7 of 48

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 membrane normal at a rate of 1 Å ns-1 for a total distance of 18 Å. In a SMD simulation, the orientation of the PH domain and the PIP2 backbone phosphate atom was restrained relative to their initial orientation or position using collective variablebased calculations. Snapshots at separation intervals of 1.2 Å were used as initial configurations to define a reaction coordinate for 16 windows of umbrella sampling (US) simulations. The first window was the initial bound state of PH domain. The first 11 windows were let to relax for 100 ns, followed by 50 ns of productive simulation; while the last 5 windows were relaxed for 50 ns, followed by 50 ns of productive simulation (the total amount of simulation per system was 2150 ns). As a result, window centers ranged from 58.3 Å to 70.3 Å were used for the free energy calculations. Further extended windows (71.5 Å to 76.3 Å) were added for orientation analysis of binding pockets. To eliminate perturbation to the PIP2 lipid by US window restraints, we have chosen the separation of protein COM and lower-layer lipid backbone phosphorus atoms COM to be restrained as the window centers (with restraint constant k=1000 kJ mol-1 nm-2). No further restraints were added to the system, except the US window restraints. Potential of mean force (PMF) profiles were calculated using the GROMACS g_wham tool58. Convergence was assessed by calculating profiles for 30 ns intervals in 10 ns sliding blocks (Fig. S2). Data before equilibrium was discarded in order to calculate the final profile. Errors were obtained from bootstrap analysis using 200 bootstraps. All profiles have been shifted so that the PMF value in bulk is 0 kcal mol-1.



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

RESULTS MD simulations grant us valuable molecular details of dynamical biological systems. In order to study interactions between ACAP1BAR-PH protein dimer and PIP2-containing multi-component phospholipid membranes, we have performed a number of atomistic MD simulations, including (1) one single ACAP1BAR-PH protein dimer in ionic solvent (2) one single ACAP1BAR-PH protein dimer interacting with PIP2-containing bilayers in ionic solvent. Combining multiple atomistic MD simulations not only improved statistics of analysis by reducing errors due to insufficient sampling, but also uncovered interesting protein dynamics under different environments. Within the availability of computational resources, all molecular systems have been performed multiple times of simulations for ensuring repeatability and improving analytical statistics (simulation details in Methods). Hinted by the molecular details of interactions between ACAP1BAR-PH protein dimer with PIP2-containing bilayers, we also identified two binding pockets on the PH domain of ACAP1BAR-PH protein dimer for PIP2 lipids. We found that although both pockets bound to PIP2 lipids consistently in multiple atomistic MD simulations, the binding pockets could indeed behave very differently. In addition to ordinary unbiased MD simulations, multiple restrained MD simulations have also been performed as umbrella sampling (US) in order to compute the potential of mean force (PMF) profiles. One single PH domain was obtained from the ACAP1BAR-PH protein dimer and placed above either a PIP2- or POPS-containing bilayers. As a result, PMF profiles of a PH domain interacting with a PIP2 or POPS lipid molecule along the normal of the membrane were obtained (simulation details in Methods). We have performed four US simulations in total in which the two binding pockets were bound to either PIP2 or POPS lipid respectively. Combining unbiased and restrained MD simulations, we examined carefully the nature of the two binding pockets.

A flexible loop linking the PH domain and the BAR domain The BAR and PH domains were reported, respectively, as conserved motifs in many different proteins23,33-35. In order to study the cooperative work carried out by these two domains upon the membrane binding events, we sought to examine the relative motion of the PH domain with respect to the BAR domain in ACAP1 proteins. By computing the 8 ACS Paragon Plus Environment

Page 8 of 48

Page 9 of 48

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


three principle axes of the BAR domain, we obtained the orientation of the BAR domain as the three reference axes (yellow arrows in Figure 1). Orientation of the PH domain was defined by one of the principle axes of the sole helix (residue 345-361) in the PH domain pointing along the length of the helix (blue arrows in Figure 1). As an ACAP1BARPH

domain is a dimeric structure of two BAR-PH monomers, the structure show two-fold

symmetry. Therefore, when computing the three reference axes for each PH domain in an ACAP1BAR-PH protein dimer, we would invert the first and the third principle axes of the BAR domain. As a result, the first principle axis of the BAR domain always points in the direction of the blue arrow of the desired PH domain; the third principle axis of the BAR domain always points from the geometry center of the BAR domain towards the desired PH domain; while the second principle axis of the BAR domain remains the same for both PH domains in an ACAP1BAR-PH protein dimer. Lastly, the orientation of PH domain with respect to the BAR domain (α angles) in ACAP1BAR-PH protein is defined as the angles between the orientation of the PH domain (blue arrows in Figure 1) and the three reference axes (yellow arrows in Figure 1) of the BAR domain (α1, α2, and α3). We first computed the α angles of the crystal structure of the ACAP1BAR-PH protein dimer (PDB ID: 4NSW). Though the dimer is composed of two identical monomers, the α angles were found to be different for two PH domains. The results are as follows: α1, 21.92⁰ and 5.58⁰; α2, 68.20⁰ and 85.04⁰; α3, 87.76⁰ and 87.46⁰. Theses angles implied that the PH domain may have flexible orientation angles relative to the BAR domain. The subsequent atomistic MD simulations would start from the crystal structures, therefore these angles also served as the initial α angles for the following dynamical evolution. We have carried out three atomistic MD simulations with a single ACAP1BAR-PH protein dimer in ionic solvent (protein-solvent simulation, details in Methods). Flexibility of the PH domain in ACAP1BAR-PH dimer was examined in the MD simulations. The temporal profiles of the α angles for each PH domain of the dimer for all the three simulations were shown in Figure 2. We note that all the α angles were disperse and indicate high dynamical flexibility of the loop linking the PH domain and the BAR domain. Although α1 angles started differently for the two PH domains in the crystal structure, they seemed to coincide in individual simulation ensembles (Table 1). However, the 9 ACS Paragon Plus Environment


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 48

equilibrated α1 angles were different for three independent simulations (ranging from 11.62⁰ to 22.58⁰), indicating that the dimer permits a variety of α1 angles. Similarly, the α2 angles also started from different angles for the two PH domains, then evolved to angles ranging from 71.82⁰ to 96.00⁰. The α2 angles generally have larger uncertainties, reflected in larger standard deviations (SD), than those of the α1 angles, also indicating flexibility of the loop. Interesting, α3 angles of the two PH domains started from same initial angles, but diverged in some simulations and caused very different residual α3 angles. All the evidence show that the PH domains are loosely linked to the BAR domains. We note that transition between states with different protein structures might require at least tens of microseconds59 MD simulations. We have instead shown that ACAP1BAR-PH protein dimer deviated from identical crystal structures and equilibrated to structures with very different orientation of the PH domain with respect to the BAR domain in multiple MD simulations.

One of the PH domains rotates when binding to the membrane asymmetrically Then, we performed two independent atomistic MD simulations of one single ACAP1BAR-PH protein dimer being placed above PIP2-containing lipid bilayers (proteinlipid simulation, details in Methods). In both simulations, we have observed only one of the PH domains binding to the underneath membrane, i.e. asymmetrical protein binding events. More importantly, with similar starting structures, different ends of the ACAP1BAR-PH protein dimer were observed to bind to the membrane in the simulations, i.e. one end binds in simulation1 and the another end binds in simulation 2. As the PH domains are loosely linked to the BAR domains, the PH domain is expected to explore a wide range of orientation angles. We note that the flexible loop between the BAR and PH domain may provide the chance for one of the PH domain to rotate to an angle facilitating the binding of protein to the lipid bilayers. We characterized the two ends into binding (P-M, P for protein, M for membrane) and nonbinding (P-W, W for water) ends. Figure 3 shows that, upon binding of the protein to the membrane, both α1 and α2 angles converged to a much smaller range of angles compared to the case of protein in water. The α3 angles, in contrast, show very different 10 ACS Paragon Plus Environment

Page 11 of 48

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


equilibrated angles (in average of two PH domains, 93.24⁰ for P-M and 86.81⁰ for P-W), clearly distinguished the binding and non-binding state of the PH domains (Table 1). Moreover, the P-M angles were seen to have smaller SD than those of the P-W angles in both protein-lipid and protein-solvent simulations. The smaller SD is intuitively caused by the interactions with the underneath membranes. As a result, the bound PH domain show consistent orientations with small SD, while the unbound PH domain was associated with a different orientation with relatively larger SD of α3 angles. To conclude, introducing PIP2-containing lipid bilayers under ACAP1BAR-PH protein dimer will cause asymmetrical binding events, which associate with protein structural changes reflected by rotation of one of the PH domains with respect to the BAR domain.

Orientation analysis revealed the process of protein-membrane binding To understand protein binding to the membrane process from its configuration in the aqueous environment, we summarized the α angle distribution of all conducted MD simulations in Figure 4. We suggest that internal structural changes might have occurred when ACAP1BAR-PH protein dimers interact with the lipids (Figure 5). First, when an ACAP1BAR-PH protein dimer stands alone under physiological environment, both of the PH domains is at the unbound (P-W) state (Figure 5A). Without neighboring lipids, the PH domain of ACAP1BAR-PH protein dimer was found to access a wide range of α angle (Figure 4A-C) with relatively large SD (up to 8⁰ in Table 1). When placing PIP2containing membranes into the simulation system, ACAP1BAR-PH protein dimer binds to the lipids asymmetrically and develops P-M and P-W states for the two PH domains (Figure 5B). Asymmetrical binding events have been observed in multiple independent protein-lipid simulations. It has been shown that both PH domains in an ACAP1BAR-PH dimer are capable for binding to PIP2-containing lipid bilayers. Comparing the α angle to protein-solvent systems, protein-lipid systems certainly have more concentrated angle distribution due to the interactions between the protein and lipids (Figure 4D-F). The presence of lipids established bound states of the PH domain, therefore confined the dynamics of the bound protein dimer. In general, the fluctuations of the α angle in the bounded state decrease to 3-5⁰ from 3-9⁰ in unbounded states (Table 1). Especially for



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 48

the α3 angle (Figure 4F), a deviation in P-M and P-W angles indicates strong interactions between the protein and nearby lipids, or even specific kind of lipids. By combining multiple all-atom MD simulations, orientation analysis has shown some molecular details of the membrane binding process of ACAP1BAR-PH protein dimers. We note that as ACAP1 protein was previously reported to tubulate lipid liposomes. Well-structured lattices of ACAP1 protein upon binding to cell membrane were observed by cryo-electron microscopy (cryo-EM)32. Due to asymmetrical binding of ACAP1BAR-PH protein dimer, formation of protein lattice would intrinsically introduce large area of curvatures, causing tubulation upon establishment of protein-lipid and proteinprotein interactions. Recent studies on other BAR proteins60,61 has suggested factors, such as protein density, behind the complicated oligomerization events. For ACAP1 proteins, our orientation analysis might be able to provide valuable molecular information for the tubulation process. While this study focuses on study of protein-lipid interactions, details of protein-protein interactions require further studies in the future.

MD simulations reveal two binding pockets of PH domains Analysis of orientation of the PH domains has pushed further our studies on the protein-lipid binding state of ACAP1BAR-PH protein dimer. In agreement with the previously reported PIP2 specificity of ACAP1 protein32, we have observed binding of PIP2 to PH domains in two protein-lipid simulations (Figure 6A). PIP2 binding events were facilitated by two binding pockets in the PH domain of ACAP1BAR-PH protein dimer. In both simulations, we observed at least one PIP2 bound to each binding pocket. Two binding pockets are defined according to amino acid loops locate near where the PIP2 lipid bound. Three key residue loops were thus noted and colored in Figure 6A. Loop 1 (residue 277-283, β1/β2 loop), reported earlier by Pang et al.32, inserts into the lipid bilayer and contributes most to the protein-lipid binding event. Loop 3 (residue 321-327, β5/β6 loop), together with Loop 1, composes the binding pocket 1 (Pocket 1) and interacts with a PIP2 lipid during the early stage of the simulations. Loop 2 (residue 300305, β3/β4 loop), together with Loop 1, composes the binding pocket 2 (Pocket 2) and interacts with another PIP2 lipid after a period of time. Trajectories of the bound PIP2 12 ACS Paragon Plus Environment

Page 13 of 48

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


were shown as the separation between the lipid backbone phosphate and the binding pockets (Figure 6B-C). Upon the binding of ACAP1BAR-PH protein dimer to lipid bilayer underneath, Pocket 1 was found to establish stable interactions between a nearby PIP2 lipid since a sharp decrease in pocket-lipid separation at ~100 ns of simulation 1 and ~20 ns of simulation 2 until the end of the simulation; while Pocket 2 gradually established such interactions throughout the simulation. Interestingly, although simulation 1 and 2 were independent simulations, the development of pocket-lipid interactions was very similar. Our orientation analysis suggested that the bound PH domain might undergo gradual rotation after inserting Loop 1 into lipid bilayers at ~100ns (Figure 2A-C). As a result, rotation of PH domain gradually decreased separation between the pocket and bilayer surface, facilitating binding of PIP2 lipid to Pocket 2. In order to identify key residues for the two pockets, Figure 6D show contacts between residues of the PH domain with nearby PIP2 lipids. Mapping of Figure 6D to molecular surface of a PH domain confirmed that the two binding pockets locate at two sides of the inserting Loop 1 (Figure 6E). Loop 1, 2 and 3 all have a number of contacts with PIP2 molecules. Specifically, residues D322, S323, E324, and R325 of loop 3 have the highest amount of contacts among the loops and contribute mainly to Pocket 1. Residues N281 and T282 of loop 1 and residues K299 and K300 of loop 2 contribute mostly to the formation of Pocket 2. These results confirm the significance of electrostatic interactions between the PH domain of ACAP1 protein and PIP2 lipids. Identification of key residues for the protein-lipid interactions not only hints for the binding interfaces, but also provides valuable information for future mutation experiments.

Nature of the two binding pockets of PH domains We examined the binding interface between the PH domains and membrane and found the PH domain rotates relative to the membrane normal direction during the PIP2 binding processes. Definition of the binding interface of the two binding pockets involves coarse-graining of the identified key residue loops. Loop 1 was coarse-grained into one bead according to the COM of the residues (Figure 7A-B). Loop 2 (in Figure 7A) and 13 ACS Paragon Plus Environment


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 48

Loop 3 (in Figure 7B) were split into halves in order to define two beads for each loop. As a result, we were able to define a binding interface by connecting the three beads for each binding pocket (the orange triangle in Figure 7A and the blue triangle in Figure 7B). Dynamics of the binding interface could then be characterized by calculating angles between the normal of the interface and the normal of the membrane surface (Figure 7C-F). In both simulations, Pocket 1 was seen to be heavily affected by the underneath membrane, reflected by the P-M angle having significantly smaller SD (3-5º) than the PW angle (9-11º) (Figure 7C, E). In contrast, Pocket 2 underwent steady changes for both P-M and P-W angles, having relatively smaller SD (both 3-6º) (Figure 7D, F). This result supports our previous observation that PIP2 lipids follow a different but consistent trajectory heading the two binding pockets (Figure 6B-C). Interestingly, we observed that the orientation of Pocket 1 changed from ~50º to 60º (Figure 7C) and that of Pocket 2 changed gradually from ~80º to 60º (Figure 7D). This observation coincided with the fact that PIP2 lipids took time to approach Pocket 2 (Figure 6B) and the PH domain rotated to complete the binding process. Requirement in orientation of the PH domain has also been reported in various studies10,12 in order to facilitate binding of lipids to different binding sites. Previous studies of GRP1 PH domain11,12 have revealed the same two binding pockets, in which a PH domain was placed alone and allowed to rotate freely, residues correspond to Pocket 2 were reported to have the highest contact frequencies with bound PIP2 ligand, namely the “canonical site”. In order to find out whether the same difference preserves in PH domain of the ACAP1BAR-PH protein dimer, we conducted multiple umbrella sampling simulations (Figure 8A) and generated PMF profiles of two pockets binding to PIP2 lipids respectively (the blue and red profile in Figure 8B). Both profiles have a global minimum at a protein-lipid separation of ~58 Å, with a well depth of -4.5 kcal mol-1 for Pocket 1 and -6.5 kcal mol-1 for Pocket 2. These profiles suggest favorable binding of both pockets to PIP2 with a preference for Pocket 2. Identification of two pockets and associated PIP2 preference is in agreement with both previous experimental and computational results9-12,62,63. In order to account for the specific binding to PIP2 lipids by PH domain, we have performed two more US simulations by replacing the PIP2 lipid at two binding pockets 14 ACS Paragon Plus Environment

Page 15 of 48

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


with POPS lipid respectively. PMF profiles of two pockets binding to POPS lipids are also shown as yellow and purple plots in Figure 8B. The PMF of both pockets binding to POPS have a well depth of -2.0 kcal mol-1 and relative shorter range of coverage (about 6 Å), we demonstrate their preference for PIP2 over POPS lipid. Different well depths of the PMF profiles could be explained partly by the charged residues located inside the pockets. Pocket 1 consists of three positively charged residues: R275, K281, and R325; while Pocket 2 consists four positively charged residues: K274, R286, K299, and K336. Pocket 1 consists of one negatively charged residue E324, while Pocket 2 does not have any. Therefore, the less negative POPS lipids (-1e) are expected to be electrostatically attracted by the charged residues less firmly, compared to PIP2 lipids (3e). As a result, Pocket 2/PIP2 possessed the lowest well depth of PMF profile, Pocket 1/PIP2 came the second. Pocket 1/POPS and Pocket 2/POPS was seen to have similar and the shallowest PMF well. In addition to the PMF profile, our US simulations also provide further insights into nature of the binding pockets. Angle distributions of orientation of the binding interfaces during US simulations were given in Figure 9, in which the restrained separation covered a distance of ~15 Å from the bound state. With a decrease in separations, Pocket 1 first spread over a wide range of orientation angles, then quickly converged to ~60° when the PIP2 lipids came near the pocket (Figure 9A). Pocket 2 proved itself as a long-range interacting binding pocket, compared to Pocket 1, as the pocket orientation evolved gradually from ~80° to ~40° with decreasing restrained separations (Figure 9B). Difference in the effective range of the pockets also reflects in the PMF profiles, in which the Pocket 2/PIP2 well lifted over a longer distance of ~2 Å than the Pocket 1/PIP2. We note that although both our energetic analysis and previous studies11,12 confirm that Pocket 2 is the preferred binding pocket for PIP2 lipids, Pocket 2 was found to bind PIP2 lipids later than Pocket 1 in both of our protein-lipid simulations (Figure 6BC). It implies that multiple factors could lie behind protein-membrane binding events. As seen in Figure 6A, existence of the BAR domain in ACAP1BAR-PH protein dimer might assist the binding, therefore confine the binding angle of the PH domain to cell membranes. As reported by Cui et al.31, residues at the tip of the BAR domain might 15 ACS Paragon Plus Environment


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 48

also involve in membrane binding. Preference binding by Pocket 2 then took effects by gradually rotating the PH domain and establishing interactions with a relatively far away PIP2 lipid. In this way, the ACAP1BAR-PH protein dimers could quickly target and bind to desire regions of the cell membrane using Pocket 1, then reinforce such binding by strongly interacting with more PIP2 lipids using Pocket 2.


Page 17 of 48

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


CONCLUSION Our results provide insight into the structural dynamics of ACAP1BAR-PH protein dimer when binding to cell membranes. Our simulations reveal that the PH domains of an ACAP1BAR-PH protein dimer would undergo orientation changes during the membrane binding process. Such process was difficult to be caught in conventional experiments yet has suggested important residues for protein-lipid interactions. The list of dominating charged residues serve for the purpose of understanding and manipulating the membrane binding process by the PH domain of ACAP1BAR-PH protein or other similar peripheral protein domains. A full understanding of membrane tubulation process will have to await explicit simulations of multiple ACAP1BAR-PH protein dimers on PIP2containing membranes. However, the present simulations already reveal the principle protein structural deformation and key residue information upon protein-lipid interactions. In addition, protein-lipid interactions between ACAP1BAR-PH protein dimer and PIP2-containing membranes have been studied by revealing two binding pockets on the PH domain. Preference for the two binding pockets is explained by PMF profiles generated by US simulations over separation of 12 Å between protein and lipids. Nature of the two binding pockets has been discussed in details with evolution of the binding process and binding interface orientations. More importantly, preference of the PH domain binding to PIP2 lipids over POPS lipids was also shown in the PMF results. Exploring the free energy surface of the protein-lipid systems revealed a much deeper well when a PIP2 lipid was in either binding pocket of the protein, compared to the case of POPS lipid.



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

Supporting Information Figure S1 displays the Root-Mean Square Deviation (RMSD) profiles of the five MD simulations (namely sim1-5) in order to show the convergence of the simulations. Convergence of the US simulations was shown in Figure S2 by calculating profiles of 30 ns intervals in 10 ns sliding blocks. Histograms of the last 50 ns interval were displayed to the right of the PMF profiles.


Page 18 of 48

Page 19 of 48

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


AUTHOR INFORMATION Corresponding author *Jun Fan Department of Physics and Materials Science, City University of Hong Kong Tat Chee Avenue, Kowloon Tong, Kowloon, Hong Kong, China Tel: +(852) 3442 9978 Email: [email protected]



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 48

ACKNOWLEDGEMENTS The authors thanks the projects supported by the National Natural Science Foundation of China (21403182), the Research Grants Council of Hong Kong (CityU 21300014), and CityU grants (7004387 & 9680136). The authors would also like to acknowledge grants to F.S. by the Chinese Ministry of Science and Technology (2014CB910700) and the Strategic Priority Research Program of Chinese Academy of Sciences (XDB08030202). This research has used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC0500OR22725. This research has also used computing resources supported by Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase). The authors have no conflicts of interest to dis-close.


Page 21 of 48

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


REFERENCES (1) (2) (3) (4) (5) (6) (7)

(8)

(9)

(10)

(11)

(12)

(13) (14)

(15)

(16)

(17)

Kuriyan, J.; Cowburn, D. Modular Peptide Recognition Domains in Eukaryotic Signaling. Annu Rev Biophys Biomol Struct 1997, 26, 259–288. Pawson, T.; Scott, J. D. Signaling Through Scaffold, Anchoring, and Adaptor Proteins. Science 1997, 278, 2075–2080. Hurley, J. H.; Misra, S. Signaling and Subcellular Targeting by MembraneBinding Domains. Annu Rev Biophys Biomol Struct 2000, 29, 49–79. Lemmon, M. A. Membrane Recognition by Phospholipid-Binding Domains. Nat. Rev. Mol. Cell Biol. 2008, 9, 99–111. Kutateladze, T. G. Translation of the Phosphoinositide Code by PI Effectors. Nature Chemical Biology 2010, 6, 507–513. Stahelin, R. V.; Scott, J. L.; Frick, C. T. Cellular and Molecular Interactions of Phosphoinositides and Peripheral Proteins. Chem. Phys. Lipids 2014, 182, 3–18. Perilla, J. R.; Goh, B. C.; Cassidy, C. K.; Liu, B.; Bernardi, R. C.; Rudack, T.; Yu, H.; Wu, Z.; Schulten, K. Molecular Dynamics Simulations of Large Macromolecular Complexes. Current Opinion in Structural Biology 2015, 31, 64– 74. Simunovic, M.; Voth, G. A.; Callan-Jones, A.; Bassereau, P. When Physics Takes Over: BAR Proteins and Membrane Curvature. Trends Cell Biol. 2015, 1– 13. Lumb, C. N.; He, J.; Xue, Y.; Stansfeld, P. J.; Stahelin, R. V.; Kutateladze, T. G.; Sansom, M. S. P. Biophysical and Computational Studies of Membrane Penetration by the GRP1 Pleckstrin Homology Domain. Structure 2011, 19, 1338–1346. Lai, C.-L.; Srivastava, A.; Pilling, C.; Chase, A. R.; Falke, J. J.; Voth, G. A. Molecular Mechanism of Membrane Binding of the GRP1 PH Domain. Journal of Molecular Biology 2013, 425, 3073–3090. Buyan, A.; Kalli, A. C.; Sansom, M. S. P. Multiscale Simulations Suggest a Mechanism for the Association of the Dok7 PH Domain with PIP-Containing Membranes. PLoS Comput Biol 2016, 12, e1005028–14. Naughton, F. B.; Kalli, A. C.; Sansom, M. S. P. Association of Peripheral Membrane Proteins with Membranes: Free Energy of Binding of GRP1 PH Domain with Phosphatidylinositol Phosphate-Containing Model Bilayers. The Journal of Physical Chemistry Letters 2016, 7, 1219–1224. Roux, B. The Calculation of the Potential of Mean Force Using ComputerSimulations. 1995, 91, 275–282. Christ, C. D.; Mark, A. E.; van Gunsteren, W. F. Basic Ingredients of Free Energy Calculations: a Review. Journal of Computational Chemistry 2009, NA– NA. Deng, Y.; Roux, B. Computations of Standard Binding Free Energies with Molecular Dynamics Simulations. The Journal of Physical Chemistry B 2009, 113, 2234–2246. Zeller, F.; Zacharias, M. Efficient Calculation of Relative Binding Free Energies by Umbrella Sampling Perturbation. Journal of Computational Chemistry 2014, 35, 2256–2262. Grasso, G.; Deriu, M. A.; Prat, M.; Rimondini, L.; Vernè, E.; Follenzi, A.; Danani, 21 ACS Paragon Plus Environment


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

(18)

(19)

(20)

(21)

(22)

(23)

(24)

(25)

(26) (27) (28)

(29) (30) (31) (32)

(33)

A. Cell Penetrating Peptide Adsorption on Magnetite and Silica Surfaces: a Computational Investigation. The Journal of Physical Chemistry B 2015, 119, 8239–8246. Sancho, M. I.; Andujar, S.; Porasso, R. D.; Enriz, R. D. Theoretical and Experimental Study of Inclusion Complexes of Β-Cyclodextrins with Chalcone and 2“,4-”Dihydroxychalcone. The Journal of Physical Chemistry B 2016, 120, 3000–3011. Jackson, T. R.; Brown, F. D.; Nie, Z.; Miura, K.; Foroni, L.; Sun, J.; Hsu, V. W.; Donaldson, J. G.; Randazzo, P. A. ACAPs Are Arf6 GTPase-Activating Proteins That Function in the Cell Periphery. J Cell Biol 2000, 151, 627–638. Dai, J.; Li, J.; Bos, E.; Porcionatto, M.; Premont, R. T.; Bourgoin, S.; Peters, P. J.; Hsu, V. W. ACAP1 Promotes Endocytic Recycling by Recognizing Recycling Sorting Signals. Developmental Cell 2004, 7, 771–776. Li, J.; Peters, P. J.; Bai, M.; Bai, M.; Dai, J.; Bos, E.; Kirchhausen, T.; Kandror, K. V.; Hsu, V. W. An ACAP1-Containing Clathrin Coat Complex for Endocytic Recycling. J Cell Biol 2007, 178, 453–464. Bai, M.; Pang, X.; Lou, J.; Zhou, Q.; Zhang, K.; Ma, J.; Li, J.; Sun, F.; Hsu, V. W. Mechanistic Insights Into Regulated Cargo Binding by ACAP1 Protein. J Biol Chem 2012, 287, 28675–28685. Peter, B. J.; Kent, H. M.; Mills, I. G.; Vallis, Y.; Butler, P. J. G.; Evans, P. R.; McMahon, H. T. BAR Domains as Sensors of Membrane Curvature: the Amphiphysin BAR Structure. Science 2004, 303, 495–499. Blood, P. D.; Voth, G. A. Direct Observation of Bin/Amphiphysin/Rvs (BAR) Domain-Induced Membrane Curvature by Means of Molecular Dynamics Simulations. PNAS 2006, 103, 15068–15072. Ayton, G. S.; Lyman, E.; Krishna, V.; Mim, C.; Unger, V. M.; Voth, G. A. New Insights Into BAR Domain-Induced Membrane Remodeling. Biophysical Journal 2009, 97, 1616–1625. Yu, H.; Schulten, K. Membrane Sculpting by F-BAR Domains Studied by Molecular Dynamics Simulations. PLoS Comput Biol 2013, 9, e1002892. Arkhipov, A.; Yin, Y.; Schulten, K. Four-Scale Description of Membrane Sculpting by BAR Domains. Biophysical Journal 2008, 95, 2806–2821. Simunovic, M.; Srivastava, A.; Voth, G. A. Linear Aggregation of Proteins on the Membrane as a Prelude to Membrane Remodeling. PNAS 2013, 110, 20396– 20401. Arkhipov, A.; Yin, Y.; Schulten, K. Membrane-Bending Mechanism of Amphiphysin N-BAR Domains. Biophysical Journal 2009, 97, 2727–2735. Yin, Y.; Arkhipov, A.; Schulten, K. Simulations of Membrane Tubulation by Lattices of Amphiphysin N-BAR Domains. Structure 2009, 17, 882–892. Cui, H.; Ayton, G. S.; Voth, G. A. Membrane Binding by the Endophilin N-BAR Domain. Biophysical Journal 2009, 97, 2746–2753. Pang, X.; Fan, J.; Zhang, Y.; Zhang, K.; Gao, B.; Ma, J.; Li, J.; Deng, Y.; Zhou, Q.; Egelman, E. H.; et al. A PH Domain in ACAP1 Possesses Key Features of the BAR Domain in Promoting Membrane Curvature. Developmental Cell 2014, 31, 73–86. Haslam, R. J.; Koide, H. B.; Hemmings, B. A. Pleckstrin Domain Homology. 22 ACS Paragon Plus Environment

Page 22 of 48

Page 23 of 48

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


(34) (35)

(36) (37) (38) (39) (40)

(41) (42) (43) (44)

(45)

(46)

(47)

(48)

(49)

(50)

Nature 1993, 363, 309–310. Mayer, B. J.; Ren, R.; Clark, K. L.; Baltimore, D. A Putative Modular Domain Present in Diverse Signaling Proteins. Cell 1993, 73, 629–630. Musacchio, A.; Gibson, T.; Rice, P.; Thompson, J.; Saraste, M. The PH Domain: a Common Piece in the Structural Patchwork of Signalling Proteins. Trends Biochem. Sci. 1993, 18, 343–348. Bottomley, M. J.; Salim, K.; Panayotou, G. Phospholipid-Binding Protein Domains. Biochim. Biophys. Acta 1998, 1436, 165–183. Lemmon, M. A.; Ferguson, K. M. Pleckstrin Homology Domains. Curr. Top. Microbiol. Immunol. 1998, 228, 39–74. Di Paolo, G.; De Camilli, P. Phosphoinositides in Cell Regulation and Membrane Dynamics. Nature 2006, 443, 651–657. Wymann, M. P.; Schneiter, R. Lipid Signalling in Disease. Nat. Rev. Mol. Cell Biol. 2008, 9, 162–176. Posor, Y.; Eichhorn-Gruenig, M.; Puchkov, D.; Schöneberg, J.; Ullrich, A.; Lampe, A.; Müller, R.; Zarbakhsh, S.; Gulluni, F.; Hirsch, E.; et al. Spatiotemporal Control of Endocytosis by Phosphatidylinositol-3,4-Bisphosphate. Nature 2013, 499, 233–237. van Meer, G.; Voelker, D. R.; Feigenson, G. W. Membrane Lipids: Where They Are and How They Behave. Nat. Rev. Mol. Cell Biol. 2008, 9, 112–124. Rebecchi, M. J.; Scarlata, S. Pleckstrin Homology Domains: a Common Fold with Diverse Functions. Annu Rev Biophys Biomol Struct 1998, 27, 503–528. Lemmon, M.; Ferguson, K. Signal-Dependent Membrane Targeting by Pleckstrin Homology (PH) Domains. Biochem. J. 2000, 1, 9–9. Miaczynska, M.; Christoforidis, S.; Giner, A.; Shevchenko, A.; UttenweilerJoseph, S.; Habermann, B.; Wilm, M.; Parton, R. G.; Zerial, M. APPL Proteins Link Rab5 to Nuclear Signal Transduction via an Endosomal Compartment. Cell 2004, 116, 445–456. Zhu, G.; Chen, J.; Liu, J.; Brunzelle, J. S.; Huang, B.; Wakeham, N.; Terzyan, S.; Li, X.; Rao, Z.; Li, G.; et al. Structure of the APPL1 BAR-PH Domain and Characterization of Its Interaction with Rab5. The EMBO Journal 2007, 26, 3484–3493. Lumb, C. N.; Sansom, M. S. P. Finding a Needle in a Haystack: the Role of Electrostatics in Target Lipid Recognition by PH Domains. PLoS Comput Biol 2012, 8, e1002617. Jo, S.; Kim, T.; Iyer, V. G.; Im, W. CHARMM-GUI: a Web-Based Graphical User Interface for CHARMM. Journal of Computational Chemistry. August 2008, pp 1859–1865. Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. The Journal of Chemical Physics 1983, 79, 926–11. Humphrey, W.; Humphrey, W.; Dalke, A.; Dalke, A.; Schulten, K.; Schulten, K. VMD: Visual Molecular Dynamics. Journal of Molecular Graphics and Modelling 1996, 14, 33–38. Phillips, J. C.; Braun, R.; Wang, W.; Gumbart, J.; Tajkhorshid, E.; Villa, E.; Chipot, C.; Skeel, R. D.; Kalé, L.; Schulten, K. Scalable Molecular Dynamics 23 ACS Paragon Plus Environment


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

(51)

(52)

(53)

(54)

(55)

(56)

(57)

(58)

(59)

(60)

(61)

(62)

Page 24 of 48

with NAMD. Journal of Computational Chemistry 2005, 26, 1781–1802. MacKerell, A. D.; Bashford, D.; Bellott, M.; Dunbrack, R. L.; Evanseck, J. D.; Field, M. J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. The Journal of Physical Chemistry B 1998, 102, 3586–3616. Mackerell, A. D., Jr.; Feig, M.; Brooks, C. L., III. Extending the Treatment of Backbone Energetics in Protein Force Fields: Limitations of Gas-Phase Quantum Mechanics in Reproducing Protein Conformational Distributions in Molecular Dynamics Simulations. Journal of Computational Chemistry 2004, 25, 1400–1415. Best, R. B.; Zhu, X.; Shim, J.; Lopes, P. E. M.; Mittal, J.; Feig, M.; Mackerell, A. D., Jr. Optimization of the Additive CHARMM All-Atom Protein Force Field Targeting Improved Sampling of the Backbone Φ, Ψ and Side-Chain Χ 1and Χ 2Dihedral Angles. Journal of Chemical Theory and Computation 2012, 8, 3257– 3273. Klauda, J. B.; Venable, R. M.; Freites, J. A.; O'Connor, J. W.; Tobias, D. J.; Mondragon-Ramirez, C.; Vorobyov, I.; Mackerell, A. D.; Pastor, R. W. Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid Types. The Journal of Physical Chemistry B 2010, 114, 7830–7843. Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: an N⋅ Log (N) Method for Ewald Sums in Large Systems. The Journal of Chemical Physics 1993, 98, 10089. Ryckaert, J.-P.; Ciccotti, G.; Berendsen, H. J. C. Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of N-Alkanes. Journal of Computational Physics 1977, 23, 327–341. Feller, S. E.; Zhang, Y.; Pastor, R. W.; Brooks, B. R. Constant Pressure Molecular Dynamics Simulation: the Langevin Piston Method. The Journal of Chemical Physics 1995, 103, 4613–4621. Hub, J. S.; de Groot, B. L.; Van Der Spoel, D. G_Wham—a Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. Journal of Chemical Theory and Computation 2010, 6, 3713–3720. Shaw, D. E.; Maragakis, P.; Lindorff-Larsen, K.; Piana, S.; Dror, R. O.; Eastwood, M. P.; Bank, J. A.; Jumper, J. M.; Salmon, J. K.; Shan, Y.; et al. Atomic-Level Characterization of the Structural Dynamics of Proteins. Science 2010, 330, 341–346. Kelley, C. F.; Messelaar, E. M.; Eskin, T. L.; Wang, S.; Song, K.; Vishnia, K.; Becalska, A. N.; Shupliakov, O.; Hagan, M. F.; Danino, D.; et al. Membrane Charge Directs the Outcome of F-BAR Domain Lipid Binding and Autoregulation. CellReports 2015, 13, 2597–2609. Isas, J. M.; Ambroso, M. R.; Hegde, P. B.; Langen, J.; Langen, R. Tubulation by Amphiphysin Requires Concentration-Dependent Switching From Wedging to Scaffolding. Structure 2015, 23, 873–881. Chen, H.-C.; Ziemba, B. P.; Landgraf, K. E.; Corbin, J. A.; Falke, J. J. Membrane Docking Geometry of GRP1 PH Domain Bound to a Target Lipid Bilayer: an EPR Site-Directed Spin-Labeling and Relaxation Study. PLoS ONE 2012, 7, e33640. 24 ACS Paragon Plus Environment

Page 25 of 48

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


(63)

He, J.; Haney, R. M.; Vora, M.; Verkhusha, V. V.; Stahelin, R. V.; Kutateladze, T. G. Molecular Mechanism of Membrane Targeting by the GRP1 PH Domain. J. Lipid Res. 2008, 49, 1807–1815.



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 26 of 48

Table 1. Orientation analysis of the PH domain in ACAP1BAR-PH protein dimer. Simulation 1-3 were the protein-solvent simulations. Simulations 4-5 were the proteinmembrane simulations. α1, α2, and α3 are the angles between the orientation of the PH domain (blue arrows in Figure 1) and the three reference axes (yellow arrows in Figure 1) of the BAR domain. P-W refers to the PH domain interacting with solvent. P-M refers to the PH domain interacting with the membrane. Error estimation from standard deviation of the data is included in brackets. Simulation 1 Simulation 2 Simulation 3 Simulation 4 Simulation 5

11.62 (4.32)

8.24 (4.43)

21.50 (4.81)

P-W

α1

12.07 (6.69)

8.11 (4.29)

89.27 (5.87)

83.08 (8.26)

87.02 (6.44)

92.03 (6.93)

96.18 (5.81)

93.19 (4.82)

86.85 (9.06)

84.28 (5.35)

81.89 (5.78)

84.28 (4.70)

87.22 (7.59)

86.39 (3.38)

92.63 (4.09)

93.85 (3.44)

70.35 (4.57)

P-W

α3

8.23 (3.76)

71.82 (5.35)

P-M 80.64 (4.53)

10.23 (4.30)

96.00 (5.99)

P-W

α2

8.55 (3.59)

22.58 (5.59)

P-M 90.32 (6.66)

10.92 (6.36)

102.37 (3.79)

P-M


Page 27 of 48

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


FIGURE LEGENDS Figure 1. Definition of the orientation of PH domain with respect to the BAR domain. An ACAP1BAR-PH protein dimer is shown with the BAR domain colored in blue and the PH domain in magenta. Three principle axes of the BAR domain are labeled shown as yellow arrows. Orientation of the PH domain corresponds to the principle axis of residue 345-361 (α5 helix) along the length of the helix and is shown as the blue arrow. The snapshot was rendered with VMD49. Figure 2. Temporal profile of the orientation analysis of three protein-solvent simulations, namely simulation 1 (sim 1), 2 (sim 2), and 3 (sim 3). (A-C) α1, α2, and α3, angles. Profiles of the two PH domains of the ACAP1BAR-PH protein dimer are shown separately (PH1 and PH2). Running averages were computed to show more clearly the evolution of the orientation angles. Figure 3. Temporal profile of the orientation analysis of two protein-lipid simulations, namely simulation 4 (sim 4) and 5 (sim 5). (A-C) α1, α2, and α3, angles. Profiles of the two PH domains of the ACAP1BAR-PH protein dimer are shown separately as P-M (blue) and P-W (green). As while one PH domain bound to the lipid bilayers, the another one remained unbound. Running averages were computed to show more clearly the evolution of the orientation angles. Figure 4. Angle distribution of all the equilibrated data from Figure 2-3. (A-C) Proteinsolvent simulations. (D-F) Protein-lipid simulations. Angles at the P-W state are shown in green, and P-M in blue. Bin sizes of 2° were used. Figure 5. A simplified model of an ACAP1BAR-PH protein dimer during membrane binding process. Snapshot of an ACAP1BAR-PH dimer from (A) MD simulation of an ACAP1BAR-PH dimer in ionic solvent (B) MD simulation of an ACAP1BAR-PH dimer on lipid bilayers. (A-B) Two monomers of the protein dimer are colored in a sense that PH domain of the green monomer interacting mostly with water and ion molecules (P-W state) and the blue one interacting with the membrane (P-M state). Similar to Figure 1, the blue arrows define the orientation of the PH domains; the three yellow arrows in (B), which were calculated



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 28 of 48

from the principle axes of the BAR domains, define the orientation of the BAR domains. All snapshots were rendered with VMD49. Figure 6. (A) Snapshot of a PH domain of the ACAP1BAR-PH protein dimer from proteinlipid simulation 4. The PH domain is shown with cartoon representations in silver, except for the key residue loops. Three key residue loops defining the binding pocket interface were colored in lighter blue (Loop 1, residue 277-283), green (Loop 2, 300305), and red (Loop 3, residue 321-327). Two bound PIP2 lipid head molecules are colored in orange (Pocket 1) and cyan (Pocket 2) to distinguish the binding pockets. Residues in contact with the PIP2 lipids are colored according to their side chain properties (magenta for negative, blue for positive, green for polar and white for nonpolar). The BAR domain is shown with cartoon representations in transparent white. (BC) Temporal profiles of separation between the pocket interface COM and the bound PIP2 lipid backbone phosphorus atoms. COM of Loop 1 and Loop 3 defines the Pocket 1 interface COM, and that of Loop 1 and Loop 2 defines the Pocket 2. (D) Number of contacts were counted every 50 ps for protein-lipid simulation 4 and then averaged over the number of frames. A contact is defined with a 12 Å cutoff between COM of the residue and the lipid backbone phosphorus atom. (E) The contact distribution in (D) projected onto the PH domain molecular surface, with blue corresponds to no contact and red corresponds to the most contacts. Loop 0, which locates between the two binding pockets, is circled. All snapshots were rendered with VMD49. Figure 7. (A) Snapshot of Pocket 1 from protein-lipid simulation 4 at 255 ns. The blue sphere corresponds to the COM of Loop 1. The red sphere corresponds to the COM of two halves of Loop 3. The binding interface is represented as orange triangle with orange arrow indicating the normal. (B) Snapshot of Pocket 2 from protein-lipid simulation 4 at 418 ns. The blue sphere corresponds to the COM of Loop 1. The green sphere corresponds to the COM of two halves of Loop 2. The binding interface is represented as cyan triangle with cyan arrow indicating the normal. (A-B) Loop 1, 2, and 3 are colored in blue, green, and red respectively. The rest residues are colored in white. Underneath lipid backbone phosphorus atoms are represented as transparent spheres. The black arrow indicates the normal of the membrane. (C-D) Orientation of Pocket 1 28 ACS Paragon Plus Environment

Page 29 of 48

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


and Pocket 2 in protein-lipid simulation 4. (E-F) Orientation of Pocket 1 and Pocket 2 in protein-lipid simulation 5. (C-F) Running averages were computed to show more clearly the evolution of the angles. Color convention refers to the previous orientation analysis of PH domains in Figure 3. Figure 8. (A) Snapshot of a PH domain at bound state with PIP2 lipids, from which we prepare all the initial structures for the US simulations. The atom coordinates were obtained from the last frame of protein-lipid simulation 4. The PH domain is shown in gray in cartoon representation; DOPC backbone phosphorus atoms of the lipid bilayers are represented by pink spheres. Two PIP2 lipids show in licorice representation indicate the position of PIP2 lipids in two independent US setups. The orange PIP2 lipid represents the initial bound state of the Pocket 1/PIP2 setup; the cyan PIP2 lipid represents that of the Pocket 2/PIP2 setup. The black arrow indicates the direction of pulling during SMD simulations for preparing the subsequent US windows. The black dotted line indicates the separation used for restraints of US windows, which is the separation between protein COM and lower-layer lipid backbone phosphorus atoms COM. (B) PMF profiles for PH domain of ACAP1BAR-PH protein dimer bound to PIP2 at Pocket 1 and Pocket 2. Protein-lipid separation was measured from the protein COM to the COM of lower layer lipid backbone phosphorus atoms along the membrane normal. Error estimates were obtained from bootstrap analysis. Figure 9. Three-dimensional histograms of orientation angle of (A) the Pocket 1/PIP2 US (B) the Pocket 2/PIP2 US throughout the 15 US windows. The angle refers to the orientation of the two binding interfaces. The restrained separations are the window centers of US simulations which measured the separation between protein COM and lower-layer lipid backbone phosphorus atoms COM. The bin size is 1⁰.



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

Figure 1.


Page 30 of 48

Page 31 of 48

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


Figure 2.



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

Figure 3.


Page 32 of 48

Page 33 of 48

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


Figure 4.



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

Figure 5.


Page 34 of 48

Page 35 of 48

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


Figure 6.



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

Figure 7.


Page 36 of 48

Page 37 of 48

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


Figure 8.



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

Figure 9.


Page 38 of 48

Page 39 of 48

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


TOC Figure



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

Figure 1. Definition of the orientation of PH domain with respect to the BAR domain. An ACAP1BAR-PH protein dimer is shown with the BAR domain colored in blue and the PH domain in magenta. Three principle axes of the BAR domain are labeled shown as yellow arrows. Orientation of the PH domain corresponds to the principle axis of residue 345-361 (5 helix) along the length of the helix and is shown as the blue arrow. 134x74mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 40 of 48

Page 41 of 48

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


Figure 2. Temporal profile of the orientation analysis of three protein-solvent simulations, namely simulation 1 (sim 1), 2 (sim 2), and 3 (sim 3). (A-C) 1, 2, and 3, angles. Profiles of the two PH domains of the ACAP1BAR-PH protein dimer are shown separately (PH1 and PH2). Running averages were computed to show more clearly the evolution of the orientation angles.

209x78mm (300 x 300 DPI)



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

Figure 3. Temporal profile of the orientation analysis of two protein-lipid simulations, namely simulation 4 (sim 4) and 5 (sim 5). (A-C) 1, 2, and 3, angles. Profiles of the two PH domains of the ACAP1BAR-PH protein dimer are shown separately as P-M (blue) and P-W (green). As while one PH domain bound to the lipid bilayers, the another one remained unbound. Running averages were computed to show more clearly the evolution of the orientation angles.

209x78mm (300 x 300 DPI)


Page 42 of 48

Page 43 of 48

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


Figure 4. Angle distribution of all the equilibrated data from Figure 2-3. (A-C) Protein-solvent simulations. (D-F) Protein-lipid simulations. Angles at the P-W state are shown in green, and P-M in blue. Bin sizes of 2 were used.

209x152mm (300 x 300 DPI)



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

Figure 5. A simplified model of an ACAP1BAR-PH protein dimer during membrane binding process. Snapshot of an ACAP1BAR-PH dimer from (A) MD simulation of an ACAP1BAR-PH dimer in ionic solvent (B) MD simulation of an ACAP1BAR-PH dimer on lipid bilayers. (A-B) Two monomers of the protein dimer are colored in a sense that PH domain of the green monomer interacting mostly with water and ion molecules (P-W state) and the blue one interacting with the membrane (P-M state). Similar to Figure 1, the blue arrows define the orientation of the PH domains; the three yellow arrows in (B), which were calculated from the principle axes of the BAR domains, define the orientation of the BAR domains. 209x86mm (300 x 300 DPI)


Page 44 of 48

Page 45 of 48

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


Figure 6. (A) Snapshot of a PH domain of the ACAP1BAR-PH protein dimer from protein-lipid simulation 4. The PH domain is shown with cartoon representations in silver, except for the key residue loops. Three key residue loops defining the binding pocket interface were colored in lighter blue (Loop 1, residue 277-283), green (Loop 2, 300-305), and red (Loop 3, residue 321-327). Two bound PIP2 lipid head molecules are colored in orange (Pocket 1) and cyan (Pocket 2) to distinguish the binding pockets. Residues in contact with the PIP2 lipids are colored according to their side chain properties (magenta for negative, blue for positive, green for polar and white for non-polar). The BAR domain is shown with cartoon representations in transparent white. (B-C) Temporal profiles of separation between the pocket interface COM and the bound PIP2 lipid backbone phosphorus atoms. COM of Loop 1 and Loop 3 defines the Pocket 1 interface COM, and that of Loop 1 and Loop 2 defines the Pocket 2. (D) Number of contacts were counted every 50 ps for protein-lipid simulation 4 and then averaged over the number of frames. A contact is defined with a 12 Å cutoff between COM of the residue and the lipid backbone phosphorus atom. (E) The contact distribution in (D) projected onto the PH domain molecular surface, with blue corresponds to no contact and red corresponds to the most contacts. Loop 0, which locates between the two binding pockets, is circled. 209x161mm (300 x 300 DPI)



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

Figure 7. (A) Snapshot of Pocket 1 from protein-lipid simulation 4 at 255 ns. The blue sphere corresponds to the COM of Loop 1. The red sphere corresponds to the COM of two halves of Loop 3. The binding interface is represented as orange triangle with orange arrow indicating the normal. (B) Snapshot of Pocket 2 from protein-lipid simulation 4 at 418 ns. The blue sphere corresponds to the COM of Loop 1. The green sphere corresponds to the COM of two halves of Loop 2. The binding interface is represented as cyan triangle with cyan arrow indicating the normal. (A-B) Loop 1, 2, and 3 are colored in blue, green, and red respectively. The rest residues are colored in white. Underneath lipid backbone phosphorus atoms are represented as transparent spheres. The black arrow indicates the normal of the membrane. (C-D) Orientation of Pocket 1 and Pocket 2 in protein-lipid simulation 4. (E-F) Orientation of Pocket 1 and Pocket 2 in protein-lipid simulation 5. (C-F) Running averages were computed to show more clearly the evolution of the angles. Color convention refers to the previous orientation analysis of PH domains in Figure 3.


Page 46 of 48

Page 47 of 48

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


209x297mm (300 x 300 DPI)



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

Figure 8. (A) Snapshot of a PH domain at bound state with PIP2 lipids, from which we prepare all the initial structures for the US simulations. The atom coordinates were obtained from the last frame of protein-lipid simulation 4. The PH domain is shown in gray in cartoon representation; DOPC backbone phosphorus atoms of the lipid bilayers are represented by pink spheres. Two PIP2 lipids show in licorice representation indicate the position of PIP2 lipids in two independent US setups. The orange PIP2 lipid represents the initial bound state of the Pocket 1/PIP2 setup; the cyan PIP2 lipid represents that of the Pocket 2/PIP2 setup. The black arrow indicates the direction of pulling during SMD simulations for preparing the subsequent US windows. The black dotted line indicates the separation used for restraints of US windows, which is the separation between protein COM and lower-layer lipid backbone phosphorus atoms COM. (B) PMF profiles for PH domain of ACAP1BAR-PH protein dimer bound to PIP2 at Pocket 1 and Pocket 2. Protein-lipid separation was measured from the protein COM to the COM of lower layer lipid backbone phosphorus atoms along the membrane normal. Error estimates were obtained from bootstrap analysis.

189x112mm (300 x 300 DPI)


Page 48 of 48

Page 49 of 48

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


Figure 9. Three-dimensional histograms of orientation angle of (A) the Pocket 1/PIP2 US (B) the Pocket 2/PIP2 US throughout the 15 US windows. The angle refers to the orientation of the two binding interfaces. The restrained separations are the window centers of US simulations which measured the separation between protein COM and lower-layer lipid backbone phosphorus atoms COM. The bin size is 1⁰.

209x177mm (300 x 300 DPI)


Molecular Details of the PH Domain of ACAP1BAR-PH Protein

Recommend Documents