A novel twist on molecular interactions between thioredoxin and nicotinamide adenine dinucleotide phosphate‐dependent thioredoxin reductase

The ubiquitous disulfide reductase thioredoxin (Trx) regulates several important biological processes such as seed germination in plants. Oxidized cytosolic Trx is regenerated by nicotinamide adenine dinucleotide phosphate (NADPH)‐dependent thioredoxin reductase (NTR) in a multistep transfer of reducing equivalents from NADPH to Trx via a tightly NTR‐bound flavin. Here, interactions between NTR and Trx are predicted by molecular modelling of the barley NTR:Trx complex (HvNTR2:HvTrxh2) and probed by site directed mutagenesis. Enzyme kinetics analysis reveals mutants in a loop of the flavin adenine dinucleotide (FAD)‐binding domain of HvNTR2 to strongly affect the interaction with Trx. In particular, Trp42 and Met43 play key roles for recognition of the endogenous HvTrxh2. Trx from Arabidopsis thaliana is also efficiently recycled by HvNTR2 but turnover in this case appears to be less dependent on these two residues, suggesting a distinct mode for NTR:Trx recognition. Comparison between the HvNTR2:HvTrxh2 model and the crystal structure of the Escherichia coli NTR:Trx complex reveals major differences in interactions involving the FAD‐ and NADPH‐binding domains as supported by our experiments. Overall, the findings suggest that NTR:Trx interactions in different biological systems are fine‐tuned by multiple intermolecular contacts. Proteins 2014; 82:607–619. © 2013 Wiley Periodicals, Inc.     

A data-driven framework to model the organism–environment system

Organisms modify their development and function in response to the environment. At the same time, the environment is modified by the activities of the organism. Despite the ubiquity of such dynamical interactions in nature, it remains challenging to develop models that accurately represent them, and that can be fitted using data. These features are desirable when modeling phenomena such as phenotypic plasticity, to generate quantitative predictions of how the system will respond to environmental signals of different magnitude or at different times, for example, during ontogeny. Here, we explain a modeling framework that represents the organism and environment as a single coupled dynamical system in terms of inputs and outputs. Inputs are external signals, and outputs are measurements of the system in time. The framework uses time-series data of inputs and outputs to fit a nonlinear black-box model that allows to predict how the system will respond to novel input signals. The framework has three key properties: it captures the dynamical nature of the organism–environment system, it can be fitted with data, and it can be applied without detailed knowledge of the system. We study phenotypic plasticity using in silico experiments and demonstrate that the framework predicts the response to novel environmental signals. The framework allows us to model plasticity as a dynamical property that changes in time during ontogeny, reflecting the well-known fact that organisms are more or less plastic at different developmental stages.     

The potential for signal integration and processing in interacting MAP kinase cascades

The cellular response to environmental stimuli requires biochemical information processing through which sensory inputs and cellular status are integrated and translated into appropriate responses by way of interacting networks of enzymes. One such network, the mitogen-activated protein (MAP) kinase cascade is a highly conserved signal transduction module that propagates signals from cell surface receptors to various cytosolic and nuclear targets by way of a phosphorylation cascade. We have investigated the potential for signal processing within a network of interacting feed-forward kinase cascades typified by the MAP kinase cascade. A genetic algorithm was used to search for sets of kinetic parameters demonstrating representative key input-output patterns of interest. We discuss two of the networks identified in our study, one implementing the exclusive-or function (XOR) and another implementing what we refer to as an in-band detector (IBD) or two-sided threshold. These examples confirm the potential for logic and amplitude-dependent signal processing in interacting MAP kinase cascades demonstrating limited cross-talk. Specifically, the XOR function allows the network to respond to either one, but not both signals simultaneously, while the IBD permits the network to respond exclusively to signals within a given range of strength, and to suppress signals below as well as above this range. The solution to the XOR problem is interesting in that it requires only two interacting pathways, crosstalk at only one layer, and no feedback or explicit inhibition. These types of responses are not only biologically relevant but constitute signal processing modules that can be combined to create other logical functions and that, in contrast to amplification, cannot be achieved with a single cascade or with two non-interacting cascades. Our computational results revealed surprising similarities between experimental data describing the JNK/MKK4/MKK7 pathway and the solution for the IBD that evolved from the genetic algorithm. The evolved IBD not only exhibited the required non-monotonic signal strength-response, but also demonstrated transient and sustained responses that properly reflected the input signal strength, dependence on both of the MAPKKs for signaling, phosphorylation site preferences by each of the MAPKKs, and both activation and inhibition resulting from the overexpression of one of the MAPKKs.     

European Roe Deer Increase Vigilance When Faced with Immediate Predation Risk by Eurasian Lynx

Prey animals must balance antipredatory behaviour with foraging behaviour. According to the risk allocation hypothesis, prey increase antipredatory behaviour and reduce foraging activity during pulses of high risk, but with continuous risk, other activities must continue and antipredatory behaviour decreases despite the risk. We studied the impact of lynx presence on the vigilance behaviour of wild roe deer under conditions of (i) a pulsed elevated risk by experimentally spreading lynx urine as an olfactory cue, and (ii) continuous risk by comparing an area where lynx was eradicated 160 years ago to an area where lynx has been re‐introduced 30 years ago. Roe deer were extremely vigilant in response to the predator olfactory cue; however, roe deer vigilance did not differ measurably among areas with or without potential lynx presence. Deer were more vigilant before sunset than during the night at both study areas, probably due to long‐term adaptation of roe deer to human hunting during daytime. Vigilance decreased from August to September even though activity of lynx increases in autumn, which may be a result either of increased foraging due to decrease in food quality in autumn, or of changes in social organization of the deer. Our results suggested that the degree of vigilance depends on environmental cues. We found that roe deer respond to lynx urine despite a long absence of lynx in the ecosystem. Our results support the risk allocation hypothesis for responses to pulses of high risk but not for responses to continuous elevated levels of risk.     

Drug block of the hERG potassium channel: insight from modeling

Many commonly used, structurally diverse, drugs block the human ether-a-go-go-related gene (hERG) K(+) channel to cause acquired long QT syndrome, which can lead to sudden death via lethal cardiac arrhythmias. This undesirable side effect is a major hurdle in the development of safe drugs. To gain insight about the structure of hERG and the nature of drug block we have produced structural models of the channel pore domain, into each of which we have docked a set of 20 hERG blockers. In the absence of an experimentally determined three-dimensional structure of hERG, each of the models was validated against site-directed mutagenesis data. First, hERG models were produced of the open and closed channel states, based on homology with the prokaryotic K(+) channel crystal structures. The modeled complexes were in partial agreement with the mutagenesis data. To improve agreement with mutagenesis data, a KcsA-based model was refined by rotating the four copies of the S6 transmembrane helix half a residue position toward the C-terminus, so as to place all residues known to be involved in drug binding in positions lining the central cavity. This model produces complexes that are consistent with mutagenesis data for smaller, but not larger, ligands. Larger ligands could be accommodated following refinement of this model by enlarging the cavity using the inherent flexibility about the glycine hinge (Gly648) in S6, to produce results consistent with the experimental data for the majority of ligands tested.     

A subset of chemosensory genes differs between two populations of a specialized leaf beetle after host plant shift

Due to its fundamental role in shaping host selection behavior, we have analyzed the chemosensory repertoire of Chrysomela lapponica. This specialized leaf beetle evolved distinct populations which shifted from the ancestral host plant, willow (Salix sp., Salicaceae), to birch (Betula rotundifolia, Betulaceae). We identified 114 chemosensory candidate genes in adult C. lapponica: 41 olfactory receptors (ORs), eight gustatory receptors, 17 ionotropic receptors, four sensory neuron membrane proteins, 32 odorant binding proteins (OBPs), and 12 chemosensory proteins (CSP) by RNA‐seq. Differential expression analyses in the antennae revealed significant upregulation of one minus‐C OBP (ClapOBP27) and one CSP (ClapCSP12) in the willow feeders. In contrast, one OR (ClapOR17), four minus‐C OBPs (ClapOBP02, 07, 13, 20), and one plus‐C OBP (ClapOBP32) were significantly upregulated in birch feeders. The differential expression pattern in the legs was more complex. To narrow down putative ligands acting as cues for host discrimination, the relative abundance and diversity of volatiles of the two host plant species were analyzed. In addition to salicylaldehyde (willow‐specific), both plant species differed mainly in their emission rate of terpenoids such as (E,E)‐α‐farnesene (high in willow) or 4,8‐dimethylnona‐1,3,7‐triene (high in birch). Qualitatively, the volatiles were similar between willow and birch leaves constituting an “olfactory bridge” for the beetles. Subsequent structural modeling of the three most differentially expressed OBPs and docking studies using 22 host volatiles indicated that ligands bind with varying affinity. We suggest that the evolution of particularly minus‐C OBPs and ORs in C. lapponica facilitated its host plant shift via chemosensation of the phytochemicals from birch as novel host plant.