General principles of secondary active transporter function

Kinetic Diagram Analysis: A Python Library for Calculating Steady-State Observables of Biochemical Systems Analytically

Kinetic diagrams are commonly used to represent biochemical systems in order to study phenomena such as free energy transduction and ion selectivity. While numerical methods are commonly used to analyze such kinetic networks, the diagram method by King, Altman and Hill makes it possible to construct exact algebraic expressions for steady-state observables in terms of the rate constants of the kinetic diagram. However, manually obtaining these expressions becomes infeasible for models of even modest complexity as the number of the required intermediate diagrams grows with the factorial of the number of states in the diagram. We developed Kinetic Diagram Analysis (KDA), a Python library that programmatically generates the relevant diagrams and expressions from a user-defined kinetic diagram. KDA outputs symbolic expressions for state probabilities and cycle fluxes at steady-state that can be symbolically manipulated and evaluated to quantify macroscopic system observables. We demonstrate the KDA approach for examples drawn from the biophysics of active secondary transmembrane transporters. For a generic 6-state antiporter model, we show how the introduction of a single leakage transition reduces transport efficiency by quantifying substrate turnover. We apply KDA to a real-world example, the 8-state free exchange model of the small multidrug resistance transporter EmrE of Hussey et al. (J. Gen. Physiol., 2020, 152, e201912437), where a change in transporter phenotype is achieved by biasing two different subsets of kinetic rates: alternating access and substrate unbinding rates. KDA is made available as open source software under the GNU General Public License version 3.

Intraventricular haemorrhage in premature infants: the role of immature neuronal salt and water transport

Intraventricular haemorrhage is a common complication of premature birth. Survivors are often left with cerebral palsy, intellectual disability and/or hydrocephalus. Animal models suggest that brain tissue shrinkage, with subsequent vascular stretch and tear, is an important step in the pathophysiology, but the cause of this shrinkage is unknown. Clinical risk factors for intraventricular haemorrhage are biomarkers of hypoxic-ischaemic stress, which causes mature neurons to swell. However, immature neuronal volume might shift in the opposite direction in these conditions. This is because immature neurons express the chloride, salt and water transporter NKCC1, which subserves regulatory volume increases in non-neural cells, whereas mature neurons express KCC2, which subserves regulatory volume decreases. When hypoxic-ischaemic conditions reduce active ion transport and increase the cytoplasmic membrane permeability, the effects of these transporters are diminished. Consequentially, mature neurons swell (cytotoxic oedema), whereas immature neurons might shrink. After hypoxic-ischaemic stress, in vivo and in vitro multi-photon imaging of perinatal transgenic mice demonstrated shrinkage of viable immature neurons, bulk tissue shrinkage and blood vessel displacement. Neuronal shrinkage was correlated with age-dependent membrane salt and water transporter expression using immunohistochemistry. Shrinkage of immature neurons was prevented by prior genetic or pharmacological inhibition of NKCC1 transport. These findings open new avenues of investigation for the detection of acute brain injury by neuroimaging, in addition to prevention of neuronal shrinkage and the ensuing intraventricular haemorrhage, in premature infants.

Kinetic Diagram Analysis: A Python Library for Calculating Steady-State Observables of Biochemical Systems Analytically

Kinetic diagrams are commonly used to represent biochemical systems in order to study phenomena such as free energy transduction and ion selectivity. While numerical methods are commonly used to analyze such kinetic networks, the diagram method by King, Altman and Hill makes it possible to construct exact algebraic expressions for steady-state observables in terms of the rate constants of the kinetic diagram. However, manually obtaining these expressions becomes infeasible for models of even modest complexity as the number of the required intermediate diagrams grows with the factorial of the number of states in the diagram. We developed Kinetic Diagram Analysis (KDA), a Python library that programmatically generates the relevant diagrams and expressions from a user-defined kinetic diagram. KDA outputs symbolic expressions for state probabilities and cycle fluxes at steady-state that can be symbolically manipulated and evaluated to quantify macroscopic system observables. We demonstrate the KDA approach for examples drawn from the biophysics of active secondary transmembrane transporters. For a generic 6-state antiporter model, we show how the introduction of a single leakage transition reduces transport efficiency by quantifying substrate turnover. We apply KDA to a real-world example, the 8-state free exchange model of the small multidrug resistance transporter EmrE of Hussey et al (J General Physiology 152 (2020), e201912437), where a change in transporter phenotype is achieved by biasing two different subsets of kinetic rates: alternating access and substrate unbinding rates. KDA is made available as open source software under the GNU General Public License version 3.

Structure and Function of Auxin Transporters

Auxins, a group of central hormones in plant growth and development, are transported by a diverse range of transporters with distinct biochemical and structural properties. This review summarizes the current knowledge on all known auxin transporters with respect to their biochemical and biophysical properties and the methods used to characterize them. In particular, we focus on the recent advances that were made concerning the PIN-FORMED family of auxin exporters. Insights derived from solving their structures have improved our understanding of the auxin export process, and we discuss the current state of the art on PIN-mediated auxin transport, including the use of biophysical methods to examine their properties. Understanding the mechanisms of auxin transport is crucial for understanding plant growth and development, as well as for the development of more effective strategies for crop production and plant biotechnology.

The Arabidopsis AtSWEET13 transporter discriminates sugars by selective facial and positional substrate recognition

Transporters are targeted by endogenous metabolites and exogenous molecules to reach cellular destinations, but it is generally not understood how different substrate classes exploit the same transporter's mechanism. Any disclosure of plasticity in transporter mechanism when treated with different substrates becomes critical for developing general selectivity principles in membrane transport catalysis. Using extensive molecular dynamics simulations with an enhanced sampling approach, we select the Arabidopsis sugar transporter AtSWEET13 as a model system to identify the basis for glucose versus sucrose molecular recognition and transport. Here we find that AtSWEET13 chemical selectivity originates from a conserved substrate facial selectivity demonstrated when committing alternate access, despite mono-/di-saccharides experiencing differing degrees of conformational and positional freedom throughout other stages of transport. However, substrate interactions with structural hallmarks associated with known functional annotations can help reinforce selective preferences in molecular transport.

19F-NMR Probing of Ion-Induced Conformational Changes in Detergent-Solubilized and Nanodisc-Reconstituted NCX_Mj

Consecutive interactions of 3Na+ or 1Ca2+ with the Na+/Ca2+ exchanger (NCX) result in an alternative exposure (access) of the cytosolic and extracellular vestibules to opposite sides of the membrane, where ion-induced transitions between the outward-facing (OF) and inward-facing (IF) conformational states drive a transport cycle. Here, we investigate sub-state populations of apo and ion-bound species in the OF and IF states by analyzing detergent-solubilized and nanodisc-reconstituted preparations of NCX_Mj with 19F-NMR. The 19F probe was covalently attached to the cysteine residues at entry locations of the cytosolic and extracellular vestibules. Multiple sub-states of apo and ion-bound species were observed in nanodisc-reconstituted (but not in detergent-solubilized) NCX_Mj, meaning that the lipid-membrane environment preconditions multiple sub-state populations toward the OF/IF swapping. Most importantly, ion-induced sub-state redistributions occur within each major (OF or IF) state, where sub-state interconversions may precondition the OF/IF swapping. In contrast with large changes in population redistributions, the sum of sub-state populations within each inherent state (OF or IF) remains nearly unchanged upon ion addition. The present findings allow the further elucidation of structure-dynamic modules underlying ion-induced conformational changes that determine a functional asymmetry of ion access/translocation at opposite sides of the membrane and ion transport rates concurring physiological demands.

Molecular mechanisms of Na+-driven bile acid transport in human NTCP

Human Na+ taurocholate co-transporting protein (hNTCP) is a key bile salt transporter to maintain enterohepatic circulation and is responsible for the recognition of hepatitis B and D viruses. Despite landmark cryoelectron microscopy studies revealing open-pore and inward-facing states of hNTCP stabilized by antibodies, the transport mechanism remains largely unknown. To address this knowledge gap, we used molecular dynamics and enhanced sampling metadynamics simulations to elucidate the intrinsic mechanism of hNTCP-mediated taurocholate acid (TCA) transport driven by Na+ binding. We uncovered three TCA-binding modes, including one that closely matched the limited cryoelectron microscopy density observed in the open-pore hNTCP. We also captured several key hNTCP conformations in the substrate transport cycle, particularly including an outward-facing, substrate-bound state. Furthermore, we provided thermodynamic evidence supporting that changes in the Na+-binding state drive the TCA transport by exploiting the amphiphilic nature of the substrate and modulating the protein environment, thereby enabling the TCA molecule to flip through. Understanding these mechanistic details of Na+-driven bile acid transport may aid in the development of hNTCP-targeted therapies for liver diseases.

Conformational free-energy landscapes of a Na+/Ca2+ exchanger explain its alternating-access mechanism and functional specificity

Secondary-active transporters catalyze the movement of myriad substances across all cellular membranes, typically against opposing concentration gradients, and without consuming any ATP. To do so, these proteins employ an intriguing structural mechanism evolved to be activated only upon recognition or release of the transported species. We examine this self-regulated mechanism using a homolog of the cardiac Na+/Ca2+ exchanger as a model system. Using advanced computer simulations, we map out the complete functional cycle of this transporter, including unknown conformations that we validate against existing experimental data. Calculated free-energy landscapes reveal why this transporter functions as an antiporter rather than a symporter, why it specifically exchanges Na+ and Ca2+, and why the stoichiometry of this exchange is exactly 3:1. We also rationalize why the protein does not exchange H+ for either Ca2+ or Na+, despite being able to bind H+ and its high similarity with H+/Ca2+ exchangers. Interestingly, the nature of this transporter is not explained by its primary structural states, known as inward- and outward-open conformations; instead, the defining factor is the feasibility of conformational intermediates between those states, wherein access pathways leading to the substrate binding sites become simultaneously occluded from both sides of the membrane. This analysis offers a physically coherent, broadly transferable route to understand the emergence of function from structure among secondary-active membrane transporters.

Oligomerization of Monoamine Transporters

Transporters of the monoamine transporter (MAT) family regulate the uptake of important neurotransmitters like dopamine, serotonin, and norepinephrine. The MAT family functions using the electrochemical gradient of ions across the membrane and comprises three transporters, dopamine transporter (DAT), serotonin transporter (SERT), and norepinephrine transporter (NET). MAT transporters have been observed to exist in monomeric states to higher-order oligomeric states. Structural features, allosteric modulation, and lipid environment regulate the oligomerization of MAT transporters. NET and SERT oligomerization are regulated by levels of PIP2 present in the membrane. The kink present in TM12 in the MAT family is crucial for dimer interface formation. Allosteric modulation in the dimer interface hinders dimer formation. Oligomerization also influences the transporters' function, trafficking, and regulation. This chapter will focus on recent studies on monoamine transporters and discuss the factors affecting their oligomerization and its impact on their function.

Mechanism of substrate binding and transport in BASS transporters

The bile acid sodium symporter (BASS) family transports a wide array of molecules across membranes, including bile acids in humans, and small metabolites in plants. These transporters, many of which are sodium-coupled, have been shown to use an elevator mechanism of transport, but exactly how substrate binding is coupled to sodium ion binding and transport is not clear. Here, we solve the crystal structure at 2.3 Å of a transporter from Neisseria meningitidis (ASBTNM) in complex with pantoate, a potential substrate of ASBTNM. The BASS family is characterised by two helices that cross-over in the centre of the protein in an arrangement that is intricately held together by two sodium ions. We observe that the pantoate binds, specifically, between the N-termini of two of the opposing helices in this cross-over region. During molecular dynamics simulations the pantoate remains in this position when sodium ions are present but is more mobile in their absence. Comparison of structures in the presence and absence of pantoate demonstrates that pantoate elicits a conformational change in one of the cross-over helices. This modifies the interface between the two domains that move relative to one another to elicit the elevator mechanism. These results have implications, not only for ASBTNM but for the BASS family as a whole and indeed other transporters that work through the elevator mechanism.

Slc11 Synapomorphy: A Conserved 3D Framework Articulating Carrier Conformation Switch

Transmembrane carriers of the Slc11 family catalyze proton (H+)-dependent uptake of divalent metal ions (Me2+) such as manganese and iron-vital elements coveted during infection. The Slc11 mechanism of high-affinity Me2+ cell import is selective and conserved between prokaryotic (MntH) and eukaryotic (Nramp) homologs, though processes coupling the use of the proton motive force to Me2+ uptake evolved repeatedly. Adding bacterial piracy of Nramp genes spread in distinct environmental niches suggests selective gain of function that may benefit opportunistic pathogens. To better understand Slc11 evolution, Alphafold (AF2)/Colabfold (CF) 3D predictions for bacterial sequences from sister clades of eukaryotic descent (MCb and MCg) were compared using both native and mutant templates. AF2/CF model an array of native MCb intermediates spanning the transition from outwardly open (OO) to inwardly open (IO) carriers. In silico mutagenesis targeting (i) a set of (evolutionarily coupled) sites that may define Slc11 function (putative synapomorphy) and (ii) residues from networked communities evolving during MCb transition indicates that Slc11 synapomorphy primarily instructs a Me2+-selective conformation switch which unlocks carrier inner gate and contributes to Me2+ binding site occlusion and outer gate locking. Inner gate opening apparently proceeds from interaction between transmembrane helix (h) h5, h8 and h1a. MCg1 xenologs revealed marked differences in carrier shape and plasticity, owing partly to an altered intramolecular H+ network. Yet, targeting Slc11 synapomorphy also converted MCg1 IO models to an OO state, apparently mobilizing the same residues to control gates. But MCg1 response to mutagenesis differed, with extensive divergence within this clade correlating with MCb-like modeling properties. Notably, MCg1 divergent epistasis marks the emergence of the genus Bordetella-Achromobacter. Slc11 synapomorphy localizes to the 3D areas that deviate least among MCb and MCg1 models (either IO or OO) implying that it constitutes a 3D network of residues articulating a Me2+-selective carrier conformation switch which is maintained in fast-evolving clades at the cost of divergent epistatic interactions impacting carrier shape and dynamics.

Thermodynamically consistent determination of free energies and rates in kinetic cycle models

Kinetic and thermodynamic models of biological systems are commonly used to connect microscopic features to system function in a bottom-up multiscale approach. The parameters of such models-free energy differences for equilibrium properties and in general rates for equilibrium and out-of-equilibrium observables-have to be measured by different experiments or calculated from multiple computer simulations. All such parameters necessarily come with uncertainties so that when they are naively combined in a full model of the process of interest, they will generally violate fundamental statistical mechanical equalities, namely detailed balance and an equality of forward/backward rate products in cycles due to Hill. If left uncorrected, such models can produce arbitrary outputs that are physically inconsistent. Here, we develop a maximum likelihood approach (named multibind) based on the so-called potential graph to combine kinetic or thermodynamic measurements to yield state-resolved models that are thermodynamically consistent while being most consistent with the provided data and their uncertainties. We demonstrate the approach with two theoretical models, a generic two-proton binding site and a simplified model of a sodium/proton antiporter. We also describe an algorithm to use the multibind approach to solve the inverse problem of determining microscopic quantities from macroscopic measurements and, as an example, we predict the microscopic p K a values and protonation states of a small organic molecule from 1D NMR data. The multibind approach is applicable to any thermodynamic or kinetic model that describes a system as transitions between well-defined states with associated free energy differences or rates between these states. A Python package multibind, which implements the approach described here, is made publicly available under the MIT Open Source license.

Mechanism of substrate binding and transport in BASS transporters

The Bile Acid Sodium Symporter (BASS) family transports a wide array of molecules across membranes, including bile acids in humans, and small metabolites in plants. These transporters, many of which are sodium-coupled, have been shown to use an elevator mechanism of transport, but exactly how substrate binding is coupled to sodium ion binding and transport is not clear. Here we solve the crystal structure at 2.3 Å of a transporter from Neisseria Meningitidis (ASBTNM) in complex with pantoate, a potential substrate of ASBTNM. The BASS family is characterised by two helices that cross-over in the centre of the protein in an arrangement that is intricately held together by two sodium ions. We observe that the pantoate binds, specifically, between the N-termini of two of the opposing helices in this cross-over region. During molecular dynamics simulations the pantoate remains in this position when sodium ions are present but is more mobile in their absence. Comparison of structures in the presence and absence of pantoate demonstrates that pantoate elicits a conformational change in one of the cross-over helices. This modifies the interface between the two domains that move relative to one another to elicit the elevator mechanism. These results have implications, not only for ASBTNM but for the BASS family as a whole and indeed other transporters that work through the elevator mechanism.

Thermodynamically consistent determination of free energies and rates in kinetic cycle models

Kinetic and thermodynamic models of biological systems are commonly used to connect microscopic features to system function in a bottom-up multiscale approach. The parameters of such models-free energy differences for equilibrium properties and in general rates for equilibrium and out-of-equilibrium observables-have to be measured by different experiments or calculated from multiple computer simulations. All such parameters necessarily come with uncertainties so that when they are naively combined in a full model of the process of interest, they will generally violate fundamental statistical mechanical equalities, namely detailed balance and an equality of forward/backward rate products in cycles due to T. Hill. If left uncorrected, such models can produce arbitrary outputs that are physically inconsistent. Here we develop a maximum likelihood approach (named multibind ) based on the so-called potential graph to combine kinetic or thermodynamic measurements to yield state resolved models that are thermodynamically consistent while being most consistent with the provided data and their uncertainties. We demonstrate the approach with two theoretical models, a generic two-proton binding site and a simplified model of a sodium/proton antiporter. We also describe an algorithm to use the multibind approach to solve the inverse problem of determining microscopic quantities from macroscopic measurements and as an example we predict the microscopic p K a s and protonation states of a small organic molecule from 1D NMR data. The multibind approach is applicable to any thermodynamic or kinetic model that describes a system as transitions between well-defined states with associated free energy differences or rates between these states. A Python package multibind , which implements the approach described here, is made publicly available under the MIT Open Source license.

WHY IT MATTERS

The increase in computational efficiency and rapid advances in methodology for quantitative free energy and rate calculations has allowed for the construction of increasingly complex thermodynamic or kinetic "bottom-up" models of chemical and biological processes. These multi-scale models serve as a framework for analyzing aspects of cellular function in terms of microscopic, molecular properties and provide an opportunity to connect molecular mechanisms to cellular function. The underlying model parameters-free energy differences or rates-are constrained by thermodynamic identities over cycles of states but these identities are not necessarily obeyed during model construction, thus potentially leading to inconsistent models. We address these inconsistencies through the use of a maximum likelihood approach for free energies and rates to adjust the model parameters in such a way that they are maximally consistent with the input parameters and exactly fulfill the thermodynamic cycle constraints. This approach enables formulation of thermodynamically consistent multi-scale models from simulated or experimental measurements.

Ins and Outs of Rocker Switch Mechanism in Major Facilitator Superfamily of Transporters

The major facilitator superfamily (MFS) of transporters consists of three classes of membrane transporters: symporters, uniporters, and antiporters. Despite such diverse functions, MFS transporters are believed to undergo similar conformational changes within their distinct transport cycles, known as the rocker-switch mechanism. While the similarities between conformational changes are noteworthy, the differences are also important since they could potentially explain the distinct functions of symporters, uniporters, and antiporters of the MFS superfamily. We reviewed a variety of experimental and computational structural data on a select number of antiporters, symporters, and uniporters from the MFS family to compare the similarities and differences of the conformational dynamics of three different classes of transporters.

The role of conformational change and key glutamic acid residues in the ClC-ec1 antiporter

The triple glutamine (Q) mutant (QQQ) structure of a Cl-/H+ antiporter from Escherichia coli (ClC-ec1) displaying a novel backbone arrangement has been used to challenge the long-held notion that Cl-/H+ antiporters do not operate through large conformational motions. The QQQ mutant substitutes the glutamine residue for an external glutamate E148, an internal glutamate E203, and a third glutamate E113 that hydrogen-bonds with E203. However, it is unknown if QQQ represents a physiologically relevant state, as well as how the protonation of the wild-type glutamates relates to the global dynamics. We herein apply continuous constant-pH molecular dynamics to investigate the H+-coupled dynamics of ClC-ec1. Although any large-scale conformational rearrangement upon acidification would be due to the accumulation of excess charge within the protein, protonation of the glutamates significantly impacts mainly the local structure and dynamics. Despite the fact that the extracellular pore enlarges at acidic pHs, an occluded ClC-ec1 within the active pH range of 3.5-7.5 requires a protonated E148 to facilitate extracellular Cl- release. E203 is also involved in the intracellular H+ transfer as an H+ acceptor. The water wire connection of E148 with the intracellular solution is regulated by the charge states of the E113/E203 dyad with coupled proton titration. However, the dynamics extracted from our simulations are not QQQ-like, indicating that the QQQ mutant does not represent the behavior of the wild-type ClC-ec1. These findings reinforce the necessity of having a protonatable residue at the E203 position in ClC-ec1 and suggest that a higher level of complexity exists for the intracellular H+ transfer in Cl-/H+ antiporters.