Skip to main content

National Institutes of Health

Eunice Kennedy Shriver National Institute of Child Health and Human Development

2020 Annual Report of the Division of Intramural Research

Structural and Chemical Biology of Membrane Proteins

Anirban Banerjee
  • Anirban Banerjee, PhD, Head, Unit on Structural and Chemical Biology of Membrane Proteins
  • Eric Christenson, PhD, Postdoctoral Fellow
  • Chul-jin Lee, PhD, Postdoctoral Fellow
  • Sayan Chakraborty, PhD, Visiting Fellow
  • Debabrata Dutta, PhD, Visiting Fellow
  • Robbins Puthenveetil, PhD, Visiting Fellow
  • Antony Yerabham, PhD, Visiting Fellow
  • Augustus Lowry, BA, Postbaccalaureate Fellow
  • Flowreen Shikwana, BA, Postbaccalaureate Fellow

Molecular mechanism of post-translational protein lipidation by DHHC palmitoyltransferases

Post-translational modifications greatly expand the structural, chemical, and functional diversity of the proteome. Of these, protein lipidation, which collectively refers to covalent modification of proteins by lipids, constitutes a centrally important class of post-translational modification. Protein palmitoylation, or more generally, protein S-acylation is a specific form of protein lipidation, whereby long-chain, typically C16, fatty acids become covalently attached to internal cysteines through a thioester linkage. Palmitoylation is one of the most pervasive and physiologically important post-translational modifications, and the targets of palmitoylation span a very wide range of proteins, ranging from ion channels to cell-surface receptors, neuronal scaffolding proteins, and small GTPases. The repertoire of palmitoylated proteins has expanded rapidly in recent years, with thousands of proteins now known to be part of the cellular “palmitoylome.” The physicochemical effect of palmitoylation is to alter the local hydrophobicity of the substrate protein. The thioester bond makes S-acylation unique in that it is a labile moiety and can be cleaved, in the cellular context, by thioesterase enzymes, which makes S-acylation one of the few dynamic post-translational modifications and unique among different forms of protein lipidation. The physiological effects of S-acylation are diverse and have critical cellular importance; for example, Ras, a small GTPase that is critical for cellular growth and differentiation and is mutated in about one-third of all human cancers, is palmitoylated at the Golgi and subsequently targeted to the plasma membrane by vesicular transport. Palmitoylated Ras localizes to cholesterol-rich domains on the plasma membrane. However, it is subsequently depalmitoylated by the thioesterase APT1, dissociates from the plasma membrane, and redistributes on endomembranes, including the Golgi. Such dynamic recycling of Ras is critical for its function.

Protein S-acylation is catalyzed by a large group of enzymes known as DHHC–palmitoyl acyltransferases (DHHC-PAT), so named because they contain a signature D-H-H-C motif (aspartate-histidine-histidine-cysteine) in a cysteine-rich domain (CRD) in an intracellular loop (Figure 1). These are low-abundance polytopic integral membrane proteins localized at a variety of cellular compartments. Humans have 23 DHHC-PATs encoded in their genome. Beyond the shared DHHC domain, DHHC-PATs vary considerably; some possess ankyrin repeats (structural protein motifs that mediate protein-protein interactions), a few have six transmembrane helices instead of the usual four, and at least one forms a functional heterodimer with a cytoplasmic auxiliary subunit (Figure 1). To date, no consensus sequences have been reported for palmitoylation. A specific DHHC-PAT can palmitoylate many substrates, and, conversely, a given substrate can be palmitoylated by many DHHC-PATs. Such redundancy has been one of the most intriguing aspects of DHHC-PATs and makes it difficult to assign substrates by overexpression/knockout strategies, given that, in the absence of one DHHC-PAT enzyme, others can take over. However, this does not necessarily reflect the true enzyme-substrate relationship. The situation has been even more confounded by the lack of specific inhibitors of DHHC-PATs. Even though 2-bromopalmitate is widely used as a global inhibitor of DHHC-PATs, it has been shown that it also broadly targets other proteins involved in lipid metabolism.

Figure 1. Organization and properties of DHHC palmitoyltransferases (PATs)

Figure 1

Click image to view.

a. The organization of three different DHHC-PATs are shown schematically. The spheres indicate protein-protein interaction domains. Erf2 associates with a cytoplasmic subunit, Erf4, to form the active enzyme.

b. The DHHC-CRD region of a few representative DHHC-PATs are aligned. The conserved amino acids are shown in red.

c. Reaction catalyzed by DHHC-PATs; the reverse reaction is catalyzed by acylprotein thioesterases (APT).

Besides its broad importance in cell biology, palmitoylation has been linked to several diseases, most notably neuropsychiatric disorders such as Huntington’s disease and various forms of cancer. Recently, it was shown that DHHC20 palmitoylates epidermal growth factor (EGFR) and is thus a potential therapeutic target for a wide range of cancers. More recently, DHHC3 has been proposed as a target for cancer treatment owing to its activity as the palmitoyltransferase for programmed cell-death ligand 1(PD-L1). However, when we started working on this family, very little was known about the molecular mechanism of DHHC palmitoyltransferases, despite their importance across a broad spectrum of biological pathways and their biomedical importance. Nothing was known about their structural organization or how they interact with substrates and the fatty acyl coenzyme A (CoA), which serves as the acyl donor.

In a major breakthrough in this field, we had earlier solved the high-resolution crystal structures of two members of the DHHC family, human DHHC20 (hDHHC20) and zebrafish DHHC15 (Figure 2a). They are the first structures of any member of this family to be characterized and reveal a tepee-like transmembrane domain organization, which splays apart towards the cytoplasmic side and harbors the active site at the membrane-aqueous interfacial region (Figure 2b), thus readily explaining why membrane-proximal cysteines are palmitoylated. We also solved the structure of hDHHC20 irreversibly modified by a covalent inhibitor, 2-bromopalmitate. The structure mimics the auto-acylated intermediate state in the enzymatic pathway and thus reveals how the acyl group of fatty acyl–CoA binds in a cavity formed in the bilayer by the transmembrane domain (Figure 2c). Residues lining the cavity contact the acyl chain, and their mutation affects enzymatic activity. By mutating two residues at the tapering end of the cavity, we also showed that we can change the acyl chain-length selectivity of the mutant enzymes (Figure 2d). Thus, the cavity functions as a molecular ruler in determining the acyl chain-length selectivity of hDHHC20, important because, although palmitate is the most prevalent fatty acid used by DHHC palmitoyltransferases, they can use fatty acyl–CoAs of other chain lengths, a property that varies between different members of the DHHC family.

Figure 2. Structure, function, and membrane deformation of DHHC palmitoyltransferases

Figure 2

Click image to view.

a. The structure of human DHHC20 and a catalytically inactive mutant of zebrafish DHHC15, shown in ribbon trace. The transmembrane domain (TM) is shown in green, the DHHC–containing cysteine-rich domain in blue and the C-terminal domain in orange. The grey spheres indicate Zn2+ ions. These are both Golgi-resident enzymes, and thus the top side faces the Golgi lumen and the active site the cytoplasm.

b. Active site of human DHHC20 showing the catalytic triad containing the active-site cysteine.

c. Structure of human DHHC20 irreversibly modified with 2-bromopalmitate, which results in the active-site cysteine linking to the alpha-carbon of palmitic acid. The acyl group of palmitic acid is shown in stick rendition. Also shown are two residues towards the top of the tapering cavity to where the palmitate binds.

d. The acyl chain-length selectivity of wild-type (WT) human DHHC20. Mutation of tyrosine181 to alanine (Y181A) expands the cavity and shifts the acyl selectivity to the longer side. On the other hand, mutation of serine29 to phenylalanine (S29F) contracts the cavity and thus shifts the acyl selectivity to the shorter side.

e. 3D density maps for the POPC alkyl-chain core (yellow surface) near the protein surface (blue), calculated from the molecular dynamics (MD) trajectories. The catalytic cysteine (Cys156) is highlighted (magenta). The four regions where the membrane is deformed are indicated. The approximate minimum and maximum widths of the alkyl-chain core of the bilayer are also indicated, showing that the maximum deformation is around the active site.

In the past year, we focused on structure-based understanding of several aspects of DHHC function. Starting from our high-resolution structure and using molecular dynamics simulations and in collaboration with José Faraldo-Gómez, we showed that the protein surface of hDHHC20 severely distorts the membrane, which causes the active site to be exposed to the aqueous side of the membrane. Otherwise, the active site would have been buried in the hydrophobic part of the membrane and would have impaired catalysis, given that this kind of catalysis involves charge separation. Using a combination of X-ray crystallography and atomistic molecular dynamics simulations, we also determined a high-resolution structure of human DHHC20 in complex with an intact palmitoyl CoA. The structure represents the pre-catalytic complex of hDHHC20 with palmitoyl CoA and reveals how the polar CoA part of the fatty acyl–CoA is recognized by hDHHC20. Mutation of residues involved in binding to CoA impair catalysis, thus lending credence to our structural model. We also showed that recognition of palmitoyl CoA by hDHHC20 depends on the coincident recognition of both the fatty acyl chain and the CoA headgroup, thus explaining how high concentrations of free intracellular CoA cannot inhibit DHHC enzymes.

Molecular mechanism of Porcupine, an integral membrane enzyme of the MBOAT family, which catalyzes Wnt lipidation

In metazoans, Wnt proteins regulate many processes, including cellular growth, differentiation, and tissue homeostasis, through the highly conserved Wnt signaling pathway. Porcupine is an endoplasmic reticulum (ER)–resident integral membrane enzyme, which catalyzes post-translational modification of Wnts with the unsaturated lipid palmitoleic acid. This unique form of lipidation with palmitoleic acid is a vital step in the biogenesis and secretion of Wnt. Porcupine inhibitors are currently in clinical trials for cancer treatment. However, Porcupine-mediated Wnt lipidation has not been reconstituted in vitro with purified enzyme. We recently reported the first successful purification of human Porcupine and confirmed, through in vitro reconstitution with the purified enzyme, that Porcupine is necessary and sufficient for Wnt acylation. By systematically examining a series of substrate variants, we showed that Porcupine intimately recognizes the local structure of Wnt around the site of acylation. Our in vitro assay enabled us to examine the activity of Porcupine with a range of fatty acyl–CoAs of varying length and unsaturation. The selectivity of human Porcupine across a spectrum of fatty acyl–CoAs suggested that the kink in the unsaturated acyl chain is a key determinant in Porcupine-mediated catalysis. We also showed that two putative Porcupine inhibitors, which were discovered with cell-based assays, indeed target human Porcupine. Together, the results provide several, high-resolution biochemical insights into the mechanism of Porcupine-mediated Wnt acylation and pave the way for further detailed biochemical and structural studies.

Molecular mechanism of iron transport across cellular membranes

The importance of iron in biology cannot be overstated. In higher organisms, mitochondria are the ‘hotspot’ for the cell biology of iron, because Fe-S clusters are biosynthesized and iron is inserted into heme there. Mitochondrial iron homeostasis plays a critical role in cellular iron homeostasis and in the overall physiology of the cell. In vertebrates, the only known major transporters of iron into mitochondria are mitoferrin-1 and mitoferrin-2, two homologous members of a large group of mitochondrial transporters known as the Mitochondrial Carrier family (Figures 3a and 3b). Mitoferrin-1 (Mfrn1) is expressed mainly in erythroid cells, while mitoferrin-2 is expressed ubiquitously. Knockout of Mfrn1 is embryonically lethal, reflecting the importance of mitoferrins in vertebrate physiology.

Mfrn1 and Mfrn2 were discovered more than 10 years ago. However, the proposed iron-transport activity had not been demonstrated using an in vitro functional reconstitution assay. Also no report about their interaction with iron or other related metal ions existed, most likely because heterologous overexpression and purification of mitoferrins were not reported in the literature. We carried out heterologous purification, in vitro functional reconstitution, and mutational dissection of a vertebrate Mfrn1 (Figure 3c), the first demonstration that Mfrn1 can indeed transport iron. We showed that Mfrn1 is a promiscuous metal ion transporter in that it also transports other first-row transition metal ions. Through mutagenesis, we discovered candidate residues that are important for metal-ion transport by Mfrn1 and those that could be involved in forming metal-ion binding sites during transport. Our studies provided the first biochemical insights into Mfrn function and form the starting point for future high-resolution structural studies of Mfrn function. Our transport assay and the purification strategy will lead to more detailed biochemical and biophysical experiments into the mechanistic basis of iron transport by Mfrn1 and Mfrn2.

Figure 3. Iron transport by Mitoferrin and MavN

Figure 3

Click image to view.

a. Iron is imported through the plasma membrane by the Transferrin/Transferrin Receptor (blue) cycle and is transported out of endosomes by the divalent metal ion transporter (DMT) (grey). Iron is delivered to Mitoferrin (yellow cylinders) by unknown means. Mitoferrin delivers iron to unknown partners in mitochondria, which become available for heme and Fe-S cluster biosynthesis.

b. Schematic depiction of the simplest transporter cycle and the different putative states of Mitoferrin involved.

c. Assay of iron-transport activity by TMfrn1, a Mitoferrin1 homolog. Representative PGSK (fluorescent reporter of iron) quenching curves upon addition of iron to TMfrn1 proteoliposome. Red traces are increasing concentrations of iron. Black dotted traces show protein-free liposomes.

d. Schematic showing Legionella entering a host cell and sequestering itself in a Legionella-containing vacuole (LCV). MavN is inserted into the membrane of the LCV and hijacks iron from the host cell.

e. Assay showing iron transport by proteoliposome-reconstituted MavN, using the identical assay as in c.

In the past year, we used our in vitro proteoliposome-reconstituted iron-transport assay, the first such assay to be reported in the literature, to dissect the iron-transport activity of MavN, another proposed iron transporter, in the bacterial pathogen Legionella pneumophila. Legionella is an intracellular pathogen that enters the host cell and sequesters itself in the Legionella-Containing Vacuole (LCV), where it survives and proliferates. This is achieved by secreting hundreds of so-called "effector proteins," which are secreted through the Type IV secretion system (T4SS) into the host cell. Of the hundreds of effector proteins, there are seven core effector proteins, and MavN is one of them. MavN is inserted into the membrane of the LCV. Using genetic and cell-based studies, the laboratory of our collaborator Ralph Isberg proposed that MavN is an iron transporter. By heterologous purification and proteoliposome reconstitution, we showed that MavN is indeed an iron transporter. Mutations in MavN that impaired iron transport in vitro also impaired Legionella growth in a cell-based assay, pointing to the importance of MavN for Legionella.

Structure and molecular mechanism of ATG9, the only transmembrane component of the core autophagy machinery

Autophagy is a critical process in both health and disease. Serving mainly to adapt organisms to a diverse range of conditions such as metabolic stress, pathogenic infection, and ageing, autophagy is conserved from yeast to humans. In addition, autophagy is involved in a large number of human diseases, including neurodegenerative disorders, various forms of cancer, and inflammatory diseases. Although autophagy was discovered more than 50 years ago, some its fundamental aspects remain rather poorly understood and are only coming to light now. An outstanding example of this lack of knowledge is the structure and function of ATG9, the only essential transmembrane protein in the pathway. ATG9 localizes to small, 30–60 nm vesicles known as “ATG9 vesicles,” which are critically important for the expansion of the autophagosome membrane. ATG9 vesicles traffic between the trans-Golgi network (TGN) and pre-autophagosomal structures in response to stimuli that initiate autophagy. Despite the clear and obvious importance of ATG9 in autophagy, its atomic structure and molecular function remain unknown. It has been speculated that ATG9 transports lipids that contribute to the growth of the double-membrane structure of the autophagosome. However, progress in this regard has been completely thwarted by the lack of an atomic structure that could serve as a starting point for a range of experiments, from cell-biological to biochemical, to elucidate its function.

In order to address this knowledge gap, we initiated a collaboration with the labs of Juan Bonifacino and Jiansen Jiang. This led to a 2.9 Å–resolution structure of the ubiquitous human ATG9A isoform by cryo-EM (Figure 4a). The structure revealed that ATG9A adopts a unique fold (Figure 4b, 4c), assembling into a domain-swapped homotrimer, with each wedge-shaped protomer comprising four transmembrane α-helices and two α-helices that are only partially embedded in the membrane. The transmembrane domain is capped by a cytosolic domain, which is mostly α-helical. Classification of the cryo-EM imaging revealed two predominant ATG9A conformers (states A and B). An important feature of the ATG9A structure is the presence of a branched network of pores through the protein, with outlets facing the cytosol, the outer leaflet lipid head groups, and the lumen (Figure 4d). Cell-based functional assays demonstrated importance of pore-lining residues, suggesting an essential function of the pore. The dimensions of the pore are sufficient to accommodate phospholipids and likely permit lipid transport from or to a lipid chaperone such as ATG2. Alternatively ATG9 may act as a lipid scramblase, distributing the lipids that arrive through ATG2 among two leaflets of the bilayer. The branched pores may also function as solvent conduits to relieve osmotic pressures in growing and shrinking ATG9 vesicles. ATG9A directly interacts with ATG2A, interactions that are mediated by the ATG9A cytosolic domain, deletion of which abrogates both ATG2A co-immunoprecipitation and autophagosome maturation. Because ATG2A has demonstrated lipid-transport capability in vitro, the plausible functions of ATG9A are to: (1) facilitate movement of lipids across the bilayer of the growing phagophore, which could be achieved by using the branched pores we see in our structure for moving polar lipid headgroups or by distorting the bilayer, lowering the energetic barrier for trans-bilayer lipid movement; (2) act as a membrane-anchored “lasso” to capture and accurately target ATG2A or other lipid chaperones. It is interesting to note that two helices, which plug the central pore of the trimer, are discernible only in state A, indicating that flexibility of the putative “lasso” could be partly dependent on the conformation of the transmembrane core.

The early stages of autophagy involve a reorganization of intracellular membranes to form the nascent autophagosome. ATG9 is essential for this process, and our structure provides a framework for interrogating its function in bilayer remodeling. Using molecular dynamics simulations in collaboration with the lab of José Faraldo-Gómez, and based on the cryo-EM data, we found that the architecture of the ATG9A trimer elicits long-range positive membrane curvature, owing to the geometric features and amino-acid make-up of the protein surface exposed to the lipid. The simulations also show that co-localization of multiple ATG9A trimers greatly amplifies this effect and induces macroscopic changes in membrane morphology (Figure 4e, 4f). The finding corresponds with observations that ATG9 is present in small (30–60 nm) vesicles and thin tubules, and, at least in yeast, localizes to the edges of the growing phagophore. Interestingly, states A and B differ in the tilt-angle of the protomers relative to the membrane normal, implying a tolerance for a range of bilayer curvatures. This leads to the intriguing possibility that ATG9A–induced curvature may modulate the energetic scale and effect lipid diffusion out of the vesicles into another membrane.

Figure 4

Click image to view.

Figure 4. Structure and membrane interactions of human ATG9A

a. Cryo-EM density map of ATG9A in state A is shown in side and top views and colored according to the protomer (cyan, magenta, and orange). Translucent surface shows disordered detergent molecules surrounding the transmembrane surface of the trimer.

b. Structure of an isolated protomer (cyan) and two alpha helices from another protomer (magenta). Numbers correspond to alpha helices.

c. Topology of the ATG9A protomer showing numbered alpha helices as rounded rectangles and beta strands as filled arrows.

d. Branched network of internal pores in the ATG9A trimer shown by electrostatic potential surface of the central pore, lateral branch, and perpendicular branch (with 180° rotation are shown). Green, yellow, and pink arrows are used to show the positions of the central pore, lateral branch, and perpendicular branch, respectively.

e. and f. Molecular dynamics simulation of a hypothetical ATG9A cluster in a POPC bilayer.

e. Configuration of the system at the beginning of the simulation; 36 ATG9A proteins (gray) are arranged around a central trimer in three concentric hexagons (inset). Note the proteins are not in contact. The system amounts to about 11,220,000 particles.

f. At the end of the calculated trajectory of 1 ms.

Additional Funding

  • Office of Aids Research (OAR) Innovation award to Anirban Banerjee (2019)
  • NICHD DIR Director's Investigator Award

Publications

  1. Stix R, Song J, Banerjee A, Faraldo-Gómez J. DHHC20 palmitoyl-transferase reshapes the membrane to foster catalysis. biorXiv 2019;https://doi.org/10.1101/808469.
  2. Christenson ET, Isaac DT, Yoshida K, Lipo E, Kim JS, Ghirlando R, Isberg RR, Banerjee A. The iron-regulated vacuolar Legionella pneumophila MavN protein is a transition-metal transporter. Proc Natl Acad Sci USA 2019;116:17775-17785.
  3. Lee CJ, Rana MS, Bae C, Li Y, Banerjee A. In vitro reconstitution of Wnt acylation reveals structural determinants of substrate recognition by the acyltransferase human Porcupine. J Biol Chem 2019;294:231-245.
  4. Rana M, Lee C, Banerjee A. The molecular mechanism of DHHC palmitoyltransferases. Biochem Soc Transactions 2019;47:157-167.
  5. Rana MS, Kumar P, Lee CJ, Verardi R, Rajashankar KR, Banerjee A. Fatty acyl recognition and transfer by an integral membrane S-acyltransferase. Science 2018;359:6372.
  6. Christenson ET, Gallegos AS, Banerjee A. In vitro reconstitution, functional dissection, and mutational analysis of metal ion transport by mitoferrin-1. J Biol Chem 2018;293:3819-3828.

Collaborators

  • José Faraldo-Gómez, PhD, Theoretical Molecular Biophysics Section, NHLBI, Bethesda, MD
  • Rodolfo Ghirlando, PhD, Laboratory of Molecular Biology, NIDDK, Bethesda, MD
  • James Inglese, PhD, Assay Development & Screening Technology Laboratory, NCATS, Bethesda, MD
  • Ralph Isberg, PhD, Tufts University School of Medicine, Boston, MA

Contact

For more information, email anirban.banerjee@nih.gov or visit http://banerjee.nichd.nih.gov.

Top of Page