Metadynamics simulations reveal mechanisms of Na+ and Ca2+ transport in two open states of the channelrhodopsin chimera, C1C2.

Cation conducting channelrhodopsins (ChRs) are a popular tool used in optogenetics to control the activity of excitable cells and tissues using light. ChRs with altered ion selectivity are in high demand for use in different cell types and for other specialized applications. However, a detailed mechanism of ion permeation in ChRs is not fully resolved. Here, we use complementary experimental and computational methods to uncover the mechanisms of cation transport and valence selectivity through the channelrhodopsin chimera, C1C2, in the high- and low-conducting open states. Electrophysiology measurements identified a single-residue substitution within the central gate, N297D, that increased Ca2+ permeability vs. Na+ by nearly two-fold at peak current, but less so at stationary current. We then developed molecular models of dimeric wild-type C1C2 and N297D mutant channels in both open states and calculated the PMF profiles for Na+ and Ca2+ permeation through each protein using well-tempered/multiple-walker metadynamics. Results of these studies agree well with experimental measurements and demonstrate that the pore entrance on the extracellular side differs from original predictions and is actually located in a gap between helices I and II. Cation transport occurs via a relay mechanism where cations are passed between flexible carboxylate sidechains lining the full length of the pore by sidechain swinging, like a monkey swinging on vines. In the mutant channel, residue D297 enhances Ca2+ permeability by mediating the handoff between the central and cytosolic binding sites via direct coordination and sidechain swinging. We also found that altered cation binding affinities at both the extracellular entrance and central binding sites underly the distinct transport properties of the low-conducting open state. This work significantly advances our understanding of ion selectivity and permeation in cation channelrhodopsins and provides the insights needed for successful development of new ion-selective optogenetic tools.

Single-molecule quantification of the oligomeric state of ZIP transporters in mammalian cells with fluorescence correlation spectroscopy.

The SLC39 family of transporters, otherwise known as ZIPs for Zrt and Irt-like Proteins, function to increase cytosolic levels of transition metals. ZIP transporters have been identified at all phylogenetic levels and are members of the SoLute Carrier (SLC) superfamily. There are fourteen ZIP transporters encoded in the human genome. ZIP transmembrane proteins are expressed in the plasma membrane or membranes of intracellular organelles and have unique expression profiles across cell types. While direct structural efforts including x-ray crystallography, NMR and ab initio approaches have been effective tools in elucidating the structure of ZIPs, direct elucidation of the oligomeric state of these proteins is essential in understanding how wild type ZIP proteins function and whether mutations alter the oligomeric state of ZIPs. Unfortunately, several tools to quantify oligomeric states of proteins require overexpression of proteins which can lead to artifacts in experimental results. In contrast, fluorescence correlation spectroscopy (FCS) is a single-molecule technique which can be used to quantify the oligomeric state of transmembrane proteins. FCS takes advantage of the observation that the molecular brightness of a cluster of fluorescent molecules is directly proportional to the number of fluorescent molecules within the protein complex. This chapter describes how to implement FCS, focused on ZIP transporters, to quantify the oligomeric state of transmembrane in vivo. Included within this chapter are procedures to design constructs for experiments, transfection of mammalian cells as well as data acquisition and analysis. Taken together, FCS is a powerful mechanism to investigate the oligomeric state of proteins embedded within membranes of cells.

The transmembrane domains mediate oligomerization of the human ZIP4 transporter in vivo.

The human (h) ZIP4 is a plasma membrane transporter that functions to increase cytosolic zinc levels. hZIP4 encodes eight transmembrane domains and a large extracellular domain (ECD). This ECD is cleaved from the holo-transporter when cells are zinc-deficient. At the same time, mutations in the ECD can result in the zinc-deficiency disease Acrodermatitis enteropathica. Previously, it was shown that hZIP4's ECD is comprised of two structurally independent subdomains where contacts between the ECD monomeric units are centered at the PAL motif. These results lead to the hypothesis that ZIP4-ECD is essential to the dimerization of the holo-transporter. To test this hypothesis, we used Fluorescence Correlation Spectroscopy (FCS) to quantify the oligomeric state of full-length hZIP4 and hZIP4 lacking the ECD domain, each tagged with eGFP. Inspection of our experimental results demonstrate that both the full-length and truncated hZIP4 is a dimer when expressed in HEK293 cells. Parallel functional experiments demonstrate that the Km and Vmax for truncated and full-length hZIP4/eGFP are similar. Determining that truncated hZIP4/eGFP forms a dimer is a crucial step for understanding the function of the hZIP4-ECD, which provides more insight into how the diseases related to hZIP4 protein.

Heterologous Expression of Full-Length and Truncated Human ZIP4 Zinc Transporter in Saccharomyces cerevisiae.

The human (h) transporter hZIP4 is the primary Zn2+ importer in the intestine. hZIP4 is also expressed in a variety of organs such as the pancreas and brain. Dysfunction of hZIP4 can result in the Zn2+ deficiency disease acrodermatitis enteropathica (AE). AE can disrupt digestive and immune system homeostasis. A limited number of hZIP4 expression strategies have hindered increasing knowledge about this essential transmembrane protein. Here, we report the heterologous expression of hZIP4 in Saccharomyces cerevisiae. Both a wild-type and a mutant S. cerevisiae strain, in which the endogenous Zn2+ transporters were deleted, were used to test the expression and localization of an hZIP4-GFP fusion protein. A full-length hZIP4-GFP and a truncated membrane-domain-only (mhZIP4-GFP) protein were observed to be present in the plasma membrane in yeast.

Channelrhodopsin C1C2: Photocycle kinetics and interactions near the central gate.

Channelrhodopsins (ChR) are light-sensitive cation channels used in optogenetics, a technique that applies light to control cells (e.g., neurons) that have been modified genetically to express those channels. Although mutations are known to affect pore kinetics, little is known about how mutations induce changes at the molecular scale. To address this issue, we first measured channel opening and closing rates of a ChR chimera (C1C2) and selected variants (N297D, N297V, and V125L). Then, we used atomistic simulations to correlate those rates with changes in pore structure, hydration, and chemical interactions among key gating residues of C1C2 in both closed and open states. Overall, the experimental results show that C1C2 and its mutants do not behave like ChR2 or its analogous variants, except V125L, making C1C2 a unique channel. Our atomistic simulations confirmed that opening of the channel and initial hydration of the gating regions between helices I, II, III, and VII of the channel occurs with 1) the presence of 13-cis retinal; 2) deprotonation of a glutamic acid gating residue, E129; and 3) subsequent weakening of the central gate hydrogen bond between the same glutamic acid E129 and asparagine N297 in the central region of the pore. Also, an aspartate (D292) is the unambiguous primary proton acceptor for the retinal Schiff base in the hydrated channel.

Characterizing Channelrhodopsin Channel Properties Via Two-Electrode Voltage Clamp and Kinetic Modeling.

Two-electrode voltage clamp (TEVC) is a preferred electrophysiological technique used to study gating kinetics and ion selectivity of light-activated channelrhodopsins (ChRs). The method uses two intracellular microelectrodes to hold, or clamp, the membrane potential at a specific value and measure the flow of ions across the plasma membrane. Here, we describe the use of TEVC and a simple solution exchange protocol to measure cation selectivity and analyze gating kinetics of the C1C2 chimera expressed in Xenopus laevis oocytes. Detailed instructions on how to process the collected data and interpret the results are also provided.

Quantifying the Oligomeric State of hZIP4 on the Surface of Cells.

The human (h) zinc transporter ZIP4 is expressed on the plasma membrane and functions to increase cytosolic zinc levels. Mutations in hZIP4 cause the disease acrodermatitis enteropathica. Dysfunction in the regulation of hZIP4 has also been indicated in solid tissue cancers, including pancreatic and prostate cancer. Although structural studies of the extracellular domain and computational modeling of the membrane domain suggest hZIP4 exists as a dimer, the oligomerization status of hZIP4 in the plasma membrane of mammalian cells has not been directly quantified in vivo. Here, the oligomeric state of hZIP4 expressed in HEK293 cells was quantified using fluorescence correlation spectroscopy. hZIP4 was tagged with eGFP, and by comparing brightness values (ε) of monomer and tandem eGFP constructs to that of an hZIP4/eGFP, we show that hZIP4 is a dimer. Determining that hZIP4 is a dimer is an important step toward understanding the function and processing of the protein, which can provide more insight into how diseases affected by hZIP4 occur and can be managed.

Concomitant disorder and high-affinity zinc binding in the human zinc- and iron-regulated transport protein 4 intracellular loop.

The human zinc- and iron-regulated transport protein 4 (hZIP4) protein is the major plasma membrane protein responsible for the uptake of zinc in the body, and as such it plays a key role in cellular zinc homeostasis. hZIP4 plasma membrane levels are regulated through post-translational modification of its large, disordered, histidine-rich cytosolic loop (ICL2) in response to intracellular zinc concentrations. Here, structural characteristics of the isolated disordered loop region, both in the absence and presence of zinc, were investigated using nuclear magnetic resonance (NMR) spectroscopy. NMR chemical shifts, coupling constants and temperature coefficients of the apoprotein, are consistent with a random coil with minor propensities for transient polyproline Type II helices and ß-strand in regions implicated in post-translational modifications. The ICL2 protein remains disordered upon zinc binding, which induces exchange broadening. Paramagnetic relaxation enhancement experiments reveal that the histidine-rich region in the apoprotein makes transient tertiary contacts with predicted post-translational modification sites. The residue-specific data presented here strengthen the relationship between hZIP4 post-translational modifications, which impact its role in cellular zinc homeostasis, and zinc sensing by the intracellular loop. Furthermore, the zinc sensing mechanism employed by the ICL2 protein demonstrates that high-affinity interactions can occur in the presence of conformational disorder.

An optimized and automated approach to quantifying channelrhodopsin photocurrent kinetics.

Channelrhodopsins are light-activated ion channels that enable targetable activation or inhibition of excitable cells with light. Ion conductance can generally be described by a four step photocycle, which includes two open and two closed states. While a complete understanding of channelrhodopsin function cannot be understood in the absence of kinetic modeling, model fitting requires manual fitting, which is laborious and technically complicated for non-experts. To enhance analysis of photocurrent data, this manuscript describes a fitting program where electrophysiology data can be automatically and quantitatively analyzed. Significant improvement in this program when compared to our previous version includes 1) the ability to automatically find the experiment start time using the derivative of the current signal, 2) utilizing the Object Oriented Programing (OPP) paradigm which is significantly more reliable if the code is used by people with little to no programming experience and 3) the distribution of the code is simplified to sharing a single MATLAB file, including rigorous comments throughout. To demonstrate the utility of this program, we show automated fitting of photocurrents from two member proteins: channelrhodopsin-2 and a chimera between channelrhodopsin-1 and channelrhodopsin-2 (C1C2).

The emerging role of zinc transporters in cellular homeostasis and cancer.

Zinc is an essential micronutrient that plays a role in the structural or enzymatic functions of many cellular proteins. Cellular zinc homeostasis involves the opposing action of two families of metal transporters: the ZnT (SLC30) family that functions to reduce cytoplasmic zinc concentrations and the ZIP (SLC39) family that functions to increase cytoplasmic zinc concentrations. Fluctuations in intracellular zinc levels mediated by these transporter families affect signaling pathways involved in normal cell development, growth, differentiation and death. Consequently, changes in zinc transporter localization and function resulting in zinc dyshomeostasis have pathophysiological effects. Zinc dyshomeostasis has been implicated in the progression of cancer. Here we review recent progress toward understanding the structural basis for zinc transport by ZnT and ZIP family proteins, as well as highlight the roles of zinc as a signaling molecule in physiological conditions and in various cancers. As zinc is emerging as an important signaling molecule in the development and progression of cancer, the ZnT and ZIP transporters that regulate cellular zinc homeostasis are promising candidates for targeted cancer therapy.

Adjacent channelrhodopsin-2 residues within transmembranes 2 and 7 regulate cation selectivity and distribution of the two open states.

Channelrhodopsin-2 (ChR2) is a light-activated channel that can conduct cations of multiple valencies down the electrochemical gradient. Under continuous light exposure, ChR2 transitions from a high-conducting open state (O1) to a low-conducting open state (O2) with differing ion selectivity. The molecular basis for the O1 → O2 transition and how ChR2 modulates selectivity between states is currently unresolved. To this end, we used steered molecular dynamics, electrophysiology, and kinetic modeling to identify residues that contribute to gating and selectivity in discrete open states. Analysis of steered molecular dynamics experiments identified three transmembrane residues (Val-86, Lys-93, and Asn-258) that form a putative barrier to ion translocation. Kinetic modeling of photocurrents generated from ChR2 proteins with conservative mutations at these positions demonstrated that these residues contribute to cation selectivity (Val-86 and Asn-258), the transition between the two open states (Val-86), open channel stability, and the hydrogen-bonding network (K93I and K93N). These results suggest that this approach can be used to identify residues that contribute to the open-state transitions and the discrete ion selectivity within these states. With the rise of ChR2 use in optogenetics, it will be critical to identify residues that contribute to O1 or O2 selectivity and gating to minimize undesirable effects.

Zn(2+) at a cellular crossroads.

Zinc is an essential micronutrient for cellular homeostasis. Initially proposed to only contribute to cellular viability through structural roles and non-redox catalysis, advances in quantifying changes in nM and pM quantities of Zn(2+) have elucidated increasing functions as an important signaling molecule. This includes Zn(2+)-mediated regulation of transcription factors and subsequent protein expression, storage and release of intracellular compartments of zinc quanta into the extracellular space which modulates plasma membrane protein function, as well as intracellular signaling pathways which contribute to the immune response. This review highlights some recent advances in our understanding of zinc signaling.

Voltage Clamp Fluorometry of P-Type ATPases.

Voltage clamp fluorometry has become a powerful tool to compare partial reactions of P-type ATPases such as the Na(+),K(+)-ATPase and H(+),K(+)-ATPase with conformational dynamics of these ion pumps. Here, we describe the methodology to heterologously express membrane proteins in X. laevis oocytes and site-specifically label these proteins with one or more fluorophores. Fluorescence changes are measured simultaneously with current measurements under two-electrode voltage clamp conditions.

A Zinc(II) Photocage Based on a Decarboxylation Metal Ion Release Mechanism for Investigating Homeostasis and Biological Signaling.

Metal ion signaling in biology has been studied extensively with ortho-nitrobenzyl photocages; however, the low quantum yields and other optical properties are not ideal for these applications. We describe the synthesis and characterization of NTAdeCage, the first member in a new class of Zn(2+) photocages that utilizes a light-driven decarboxylation reaction in the metal ion release mechanism. NTAdeCage binds Zn(2+) with sub-pM affinity using a modified nitrilotriacetate chelator and exhibits an almost 6 order of magnitude decrease in metal binding affinity upon uncaging. In contrast to other metal ion photocages, NTAdeCage and the corresponding Zn(2+) complex undergo efficient photolysis with quantum yields approaching 30 %. The ability of NTAdeCage to mediate the uptake of (65) Zn(2+) by Xenopus laevis oocytes expressing hZIP4 demonstrates the viability of this photocaging strategy to execute biological assays.

Cysteine Substitution and Labeling Provide Insight into Channelrhodopsin-2 Ion Conductance.

Channelrhodopsin-2 is a light-activated cation channel. However, the mechanism of ion conductance is unresolved. Here, we performed cysteine scanning mutagenesis on transmembrane domain 7 followed by labeling with a methanethiosulfonate compound. Analysis of our results shows that residues that line the putative pore and interface with adjacent transmembrane domains 1 and 3, as proposed by our channelrhodopsin-2 homology model, affect ion conductance, decay kinetics, and/or off kinetics. Combined, these results suggest that negative charges at the extracellular side of transmembrane domain 7 funnel cations into the pore.

Computation and Functional Studies Provide a Model for the Structure of the Zinc Transporter hZIP4.

Members of the Zrt and Irt protein (ZIP) family are a central participant in transition metal homeostasis as they function to increase the cytosolic concentration of zinc and/or iron. However, the lack of a crystal structure hinders elucidation of the molecular mechanism of ZIP proteins. Here, we employed GREMLIN, a co-evolution-based contact prediction approach in conjunction with the Rosetta structure prediction program to construct a structural model of the human (h) ZIP4 transporter. The predicted contact data are best fit by modeling hZIP4 as a dimer. Mutagenesis of residues that comprise a central putative hZIP4 transmembrane transition metal coordination site in the structural model alter the kinetics and specificity of hZIP4. Comparison of the hZIP4 dimer model to all known membrane protein structures identifies the 12-transmembrane monomeric Piriformospora indica phosphate transporter (PiPT), a member of the major facilitator superfamily (MFS), as a likely structural homolog.

The large intracellular loop of hZIP4 is an intrinsically disordered zinc binding domain.

The human (h) ZIP4 transporter is a plasma membrane protein which functions to increase the cytosolic concentration of zinc. hZIP4 transports zinc into intestinal cells and therefore has a central role in the absorption of dietary zinc. hZIP4 has eight transmembrane domains and encodes a large intracellular loop between transmembrane domains III and IV, M3M4. Previously, it has been postulated that this domain regulates hZIP4 levels in the plasma membrane in a zinc-dependent manner. The objective of this research was to examine the zinc binding properties of the large intracellular loop of hZIP4. Therefore, we have recombinantly expressed and purified M3M4 and showed that this domain binds two zinc ions. Using a combination of site-directed mutagenesis, metal binding affinity assays, and X-ray absorption spectroscopy, we demonstrated that the two Zn(2+) ions bind sequentially, with the first Zn(2+) binding to a CysHis3 site with a nanomolar binding affinity, and the second Zn(2+) binding to a His4 site with a weaker affinity. Circular dichroism spectroscopy revealed that the M3M4 domain is intrinsically disordered, with only a small structural change induced upon Zn(2+) coordination. Our data supports a model in which the intracellular M3M4 domain senses high cytosolic Zn(2+) concentrations and regulates the plasma membrane levels of the hZIP4 transporter in response to Zn(2+) binding.

Transmembrane domain three contributes to the ion conductance pathway of channelrhodopsin-2.

Channelrhodopsin-2 (ChR2) is a light-activated nonselective cation channel that is found in the eyespot of the unicellular green alga Chlamydomonas reinhardtii. Despite the wide employment of this protein to control the membrane potential of excitable membranes, the molecular determinants that define the unique ion conductance properties of this protein are not well understood. To elucidate the cation permeability pathway of ion conductance, we performed cysteine scanning mutagenesis of transmembrane domain three followed by labeling with methanethiosulfonate derivatives. An analysis of our experimental results as modeled onto the crystal structure of the C1C2 chimera demonstrate that the ion permeation pathway includes residues on one face of transmembrane domain three at the extracellular side of the channel that face the center of ChR2. Furthermore, we examined the role of a residue at the extracellular side of transmembrane domain three in ion conductance. We show that ion conductance is mediated, in part, by hydrogen bonding at the extracellular side of transmembrane domain three. These results provide a starting point for examining the cation permeability pathway for ChR2.

Re-introduction of transmembrane serine residues reduce the minimum pore diameter of channelrhodopsin-2.

Channelrhodopsin-2 (ChR2) is a microbial-type rhodopsin found in the green algae Chlamydomonas reinhardtii. Under physiological conditions, ChR2 is an inwardly rectifying cation channel that permeates a wide range of mono- and divalent cations. Although this protein shares a high sequence homology with other microbial-type rhodopsins, which are ion pumps, ChR2 is an ion channel. A sequence alignment of ChR2 with bacteriorhodopsin, a proton pump, reveals that ChR2 lacks specific motifs and residues, such as serine and threonine, known to contribute to non-covalent interactions within transmembrane domains. We hypothesized that reintroduction of the eight transmembrane serine residues present in bacteriorhodopsin, but not in ChR2, will restrict the conformational flexibility and reduce the pore diameter of ChR2. In this work, eight single serine mutations were created at homologous positions in ChR2. Additionally, an endogenous transmembrane serine was replaced with alanine. We measured kinetics, changes in reversal potential, and permeability ratios in different alkali metal solutions using two-electrode voltage clamp. Applying excluded volume theory, we calculated the minimum pore diameter of ChR2 constructs. An analysis of the results from our experiments show that reintroducing serine residues into the transmembrane domain of ChR2 can restrict the minimum pore diameter through inter- and intrahelical hydrogen bonds while the removal of a transmembrane serine results in a larger pore diameter. Therefore, multiple positions along the intracellular side of the transmembrane domains contribute to the cation permeability of ChR2.

The cation selectivity of the ZIP transporters.

Zinc is a required micronutrient for cellular homeostasis and is essential for the structure and/or function of 100s of biological processes. Despite the central importance of zinc in physiology, the mechanism by which this transition metal is transported into cells is not well understood. The first human zinc importer was identified in 2000. Since that time, a new family of proteins, named ZIP, for Zrt-, Irt-like proteins, has been shown to be expressed in a tissue-specific manner in humans. There are a total of 14 members of this family, which can be further divided into four subfamilies based on sequence similarities: ZIPI, ZIPII, gufA and LIV-1. It has been shown that these proteins are expressed on the plasma membrane as well as on intracellular organelles and each function or are proposed to function to increase the cytosolic concentration of zinc. While the subcellular localization of most of the ZIP family of proteins has been elucidated, the mechanistic details of these proteins including the driving force of cation translocation, residues essential for transport and the molecular determinants which define the differing cation selectivity between the ZIP family of proteins have not been well resolved. The objective of this review is to describe the current understanding of cation transport, mediated by the mammalian family of ZIP proteins and to present some molecular determinants that may contribute to the differing substrate specificity for this family of proteins.