Transcending Markov: non-Markovian rate processes of thermosensitive TRP ion channels.

The Markov state model (MSM) is a popular theoretical tool for describing the hierarchy of time scales involved in the function of many proteins especially ion channel gating. An MSM is a particular case of the general non-Markovian model, where the rate of transition from one state to another does not depend on the history of state occupancy within the system, i.e. it only includes reversible, non-dissipative processes. However, an MSM requires knowledge of the precise conformational state of the protein and is not predictive when those details are not known. In the case of ion channels, this simple description fails in real (non-equilibrium) situations, for example when local temperature changes, or when energy losses occur during channel gating. Here, we show it is possible to use non-Markovian equations (i.e. offer a general description that includes the MSM as a particular case) to develop a relatively simple analytical model that describes the non-equilibrium behaviour of the temperature-sensitive transient receptor potential (TRP) ion channels, TRPV1 and TRPM8. This model accurately predicts asymmetrical opening and closing rates, infinite processes and the creation of new states, as well as the effect of temperature changes throughout the process. This approach therefore overcomes the limitations of the MSM and allows us to go beyond a mere phenomenological description of the dynamics of ion channel gating towards a better understanding of the physics underlying these processes.

Quantum mechanics in ion channel.

Many attempts have been made to associate physical systems with system biology. Most of them are showing very impressive analogs. However, in quantum mechanics, very limited data and knowledge are available regarding analogs. The ion movement is approximately constant in velocity since the time periods are so small and are negligible. In this study, I focused on a new aspect of ion channels using quantum mechanics. Suppose that the time required for the ion to cross a distance between two gates.

Metabolome dynamics during wheat domestication.

One of the most important crops worldwide is wheat. Wheat domestication took place about 10,000 years ago. Not only that its wild progenitors have been discovered and phenotypically characterized, but their genomes were also sequenced and compared to modern wheat. While comparative genomics is essential to track genes that contribute to improvement in crop yield, comparative analyses of functional biological end-products, such as metabolites, are still lacking. With the advent of rigorous mass-spectrometry technologies, it is now possible to address that problem on a big-data scale. In attempt to reveal classes of metabolites, which are associated with wheat domestication, we analyzed the metabolomes of wheat kernel samples from various wheat lines. These wheat lines represented subspecies of tetraploid wheat along primary and secondary domestications, including wild emmer, domesticated emmer, landraces durum, and modern durum. We detected that the groups of plant metabolites such as plant-defense metabolites, antioxidants and plant hormones underwent significant changes during wheat domestication. Our data suggest that these metabolites may have contributed to the improvement in the agricultural fitness of wheat. Closer evaluation of specific metabolic pathways may result in the future in genetically-engineered high-yield crops.

Changes in "natural antibiotic" metabolite composition during tetraploid wheat domestication.

Gramineous plants protect their seeds from a variety of biotic stresses by producing toxic and deterrent secondary metabolites such as benzoxazinoids. It is unclear how the composition and abundance of these natural toxins has changed over the course of crop-plant domestication. To address this uncertainty, we characterized differences in metabolic levels of benzoxazinoids and their derivatives, between four lines of tetraploid wheat: wild emmer wheat (WEW), the direct progenitor of modern wheat; non-fragile domesticated emmer wheat (DEW), which was first domesticated about 11,000 years ago; the subsequently developed non-fragile and free-threshing durum landraces (LD); and modern durum (MD) varieties. Three-dimensional principal component analysis of mass spectrometry data of wheat metabolites showed with high resolution clear differences between metabolic profiles of WEW, DEW, and durum (LD + MD) and similarity in the metabolic profiles of the two durum lines (LD and MD) that is coherent with the phylogenetic relationship between the corresponding wheat lines. Moreover, our results indicated that some secondary metabolites involved in plant defense mechanisms became significantly more abundant during wheat domestication, while other defensive metabolites decreased or were lost. These metabolic changes reflect the beneficial or detrimental roles the corresponding metabolites might play during the domestication of three taxonomic subspecies of tetraploid wheat (Triticum turgidum).

Neuromodulation of Astrocytic K+ Clearance.

Potassium homeostasis is fundamental for brain function. Therefore, effective removal of excessive K+ from the synaptic cleft during neuronal activity is paramount. Astrocytes play a key role in K+ clearance from the extracellular milieu using various mechanisms, including uptake via Kir channels and the Na+-K+ ATPase, and spatial buffering through the astrocytic gap-junction coupled network. Recently we showed that alterations in the concentrations of extracellular potassium ([K+]o) or impairments of the astrocytic clearance mechanism affect the resonance and oscillatory behavior of both the individual and networks of neurons. These results indicate that astrocytes have the potential to modulate neuronal network activity, however, the cellular effectors that may affect the astrocytic K+ clearance process are still unknown. In this study, we have investigated the impact of neuromodulators, which are known to mediate changes in network oscillatory behavior, on the astrocytic clearance process. Our results suggest that while some neuromodulators (5-HT; NA) might affect astrocytic spatial buffering via gap-junctions, others (DA; Histamine) primarily affect the uptake mechanism via Kir channels. These results suggest that neuromodulators can affect network oscillatory activity through parallel activation of both neurons and astrocytes, establishing a synergistic mechanism to maximize the synchronous network activity.

The dynamics of K+ channel gates as a biological transistor.

Potassium channels are pore-forming membrane proteins that open and close in response to changes in a chemical or electrical potential, thereby regulating the flow of potassium ions across biological membranes. Two regions of the same channels are acting in tandem and enable ion flow through the channel pore. I refer to this coupled action as a "gate linker". To closely examine the role of the gate linker in the channel function, I mutated the amino acids in the cDNA of this region, and used from knowen mutaion, either alone or together with the amino acids of adjacent regions. I have emphasized the importance of the linker between these two gates - mutations in this region may cause conformational changes that play a fundamental role in mediating the coupling between the voltage sensor, activation gate and selectivity filter elements of Kv channels. I observe that free energy considerations show the significance of the coupling between the activation and inactivation gates. Moreover, a symmetry between the coupling and sensor spring strength leads to the destruction of ion conductivity. I present a thermodynamic framework for the possible study of multiple channel blocks. The arising physical perspective of the gating process gives rise to new research avenues of the coupling mode of potassium channels and may assist in explaining the centrality of the "gate linker" to the channel function.

Adaptive evolution of benzoxazinoids in wild emmer wheat, Triticum dicoccoides, at "Evolution Canyon", Mount Carmel, Israel.

BACKGROUND: "Evolution Canyon" (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unraveling evolution-in-action, highlighting the evolutionary processes of biodiversity evolution, adaptation, and incipient sympatric speciation. A major model organism in ECI is the tetraploid wild emmer wheat, Triticum dicoccoides (TD), the progenitor of cultivated emmer and durum wheat. TD displays dramatic interslope adaptive evolutionary divergence on the tropical, savannoid-hot and dry south-facing, "African" slope (AS), and on the temperate, forested, cool and humid, north-facing, "European" slope (ES), separated on average by 250 m. From the perspective of chemical evolution and metabolomics, it is important to unravel interslope divergence in biologically relevant secondary metabolites between the abutting slope populations. Here, in TD we examined hydroxamic acid (Hx), which is a family of secondary cereal metabolites, and plays a major role in defending the plant against fungi, insects and weeds. RESULTS: Our examination revealed that higher concentrations of DIBOA and DIMBOA were found in seedlings growing in the same greenhouse from seeds collected from the cool and humid forested ES, whereas the seedlings of seeds collected from the savannoid AS (both in root and shoot tissues), showed no DIMBOA. Remarkably, only DIBOA appears in both shoots and roots of the AS seedlings. It rises to a peak and then decreases in both organs and in seedlings from both slopes. The DIMBOA, which appears only in the ES seedlings, rises to a peak and decreases in the shoot, but increased and remained in a plateau in the root, till the end of the experiment. CONCULSIONS/SIGNIFICANCE: The results suggest stronger genetic resistance of defense compounds DIBOA and DIMBOA against biotic stresses (fungi and other pathogens) by ES seedlings. However, AS seedlings responded earlier but were to the same biotic stresses. The genetic difference found in AS seedlings was caused by the main adaptive selection in AS, which was against climatic, abiotic stresses, and was weaker, or not at all, against biotic stresses. The distinct genetic interslope differences appear important and is very significant and are elaborated in the discussion.

Adaptive methylation regulation of p53 pathway in sympatric speciation of blind mole rats, Spalax.

Epigenetic modifications play significant roles in adaptive evolution. The tumor suppressor p53, well known for controlling cell fate and maintaining genomic stability, is much less known as a master gene in environmental adaptation involving methylation modifications. The blind subterranean mole rat Spalax eherenbergi superspecies in Israel consists of four species that speciated peripatrically. Remarkably, the northern Galilee species Spalax galili (2n = 52) underwent adaptive ecological sympatric speciation, caused by the sharply divergent chalk and basalt ecologies. This was demonstrated by mitochondrial and nuclear genomic evidence. Here we show that the expression patterns of the p53 regulatory pathway diversified between the abutting sympatric populations of S. galili in sharply divergent chalk-basalt ecologies. We identified higher methylation on several sites of the p53 promoter in the population living in chalk soil (chalk population). Site mutagenesis showed that methylation on these sites linked to the transcriptional repression of p53 involving Cut-Like Homeobox 1 (Cux1), paired box 4 (Pax 4), Pax 6, and activator protein 1 (AP-1). Diverse expression levels of p53 between the incipiently sympatrically speciating chalk-basalt abutting populations of S. galili selectively affected cell-cycle arrest but not apoptosis. We hypothesize that methylation modification of p53 has adaptively shifted in supervising its target genes during sympatric speciation of S. galili to cope with the contrasting environmental stresses of the abutting divergent chalk-basalt ecologies.

Putative adaptive inter-slope divergence of transposon frequency in fruit flies (Drosophila melanogaster) at "Evolution Canyon", Mount Carmel, Israel.

BACKGROUND: The current analysis of transposon elements (TE) in Drosophila melanogaster at Evolution Canyon, (EC), Israel, is based on data and analysis done by our collaborators (Drs. J. Gonzalez, J. Martinez and W. Makalowski, this issue). They estimated the frequencies of 28 TEs (transposon elements) in fruit flies (D. melanogaster) from the ecologically tropic, hot, and dry south-facing slope (SFS) or "African" slope (AS) of EC and compared it with the TE frequencies on the temperate-cool and humid north-facing slope (NFS) or "European" slope (ES), separated, on average, by 250 m. The flies were sampled from two stations on each slope. We received their results, including the frequencies of each TE on each slope, and the probabilities of the statistical analyses (G-tests) of each TE separately. We continued the analysis of the inter-slope differences of the frequencies of the TEs, and based our different conclusions on that analysis and on the difference between micro (=EC) and macro (2000 km.) comparisons [Gonzalez et al. 2015 doi: 10.1186/s13062-015-0075-4 ]. RESULTS: Our collaborators based all their conclusions on the non-significant results of each of the individual tests of the 28 TEs. We analysed also the distribution of the TE differences between the slopes, based on their results. Thirteen TEs were more frequent on the SFS, 11 were more frequent on the NFS, and four had equal frequencies. Because of the equalizing effect of the ongoing migration, only small and temporary differences between the slopes (0 - 0.06) were regarded by us as random fluctuations (drift). Three TEs were intermediate (0.08-0.09) and await additional research. The 11 TEs with large frequency differences (0.12 - 0.22) were regarded by us as putative adaptive TEs, because the equalizing power of ongoing migration will eliminate random large differences. Five of them were higher on the SFS and six were higher on the NFS. Gaps in the distribution of the differences distinguished between the large and small differences. The large gap among the 11 TEs favored on the NFS was significant and supports our rejection of drift as the only explanation of the distribution of the slope differences. The gaps in the distribution of the differences separated the putative TEs with strong enough selection from those TEs that couldn't overrule the migration. The results are compared and contrasted with the directional effect of the frequencies of the same TEs in the study of global climatic comparisons across thousands of kilometers. From the 11 putative adaptive TEs in the local "Evolution Canyon," six differentiate in the same direction as in the continental comparisons and four in the opposite direction. One TE, FBti0019144, differentiated in EC in the same direction as in Australia and in the opposite direction to that of North America. CONCLUSIONS: We presume that the major divergent evolutionary driving force at the local EC microsite is natural selection overruling gene flow. Therefore, after we rejected drift as an explanation of all the large slope differences, we regarded them as putatively adaptive. In order to substantiate the individual TE adaptation, we need to increase the sample sizes and reveal the significant adaptive TEs. The comparison of local and global studies show only partial similarity in the adaptation of the TEs, because of the dryness of the ecologically tropical climate in EC, in contrast to the wet tropical climate in the global compared climates. Moreover, adaptation of a TE may be expressed only in part of the time and specific localities.

Natural selection causes adaptive genetic resistance in wild emmer wheat against powdery mildew at "Evolution Canyon" microsite, Mt. Carmel, Israel.

BACKGROUND: "Evolution Canyon" (ECI) at Lower Nahal Oren, Mount Carmel, Israel, is an optimal natural microscale model for unraveling evolution in action highlighting the basic evolutionary processes of adaptation and speciation. A major model organism in ECI is wild emmer, Triticum dicoccoides, the progenitor of cultivated wheat, which displays dramatic interslope adaptive and speciational divergence on the tropical-xeric "African" slope (AS) and the temperate-mesic "European" slope (ES), separated on average by 250 m. METHODS: We examined 278 single sequence repeats (SSRs) and the phenotype diversity of the resistance to powdery mildew between the opposite slopes. Furthermore, 18 phenotypes on the AS and 20 phenotypes on the ES, were inoculated by both Bgt E09 and a mixture of powdery mildew races. RESULTS: In the experiment of genetic diversity, very little polymorphism was identified intra-slope in the accessions from both the AS or ES. By contrast, 148 pairs of SSR primers (53.23%) amplified polymorphic products between the phenotypes of AS and ES. There are some differences between the two wild emmer wheat genomes and the inter-slope SSR polymorphic products between genome A and B. Interestingly, all wild emmer types growing on the south-facing slope (SFS=AS) were susceptible to a composite of Blumeria graminis, while the ones growing on the north-facing slope (NFS=ES) were highly resistant to Blumeria graminis at both seedling and adult stages. CONCLUSION/SIGNIFICANCE: Remarkable inter-slope evolutionary divergent processes occur in wild emmer wheat, T. dicoccoides at EC I, despite the shot average distance of 250 meters. The AS, a dry and hot slope, did not develop resistance to powdery mildew, whereas the ES, a cool and humid slope, did develop resistance since the disease stress was strong there. This is a remarkable demonstration in host-pathogen interaction on how resistance develops when stress causes an adaptive result at a micro-scale distance.

Probing the transition state of the allosteric pathway of the Shaker Kv channel pore by linear free-energy relations.

Long-range coupling between distant functional elements of proteins may rely on allosteric communication trajectories lying along the protein structure, as described in the case of the Shaker voltage-activated potassium (Kv) channel model allosteric system. Communication between the distant Kv channel activation and slow inactivation pore gates was suggested to be mediated by a network of local pairwise and higher-order interactions among the functionally unique residues that constitute the allosteric trajectory. The mechanism by which conformational changes propagate along the Kv channel allosteric trajectory to achieve pore opening, however, remains unclear. Such conformational changes may propagate in either a concerted or a sequential manner during the reaction coordinate of channel opening. Residue-level structural information on the transition state of channel gating is required to discriminate between these possibilities. Here, we combine patch-clamp electrophysiology recordings of Kv channel gating and analysis using linear free-energy relations, focusing on a select set of residues spanning the allosteric trajectory of the Kv channel pore. We show that all allosteric trajectory residues tested exhibit an open-like conformation in the transition state of channel opening, implying that coupling interactions occur along the trajectory break in a concerted manner upon moving from the closed to the open state. Energetic coupling between the Kv channel gates thus occurs in a concerted fashion in both the spatial and the temporal dimensions, strengthening the notion that such trajectories correspond to pathways of mechanical deformation along which conformational changes propagate.

Gating the pore of potassium leak channels.

A key feature of potassium channel function is the ability to switch between conducting and non-conducting states by undergoing conformational changes in response to cellular or extracellular signals. Such switching is facilitated by the mechanical coupling of gating domain movements to pore opening and closing. Two-pore domain potassium channels (K(2P)) conduct leak or background potassium-selective currents that are mostly time- and voltage-independent. These channels play a significant role in setting the cell resting membrane potential and, therefore modulate cell responsiveness and excitability. Thus, K(2P) channels are key players in numerous physiological processes and were recently shown to also be involved in human pathologies. It is well established that K(2P) channel conductance, open probability and cell surface expression are significantly modulated by various physical and chemical stimuli. However, in understanding how such signals are translated into conformational changes that open or close the channels gate, there remain more open questions than answers. A growing line of evidence suggests that the outer pore area assumes a critical role in gating K(2P) channels, in a manner reminiscent of C-type inactivation of voltage-gated potassium channels. In some K(2P) channels, this gating mechanism is facilitated in response to external pH levels. Recently, it was suggested that K(2P) channels also possess a lower activation gate that is positively coupled to the outer pore gate. The purpose of this review is to present an up-to-date summary of research describing the conformational changes and gating events that take place at the K(2P) channel ion-conducting pathway during the channel regulation.

Inverse coupling in leak and voltage-activated K+ channel gates underlies distinct roles in electrical signaling.

Voltage-activated (Kv) and leak (K(2P)) K(+) channels have key, yet distinct, roles in electrical signaling in the nervous system. Here we examine how differences in the operation of the activation and slow inactivation pore gates of Kv and K(2P) channels underlie their unique roles in electrical signaling. We report that (i) leak K(+) channels possess a lower activation gate, (ii) the activation gate is an important determinant controlling the conformational stability of the K(+) channel pore, (iii) the lower activation and upper slow inactivation gates of leak channels cross-talk and (iv) unlike Kv channels, where the two gates are negatively coupled, these two gates are positively coupled in K(2P) channels. Our results demonstrate how basic thermodynamic properties of the K(+) channel pore, particularly conformational stability and coupling between gates, underlie the specialized roles of Kv and K(2P) channel families in electrical signaling.

A novel mechanism for human K2P2.1 channel gating. Facilitation of C-type gating by protonation of extracellular histidine residues.

The mammalian K2P2.1 potassium channel (TREK-1, KCNK2) is highly expressed in excitable tissues, where it plays a key role in the cellular mechanisms of neuroprotection, anesthesia, pain perception, and depression. Here, we report that external acidification, within the physiological range, strongly inhibits the human K2P2.1 channel by inducing "C-type" closure. We have identified two histidine residues (i.e. His-87 and His-141), located in the first external loop of the channel, that govern the response of the channel to external pH. We demonstrate that these residues are within physical proximity to glutamate 84, homologous to Shaker Glu-418, KcsA Glu-51, and KCNK0 Glu-28 residues, all previously argued to stabilize the outer pore gate in the open conformation by forming hydrogen bonds with pore-adjacent residues. We thus propose a novel mechanism for pH sensing in which protonation of His-141 and His-87 generates a local positive charge that serves to draw Glu-84 away from its natural interactions, facilitating the collapse of the selectivity filter region. In accordance with this proposed mechanism, low pH modified K2P2.1 selectivity toward potassium. Moreover, the proton-mediated effect was inhibited by external potassium ions and was enhanced by a mutation (S164Y) known to accelerate C-type gating. Furthermore, proton-induced current inhibition was more pronounced at negative potentials. Thus, voltage-dependent C-type gating acceleration by protons represents a novel mechanism for K2P2.1 outward rectification.