Saltatory Conduction along Myelinated Axons Involves a Periaxonal Nanocircuit.

The propagation of electrical impulses along axons is highly accelerated by the myelin sheath and produces saltating or "jumping" action potentials across internodes, from one node of Ranvier to the next. The underlying electrical circuit, as well as the existence and role of submyelin conduction in saltatory conduction remain, however, elusive. Here, we made patch-clamp and high-speed voltage-calibrated optical recordings of potentials across the nodal and internodal axolemma of myelinated neocortical pyramidal axons combined with electron microscopy and experimentally constrained cable modeling. Our results reveal a nanoscale yet conductive periaxonal space, incompletely sealed at the paranodes, which separates the potentials across the low-capacitance myelin sheath and internodal axolemma. The emerging double-cable model reproduces the recorded evolution of voltage waveforms across nodes and internodes, including rapid nodal potentials traveling in advance of attenuated waves in the internodal axolemma, revealing a mechanism for saltation across time and space.

A Growth-Based Framework for Leaf Shape Development and Diversity.

How do genes modify cellular growth to create morphological diversity? We study this problem in two related plants with differently shaped leaves: Arabidopsis thaliana (simple leaf shape) and Cardamine hirsuta (complex shape with leaflets). We use live imaging, modeling, and genetics to deconstruct these organ-level differences into their cell-level constituents: growth amount, direction, and differentiation. We show that leaf shape depends on the interplay of two growth modes: a conserved organ-wide growth mode that reflects differentiation; and a local, directional mode that involves the patterning of growth foci along the leaf edge. Shape diversity results from the distinct effects of two homeobox genes on these growth modes: SHOOTMERISTEMLESS broadens organ-wide growth relative to edge-patterning, enabling leaflet emergence, while REDUCED COMPLEXITY inhibits growth locally around emerging leaflets, accentuating shape differences created by patterning. We demonstrate the predictivity of our findings by reconstructing key features of C. hirsuta leaf morphology in A. thaliana. VIDEO ABSTRACT.

Local and global changes in cell density induce reorganisation of 3D packing in a proliferating epithelium.

Tissue morphogenesis is intimately linked to the changes in shape and organisation of individual cells. In curved epithelia, cells can intercalate along their own apicobasal axes, adopting a shape named 'scutoid' that allows energy minimization in the tissue. Although several geometric and biophysical factors have been associated with this 3D reorganisation, the dynamic changes underlying scutoid formation in 3D epithelial packing remain poorly understood. Here, we use live imaging of the sea star embryo coupled with deep learning-based segmentation to dissect the relative contributions of cell density, tissue compaction and cell proliferation on epithelial architecture. We find that tissue compaction, which naturally occurs in the embryo, is necessary for the appearance of scutoids. Physical compression experiments identify cell density as the factor promoting scutoid formation at a global level. Finally, the comparison of the developing embryo with computational models indicates that the increase in the proportion of scutoids is directly associated with cell divisions. Our results suggest that apico-basal intercalations appearing immediately after mitosis may help accommodate the new cells within the tissue. We propose that proliferation in a compact epithelium induces 3D cell rearrangements during development.

Impaired value-based decision-making in Parkinson's disease apathy.

Apathy is a common and disabling complication of Parkinson's disease characterized by reduced goal-directed behaviour. Several studies have reported dysfunction within prefrontal cortical regions and projections from brainstem nuclei whose neuromodulators include dopamine, serotonin and noradrenaline. Work in animal and human neuroscience have confirmed contributions of these neuromodulators on aspects of motivated decision-making. Specifically, these neuromodulators have overlapping contributions to encoding the value of decisions, and influence whether to explore alternative courses of action or persist in an existing strategy to achieve a rewarding goal. Building upon this work, we hypothesized that apathy in Parkinson's disease should be associated with an impairment in value-based learning. Using a four-armed restless bandit reinforcement learning task, we studied decision-making in 75 volunteers; 53 patients with Parkinson's disease, with and without clinical apathy, and 22 age-matched healthy control subjects. Patients with apathy exhibited impaired ability to choose the highest value bandit. Task performance predicted an individual patient's apathy severity measured using the Lille Apathy Rating Scale (R = -0.46, P < 0.001). Computational modelling of the patient's choices confirmed the apathy group made decisions that were indifferent to the learnt value of the options, consistent with previous reports of reward insensitivity. Further analysis demonstrated a shift away from exploiting the highest value option and a reduction in perseveration, which also correlated with apathy scores (R = -0.5, P < 0.001). We went on to acquire functional MRI in 59 volunteers; a group of 19 patients with and 20 without apathy and 20 age-matched controls performing the Restless Bandit Task. Analysis of the functional MRI signal at the point of reward feedback confirmed diminished signal within ventromedial prefrontal cortex in Parkinson's disease, which was more marked in apathy, but not predictive of their individual apathy severity. Using a model-based categorization of choice type, decisions to explore lower value bandits in the apathy group activated prefrontal cortex to a similar degree to the age-matched controls. In contrast, Parkinson's patients without apathy demonstrated significantly increased activation across a distributed thalamo-cortical network. Enhanced activity in the thalamus predicted individual apathy severity across both patient groups and exhibited functional connectivity with dorsal anterior cingulate cortex and anterior insula. Given that task performance in patients without apathy was no different to the age-matched control subjects, we interpret the recruitment of this network as a possible compensatory mechanism, which compensates against symptomatic manifestation of apathy in Parkinson's disease.

Interactions between calcium-induced arrhythmia triggers and the electrophysiological-anatomical substrate underlying the induction of atrial fibrillation.

Atrial fibrillation (AF) is the most common cardiac arrhythmia and is sustained by spontaneous focal excitations and re-entry. Spontaneous electrical firing in the pulmonary vein (PV) sleeves is implicated in AF generation. The aim of this simulation study was to identify the mechanisms determining the localisation of AF triggers in the PVs and their contribution to the genesis of AF. A novel biophysical model of the canine atria was used that integrates stochastic, spontaneous subcellular Ca2+ release events (SCRE) with regional electrophysiological heterogeneity in ionic properties and a detailed three-dimensional model of atrial anatomy, microarchitecture and patchy fibrosis. Simulations highlighted the importance of the smaller inward rectifier potassium current (IK1 ) in PV cells compared to the surrounding atria, which enabled SCRE more readily to result in delayed-afterdepolarisations that induced triggered activity. There was a leftward shift in the dependence of the probability of triggered activity on sarcoplasmic reticulum Ca2+ load. This feature was accentuated in 3D tissue compared to single cells (Δ half-maximal [Ca2+ ]SR  = 58 µM vs. 22 µM). In 3D atria incorporating electrical heterogeneity, excitations preferentially emerged from the PV region. These triggered focal excitations resulted in transient re-entry in the left atrium. Addition of fibrotic patches promoted localised emergence of focal excitations and wavebreaks that had a more substantial impact on generating AF-like patterns than the PVs. Thus, a reduced IK1 , less negative resting membrane potential, and fibrosis-induced changes of the electrotonic load all contribute to the emergence of complex excitation patterns from spontaneous focal triggers. KEY POINTS: Focal excitations in the atria are most commonly associated with the pulmonary veins, but the mechanisms for this localisation are yet to be elucidated. We applied a multi-scale computational modelling approach to elucidate the mechanisms underlying such localisations. Myocytes in the pulmonary vein region of the atria have a less negative resting membrane potential and reduced time-independent potassium current; we demonstrate that both of these factors promote triggered activity in single cells and tissues. The less negative resting membrane potential also contributes to heterogeneous inactivation of the fast sodium current, which can enable re-entrant-like excitation patterns to emerge without traditional conduction block.

Pre-supplementary motor area strengthens reward sensitivity in intertemporal choice.

Previous investigations on the causal neural mechanisms underlying intertemporal decision making focused on the dorsolateral prefrontal cortex as neural substrate of cognitive control. However, little is known, about the causal contributions of further parts of the frontoparietal control network to delaying gratification, including the pre-supplementary motor area (pre-SMA) and posterior parietal cortex (PPC). Conflicting previous evidence related pre-SMA and PPC either to evidence accumulation processes, choice biases, or response caution. To disentangle between these alternatives, we combined drift diffusion models of decision making with online transcranial magnetic stimulation (TMS) over pre-SMA and PPC during an intertemporal decision task. While we observed no robust effects of PPC TMS, perturbation of pre-SMA activity reduced preferences for larger over smaller rewards. A drift diffusion model of decision making suggests that pre-SMA increases the weight assigned to reward magnitudes during the evidence accumulation process without affecting choice biases or response caution. Taken together, the current findings reveal the computational role of the pre-SMA in value-based decision making, showing that pre-SMA promotes choices of larger, costly rewards by strengthening the sensitivity to reward magnitudes.

Evoked resonant neural activity long-term dynamics can be reproduced by a computational model with vesicle depletion.

Subthalamic deep brain stimulation (DBS) robustly generates high-frequency oscillations known as evoked resonant neural activity (ERNA). Recently the importance of ERNA has been demonstrated through its ability to predict the optimal DBS contact in the subthalamic nucleus in patients with Parkinson's disease. However, the underlying mechanisms of ERNA are not well understood, and previous modelling efforts have not managed to reproduce the wealth of published data describing the dynamics of ERNA. Here, we aim to present a minimal model capable of reproducing the characteristics of the slow ERNA dynamics published to date. We make biophysically-motivated modifications to the Kuramoto model and fit its parameters to the slow dynamics of ERNA obtained from data. Our results demonstrate that it is possible to reproduce the slow dynamics of ERNA (over hundreds of seconds) with a single neuronal population, and, crucially, with vesicle depletion as one of the key mechanisms behind the ERNA frequency decay in our model. We further validate the proposed model against experimental data from Parkinson's disease patients, where it captures the variations in ERNA frequency and amplitude in response to variable stimulation frequency, amplitude, and to stimulation pulse bursting. We provide a series of predictions from the model that could be the subject of future studies for further validation.

Histamine in the neocortex: Towards integrating multiscale effectors.

Histamine is a modulatory neurotransmitter, which has received relatively less attention in the central nervous system than other neurotransmitters. The functional role of histamine in the neocortex, the brain region that controls higher-order cognitive functions such as attention, learning and memory, remains largely unknown. This article focuses on the emerging roles and mechanisms of histamine release in the neocortex. We describe gaps in current knowledge and propose the application of interdisciplinary tools to dissect the detailed multiscale functional logic of histaminergic action in the neocortex ranging from sub-cellular, cellular, dendritic and synaptic levels to microcircuits and mesoscale effects.

Bayesian reduced rank regression models generalizable neural fingerprints that differentiate between individuals in magnetoencephalography data.

Recent magnetoencephalography (MEG) studies have reported that functional connectivity (FC) and power spectra can be used as neural fingerprints in differentiating individuals. Such studies have mainly used correlations between measurement sessions to distinguish individuals from each other. However, it has remained unclear whether such correlations might reflect a more generalizable principle of individually distinctive brain patterns. Here, we evaluated a machine-learning based approach, termed latent-noise Bayesian reduced rank regression (BRRR) as a means of modelling individual differences in the resting-state MEG data of the Human Connectome Project (HCP), using FC and power spectra as neural features. First, we verified that BRRR could model and reproduce the differences between metrics that correlation-based fingerprinting yields. We trained BRRR models to distinguish individuals based on data from one measurement and used the models to identify subsequent measurement sessions of those same individuals. The best performing BRRR models, using only 20 spatiospectral components, were able to identify subjects across measurement sessions with over 90% accuracy, approaching the highest correlation-based accuracies. Using cross-validation, we then determined whether that BRRR model could generalize to unseen subjects, successfully classifying the measurement sessions of novel individuals with over 80% accuracy. The results demonstrate that individual neurofunctional differences can be reliably extracted from MEG data with a low-dimensional predictive model and that the model is able to classify novel subjects.

Molecular Basis for Chemoselectivity Control in Oxidations of Internal Aryl-Alkenes Catalyzed by Laboratory Evolved P450s.

P450 enzymes naturally perform selective hydroxylations and epoxidations of unfunctionalized hydrocarbon substrates, among other reactions. The adaptation of P450 enzymes to a particular oxidative reaction involving alkenes is of great interest for the design of new synthetically useful biocatalysts. However, the mechanism that these enzymes utilize to precisely modulate the chemoselectivity and distinguishing between competing alkene double bond epoxidations and allylic C-H hydroxylations is sometimes not clear, which hampers the rational design of specific biocatalysts. In a previous work, a P450 from Labrenzia aggregata (P450LA1) was engineered in the laboratory using directed evolution to catalyze the direct oxidation of trans-ß-methylstyrene to phenylacetone. The final variant, KS, was able to overcome the intrinsic preference for alkene epoxidation to directly generate a ketone product via the formation of a highly reactive carbocation intermediate. Here, additional library screening along this evolutionary lineage permitted to serendipitously detect a mutation that overcomes epoxidation and carbonyl formation by exhibiting a large selectivity of 94 % towards allylic C-H hydroxylation. A multiscalar computational methodology was applied to reveal the molecular basis towards this hydroxylation preference. Enzyme modelling suggests that introduction of a bulky substitution dramatically changes the accessible conformations of the substrate in the active site, thus modifying the enzymatic selectivity towards terminal hydroxylation and avoiding the competing epoxidation pathway, which is sterically hindered.

Addressing barriers in comprehensiveness, accessibility, reusability, interoperability and reproducibility of computational models in systems biology.

Computational models are often employed in systems biology to study the dynamic behaviours of complex systems. With the rise in the number of computational models, finding ways to improve the reusability of these models and their ability to reproduce virtual experiments becomes critical. Correct and effective model annotation in community-supported and standardised formats is necessary for this improvement. Here, we present recent efforts toward a common framework for annotated, accessible, reproducible and interoperable computational models in biology, and discuss key challenges of the field.

Computational modeling of light processing in the habenula and dorsal raphe based on laser ablation of functionally-defined cells.

BACKGROUND: The habenula is a major regulator of serotonergic neurons in the dorsal raphe, and thus of brain state. The functional connectivity between these regions is incompletely characterized. Here, we use the ability of changes in irradiance to trigger reproducible changes in activity in the habenula and dorsal raphe of zebrafish larvae, combined with two-photon laser ablation of specific neurons, to establish causal relationships. RESULTS: Neurons in the habenula can show an excitatory response to the onset or offset of light, while neurons in the anterior dorsal raphe display an inhibitory response to light, as assessed by calcium imaging. The raphe response changed in a complex way following ablations in the dorsal habenula (dHb) and ventral habenula (vHb). After ablation of the ON cells in the vHb (V-ON), the raphe displayed no response to light. After ablation of the OFF cells in the vHb (V-OFF), the raphe displayed an excitatory response to darkness. After ablation of the ON cells in the dHb (D-ON), the raphe displayed an excitatory response to light. We sought to develop in silico models that could recapitulate the response of raphe neurons as a function of the ON and OFF cells of the habenula. Early attempts at mechanistic modeling using ordinary differential equation (ODE) failed to capture observed raphe responses accurately. However, a simple two-layer fully connected neural network (NN) model was successful at recapitulating the diversity of observed phenotypes with root-mean-squared error values ranging from 0.012 to 0.043. The NN model also estimated the raphe response to ablation of D-off cells, which can be verified via future experiments. CONCLUSION: Lesioning specific cells in different regions of habenula led to qualitatively different responses to light in the dorsal raphe. A simple neural network is capable of mimicking experimental observations. This work illustrates the ability of computational modeling to integrate complex observations into a simple compact formalism for generating testable hypotheses, and for guiding the design of biological experiments.

Defining the brain control of physiological stability.

The last few decades have seen major advances in neurobiology and uncovered novel genetic and cellular substrates involved in the control of physiological set points. In this Review, I discuss the limitations in the definition of homeostatic set points established by Walter B Canon and highlight evidence that two other physiological systems, namely rheostasis and allostasis provide distinct inputs to independently modify set-point levels. Using data collected over the past decade, the hypothalamic and genetic basis of regulated changes in set-point values by rheostatic mechanisms are described. Then, the role of higher-order brain regions, such as hippocampal circuits, for experience-dependent, allostatic induced changes in set-points are outlined. I propose that these systems provide a hierarchical organization of physiological stability that exists to maintain set-point values. The hierarchical organization of physiology has direct implications for basic and medical research, and clinical practice.

Effort avoidance as a core mechanism of apathy in frontotemporal dementia.

Apathy is a core symptom in patients with behavioural variant frontotemporal dementia (bvFTD). It is defined by the observable reduction in goal-directed behaviour, but the underlying mechanisms are poorly understood. According to decision theory, engagement in goal-directed behaviour depends on a cost-benefit optimization trading off the estimated effort (related to the behaviour) against the expected reward (related to the goal). In this framework, apathy would thus result from either a decreased appetence for reward, or from an increased aversion to effort. Here, we phenotyped the motivational state of 21 patients with bvFTD and 40 matched healthy controls using computational analyses of behavioural responses in a comprehensive series of behavioural tasks, involving both expression of preference (comparing reward value and effort cost) and optimization of performance (adjusting effort production to the reward at stake). The primary finding was an elevated aversion to effort, consistent across preference and performance tasks in patients with bvFTD compared to controls. Within the bvFTD group, effort avoidance was correlated to cortical atrophy in the dorsal anterior cingulate cortex and to apathy score measured on a clinical scale. Thus, our results highlight elevated effort aversion (not reduced reward appetence) as a core dysfunction that might generate apathy in patients with bvFTD. More broadly, they provide novel behavioural tests and computational tools to identify the dysfunctional mechanisms producing motivation deficits in patients with brain damage.

Current developments in elastic and acoustic metamaterials science.

The concept of metamaterial recently emerged as a new frontier of scientific research, encompassing physics, materials science and engineering. In a broad sense, a metamaterial indicates an engineered material with exotic properties not found in nature, obtained by appropriate architecture either at macro-scale or at micro-/nano-scales. The architecture of metamaterials can be tailored to open unforeseen opportunities for mechanical and acoustic applications, as demonstrated by an impressive and increasing number of studies. Building on this knowledge, this theme issue aims to gather cutting-edge theoretical, computational and experimental studies on elastic and acoustic metamaterials, with the purpose of offering a wide perspective on recent achievements and future challenges. This article is part of the theme issue 'Current developments in elastic and acoustic metamaterials science (Part 1)'.

Current developments in elastic and acoustic metamaterials science.

The concept of metamaterial recently emerged as a new frontier of scientific research, encompassing physics, materials science and engineering. In a broad sense, a metamaterial indicates an engineered material with exotic properties not found in nature, obtained by appropriate architecture either at macro-scale or at micro-/nano-scales. The architecture of metamaterials can be tailored to open unforeseen opportunities for mechanical and acoustic applications, as demonstrated by an impressive and increasing number of studies. Building on this knowledge, this theme issue aims to gather cutting-edge theoretical, computational and experimental studies on elastic and acoustic metamaterials, with the purpose of offering a wide perspective on recent achievements and future challenges.This article is part of the theme issue, 'Current developments in elastic and acoustic metamaterials science (Part 2)'.

Antimicrobial activity of compounds identified by artificial intelligence discovery engine targeting enzymes involved in Neisseria gonorrhoeae peptidoglycan metabolism.

BACKGROUND: Neisseria gonorrhoeae (Ng) causes the sexually transmitted disease gonorrhoea. There are no vaccines and infections are treated principally with antibiotics. However, gonococci rapidly develop resistance to every antibiotic class used and there is a need for developing new antimicrobial treatments. In this study we focused on two gonococcal enzymes as potential antimicrobial targets, namely the serine protease L,D-carboxypeptidase LdcA (NgO1274/NEIS1546) and the lytic transglycosylase LtgD (NgO0626/NEIS1212). To identify compounds that could interact with these enzymes as potential antimicrobials, we used the AtomNet virtual high-throughput screening technology. We then did a computational modelling study to examine the interactions of the most bioactive compounds with their target enzymes. The identified compounds were tested against gonococci to determine minimum inhibitory and bactericidal concentrations (MIC/MBC), specificity, and compound toxicity in vitro. RESULTS: AtomNet identified 74 compounds that could potentially interact with Ng-LdcA and 84 compounds that could potentially interact with Ng-LtgD. Through MIC and MBC assays, we selected the three best performing compounds for both enzymes. Compound 16 was the most active against Ng-LdcA, with a MIC50 value < 1.56 µM and MBC50/90 values between 0.195 and 0.39 µM. In general, the Ng-LdcA compounds showed higher activity than the compounds directed against Ng-LtgD, of which compound 45 had MIC50 values of 1.56-3.125 µM and MBC50/90 values between 3.125 and 6.25 µM. The compounds were specific for gonococci and did not kill other bacteria. They were also non-toxic for human conjunctival epithelial cells as judged by a resazurin assay. To support our biological data, in-depth computational modelling study detailed the interactions of the compounds with their target enzymes. Protein models were generated in silico and validated, the active binding sites and amino acids involved elucidated, and the interactions of the compounds interacting with the enzymes visualised through molecular docking and Molecular Dynamics Simulations for 50 ns and Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA). CONCLUSIONS: We have identified bioactive compounds that appear to target the N. gonorrhoeae LdcA and LtgD enzymes. By using a reductionist approach involving biological and computational data, we propose that compound Ng-LdcA-16 and Ng-LtgD-45 are promising anti-gonococcal compounds for further development.

Systems Biology of Ageing.

The ageing process is highly complex involving multiple processes operating at different biological levels. Systems Biology presents an approach using integrative computational and laboratory study that allows us to address such complexity. The approach relies on the computational analysis of knowledge and data to generate predictive models that may be validated with further laboratory experimentation. Our understanding of ageing is such that translational opportunities are within reach and systems biology offers a means to ensure that optimal decisions are made. We present an overview of the methods employed from bioinformatics and computational modelling and describe some of the insights into ageing that have been gained.

A neuromechanical model characterizing the motor planning and posture control in the voluntary lean in Parkinson's disease.

Parkinson's disease targets patients' cognitive and motor abilities, including postural control. Many studies have been carried out to introduce mathematical models for a better understanding of postural control in such patients and the relation between the model parameters and the clinical assessments. So far, these studies have addressed this connection merely in static tests, such as quiet stance. The aim of this study is to develop a model for voluntary lean, and as such, identify the model parameters for both PD patients and healthy subjects from experimental data. The proposed model comprises planning and control sections. The model parameters for the planning section were extracted from the time response characteristics. Parameters for the control section were identified based on the spatial characteristics of the center-of-pressure (COP) response using an optimization process. 24 PD patients along with 24 matched healthy subjects participated in the study. The results showed a significant difference between the two groups in terms of temporal parameters for the planning section. This difference emphasizes bradykinesia as an essential symptom of PD. Also, differences were found for the postural control section. In all directions, the proportional gain of the feedback controller was significantly larger in PD patients; however, the gain of the feedforward controller was significantly smaller in PD patients. Furthermore, the control gains were strongly correlated with the clinical scales (Functional Reach Test and Unified Parkinson's Disease Rating Scale) in certain directions. In conclusion, the new model helps to better understand and quantify some PD symptoms in voluntary lean tasks.

The influence of emotion on temporal context models.

Temporal context models (TCMs) have been influential in understanding episodic memory and its neural underpinnings. Recently, TCMs have been extended to explain emotional memory effects, one of the most clinically important findings in the field of memory research. This review covers recent advances in hypotheses for the neural representation of spatiotemporal context through the lens of TCMs, including their ability to explain the influence of emotion on episodic and temporal memory. In recent years, simplifying assumptions of "classical" TCMs - with exponential trace decay and the mechanism by which temporal context is recovered - have become increasingly clear. The review also outlines how recent advances could be incorporated into a future TCM, beyond classical assumptions, to integrate emotional modulation.