Energetics and mechanism of anion permeation across formate-nitrite transporters

Single-cell transcriptomics across 2,534 microbial species reveals functional heterogeneity in the rumen microbiome

Deciphering the activity of individual microbes within complex communities and environments remains a challenge. Here we describe the development of microbiome single-cell transcriptomics using droplet-based single-cell RNA sequencing and pangenome-based computational analysis to characterize the functional heterogeneity of the rumen microbiome. We generated a microbial genome database (the Bovine Gastro Microbial Genome Map) as a functional reference map for the construction of a single-cell transcriptomic atlas of the rumen microbiome. The atlas includes 174,531 microbial cells and 2,534 species, of which 172 are core active species grouped into 12 functional clusters. We detected single-cell-level functional roles, including a key role for Basfia succiniciproducens in the carbohydrate metabolic niche of the rumen microbiome. Furthermore, we explored functional heterogeneity and reveal metabolic niche trajectories driven by biofilm formation pathway genes within B. succiniciproducens. Our results provide a resource for studying the rumen microbiome and illustrate the diverse functions of individual microbial cells that drive their ecological niche stability or adaptation within the ecosystem.

The formate-hydrogen axis and its impact on the physiology of enterobacterial fermentation

Formic acid (HCOOH) and dihydrogen (H2) are characteristic products of enterobacterial mixed-acid fermentation, with H2 generation increasing in conjunction with a decrease in extracellular pH. Formate and acetyl-CoA are generated by radical-based and coenzyme A-dependent cleavage of pyruvate catalysed by pyruvate formate-lyase (PflB). Formate is also the source of H2, which is generated along with carbon dioxide through the action of the membrane-associated, cytoplasmically-oriented formate hydrogenlyase (FHL-1) complex. Synthesis of the FHL-1 complex is completely dependent on the cytoplasmic accumulation of formate. Consequently, formate determines its own disproportionation into H2 and CO2 by the FHL-1 complex. Cytoplasmic formate levels are controlled by FocA, a pentameric channel that translocates formic acid/formate bidirectionally between the cytoplasm and periplasm. Each protomer of FocA has a narrow hydrophobic pore through which neutral formic acid can pass. Two conserved amino acid residues, a histidine and a threonine, at the center of the pore control directionality of translocation. The histidine residue is essential for pH-dependent influx of formic acid. Studies with the formate analogue hypophosphite and amino acid variants of FocA suggest that the mechanisms of formic acid efflux and influx differ. Indeed, current data suggest, depending on extracellular formate levels, two separate uptake mechanisms exist, both likely contributing to maintain pH homeostasis. Bidirectional formate/formic acid translocation is dependent on PflB and influx requires an active FHL-1 complex. This review describes the coupling of formate and H2 production in enterobacteria.

Evidence for bidirectional formic acid translocation in vivo via the Escherichia coli formate channel FocA

Pentameric FocA permeates either formate or formic acid bidirectionally across the cytoplasmic membrane of anaerobically growing Escherichia coli. Each protomer of FocA has its own hydrophobic pore, but it is unclear whether formate or neutral formic acid is translocated in vivo. Here, we measured total and dicyclohexylcarbodiimide (DCCD)-inhibited proton flux out of resting, fermentatively grown, stationary-phase E. coli cells in dependence on FocA. Using a wild-type strain synthesizing native FocA, it was shown that using glucose as a source of formate, DCCD-independent proton efflux was ∼2.5 mmol min-1, while a mutant lacking FocA showed only DCCD-inhibited, FOF1-ATPase-dependent proton-efflux. A strain synthesizing a chromosomally-encoded FocAH209N variant that functions exclusively to translocate formic acid out of the cell, showed a further 20 % increase in FocA-dependent proton efflux relative to the parental strain. Cells synthesizing a FocAT91A variant, which is unable to translocate formic acid out of the cell, showed only DCCD-inhibited proton efflux. When exogenous formate was added, formic acid uptake was shown to be both FocA- and proton motive force-dependent. By measuring rates of H2 production, potassium ion flux and ATPase activity, these data support a role for coupling between formate, proton and K+ ion translocation in maintaining pH and ion gradient homeostasis during fermentation. FocA thus plays a key role in maintaining this homeostatic balance in fermenting cells by bidirectionally translocating formic acid.

Investigation of molecular details of a bacterial cationic amino acid transporter (GkApcT) during arginine transportation using molecular dynamics simulation and umbrella sampling techniques

CONTEXT

Cationic amino acid transporters (CATs) facilitate arginine transport across membranes and maintain its levels in various tissues and organs, but their overexpression has been associated with severe cancers. A recent study identified the alternating access mechanism and critical residues involved in arginine transportation in a cationic amino acid transporter from Geobacillus kaustophilus (GkApcT). Here, we used molecular dynamics (MD) simulation methods to investigate the transportation mechanism of arginine (Arg) through GkApcT. The results revealed that arginine strongly interacts with specific binding site residues (Thr43, Asp111, Glu115, Lys191, Phe231, Ile234, and Asp237). Based on the umbrella sampling, the main driving force for arginine transport is the polar interactions of the arginine with channel-lining residues. An in-depth description of the dissociation mechanism and binding energy analysis brings valuable insight into the interactions between arginine and transporter residues, facilitating the design of effective CAT inhibitors in cancer cells.

METHODS

The membrane-protein system was constructed by uploading the prokaryotic CAT (PDB ID: 6F34) to the CHARMM-GUI web server. Molecular dynamics simulations were done using the GROMACS package, version 5.1.4, with the CHARMM36 force field and TIP3P water model. The MM-PBSA approach was performed for determining the arginine binding free energy. Furthermore, the hotspot residues were identified through per-residue decomposition analysis. The characteristics of the channel such as bottleneck radius and channel length were analyzed using the CaverWeb 1.1 web server. The proton wire inside the transporter was investigated based on the classic Grotthuss mechanism. We also investigated the atomistic details of arginine transportation using the path-based free energy umbrella sampling technique (US).

Insights into the physiological and genomic characterization of three bacterial isolates from a highly alkaline, terrestrial serpentinizing system

The terrestrial serpentinite-hosted ecosystem known as "The Cedars" is home to a diverse microbial community persisting under highly alkaline (pH ~ 12) and reducing (Eh < -550 mV) conditions. This extreme environment presents particular difficulties for microbial life, and efforts to isolate microorganisms from The Cedars over the past decade have remained challenging. Herein, we report the initial physiological assessment and/or full genomic characterization of three isolates: Paenibacillus sp. Cedars ('Paeni-Cedars'), Alishewanella sp. BS5-314 ('Ali-BS5-314'), and Anaerobacillus sp. CMMVII ('Anaero-CMMVII'). Paeni-Cedars is a Gram-positive, rod-shaped, mesophilic facultative anaerobe that grows between pH 7-10 (minimum pH tested was 7), temperatures 20-40°C, and 0-3% NaCl concentration. The addition of 10-20 mM CaCl₂ enhanced growth, and iron reduction was observed in the following order, 2-line ferrihydrite > magnetite > serpentinite ~ chromite ~ hematite. Genome analysis identified genes for flavin-mediated iron reduction and synthesis of a bacillibactin-like, catechol-type siderophore. Ali-BS5-314 is a Gram-negative, rod-shaped, mesophilic, facultative anaerobic alkaliphile that grows between pH 10-12 and temperatures 10-40°C, with limited growth observed 1-5% NaCl. Nitrate is used as a terminal electron acceptor under anaerobic conditions, which was corroborated by genome analysis. The Ali-BS5-314 genome also includes genes for benzoate-like compound metabolism. Anaero-CMMVII remained difficult to cultivate for physiological studies; however, growth was observed between pH 9-12, with the addition of 0.01-1% yeast extract. Anaero-CMMVII is a probable oxygen-tolerant anaerobic alkaliphile with hydrogenotrophic respiration coupled with nitrate reduction, as determined by genome analysis. Based on single-copy genes, ANI, AAI and dDDH analyses, Paeni-Cedars and Ali-BS5-314 are related to other species (P. glucanolyticus and A. aestuarii, respectively), and Anaero-CMMVII represents a new species. The characterization of these three isolates demonstrate the range of ecophysiological adaptations and metabolisms present in serpentinite-hosted ecosystems, including mineral reduction, alkaliphily, and siderophore production.

FocA and its central role in fine-tuning pH homeostasis of enterobacterial formate metabolism

During enterobacterial mixed-acid fermentation, formate is generated from pyruvate by the glycyl-radical enzyme pyruvate formate-lyase (PflB). In Escherichia coli, especially at low pH, formate is then disproportionated to CO₂ and H₂ by the cytoplasmically oriented, membrane-associated formate hydrogenlyase (FHL) complex. If electron acceptors are available, however, formate is oxidized by periplasmically oriented, respiratory formate dehydrogenases. Formate translocation across the cytoplasmic membrane is controlled by the formate channel, FocA, a member of the formate-nitrite transporter (FNT) family of homopentameric anion channels. This review highlights recent advances in our understanding of how FocA helps to maintain intracellular formate and pH homeostasis during fermentation. Efflux and influx of formate/formic acid are distinct processes performed by FocA and both are controlled through protein interaction between FocA's N-terminal domain with PflB. Formic acid efflux by FocA helps to maintain cytoplasmic pH balance during exponential-phase growth. Uptake of formate against the electrochemical gradient (inside negative) is energetically and mechanistically challenging for a fermenting bacterium unless coupled with proton/cation symport. Translocation of formate/formic acid into the cytoplasm necessitates an active FHL complex, whose synthesis also depends on formate. Thus, FocA, FHL and PflB function together to govern formate homeostasis. We explain how FocA achieves efflux of formic acid and propose mechanisms for pH-dependent uptake of formate both with and without proton symport. We propose that FocA displays both channel- and transporter-like behaviour. Whether this translocation behaviour is shared by other members of the FNT family is also discussed.

Exploring the Energy Landscape of Riboswitches Using Collective Variables Based on Tertiary Contacts

Messenger RNA regulatory elements, such as riboswitches, can display a high degree of flexibility. By characterizing their energy landscapes, and corresponding distributions of 3D configurations, structure-function relationships can be elucidated. Molecular dynamics simulation with enhanced sampling is an important strategy used to computationally access free energy landscapes characterizing the accessible 3D conformations of RNAs. While tertiary contacts are thought to play important roles in RNA dynamics, it is difficult, in explicit solvent, to sample the formation and breakage of tertiary contacts, such as helix-helix interactions, pseudoknot interactions, and junction interactions, while maintaining intact secondary structure elements. To this end, we extend previously developed collective variables and metadynamics efforts, to establish a simple metadynamics protocol, which utilizes only one collective variable, based on multiple tertiary contacts, to characterize the underlying free energy landscape of any RNA molecule. We develop a modified collective variable, the tertiary contacts distance (Q_TC), which can probe the formation and breakage of all or selectively chosen tertiary contacts of the RNA. The SAM-I riboswitch in the presence of three ionic and substrate conditions was investigated and validated against the structure ensemble previously generated using SAXS experiments. This efficient and easy to implement all-atom MD simulation based approach incorporating metadynamics to study RNA conformational dynamics can also be transferred to any other type of biomolecule.

Distinguishing functional from structural roles of conserved pore residues during formate translocation by the FocA anion channel

The formate‐specific anion channel FocA of Escherichia coli belongs to the superfamily of homopentameric formate‐nitrite transporters (FNT). Minimally nine amino acid residues are conserved in the formate translocation pore of each protomer of the pentamer, including a histidine (H209) and a threonine (T91), both of which are crucial for bidirectional formate translocation through the pore. Information regarding in vivo functional or structural roles for the other seven conserved residues is limited, or nonexistent. Here, we conducted an amino acid‐exchange analysis of these seven conserved residues. Using an established formate‐responsive lacZ‐based assay to monitor changes in intracellular formate levels and anaerobic growth rate due to the inhibitory formate analog hypophosphite, we identified five of the seven residues analyzed to be important for the structural integrity of the pentamer, in particular, two highly conserved asparagine residues, N213 and N262. The remaining two conserved residues, K156 and N172, were essential for formate/hypophosphite translocation. K156 is located on the periplasmic fringe of the pore and aids the attraction of formate to the channel. Here, we show that this residue is also important for formate efflux from the cytoplasm to the periplasm, suggesting a role in formate release from the pore. N172 could be replaced by alanine with retention of low‐level bidirectional anion translocation function; however, exchange for threonine abolished anion translocation. N172 is, therefore, crucial for bidirectional formate translocation, possibly through its interaction with the conserved pore residue, T91.

Interplay between the Conserved Pore Residues Thr-91 and His-209 Controls Formate Translocation through the FocA Channel

The formate channel A (FocA) belongs to the formate-nitrite transporter (FNT) family, members of which permeate small monovalent anions. FocA from Escherichia coli translocates formate/formic acid bi-directionally across the cytoplasmic membrane during fermentative growth. Two residues are particularly well conserved within the translocation pores of FNTs: threonine-91 and histidine-209, based on E. coli FocA numbering. These residues are located at the tips of two broken transmembrane helices and control anion passage. H209 is the only charged residue within the pore and interacts with T91. Here, we addressed the role of the T91-H209 interaction network in the permeation of formate in vivo through FocA by performing an extensive amino acid-exchange study. Monitoring changes in intracellular formate using a formate-responsive fdhFP::lacZ reporter system revealed that T91 is essential for the ability of FocA to translocate formate bi-directionally. Only exchange for serine was partially tolerated, indicating that the hydroxyl group of T91 is mechanistically important. Substitution of H209 with N or Q was previously shown to convert FocA into a formate efflux channel. We show here that residue-exchanges A, I and T at this position resulted in a similar phenotype. Moreover, efflux function was confirmed for these FocA variants by measuring excreted formate in the culture medium. Substitution of bulky or charged residues for H209 prevented bi-directional formate passage. Studies using hypophosphite, a toxic analogue of formate taken up by FocA, and which causes impaired growth, confirmed that T91 and H209 substitutions essentially abolished, or drastically reduced, FocA´s translocation activity, as shown by reduction of growth rate. The exceptions were T91S- and T91Y-exchange variants that retained partial ability to take up inhibitory hypophosphite. Together, our findings indicate that T91 is essential for formate permeation in both directions; however, it is particularly important to allow anion efflux. Moreover, H209 is essential for formate uptake by FocA, strongly suggesting that protonation-deprotonation of this residue plays a role in formate uptake. Finally, our results substantiate the premise that efflux and influx of formate by FocA are mechanistically distinct processes that are controlled by the interplay between T91 and H209.

The FocA channel functions to maintain intracellular formate homeostasis during Escherichia coli fermentation

FocA translocates formate/formic acid bi-directionally across the cytoplasmic membrane when Escherichia coli grows by fermentation. It remains unclear, however, what physiological benefit is imparted by FocA, because formic acid (pK _a=3.75) can diffuse passively across the membrane, especially at low pH. Here, we monitored changes in intra- and extracellular formate levels during batch-culture fermentation, comparing a parental E. coli K-12 strain with its isogenic focA mutant. Our results show that, regardless of the initial pH in the culture, FocA functions to maintain relatively constant intracellular formate levels during growth. Analysis of a strain synthesizing a FocA_T91A variant with an exchange in a conserved threonine residue within the translocation pore revealed the strain accumulated formate intracellularly and imported formate poorly, but in a pH-dependent manner, which was different to uptake by native FocA. We conclude that FocA maintains formate homeostasis, using different mechanisms for efflux and uptake of the anion.

Hydrogen production in the presence of oxygen by Escherichia coli K-12

Escherichia coli is a facultative anaerobe that can grow in a variety of environmental conditions. In the complete absence of O₂, E. coli can perform a mixed-acid fermentation that contains within it an elaborate metabolism of formic acid. In this study, we use cavity-enhanced Raman spectroscopy (CERS), FTIR, liquid Raman spectroscopy, isotopic labelling and molecular genetics to make advances in the understanding of bacterial formate and H₂ metabolism. It is shown that, under anaerobic (anoxic) conditions, formic acid is generated endogenously, excreted briefly from the cell, and then taken up again to be disproportionated to H₂ and CO₂ by formate hydrogenlyase (FHL-1). However, exogenously added D-labelled formate behaves quite differently from the endogenous formate and is taken up immediately, independently, and possibly by a different mechanism, by the cell and converted to H₂ and CO₂. Our data support an anion-proton symport model for formic acid transport. In addition, when E. coli was grown in a micro-aerobic (micro-oxic) environment it was possible to analyse aspects of formate and O₂ respiration occurring alongside anaerobic metabolism. While cells growing under micro-aerobic conditions generated endogenous formic acid, no H₂ was produced. However, addition of exogenous formate at the outset of cell growth did induce FHL-1 biosynthesis and resulted in formate-dependent H₂ production in the presence of O₂.

Permeation of Fosfomycin through the Phosphate-Specific Channels OprP and OprO of Pseudomonas aeruginosa

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen responsible for many nosocomial infections. It is quite resistant to various antibiotics, caused by the absence of general diffusion pores in the outer membrane. Instead, it contains many substrate-specific channels. Among them are the two phosphate- and pyrophosphate-specific porins OprP and OprO. Phosphonic acid antibiotics such as fosfomycin and fosmidomycin seem to be good candidates for using these channels to enter P. aeruginosa bacteria. Here, we investigated the permeation of fosfomycin through OprP and OprO using electrophysiology and molecular dynamics (MD) simulations. The results were compared to those of the fosmidomycin translocation, for which additional MD simulations were performed. In the electrophysiological approach, we noticed a higher binding affinity of fosfomycin than of fosmidomycin to OprP and OprO. In MD simulations, the ladder of arginine residues and the cluster of lysine residues play an important role in the permeation of fosfomycin through the OprP and OprO channels. Molecular details on the permeation of fosfomycin through OprP and OprO channels were derived from MD simulations and compared to those of fosmidomycin translocation. In summary, this study demonstrates that the selectivity of membrane channels can be employed to improve the permeation of antibiotics into Gram-negative bacteria and especially into resistant P. aeruginosa strains.

A single amino acid exchange converts FocA into a unidirectional efflux channel for formate

During mixed-acid fermentation, Escherichia coli initially translocates formate out of the cell, but re-imports it at lower pH. This is performed by FocA, the archetype of the formate-nitrite transporter (FNT) family of pentameric anion channels. Each protomer of FocA has a hydrophobic pore through which formate/formic acid is bidirectionally translocated. It is not understood how the direction of formate/formic acid passage through FocA is controlled by pH. A conserved histidine residue (H209) is located within the translocation pore, suggesting that protonation/deprotonation might be linked to the direction of formate translocation. Using a formate-responsive lacZ-based reporter system we monitored changes in formate levels in vivo when H209 in FocA was exchanged for either of the non-protonatable amino acids asparagine or glutamine, which occur naturally in some FNTs. These FocA variants (with N or Q) functioned as highly efficient formate efflux channels and the bacteria could neither accumulate formate nor produce hydrogen gas. Therefore, the data in this study suggest that this central histidine residue within the FocA pore is required for pH-dependent formate uptake into E. coli cells. We also address why H209 is evolutionarily conserved and provide a physiological rationale for the natural occurrence of N/Q variants of FNT channels.

Mutational Widening of Constrictions in a Formate-Nitrite/H+ Transporter Enables Aquaporin-Like Water Permeability and Proton Conductance

The unrelated protein families of the microbial formate–nitrite transporters (FNTs) and aquaporins (AQP) likely adapted the same protein fold through convergent evolution. FNTs facilitate weak acid anion/H⁺ cotransport, whereas AQP water channels strictly exclude charged substrates including protons. The FNT channel–like transduction pathway bears two lipophilic constriction sites that sandwich a highly conserved histidine residue. Because of lacking experiments, the function of these constrictions is unclear, and the protonation status of the central histidine during substrate transport remains a matter of debate. Here, we introduced constriction-widening mutations into the prototypical FNT from Escherichia coli, FocA, and assayed formate/H⁺ transport properties, water/solute permeability, and proton conductance. We found that enlargement of these constrictions concomitantly decreased formate/formic acid transport. In contrast to wildtype FocA, the mutants were unable to make use of a transmembrane proton gradient as a driving force. A construct in which both constrictions were eliminated exhibited water permeability, similar to AQPs, although accompanied by a proton conductance. Our data indicate that the lipophilic constrictions mainly act as barriers to isolate the central histidine from the aqueous bulk preventing protonation via proton wires. These results are supportive of an FNT transport model in which the central histidine is uncharged, and weak acid substrate anion protonation occurs in the vestibule regions of the transporter before passing the constrictions.

Improved Sampling and Free Energy Estimates for Antibiotic Permeation through Bacterial Porins

Antibiotics enter into bacterial cells via protein channels that serve as low-energy pathways through the outer membrane, which is otherwise impenetrable. Insights into the molecular mechanisms underlying the transport processes are vital for the development of effective antibacterials. A much-desired prerequisite is an accurate and reproducible determination of free energy surfaces for antibiotic translocation, enabling quantitative and meaningful comparisons of permeation mechanisms for different classes of antibiotics. Inefficient sampling along the orthogonal degrees of freedom, for example, in umbrella sampling and metadynamics approaches, is however a key limitation affecting the accuracy and the convergence of free energy estimates. To overcome this limitation, two sampling methods have been employed in the present study that, respectively, combine umbrella sampling and metadynamics-style biasing schemes with temperature acceleration for improved sampling along orthogonal degrees of freedom. As a model for the transport of bulky solutes, the ciprofloxacin-OmpF system has been selected. The well-tempered metadynamics approach with multiple walkers is compared with its "temperature-accelerated" variant in terms of improvements in sampling and convergence of free energy estimates. We find that the inclusion of collective variables governing solute degrees of freedom and solute-water interactions within the sampling scheme largely alleviates sampling issues. Concerning improved sampling and convergence of free energy estimates from independent simulations, the temperature-accelerated sliced sampling approach that combines umbrella sampling with temperature-accelerated molecular dynamics performs even better as shown for the ciprofloxacin-OmpF system.

Millisecond-Long Simulations of Antibiotics Transport through Outer Membrane Channels

To reach their target site inside Gram-negative bacteria, almost all antibiotics need to cross the outer membrane. Computational modeling of such processes can be numerically demanding due to the size of the systems and especially due to the timescales involved. Recently, a hybrid Brownian and molecular dynamics approach, i.e., Brownian dynamics including explicit atoms (BRODEA), has been developed and evaluated for studying the transport of monoatomic ions through membrane channels. Later on, this numerically efficient scheme has been applied to determine the free energy surfaces of the ciprofloxacin and enrofloxacin translocation through the porin OmpC using temperature-accelerated simulations. To improve the usability and accuracy of the approach, schemes to approximate the position-dependent diffusion constant of the molecule while traversing the pore had to be established. To this end, we have studied the translocation of the charged phosphonic acid antibiotic fosfomycin through the porin OmpF from Escherichia coli devising and benchmarking several diffusion models. To test the efficiency and sensitivity of these models, the effect of OmpF mutations on the permeation of fosfomycin was analyzed. Permeation events have been recorded over millisecond-long biased and unbiased simulations, from which thermodynamics and kinetics quantities of the translocation processes were determined. As a result, the use of the BRODEA approach, together with the appropriate diffusion model, was seen to accurately reproduce the findings observed in electrophysiology experiments and all-atom molecular dynamics simulations. These results suggest that the BRODEA approach can become a valuable tool for screening numerous compounds to evaluate their outer membrane permeability, a property important in the development of new antibiotics.

The soluble cytoplasmic N-terminal domain of the FocA channel gates bidirectional formate translocation

FocA belongs to the pentameric FNT (formate-nitrite transporter) superfamily of anion channels, translocating formate bidirectionally across the cytoplasmic membrane of Escherichia coli and other microorganisms. While the membrane-integral core of FocA shares considerable amino acid sequence conservation with other FNT family members, the soluble cytoplasmic N-terminal domain does not. To analyze the potential biochemical function of FocA's N-terminal domain in vivo, we constructed truncation derivatives and amino acid-exchange variants, and determined their ability to translocate formate across the membrane of E. coli cells by monitoring intracellular formate levels using a formate-sensitive reporter system. Analysis of strains synthesizing these FocA variants provided insights into formate efflux. Strains lacking the ability to generate formate intracellularly allowed us to determine whether these variants could import formate or its toxic chemical analog hypophosphite. Our findings reveal that the N-terminal domain of FocA is crucial for bidirectional FocA-dependent permeation of formate across the membrane. Moreover, we show that an amino acid sequence motif and secondary structural features of the flexible N-terminal domain are important for formate translocation, and efflux/influx is influenced by pyruvate formate-lyase. The soluble N-terminal domain is, therefore, essential for bidirectional formate translocation by FocA, suggesting a "gate-keeper" function controlling anion accessibility.

Voltage-Dependent Transport of Neutral Solutes through Nanopores: A Molecular View

The permeation of (neutral) molecules through nanopores in the presence of external voltages depends on several factors including pore electrostatics, electrophoretic force, and electro-osmotic drag. In earlier single-channel electrophysiology experiments, voltage-dependent asymmetric transport of neutral α-cyclodextrin (α-CD) molecules through the biological nanopore ΔCymA was observed. The voltage-dependent ion-associated flow of water, the so-called electro-osmotic flow, has been suggested to be the key factor behind the observed asymmetric behavior. The influence of pore electrostatics and electrophoretic force and their interplay with the electro-osmotic drag with varying buffers and voltages has not yet been analyzed at the molecular level. Hence, the detailed physical mechanism behind this intriguing permeation process is in part still unclear. In the present study, we have performed 36 μs all-atom free energy calculations by combining applied-field molecular dynamics simulations with metadynamics techniques. The influence of several ionic conditions as well as external voltages on the permeation of α-CD molecules across the ΔCymA pore has been investigated. To decipher the thermodynamic and kinetic details, the lowest energy paths and the permeation times for α-CD translocation have been estimated. In the presence of KCl or MgCl₂ salts, the charge of the cations is found to control the direction and magnitude of the electro-osmotic flow, which in turn strongly affects α-CD permeation. Overall, the present findings significantly improve the fundamental understanding of the voltage-dependent transport of neutral solutes across nanopores.

Water in Nanopores and Biological Channels: A Molecular Simulation Perspective

This Review explores the dynamic behavior of water within nanopores and biological channels in lipid bilayer membranes. We focus on molecular simulation studies, alongside selected structural and other experimental investigations. Structures of biological nanopores and channels are reviewed, emphasizing those high-resolution crystal structures, which reveal water molecules within the transmembrane pores, which can be used to aid the interpretation of simulation studies. Different levels of molecular simulations of water within nanopores are described, with a focus on molecular dynamics (MD). In particular, models of water for MD simulations are discussed in detail to provide an evaluation of their use in simulations of water in nanopores. Simulation studies of the behavior of water in idealized models of nanopores have revealed aspects of the organization and dynamics of nanoconfined water, including wetting/dewetting in narrow hydrophobic nanopores. A survey of simulation studies in a range of nonbiological nanopores is presented, including carbon nanotubes, synthetic nanopores, model peptide nanopores, track-etched nanopores in polymer membranes, and hydroxylated and functionalized nanoporous silica. These reveal a complex relationship between pore size/geometry, the nature of the pore lining, and rates of water transport. Wider nanopores with hydrophobic linings favor water flow whereas narrower hydrophobic pores may show dewetting. Simulation studies over the past decade of the behavior of water in a range of biological nanopores are described, including porins and β-barrel protein nanopores, aquaporins and related polar solute pores, and a number of different classes of ion channels. Water is shown to play a key role in proton transport in biological channels and in hydrophobic gating of ion channels. An overall picture emerges, whereby the behavior of water in a nanopore may be predicted as a function of its hydrophobicity and radius. This informs our understanding of the functions of diverse channel structures and will aid the design of novel nanopores. Thus, our current level of understanding allows for the design of a nanopore which promotes wetting over dewetting or vice versa. However, to design a novel nanopore, which enables fast, selective, and gated flow of water de novo would remain challenging, suggesting a need for further detailed simulations alongside experimental evaluation of more complex nanopore systems.

Transmembrane Facilitation of Lactate/H+ Instead of Lactic Acid Is Not a Question of Semantics but of Cell Viability

Transmembrane transport of monocarboxylates is conferred by structurally diverse membrane proteins. Here, we describe the pH dependence of lactic acid/lactate facilitation of an aquaporin (AQP9), a monocarboxylate transporter (MCT1, SLC16A1), and a formate–nitrite transporter (plasmodium falciparum FNT, PfFNT) in the equilibrium transport state. FNTs exhibit a channel-like structure mimicking the aquaporin-fold, yet act as secondary active transporters. We used radiolabeled lactate to monitor uptake via yeast-expressed AQP9, MCT1, and PfFNT for long enough time periods to reach the equilibrium state in which import and export rates are balanced. We confirmed that AQP9 behaved perfectly equilibrative for lactic acid, i.e., the neutral lactic acid molecule enters and passes the channel. MCT1, in turn, actively used the transmembrane proton gradient and acted as a lactate/H⁺ co-transporter. PfFNT behaved highly similar to the MCT in terms of transport properties, although it does not adhere to the classical alternating access transporter model. Instead, the FNT appears to use the proton gradient to neutralize the lactate anion in the protein’s vestibule to generate lactic acid in a place that traverses the central hydrophobic transport path. In conclusion, we propose to include FNT-type proteins into a more generalized, function-based transporter definition.

Nutrient Limitation Causes Differential Expression of Transport- and Metabolism Genes in the Compartmentalized Anammox Bacterium Kuenenia stuttgartiensis

Anaerobic ammonium-oxidizing (anammox) bacteria, members of the "Candidatus Brocadiaceae" family, play an important role in the nitrogen cycle and are estimated to be responsible for about half of the oceanic nitrogen loss to the atmosphere. Anammox bacteria combine ammonium with nitrite and produce dinitrogen gas via the intermediates nitric oxide and hydrazine (anammox reaction) while nitrate is formed as a by-product. These reactions take place in a specialized, membrane-enclosed compartment called the anammoxosome. Therefore, the substrates ammonium, nitrite and product nitrate have to cross the outer-, cytoplasmic-, and anammoxosome membranes to enter or exit the anammoxosome. The genomes of all anammox species harbor multiple copies of ammonium-, nitrite-, and nitrate transporter genes. Here we investigated how the distinct genes for ammonium-, nitrite-, and nitrate- transport were expressed during substrate limitation in membrane bioreactors. Transcriptome analysis of Kuenenia stuttgartiensis planktonic cells showed that four of the seven ammonium transporter homologs and two of the nine nitrite transporter homologs were significantly upregulated during ammonium-limited growth, while another ammonium transporter- and four nitrite transporter homologs were upregulated in nitrite limited growth conditions. The two nitrate transporters were expressed to similar levels in both conditions. In addition, genes encoding enzymes involved in the anammox reaction were differentially expressed, with those using nitrite as a substrate being upregulated under nitrite limited growth and those using ammonium as a substrate being upregulated during ammonium limitation. Taken together, these results give a first insight in the potential role of the multiple nutrient transporters in regulating transport of substrates and products in and out of the compartmentalized anammox cell.

Exploration of Free Energy Surfaces Across a Membrane Channel Using Metadynamics and Umbrella Sampling

To reach their site of action, it is essential for antibiotic molecules to cross the bacterial outer membrane. The progress of enhanced sampling techniques in molecular dynamics simulations enables us to understand these translocations at an atomic level. To this end, calculations of free energy surfaces for these permeation processes are of key importance. Herein, we investigate the translocation of a variety of anionic solutes through the outer membrane pore OprO of the Gram-negative bacterium Pseudomonas aeruginosa using the metadynamics and umbrella sampling techniques at the all-atom level. Free energy calculations have been performed employing these two distinct methods in order to illustrate the difference in computed free energies, if any. The investigated solutes range from a single atomic chloride ion over a multiatomic monophosphate ion to a more bulky fosmidomycin antibiotic. The role of complexity of the permeating solutes in estimating accurate free energy profiles is demonstrated by performing extensive convergence analysis. For simple monatomic ions, good agreement between the well-tempered metadynamics and the umbrella sampling approaches is achieved, while for the permeation of the monophosphate ion differences start to appear. In the case of larger molecules such as fosmidomycin it is a tough challenge to achieve converged free energy profiles. This issue is mainly due to neglecting orthogonal degrees of freedom during the free energy calculations. Nevertheless, the freely driven metadynamics approach leads to clearly advantageous results. Additionally, atomistic insights of the translocation mechanisms of all three solutes are discussed.

Is the E. coli Homolog of the Formate/Nitrite Transporter Family an Anion Channel? A Computational Study

Formate/nitrite transporters (FNTs) selectively transport monovalent anions and are found in prokaryotes and lower eukaryotes. They play a significant role in bacterial growth and act against the defense mechanism of infected hosts. Because FNTs do not occur in higher animals, they are attractive drug targets for many bacterial diseases. Phylogenetic analysis revealed that they can be classified into eight subgroups, two of which belong to the uncharacterized YfdC-α and YfdC-β groups. Experimentally determined structures of FNTs belonging to different phylogenetic groups adopt the unique aquaporin-like hourglass helical fold. We considered the formate channel from Vibrio cholerae, the hydrosulphide channel from Clostridium difficile, and the uncharacterized channel from Escherichia coli (EcYfdC) to investigate the mechanism of transport and selectivity. Using equilibrium molecular dynamics and umbrella sampling studies, we determined temporal channel radius profiles, permeation events, and potential of mean force profiles of different substrates with the conserved central histidine residue in protonated or neutral form. Unlike the formate channel from V. cholerae and the hydrosulphide channel from C. difficile, molecular dynamics studies showed that the formate substrate was unable to enter the vestibule region of EcYfdC. Absence of a conserved basic residue and presence of acidic residues in the vestibule regions, conserved only in YfdC-α, were found to be responsible for high energy barriers for the anions to enter EcYfdC. Potential of mean force profiles generated for ammonia and ammonium ion revealed that EcYfdC can transport neutral solutes and could possibly be involved in the transport of cations analogous to the mechanism proposed for ammonium transporters. Although YfdC members belong to the FNT family, our studies strongly suggest that EcYfdC is not an anion channel. Absence or presence of specific charged residues at particular positions makes EcYfdC selective for neutral or possibly cationic substrates. Further experimental studies are needed to get a definitive answer to the question of the substrate selectivity of EcYfdC. This provides an example of membrane proteins from the same family transporting substrates of different chemical nature.

APH-1A Component of γ-Secretase Forms an Internal Water and Ion-Containing Cavity

Anterior pharynx-defective 1A (APH-1A) is a seven transmembrane component of γ-secretase (GS), an aspartyl protease enzyme involved in the production of toxic amyloid-β peptides in Alzheimer's disease patients. Cryo-electron microscopy structures of the enzyme complex revealed a central cavity in its APH-1A component, similar to water-containing cavities in G-protein coupled receptors (GPCRs). In this work, we performed molecular dynamics and umbrella sampling simulations to understand the role of the APH-1A cavity in the GS complex. Our results suggest that APH-1A is able to store water molecules in its inner cavity and transport some of them between cell spaces. Additionally, APH-1A allows the influx of extracellular cations into a central hydrophilic cavity but cannot transport them into the intracellular space. Overall, this study seeks to describe an alternative APH-1A function in GS besides its complex stabilization role and provide novel approaches to understand the functioning of the GS enzyme.

Formate-nitrite transporters carrying nonprotonatable amide amino acids instead of a central histidine maintain pH-dependent transport

Microbial formate-nitrite transporter-type proteins (FNT) exhibit dual transport functionality. At neutral pH, electrogenic anion currents are detectable, whereas upon acidification transport of the neutral, protonated monoacid predominates. Physiologically, FNT-mediated proton co-transport is vital when monocarboxylic acid products of the energy metabolism, such as l-lactate, are released from the cell. Accordingly, Plasmodium falciparum malaria parasites can be killed by small-molecule inhibitors of PfFNT. Two opposing hypotheses on the site of substrate protonation are plausible. The proton relay mechanism postulates proton transfer from a highly conserved histidine centrally positioned in the transport path. The dielectric slide mechanism assumes decreasing acidity of substrates entering the lipophilic vestibules and protonation via the bulk water. Here, we defined the transport mechanism of the FNT from the amoebiasis parasite Entamoeba histolytica, EhFNT, and also show that BtFdhC from Bacillus thuringiensis is a functional formate transporter. Both FNTs carry a nonprotonatable amide amino acid, asparagine or glutamine, respectively, at the central histidine position. Despite having a nonprotonatable residue, EhFNT displayed the same substrate selectivity for larger monocarboxylates including l-lactate, a low substrate affinity as is typical for FNTs, and, strikingly, proton motive force-dependent transport as observed for PfFNT harboring a central histidine. These results argue against a proton relay mechanism, indicating that substrate protonation must occur outside of the central histidine region, most likely in the vestibules. Furthermore, EhFNT is the sole annotated FNT in the Entamoeba genome suggesting that it could be a putative new drug target with similar utility as that of the malarial PfFNT.