Knowing the enemy: homoacetogens in hydrogen production reactors.

One of the bottlenecks of the hydrogen production by dark fermentation is the low yields obtained because of the homoacetogenesis persistence, a metabolic pathway where H2 and CO2 are consumed to produce acetate. The central reactions of H2 production and homoacetogenesis are catalyzed by enzyme hydrogenase and the formyltetrahydrofolate synthetase, respectively. In this work, genes encoding for the formyltetrahydrofolate synthetase (fthfs) and hydrogenase (hydA) were used to investigate the diversity of homoacetogens as well as their phylogenetic relationships through quantitative PCR (qPCR) and next-generation amplicon sequencing. A total of 70 samples from 19 different H2-producing bioreactors with different configurations and operating conditions were analyzed. Quantification through qPCR showed that the abundance of fthfs and hydA was strongly associated with the type of substrate, organic loading rate, and H2 production performance. In particular, fthfs sequencing revealed that homoacetogens diversity was low with one or two dominant homoacetogens in each sample. Clostridium carboxivorans was detected in the reactors fed with agave hydrolisates; Acetobacterium woodii dominated in systems fed with glucose; Blautia coccoides and unclassified Sporoanaerobacter species were present in reactors fed with cheese whey; finally, Eubacterium limosum and Selenomonas sp. were co-dominant in reactors fed with glycerol. Altogether, quantification and sequencing analysis revealed that the occurrence of homoacetogenesis could take place due to (1) metabolic changes of H2-producing bacteria towards homoacetogenesis or (2) the displacement of H2-producing bacteria by homoacetogens. Overall, it was demonstrated that the fthfs gene was a suitable marker to investigate homoacetogens in H2-producing reactors. KEY POINTS: • qPCR and sequencing analysis revealed two homoacetogenesis phenomena. • fthfs gene was a suitable marker to investigate homoacetogens in H2 reactors.

Continuous thermophilic hydrogen production from an enzymatic hydrolysate of agave bagasse: Inoculum origin, homoacetogenesis and microbial community analysis.

In this research, the performance of two thermophilic inocula of different origin on continuous hydrogen production from an enzymatic hydrolysate of agave bagasse were compared; one of them was obtained from a thermophilic reactor and the second one was taken from a mesophilic reactor and acclimated to thermophilic conditions. The acclimation process in one-step quickly established a high-performance hydrogen producing community, obtaining a volumetric hydrogen production rate of 3811 ± 19 mL H2/L-d with an hydrogen yield of 121 L H2/kg bagasse compared to 1473 ± 6 mL H2/L-d and 26.6 L H2/kg obtained with the thermophilic-origin inoculum. The differences in the performance of both inocula were closely linked to the profile of volatile fatty acids produced, the homoacetogenic pathway and the microbial community, the latter being the determining factor. The use of mesophilic-origin inoculum acclimated to thermophilic conditions can significantly improve the hydrogen production from lignocellulosic bagasse.

Standardized protocol for determination of biohydrogen potential.

Biohydrogen production potential (BHP) depends on several factors like inoculum source, substrate, pH, among many others. Batch assays are the most common strategy to evaluate such parameters, where the comparison is a challenging task due to the different procedures used. The present method introduces the first internationally validated protocol, evaluated by 8 independent laboratories from 5 different countries, to assess the biohydrogen potential. As quality criteria, a coefficient of variation of the cumulative hydrogen production (H max) was defined to be <15 %. Two options to run BHP batch tests were proposed; a manual protocol with periodic measurements of biogas production, needing conventional laboratory materials and analytical equipment for biogas characterization; and an automatic protocol, which is run in a device developed for online measurements of low biogas production. The detailed procedures for both protocol options are presented, as well as data validating them. The validation showed acceptable repeatability and reproducibility, measured as intra- and inter-laboratory coefficient of variation, which can be reduced up to 9 %.

Improvement of continuous hydrogen production using individual and binary enzymatic hydrolysates of agave bagasse in suspended-culture and biofilm reactors.

Continuous hydrogen (H2) production from individual (Stonezyme, IH) and binary (Celluclast-Viscozyme, BH) enzymatic hydrolysates of agave bagasse was evaluated in continuous stirred-tank reactors (CSTR) and trickling bed reactors (TBR). The volumetric H2 production rates (VHPR) in CSTR were 13 and 2.25 L H2/L-d with BH and IH, respectively. Meanwhile, VHPR of 5.76 and 2.0 L H2/L-d were obtained in the TBR configuration using BH and IH, respectively. Differences on VHPR between reactors could be explained by substrate availability, which is intrinsic to the growth mode of each reactor configuration; while differences of VHPR between hydrolysates were possibly related to the composition of enzymatic hydrolysates. Furthermore, homoacetogenesis was strongly influenced by H2 and substrate transfer conditions. Considering VHPR, H2 yields, and costs of hydrolysis, hydrogen production from binary hydrolysates of agave bagasse was identified as the most promising alternative evaluated with scale-up potential for the production of energy biofuels.

Microbial contamination in methanol biofilters inoculated with a pure strain of Pichia pastoris: A potential limitation for waste revalorization.

Novel biotechnologies to valorize waste emissions are based on the use of specialized microbial groups that produce different compounds of industrial interest. On this scenario, the retention of such specific microorganisms in the system is of critical interest; however, the potential limitations of working with simplified cultures in a competitive open environment are neither fully explored nor well understood. In this work, a series of biofilters treating methanol vapors coupled with heterologous endochitinase production were used to evaluate the performance of a specialized microbial population during a typical open-to-environment operation. The biofilters were inoculated with a transformed strain of Pichia pastoris and were operated identically for about 90 days. The results showed that the biofiltration performance became diverse with time in terms of the elimination capacity (EC) shifting from a variation coefficient of 1.5% (EC = 274 ± 24, 279 ± 5, and 281.9 ± 25 g/[m3 h]) at the beginning of the operation to 33% (EC = 297 ± 9, 338 ± 7, and 341 ± 2 g/[m3 h]) at the end of operation. Epifluorescence analysis and cloning-sequencing suggested that P. pastoris remained as the dominant microorganism of methanol degradation, whereas diverse airborne bacteria, including Ochrobactrum spp. and Klebsiella oxytoca, played a secondary role possibly associated with the consumption of intermediates. Overall, this study found that low diversity systems operated under non-sterile conditions could be susceptible to contamination with external microorganisms causing a diversifying behavior at the performance and microbial community levels. © 2018 American Institute of Chemical Engineers Biotechnol. Prog., 35: e2715, 2019.

GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics Simulations.

We introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged density maps from molecular dynamics simulations. GROmaρs efficiently computes density maps by fast multi-Gaussian spreading of atomic densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of atomic average localization during the simulations. Most importantly, it allows the comparison between computed and reference maps (e.g., experimental) through calculation of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quantitatively contrast perturbed and control simulation data sets and to examine how much biomolecular systems resemble both synthetic and experimental density maps. This was especially advantageous for multimolecule systems in which standard comparisons like RMSDs are difficult to compute. In addition, GROmaρs incorporates absolute and relative spatial free-energy estimates to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the density calculation. Furthermore, masking of regions was implemented to speed up calculations and to facilitate the comparison with experimental maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and average charge distribution in biomolecular systems. We employed all these functionalities to inspect the localization of lipid and water molecules in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, we anticipate a high applicability of GROmaρs for the analysis of molecular dynamics simulations and their comparison with experimentally determined densities.

Hydrogen metabolic patterns driven by Clostridium-Streptococcus community shifts in a continuous stirred tank reactor.

The hydrogen (H2) production efficiency in dark fermentation systems is strongly dependent on the occurrence of metabolic pathways derived from the selection of microbial species that either consume molecular H2 or outcompete hydrogenogenic bacteria for the organic substrate. In this study, the effect of organic loading rate (OLR) on the H2 production performance, the metabolic pathways, and the microbial community composition in a continuous system was evaluated. Two bacterial genera, Clostridium and Streptococcus, were dominant in the microbial community depending on the OLR applied. At low OLR (14.7-44.1 gLactose/L-d), Clostridium sp. was dominant and directed the system towards the acetate-butyrate fermentation pathway, with a maximum H2 yield of 2.14 molH2/molHexose obtained at 29.4 gLactose/L-d. Under such conditions, the volumetric hydrogen production rate (VHPR) was between 3.2 and 11.6 LH2/L-d. In contrast, relatively high OLR (58.8 and 88.2 gLactose/L-d) favored the dominance of Streptococcus sp. as co-dominant microorganism leading to lactate production. Under these conditions, the formate production was also stimulated serving as a strategy to dispose the surplus of reduced molecules (e.g., NADH2+), which theoretically consumed up to 5.72 LH2/L-d. In such scenario, the VHPR was enhanced (13.7-14.5 LH2/L-d) but the H2 yield dropped to a minimum of 0.74 molH2/molHexose at OLR = 58.8 gLactose/L-d. Overall, this research brings clear evidence of the intrinsic occurrence of metabolic pathways detrimental for biohydrogen production, i.e., lactic acid fermentation and formate production, suggesting the use of low OLR as a strategy to control them.

Gating Charge Calculations by Computational Electrophysiology Simulations.

Electrical cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the molecular details of voltage sensing remain insufficiently understood. Here, we present a molecular dynamics simulation-based method to determine the underlying charge movement across the membrane-the gating charge-by measuring electrical capacitor properties of membrane-embedded proteins. We illustrate the approach by calculating the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K+ channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against experimental data. We then use the gating charge analysis to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na+ ions in Na+-coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with experiments, and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions.

Localization and Ordering of Lipids Around Aquaporin-0: Protein and Lipid Mobility Effects.

Hydrophobic matching, lipid sorting, and protein oligomerization are key principles by which lipids and proteins organize in biological membranes. The Aquaporin-0 channel (AQP0), solved by electron crystallography (EC) at cryogenic temperatures, is one of the few protein-lipid complexes of which the structure is available in atomic detail. EC and room-temperature molecular dynamics (MD) of dimyristoylglycerophosphocholine (DMPC) annular lipids around AQP0 show similarities, however, crystal-packing and temperature might affect the protein surface or the lipids distribution. To understand the role of temperature, lipid phase, and protein mobility in the localization and ordering of AQP0-lipids, we used MD simulations of an AQP0-DMPC bilayer system. Simulations were performed at physiological and at DMPC gel-phase temperatures. To decouple the protein and lipid mobility effects, we induced gel-phase in the lipids or restrained the protein. We monitored the lipid ordering effects around the protein. Reducing the system temperature or inducing lipid gel-phase had a marginal effect on the annular lipid localization. However, restraining the protein mobility increased the annular lipid localization around the whole AQP0 surface, resembling EC. The distribution of the inter-phosphate and hydrophobic thicknesses showed that stretching of the DMPC annular layer around AQP0 surface is the mechanism that compensates the hydrophobic mismatch in this system. The distribution of the local area-per-lipid and the acyl-chain order parameters showed particular fluid- and gel-like areas that involved several lipid layers. These areas were in contact with the surfaces of higher and lower protein mobility, respectively. We conclude that the AQP0 surfaces induce specific fluid- and gel-phase prone areas. The presence of these areas might guide the AQP0 lipid sorting interactions with other membrane components, and is compatible with the squared array oligomerization of AQP0 tetramers separated by a layer of annular lipids.

Voltage Dependence of Conformational Dynamics and Subconducting States of VDAC-1.

The voltage-dependent anion channel 1 (VDAC-1) is an important protein of the outer mitochondrial membrane that transports energy metabolites and is involved in apoptosis. The available structures of VDAC proteins show a wide ß-stranded barrel pore, with its N-terminal α-helix (N-α) bound to its interior. Electrophysiology experiments revealed that voltage, its polarity, and membrane composition modulate VDAC currents. Experiments with VDAC-1 mutants identified amino acids that regulate the gating process. However, the mechanisms for how these factors regulate VDAC-1, and which changes they trigger in the channel, are still unknown. In this study, molecular dynamics simulations and single-channel experiments of VDAC-1 show agreement for the current-voltage relationships of an "open" channel and they also show several subconducting transient states that are more cation selective in the simulations. We observed voltage-dependent asymmetric distortions of the VDAC-1 barrel and the displacement of particular charged amino acids. We constructed conformational models of the protein voltage response and the pore changes that consistently explain the protein conformations observed at opposite voltage polarities, either in phosphatidylethanolamine or phosphatidylcholine membranes. The submicrosecond VDAC-1 voltage response shows intrinsic structural changes that explain the role of key gating amino acids and support some of the current gating hypotheses. These voltage-dependent protein changes include asymmetric barrel distortion, its interaction with the membrane, and significant displacement of N-α amino acids.

The membrane anchor of the transcriptional activator SREBP is characterized by intrinsic conformational flexibility.

Regulated intramembrane proteolysis (RIP) is a conserved mechanism crucial for numerous cellular processes, including signaling, transcriptional regulation, axon guidance, cell adhesion, cellular stress responses, and transmembrane protein fragment degradation. Importantly, it is relevant in various diseases including Alzheimer's disease, cardiovascular diseases, and cancers. Even though a number of structures of different intramembrane proteases have been solved recently, fundamental questions concerning mechanistic underpinnings of RIP and therapeutic interventions remain. In particular, this includes substrate recognition, what properties render a given substrate amenable for RIP, and how the lipid environment affects the substrate cleavage. Members of the sterol regulatory element-binding protein (SREBP) family of transcription factors are critical regulators of genes involved in cholesterol/lipid homeostasis. After site-1 protease cleavage of the inactive SREBP transmembrane precursor protein, RIP of the anchor intermediate by site-2 protease generates the mature transcription factor. In this work, we have investigated the labile anchor intermediate of SREBP-1 using NMR spectroscopy. Surprisingly, NMR chemical shifts, site-resolved solvent exposure, and relaxation studies show that the cleavage site of the lipid-signaling protein intermediate bears rigid α-helical topology. An evolutionary conserved motif, by contrast, interrupts the secondary structure ∼9-10 residues C-terminal of the scissile bond and acts as an inducer of conformational flexibility within the carboxyl-terminal transmembrane region. These results are consistent with molecular dynamics simulations. Topology, stability, and site-resolved dynamics data suggest that the cleavage of the α-helical substrate in the case of RIP may be associated with a hinge motion triggered by the molecular environment.

ClC-K chloride channels: emerging pathophysiology of Bartter syndrome type 3.

The mutations in the CLCNKB gene encoding the ClC-Kb chloride channel are responsible for Bartter syndrome type 3, one of the four variants of Bartter syndrome in the genetically based nomenclature. All forms of Bartter syndrome are characterized by hypokalemia, metabolic alkalosis, and secondary hyperaldosteronism, but Bartter syndrome type 3 has the most heterogeneous presentation, extending from severe to very mild. A relatively large number of CLCNKB mutations have been reported, including gene deletions and nonsense or missense mutations. However, only 20 CLCNKB mutations have been functionally analyzed, due to technical difficulties regarding ClC-Kb functional expression in heterologous systems. This review provides an overview of recent progress in the functional consequences of CLCNKB mutations on ClC-Kb chloride channel activity. It has been observed that 1) all ClC-Kb mutants have an impaired expression at the membrane; and 2) a minority of the mutants combines reduced membrane expression with altered pH-dependent channel gating. Although further investigation is needed to fully characterize disease pathogenesis, Bartter syndrome type 3 probably belongs to the large family of conformational diseases, in which the mutations destabilize channel structure, inducing ClC-Kb retention in the endoplasmic reticulum and accelerated channel degradation.

Mechanisms of anion conduction by coupled glutamate transporters.

Excitatory amino acid transporters (EAATs) are essential for terminating glutamatergic synaptic transmission. They are not only coupled glutamate/Na(+)/H(+)/K(+) transporters but also function as anion-selective channels. EAAT anion channels regulate neuronal excitability, and gain-of-function mutations in these proteins result in ataxia and epilepsy. We have combined molecular dynamics simulations with fluorescence spectroscopy of the prokaryotic homolog GltPh and patch-clamp recordings of mammalian EAATs to determine how these transporters conduct anions. Whereas outward- and inward-facing GltPh conformations are nonconductive, lateral movement of the glutamate transport domain from intermediate transporter conformations results in formation of an anion-selective conduction pathway. Fluorescence quenching of inserted tryptophan residues indicated the entry of anions into this pathway, and mutations of homologous pore-forming residues had analogous effects on GltPh simulations and EAAT2/EAAT4 measurements of single-channel currents and anion/cation selectivities. These findings provide a mechanistic framework of how neurotransmitter transporters can operate as anion-selective and ligand-gated ion channels.

CLCNKB mutations causing mild Bartter syndrome profoundly alter the pH and Ca2+ dependence of ClC-Kb channels.

ClC-Kb, a member of the ClC family of Cl(-) channels/transporters, plays a major role in the absorption of NaCl in the distal nephron. CLCNKB mutations cause Bartter syndrome type 3, a hereditary renal salt-wasting tubulopathy. Here, we investigate the functional consequences of a Val to Met substitution at position 170 (V170M, α helix F), which was detected in eight patients displaying a mild phenotype. Conductance and surface expression were reduced by ~40-50 %. The regulation of channel activity by external H(+) and Ca(2+) is a characteristic property of ClC-Kb. Inhibition by external H(+) was dramatically altered, with pKH shifting from 7.6 to 6.0. Stimulation by external Ca(2+) on the other hand was no longer detectable at pH 7.4, but was still present at acidic pH values. Functionally, these regulatory modifications partly counterbalance the reduced surface expression by rendering V170M hyperactive. Pathogenic Met170 seems to interact with another methionine on α helix H (Met227) since diverse mutations at this site partly removed pH sensitivity alterations of V170M ClC-Kb. Exploring other disease-associated mutations, we found that a Pro to Leu substitution at position 124 (α helix D, Simon et al., Nat Genet 1997, 17:171-178) had functional consequences similar to those of V170M. In conclusion, we report here for the first time that ClC-Kb disease-causing mutations located around the selectivity filter can result in both reduced surface expression and hyperactivity in heterologous expression systems. This interplay must be considered when analyzing the mild phenotype of patients with type 3 Bartter syndrome.

Computational analysis of local membrane properties.

In the field of biomolecular simulations, dynamics of phospholipid membranes is of special interest. A number of proteins, including channels, transporters, receptors and short peptides are embedded in lipid bilayers and tightly interact with phospholipids. While the experimental measurements report on the spatial and/or temporal average membrane properties, simulation results are not restricted to the average properties. In the current study, we present a collection of methods for an efficient local membrane property calculation, comprising bilayer thickness, area per lipid, deuterium order parameters, Gaussian and mean curvature. The local membrane property calculation allows for a direct mapping of the membrane features, which subsequently can be used for further analysis and visualization of the processes of interest. The main features of the described methods are highlighted in a number of membrane systems, namely: a pure dimyristoyl-phosphatidyl-choline (DMPC) bilayer, a fusion peptide interacting with a membrane, voltage-dependent anion channel protein embedded in a DMPC bilayer, cholesterol enriched bilayer and a coarse grained simulation of a curved palmitoyl-oleoyl-phosphatidyl-choline lipid membrane. The local membrane property analysis proves to provide an intuitive and detailed view on the observables that are otherwise interpreted as averaged bilayer properties.

Partial least-squares functional mode analysis: application to the membrane proteins AQP1, Aqy1, and CLC-ec1.

We introduce an approach based on the recently introduced functional mode analysis to identify collective modes of internal dynamics that maximally correlate to an external order parameter of functional interest. Input structural data can be either experimentally determined structure ensembles or simulated ensembles, such as molecular dynamics trajectories. Partial least-squares regression is shown to yield a robust solution to the multidimensional optimization problem, with a minimal and controllable risk of overfitting, as shown by extensive cross-validation. Several examples illustrate that the partial least-squares-based functional mode analysis successfully reveals the collective dynamics underlying the fluctuations in selected functional order parameters. Applications to T4 lysozyme, the Trp-cage, the aquaporin channels Aqy1 and hAQP1, and the CLC-ec1 chloride antiporter are presented in which the active site geometry, the hydrophobic solvent-accessible surface, channel gating dynamics, water permeability (p(f)), and a dihedral angle are defined as functional order parameters. The Aqy1 case reveals a gating mechanism that connects the inner channel gating residues with the protein surface, thereby providing an explanation of how the membrane may affect the channel. hAQP1 shows how the p(f) correlates with structural changes around the aromatic/arginine region of the pore. The CLC-ec1 application shows how local motions of the gating Glu(148) couple to a collective motion that affects ion affinity in the pore.

ß-Barrel mobility underlies closure of the voltage-dependent anion channel.

The voltage-dependent anion channel (VDAC) is the major protein in the outer mitochondrial membrane, where it mediates transport of ATP and ADP. Changes in its permeability, induced by voltage or apoptosis-related proteins, have been implicated in apoptotic pathways. The three-dimensional structure of VDAC has recently been determined as a 19-stranded ß-barrel with an in-lying N-terminal helix. However, its gating mechanism is still unclear. Using solid-state NMR spectroscopy, molecular dynamics simulations, and electrophysiology, we show that deletion of the rigid N-terminal helix sharply increases overall motion in VDAC's ß-barrel, resulting in elliptic, semicollapsed barrel shapes. These states quantitatively reproduce conductance and selectivity of the closed VDAC conformation. Mutation of the N-terminal helix leads to a phenotype intermediate to the open and closed states. These data suggest that the N-terminal helix controls entry into elliptic ß-barrel states which underlie VDAC closure. Our results also indicate that ß-barrel channels are intrinsically flexible.

Molecular driving forces defining lipid positions around aquaporin-0.

Lipid-protein interactions play pivotal roles in biological membranes. Electron crystallographic studies of the lens-specific water channel aquaporin-0 (AQP0) revealed atomistic views of such interactions, by providing high-resolution structures of annular lipids surrounding AQP0. It remained unclear, however, whether these lipid structures are representative of the positions of unconstrained lipids surrounding an individual protein, and what molecular determinants define the lipid positions around AQP0. We addressed these questions by using molecular dynamics simulations and crystallographic refinement, and calculated time-averaged densities of dimyristoyl-phosphatidylcholine lipids around AQP0. Our simulations demonstrate that, although the experimentally determined crystallographic lipid positions are constrained by the crystal packing, they appropriately describe the behavior of unconstrained lipids around an individual AQP0 tetramer, and thus likely represent physiologically relevant lipid positions.While the acyl chains were well localized, the lipid head groups were not. Furthermore, in silico mutations showed that electrostatic interactions do not play a major role attracting these phospholipids towards AQP0. Instead, the mobility of the protein crucially modulates the lipid localization and explains the difference in lipid density between extracellular and cytoplasmic leaflets. Moreover, our simulations support a general mechanism in which membrane proteins laterally diffuse accompanied by several layers of localized lipids, with the positions of the annular lipids being influenced the most by the protein surface. We conclude that the acyl chains rather than the head groups define the positions of dimyristoyl-phosphatidylcholine lipids around AQP0. Lipid localization is largely determined by the mobility of the protein surface, whereas hydrogen bonds play an important but secondary role.

Functional dynamics in the voltage-dependent anion channel.

The voltage-dependent anion channel (VDAC), located in the outer mitochondrial membrane, acts as a gatekeeper for the entry and exit of mitochondrial metabolites. Here we reveal functional dynamics of isoform one of VDAC (VDAC1) by a combination of solution NMR spectroscopy, Gaussian network model analysis, and molecular dynamics simulation. Micro- to millisecond dynamics are significantly increased for the N-terminal six ß-strands of VDAC1 in micellar solution, in agreement with increased B-factors observed in the same region in the bicellar crystal structure of VDAC1. Molecular dynamics simulations reveal that a charge on the membrane-facing glutamic acid 73 (E73) accounts for the elevation of N-terminal protein dynamics as well as a thinning of the nearby membrane. Mutation or chemical modification of E73 strongly reduces the micro- to millisecond dynamics in solution. Because E73 is necessary for hexokinase-I-induced VDAC channel closure and inhibition of apoptosis, our results imply that micro- to millisecond dynamics in the N-terminal part of the barrel are essential for VDAC interaction and gating.

Voltage-dependent and -independent titration of specific residues accounts for complex gating of a ClC chloride channel by extracellular protons.

The ClC transport protein family comprises both Cl(-) ion channel and H(+)/Cl(-) and H(+)/NO(3)(-) exchanger members. Structural studies on a bacterial ClC transporter reveal a pore obstructed at its external opening by a glutamate side-chain which acts as a gate for Cl(-) passage and in addition serves as a staging post for H(+) exchange. This same conserved glutamate acts as a gate to regulate Cl(-) flow in ClC channels. The activity of ClC-2, a genuine Cl(-) channel, has a biphasic response to extracellular pH with activation by moderate acidification followed by abrupt channel closure at pH values lower than approximately 7. We have now investigated the molecular basis of this complex gating behaviour. First, we identify a sensor that couples extracellular acidification to complete closure of the channel. This is extracellularly-facing histidine 532 at the N-terminus of transmembrane helix Q whose neutralisation leads to channel closure in a cooperative manner. We go on to show that acidification-dependent activation of ClC-2 is voltage dependent and probably mediated by protonation of pore gate glutamate 207. Intracellular Cl(-) acts as a voltage-independent modulator, as though regulating the pK(a) of the protonatable residue. Our results suggest that voltage dependence of ClC-2 is given by hyperpolarisation-dependent penetration of protons from the extracellular side to neutralise the glutamate gate deep within the channel, which allows Cl(-) efflux. This is reminiscent of a partial exchanger cycle, suggesting that the ClC-2 channel evolved from its transporter counterparts.