Tutorials › Quantitative Systems Pharmacology (QSP) › Tracing Causal Paths to Set the Mechanistic Boundary

Defining Model Scope and Boundaries · Tutorial 17 of 799

Tracing Causal Paths to Set the Mechanistic Boundary

Trace causal pathways from intervention to outcome to identify the mechanisms a decision-focused QSP model must represent—and biology it may leave outside its boundary.

Intermediate 12 min read

What You'll Learn

  • Define a causal path from an intervention to a decision-relevant outcome.
  • Distinguish direct routes, mediated routes, and branches that do not reach the outcome.
  • Use a causal diagram to propose a mechanistic model boundary.
  • Calculate intervention contrasts along a simple pathway model.
  • Recognize when context variables or a different outcome change the boundary.
  • Identify assumptions and failure modes in path-based boundary setting.

Why Trace Causal Paths?

A biological system can contain far more processes than a particular QSP decision requires. The model boundary is not a catalogue of everything known about the disease or drug. It is a proposed scope for representing the causal processes needed to connect an intervention to the outcome that informs a decision. Causal path tracing makes that scope explicit: start at the intervention, follow plausible causal links, and ask which processes can transmit its effect to the decision-relevant outcome.

The outcome and decision context must be stated before drawing the boundary. For example, a pathway irrelevant to an efficacy endpoint may matter to a safety endpoint, and a mechanism important for one patient population may be inactive in another. The earlier material on decision questions and model output boundaries establishes why that context matters. Here, the focus is the causal diagram used to trace routes within it.

Definition: A causal path from an intervention \(D\) to an outcome \(Y\) is a sequence of directed causal links beginning at \(D\) and ending at \(Y\). Variables or processes along such a sequence are on a path from \(D\) to \(Y\). A mechanistic boundary is the set of processes represented to describe the causal routes needed for the specified intervention-to-outcome question, together with any additional variables needed to define those routes in the stated context.

A causal diagram uses nodes for variables or processes and arrows for hypothesized causal influence. A direct path is a single arrow, such as \(D \rightarrow Y\). A mediated path includes one or more intervening nodes, such as \(D \rightarrow M \rightarrow Y\). Several paths may run in parallel. A node affected by the intervention but with no directed route onward to the specified outcome is a branch that does not transmit an effect to that outcome in the diagram.

The arrows are causal claims, not merely observed associations. A diagram organizes a proposed explanation; it does not establish that the links are true. Its use assumes that the intervention is sufficiently well defined, the relevant context and outcome are specified, the temporal ordering of causes is coherent, and the diagram includes the plausible causal routes needed for the question. If important causes or routes are omitted, a boundary drawn from the diagram may be inadequate.

From a Diagram to a Boundary

Begin at the intervention and trace forward along every plausible directed route that reaches the decision-relevant outcome. The mechanisms on those routes are candidates for representation. Then consider whether other variables are needed to specify how those routes operate—for example, a context variable that changes a mediator’s level or a cause of both a mediator and the outcome. Such a variable is not necessarily a drug mechanism, but omitting it may prevent the model from representing the intervention contrast appropriately.

Do not equate “not on a directed path from the intervention” with “always irrelevant.” A variable may be a common cause of a pathway component and the outcome, may define the population, or may be relevant to another decision outcome. Conversely, a biologically interesting process is not automatically required merely because it changes after treatment. The question is whether it can transmit an effect to the stated outcome under the causal structure and context being modeled.

1
Fix the question.
State the intervention or contrast, the population and context, and the outcome used for the decision.
2
Draw plausible causal links.
Separate causal hypotheses from associations and order the nodes consistently with the proposed biology.
3
Trace routes to the outcome.
Mark every plausible directed route from the intervention to the outcome, including parallel routes and intermediate processes.
4
Check context and omitted causes.
Ask what other variables are needed to specify or interpret the routes in the stated population and intervention setting.
5
Set a provisional boundary.
Include the required route processes; document branches excluded for this outcome and context rather than treating them as universally unimportant.

Why Path Coefficients Matter

A simple linear structural model illustrates how paths carry an intervention effect. Suppose \(M\) is a mediator and \(Y\) is the outcome, with equations \(M=aD+U\) and \(Y=bM+cD+V\), where \(U\) and \(V\) are fixed for the intervention contrast. Substituting the first equation into the second gives \(Y=abD+bU+cD+V\). Therefore, changing \(D\) by one unit changes the predicted \(Y\) by \(ab+c\): \(ab\) is the contribution along \(D \rightarrow M \rightarrow Y\), and \(c\) is the contribution along the direct arrow \(D \rightarrow Y\). This calculation relies on the stated linear equations and fixed context; it is not a general formula for nonlinear or confounded systems.

More generally, in a linear acyclic structural model, a directed path contributes the product of its link coefficients to the total intervention effect, and the total effect is the sum of contributions across directed paths. This follows by repeatedly substituting each node’s equation into the outcome equation: each substitution generates terms for the routes through that node, and terms reaching the outcome correspond to directed paths. The result is a way to check a specified model’s implications, not evidence that the diagram or coefficient values are correct.

Worked Example: A Mediated Effect With a Context Variable

Consider a hypothetical model of an intervention \(D\), a biological mediator \(M\), a decision-relevant score \(Y\), and a context variable \(U\). Let \(D=0\) denote control and \(D=1\) denote intervention. The proposed structural equations are \(M=8-2D+U\) and \(Y=100+3M-5U\). The diagram has \(D \rightarrow M \rightarrow Y\), \(U \rightarrow M\), and \(U \rightarrow Y\). Calculate the predicted scores and intervention contrast at each context value \(U=0\) and \(U=1\).

At \(U=0\), control gives \(M=8-2(0)+0=8\), so \(Y=100+3(8)-5(0)=124\). Intervention gives \(M=8-2(1)+0=6\), so \(Y=100+3(6)-5(0)=118\). The contrast is \(118-124=-6\) score units.

At \(U=1\), control gives \(M=8-2(0)+1=9\), so \(Y=100+3(9)-5(1)=100+27-5=122\). Intervention gives \(M=8-2(1)+1=7\), so \(Y=100+3(7)-5(1)=100+21-5=116\). The contrast is \(116-122=-6\) score units.

In both contexts the intervention contrast is \((-2)(3)=-6\) score units: the intervention changes \(M\) by \(-2\) units, and each unit of \(M\) changes \(Y\) by \(3\) score units. The direct \(U\)-to-\(Y\) term shifts the scores with context but cancels in the within-context contrast. In this example \(U\) is needed to predict absolute scores, though it does not alter the contrast. The equations and values are hypothetical and do not establish a biological mechanism.

Parallel Routes and Excluded Branches

Path tracing is especially useful when a drug may affect more than one process. If two plausible routes reach the outcome, excluding either route can change the predicted intervention effect. A branch that changes under treatment but does not lead to the selected outcome may be outside the boundary for that question. The conclusion is specific: if the branch becomes relevant to another endpoint, the boundary must be reconsidered.

Worked Example: Two Routes to an Outcome

In a hypothetical linear model, intervention \(D\) affects a target state \(T\) and a second process \(O\), with both influencing an inflammation score \(Y\). Let \(T=10-2D\), \(O=4-D\), and \(Y=50-3T-2O\). Assume these equations correctly describe the specified context and that no other routes affect the intervention contrast. Calculate the contrast and identify the route contributions.

Under control, \(D=0\), so \(T=10\), \(O=4\), and \(Y=50-3(10)-2(4)=12\). Under intervention, \(D=1\), so \(T=8\), \(O=3\), and \(Y=50-3(8)-2(3)=20\). Thus the intervention contrast is \(20-12=8\) inflammation-score units.

Along the target route, \(D\) changes \(T\) by \(-2\), and the coefficient from \(T\) to \(Y\) is \(-3\). Its contribution is \((-2)(-3)=6\) score units. Along the second route, \(D\) changes \(O\) by \(-1\), and the coefficient from \(O\) to \(Y\) is \(-2\). Its contribution is \((-1)(-2)=2\) score units. The contributions sum to \(6+2=8\), matching the direct calculation. A boundary that included only \(T\) would predict a contrast of \(6\), not \(8\), under these assumed equations.

Worked Example: A Treatment-Affected Branch Outside the Outcome Boundary

Suppose a hypothetical intervention changes an immune-cell marker \(B\), while a separate mediator \(M\) changes the selected efficacy outcome \(Y\). The proposed equations are \(B=5+2D\), \(M=10-3D\), and \(Y=40-2M\). The diagram contains \(D \rightarrow B\), \(D \rightarrow M \rightarrow Y\), and no arrow from \(B\) to \(Y\) or to \(M\). Compare the predicted efficacy outcome under control and intervention.

For \(D=0\), \(M=10\), so \(Y=40-2(10)=20\). For \(D=1\), \(M=7\), so \(Y=40-2(7)=26\). The intervention contrast is \(26-20=6\) outcome units. Although \(B\) changes from \(5\) to \(7\), it does not enter the equation for \(Y\); including it would not change this efficacy prediction under the stated diagram.

This exclusion is conditional on the question and causal structure. If \(B\) were itself a safety outcome, or if evidence supported \(B \rightarrow Y\), it would no longer be an irrelevant branch for the corresponding model purpose. A missing arrow cannot be justified merely because the available data do not show an association.

Diagnostics, Alternatives, and Limitations

A path-based boundary is a hypothesis to scrutinize, not a guarantee of adequacy. Warning signs include a claimed intervention effect with no route to the modeled outcome, a route that depends on an unsupported causal link, or a process that changes after treatment but has been excluded without considering whether it affects the chosen endpoint. Also check the direction and timing of each proposed link. An observed association between two nodes does not establish which causes which, and a biomarker can track an outcome without mediating the intervention’s effect.

The boundary should also be checked against the model’s intended use. A single-outcome boundary may be narrower than one needed for a joint efficacy-and-safety decision. A population-specific route may fail in another population. If feedback is important, a simple acyclic diagram at one time point may not represent it; time-indexed states or another explicit representation may be needed. These are reasons to revisit the proposed structure, not reasons to include all biology by default.

A broad mechanistic network can represent many candidate routes, but it increases the number of assumptions and quantities that must be supported. A reduced mechanistic model focuses on routes needed for the decision, while retaining causal interpretation for those routes. An empirical exposure–response fit may summarize a relationship in the data without representing the intervening biological pathways; it can be useful for interpolation within its supported setting, but may not answer a counterfactual that depends on changing a mechanism. These approaches answer different questions, and no diagram alone determines which is adequate.

Common pitfall: “The intervention changes this biomarker, so the biomarker must be inside the model.” A change after intervention establishes neither that the biomarker lies on a causal path to the selected outcome nor that it is needed to predict the decision. Trace a supported route to the specified outcome; if none exists, the biomarker may be outside that outcome’s boundary, while remaining relevant to another endpoint or question.

Path tracing has limits. The diagram may omit a route, misstate a causal direction, or oversimplify interactions and context dependence. A process can influence a path without appearing as a mediator on the path itself, and nonlinear effects can make a single set of path contributions inadequate. Thus, use the diagram to make scope and assumptions inspectable, then compare its implied predictions with biological knowledge and available evidence. A defensible boundary is explicit about what is represented, what is excluded, and for which outcome and context that choice is being made.

Check Your Understanding

Answer each question, then compare with the solutions.

  1. A diagram contains \(D \rightarrow A \rightarrow Y\) and \(D \rightarrow B\), with no route from \(B\) to \(Y\). Is \(B\) on a directed path from \(D\) to \(Y\)?
  2. In the linear model \(M=5- D\), \(Y=20+4M+2D\), calculate the predicted \(Y\) at \(D=0\) and \(D=1\), and give the intervention contrast.
  3. Why does an observed association between a biomarker and an outcome not by itself establish that the biomarker mediates the intervention effect?
  4. A process has no route to the efficacy outcome but is itself a decision-relevant safety endpoint. Should it automatically be excluded from all models? Explain.
  5. In a diagram with two directed routes from \(D\) to \(Y\), what is the risk of omitting one route from a model intended to predict the intervention contrast?

Solutions. (1) No. The route ends at \(B\) and does not continue to \(Y\) in the stated diagram. (2) At \(D=0\), \(M=5\) and \(Y=20+4(5)+2(0)=40\). At \(D=1\), \(M=4\) and \(Y=20+4(4)+2(1)=38\). The contrast is \(38-40=-2\) outcome units; the mediated contribution is \((-1)(4)=-4\), and the direct contribution is \(+2\), summing to \(-2\). (3) Association alone does not establish causal direction or rule out common causes; mediation requires a supported causal route through the biomarker. (4) No. Its lack of a path to efficacy does not make it irrelevant to a safety decision, so the model boundary depends on the endpoint and use. (5) The omitted route may contribute to the total intervention effect, leading to a biased or incomplete prediction under the proposed causal structure.

Clinical Biostats