A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [13, 12], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [13] to further simplify the kinetic theory.

The spatio-temporal brain dynamics of processing and integrating sound localization cues in humans.

Interaural intensity and time differences (IID and ITD) are two binaural auditory cues for localizing sounds in space. This study investigated the spatio-temporal brain mechanisms for processing and integrating IID and ITD cues in humans. Auditory-evoked potentials were recorded, while subjects passively listened to noise bursts lateralized with IID, ITD or both cues simultaneously, as well as a more frequent centrally presented noise. In a separate psychophysical experiment, subjects actively discriminated lateralized from centrally presented stimuli. IID and ITD cues elicited different electric field topographies starting at approximately 75 ms post-stimulus onset, indicative of the engagement of distinct cortical networks. By contrast, no performance differences were observed between IID and ITD cues during the psychophysical experiment. Subjects did, however, respond significantly faster and more accurately when both cues were presented simultaneously. This performance facilitation exceeded predictions from probability summation, suggestive of interactions in neural processing of IID and ITD cues. Supra-additive neural response interactions as well as topographic modulations were indeed observed approximately 200 ms post-stimulus for the comparison of responses to the simultaneous presentation of both cues with the mean of those to separate IID and ITD cues. Source estimations revealed differential processing of IID and ITD cues initially within superior temporal cortices and also at later stages within temporo-parietal and inferior frontal cortices. Differences were principally in terms of hemispheric lateralization. The collective psychophysical and electrophysiological results support the hypothesis that IID and ITD cues are processed by distinct, but interacting, cortical networks that can in turn facilitate auditory localization.

Mechanical external work and recovery at preferred walking speed in obese subjects

PURPOSE: The aim of this study was to compare the mechanical external work (per kg) and pendular energy transduction at preferred walking speed (PWS) in obese versus normal body mass subjects to investigate whether obese adults adopt energy conserving gait mechanics. METHODS: The mechanical external work (Wext) and the fraction of mechanical energy recovered by the pendular mechanism (Rstep) were computed using kinematic data acquired by an optoelectronic system and were compared in 30 obese (OG; body mass index [BMI] = 39.6 +/- 0.6 kg m(-2); 29.5 +/- 1.3 yr) and 19 normal body mass adults (NG; BMI = 21.4 +/- 0.5 kg m(-2); 31.2 +/- 1.2 yr) walking at PWS. RESULTS: PWS was significantly lower in OG (1.18 +/- 0.02 m s(-1)) than in NG (1.33 +/- 0.02 m s(-1); P <or= 0.001). There was no significant difference in Wext per unit mass between groups (OG: 0.36 +/- 0.03 J kg(-1) m(-1); NG: 0.31 +/- 0.02 J kg(-1) m(-1); P = 0.12). Rstep was significantly lower in OG (68.4% +/- 2.0%) compared with NG (74.4% +/- 1.0%; P = 0.01). In OG only, Wext per unit mass was positively correlated with PWS (r = 0.57; P < 0.001). CONCLUSION: Obese adults do not appear to alter their gait to improve pendular energy transduction and may select slower PWS to reduce mechanical and metabolic work.

Dynamics of phonological-phonetic encoding in word production: Evidence from diverging ERPs between stroke patients and controls.

While the dynamics of lexical-semantic and lexical-phonological encoding in word production have been investigated in several event-related potential (ERP) studies, the estimated time course of phonological-phonetic encoding is the result of rather indirect evidence. We investigated the dynamics of phonological-phonetic encoding combining ERP analyses covering the entire encoding process in picture naming and word reading tasks by comparing ERP modulations in eight brain-damaged speakers presenting impaired phonological-phonetic encoding relative to 16 healthy controls. ERPs diverged between groups in terms of local waveform amplitude and global topography at ∼400ms after stimulus onset in the picture naming task and at ∼320-350ms in word reading and sustained until 100ms before articulation onset. These divergences appeared in later time windows than those found in patients with underlying lexical-semantic and lexical-phonological impairment in previous studies, providing evidence that phonological-phonetic encoding is engaged around 400ms in picture naming and around 330ms in word reading.

Banning Obsolete Weapons or Reshaping Perceptions of Military Utility: Discursive Dynamics in Weapons Prohibitions

The paper focuses on the argumentative process through which new international norms prohibiting the use of weapons causing severe civilian harm emerge. It examines the debate surrounding the use and usefulness of landmines and cluster munitions and traces the process through which NGOs change conceptions of military utility and effectiveness of certain weapons by highlighting their humanitarian problems and questioning their military value. By challenging military thinking on these issues, NGOs redefine the terms of the debate – from a commonplace practice, the use of such weapons becomes controversial and military decisions need to be justified. The argument-counterargument dynamic shifts the burden of proof of the necessity and safety of the weapons to the users. The process witnesses the ability of NGOs to influence debates on military issues despite their disadvantaged position in hard security issue areas. It also challenges realist assumptions that only weapons that are obsolete or low-cost force equalizers for weak actors can be banned. To the contrary, the paper shows that in the case of landmines and cluster munitions, defining the military (in)effectiveness of the weapons is part and parcel of the struggle for their prohibition.

Understanding the low-temperature limitations to forest growth through calibration of a forest dynamics model with tree-ring data

The sensitivity of altitudinal and latitudinal tree-line ecotones to climate change, particularly that of temperature, has received much attention. To improve our understanding of the factors affecting tree-line position, we used the spatially explicit dynamic forest model TreeMig. Although well-suited because of its landscape dynamics functions, TreeMig features a parabolic temperature growth response curve, which has recently been questioned. and the species parameters are not specifically calibrated for cold temperatures. Our main goals were to improve the theoretical basis of the temperature growth response curve in the model and develop a method for deriving that curve's parameters from tree-ring data. We replaced the parabola with an asymptotic curve, calibrated for the main species at the subalpine (Swiss Alps: Pinus cembra, Larix decidua, Picea abies) and boreal (Fennoscandia: Pinus sylvestris, Betula pubescens, P. abies) tree-lines. After fitting new parameters, the growth curve matched observed tree-ring widths better. For the subalpine species, the minimum degree-day sum allowing, growth (kDDMin) was lowered by around 100 degree-days; in the case of Larix, the maximum potential ring-width was increased to 5.19 mm. At the boreal tree-line, the kDDMin for P. sylvestris was lowered by 210 degree-days and its maximum ring-width increased to 2.943 mm; for Betula (new in the model) kDDMin was set to 325 degree-days and the maximum ring-width to 2.51 mm; the values from the only boreal sample site for Picea were similar to the subalpine ones, so the same parameters were used. However, adjusting the growth response alone did not improve the model's output concerning species' distributions and their relative importance at tree-line. Minimum winter temperature (MinWiT, mean of the coldest winter month), which controls seedling establishment in TreeMig, proved more important for determining distribution. Picea, P. sylvestris and Betula did not previously have minimum winter temperature limits, so these values were set to the 95th percentile of each species' coldest MinWiT site (respectively -7, -11, -13). In a case study for the Alps, the original and newly calibrated versions of TreeMig were compared with biomass data from the National Forest Inventor), (NFI). Both models gave similar, reasonably realistic results. In conclusion, this method of deriving temperature responses from tree-rings works well. However, regeneration and its underlying factors seem more important for controlling species' distributions than previously thought. More research on regeneration ecology, especially at the upper limit of forests. is needed to improve predictions of tree-line responses to climate change further.

Dynamics of experimental populations of native and introduced blowflies (Diptera: Calliphoridae): mathematical modelling and the transition from asymptotic equilibrium to bounded oscillations

The equilibrium dynamics of native and introduced blowflies is modelled using a density-dependent model of population growth that takes into account important features of the life-history in these flies. A theoretical analysis indicates that the product of maximum fecundity and survival is the primary determinant of the dynamics. Cochliomyia macellaria, a blowfly native to the Americas and the introduced Chrysomya megacephala and Chrysomya putoria, differ in their dynamics in that the first species shows a damping oscillatory behavior leading to a one-point equilibrium, whereas in the last two species population numbers show a two-point limit cycle. Simulations showed that variation in fecundity has a marked effect on the dynamics and indicates the possibility of transitions from one-point equilibrium to bounded oscillations and aperiodic behavior. Variation in survival has much less influence on the dynamics.

The new direct oral anticoagulants in special indications: rationale and preliminary data in cancer, mechanical heart valves, anti-phospholipid syndrome, and heparin-induced thrombocytopenia and beyond.

The present review will briefly summarize the interplay between coagulation and inflammation, highlighting possible effects of direct inhibition of factor Xa and thrombin beyond anticoagulation. Additionally, the rationale for the use of the new direct oral anticoagulants (DOACs) for indications such as cancer-associated venous thromboembolism (CAT), mechanical heart valves, thrombotic anti-phospholipid syndrome (APS), and heparin-induced thrombocytopenia (HIT) will be explored. Published data on patients with cancer or mechanical heart valves treated with DOAC will be discussed, as well as planned studies in APS and HIT. Although at the present time published evidence is insufficient for recommending DOAC in the above-mentioned indications, there are good arguments in favor of clinical trials investigating their efficacy in these contexts. Direct inhibition of factor Xa or thrombin may reveal interesting effects beyond anticoagulation as well.

Direct Determination of mechanical properties in shallow formations by (ultra) sonic methods

For several decades mechanical properties of shallow formations (soil) obtained by sonic to ultrasonic wave testing were reported to be greater than those based on mechanical tests. The present article relying on a statistical analysis of more than 300 tests shows that elastic moduli of the soil can indeed be obtained from (ultra)sonic tests and that they are identical to those resulting from mechanical tests.

Structural prediction of peptides bound to MHC class I.

An ab initio structure prediction approach adapted to the peptide-major histocompatibility complex (MHC) class I system is presented. Based on structure comparisons of a large set of peptide-MHC class I complexes, a molecular dynamics protocol is proposed using simulated annealing (SA) cycles to sample the conformational space of the peptide in its fixed MHC environment. A set of 14 peptide-human leukocyte antigen (HLA) A0201 and 27 peptide-non-HLA A0201 complexes for which X-ray structures are available is used to test the accuracy of the prediction method. For each complex, 1000 peptide conformers are obtained from the SA sampling. A graph theory clustering algorithm based on heavy atom root-mean-square deviation (RMSD) values is applied to the sampled conformers. The clusters are ranked using cluster size, mean effective or conformational free energies, with solvation free energies computed using Generalized Born MV 2 (GB-MV2) and Poisson-Boltzmann (PB) continuum models. The final conformation is chosen as the center of the best-ranked cluster. With conformational free energies, the overall prediction success is 83% using a 1.00 Angstroms crystal RMSD criterion for main-chain atoms, and 76% using a 1.50 Angstroms RMSD criterion for heavy atoms. The prediction success is even higher for the set of 14 peptide-HLA A0201 complexes: 100% of the peptides have main-chain RMSD values &lt; or =1.00 Angstroms and 93% of the peptides have heavy atom RMSD values &lt; or =1.50 Angstroms. This structure prediction method can be applied to complexes of natural or modified antigenic peptides in their MHC environment with the aim to perform rational structure-based optimizations of tumor vaccines.

Regulation of microtubule dynamics by the neuronal growth-associated protein SCG10.

Dynamic assembly and disassembly of microtubules is essential for cell division, cell movements, and intracellular transport. In the developing nervous system, microtubule dynamics play a fundamental role during neurite outgrowth, elongation, and branching, but the molecular mechanisms involved are unknown. SCG10 is a neuron-specific protein that is membrane-associated and highly enriched in growth cones. Here we show that SCG10 binds to microtubules, inhibits their assembly, and can induce microtubule disassembly. We also show that SCG10 overexpression enhances neurite outgrowth in a stably transfected neuronal cell line. These data identify SCG10 as a key regulator of neurite extension through regulation of microtubule instability.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex.

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.