Conservation of knotted and slipknotted topology in transmembrane transporters.

The existence of nontrivial topology is well accepted in globular proteins but not in membrane proteins. Our comprehensive topological analysis of the Protein Data Bank structures reveals 18 families of transmembrane proteins with nontrivial topology, showing that they constitute a significant number of membrane proteins. Moreover, we found that they comprise one of the largest groups of secondary active transporters. We classified them based on their knotted fingerprint into four groups: three slipknotted and one knotted. Unexpectedly, we found that the same protein can possess two distinct slipknot motifs that correspond to its outward- and inward-open conformational state. Based on the analysis of structures and knotted fingerprints, we show that slipknot topology is directly involved in the conformational transition and substrate transfer. Therefore, entanglement can be used to classify proteins and to find their structure-function relationship. Furthermore, based on the topological analysis of the transmembrane protein structures predicted by AlphaFold, we identified new potentially slipknotted protein families.

Slipknotted and unknotted monovalent cation-proton antiporters evolved from a common ancestor.

While the slipknot topology in proteins has been known for over a decade, its evolutionary origin is still a mystery. We have identified a previously overlooked slipknot motif in a family of two-domain membrane transporters. Moreover, we found that these proteins are homologous to several families of unknotted membrane proteins. This allows us to directly investigate the evolution of the slipknot motif. Based on our comprehensive analysis of 17 distantly related protein families, we have found that slipknotted and unknotted proteins share a common structural motif. Furthermore, this motif is conserved on the sequential level as well. Our results suggest that, regardless of topology, the proteins we studied evolved from a common unknotted ancestor single domain protein. Our phylogenetic analysis suggests the presence of at least seven parallel evolutionary scenarios that led to the current diversity of proteins in question. The tools we have developed in the process can now be used to investigate the evolution of other repeated-domain proteins.

Dimerisation of the Yeast K+ Translocation Protein Trk1 Depends on the K+ Concentration.

In baker's yeast (Saccharomyces cerevisiae), Trk1, a member of the superfamily of K-transporters (SKT), is the main K+ uptake system under conditions when its concentration in the environment is low. Structurally, Trk1 is made up of four domains, each similar and homologous to a K-channel α subunit. Because most K-channels are proteins containing four channel-building α subunits, Trk1 could be functional as a monomer. However, related SKT proteins TrkH and KtrB were crystallised as dimers, and for Trk1, a tetrameric arrangement has been proposed based on molecular modelling. Here, based on Bimolecular Fluorescence Complementation experiments and single-molecule fluorescence microscopy combined with molecular modelling; we provide evidence that Trk1 can exist in the yeast plasma membrane as a monomer as well as a dimer. The association of monomers to dimers is regulated by the K+ concentration.

Communication between N terminus and loop2 tunes Orai activation.

Ca2+ release-activated Ca2+ (CRAC) channels constitute the major Ca2+ entry pathway into the cell. They are fully reconstituted via intermembrane coupling of the Ca2+-selective Orai channel and the Ca2+-sensing protein STIM1. In addition to the Orai C terminus, the main coupling site for STIM1, the Orai N terminus is indispensable for Orai channel gating. Although the extended transmembrane Orai N-terminal region (Orai1 amino acids 73-91; Orai3 amino acids 48-65) is fully conserved in the Orai1 and Orai3 isoforms, Orai3 tolerates larger N-terminal truncations than Orai1 in retaining store-operated activation. In an attempt to uncover the reason for these isoform-specific structural requirements, we analyzed a series of Orai mutants and chimeras. We discovered that it was not the N termini, but the loop2 regions connecting TM2 and TM3 of Orai1 and Orai3 that featured distinct properties, which explained the different, isoform-specific behavior of Orai N-truncation mutants. Atomic force microscopy studies and MD simulations suggested that the remaining N-terminal portion in the non-functional Orai1 N-truncation mutants formed new, inhibitory interactions with the Orai1-loop2 regions, but not with Orai3-loop2. Such a loop2 swap restored activation of the N-truncation Orai1 mutants. To mimic interactions between the N terminus and loop2 in full-length Orai1 channels, we induced close proximity of the N terminus and loop2 via cysteine cross-linking, which actually caused significant inhibition of STIM1-mediated Orai currents. In aggregate, maintenance of Orai activation required not only the conserved N-terminal region but also permissive communication of the Orai N terminus and loop2 in an isoform-specific manner.

A refined atomic scale model of the Saccharomyces cerevisiae K+-translocation protein Trk1p combined with experimental evidence confirms the role of selectivity filter glycines and other key residues.

Potassium ion (K+) uptake in yeast is mediated mainly by the Trk1/2 proteins that enable cells to survive on external K+ concentration as low as a few µM. Fungal Trks are related to prokaryotic TRK and Ktr and plant HKT K+ transport systems. Overall sequence similarity is very low, thus requiring experimental verification of homology models. Here a refined structural model of the Saccharomyces cerevisiae Trk1 is presented that was obtained by combining homology modeling, molecular dynamics simulation and experimental verification through functional analysis of mutants. Structural models and experimental results showed that glycines within the selectivity filter, conserved among the K-channel/transporter family, are not only important for protein function, but are also required for correct folding/membrane targeting. A conserved aspartic acid in the PA helix (D79) and a lysine in the M2D helix (K1147) were proposed earlier to interact. Our results suggested individual roles of these residues in folding, structural integrity and function. While mutations of D79 completely abolished protein folding, mutations at position 1147 were tolerated to some extent. Intriguingly, a secondary interaction of D79 with R76 could enhance folding/stability of Trk1 and enable a fraction of Trk1[K1147A] to fold. The part of the ion permeation path containing the selectivity filter is shaped similar to that of ion channels. However below the selectivity filter it is obstructed or regulated by a proline containing loop. The presented model could provide the structural basis for addressing the long standing question if Trk1 is a passive or active ion-translocation system.

MPM motifs of the yeast SKT protein Trk1 can assemble to form a functional K+-translocation system.

The yeast Trk1 polypeptide, like other members of the Superfamily of K Transporters (SKT proteins) consists of four Membrane-Pore-Membrane motifs (MPMs A-D) each of which is homologous to a single K-channel subunit. SKT proteins are thought to have evolved from ancestral K-channels via two gene duplications and thus single MPMs might be able to assemble when located on different polypeptides. To test this hypothesis experimentally we generated a set of partial gene deletions to create alleles encoding one, two, or three MPMs, and analysed the cellular localisation and interactions of these Trk1 fragments using GFP tags and Bimolecular Fluorescence Complementation (BiFC). The function of these partial Trk1 proteins either alone or in combinations was assessed by expressing the encoding genes in a K+-uptake deficient strain lacking also the K-channel Tok1 (trk1,trk2,tok1Δ) and (i) analysing their ability to promote growth in low [K+] media and (ii) by ion flux measurements using "microelectrode based ion flux estimation" (MIFE). We found that proteins containing only one or two MPM motifs can interact with each other and assemble with a polypeptide consisting of the rest of the Trk system to form a functional K+-translocation system.

The S. cerevisiae cation translocation protein Trk1 is functional without its "long hydrophilic loop" but LHL regulates cation translocation activity and selectivity.

In Saccharomyces cerevisiae, K+-uptake under K+-limiting conditions is largely mediated by the cation translocation systems Trk1 and Trk2 belonging to the family of SKT proteins. They are related to two-transmembrane-domain (inward rectifying K-) channels but unlike the symmetrical tetrameric structure of K-channels, a single Trk contains four pore-forming domains (A-D) encoded on one polypeptide chain. Between domains A and B Trks contain large cytosolic regions dubbed "long hydrophilic loop" (LHL). LHLs are not homologous/similar to any other identified protein (domain) and also show little similarity between Trk1 and Trk2. Here we demonstrate that Trk1 is functional without LHL. However, in growth experiments NaCl sensitivity of Trk1[ΔLHL] expressing cells is increased under K+-limiting conditions compared to full-length Trk1. Non-invasive ion flux measurements showed that K+-influx through Trk1 and Trk1[ΔLHL] is decreased in the presence of surplus Na+ due to permeability of the proteins for both cations and competition between them. Trk1[ΔLHL] is less affected than full-length Trk1 because it is more selective for K+ over Na+. Furthermore, K+ re-uptake after starvation is delayed and decreased in Trk1[ΔLHL]. Thus, a role of LHL is regulating cation fluxes through Trk1 by (i) allowing also Na+ to pass if monovalent cations (mainly K+) are limiting and (ii) by accelerating/enhancing a switch from low to high affinity ion translocation. We propose that LHL could modulate Trk1 transport properties via direct influence on a transmembrane helix (M2A) which can switch between bent and straight conformation, thereby directly modifying the radius of the pore and selectivity filter.

Transmembrane helix connectivity in Orai1 controls two gates for calcium-dependent transcription.

The channel Orai1 requires Ca2+ store depletion in the endoplasmic reticulum and an interaction with the Ca2+ sensor STIM1 to mediate Ca2+ signaling. Alterations in Orai1-mediated Ca2+ influx have been linked to several pathological conditions including immunodeficiency, tubular myopathy, and cancer. We screened large-scale cancer genomics data sets for dysfunctional Orai1 mutants. Five of the identified Orai1 mutations resulted in constitutively active gating and transcriptional activation. Our analysis showed that certain Orai1 mutations were clustered in the transmembrane 2 helix surrounding the pore, which is a trigger site for Orai1 channel gating. Analysis of the constitutively open Orai1 mutant channels revealed two fundamental gates that enabled Ca2+ influx: Arginine side chains were displaced so they no longer blocked the pore, and a chain of water molecules formed in the hydrophobic pore region. Together, these results enabled us to identify a cluster of Orai1 mutations that trigger Ca2+ permeation associated with gene transcription and provide a gating mechanism for Orai1.

Cholesterol modulates Orai1 channel function.

STIM1 (stromal interaction molecule 1) and Orai proteins are the essential components of Ca(2+) release-activated Ca(2+) (CRAC) channels. We focused on the role of cholesterol in the regulation of STIM1-mediated Orai1 currents. Chemically induced cholesterol depletion enhanced store-operated Ca(2+) entry (SOCE) and Orai1 currents. Furthermore, cholesterol depletion in mucosal-type mast cells augmented endogenous CRAC currents, which were associated with increased degranulation, a process that requires calcium influx. Single point mutations in the Orai1 amino terminus that would be expected to abolish cholesterol binding enhanced SOCE to a similar extent as did cholesterol depletion. The increase in Orai1 activity in cells expressing these cholesterol-binding-deficient mutants occurred without affecting the amount in the plasma membrane or the coupling of STIM1 to Orai1. We detected cholesterol binding to an Orai1 amino-terminal fragment in vitro and to full-length Orai1 in cells. Thus, our data showed that Orai1 senses the amount of cholesterol in the plasma membrane and that the interaction of Orai1 with cholesterol inhibits its activity, thereby limiting SOCE.

A calcium-accumulating region, CAR, in the channel Orai1 enhances Ca(2+) permeation and SOCE-induced gene transcription.

The Ca(2+) release-activated Ca(2+) channel mediates Ca(2+) influx in a plethora of cell types, thereby controlling diverse cellular functions. The channel complex is composed of stromal interaction molecule 1 (STIM1), an endoplasmic reticulum Ca(2+)-sensing protein, and Orai1, a plasma membrane Ca(2+) channel. Channels composed of STIM1 and Orai1 mediate Ca(2+) influx even at low extracellular Ca(2+) concentrations. We investigated whether the activity of Orai1 adapted to different environmental Ca(2+) concentrations. We used homology modeling and molecular dynamics simulations to predict the presence of an extracellular Ca(2+)-accumulating region (CAR) at the pore entrance of Orai1. Furthermore, simulations of Orai1 proteins with mutations in CAR, along with live-cell experiments, or simulations and electrophysiological recordings of the channel with transient, electrostatic loop3 interacting with loop1 (the site of CAR) determined that CAR enhanced Ca(2+) permeation most efficiently at low external Ca(2+) concentrations. Consistent with these results, cells expressing Orai1 CAR mutants exhibited impaired gene expression stimulated by the Ca(2+)-activated transcription factor nuclear factor of activated T cells (NFAT). We propose that the Orai1 channel architecture with a close proximity of CAR to the selectivity filter, which enables Ca(2+)-selective ion permeation, enhances the local extracellular Ca(2+) concentration to maintain Ca(2+)-dependent gene regulation even in environments with relatively low Ca(2+)concentrations.

Binding-competent states for L-arginine in E. coli arginine repressor apoprotein.

Arginine repressor of E. coli is a multifunctional hexameric protein that provides feedback regulation of arginine metabolism upon activation by the negatively cooperative binding of L-arginine. Interpretation of this complex system requires an understanding of the protein's conformational landscape. The ~50 kDa hexameric C-terminal domain was studied by 100 ns molecular dynamics simulations in the presence and absence of the six L-arg ligands that bind at the trimer-trimer interface. A rotational shift between trimers followed by rotational oscillation occurs in the production phase of the simulations only when L-arg is absent. Analysis of the system reveals that the degree of rotation is correlated with the number of hydrogen bonds across the trimer interface. The trajectory presents frames with one or more apparently open binding sites into which one L-arg could be docked successfully in three different instances, indicating that a binding-competent state of the system is occasionally sampled. Simulations of the resulting singly-liganded systems reveal for the first time that the binding of one L-arg results in a holoprotein-like conformational distribution.

Regulation of the transient receptor potential channel TRPA1 by its N-terminal ankyrin repeat domain.

The transient receptor potential channel A1 (TRPA1) is unique among ion channels of higher vertebrates in that it harbors a large ankyrin repeat domain. The TRPA1 channel is expressed in the inner ear and in nociceptive neurons. It is involved in hearing as well as in the perception of pungent and irritant chemicals. The ankyrin repeat domain has special mechanical properties, which allows it to function as a soft spring that can be extended over a large range while maintaining structural integrity. A calcium-binding site has been experimentally identified within the ankyrin repeats. We built a model of the N-terminal 17 ankyrin repeat structure, including the calcium-binding EF-hand. In our simulations we find the calcium-bound state to be rigid as compared to the calcium-free state. While the end-to-end distance can change by almost 50% in the apo form, these fluctuations are strongly reduced by calcium binding. This increase in stiffness that constraints the end-to-end distance in the holo form is predicted to affect the force acting on the gate of the TRPA1 channel, thereby changing its open probability. Simulations of the transmembrane domain of TRPA1 show that residue N855, which has been associated with familial episodic pain syndrome, forms a strong link between the S4-S5 connecting helix and S1, thereby creating a direct force link between the N-terminus and the gate. The N855S mutation weakens this interaction, thereby reducing the communication between the N-terminus and the transmembrane part of TRPA1.

Biphasic kinetic behavior of E. coli WrbA, an FMN-dependent NAD(P)H:quinone oxidoreductase.

The E. coli protein WrbA is an FMN-dependent NAD(P)H:quinone oxidoreductase that has been implicated in oxidative defense. Three subunits of the tetrameric enzyme contribute to each of four identical, cavernous active sites that appear to accommodate NAD(P)H or various quinones, but not simultaneously, suggesting an obligate tetramer with a ping-pong mechanism in which NAD departs before oxidized quinone binds. The present work was undertaken to evaluate these suggestions and to characterize the kinetic behavior of WrbA. Steady-state kinetics results reveal that WrbA conforms to a ping-pong mechanism with respect to the constancy of the apparent Vmax to Km ratio with substrate concentration. However, the competitive/non-competitive patterns of product inhibition, though consistent with the general class of bi-substrate reactions, do not exclude a minor contribution from additional forms of the enzyme. NMR results support the presence of additional enzyme forms. Docking and energy calculations find that electron-transfer-competent binding sites for NADH and benzoquinone present severe steric overlap, consistent with the ping-pong mechanism. Unexpectedly, plots of initial velocity as a function of either NADH or benzoquinone concentration present one or two Michaelis-Menten phases depending on the temperature at which the enzyme is held prior to assay. The effect of temperature is reversible, suggesting an intramolecular conformational process. WrbA shares these and other details of its kinetic behavior with mammalian DT-diaphorase, an FAD-dependent NAD(P)H:quinone oxidoreductase. An extensive literature review reveals several other enzymes with two-plateau kinetic plots, but in no case has a molecular explanation been elucidated. Preliminary sedimentation velocity analysis of WrbA indicates a large shift in size of the multimer with temperature, suggesting that subunit assembly coupled to substrate binding may underlie the two-plateau behavior. An additional aim of this report is to bring under wider attention the apparently widespread phenomenon of two-plateau Michaelis-Menten plots.