A dynamic model of the hypoxia-inducible factor 1α (HIF-1α) network

Activation of the hypoxia-inducible factor (HIF) pathway is a critical step in the transcriptional response to hypoxia. Although many of the key proteins involved have been characterised, the dynamics of their interactions in generating this response remain unclear. In the present study, we have generated a comprehensive mathematical model of the HIF-1a pathway based on core validated components and dynamic experimental data, and confirm the previously described connections within the predicted network topology. Our model confirms previous work demonstrating that the steps leading to optimal HIF-1a transcriptional activity require sequential inhibition of both prolyl- and asparaginyl-hydroxylases. We predict from our model (and confirm experimentally) that there is residual activity of the asparaginyl-hydroxylase FIH (factor inhibiting HIF) at low oxygen tension. Furthermore, silencing FIH under conditions where prolyl-hydroxylases are inhibited results in increased HIF-1a transcriptional activity, but paradoxically decreases HIF-1a stability. Using a core module of the HIF network and mathematical proof supported by experimental data, we propose that asparaginyl hydroxylation confers a degree of resistance upon HIF-1a to proteosomal degradation. Thus, through in vitro experimental data and in silico predictions, we provide a comprehensive model of the dynamic regulation of HIF-1a transcriptional activity by hydroxylases and use its predictive and adaptive properties to explain counter-intuitive biological observations.

Hypoxia-inducible factor (HIF) network:insights from mathematical models

Oxygen is a crucial molecule for cellular function. When oxygen demand exceeds supply, the oxygen sensing pathway centred on the hypoxia inducible factor (HIF) is switched on and promotes adaptation to hypoxia by up-regulating genes involved in angiogenesis, erythropoiesis and glycolysis. The regulation of HIF is tightly modulated through intricate regulatory mechanisms. Notably, its protein stability is controlled by the oxygen sensing prolyl hydroxylase domain (PHD) enzymes and its transcriptional activity is controlled by the asparaginyl hydroxylase FIH (factor inhibiting HIF-1).To probe the complexity of hypoxia-induced HIF signalling, efforts in mathematical modelling of the pathway have been underway for around a decade. In this paper, we review the existing mathematical models developed to describe and explain specific behaviours of the HIF pathway and how they have contributed new insights into our understanding of the network. Topics for modelling included the switch-like response to decreased oxygen gradient, the role of micro environmental factors, the regulation by FIH and the temporal dynamics of the HIF response. We will also discuss the technical aspects, extent and limitations of these models. Recently, HIF pathway has been implicated in other disease contexts such as hypoxic inflammation and cancer through crosstalking with pathways like NF?B and mTOR. We will examine how future mathematical modelling and simulation of interlinked networks can aid in understanding HIF behaviour in complex pathophysiological situations. Ultimately this would allow the identification of new pharmacological targets in different disease settings.

Cortical dynamics and synchronization related to multiple target consolidation under rapid-serial-visual-presentation conditions

The present report reviews behavioural, electroencephalographic, and especially magnetoencephalographic findings on the cortical mechanisms underlying attentional processes that separate targets from distractors and that ensure durable target representations for goal-directed action. A common way of investigation is to observe the system’s overt and covert behaviour when capacity limitations are reached. Here we focus on the aspect of temporally enhanced processing load, namely on performance deficits occurring under rapid-serial-visual-presentation (RSVP) conditions. The most prominent of these deficits is the so-called “attentional blink” (AB) effect. We first report MEG findings with respect to the time course of activation that shows modulations around 300 ms after target onset which reflect demands and success of target consolidation. Then, findings regarding long-range inter-area phase synchronization are reported that are hypothesized to mediate communication within the attentional network. Changes in synchronization reflect changes in the attentional demands of the task and are directly related to behavioural performance. Furthermore, enhanced vigilance of the system elicits systematically increased synchronization indices. A hypothetical framework is sketched out that aims at explaining limitations in multiple target consolidation under RSVP conditions.

Novel neuronal and astrocytic mechanisms in thalamocortical loop dynamics

In this review, we summarize three sets of findings that have recently been observed in thalamic astrocytes and neurons, and discuss their significance for thalamocortical loop dynamics. (i) A physiologically relevant ‘window’ component of the low–voltage–activated, T–type Ca2+ current (ITwindow) plays an essential part in the slow (less than 1 Hz) sleep oscillation in adult thalamocortical (TC) neurons, indicating that the expression of this fundamental sleep rhythm in these neurons is not a simple reflection of cortical network activity. It is also likely that ITwindow underlies one of the cellular mechanisms enabling TC neurons to produce burst firing in response to novel sensory stimuli. (ii) Both electrophysiological and dye–injection experiments support the existence of gap junction–mediated coupling among young and adult TC neurons. This finding indicates that electrical coupling–mediated synchronization might be implicated in the high and low frequency oscillatory activities expressed by this type of thalamic neuron. (iii) Spontaneous intracellular Ca2+ ([Ca2+]i) waves propagating among thalamic astrocytes are able to elicit large and long–lasting N–methyl–D–aspartate–mediated currents in TC neurons. The peculiar developmental profile within the first two postnatal weeks of these astrocytic [Ca2+]i transients and the selective activation of these glutamate receptors point to a role for this astrocyte–to–neuron signalling mechanism in the topographic wiring of the thalamocortical loop. As some of these novel cellular and intracellular properties are not restricted to thalamic astrocytes and neurons, their significance may well apply to (patho)physiological functions of glial and neuronal elements in other brain areas.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

Indirect ties in knowledge networks:a social network analysis with ordered weighted averaging operators

This PhD thesis analyses networks of knowledge flows, focusing on the role of indirect ties in the knowledge transfer, knowledge accumulation and knowledge creation process. It extends and improves existing methods for mapping networks of knowledge flows in two different applications and contributes to two stream of research. To support the underlying idea of this thesis, which is finding an alternative method to rank indirect network ties to shed a new light on the dynamics of knowledge transfer, we apply Ordered Weighted Averaging (OWA) to two different network contexts. Knowledge flows in patent citation networks and a company supply chain network are analysed using Social Network Analysis (SNA) and the OWA operator. The OWA is used here for the first time (i) to rank indirect citations in patent networks, providing new insight into their role in transferring knowledge among network nodes; and to analyse a long chain of patent generations along 13 years; (ii) to rank indirect relations in a company supply chain network, to shed light on the role of indirectly connected individuals involved in the knowledge transfer and creation processes and to contribute to the literature on knowledge management in a supply chain. In doing so, indirect ties are measured and their role as means of knowledge transfer is shown. Thus, this thesis represents a first attempt to bridge the OWA and SNA fields and to show that the two methods can be used together to enrich the understanding of the role of indirectly connected nodes in a network. More specifically, the OWA scores enrich our understanding of knowledge evolution over time within complex networks. Future research can show the usefulness of OWA operator in different complex networks, such as the on-line social networks that consists of thousand of nodes.

Cortical oscillatory dynamics and benzodiazepine-site modulation of tonic inhibition in fast spiking interneurons

Tonic conductance mediated by extrasynaptic GABAA receptors has been implicated in the modulation of network oscillatory activity. Using an in vitro brain slice to produce oscillatory activity and a kinetic model of GABAA receptor dynamics, we show that changes in tonic inhibitory input to fast spiking interneurons underlie benzodiazepine-site mediated modulation of neuronal network synchrony in rat primary motor cortex. We found that low concentrations (10 nM) of the benzodiazepine site agonist, zolpidem, reduced the power of pharmacologically-induced beta-frequency (15–30 Hz) oscillatory activity. By contrast, higher doses augmented beta power. Application of the antagonist, flumazenil, also increased beta power suggesting endogenous modulation of the benzodiazepine binding site. Voltage-clamp experiments revealed that pharmacologically-induced rhythmic inhibitory postsynaptic currents were reduced by 10 nM zolpidem, suggesting an action on inhibitory interneurons. Further voltage -clamp studies of fast spiking cells showed that 10 nM zolpidem augmented a tonic inhibitory GABAA receptor mediated current in fast spiking cells whilst higher concentrations of zolpidem reduced the tonic current. A kinetic model of zolpidem-sensitive GABAA receptors suggested that incubation with 10 nM zolpidem resulted in a high proportion of GABAA receptors locked in a kinetically slow desensitized state whilst 30 nM zolpidem favoured rapid transition into and out of desensitized states. This was confirmed experimentally using a challenge with saturating concentrations of GABA. Selective modulation of an interneuron-specific tonic current may underlie the reversal of cognitive and motor deficits afforded by low-dose zolpidem in neuropathological states.

Anisotropic opinion dynamics

We consider the process of opinion formation in a society of interacting agents, where there is a set B of socially accepted rules. In this scenario, we observed that agents, represented by simple feed-forward, adaptive neural networks, may have a conservative attitude (mostly in agreement with B) or liberal attitude (mostly in agreement with neighboring agents) depending on how much their opinions are influenced by their peers. The topology of the network representing the interaction of the society's members is determined by a graph, where the agents' properties are defined over the vertexes and the interagent interactions are defined over the bonds. The adaptability of the agents allows us to model the formation of opinions as an online learning process, where agents learn continuously as new information becomes available to the whole society (online learning). Through the application of statistical mechanics techniques we deduced a set of differential equations describing the dynamics of the system. We observed that by slowly varying the average peer influence in such a way that the agents attitude changes from conservative to liberal and back, the average social opinion develops a hysteresis cycle. Such hysteretic behavior disappears when the variance of the social influence distribution is large enough. In all the cases studied, the change from conservative to liberal behavior is characterized by the emergence of conservative clusters, i.e., a closed knitted set of society members that follow a leader who agrees with the social status quo when the rule B is challenged.

Conformational dynamics associated with ligand binding to vertebrate hexa-coordinate hemoglobins

Neuroglobin (Ngb) and cytoglobin (Cygb) are two new additions to the globin family, exhibiting heme iron hexa-coordination, a disulfide bond and large internal cavities. These proteins are implicated in cytoprotection under hypoxic-ischemic conditions, but the molecular basis of their cytoprotective function is unclear. Herein, a photothermal and spectroscopic study of the interactions of diatomic ligands with Ngb, Cygb, myoglobin and hemoglobin is presented. The impact of the disulfide bond in Ngb and Cygb and role of conserved residues in Ngb His64, Val68, Cys55, Cys120 and Tyr44 on conformational dynamics associated with ligand binding/dissociation were investigated. Transient absorption and photoacoustic calorimetry studies indicate that CO photo-dissociation from Ngb leads to a volume expansion (13.4±0.9 mL mol-1), whereas a smaller volume change was determined for Ngb with reduced Cys (ΔV=4.6±0.3 mL mol-1). Furthermore, Val68 side chain regulates ligand migration between the distal pocket and internal hydrophobic cavities since Val68Phe geminate quantum yield is ∼2.7 times larger than that of WT Ngb. His64Gln and Tyr44Phe mutations alter the thermodynamic parameters associated with CO photo-release indicating that electrostatic/hydrogen binding network that includes heme propionate groups, Lys 67, His64, and Tyr 44 in Ngb modulates the energetics of CO photo-dissociation. In Cygb, CO escape from the protein matrix is fast (< 40 ns) with a ΔH of 18±2 kcal mol-1 in Cygbred, whereas disulfide bridge formation promotes a biphasic ligand escape associated with an overall enthalpy change of 9±4 kcal mol-1. Therefore, the disulfide bond modulates conformational dynamics in Ngb and Cygb. I propose that in Cygb with reduced Cys the photo-dissociated ligand escapes through the hydrophobic tunnel as occurs in Ngb, whereas the CO preferentially migrates through the His64 gate in Cygbox. To characterize Cygb surface 1,8-ANS interactions with Cygb were investigated employing fluorescence spectroscopy, ITC and docking simulations. Two 1,8-ANS binding sites were identified. One binding site is located close to the extended N-terminus of Cygb and was also identified as a binding site for oleate. Furthermore, guanidinium hydrochloride-induced unfolding studies of Cygb reveal that the disulfide bond does not impact Cygb stability, whereas binding of cyanide slightly increases the protein stability.

Basic measurements of radiation from the Baseline Surface Radiation Network (BSRN) in the years 1992 to 2012, reference list of 6378 datasets

In the framework of the global energy balance, the radiative energy exchanges between Sun, Earth and space are now accurately quantified from new satellite missions. Much less is known about the magnitude of the energy flows within the climate system and at the Earth surface, which cannot be directly measured by satellites. In addition to satellite observations, here we make extensive use of the growing number of surface observations to constrain the global energy balance not only from space, but also from the surface. We combine these observations with the latest modeling efforts performed for the 5th IPCC assessment report to infer best estimates for the global mean surface radiative components. Our analyses favor global mean downward surface solar and thermal radiation values near 185 and 342 Wm**-2, respectively, which are most compatible with surface observations. Combined with an estimated surface absorbed solar radiation and thermal emission of 161 Wm**-2 and 397 Wm**-2, respectively, this leaves 106 Wm**-2 of surface net radiation available for distribution amongst the non-radiative surface energy balance components. The climate models overestimate the downward solar and underestimate the downward thermal radiation, thereby simulating nevertheless an adequate global mean surface net radiation by error compensation. This also suggests that, globally, the simulated surface sensible and latent heat fluxes, around 20 and 85 Wm**-2 on average, state realistic values. The findings of this study are compiled into a new global energy balance diagram, which may be able to reconcile currently disputed inconsistencies between energy and water cycle estimates.

Measuring the Diversity and the Integration of Film Production in Mercosur: A Social Network Analysis

Based on an original and comprehensive database of all feature fiction films produced in Mercosur between 2004 and 2012, the paper analyses whether the Mercosur film industry has evolved towards an integrated and culturally more diverse market. It provides a summary of policy opportunities in terms of integration and diversity, emphasizing the limiter role played by regional policies. It then shows that although the Mercosur film industry remains rather disintegrated, it tends to become more integrated and culturally more diverse. From a methodological point of view, the combination of Social Network Analysis and the Stirling Model opens up interesting research tracks to analyse creative industries in terms of their market integration and their cultural diversity.

Changes in biodiversity and trade-offs among ecosystem services, stakeholders, and components of well-being: the contribution of the International Long-Term Ecological Research network (ILTER) to Programme on Ecosystem Change and Society (PECS)

The International Long-Term Ecological Research (ILTER) network comprises > 600 scientific groups conducting site-based research within 40 countries. Its mission includes improving the understanding of global ecosystems and informs solutions to current and future environmental problems at the global scales. The ILTER network covers a wide range of social-ecological conditions and is aligned with the Programme on Ecosystem Change and Society (PECS) goals and approach. Our aim is to examine and develop the conceptual basis for proposed collaboration between ILTER and PECS. We describe how a coordinated effort of several contrasting LTER site-based research groups contributes to the understanding of how policies and technologies drive either toward or away from the sustainable delivery of ecosystem services. This effort is based on three tenets: transdisciplinary research; cross-scale interactions and subsequent dynamics; and an ecological stewardship orientation. The overarching goal is to design management practices taking into account trade-offs between using and conserving ecosystems toward more sustainable solutions. To that end, we propose a conceptual approach linking ecosystem integrity, ecosystem services, and stakeholder well-being, and as a way to analyze trade-offs among ecosystem services inherent in diverse management options. We also outline our methodological approach that includes: (i) monitoring and synthesis activities following spatial and temporal trends and changes on each site and by documenting cross-scale interactions; (ii) developing analytical tools for integration; (iii) promoting trans-site comparison; and (iv) developing conceptual tools to design adequate policies and management interventions to deal with trade-offs. Finally, we highlight the heterogeneity in the social-ecological setting encountered in a subset of 15 ILTER sites. These study cases are diverse enough to provide a broad cross-section of contrasting ecosystems with different policy and management drivers of ecosystem conversion; distinct trends of biodiversity change; different stakeholders’ preferences for ecosystem services; and diverse components of well-being issues.

Apicomplexans pulling the strings: manipulation of the host cell cytoskeleton dynamics

Invasive stages of apicomplexan parasites require a host cell to survive, proliferate and advance to the next life cycle stage. Once invasion is achieved, apicomplexans interact closely with the host cell cytoskeleton, but in many cases the different species have evolved distinct mechanisms and pathways to modulate the structural organization of cytoskeletal filaments. The host cell cytoskeleton is a complex network, largely, but not exclusively, composed of microtubules, actin microfilaments and intermediate filaments, all of which are modulated by associated proteins, and it is involved in diverse functions including maintenance of cell morphology and mechanical support, migration, signal transduction, nutrient uptake, membrane and organelle trafficking and cell division. The ability of apicomplexans to modulate the cytoskeleton to their own advantage is clearly beneficial. We here review different aspects of the interactions of apicomplexans with the three main cytoskeletal filament types, provide information on the currently known parasite effector proteins and respective host cell targets involved, and how these interactions modulate the host cell physiology. Some of these findings could provide novel targets that could be exploited for the development of preventive and/or therapeutic strategies.

Large-scale dynamic hydrofracturing, healing and fracture network characterization

Permeability of a rock is a dynamic property that varies spatially and temporally. Fractures provide the most efficient channels for fluid flow and thus directly contribute to the permeability of the system. Fractures usually form as a result of a combination of tectonic stresses, gravity (i.e. lithostatic pressure) and fluid pressures. High pressure gradients alone can cause fracturing, the process which is termed as hydrofracturing that can determine caprock (seal) stability or reservoir integrity. Fluids also transport mass and heat, and are responsible for the formation of veins by precipitating minerals within open fractures. Veining (healing) thus directly influences the rock’s permeability. Upon deformation these closed factures (veins) can refracture and the cycle starts again. This fracturing-healing-refacturing cycle is a fundamental part in studying the deformation dynamics and permeability evolution of rock systems. This is generally accompanied by fracture network characterization focusing on network topology that determines network connectivity. Fracture characterization allows to acquire quantitative and qualitative data on fractures and forms an important part of reservoir modeling. This thesis highlights the importance of fracture-healing and veins’ mechanical properties on the deformation dynamics. It shows that permeability varies spatially and temporally, and that healed systems (veined rocks) should not be treated as fractured systems (rocks without veins). Field observations also demonstrate the influence of contrasting mechanical properties, in addition to the complexities of vein microstructures that can form in low-porosity and permeability layered sequences. The thesis also presents graph theory as a characterization method to obtain statistical measures on evolving network connectivity. It also proposes what measures a good reservoir should have to exhibit potentially large permeability and robustness against healing. The results presented in the thesis can have applications for hydrocarbon and geothermal reservoir exploration, mining industry, underground waste disposal, CO2 injection or groundwater modeling.

Continuum tensor network field states, path integral representations and spatial symmetries

A natural way to generalize tensor network variational classes to quantum field systems is via a continuous tensor contraction. This approach is first illustrated for the class of quantum field states known as continuous matrix-product states (cMPS). As a simple example of the path-integral representation we show that the state of a dynamically evolving quantum field admits a natural representation as a cMPS. A completeness argument is also provided that shows that all states in Fock space admit a cMPS representation when the number of variational parameters tends to infinity. Beyond this, we obtain a well-behaved field limit of projected entangled-pair states (PEPS) in two dimensions that provide an abstract class of quantum field states with natural symmetries. We demonstrate how symmetries of the physical field state are encoded within the dynamics of an auxiliary field system of one dimension less. In particular, the imposition of Euclidean symmetries on the physical system requires that the auxiliary system involved in the class' definition must be Lorentz-invariant. The physical field states automatically inherit entropy area laws from the PEPS class, and are fully described by the dissipative dynamics of a lower dimensional virtual field system. Our results lie at the intersection many-body physics, quantum field theory and quantum information theory, and facilitate future exchanges of ideas and insights between these disciplines.