CRISMA Laboratory UPMC Critical Care www.ccm.pitt.edu Exemplary - - PowerPoint PPT Presentation

crisma laboratory upmc critical care
SMART_READER_LITE
LIVE PREVIEW

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


slide-1
SLIDE 1

 Exemplary Care  Cut t ing-edge Research  World-class Education 

CRISMA Laboratory UPMC Critical Care

www.ccm.pitt.edu

slide-2
SLIDE 2

 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

slide-3
SLIDE 3

 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)

slide-4
SLIDE 4

 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)

slide-5
SLIDE 5

 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

slide-6
SLIDE 6

 Exemplary Care  Cut t ing-edge Research  World-class Education 

slide-7
SLIDE 7

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Influenza A virus

slide-8
SLIDE 8

 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

slide-9
SLIDE 9

 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

slide-10
SLIDE 10

 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

slide-11
SLIDE 11

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Multiscale modeling – MIDAS

slide-12
SLIDE 12

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Multiscale modeling – Population

slide-13
SLIDE 13

 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
slide-14
SLIDE 14

 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

slide-15
SLIDE 15

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Multiscale modeling – Viral evolution

slide-16
SLIDE 16

 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
slide-17
SLIDE 17

 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 β β δ = − = − = −

slide-18
SLIDE 18

 Exemplary Care  Cut t ing-edge Research  World-class Education 

IVA-host model 2.0

Hancioglu et al, J Theor Biol 2007

slide-19
SLIDE 19

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Host-level model 3.0

Price et al. JCC 2008

slide-20
SLIDE 20

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Looking at individuals

Baccam et al. Virology 2007

slide-21
SLIDE 21

 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 β β δ = − = − = −

slide-22
SLIDE 22

 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

slide-23
SLIDE 23

 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

slide-24
SLIDE 24

 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

slide-25
SLIDE 25

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Chain mixing

1

β

2

β

3

β

4

β

5

β

tests

slide-26
SLIDE 26

 Exemplary Care  Cut t ing-edge Research  World-class Education 

The population ensemble

slide-27
SLIDE 27

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Looking at a population

Baccam et al. Virology 2007

slide-28
SLIDE 28

 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)

slide-29
SLIDE 29

 Exemplary Care  Cut t ing-edge Research  World-class Education 

A population ensemble

Patient 1 Patient 4

slide-30
SLIDE 30

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Host-level models

slide-31
SLIDE 31

 Exemplary Care  Cut t ing-edge Research  World-class Education 

V3.0 ensemble

slide-32
SLIDE 32

 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)

slide-33
SLIDE 33

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Multiobjective estimation

2 2 a b

E E E = +

a b

E E E

α β

=

slide-34
SLIDE 34

 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)

slide-35
SLIDE 35

 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

slide-36
SLIDE 36

 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

slide-37
SLIDE 37

 Exemplary Care  Cut t ing-edge Research  World-class Education 

Predicted impact of therapy

Virulent pathogen Si ifi i i i i

slide-38
SLIDE 38

 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

slide-39
SLIDE 39

 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

slide-40
SLIDE 40

 Exemplary Care  Cut t ing-edge Research  World-class Education 

UPMC Critical Care

www.ccm.pitt.edu

slide-41
SLIDE 41

 Exemplary Care  Cut t ing-edge Research  World-class Education 