A quantitative system pharmacology computer model for cognitive deficits in schizophrenia

Mechanistic modeling as an explanatory tool for clinical treatment of chronic catatonia

Mathematical modeling of neural systems is an effective means to integrate complex information about the brain into a numerical tool that can help explain observations. However, the use of neural models to inform clinical decisions has been limited. In this study, we use a simple model of brain circuitry, the Wilson-Cowan model, to predict changes in a clinical measure for catatonia, the Bush-Francis Catatonia Rating Scale, for use in clinical treatment of schizophrenia. This computational tool can then be used to better understand mechanisms of action for pharmaceutical treatments, and to fine-tune dosage in individual cases. We present the conditions of clinical care for a residential patient cohort, and describe methods for synthesizing data to demonstrated the functioning of the model. We then show that the model can be used to explain effect sizes of treatments and estimate outcomes for combinations of medications. We conclude with a demonstration of how this model could be personalized for individual patients to inform ongoing treatment protocols.

Hallmarks of neurodegenerative disease: A systems pharmacology perspective

Age-related central neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, are a rising public health concern and have been plagued by repeated drug development failures. The complex nature and poor mechanistic understanding of the etiology of neurodegenerative diseases has hindered the discovery and development of effective disease-modifying therapeutics. Quantitative systems pharmacology models of neurodegeneration diseases may be useful tools to enhance the understanding of pharmacological intervention strategies and to reduce drug attrition rates. Due to the similarities in pathophysiological mechanisms across neurodegenerative diseases, especially at the cellular and molecular levels, we envision the possibility of structural components that are conserved across models of neurodegenerative diseases. Conserved structural submodels can be viewed as building blocks that are pieced together alongside unique disease components to construct quantitative systems pharmacology (QSP) models of neurodegenerative diseases. Model parameterization would likely be different between the different types of neurodegenerative diseases as well as individual patients. Formulating our mechanistic understanding of neurodegenerative pathophysiology as a mathematical model could aid in the identification and prioritization of drug targets and combinatorial treatment strategies, evaluate the role of patient characteristics on disease progression and therapeutic response, and serve as a central repository of knowledge. Here, we provide a background on neurodegenerative diseases, highlight hallmarks of neurodegeneration, and summarize previous QSP models of neurodegenerative diseases.

Toward Predicting Impact of Common Genetic Variants on Schizophrenia Clinical Responses With Antipsychotics: A Quantitative System Pharmacology Study

CNS disorders are lagging behind other indications in implementing genotype-dependent treatment algorithms for personalized medicine. This report uses a biophysically realistic computer model of an associative and dorsal motor cortico-striatal-thalamo-cortical loop and a working memory cortical model to investigate the pharmacodynamic effects of COMTVal158Met rs4680, 5-HTTLPR rs 25531 s/L and D2DRTaq1A1 genotypes on the clinical response of 7 antipsychotics. The effect of the genotypes on dopamine and serotonin dynamics and the level of target exposure for the drugs was calibrated from PET displacement studies. The simulations suggest strong gene-gene pharmacodynamic interactions unique to each antipsychotic. For PANSS Total, the D2DRTaq1 allele has the biggest impact, followed by the 5-HTTLPR rs25531. The A2A2 genotype improved efficacy for all drugs, with a more complex outcome for the 5-HTTLPR rs25531 genotype. Maximal range in PANSS Total for all 27 individual combinations is 3 (aripiprazole) to 5 points (clozapine). The 5-HTTLPR L/L with aripiprazole and risperidone and the D2DRTaq1A2A2 allele with haloperidol, clozapine and quetiapine reduce the motor side-effects with opposite effects for the s/s genotype. The COMT genotype has a limited effect on antipsychotic effect and EPS. For cognition, the COMT MM 5-HTTLPR L/L genotype combination has the best performance for all antipsychotics, except clozapine. Maximal difference is 25% of the total dynamic range in a 2-back working memory task. Aripiprazole is the medication that is best suited for the largest number of genotype combinations (10) followed by Clozapine and risperidone (6), haloperidol and olanzapine (3) and quetiapine and paliperidone for one genotype. In principle, the platform could identify the best antipsychotic treatment balancing efficacy and side-effects for a specific individual genotype. Once the predictions of this platform are validated in a clinical setting the platform has potential to support rational personalized treatment guidance in clinical practice.

Exploring the relation between BOLD fMRI and cognitive performance using a computer-based quantitative systems pharmacology model: Applications to the COMTVAL158MET genotype and ketamine

BOLD fMRI is increasingly used mostly in an observational way to probe the effect of genotypes or therapeutic intervention in normal and diseased subjects. We use a mechanism-based quantitative systems pharmacology computer model of a human cortical microcircuit, previously calibrated for the 2-back working memory paradigm, adding established biophysical principles, of glucose metabolism, oxygen consumption, neurovascular effects and the paramagnetic impact on blood oxygen levels to calculate a readout for the voxel-based BOLD fMRI signal. The objective was to study the effect of the Catechol-O-methyl Transferase Val158Met (COMT) genotype on performance and BOLD fMRI. While the simulation suggests that on average virtual COMTVV genotype subjects perform worse, subjects with lower GABA, lower 5-HT₃ and higher 5-HT_1A activation can improve cognitive performance to the level of COMTMM subjects but at the expense of higher BOLD fMRI signal. In a schizophrenia condition, increased NMDA, GABA tone and noise levels, and lower D₁R activity can improve cognitive outcome with greater BOLD fMRI signal in COMT Val-carriers. We further generate hypotheses about why ketamine in healthy controls increases the BOLD fMRI signal but reduces cognitive performance. These simulations suggest a strong non-linear relationship between BOLD fMRI signal and cognitive performance. When validated, this mechanistic approach can be useful for moving beyond the descriptive nature of BOLD fMRI imaging and supporting the proper interpretation of imaging biomarkers in CNS disorders.

Simulating the Effects of Common Comedications and Genotypes on Alzheimer's Cognitive Trajectory Using a Quantitative Systems Pharmacology Approach

BACKGROUND

Many Alzheimer's disease patients in clinical practice are on polypharmacy for treatment of comorbidities.

OBJECTIVE

While pharmacokinetic interactions between drugs have been relatively well established with corresponding treatment guidelines, many medications and common genotype variants also affect central brain circuits involved in cognitive trajectory, leading to complex pharmacodynamic interactions and a large variability in clinical trials.

METHODS

We applied a mechanism-based and ADAS-Cog calibrated Quantitative Systems Pharmacology biophysical model of neuronal circuits relevant for cognition in Alzheimer's disease, to standard-of-care cholinergic therapy with COMTVal158Met, 5-HTTLPR rs25531, and APOE genotypes and with benzodiazepines, antidepressants, and antipsychotics, all together 9,585 combinations.

RESULTS

The model predicts a variability of up to 14 points on ADAS-Cog at baseline (COMTVV 5-HTTLPRss APOE 4/4 combination is worst) and a four-fold range for the rate of progression. The progression rate is inversely proportional to baseline ADAS-Cog. Antidepressants, benzodiazepines, first-generation more than second generation, and most antipsychotics with the exception of aripiprazole worsen the outcome when added to standard-of-care in mild cases. Low dose second-generation benzodiazepines revert the negative effects of risperidone and olanzapine, but only in mild stages. Non APOE4 carriers with a COMTMM and 5HTTLPRLL are predicted to have the best cognitive performance at baseline but deteriorate somewhat faster over time. However, this effect is significantly modulated by comedications.

CONCLUSION

Once these simulations are validated, the platform can in principle provide optimal treatment guidance in clinical practice at an individual patient level, identify negative pharmacodynamic interactions with novel targets and address protocol amendments in clinical trials.

Neuronal Circuit-Based Computer Modeling as a Phenotypic Strategy for CNS R&D

With the success rate of drugs for CNS indications at an all-time low, new approaches are needed to turn the tide of failed clinical trials. This paper reviews the history of CNS drug Discovery over the last 60 years and proposes a new paradigm based on the lessons learned. The initial wave of successful therapeutics discovered using careful clinical observations was followed by an emphasis on a phenotypic target-agnostic approach, often leading to successful drugs with a rich pharmacology. The subsequent introduction of molecular biology and the focus on a target-driven strategy has largely dominated drug discovery efforts over the last 30 years, but has not increased the probability of success, because these highly selective molecules are unlikely to address the complex pathological phenotypes of most CNS disorders. In many cases, reliance on preclinical animal models has lacked robust translational power. We argue that Quantitative Systems Pharmacology (QSP), a mechanism-based computer model of biological processes informed by preclinical knowledge and enhanced by neuroimaging and clinical data could be a new powerful knowledge generator engine and paradigm for rational polypharmacy. Progress in the academic discipline of computational neurosciences, allows one to model the effect of pathology and therapeutic interventions on neuronal circuit firing activity that can relate to clinical phenotypes, driven by complex properties of specific brain region activation states. The model is validated by optimizing the correlation between relevant emergent properties of these neuronal circuits and historical clinical and imaging datasets. A rationally designed polypharmacy target profile will be discovered using reverse engineering and sensitivity analysis. Small molecules will be identified using a combination of Artificial Intelligence methods and computational modeling, tested subsequently in heterologous cellular systems with human targets. Animal models will be used to establish target engagement and for ADME-Tox, with the QSP approach complemented by in vivo preclinical models that can be further refined to increase predictive validity. The QSP platform can also mitigate the variability in clinical trials with the concept of virtual patients. Because the QSP platform integrates knowledge from a wide variety of sources in an actionable simulation, it offers the possibility of substantially improving the success rate of CNS R&D programs while, at the same time, reducing both cost and the number of animals.

Predicting parkinsonism side-effects of antipsychotic polypharmacy prescribed in secondary mental healthcare

BACKGROUND

Computer-modelling approaches have the potential to predict the interactions between different antipsychotics and provide guidance for polypharmacy.

AIMS

To evaluate the accuracy of the quantitative systems pharmacology platform to predict parkinsonism side-effects in patients prescribed antipsychotic polypharmacy.

METHODS

Using anonymized data from South London and Maudsley NHS Foundation Trust electronic health records we applied quantitative systems pharmacology, a neurophysiology-based computer model of humanized neuronal circuits, to predict the risk for parkinsonism symptoms in patients with schizophrenia prescribed two concomitant antipsychotics. The performance of the quantitative systems pharmacology model was compared with the performance of simple parameters such as: combination of affinity constants (1/K_sum); sum of D₂R occupancies (D₂R) and chlorpromazine equivalent dose.

RESULTS

We identified 832 patients with schizophrenia who were receiving two antipsychotics for six or more months, between 1 January 2007 and 31 December 2014. The area under the receiver operating characteristic (AUROC) curve for the quantitative systems pharmacology model was 0.66 ( p = 0.01), while AUROCs for D₂R, 1/K_sum and chlorpromazine equivalent dose were 0.52 ( p = 0.350), 0.53 ( p = 0.347) and 0.52 ( p = 0.330) respectively.

CONCLUSION

Our results indicate that quantitative systems pharmacology has the potential to predict the risk of parkinsonism associated with antipsychotic polypharmacy from minimal source information, and thus might have potential decision-support applicability in clinical settings.

Impact of amyloid-beta changes on cognitive outcomes in Alzheimer's disease: analysis of clinical trials using a quantitative systems pharmacology model

BACKGROUND

Despite a tremendous amount of information on the role of amyloid in Alzheimer's disease (AD), almost all clinical trials testing this hypothesis have failed to generate clinically relevant cognitive effects.

METHODS

We present an advanced mechanism-based and biophysically realistic quantitative systems pharmacology computer model of an Alzheimer-type neuronal cortical network that has been calibrated with Alzheimer Disease Assessment Scale, cognitive subscale (ADAS-Cog) readouts from historical clinical trials and simulated the differential impact of amyloid-beta (Aβ40 and Aβ42) oligomers on glutamate and nicotinic neurotransmission.

RESULTS

Preclinical data suggest a beneficial effect of shorter Aβ forms within a limited dose range. Such a beneficial effect of Aβ40 on glutamate neurotransmission in human patients is absolutely necessary to reproduce clinical data on the ADAS-Cog in minimal cognitive impairment (MCI) patients with and without amyloid load, the effect of APOE genotype effect on the slope of the cognitive trajectory over time in placebo AD patients and higher sensitivity to cholinergic manipulation with scopolamine associated with higher Aβ in MCI subjects. We further derive a relationship between units of Aβ load in our model and the standard uptake value ratio from amyloid imaging. When introducing the documented clinical pharmacodynamic effects on Aβ levels for various amyloid-related clinical interventions in patients with low Aβ baseline, the platform predicts an overall significant worsening for passive vaccination with solanezumab, beta-secretase inhibitor verubecestat and gamma-secretase inhibitor semagacestat. In contrast, all three interventions improved cognition in subjects with moderate to high baseline Aβ levels, with verubecestat anticipated to have the greatest effect (around ADAS-Cog value 1.5 points), solanezumab the lowest (0.8 ADAS-Cog value points) and semagacestat in between. This could explain the success of many amyloid interventions in transgene animals with an artificial high level of Aβ, but not in AD patients with a large variability of amyloid loads.

CONCLUSIONS

If these predictions are confirmed in post-hoc analyses of failed clinical amyloid-modulating trials, one should question the rationale behind testing these interventions in early and prodromal subjects with low or zero amyloid load.

Towards the virtual human patient. Quantitative Systems Pharmacology in Alzheimer's disease

Development of successful therapeutic interventions in Central Nervous Systems (CNS) disorders is a daunting challenge with a low success rate. Probable reasons include the lack of translation from preclinical animal models, the individual variability of many pathological processes converging upon the same clinical phenotype, the pharmacodynamical interaction of various comedications and last but not least the complexity of the human brain. This paper argues for a re-engineering of the pharmaceutical CNS Research & Development strategy using ideas focused on advanced computer modeling and simulation from adjacent engineering-based industries. We provide examples that such a Quantitative Systems Pharmacology approach based on computer simulation of biological processes and that combines the best of preclinical research with actual clinical outcomes can enhance translation to the clinical situation. We will expand upon (1) the need to go from Big Data to Smart Data and develop predictive and quantitative algorithms that are actionable for the pharma industry, (2) using this platform as a "knowledge machine" that captures community-wide expertise in an active hypothesis-testing approach, (3) learning from failed clinical trials and (4) the need to go beyond simple linear hypotheses and embrace complex non-linear hypotheses. We will propose a strategy for applying these concepts to the substantial individual variability of AD patient subgroups and the treatment of neuropsychiatric problems in AD. Quantitative Systems Pharmacology is a new 'humanized' tool for supporting drug discovery and development in general and CNS disorders in particular.

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.

A Humanized Clinically Calibrated Quantitative Systems Pharmacology Model for Hypokinetic Motor Symptoms in Parkinson's Disease

The current treatment of Parkinson’s disease with dopamine-centric approaches such as L-DOPA and dopamine agonists, although very successful, is in need of alternative treatment strategies, both in terms of disease modification and symptom management. Various non-dopaminergic treatment approaches did not result in a clear clinical benefit, despite showing a clear effect in preclinical animal models. In addition, polypharmacy is common, sometimes leading to unintended effects on non-motor cognitive and psychiatric symptoms. To explore novel targets for symptomatic treatment and possible synergistic pharmacodynamic effects between different drugs, we developed a computer-based Quantitative Systems Pharmacology (QSP) platform of the closed cortico-striatal-thalamic-cortical basal ganglia loop of the dorsal motor circuit. This mechanism-based simulation platform is based on the known neuro-anatomy and neurophysiology of the basal ganglia and explicitly incorporates domain expertise in a formalized way. The calculated beta/gamma power ratio of the local field potential in the subthalamic nucleus correlates well (R² = 0.71) with clinically observed extra-pyramidal symptoms triggered by antipsychotics during schizophrenia treatment (43 drug-dose combinations). When incorporating Parkinsonian (PD) pathology and reported compensatory changes, the computer model suggests a major increase in b/g ratio (corresponding to bradykinesia and rigidity) from a dopamine depletion of 70% onward. The correlation between the outcome of the QSP model and the reported changes in UPDRS III Motor Part for 22 placebo-normalized drug-dose combinations is R² = 0.84. The model also correctly recapitulates the lack of clinical benefit for perampanel, MK-0567 and flupirtine and offers a hypothesis for the translational disconnect. Finally, using human PET imaging studies with placebo response, the computer model predicts well the placebo response for chronic treatment, but not for acute treatment in PD.

Assessing the synergy between cholinomimetics and memantine as augmentation therapy in cognitive impairment in schizophrenia. A virtual human patient trial using quantitative systems pharmacology

While many drug discovery research programs aim to develop highly selective clinical candidates, their clinical success is limited because of the complex non-linear interactions of human brain neuronal circuits. Therefore, a rational approach for identifying appropriate synergistic multipharmacology and validating optimal target combinations is desperately needed. A mechanism-based Quantitative Systems Pharmacology (QSP) computer-based modeling platform that combines biophysically realistic preclinical neurophysiology and neuropharmacology with clinical information is a possible solution. This paper reports the application of such a model for Cognitive Impairment In Schizophrenia (CIAS), where the cholinomimetics galantamine and donepezil are combined with memantine and with different antipsychotics and smoking in a virtual human patient experiment. The results suggest that cholinomimetics added to antipsychotics have a modest effect on cognition in CIAS in non-smoking patients with haloperidol and risperidone and to a lesser extent with olanzapine and aripiprazole. Smoking reduces the effect of cholinomimetics with aripiprazole and olanzapine, but enhances the effect in haloperidol and risperidone. Adding memantine to antipsychotics improves cognition except with quetiapine, an effect enhanced with smoking. Combining cholinomimetics, antipsychotics and memantine in general shows an additive effect, except for a negative interaction with aripiprazole and quetiapine and a synergistic effect with olanzapine and haloperidol in non-smokers and haloperidol in smokers. The complex interaction of cholinomimetics with memantine, antipsychotics and smoking can be quantitatively studied using mechanism-based advanced computer modeling. QSP modeling of virtual human patients can possibly generate useful insights on the non-linear interactions of multipharmacology drugs and support complex CNS R&D projects in cognition in search of synergistic polypharmacy.

Multitarget drug discovery projects in CNS diseases: quantitative systems pharmacology as a possible path forward

Clinical development in brain diseases has one of the lowest success rates in the pharmaceutical industry, and many promising rationally designed single-target R&D projects fail in expensive Phase III trials. By contrast, successful older CNS drugs do have a rich pharmacology. This article will provide arguments suggesting that highly selective single-target drugs are not sufficiently powerful to restore complex neuronal circuit homeostasis. A rationally designed multitarget project can be derisked by dialing in an additional symptomatic treatment effect on top of a disease modification target. Alternatively, we expand upon a hypothetical workflow example using a humanized computer-based quantitative systems pharmacology platform. The hope is that incorporating rationally multipharmacology drug discovery could potentially lead to more impactful polypharmacy drugs.

Understanding responder neurobiology in schizophrenia using a quantitative systems pharmacology model: application to iloperidone

The concept of targeted therapies remains a holy grail for the pharmaceutical drug industry for identifying responder populations or new drug targets. Here we provide quantitative systems pharmacology as an alternative to the more traditional approach of retrospective responder pharmacogenomics analysis and applied this to the case of iloperidone in schizophrenia. This approach implements the actual neurophysiological effect of genotypes in a computer-based biophysically realistic model of human neuronal circuits, is parameterized with human imaging and pathology, and is calibrated by clinical data. We keep the drug pharmacology constant, but allowed the biological model coupling values to fluctuate in a restricted range around their calibrated values, thereby simulating random genetic mutations and representing variability in patient response. Using hypothesis-free Design of Experiments methods the dopamine D4 R-AMPA (receptor-alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid) receptor coupling in cortical neurons was found to drive the beneficial effect of iloperidone, likely corresponding to the rs2513265 upstream of the GRIA4 gene identified in a traditional pharmacogenomics analysis. The serotonin 5-HT3 receptor-mediated effect on interneuron gamma-aminobutyric acid conductance was identified as the process that moderately drove the differentiation of iloperidone versus ziprasidone. This paper suggests that reverse-engineered quantitative systems pharmacology is a powerful alternative tool to characterize the underlying neurobiology of a responder population and possibly identifying new targets.

Quantitative Systems Pharmacology can reduce attrition and improve productivity in pharmaceutical research and development

The empirical hypothesis generation and testing approach to pharmaceutical research and development (R&D), and biomedical research has proven very effective over the last half-century; resulting in tremendous increases productivity and the rates of approval for new drug applications at the Food and Drug Administration (FDA). However, as discovery of new therapeutic approaches for diseases with unmet medical need becomes more challenging, the productivity and efficiency of the traditional approach to drug discovery and development is diminishing. Innovative approaches are needed, such as those offered by Quantitative Systems Pharmacology (QSP) modeling and simulation. This “systems” approach to modeling and simulation can be used to guide the hypothesis generation and testing process in pharmaceutical R&D, in a manner similar to its adoption in other industries in the past. Embedding QSP into the existing processes within pharmaceutical discovery and development will be required in order to realize the full beneficial impact of this innovative approach.

A computer-based quantitative systems pharmacology model of negative symptoms in schizophrenia: exploring glycine modulation of excitation-inhibition balance

Although many antipsychotics can reasonably control positive symptoms in schizophrenia, patients' return to society is often hindered by negative symptoms and cognitive deficits. As an alternative to animal rodent models that are often not very predictive for the clinical situation, we developed a new computer-based mechanistic modeling approach. This Quantitative Systems Pharmacology approach combines preclinical basic neurophysiology of a biophysically realistic neuronal ventromedial cortical-ventral striatal network identified from human imaging studies that are associated with negative symptoms. Calibration of a few biological coupling parameters using a retrospective clinical database of 34 drug-dose combinations resulted in correlation coefficients greater than 0.60, while a robust quantitative prediction of a number of independent trials was observed. We then simulated the effect of glycine modulation on the anticipated clinical outcomes. The quantitative biochemistry of glycine interaction with the different NMDA-NR₂ subunits, neurodevelopmental trajectory of the NMDA-NR_2B in the human schizophrenia pathology, their specific localization on excitatory vs. inhibitory interneurons and the electrogenic nature of the glycine transporter resulted in an inverse U-shape dose-response with an optimum in the low micromolar glycine concentration. Quantitative systems pharmacology based computer modeling of complex humanized brain circuits is a powerful alternative approach to explain the non-monotonic dose-response observed in past clinical trial outcomes with sarcosine, D-cycloserine, glycine, or D-serine or with glycine transporter inhibitors. In general it can be helpful to better understand the human neurophysiology of negative symptoms, especially with targets that show non-monotonic dose-responses.

Dopaminergic contributions to hippocampal pathophysiology in schizophrenia: a computational study

Since the original formulation of the dopamine hypothesis, a number of other cellular-level abnormalities--eg, NMDA receptor hypofunction, GABA system dysfunction, neural connectivity disturbances--have been identified in schizophrenia, but the manner in which these potentially interact with hyperdopaminergia to lead to schizophrenic symptomatology remains uncertain. Previously, we created a neuroanatomically detailed, biophysically realistic computational model of hippocampus in the control (unaffected) and schizophrenic conditions, implemented on a 72-processor supercomputer platform. In the current study, we apply the effects of dopamine (DA), dose-dependently, to both models on the basis of an exhaustive review of the neurophysiologic literature on DA's ion channel and synaptic level effects. To index schizophrenic behavior, we use the specific inability of the model to attune to the 40 Hz (gamma band) frequency, a finding that has been well replicated in the clinical electroencephalography (EEG) and magnetoencephalography literature. In trials using 20 'simulated patients', we find that DA applied to the control model produces modest increases in 40 Hz activity, similar to experimental studies. However, in the schizophrenic model, increasing DA induces a decrement in 40 Hz resonance. This modeling work is significant in that it suggests that DA's effects may vary based on the neural substrate on which it acts, and--via simulated EEG recordings-points to the neurophysiologic mechanisms by which this may occur. We also feel that it makes a methodological contribution, as it exhibits a process by which a large amount of neurobiological data can be integrated to run pharmacologically relevant in silico experiments, using a systems biology approach.