In silico prediction of amyloid-related imaging abnormalities

Abstract

The dosage of the active ingredient in a treatment provided to or under consideration for providing to a subject is specified. The level of local amyloid beta is predicted based on the dosage of the active ingredient. T 2 The severity of amyloid-related imaging abnormalities (ARIA) (ARIA-E) that appear as high signal in fluid-attenuated inversion recovery (FLAIR) images is predicted based on the predicted amount of removal of local amyloid beta, and predicting the severity of ARIA-E includes predicting the degree of blood vessel wall damage. A result corresponding to the predicted ARIA-E severity is output.

Description

[Technical Field]

[0001] Cross-reference of related applications This application claims the benefit and priority of U.S. Provisional Application No. 63 / 337,992, filed on 3 May 2022, which is incorporated herein by reference in its entirety for all purposes. [Background technology]

[0002] Over the past several decades, amyloid-beta (Aβ) aggregates have been the most common target in clinical trials investigating disease-modifying therapies in Alzheimer's disease. Several therapies based on monoclonal antibodies against Aβ aggregates have reported amyloid-associated imaging abnormalities (ARIAs). Two types of ARIA have been identified by magnetic resonance imaging (MRI) of the brain. ARIA-E, appearing as high signal intensity on T2-weighted fluid-suppressed inversion-recovery (FLAIR) images, suggests vasogenic edema, sulcal effusion, and gyral swelling. ARIA-H, appearing as low signal intensity on T2*-weighted gradient echo sequences, is thought to represent hemosiderosis and microbleeds. ARIA-E is transient and typically resolves within a few months, while ARIA-H remains visible on subsequent MRIs. Most ARIA-E events are mild to moderate in severity and usually unrelated to symptoms, as assessed by a variety of radiological measures that show good correlation.

[0003] The risk of treatment-related ARIA is dose-dependent and appears to increase in APOEε4 carriers. To mitigate the risk of ARIA, several anti-amyloid clinical trials have employed dose-escalation schemes relative to the target dose, as well as routine brain MRI monitoring. The incidence of treatment-related ARIA events across various anti-amyloid clinical trials has been reported to range from 10% to 42%, with the highest proportion occurring in APOEε4 allele carriers. Spontaneous ARIA events were rarely observed in placebo groups (e.g., <3% ARIA-E) and outside the environment of anti-amyloid clinical trials. An example of the latter is ARIA-like events occurring in cerebral amyloid angiopathy-associated inflammation (CAA-ri), a rare autoimmune response to vascular Aβ aggregates. In clinical practice, ARIA can be a significant burden for Alzheimer's patients and healthcare systems due to the MRI monitoring requirements. Furthermore, ARIA can limit patients' ability to reach their target dose and benefit from amyloid-reducing disease-modifying therapy.

[0004] Therefore, it would be advantageous to be able to generate predictions of treatment-related ARIA events and their severity. [Brief explanation of the drawing]

[0005] This disclosure will be explained in conjunction with the attached drawings. [Figure 1] Figure 1 shows the main pathophysiological pathways linking drug-mediated elimination of Aβ aggregates to ARIA. [Figure 2] Figure 2 shows an exemplary ARIA-E prediction network, including the ARIA-E prediction system and user devices. [Figure 3] Figure 3 shows a vascular wall damage (VWD) model and related case studies, along with plots of changes over time and variable relationships. [Figure 4] Figure 4 shows a flowchart of an exemplary process 400 for predicting ARIA-E according to several embodiments of the present invention. [Figure 5]Figure 5 shows an exemplary scatter plot of residuals generated by a PK model used to predict plasma drug concentrations for ARIA-E cases from SR / MR OLE studies. Individually weighted residuals (IWRES) are shown. Solid lines represent empirical (10th, 50th, and 90th) percentiles, and dotted lines represent predicted (10th, 50th, and 90th) percentiles. [Figure 6] Figure 6 shows a scatter plot of exemplary residuals generated by Part 1 (upper panel) and Part 2 (lower panel) of a two-part PD model used to predict the probability of ARIA-E occurrence and the magnitude of ARIA-E, respectively. Normalized predictive distribution error (NPDE) and IWRES are also shown. Solid lines represent empirical (10th, 50th, and 90th) percentiles, and dotted lines represent predicted (10th, 50th, and 90th) percentiles. [Figure 7] Figure 7 shows exemplary model performance data at the target level. [Figure 8A] Figure 8A shows data illustrating the effect of model parameters on amyloid-related imaging abnormalities (ARIA-E) due to edema / exudation dynamics. [Figure 8B] Figure 8B shows data illustrating the effect of model parameters on amyloid-related imaging abnormalities (ARIA-E) due to edema / exudation dynamics. [Figure 8C] Figure 8C shows data illustrating the effect of model parameters on amyloid-related imaging abnormalities (ARIA-E) due to edema / exudation dynamics. [Overview of the project]

[0006] In some embodiments, a computer implementation method is provided. This method includes the dosage of an active ingredient in a treatment provided to or under consideration for provision to a subject. The method further includes predicting the level of local amyloid-beta based on the dosage of the active ingredient and predicting the severity of amyloid-related imaging abnormalities (ARIA) (ARIA-E) that appear as high signals on T2-weighted fluid-suppressed inversion recovery (FLAIR) images based on the predicted amount of local amyloid-beta removed, wherein predicting the severity of ARIA-E includes predicting the degree of vascular wall damage. The method also includes outputting results corresponding to the predicted ARIA-E severity.

[0007] The active ingredient may include an anti-amyloid monoclonal antibody.

[0008] This method may include using a pharmacokinetic model to predict the time course of the active ingredient concentration in plasma during treatment, and the prediction of local amyloid beta is based on at least one predicted concentration of the active ingredient over the predicted time course.

[0009] Predicting the level of topical amyloid beta may include estimating a baseline level of topical amyloid beta, calculating the rate of topical amyloid beta removal based on the dose of the active ingredient and the baseline level of topical amyloid beta, and predicting the said level of topical amyloid beta based on the calculated rate of topical amyloid beta removal.

[0010] Predicting the aforementioned levels of local amyloid beta may involve solving a pharmacodynamic differential equation that assumes the rate of change of amyloid beta is proportional to the product of the concentration of the active ingredient in the subject and the local amyloid beta level.

[0011] Predicting the level of vascular wall damage may involve solving differential equations that include the level of focal amyloid-beta.

[0012] Predicting the severity of ARIA-E may involve solving an algebraic equation that assumes a non-linear relationship between the level of vascular wall damage and the Barkhof Global Total Score (BGTS).

[0013] The method may include identifying a potential schedule for monitoring ARIA events based on the predicted severity of ARIA, and the results characterize the potential schedule.

[0014] The method may include identifying a potential recommended dosage of an active ingredient for a subject based on the predicted severity of ARIA, and the results characterize the potential recommended dosage.

[0015] In some embodiments, the system includes a non-transitory computer-readable storage medium that, when executed by one or more data processors, includes instructions that cause the one or more data processors to execute some or all of one or more of the methods and / or some or all of one or more of the processes disclosed herein. Some embodiments of the present disclosure include a computer program product tangibly embodied in a non-transitory machine-readable storage medium that includes instructions configured to cause one or more data processors to execute some or all of one or more of the methods and / or some or all of one or more of the processes disclosed herein.

[0016] The terms and expressions used are used as terms of description and not of limitation, and there is no intention to exclude equivalents of the features shown and described or portions thereof. However, it is recognized that various modifications are possible within the scope of the claimed systems and methods. Accordingly, although the systems and methods are specifically disclosed by way of example and any features, modifications and variations of the concepts disclosed herein should be recognized by those skilled in the art, and such modifications and variations are considered to be within the scope of the systems and methods defined by the appended claims.

[0017] This summary is not intended to identify any significant or essential features of the claimed subject matter, nor is it intended to be used solely to determine the scope of the claimed subject matter. These exemplary examples are mentioned not to limit or define the disclosure, but to provide examples to aid its understanding. Further embodiments and examples are discussed in the detailed description, where further explanation is provided. The subject matter should be understood by referring to the appropriate parts of the entire specification of this disclosure, any or all of the drawings, and each claim.

[0018] The above, along with other features and embodiments, will become more apparent upon reference to the following specification, claims, and accompanying drawings. [Modes for carrying out the invention]

[0019] The precise pathophysiological mechanisms of ARIA are not fully understood. It has been proposed that drug-mediated removal of Aβ aggregates increases the permeability of the cerebral vascular walls to the entry of fluids and blood products into the brain, leading to ARIA-E and ARIA-H, respectively. In the brains of Alzheimer's patients, Aβ aggregates are found in the brain parenchyma as well as in the cerebral vascular walls in the form of cerebral amyloid angiopathy (CAA). Anti-amyloid antibodies associated with ARIA bind to various forms of Aβ aggregates (e.g., oligomers, fibrils, and plaques) and remove them by degradation and / or effector cell-mediated phagocytosis.

[0020] Figure 1 shows the major pathophysiological pathways thought to link drug-mediated removal of Aβ aggregates to ARIA. Overall, all hypothetical pathogenesis of ARIA appears to involve some degree of damage to the cerebral vascular wall. The extent to which the degree of amyloid loading, or the location and rate of amyloid clearance, influences the onset, time course, and severity of ARIA remains unclear. A better understanding of ARIA etiology could reveal additional risk factors that could accelerate pathways toward patient-centered strategies for medication and monitoring in individuals with Alzheimer's disease.

[0021] Some embodiments of the present invention relate to predicting the temporal dynamics and dose-dependence of ARIA-E severity by using a semi-mechanical in silico model of ARIA-E called a vascular wall injury (VWD) model. For example, the model may receive input data containing dose information indicating the amount of anti-amyloid monoclonal antibody (e.g., gantenerumab) administered in a simulation according to a given schedule (e.g., the model simulates the administration of doses every four weeks), and the model may generate results corresponding to the predicted occurrence and / or severity of ARIA (or ARIA-E) at each of a set of time points relative to a reference point in the schedule (e.g., administration of a first dose). The predicted severity can be used to generate potential monitoring schedules, administration schedules, and / or dose adjustment criteria. For example, a schedule for monitoring and safety data collection may be proposed based on an initial dosing plan, and dose adjustment criteria may be defined to propose reducing the dose by a specified amount only if MRI shows at least a specified degree of ARIA-E.

[0022] The VWD model framework, consistent with the ARIA-E pathological mechanism proposed in Figure 1, allows for different biological interpretations of the modeled mechanism. The VWD model links several pharmacological and biological factors to the observed pharmacokinetics and magnitude of ARIA-E. One biological factor considered in the model is focal amyloid, either vascular or parenchymal, whose drug-mediated removal can trigger a cascade of events leading to ARIA-E in one or more regions of the brain. Another biological variable in the model is the level of VWD, which is a hypothetical measure of susceptibility to fluid leakage into the brain resulting from the destruction of vascular integrity and / or perivascular inflammation induced by amyloid-drug interactions. The VWD model configuration allows for the assessment of the potential interaction between drug-mediated removal of focal amyloid and intrinsic vascular repair processes, which can determine or influence the level of VWD, which in turn influences the magnitude of ARIA-E. ARIA-E prediction network

[0023] Figure 2 shows the ARIA-E prediction network including the ARIA-E prediction system 205 and the user device 210. The ARIA-E prediction system 205 consists of pharmacokinetic model components (controlled by the pharmacokinetic model controller 215) and ( Pharmacodynamics The VWD model, which includes both pharmacokinetic model components (controlled by the model controller 220), is controlled by the model controller 220. The VWD model (as shown on the left side of Figure 3) may be a semi-mechanical pharmacokinetic pharmacodynamic model that integrates known and / or estimated processes and / or (for example, simultaneously) accounts for inter-individual variability of model parameters.

[0024] Each of the pharmacokinetic model controllers 215 and 220 may be configured to specify the formulas used in the model, retrieve stored data (e.g., parameter values) from memory, and execute their respective model components. The aggregate model controller 225 may call the pharmacokinetic model controller 215 and / or the pharmacodynamic model controller 220, send input data, and / or receive output. Thus, the aggregate model controller 225 may be configured to control the definition and execution of the VWD model by interacting with the pharmacokinetic model controller 215 and / or the pharmacodynamic model controller 220. The aggregate model controller 225 may further perform pre-processing and / or post-processing. For example, such pre-processing and / or post-processing may adjust for subject-specific characteristics such as age, APOE e4 genotype, baseline ARIA-H load, and vascular risk factors (which may alternatively or additionally be factors used by the pharmacokinetic and / or pharmacodynamic models to generate predictions).

[0025] The pharmacokinetic model components are configured to predict the time course of the concentration of the active ingredient in the target plasma based on administration history and long-term observation of plasma concentrations. An exemplary model, including a pharmacokinetic model, may be disclosed in Retout S et al., "Disease Modeling and Model-Based Meta-Analyses to Define a New Direction for a Phase III Program of Gantenerumab in Alzheimer's Disease," Clin Pharmacol Ther. 2022 Apr;111(4):857-866, which is incorporated herein by reference in its entirety for all purposes.

[0026] The top graph in the central panel of Figure 3 shows an exemplary time course of the concentration of the active ingredient in the target plasma (left axis) as estimated using a pharmacokinetic model, when the active ingredient is administered at the dose indicated by the vertical position of the right axis plus sign and at the time indicated by the horizontal position of the plus sign plus sign.

[0027] The pharmacodynamic model components are configured to predict ARIA-E findings based on the concentration of the active ingredient, based on intermediate calculations that predict changes in the levels of local amyloid-β and vascular wall damage caused by the active ingredient. The intermediate calculations may include solving one or more differential equations and / or performing one or more algebraic calculations. Exemplary equations that may be used in the pharmacodynamic model components are listed below.

[0028] The first exemplary differential equation states that the rate of drug-mediated removal of local Aβ is in plasma

number

number

number

[0029] In particular, this formula refers to the localized level of amyloid, which represents a hypothetical measure of either vascular amyloid or parenchymal amyloid from brain regions affected by ARIA-E. The rate of removal of this “localized amyloid” is thought to be proportional to both the drug concentration in plasma and the existing localized amyloid level. In this specification, α 除去The proportionality constant, known as the removal constant, is shown and can be interpreted as the rate constant of the biochemical reaction between the anti-amyloid antibody and the amyloid aggregates. Such a reaction corresponds to the drug-mediated amyloid removal that occurs, and its magnitude may vary between individuals. The baseline level of local amyloid (indicated herein as amyloid 0) affects the initial rate of amyloid removal and may also vary between individuals. Other processes that can alter local amyloid levels, such as physiological amyloid production and endogenous amyloid clearance, are not included in the model.

[0030] The second graph from the top in the central panel of Figure 3 shows the estimated time course of local Aβ, estimated using equation (1). In this particular example, local amyloid is estimated to decrease over time.

[0031] The second exemplary differential equation assumes that the removal of local Aβ drives the accumulation of VWD(t) [arbitrary units], where this is the repair rate k 修復 [day -1 (This process is offset by the primary vascular repair process characterized by...)

number

[0032] Parameters: amyloid d0, α 除去 and k 修復 As shown in Figure 7, this value is allowed to vary depending on the subject, and its statistical distribution (typical value and range) is estimated. One, several, or all of the remaining model parameters may be constant across individuals.

[0033] The third graph from the top in the central panel of Figure 3 shows the estimated time course of VWD, estimated using equation (2). In this particular case, VWD is estimated to increase first and then decrease overall.

[0034] The "VWD" variable corresponds to a hypothetical concept describing susceptibility to fluid leakage into the brain due to treatment-induced defects in cerebral blood vessels and perivascular structures, such as impaired blood-brain barrier and burdened intramural periarterial drainage pathways. In the absence of interference, VWD is equal to 0. The accumulation of VWD is promoted by drug-mediated amyloid removal and offset by the presumed intrinsic vascular repair process. With respect to equation (2), vascular repair is given by the rate constant (k 修復 It is assumed that the process is a primary process with a half-life independent of the size of the VWD. The size of the VWD may be related to the predicted ARIA-E score (BGTS) by a sigmoid function (for example). This relationship can guarantee zero or small BGTS at low levels of VWD, a sharp increase in the ARIA-E score at intermediate levels, and impose an upper limit on the size of the ARIA-E (BGTS = 60) at high levels of VWD.

[0035] Two additional exemplary algebraic equations relate VWD to the magnitude of observed ARIA-E. Using two equations for this relationship addresses the complexity of having a large number of zero BGTS values ​​in the ARIA-E score dataset; the first equation predicts the probability of an ARIA-E event (i.e., a non-zero BGTS) for each BGTS observation, and the second equation quantifies the magnitude of the ARIA-E when the event is predicted to occur with a high probability.

[0036] The first exemplary algebraic equation deals with the probability p(t) of a positive (non-zero) BGTS under the assumption that the logarithm of the odds ratio is related linearly to VWD(t) as follows:

number

number

[0037] The graph below the central panel in Figure 3 shows the estimated time course of BGTS for a typical case related to the other graphs in the central panel. The estimated time course of BGTS is determined using equations (3) and (4). In this typical case, BGTS increases immediately after the start of the treatment regimen and then decreases overall.

[0038] The aggregate model controller 225 can use the predicted BGTS to facilitate the determination of proposed dosing schedules and / or proposed monitoring schedules for individuals experiencing ARIA-E. For example, the aggregate model controller can adjust the generation of predicted BGTS time courses for various prescribed schedules for dosing, various prescribed schedules for monitoring, and / or various different BGTS cutoff values ​​(a threshold for ARIA-E severity that leads to treatment interruption in the model). Thus, for each prescribed dosing schedule, monitoring schedule, and / or BGTS cutoff value, a predicted progression of ARIA-E severity (BGTS time course) can be generated, which can then be output (e.g., via interface controller 235) to assist the user's (e.g., clinician) decision regarding which dosing schedule, monitoring schedule, and / or BGTS cutoff value to use in clinical trials and / or real-world use cases.

[0039] Therefore, the VWD model can interconnect the following factors: (i) the dose of gantenerumab, (ii) the drug concentration in plasma, (iii) drug-induced removal of local amyloid, (iv) VWD estimated to result in fluid leakage into the brain, (v) the intrinsic vascular repair process, and (vi) BGTS, which quantifies the severity of ARIA-E. Individual-specific information regarding a given dose and administration time can be provided as input, while parameters for some or all of the other factors can be estimated.

[0040] The parameter collector 230 first identifies the parameters of each of the pharmacokinetic and pharmacodynamic model components. In some examples, the parameter collector 230 identifies the values ​​for one or more parameters by querying an external source for the parameters. In some examples, the parameters collector230 learns the values ​​of one or more parameters by fitting the model components or VWD model using one or more training datasets. For example, the training dataset may include results from studies in which a set of subjects (e.g., each diagnosed with Alzheimer's disease) are administered a given active ingredient (or an active ingredient corresponding to a particular class of active ingredients), e.g., gantenerumab, in one or more specific doses according to one or more dosing schedules. The dosing schedule may represent the relative times at which various doses of the active ingredient are administered, e.g., an initial dose, a recent dose, or another baseline. In some examples, the dose of the active ingredient administered at each specific dosing time in the schedule is the same as for each of the other dosing times. In some examples, the dose of the active ingredient administered at at least two different dosing times. For example, the dose may increase over at least two doses. The study may have been conducted in which MRI was performed periodically and evaluated to monitor ARIA-E, and samples (e.g., blood samples) were collected periodically and evaluated to determine the concentration of the active ingredient. The assessment may include assigning a severity level to any ARIA-E detected using scoring techniques such as the Barkhof total score.

[0041] Research data from data sources may be further filtered to selectively include data corresponding to a subset of subjects, so that the training dataset includes selective data. The subset of subjects may include subjects in whom an occurrence of ARIA-E (e.g., corresponding to a severity greater than a given threshold, such as a severity greater than zero) was observed during the study or within one year of the first dose of the study.

[0042] In some examples, the parameter collector 230 identifies the parameters of a single model component in the isolation of other model components. For example, the parameters of a pharmacokinetic model component may be defined by performing a Bayesian analysis on the study data.

[0043] In some examples, parameter collector 230 fixes some parameters while adapting other parameters. For example, parameter collector 230 may fix the parameters of the pharmacokinetic model component while adapting the pharmacodynamic model component. Parameter collector 230 may estimate parameter values using a non-linear mixed effects modeling technique.

[0044] Parameter collector 230 can use assumptions regarding the distribution of variables (e.g., the shape of the distribution), maximum values, whether various variables are dependent on each other, etc. For example, the parameter Collector 230 may define the assumption that parameters amyloid0, α 除去 and k 修復 have log-normal inter-individual variability, estimate the fixed effects of their typical values, and estimate the width of the distribution (random effect). As another example, parameter collector 230 can define the estimated values of parameters β1, VWD 50 and pow without inter-individual variability (i.e., as only fixed effects). It can be shown that VWD(t) is proportional to amyloid0, meaning that parameters amyloid0 and EG 50 appear effectively only as the ratio of the expressions of BGTS(t). To ensure structural identifiability, parameter collector 230 can set EG 50 to 1 without loss of generality. Parameter collector 230 can fix the maximum value of BGTS to 60 along the limits of the 60-point BGTS severity scale. The residual (unexplained) variability that describes the difference between the observed BGTS and the predicted BGTS can be assumed to be independent of the prediction and normally distributed.

[0045] The ARIA-E prediction system 205 includes an interface controller 235 that facilitates communication between the ARIA-E prediction system 205 and the user device 210. The interface controller 235 may be configured to process incoming communications from the user device 210 to detect a request for prediction generated by (e.g.) an aggregate model controller (e.g., based on the output from the pharmacokinetic model controller 215 and / or pharmacodynamic model controller 220), and / or to detect information to be used as input to the prediction (e.g., identifying the dosage and / or administration schedule used to treat a given subject or being considered for the treatment of a given subject).

[0046] The interface controller 235 may further, or alternatively, be configured to generate code that can initiate the generation of a presentation or interface (e.g., a web page) that can be read by a software application and may contain or represent one or more predictions generated by the ARIA-E prediction system 205. For example, the interface controller 235 may generate web page code that, when executed by a browser, presents a web page (on a website) that identifies one or more predictions. The website and / or web page may initially operate to collect input on the user device 210 used to run the model.

[0047] When the user device 210 receives information and / or codes from the ARIA-E prediction system 205, the user device 210 may generate and present a presentation containing information and / or data that specifies that a code should be presented. Thus, the user device 210 may present data corresponding to a prediction about whether or to what extent a given treatment regimen can result in ARIA-E. ARIA-E prediction process

[0048] Figure 4 shows a flowchart of an exemplary process 400 for predicting ARIA-E according to several embodiments of the present invention. In block 402, the aggregate model controller 225 identifies the actual or potential dose of the active ingredient. Identification may include identifying the actual or potential dose that has been or should be administered at a given time in the treatment regimen (e.g., 0 hours, 8 weeks after the start of the regimen, at the point when a certain dose is reached). The dose may be (e.g.) the dose indicated in a communication received from the user device 210 (e.g., as identified by the interface controller 235), a default dose (e.g., as identified and / or retrieved by the aggregate model controller 225), one or more possible doses (e.g., as identified and / or retrieved by the aggregate model controller 225), etc.

[0049] In block 404, the pharmacodynamic model controller 220 predicts the level of local amyloid. More specifically, the pharmacokinetic model controller 215 The pharmacodynamic model controller 220 can predict the time course of the active ingredient concentration based on the dose of the active ingredient, and predicts the time course of local amyloid levels based on the predicted active ingredient concentration. The predicted level of local amyloid can be predicted simultaneously with the administration of the dose, predicted at a specific date compared to when a first (or given) dose was administered, and can be predicted for each of several time points compared to when a first (or given) dose was administered, and so on.

[0050] In block 406, the pharmacodynamic model controller 220 predicts the level of vascular wall damage based on the predicted level of local amyloid and predicted vascular repair. The predicted level of vascular wall damage can be generated by solving differential equations such as equation (2). The predicted level of vascular wall damage may further depend on the predicted repair rate (for example).

[0051] In block 408, the pharmacodynamic model controller 220 predicts ARIA-E severity once or multiple times (e.g., over time of ARIA-E severity). Block 408 may include predicting the level of local amyloid beta at one point in time or multiple time points (e.g., over time of level), and / or using the predicted level. Block 406 may include predicting the predicted level of vascular wall damage at one point in time or multiple time points (e.g., over time of vascular wall damage), and / or using the predicted level. The local amyloid beta level may be predicted based on the estimated active ingredient concentration level using one or more differential equations. The degree of vascular wall damage may be predicted based on the predicted local amyloid beta level (e.g., using a method such as shown with reference to block 406), using one or more differential equations.

[0052] In block 410, results corresponding to the predicted ARIA-E severity level can be output. For example, an ARIA-E prediction system 205 The results can be transmitted to the user device 210. The results may include one or more predicted ARIA levels (e.g., time course of ARIA-E levels), a potential schedule for MRI-based monitoring of ARIA-E (e.g., by BGTS scoring), or a potential revised dosing schedule (e.g., showing one or more proposed doses of the active ingredient and / or the relative timing between one or more consecutive doses). Examples Dynamic behavior of illustrative examples

[0053] The data used for model development consist of long-term measurements of plasma drug concentrations and ARIA-E size from a subset (N=112) of individuals who developed ARIA-E during the open-label extension (OLE) of the Scarlet Road (SR; NCT01224106) and Marguerite Road (MR; NCT02051608) studies of gantenerumab, including participants from the studies for prodromal Alzheimer's disease and mild Alzheimer's disease, respectively. In short, SR / MR trial participants who underwent double-blind treatment and had at least one follow-up visit were eligible to participate in the open-label extension (OLE). During the SR / MR OLE trial, all participants received subcutaneous gantenerumab every four weeks, with dose escalation towards a target dose of 1200 mg. Each participant was assigned to one of five dose-setting regimens based on their APOE4 carrier status and last dose during double-blind treatment (Table 1). Protocol-defined routine MRI monitoring for ARIA-E detection was performed at regular intervals. If ARIA-E was detected, the dosage was adjusted based on its radiological severity and / or the presence of symptoms. The size of ARIA-E on MRI scans was assessed using the Barkhof Total Gradual Score (BGTS), a 60-point severity scale where higher scores indicate greater severity. The reported BGTS described the number and size of parenchymal hyperintensity, sulcal hyperintensity, and gyral swellings present on FLAIR images. [Table 1]

[0054] The four ARIA-E cases shown in Figures 3 and 7 were selected to highlight diverse ARIA-E features (size, resolution, relapse) that are mechanistically insightful but do not reflect the typical ARIA-E profile of gantenerumab. Each trial participant was assigned to one of the five dose-setting regimens listed in Table 1. Because the model accounted for inter-individual variability, the modeled time-profile for each ARIA-E case was generated based on empirical Bayesian estimates, i.e., the most likely values ​​of the individual parameters.

[0055] VWD parameters were estimated by fitting the entire VWD model to a set of long-term BGTS observations while keeping individual pharmacokinetic parameters fixed. Boolean data from the first part (e.g., yes / no for ARIA-E occurrence) and continuous observations from the second part (e.g., BGTS values ​​within the range of ARIA-E magnitude) were implemented in the data file as two different observation types and modeled simultaneously. Parameter estimation was performed using the nonlinear mixed-effects method in Monolix (version 2020R1).

[0056] The data shown in Figure 5 demonstrate that this previously constructed population pharmacokinetic model was able to describe individual PK data from the SR / MR OLE trial, and the residuals are uniformly distributed around zero, with most values ​​within -2 and +2 standard deviations, indicating the absence of significant systematic bias.

[0057] Parameters: Amyloid 0, α 除去 and k 修復 Assuming log-normal inter-individual variability, the typical values ​​(fixed effects) and distribution width (random effects) were estimated. Parameters β1, VWD 50 And pow was estimated without inter-individual variability (i.e., as fixed effects only). It can be shown that VWD(t) is proportional to A amyloid 0, and the parameters amyloid 0 and EG 50 This means that it appears only as a ratio in the formula BGTS(t). Therefore, in order to ensure structural identifiability, without loss of generality, EG 50 I set it to 1. Furthermore, BGTS 最大The value was fixed at 60, along the limits of the 60-point BGTS severity scale. Residual (unexplained) variability describing the difference between observed BGTS and predicted BGTS was assumed to be independent of prediction and normally distributed. Estimated population parameters (representative and inter-individual variability) are shown in Table 2. Diagnostic plots from Figure 6 demonstrate the model's ability to predict the probability of ARIA-E events at the population level, as well as the magnitude of ARIA-E at the individual level in particular, with residuals over time evenly concentrated around zero, no significant systematic bias, and most values ​​within -2 and +2 standard deviations. [Table 2]

[0058] The estimated time course of active ingredient concentration was calculated using the pharmacokinetic model disclosed by Retout S et al., "Disease Modeling and Model-Based Meta-Analyses to Define a New Direction for a Phase III Program of Gantenerumab in Alzheimer's Disease," Clin Pharmacol Ther. 2022 Apr;111(4):857-866. The estimated time course of local amyloid beta was calculated using equation (1). The estimated time course of VWD was calculated using equation (2). The estimated time course of BGTS was calculated using equations (3) and (4).

[0059] The estimated time course of concentration, local amyloid beta, VWD, and BGTS is shown in the central panel of Figure 3 for Case 1, and on the left side of Figure 7 for Cases 2-4.

[0060] Cases 1 and 2 were assigned to the same rapid escalation regimen (regmen 4), yet they developed ARIA-E at significantly different times. (See the central panel of Figure 3 and the top set of graphs in Figure 7). Case 1 reached the target dose by week 16 and experienced its first ARIA-E event at week 20, while Case 2 developed ARIA-E at week 12 before reaching the target dose. These two cases appear to experience different rates of amyloid removal around week 12 (as can be seen from the different slopes of the focal amyloid curves). This is due to different estimates of the removal constant (α in the upper right plot of Figure 7). 除去 This is explained by the model (see values). This difference arises despite similar drug concentration levels and comparable estimates of baseline topical amyloid levels during the first 12 weeks.

[0061] Case 3, assigned to the slowest escalation regimen (Regimen 1), experienced two treatment interruptions due to ARIA-E (see the plus sign in the top graph of the central plot set in Figure 7). The onset of ARIA-E during slow escalation may be attributable to the relatively high baseline levels of local amyloid estimated for this subject (see the amyloid 0 value in the central right distribution plot in Figure 7). Model results suggest a recurrence of ARIA-E during readmission and dose-increase setting due to the combination of persistently high levels of local amyloid and relatively high drug concentrations, which drives another increase in VWD.

[0062] Case 4, assigned to the second slowest escalation regimen (regmen 2), developed an initial ARIA-E that was far larger than any of the ARIA-E events experienced by the other three cases (see the bottom graph in the bottom plot set of Figure 7). This model did not estimate significantly different levels of baseline topical amyloid or drug concentration that could explain the observed differences in the magnitude of ARIA-E. However, both the removal and repair constants estimated for Case 4 were significantly different from those estimated for the other three cases. Specifically, the significantly larger value of the removal constant corresponds to the smaller repair rate constant (α in Figure 7). 除去 and k修復 Local amyloid removal is driven at a high rate that cannot match the slow rate of vascular repair (as determined by the value). This imbalance results in large VWD levels and ultimately high BGTS.

[0063] The model predictions for the four ARIA-E cases share several features. The model predicted the sharpest decrease in focal amyloid levels around the time of the observed ARIA-E event, as well as the subsequent increase in the size of VWD and BGTS. This pattern is particularly pronounced at week 20 for case 1, week 12 for case 2, weeks 28 and 60 for case 3, and weeks 12 and 52 for case 4 (Figures 3 and 7). Throughout the entire treatment interruption period, the model predicted a plateau in focal amyloid levels and a gradual decrease in VWD. The latter also reflects the repair process, transitioning to a gradual decrease in the size of the ARIA-E that closely follows the observed ARIA-E resolution, which may be observed between weeks 20–36 for case 1, weeks 12–24 for case 2, weeks 28–44 for case 3, and weeks 12–40 for case 4 (Figures 3 and 7). Exploring the influence of key model parameters on ARIA-E dynamics

[0064] Considering that both baseline levels and removal constants of local amyloid affect the rate of local amyloid removal, their roles in ARIA-E kinetics were investigated separately in a series of simulations shown in Figures 8A and 8B. The control plot is shown in the central panel of Figure 3 and corresponds to Case 1. The other plots are simulations with a single variable modified (removal constant in Figure 8A, baseline local amyloid beta in Figure 8B, and repair rate constant in Figure 8C).

[0065] A larger removal constant was associated with an earlier increase in both VWD and ARIA-E, while a smaller removal constant delayed the onset of ARIA-E, but also delayed local amyloid clearance. Data points at weeks 12 and 24 provide comparison. For higher baseline levels of local amyloid, the model predicted early damage and recurrent ARIA-E events in the vascular wall due to persistently high levels of local amyloid. In contrast, for lower baseline levels of local amyloid, the model predicted sequential removal of local amyloid without significant ARIA-E events. Figure 8C shows the evolution of ARIA-E as the estimated vascular repair process is reduced or enhanced by variations in the repair rate constant. When the repair rate constant of the intrinsic vascular repair process is smaller, local amyloid removal leads to repeated damage to the vascular wall and ultimately to ARIA-E recurrence. Conversely, in this model, a vascular repair process with a larger rate constant can counteract the predicted amyloid removal rate at the end of the escalation and prevent significant ARIA-E. interpretation

[0066] This semi-mechanical in silico model addressed the ambiguity and complexity of the ARIA-E pathology mechanism in the context of long-term observation of ARIA-E size from anti-amyloid clinical trials. Previous ARIA-E models developed using bapineuzumab and gantenerumab data are event hazard models capable of predicting the probability of ARIA-E occurrence for various dosing regimens. The mathematical framework used in this embodiment allows for different biological interpretations of the modeled mechanism based on the following key factors: local amyloid loading, drug exposure (dose and drug concentration), amyloid removal, and vascular wall damage and repair. For the various ARIA-E scenarios shown in Figure 1, the same mathematical relationships describe drug-mediated removal of local amyloid, whether it is parenchymal or vascular Aβ loading being removed via direct degradation or cell-mediated phagocytosis. While the precise biological details underlying ARIA-E remain unclear, a common feature across the different mechanistic pathways that can trigger ARIA-E events appears to be damaged vascular walls. A model linking drug-mediated removal of local amyloid to VWD levels and ultimately to ARIA-E severity (as measured by BGTS) performed well in capturing inter-subject variability. Parameters were estimated using a nonlinear mixed-effects method by fitting the model to the BGTS dataset from the gantenerumab SR / MR OLE trial. The resulting VWD model provides a good explanation for the subject levels of ARIA-E reported in the gantenerumab SR / MR OLE trial, thereby validating the mechanism employed.

[0067] Beyond fitting the individual time courses of ARIA-E events, the VWD model offers several insights. According to this model, (i) high exposure (dose and drug concentration), (ii) high local amyloid and / or (iii) high efficiency in local amyloid removal (removal constant α(x) 除去Due to capture by , high rates of local amyloid removal may promote the development of ARIA-E. In the case of anti-amyloid antibodies such as aducanumab, donanemab, and gantenerumab, which are highly dependent on Fc receptor-mediated phagocytosis of amyloid aggregates by immune cells, higher estimates of the removal constant may reflect a larger, possibly excessive, response of immune cells activated by the amyloid-drug complex. Image-based analysis of a small number of ARIA-E cases from the gantenerumab study revealed that areas with FLAIR high signal often show a significant decrease in amyloid PET signal, thereby supporting the idea of ​​more efficient cell-mediated phagocytosis of Aβ aggregates at ARIA-E sites. Notably, significant amyloid reduction was also observed in areas without FLAIR abnormalities. Based on the VWD model, local amyloid reduction can occur without ARIA-E, insofar as the local rate of amyloid removal is offset by the rate of intrinsic vascular repair processes. Furthermore, according to this model, ARIA-E can completely dissipate before VWD levels return to zero (BGTS=0). This behavior is a result of the estimated nonlinear VWD-BGTS relationship, which ensures that minor damage to the vascular wall does not result in a detectable ARIA-E event, as shown in Figure 3. Finally, the VWD model suggests that recurrence of ARIA-E may be a combined effect of high drug concentrations and persistent high levels of local amyloid, and that amyloid depletion ultimately reduces the risk of ARIA-E over time, even at high drug concentrations.

[0068] When comparing model results with amyloid reduction findings in anti-amyloid clinical trials that reported ARIA, it is useful to distinguish between the global and local properties of the modeled / measured amyloid. For example, studies using comprehensive PET assessments of baseline amyloid in the bapineuzumab and gantenerumab studies showed no significant difference between the ARIA-E and non-ARIA-E groups. Instead, when amyloid PET signals were partially assessed, baseline amyloid load from the occipital brain region was found to be significantly higher in the ARIA-E group than in the non-ARIA-E group. Considering the eccentric distribution of CAA in the occipital vascular structure, higher baseline levels of local amyloid estimated by the VWD model may be associated with a higher local load of vascular amyloid.

[0069] Based solely on PET imaging, it is not possible to definitively say whether a decrease in amyloid PET signaling is due to the removal of vascular and / or parenchymal amyloid. The relative contribution of the two forms of Aβ aggregates to amyloid PET signaling is a subject of debate. While some studies have concluded that amyloid PET can detect vascular Aβ in CAA, other recent studies have reported that amyloid PET signaling appears to be primarily driven by parenchymal amyloid plaques and not significantly confounded by CAA. The latter finding is limited by the fact that, despite pathological confirmation of CAA in all subjects, only half of the subjects met the clinical diagnosis of high-probability CAA with a median number of two cerebral microhemorrhages. Notably, ten years ago, Alzheimer's Association Research Roundtable Workgroup recommended that individuals with more than four cerebral microhemorrhages at baseline should not be included in anti-amyloid clinical trials. Therefore, the contribution of vascular Aβ aggregates to PET signaling in clinically more severe CAA cases, particularly those enrolled in more recent anti-amyloid studies, remains undetermined. To further elucidate the relative risk of ARIA-E associated with local vascular amyloid burden and its removal, more analyses are needed to integrate long-term PET and MRI data from anti-amyloid studies and quantify the local distribution of amyloid signals, FLAIR hyperintensity, and CAA imaging markers. The VWD model can be improved using future biomarkers that can better distinguish between vascular and parenchymal amyloid burden, thereby predicting the evolution of ARIA-E in patients with varying degrees of AD pathology (based on parenchymal amyloid burden) and CAA (based on vascular amyloid burden).

[0070] The semi-mechanical VWD model provides a mathematical and semi-mechanical explanation for the temporal patterns observed in BGTS data from the gantenerumab study. The VWD model enables in silico individual-level investigation of biologically plausible variables and processes that influence the onset and severity of ARIA-E. Upon validation, the VWD model would be useful in generating hypotheses that can be tested in clinical studies for ARIA management. Using the VWD model, the effects of the rate of local amyloid removal and vascular repair on the temporal course of ARIA-E in individuals with a history of ARIA-E can be simulated. Such model-based predictions can assist in decisions regarding the continuation or reintroduction of treatment to minimize the risk of ARIA-E progression.

[0071] While the subject matter has been described in detail with respect to its specific embodiments, those skilled in the art will understand that, upon achieving the above understanding, modifications, variations, and equivalents to such embodiments will readily be generated. Therefore, it should be understood that this disclosure is presented for illustrative purposes only, and not as an limitation, and does not preclude such modifications, variations, and / or additional inclusions to the subject matter that would readily become apparent to those skilled in the art. Indeed, the methods and systems described herein may be embodied in various other forms. Furthermore, various omissions, substitutions, and modifications may be made in the forms of the methods and systems described herein without departing from the spirit of this disclosure. The appended claims and their equivalents are intended to encompass forms or modifications that fall within the scope and spirit of this disclosure.

[0072] Unless otherwise specified, discussions throughout this specification using terms such as “processing,” “computing,” “calculating,” “determining,” and “identifying” are understood to refer to the operation or process of one or more computers or similar electronic computing devices or devices that manipulate or transform data represented as physical electronic or magnetic quantities in the memory, registers, or other information storage devices, transmission devices, or display devices of a computing platform.

[0073] The one or more systems described herein are not limited to any particular hardware architecture or configuration. A computing device may include any suitable arrangement of components that provide a result conditional on one or more inputs. Suitable computing devices range from general-purpose computing devices to dedicated computing devices implementing one or more embodiments of this subject matter, and include multipurpose microprocessor-based computing systems that access stored software for programming or configuring computing systems. Any suitable programming, scripting, or other type of language or combination of languages ​​may be used to implement the teachings contained herein in software used for programming or configuring computing devices.

[0074] Embodiments of the methods disclosed herein may be implemented in the operation of such computing devices. The order of the blocks shown in the above examples may be changed, and for example, blocks may be rearranged, combined, and / or divided into subblocks. Certain blocks or processes may be executed in parallel.

[0075] The conditional language used herein, such as in particular "can," "could," "might," "may," and "for example," is generally intended to convey that certain examples include certain features, elements, and / or steps, but other examples do not, unless otherwise specified or understood in the context in which they are used. Therefore, such conditional language is not generally intended to imply that features, elements, and / or steps are somehow required in one or more examples, or that one or more examples necessarily include logic for determining whether these features, elements, and / or steps are included in or performed in any particular example, with or without input or prompting from the author.

[0076] Terms such as “comprising,” “including,” and “having” are synonyms and are used comprehensively and in an open-ended manner, without excluding additional elements, features, actions, etc. Similarly, the term “or” is used in a comprehensive sense (rather than an exclusive sense), for example, when used to connect a list of elements, “or” means one, some, or all of the elements in the list. The use of “adapted to” or “configured to” herein means an open and comprehensive language that does not exclude devices adapted or configured to perform additional tasks or steps. Furthermore, the use of “based on” means open and comprehensive in that a process, step, calculation, or other action “based” on one or more enumerated conditions or values ​​may actually be based on additional conditions or values ​​beyond those enumerated. Similarly, the use of “based at least in part on” means that the process, step, calculation, or other action may actually be based on additional conditions or values ​​beyond those enumerated, “at least in part on” one or more enumerated conditions or values. The headings, lists, and numbering included herein are for illustrative purposes only and are not intended to limit the scope of the text.

[0077] The various features and processes described above may be used independently of each other or combined in various ways. All possible combinations and partial combinations are intended to fall within the scope of this disclosure. Furthermore, in some implementations, certain methods or process blocks may be omitted. The methods and processes described herein are also not limited to any particular sequence, and the blocks or states associated therewith may be executed in other appropriate sequences. For example, the described blocks or states may be executed in an order other than the specifically disclosed order, or multiple blocks or states may be combined into a single block or state. The exemplary blocks or states may be executed in series, in parallel, or in any other way. Blocks or states may be added to or removed from the disclosed examples. Similarly, the exemplary systems and components described herein may be configured differently from those described. For example, elements may be added, removed, or rearranged compared to the disclosed examples.

Claims

1. Identifying the dosage of the active ingredient in treatments provided to or being considered for provision to the target population, To predict the level of local amyloid beta based on the dose of the active ingredient, Based on the predicted level of local amyloid beta, T ２ Predicting the severity of amyloid-related imaging abnormalities (ARIA) (ARIA-E) that appear as high signals in enhanced fluid-suppressed inversion recovery (FLAIR) images, wherein predicting the severity of ARIA-E includes predicting the level of vascular wall damage, and Output the results of the prediction regarding ARIA-E severity. Computer implementation methods, including those mentioned above.

2. The method according to claim 1, wherein the active ingredient comprises an anti-amyloid monoclonal antibody.

3. The method according to claim 1, further comprising using a pharmacokinetic model to predict the time course of the active ingredient concentration in plasma during the treatment, wherein the prediction of local amyloid beta is based on at least one predicted concentration of the active ingredient over the predicted time course.

4. Predicting the aforementioned levels of local amyloid beta is possible. To estimate the baseline level of local amyloid beta, The removal rate of local amyloid beta is calculated based on the dose of the active ingredient and the baseline level of local amyloid beta, Based on the calculated removal rate of local amyloid beta, predict the level of local amyloid beta. The method according to claim 1, including the method described in claim 1.

5. The method according to claim 1, wherein predicting the level of local amyloid beta involves solving a pharmacodynamic differential equation that assumes the rate of change of amyloid beta is proportional to the product of the concentration of the active ingredient and the level of local amyloid beta in the subject.

6. The method according to claim 1, wherein predicting the level of vascular wall damage includes solving a differential equation involving the level of local amyloid beta.

7. The method according to claim 1, wherein predicting the severity of ARIA-E involves solving an algebraic equation that assumes a nonlinear relationship between the level of vascular wall damage and the Barkhof total score (BGTS).

8. The method according to claim 1, comprising identifying a schedule for monitoring ARIA events based on the predicted severity of ARIA, wherein the results characterize the schedule.

9. The method according to claim 1, comprising identifying a recommended dose of the active ingredient for the subject based on the predicted severity of ARIA, wherein the result characterizes the recommended dose.

10. One or more data processors, A non-temporary computer-readable storage medium that includes an instruction causing the one or more data processors to execute the method according to any one of claims 1 to 9 when executed by the one or more data processors, A system equipped with these features.

11. A computer program product tangibly embodied in a non-temporary machine-readable storage medium, comprising instructions configured to cause one or more data processors to perform the method according to any one of claims 1 to 9.

Images