Computational modeling of psychiatric illnesses via well-defined neurophysiological and neurocognitive biomarkers

The beta component of gamma-band auditory steady-state responses in patients with schizophrenia

The mechanisms underlying circuit dysfunctions in schizophrenia (SCZ) remain poorly understood. Auditory steady-state responses (ASSRs), especially in the gamma and beta band, have been suggested as a potential biomarker for SCZ. While the reduction of 40 Hz power for 40 Hz drive has been well established and replicated in SCZ patients, studies are inconclusive when it comes to an increase in 20 Hz power during 40 Hz drive. There might be several factors explaining the inconsistencies, including differences in the sensitivity of the recording modality (EEG vs MEG), differences in stimuli (click-trains vs amplitude-modulated tones) and large differences in the amplitude of the stimuli. Here, we used a computational model of ASSR deficits in SCZ and explored the effect of three SCZ-associated microcircuit alterations: reduced GABA activity, increased GABA decay times and NMDA receptor hypofunction. We investigated the effect of input strength on gamma (40 Hz) and beta (20 Hz) band power during gamma ASSR stimulation and saw that the pronounced increase in beta power during gamma stimulation seen experimentally could only be reproduced in the model when GABA decay times were increased and only for a specific range of input strengths. More specifically, when the input was in this specific range, the rhythmic drive at 40 Hz produced a strong 40 Hz rhythm in the control network; however, in the 'SCZ-like' network, the prolonged inhibition led to a so-called 'beat-skipping', where the network would only strongly respond to every other input. This mechanism was responsible for the emergence of the pronounced 20 Hz beta peak in the power spectrum. The other two microcircuit alterations were not able to produce a substantial 20 Hz component but they further narrowed the input strength range for which the network produced a beta component when combined with increased GABAergic decay times. Our finding that the beta component only existed for a specific range of input strengths might explain the seemingly inconsistent reporting in experimental studies and suggests that future ASSR studies should systematically explore different amplitudes of their stimuli. Furthermore, we provide a mechanistic link between a microcircuit alteration and an electrophysiological marker in schizophrenia and argue that more complex ASSR stimuli are needed to disentangle the nonlinear interactions of microcircuit alterations. The computational modelling approach put forward here is ideally suited to facilitate the development of such stimuli in a theory-based fashion.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

A Computational Model of Major Depression: the Role of Glutamate Dysfunction on Cingulo-Frontal Network Dynamics

Major depression disease (MDD) is associated with the dysfunction of multinode brain networks. However, converging evidence implicates the reciprocal interaction between midline limbic regions (typified by the ventral anterior cingulate cortex, vACC) and the dorso-lateral prefrontal cortex (dlPFC), reflecting interactions between emotions and cognition. Furthermore, growing evidence suggests a role for abnormal glutamate metabolism in the vACC, while serotonergic treatments (selective serotonin reuptake inhibitor, SSRI) effective for many patients implicate the serotonin system. Currently, no mechanistic framework describes how network dynamics, glutamate, and serotonin interact to explain MDD symptoms and treatments. Here, we built a biophysical computational model of 2 areas (vACC and dlPFC) that can switch between emotional and cognitive processing. MDD networks were simulated by slowing glutamate decay in vACC and demonstrated sustained vACC activation. This hyperactivity was not suppressed by concurrent dlPFC activation and interfered with expected dlPFC responses to cognitive signals, mimicking cognitive dysfunction seen in MDD. Simulation of clinical treatments (SSRI or deep brain stimulation) counteracted this aberrant vACC activity. Theta and beta/gamma oscillations correlated with network function, representing markers of switch-like operation in the network. The model shows how glutamate dysregulation can cause aberrant brain dynamics, respond to treatments, and be reflected in EEG rhythms as biomarkers of MDD.

An in Silico, Biomarker-Based Method for the Evaluation of Virtual Neuropsychiatric Drug Effects

The recent explosion in neuroscience research has markedly increased our understanding of the neurobiological correlates of many psychiatric illnesses, but this has unfortunately not translated into more effective pharmacologic treatments for these conditions. At the same time, researchers have increasingly sought out biological markers, or biomarkers, as a way to categorize psychiatric illness, as these are felt to be closer to underlying genetic and neurobiological vulnerabilities. While biomarker-based drug discovery approaches have tended to employ in vivo (e.g., rodent) or in vitro test systems, relatively little attention has been paid to the potential of computational, or in silico, methodologies. Here we describe such a methodology, using as an example a biophysically detailed computational model of hippocampus that is made to generate putative schizophrenia biomarkers by the inclusion of a number of neuropathological changes that have been associated with the illness (NMDA system deficit, decreased neural connectivity, hyperdopaminergia). We use the specific inability to attune to gamma band (40 Hz) auditory stimulus as our illness biomarker. We expose this system to a large number of virtual medications, defined by systematic variation of model parameters corresponding to five cellular-level effects. The potential efficacy of virtual medications is determined by a wellness metric (WM) that we have developed. We identify a number of virtual agents that consist of combinations of mechanisms, which are not simply reversals of the causative lesions. The manner in which this methodology could be extended to other neuropsychiatric conditions, such as Alzheimer's disease, autism, and fragile X syndrome, is discussed.

Semi-mechanistic computer simulation of psychotic symptoms in schizophrenia with a model of a humanized cortico-striatal-thalamocortical loop

Despite new insights into the pathophysiology of schizophrenia and clinical trials with highly selective drugs, no new therapeutic breakthroughs have been identified. We present a semi-mechanistic Quantitative Systems Pharmacology (QSP) computer model of a biophysically realistic cortical-striatal-thalamo-cortical loop. The model incorporates the direct, indirect and hyperdirect pathway of the basal ganglia and CNS drug targets that modulate neuronal firing, based on preclinical data about their localization and coupling to voltage-gated ion channels. Schizophrenia pathology is introduced using quantitative human imaging data on striatal hyperdopaminergic activity and cortical dysfunction. We identified an entropy measure of neuronal firing in the thalamus, related to the bandwidth of information processing that correlates well with reported historical clinical changes on PANSS Total with antipsychotics after introduction of their pharmacology (42 drug-dose combinations, r²=0.62). This entropy measure is further validated by predicting the clinical outcome of 28 other novel stand-alone interventions, 14 of them with non-dopamine D₂R pharmacology, in addition to 8 augmentation trials (correlation between actual and predicted clinical scores r²=0.61). The platform predicts that most combinations of antipsychotics have a lower efficacy over what can be achieved by either one; negative pharmacodynamical interactions are prominent for aripiprazole added to risperidone, haloperidol, quetiapine and paliperidone. The model also recapitulates the increased probability for psychotic breakdown in a supersensitive environment and the effect of ketamine in healthy volunteers. This QSP platform, combined with similar readouts for motor symptoms, negative symptoms and cognitive impairment has the potential to improve our understanding of drug effects in schizophrenia patients.

Multifactorial Modeling of Impairment of Evoked Gamma Range Oscillations in Schizophrenia

Despite a significant increase in efforts to identify biomarkers and endophenotypic measures of psychiatric illnesses, only a very limited amount of computational models of these markers and measures has been implemented so far. Moreover, existing computational models dealing with biomarkers typically only examine one possible mechanism in isolation, disregarding the possibility that other combinations of model parameters might produce the same network behavior (what has been termed "multifactoriality"). In this study we describe a step toward a computational instantiation of an endophenotypic finding for schizophrenia, namely the impairment of evoked auditory gamma and beta oscillations in schizophrenia. We explore the multifactorial nature of this impairment using an established model of primary auditory cortex, by performing an extensive search of the parameter space. We find that single network parameters contain only little information about whether the network will show impaired gamma entrainment and that different regions in the parameter space yield similar network level oscillation abnormalities. These regions in the parameter space, however, show strong differences in the underlying network dynamics. To sum up, we present a first step toward an in silico instantiation of an important biomarker of schizophrenia, which has great potential for the identification and study of disease mechanisms and for understanding of existing treatments and development of novel ones.