Analysis of Quantitative High Throughput Screening Data

数千种广泛商业使用的化学品尚未测试对人类的不利影响，但存在于环境中。因此，需要改进用于体内毒性测试的化学品优先级，并且最终需要找到用于评估大量潜在有害物质的基于细胞的替代品。定量高通量筛选（qHTS）试验是多浓度实验，在国家毒理学计划的努力中发挥着重要作用，以应对这些测试挑战，并将毒理学从主要的观察性科学推进到主要的预测性科学。qHTS可以在一个广阔的化学空间内同时分析数千种化学物质，降低每种物质的成本。 先前用于从qHTS数据进行活动调用的方法是基于寻求最小化假阳性的制药应用，并且通常依赖于化学分析而不是统计测试来进行活动调用。例如，我们开发了一种三阶段算法，将来自qHTS数据的物质分类为与毒理学评价相关的统计学支持的活性类别，寻求提高灵敏度，同时最大限度地降低I类错误率（Shockley，2012）。我们的方法的第一阶段拟合四参数Hill方程，以找到在测试浓度范围内具有稳健浓度响应曲线的活性物质。第二阶段寻找在最低测试浓度下具有实质性活性的相对强效物质，即第一阶段未捕获的物质。该算法的第三阶段也是最后一阶段将统计上显著的特征与缺乏统计上令人信服的支持或不活跃的响应分开。该框架可容纳大量的qHTS数据，容忍丢失的数据，并且不需要重复测量。 上述三阶段算法基于希尔方程模型。然而，浓度-响应数据可能很复杂，并且在数据中找到不基于S形曲线拟合的替代模式可能更具信息性。由于非线性回归模型拟合qHTS实验中生成的数据得出的参数估计值存在很大的不确定性，因此通常不清楚如何为后续研究确定化学品的优先顺序。因此，我们开发了一种加权熵分数（WES）作为平均活性水平的度量，以便在qHTS实验中对化学品进行排名（Shockley，2014）。WES评分可用于在没有预先指定的模型结构的情况下对测试库中的所有化学品进行排名，或者WES可用于通过对返回的“命中”进行排名来补充现有方法。WES的性能进行了评估，使用希尔模型模拟的数据。在qHTS研究的典型条件下，WES优于基于AC 50（半数最大反应的估计浓度）的排名。

Thousands of chemicals in wide commercial use have not been tested for adverse effects on humans, but are present in the environment. Accordingly, there is a need to improve chemical prioritization for in vivo toxicity testing and, ultimately, to find cell-based alternatives for evaluating the large inventory of potentially harmful substances. Quantitative high throughput screening (qHTS) assays are multiple-concentration experiments with an important role in the efforts of the National Toxicology Program to meet these testing challenges and advance toxicology from a predominantly observational science to a predominantly predictive science. qHTS can simultaneously assay thousands of chemicals over a wide chemical space with reduced cost per substance. Previous approaches for making activity calls from qHTS data were based on pharmaceutical applications seeking to minimize false positives and usually relied on heuristics rather than statistical tests to make activity calls. For instance, we developed a three-stage algorithm to classify substances from qHTS data into statistically supported activity categories relevant to toxicological evaluation, seeking to improve sensitivity while minimizing Type I error rate (Shockley, 2012). The first stage of our approach fits a four-parameter Hill equation to find active substances with a robust concentration-response profile within the tested concentration range. The second stage finds relatively potent substances with substantial activity at the lowest tested concentration, substances not captured in the first stage. The third and final stage of the algorithm separates statistically significant profiles from responses that lack statistically compelling support, or inactives. This framework accommodates large volumes of qHTS data, tolerates missing data, and does not require replicate measurements. The three-stage algorithm described above is based on the Hill equation model. However, concentration-response data can be complex and it may be more informative to find alternative patterns in the data not based on fits to sigmoidal curves. It is often unclear how to prioritize chemicals for follow-up studies due to the large uncertainties that accompany parameter estimates derived from nonlinear regression model fits to data generated in qHTS experiments. Therefore, we developed a weighted entropy score (WES) as a measure of average activity level in order to rank chemical in qHTS experiments (Shockley, 2014). WES scores can be used to rank all chemicals in a tested library without a pre-specified model structure, or WES can be used to complement existing approaches by ranking returned "hits". The performance of WES has been evaluated using data simulated from a Hill model. WES outperforms rankings based on AC50 (estimated concentration of half-maximal response) across the full range of conditions that are typical of qHTS studies.