The goal of this method is to characterize the drug sensitivity of multiple myeloma (MM) primary cells to a panel of agents ex vivo, as close as possible to physiological conditions, and with sufficient precision that these results can be used to identify chemoresistant sub-populations within the tumor burden and, ultimately, to parameterize computational models designed to estimate clinical response.
Computational models are powerful tools to analyze complex systems, such as cancer-host-therapy interactions. However, models are only as good as the data used to parameterize them. Unfortunately, most data available in literature cannot be used in its current form to parameterize such models, as they often are obtained in mutually exclusive experimental conditions. Thus, new experiments are often required with the goal of obtaining the set of experimental parameters needed. Most viability assays, however, are destructive and thus limited to a small number of time points, often only one. This greatly handicaps the usage of chemosensitivity assays to parameterize computational models, since such experiments fail to provide information on the temporal dynamics of the system. This is especially true with primary cancer cells, which are often limited to a few millions per sample, and have a short lifespan after biopsy, thus limiting the number of experimental conditions available. In addition, most viability assays require the separation of the cancer from the non-cancer (stroma) cells at the time of quantification, which adds extra work to the protocol, further perturbs the cells, and limits the number of experiments that can be performed.
40 years ago, Salmon and collaborators1 proposed their famous in vitro colony formation assay for assessment of chemosensitivity of tumor cells. Unfortunately this assay found mixed success due to the small number of multiple myeloma (MM) patient samples that were capable of forming colonies under control conditions: It was estimated that only a small subgroup of 0.001% to 0.1% of MM cells were capable of replication, being termed as multiple myeloma stem cells2. This significantly limited the success rate of the assay and the number of drugs that could be tested at one time, even in more recent models3. This issue is much more important now when MM agents (standard or experimental) are more numerous than four decades ago.
A major limitation of these early assays is their dichotomized output: either a patient is “sensitive” or “resistant” to a particular drug. No information is provided regarding the degree of sensitivity and heterogeneity of the tumor population. Thus, patients with small sub-populations of fast-growing resistant cells would be classified as “sensitive” to therapy, but would relapse shortly, and thus not benefit from treatment. Given the importance of duration of response in overall survival (OS) of cancer patients4,5, it is clear how these assays were not able to properly estimate OS consistently.
In order to circumvent these limitations, and to be able to generate patient-specific computational models that would estimate personalized clinical response to a panel of drugs, we have developed a method for non-destructively testing drug sensitivity of multiple myeloma (MM) cell lines and primary MM cells in an ex vivo reconstruction of the bone marrow microenvironment, including extracellular matrix and stroma6. This assay, however, had the limitation of relying on a particular commercial microfluidic slide, whose dimensions and cost restricted the number of experiments or drugs that could be performed at once in a given piece of equipment.
The here described system extends this original assay into a high-throughput organotypic dose-response platform, for in vitro screening of drugs, based on a digital image analysis algorithm to non-destructively quantify cell viability. Each well in 384 or 1,536-well plates is a 3D reconstruction of the bone marrow microenvironment, including primary MM cells, extracellular matrix, and patient-derived stroma and growth factors. Live microscopy and digital image analysis are used to detect cell death events in different drug concentrations, which are used to generate dose-response surfaces. From the in vitro data, a mathematical model identifies the size and chemosensitivity of sub-populations within the patient’s tumor burden, and can be used to simulate how the tumor would respond to the drug(s) in physiological conditions in a clinical regimen6 (Figure 1).
The main innovations of this platform are: (a) small number of cancer cells required (1,000-10,000 per drug concentration); (b) assessment of drug efficacy in physiological conditions (extracellular matrix, stroma, patient-derived growth factors); (c) No toxicity from viability markers7 since only bright field imaging is used, thus no need to transfect cells with fluorescence8 or bioluminescence9; (d) continuous imaging provides drug effect as a function of concentration and exposure time (pharmacodynamics); and (e) the integration between in vitro and computational evolutionary models, to estimate clinical outcome6 (Silva et al., in preparation).
The key factor that allows the digital image analysis algorithm to discern cancer from stromal cells is their different affinity to the bottom of the well: MM cells do not adhere, and remain in suspension in the matrix, while stromal cells adhere to the bottom of the well and then stretch.
This separation occurs even though MM and stroma are seeded at the same time, and occurs during the O/N process of incubation. The spinning down in the centrifuge following seeding of the plate further accelerates this process. As a consequence, MM cells appear in the image as round bright disks surrounded by a dark ring, while the stroma maintains a low profile and a darker shade (Figure 2). Thus, the assay in its current form is best suited for non-adherent cancer cells, although adjustments in the algorithm could be done to separate cancer and stroma based on other morphologic features, such as size and shape. In this protocol we describe two types of extracellular matrix (collagen and basement membrane matrix) as well as two types of co-cultures, with (bone marrow derived mesenchymal stem cells) BMSC and HUVECS, representing the interstitial and perivascular niches of the bone marrow, respectively.
Using a 384-well plate, 5 concentrations per drug and two replicates, we were able to test up to 31 different drugs with the sample of a single patient. In 1,536-well plates this number is 127. Figure 3 depicts the layout currently used for manual cell seeding, while Supplemental Figure 1 represents the optimum distribution using a robotic pipettor. In the design from Figure 3, each 384-well plate can carry 3 patient samples with 7 different drugs, or one patient sample tested with 21 drugs. Each drug is represented by 5 concentrations, and each condition is seeded in duplicates. 4 control wells (no drug added) are seeded for each set of 7 drugs (e.g., wells 36-39, 126-129 and 200-203). To ensure the potency of the drugs used, a human myeloma cell line (e.g., H929 or MM1.S) is also seeded, and drugged in duplicates at the highest concentration of each of the drugs used (wells 76-89). The positive cell line control ensures that the drugs used have adequate potency, since their dose response curves are comparable across experiments. Two wells are seeded with a myeloma cell line as control for the environmental conditions, in other words, to detect possible problems in the bench top incubator during imaging (wells 75 and 90). The enumeration of the wells follows a zigzag pattern in order to reduce the distance covered by the motorized stage of the microscope, which reduces acquisition time and reduces wear of the equipment.