Exemplary Care Cut t ing-edge Research World-class Education
CRISMA Laboratory UPMC Critical Care www.ccm.pitt.edu Exemplary - - PowerPoint PPT Presentation
CRISMA Laboratory UPMC Critical Care www.ccm.pitt.edu Exemplary - - PowerPoint PPT Presentation
CRISMA Laboratory UPMC Critical Care www.ccm.pitt.edu Exemplary Care Cut t ing-edge Research World-class Education Multi-scale modeling of Influenza a virus and its containment Towards System Biology Workshop Grenoble, May 2011
Exemplary Care Cut t ing-edge Research World-class Education
Multi-scale modeling of Influenza a virus and its containment
Gilles Clermont, MD
CRISMA Center Center for Inflammation and Regenerative Modeling Critical Care Medicine, Mathematics, Industrial Engineering University of Pittsburgh
Towards System Biology Workshop Grenoble, May 2011
Exemplary Care Cut t ing-edge Research World-class Education
Acknowledgments
Colleagues Trainees
Baris Hancioglu Ian Price David Dreisigmeyer Sarah Lukens Jay Depasse
NIH(GM)/NSF(DMS)
Exemplary Care Cut t ing-edge Research World-class Education
Epidemic Influenza
In a “normal” season
Up to 30 million infections in the US 100,000 attributable hospitalizations 35,000 attributable deaths 3B (direct) -15B (total) in societal costs
Most deaths are due to secondary infections (pneumonia) Target population (old, sick and the very young)
Exemplary Care Cut t ing-edge Research World-class Education
Pandemic 2009 H1N1
In the US
61M cases, 12,700 deaths
Severe illness
The Australia-New Zealand experience (NEJM)
722 cases (29 cases/M) 14.3% mortality
The Canadian experience (JAMA)
215 cases 14.3% mortality
Our own experience Pittsburgh/UPMC
24 admissions to adult ICU starting in august 2009 75% mortality
Target population
Young adult, pregnant women
Exemplary Care Cut t ing-edge Research World-class Education
Exemplary Care Cut t ing-edge Research World-class Education
Influenza A virus
Exemplary Care Cut t ing-edge Research World-class Education
Pandemic influenza is an inflammatory disease
Perrone et al., PLoS Pathogens 2008
1918 H5N1 Tx/91 H5N1
Exemplary Care Cut t ing-edge Research World-class Education
Influenza modeling in Pittsburgh
Population (Burke) Host (Clermont/Rosenfeld) Biodefense (Morel) Epidemiology Data Center JHU Genetics (Ghedin-Virology) Vaccines-In vivo (Ross - CVR) PSC Social Sciences Allegheny county HD
Exemplary Care Cut t ing-edge Research World-class Education
Multiscale modeling
IVA Global pool Human IVA Population level Human IVA Host level Containment strategies
Exemplary Care Cut t ing-edge Research World-class Education
Multiscale modeling – MIDAS
Exemplary Care Cut t ing-edge Research World-class Education
Multiscale modeling – Population
Exemplary Care Cut t ing-edge Research World-class Education
Population-level models (SIR)
Susceptibles
P(Infection|Environment)
Infected
P(Duration) ?? P(Death)
Recovered
P(Susceptible)?=0
Resistant
P(immune)
- Human activity (travel)
- Geography/weather
Exemplary Care Cut t ing-edge Research World-class Education
Population-level models
Agent-based model (FRED) Sophisticated cellular automata Stochastic - Distribution driven Age, household size, school sizes… Transmission, duration of disease… Computational requirements 32 nodes 75GB (3e8 agents) 5 hours
Ferguson et al. Nature 2005
Exemplary Care Cut t ing-edge Research World-class Education
Multiscale modeling – Viral evolution
Exemplary Care Cut t ing-edge Research World-class Education
Multiscale modeling - Host
Host Lung Cells Pathways Tissues Virus Therapy
Vaccines Anti-virals
- M2 inhibitors
- NA inhibitors
Exemplary Care Cut t ing-edge Research World-class Education
The simplest viral model – ever ! (v1.0)
Virus Infected cell Healthy cell Killer cell Λ=0 D=0 δ - where active cell killing is implicit in this death rate p β c - virus death
- r shedding
dH HV dt dI HV I dt dV pI cV dt β β δ = − = − = −
Exemplary Care Cut t ing-edge Research World-class Education
IVA-host model 2.0
Hancioglu et al, J Theor Biol 2007
Exemplary Care Cut t ing-edge Research World-class Education
Host-level model 3.0
Price et al. JCC 2008
Exemplary Care Cut t ing-edge Research World-class Education
Looking at individuals
Baccam et al. Virology 2007
Exemplary Care Cut t ing-edge Research World-class Education
The estimation problem
Virus Infected cell Healthy cell Killer cell Λ=0 D=0 δ - where active cell killing is implicit in this death rate p β c - virus death
- r shedding
dH HV dt dI HV I dt dV pI cV dt β β δ = − = − = −
Exemplary Care Cut t ing-edge Research World-class Education
Population ensemble models
E(Mn)
E(Mn) ≡ Metamodel or Ensemble
Many many more models than individuals
Exemplary Care Cut t ing-edge Research World-class Education
Creating the ensemble
Error Probability (density) Parameter axis pA pB pC pD pE
( )
2
( | )
i D H H
L L
M D e
β − Γ +Γ
∝ L
Exemplary Care Cut t ing-edge Research World-class Education
Calculating the Ensemble
Earl, Deem 2006
Local minima Parameter space
MH-type algorithm Parallel tempering/ Replica exchange as mixing method
Exemplary Care Cut t ing-edge Research World-class Education
Chain mixing
1
β
2
β
3
β
4
β
5
β
tests
Exemplary Care Cut t ing-edge Research World-class Education
The population ensemble
Exemplary Care Cut t ing-edge Research World-class Education
Looking at a population
Baccam et al. Virology 2007
Exemplary Care Cut t ing-edge Research World-class Education
Probabilistic ensembles – for subpopulations
Sparsity -> pooling “similar” patients may be good enough The clinician does this
E(M)
Exemplary Care Cut t ing-edge Research World-class Education
A population ensemble
Patient 1 Patient 4
Exemplary Care Cut t ing-edge Research World-class Education
Host-level models
Exemplary Care Cut t ing-edge Research World-class Education
V3.0 ensemble
Exemplary Care Cut t ing-edge Research World-class Education
Probabilistic ensembles – for subpopulations
Sparsity -> pooling “similar” patients may be good enough The clinician does this
E(M)
Exemplary Care Cut t ing-edge Research World-class Education
Multiobjective estimation
2 2 a b
E E E = +
a b
E E E
α β
=
Exemplary Care Cut t ing-edge Research World-class Education
Linking scales of description
7 14 Threshold of upper respiratory symptoms
% Infected cells
10 50
B
Dmax
Viral load (log particles/cc) 7 14 12 8 Threshold of infectivity
A
Vtdur Vtmin Vmax
Time (days) 7 14 Threshold of systemic symptoms
IFN-alpha level
100
500 C
Fmin Fdur Ddur
Time (days)
Exemplary Care Cut t ing-edge Research World-class Education
Linking scales of description
Transmission = infectivity x upper respiratory symptoms Behavior (stay home or not)= systemic symptoms
Viral load (log particles/cc) 7 14 12 8 Threshold of infectivity
A
Vtdur Vtmin Vtmin Vtdur
Exemplary Care Cut t ing-edge Research World-class Education
Linking scales of description
Integrated run (explicit hybrid) = not practical Look-up tables = fast but memory hungry Response surfaces = forward simulations, polynomial fit, algebraic computation
{ }
{ }
Transmissibility Behavior Resistance Transmissibility Behavior Resistance
(Inputs) , , , , , Inputs Viral load, burst size, viral adhesion, existing antibodies, .. R f R µ µ µ σ σ σ = ∈ ∈ Inputs R
Exemplary Care Cut t ing-edge Research World-class Education
Predicted impact of therapy
Virulent pathogen Si ifi i i i i
Exemplary Care Cut t ing-edge Research World-class Education
A multiscale lung model
Inflammation occurs in the tissue barrier between air and blood. Tissue swelling impairs gas
- diffusion. Extreme inflammation
- f a respiratory unit (~25 alveoli)
can completely stop gas exchange (shunt). The global impact of inflammation depends on the combined contribution of respiratory units (RU) with diverse anatomical and physiologic properties.
Reynolds et al, J Theor Biol 2009
Exemplary Care Cut t ing-edge Research World-class Education
Knowledge and successful translation
Empirical Epicycles
- to XVII century
Kepler/Newton Ellipses First interaction between a physical law and empiric
- bservation
GR - Einstein Precessing ellipses Discrepancy between predictions and empiric observation
Depth of knowledge
Discrepancy between gravity and other forces of nature
QG - ?? Black hole physics
Exemplary Care Cut t ing-edge Research World-class Education
UPMC Critical Care
www.ccm.pitt.edu
Exemplary Care Cut t ing-edge Research World-class Education