Mechanistic basis of the dynamic response of TWIK1 ionic selectivity to pH.

Highly selective for K+ at neutral pH, the TWIK1 channel becomes permeable to Na+ upon acidification. Using molecular dynamics simulations, we identify a network of residues involved in this unique property. Between the open and closed states previously observed by electron microscopy, molecular dynamics simulations show that the channel undergoes conformational changes between pH 7.5-6 involving residues His122, Glu235, Lys246 and Phe109. A complex network of interactions surrounding the selectivity filter at high pH transforms into a simple set of stronger interactions at low pH. In particular, His122 protonated by acidification moves away from Lys246 and engages in a salt bridge with Glu235. In addition, stacking interactions between Phe109 and His122, which stabilize the selectivity filter in its K+-selective state at high pH, disappear upon acidification. This leads to dissociation of the Phe109 aromatic side chain from this network, resulting in the Na+-permeable conformation of the channel.

Case report: 'Atypical Richter transformation from CLL-type monoclonal B-cell lymphocytosis into Burkitt lymphoma in a treatment naïve patient'.

Background: Richter transformation refers to the progression of an initially slow-growing small lymphocytic lymphoma/chronic lymphocytic leukemia (SLL/CLL) into an aggressive lymphoma, typically diffuse large B-cell lymphoma (DLBCL) or Hodgkin lymphoma. Case presentation: The patient presented with a rapid onset of localized cervical swelling, accompanied by monoclonal B-cell lymphocytosis displaying a CLL immunophenotype. The histopathological analysis identified a Burkitt lymphoma (BL) located in the submandibular gland and adjacent lymph node. The patient's bone marrow displayed a minor infiltration of monoclonal B-cells with a CLL immunophenotype (< 10%). Molecular analysis demonstrated the presence of the same monoclonal rearrangement in the framework region (FR3 region) of the immunoglobulin heavy chain (IGH) locus. High-throughput sequencing of the immunoglobulin heavy and light chains also confirmed the presence of the same rearrangement in SLL/CLL and in the Burkitt lymphoma sample, but also highlighted the presence of a second rearrangement in the Burkitt lymphoma cells, not shared with the SLL/CLL cells in the bone marrow. The patient was treated with DA-EPOCH-R, which lead to a complete metabolic response. Conclusion: This report provides an exceptionally rare description of a CLL-type monoclonal B-cell lymphocytosis transforming into a very aggressive Burkitt lymphoma in a treatment naïve patient.

Autoimmunity and immunodeficiency associated with monoallelic LIG4 mutations via haploinsufficiency.

BACKGROUND: Biallelic mutations in LIG4 encoding DNA-ligase 4 cause a rare immunodeficiency syndrome manifesting as infant-onset life-threatening and/or opportunistic infections, skeletal malformations, radiosensitivity and neoplasia. LIG4 is pivotal during DNA repair and during V(D)J recombination as it performs the final DNA-break sealing step. OBJECTIVES: This study explored whether monoallelic LIG4 missense mutations may underlie immunodeficiency and autoimmunity with autosomal dominant inheritance. METHODS: Extensive flow-cytometric immune-phenotyping was performed. Rare variants of immune system genes were analyzed by whole exome sequencing. DNA repair functionality and T-cell-intrinsic DNA damage tolerance was tested with an ensemble of in vitro and in silico tools. Antigen-receptor diversity and autoimmune features were characterized by high-throughput sequencing and autoantibody arrays. Reconstitution of wild-type versus mutant LIG4 were performed in LIG4 knockout Jurkat T cells, and DNA damage tolerance was subsequently assessed. RESULTS: A novel heterozygous LIG4 loss-of-function mutation (p.R580Q), associated with a dominantly inherited familial immune-dysregulation consisting of autoimmune cytopenias, and in the index patient with lymphoproliferation, agammaglobulinemia, and adaptive immune cell infiltration into nonlymphoid organs. Immunophenotyping revealed reduced naive CD4+ T cells and low TCR-Vα7.2+ T cells, while T-/B-cell receptor repertoires showed only mild alterations. Cohort screening identified 2 other nonrelated patients with the monoallelic LIG4 mutation p.A842D recapitulating clinical and immune-phenotypic dysregulations observed in the index family and displaying T-cell-intrinsic DNA damage intolerance. Reconstitution experiments and molecular dynamics simulations categorize both missense mutations as loss-of-function and haploinsufficient. CONCLUSIONS: This study provides evidence that certain monoallelic LIG4 mutations may cause human immune dysregulation via haploinsufficiency.

Calcium regulates acid-sensing ion channel 3 activation by competing with protons in the channel pore and at an allosteric binding site.

The extracellular Ca2+ concentration changes locally under certain physiological and pathological conditions. Such variations affect the function of ion channels of the nervous system and consequently also neuronal signalling. We investigated here the mechanisms by which Ca2+ controls the activity of acid-sensing ion channel (ASIC) 3. ASICs are neuronal, H+-gated Na+ channels involved in several physiological and pathological processes, including the expression of fear, learning, pain sensation and neurodegeneration after ischaemic stroke. It was previously shown that Ca2+ negatively modulates the ASIC pH dependence. While protons are default activators of ASIC3, this channel can also be activated at pH7.4 by the removal of the extracellular Ca2+. Two previous studies concluded that low pH opens ASIC3 by displacing Ca2+ ions that block the channel pore at physiological pH. We show here that an acidic residue, distant from the pore, together with pore residues, controls the modulation of ASIC3 by Ca2+. Our study identifies a new regulatory site in ASIC3 and demonstrates that ASIC3 activation involves an allosteric mechanism together with Ca2+ unbinding from the channel pore. We provide a molecular analysis of a regulatory mechanism found in many ion channels.

Galvani Offset Potential and Constant-pH Simulations of Membrane Proteins.

A central problem in computational biophysics is the treatment of titratable residues in molecular dynamics simulations of large biological macromolecular systems. Conventional simulation methods ascribe a fixed ionization state to titratable residues in accordance with their pKa and the pH of the system, assuming that an effective average model will be able to capture the predominant behavior of the system. While this assumption may be justifiable in many cases, it is certainly limited, and it is important to design alternative methodologies allowing a more realistic treatment. Constant-pH simulation methods provide powerful approaches to handle titratable residues more realistically by allowing the ionization state to vary statistically during the simulation. Extending the molecular mechanical (MM) potential energy function to a family of potential functions accounting for different ionization states, constant-pH simulations are designed to sample all accessible configurations and ionization states, properly weighted according to their Boltzmann factor. Because protonation and deprotonation events correspond to a change in the total charge, difficulties arise when the long-range Coulomb interaction is treated on the basis of an idealized infinite simulation model and periodic boundary conditions with particle-mesh Ewald lattice sums. Charging free-energy calculations performed under these conditions in aqueous solution depend on the Galvani potential of the bulk water phase. This has important implications for the equilibrium and nonequilibrium constant-pH simulation methods grounded in the relative free-energy difference corresponding to the protonated and unprotonated residues. Here, the effect of the Galvani potential is clarified, and a simple practical solution is introduced to address this issue in constant-pH simulations of the acid-sensing ion channel (ASIC).

Kinetic analysis of ASIC1a delineates conformational signaling from proton-sensing domains to the channel gate.

Acid-sensing ion channels (ASICs) are neuronal Na+ channels that are activated by a drop in pH. Their established physiological and pathological roles, involving fear behaviors, learning, pain sensation, and neurodegeneration after stroke, make them promising targets for future drugs. Currently, the ASIC activation mechanism is not understood. Here, we used voltage-clamp fluorometry (VCF) combined with fluorophore-quencher pairing to determine the kinetics and direction of movements. We show that conformational changes with the speed of channel activation occur close to the gate and in more distant extracellular sites, where they may be driven by local protonation events. Further, we provide evidence for fast conformational changes in a pathway linking protonation sites to the channel pore, in which an extracellular interdomain loop interacts via aromatic residue interactions with the upper end of a transmembrane helix and would thereby open the gate.

The Voltage-Dependent Deactivation of the KvAP Channel Involves the Breakage of Its S4 Helix.

Voltage-gated potassium channels (Kv) allow ion permeation upon changes of the membrane electrostatic potential (Vm). Each subunit of these tetrameric channels is composed of six transmembrane helices, of which the anti-parallel helix bundle S1-S4 constitutes the voltage-sensor domain (VSD) and S5-S6 forms the pore domain. Here, using 82 molecular dynamics (MD) simulations involving 266 replicated VSDs, we report novel responses of the archaebacterial potassium channel KvAP to membrane polarization. We show that the S4 α-helix, which is straight in the experimental crystal structure solved under depolarized conditions (Vm ∼ 0), breaks into two segments when the cell membrane is hyperpolarized (Vm << 0), and reversibly forms a single straight helix following depolarization (Vm = 0). The outermost segment of S4 translates along the normal to the membrane, bringing new perspective to previously paradoxical accessibility experiments that were initially thought to imply the displacement of the whole VSD across the membrane. The novel model is applied through steered and unbiased MD simulations to the recently solved whole structure of KvAP. The simulations show that the resting state involves a re-orientation of the S5 α-helix by ∼ 5-6 degrees in respect to the normal of the bilayer, which could result in the constriction and closure of the selectivity filter. Our findings support the idea that the breakage of S4 under (hyper)polarization is a general feature of Kv channels with a non-swapped topology.

Structural and Functional Analysis of Gly212 Mutants Reveals the Importance of Intersubunit Interactions in ASIC1a Channel Function.

Acid-sensing ion channels (ASICs) act as pH sensors in neurons. ASICs contribute to pain sensation, learning, fear behavior and to neuronal death after ischemic stroke. Extracellular acidification induces a transient activation and subsequent desensitization of these Na+-selective channels. ASICs are trimeric channels made of identical or homologous subunits. We have previously shown that mutation of the highly conserved Gly212 residue of human ASIC1a to Asp affects the channel function. Gly212 is located in the proximity of a predicted Cl- binding site at a subunit interface. Here, we have measured the function of a series of Gly212 mutants. We show that substitution of Gly212 affects the ASIC1a pH dependence and current decay kinetics. Intriguingly, the mutations to the acidic residues Asp and Glu have opposing effects on the pH dependence and the current decay kinetics. Analysis of molecular dynamics simulation trajectories started with the coordinates of the closed conformation indicates that the immediate environment of residue 212 in G212E, which shifts the pH dependence to more alkaline values, adopts a conformation closer to the open state. The G212D and G212E mutants have a different pattern of intersubunit salt bridges, that, in the case of G212E, leads to an approaching of neighboring subunits. Based on the comparison of crystal structures, the conformational changes in this zone appear to be smaller during the open-desensitized transition. Nevertheless, MD simulations highlight differences between mutants, suggesting that the changed function upon substitution of residue 212 is due to differences in intra- and intersubunit interactions in its proximity.

Slowing of the Time Course of Acidification Decreases the Acid-Sensing Ion Channel 1a Current Amplitude and Modulates Action Potential Firing in Neurons.

Acid-sensing ion channels (ASICs) are H+-activated neuronal Na+ channels. They are involved in fear behavior, learning, neurodegeneration after ischemic stroke and in pain sensation. ASIC activation has so far been studied only with fast pH changes, although the pH changes associated with many roles of ASICs are slow. It is currently not known whether slow pH changes can open ASICs at all. Here, we investigated to which extent slow pH changes can activate ASIC1a channels and induce action potential signaling. To this end, ASIC1a current amplitudes and charge transport in transfected Chinese hamster ovary cells, and ASIC-mediated action potential signaling in cultured cortical neurons were measured in response to defined pH ramps of 1-40 s duration from pH 7.4 to pH 6.6 or 6.0. A kinetic model of the ASIC1a current was developed and integrated into the Hodgkin-Huxley action potential model. Interestingly, whereas the ASIC1a current amplitude decreased with slower pH ramps, action potential firing was higher upon intermediate than fast acidification in cortical neurons. Indeed, fast pH changes (<4 s) induced short action potential bursts, while pH changes of intermediate speed (4-10 s) induced longer bursts. Slower pH changes (>10 s) did in many experiments not generate action potentials. Computer simulations corroborated these observations. We provide here the first description of ASIC function in response to defined slow pH changes. Our study shows that ASIC1a currents, and neuronal activity induced by ASIC1a currents, strongly depend on the speed of pH changes. Importantly, with pH changes that take >10 s to complete, ASIC1a activation is inefficient. Therefore, it is likely that currently unknown modulatory mechanisms allow ASIC activity in situations such as ischemia and inflammation.

SDHA gain-of-function engages inflammatory mitochondrial retrograde signaling via KEAP1-Nrf2.

Whether screening the metabolic activity of immune cells facilitates discovery of molecular pathology remains unknown. Here we prospectively screened the extracellular acidification rate as a measure of glycolysis and the oxygen consumption rate as a measure of mitochondrial respiration in B cells from patients with primary antibody deficiency. The highest oxygen consumption rate values were detected in three study participants with persistent polyclonal B cell lymphocytosis (PPBL). Exome sequencing identified germline mutations in SDHA, which encodes succinate dehydrogenase subunit A, in all three patients with PPBL. SDHA gain-of-function led to an accumulation of fumarate in PPBL B cells, which engaged the KEAP1-Nrf2 system to drive the transcription of genes encoding inflammatory cytokines. In a single patient trial, blocking the activity of the cytokine interleukin-6 in vivo prevented systemic inflammation and ameliorated clinical disease. Overall, our study has identified pathological mitochondrial retrograde signaling as a disease modifier in primary antibody deficiency.

A Rigorous and Efficient Method To Reweight Very Large Conformational Ensembles Using Average Experimental Data and To Determine Their Relative Information Content.

Flexible polypeptides such as unfolded proteins may access an astronomical number of conformations. The most advanced simulations of such states usually comprise tens of thousands of individual structures. In principle, a comparison of parameters predicted from such ensembles to experimental data provides a measure of their quality. In practice, analyses that go beyond the comparison of unbiased average data have been impossible to carry out on the entirety of such very large ensembles and have, therefore, been restricted to much smaller subensembles and/or nondeterministic algorithms. Here, we show that such very large ensembles, on the order of 10(4) to 10(5) conformations, can be analyzed in full by a maximum entropy fit to experimental average data. Maximizing the entropy of the population weights of individual conformations under experimental χ(2) constraints is a convex optimization problem, which can be solved in a very efficient and robust manner to a unique global solution even for very large ensembles. Since the population weights can be determined reliably, the reweighted full ensemble presents the best model of the combined information from simulation and experiment. Furthermore, since the reduction of entropy due to the experimental constraints is well-defined, its value provides a robust measure of the information content of the experimental data relative to the simulated ensemble and an indication for the density of the sampling of conformational space. The method is applied to the reweighting of a 35,000 frame molecular dynamics trajectory of the nonapeptide EGAAWAASS by extensive NMR (3)J coupling and RDC data. The analysis shows that RDCs provide significantly more information than (3)J couplings and that a discontinuity in the RDC pattern at the central tryptophan is caused by a cluster of helical conformations. Reweighting factors are moderate and consistent with errors in MD force fields of less than 3kT. The required reweighting is larger for an ensemble derived from a statistical coil model, consistent with its coarser nature. We call the method COPER, for convex optimization for ensemble reweighting. Similar advantages of large-scale efficiency and robustness can be obtained for other ensemble analysis methods with convex targets and constraints, such as constrained χ(2) minimization and the maximum occurrence method.

Backbone hydration determines the folding signature of amino acid residues.

The relation between the sequence of a protein and its three-dimensional structure remains largely unknown. A lasting dream is to elucidate the side-chain-dependent driving forces that govern the folding process. Different structural data suggest that aromatic amino acids play a particular role in the stabilization of protein structures. To better understand the underlying mechanism, we studied peptides of the sequence EGAAXAASS (X = Gly, Ile, Tyr, Trp) through comparison of molecular dynamics (MD) trajectories and NMR residual dipolar coupling (RDC) measurements. The RDC data for aromatic substitutions provide evidence for a kink in the peptide backbone. Analysis of the MD simulations shows that the formation of internal hydrogen bonds underlying a helical turn is key to reproduce the experimental RDC values. The simulations further reveal that the driving force leading to such helical-turn conformations arises from the lack of hydration of the peptide chain on either side of the bulky aromatic side chain, which can potentially act as a nucleation point initiating the folding process.

Carbon flux and growth in mature deciduous forest trees exposed to elevated CO2.

Whether rising atmospheric carbon dioxide (CO2) concentrations will cause forests to grow faster and store more carbon is an open question. Using free air CO2 release in combination with a canopy crane, we found an immediate and sustained enhancement of carbon flux through 35-meter-tall temperate forest trees when exposed to elevated CO2. However, there was no overall stimulation in stem growth and leaf litter production after 4 years. Photosynthetic capacity was not reduced, leaf chemistry changes were minor, and tree species differed in their responses. Although growing vigorously, these trees did not accrete more biomass carbon in stems in response to elevated CO2, thus challenging projections of growth responses derived from tests with smaller trees.