Impact of protein conformations on binding free energy calculations in the beta-secretase 1 system.

In binding free energy calculations, simulations must sample all relevant conformations of the system in order to obtain unbiased results. For instance, different ligands can bind to different metastable states of a protein, and if these protein conformational changes are not sampled in relative binding free energy calculations, the contribution of these states to binding is not accounted for and thus calculated binding free energies are inaccurate. In this work, we investigate the impact of different beta-sectretase 1 (BACE1) protein conformations obtained from x-ray crystallography on the binding of BACE1 inhibitors. We highlight how these conformational changes are not adequately sampled in typical molecular dynamics simulations. Furthermore, we show that insufficient sampling of relevant conformations induces substantial error in relative binding free energy calculations, as judged by a variation in calculated relative binding free energies up to 2 kcal/mol depending on the starting protein conformation. These results emphasize the importance of protein conformational sampling and pose this BACE1 system as a challenge case for further method development in the area of enhanced protein conformational sampling, either in combination with binding calculations or as an endpoint correction.

Discovery and Optimization of Selective Brain-Penetrant EBP Inhibitors that Enhance Oligodendrocyte Formation.

The inhibition of emopamil binding protein (EBP), a sterol isomerase within the cholesterol biosynthesis pathway, promotes oligodendrocyte formation, which has been proposed as a potential therapeutic approach for treating multiple sclerosis. Herein, we describe the discovery and optimization of brain-penetrant, orally bioavailable inhibitors of EBP. A structure-based drug design approach from literature compound 1 led to the discovery of a hydantoin-based scaffold, which provided balanced physicochemical properties and potency and an improved in vitro safety profile. The long half-lives of early hydantoin-based EBP inhibitors in rodents prompted an unconventional optimization strategy, focused on increasing metabolic turnover while maintaining potency and a brain-penetrant profile. The resulting EBP inhibitor 11 demonstrated strong in vivo target engagement in the brain, as illustrated by the accumulation of EBP substrate zymostenol after repeated dosing. Furthermore, compound 11 enhanced the formation of oligodendrocytes in human cortical organoids, providing additional support for our therapeutic hypothesis.

Kartograf: A Geometrically Accurate Atom Mapper for Hybrid-Topology Relative Free Energy Calculations.

Relative binding free energy (RBFE) calculations have emerged as a powerful tool that supports ligand optimization in drug discovery. Despite many successes, the use of RBFEs can often be limited by automation problems, in particular, the setup of such calculations. Atom mapping algorithms are an essential component in setting up automatic large-scale hybrid-topology RBFE calculation campaigns. Traditional algorithms typically employ a 2D subgraph isomorphism solver (SIS) in order to estimate the maximum common substructure. SIS-based approaches can be limited by time-intensive operations and issues with capturing geometry-linked chemical properties, potentially leading to suboptimal solutions. To overcome these limitations, we have developed Kartograf, a geometric-graph-based algorithm that uses primarily the 3D coordinates of atoms to find a mapping between two ligands. In free energy approaches, the ligand conformations are usually derived from docking or other previous modeling approaches, giving the coordinates a certain importance. By considering the spatial relationships between atoms related to the molecule coordinates, our algorithm bypasses the computationally complex subgraph matching of SIS-based approaches and reduces the problem to a much simpler bipartite graph matching problem. Moreover, Kartograf effectively circumvents typical mapping issues induced by molecule symmetry and stereoisomerism, making it a more robust approach for atom mapping from a geometric perspective. To validate our method, we calculated mappings with our novel approach using a diverse set of small molecules and used the mappings in relative hydration and binding free energy calculations. The comparison with two SIS-based algorithms showed that Kartograf offers a fast alternative approach. The code for Kartograf is freely available on GitHub (https://github.com/OpenFreeEnergy/kartograf). While developed for the OpenFE ecosystem, Kartograf can also be utilized as a standalone Python package.

The DECON pilot project investigates predictive markers for successful bariatric surgery.

Obesity is a chronic, multifactorial disease which is linked to a number of adverse endocrinological and metabolic conditions. Currently, bariatric surgery is one of the most effective treatments for individuals diagnosed with severe obesity. However, the current indications for bariatric surgery are based on inadequate metrics (i.e., BMI) which do not account for the complexity of the disease, nor the heterogeneity among the patient population. Moreover, there is a lack of understanding with respect to the biological underpinnings that influence successful and sustained weight loss post-bariatric surgery. Studies have implicated age and pre-surgery body weight as two factors that are associated with favorable patient outcomes. Still, there is an urgent medical need to identify other potential factors that could improve the specificity of candidate selection and better inform the treatment plan of patients with obesity. In this report, we present and describe the cohort of the DECON pilot project, a multicenter study which aims to identify predictive biomarkers of successful weight loss after bariatric surgery.

Broadening the Scope of Binding Free Energy Calculations Using a Separated Topologies Approach.

Binding free energy calculations predict the potency of compounds to protein binding sites in a physically rigorous manner and see broad application in prioritizing the synthesis of novel drug candidates. Relative binding free energy (RBFE) calculations have emerged as an industry-standard approach to achieve highly accurate rank-order predictions of the potency of related compounds; however, this approach requires that the ligands share a common scaffold and a common binding mode, restricting the methods' domain of applicability. This is a critical limitation since complex modifications to the ligands, especially core hopping, are very common in drug design. Absolute binding free energy (ABFE) calculations are an alternate method that can be used for ligands that are not congeneric. However, ABFE suffers from a known problem of long convergence times due to the need to sample additional degrees of freedom within each system, such as sampling rearrangements necessary to open and close the binding site. Here, we report on an alternative method for RBFE, called Separated Topologies (SepTop), which overcomes the issues in both of the aforementioned methods by enabling large scaffold changes between ligands with a convergence time comparable to traditional RBFE. Instead of only mutating atoms that vary between two ligands, this approach performs two absolute free energy calculations at the same time in opposite directions, one for each ligand. Defining the two ligands independently allows the comparison of the binding of diverse ligands without the artificial constraints of identical poses or a suitable atom-atom mapping. This approach also avoids the need to sample the unbound state of the protein, making it more efficient than absolute binding free energy calculations. Here, we introduce an implementation of SepTop. We developed a general and efficient protocol for running SepTop, and we demonstrated the method on four diverse, pharmaceutically relevant systems. We report the performance of the method, as well as our practical insights into the strengths, weaknesses, and challenges of applying this method in an industrial drug design setting. We find that the accuracy of the approach is sufficiently high to rank order ligands with an accuracy comparable to traditional RBFE calculations while maintaining the additional flexibility of SepTop.

Absolute Binding Free Energy Calculations for Buried Water Molecules.

Water often plays a key role in mediating protein-ligand interactions. Understanding contributions from active-site water molecules to binding thermodynamics of a ligand is important in predicting binding free energies for ligand optimization. In this work, we tested a non-equilibrium switching method for absolute binding free energy calculations on water molecules in binding sites of 13 systems. We discuss the lessons we learned about identified issues that affected our calculations and ways to address them. This work fits with our larger focus on how to do accurate ligand binding free energy calculations when water rearrangements are very slow, such as rearrangements due to ligand modification (as in relative free energy calculations) or ligand binding (as in absolute free energy calculations). The method studied in this work can potentially be used to account for limited water sampling via providing endpoint corrections to free energy calculations using our calculated binding free energy of water.

The Unexpected Selectivity Switching from Mitochondria to Lysosome in a D-π-A Cyanine Dye.

Two interesting benzothizolium-based D-π-A type hemicyanine dyes (3a-3b) with a diphenylamine (-NPh2) donor group were evaluated for fluorescence confocal microscopy imaging ability in live cells (MO3.13, NHLF). In sharp contrast to previously reported D-π-A dyes with alkyl amine donor (-NR2) groups (1), 3a and 3b exhibited significantly different photophysical properties and organelle selectivity. Probes 3a and 3b were nearly non-fluorescent in many polar and non-polar solvents but exhibited a bright red fluorescence (λem ≈ 630-640 nm) in stained MO3.13 and NHLF with very low probe concentrations (i.e., 200 nM). Fluorescence confocal microscopy-based co-localization studies revealed excellent lysosome selectivity from the probes 3a-3b, which is in sharp contrast to previously reported D-π-A type benzothiazolium dyes (1) with an alkyl amine donor group (-NR2) (exhibiting selectivity towards cellular mitochondria). The photostability of probe 3 was found to be dependent on the substituent (R') attached to the quaternary nitrogen atom in the cyanine dye structure. The observed donor-dependent selectivity switching phenomenon can be highly useful in designing novel organelle-targeted fluorescent probes for live-cell imaging applications.

Investigating Mechanisms of Subcutaneous Preconditioning Incubation for Neural Stem Cell Embedded Hydrogels.

Stem cells are a vital component of regenerative medicine therapies, however, only a fraction of stem cells delivered to the central nervous system following injury survive the inflammatory environment. Previously, we showed that subcutaneous preconditioning of neural stem cell (NSC) embedded hydrogels for 28 days improved spinal cord injury (SCI) functional outcomes over controls. Here, we investigated the mechanism of subcutaneous preconditioning of NSC-embedded hydrogels, with and without the known neurogenic cue, interferon gamma (IFN-γ), for 3, 14, or 28 days to refine and identify subcutaneous preconditioning conditions by measurement of neurogenic markers and cytokines. Studying the preconditioning mechanism, we found that subcutaneous foreign body response (FBR) associated cytokines infiltrated the scaffold in groups with and without NSCs, with time point effects. A pro-inflammatory environment with upregulated interleukin (IL)-6, IL-10, macrophage inflammatory protein (MIP)-1, MIP-2, IFN-γ-inducible protein 10 (IP-10), tumor necrosis factor-α (TNF-α), and IL-12p70 was observed on day 3. By 14 and 28 days, there was an increase in pro-regenerative cytokines (IL-13, IL-4) along with pro-inflammatory markers IL-1ß, IP-10, and RANTES (regulated on activation, normal T cell expressed, and secreted) potentially part of the mechanism that had an increased functional outcome in SCI. Coinciding with changes in cytokines, the macrophage population increased over time from 3 to 28 days, whereas neutrophils peaked at 3 days with a significant decrease at later time points. Expression of the neuronal marker ßIII tubulin in differentiating NSCs was supported at 3 days in the presence of soluble and immobilized IFN-γ and at 14 days by immobilized IFN-γ only, but it was greatly attenuated in all conditions at 28 days, partially because of dilution via host cell infiltration. We conclude that subcutaneously incubating NSC seeded scaffolds for 3 or 14 days could act as host specific preconditioning through exposure to FBR while retaining ßIII tubulin expression of NSCs to further improve the SCI functional outcome observed with 28 day subcutaneous incubation.

Therapeutic Targeting of Stromal-Tumor HGF-MET Signaling in an Organotypic Triple-Negative Breast Tumor Model.

The tumor microenvironment (TME) promotes proliferation, drug resistance, and invasiveness of cancer cells. Therapeutic targeting of the TME is an attractive strategy to improve outcomes for patients, particularly in aggressive cancers such as triple-negative breast cancer (TNBC) that have a rich stroma and limited targeted therapies. However, lack of preclinical human tumor models for mechanistic understanding of tumor-stromal interactions has been an impediment to identify effective treatments against the TME. To address this need, we developed a three-dimensional organotypic tumor model to study interactions of patient-derived cancer-associated fibroblasts (CAF) with TNBC cells and explore potential therapy targets. We found that CAFs predominantly secreted hepatocyte growth factor (HGF) and activated MET receptor tyrosine kinase in TNBC cells. This tumor-stromal interaction promoted invasiveness, epithelial-to-mesenchymal transition, and activities of multiple oncogenic pathways in TNBC cells. Importantly, we established that TNBC cells become resistant to monotherapy and demonstrated a design-driven approach to select drug combinations that effectively inhibit prometastatic functions of TNBC cells. Our study also showed that HGF from lung fibroblasts promotes colony formation by TNBC cells, suggesting that blocking HGF-MET signaling potentially could target both primary TNBC tumorigenesis and lung metastasis. Overall, we established the utility of our organotypic tumor model to identify and therapeutically target specific mechanisms of tumor-stromal interactions in TNBC toward the goal of developing targeted therapies against the TME. IMPLICATIONS: Leveraging a state-of-the-art organotypic tumor model, we demonstrated that CAFs-mediated HGF-MET signaling drive tumorigenic activities in TNBC and presents a therapeutic target.

Osmoregulatory Role of Betaine and Betaine/γ-Aminobutyric Acid Transporter 1 in Post-Traumatic Syringomyelia.

Syringomyelia (SM) is primarily characterized by the formation of a fluid-filled cyst that forms in the parenchyma of the spinal cord following injury or other pathology. Recent omics studies in animal models have identified dysregulation of solute carriers, channels, transporters, and small molecules associated with osmolyte regulation during syrinx formation/expansion in the spinal cord. However, their connections to syringomyelia etiology are poorly understood. In this study, the biological functions of the potent osmolyte betaine and its associated solute carrier betaine/γ-aminobutyric acid (GABA) transporter 1 (BGT1) were studied in SM. First, a rat post-traumatic SM model was used to demonstrate that the BGT1 was primarily expressed in astrocytes in the vicinity of syrinxes. In an in vitro system, we found that astrocytes uptake betaine through BGT1 to regulate cell size under hypertonic conditions. Treatment with BGT1 inhibitors, especially NNC 05-2090, demonstrated midhigh micromolar range potency in vitro that reversed the osmoprotective effects of betaine. Finally, the specificity of these BGT1 inhibitors in the CNS was demonstrated in vivo, suggesting feasibility for targeting betaine transport in SM. In summary, these data provide an enhanced understanding of the role of betaine and its associated solute carrier BGT1 in cell osmoregulation and implicates the active role of betaine and BGT1 in syringomyelia progression.

Simultaneous Visualization of Mitochondria and Lysosome by a Single Cyanine Dye: The Impact of the Donor Group (-NR2) Towards Organelle Selectivity.

A benzothiazolium-based hemicyanine dye (probe 3) has been synthesized by attaching a morpholine group into a phenyl benzothiazolium skeleton. Probe 3 exhibited interesting photophysical characteristics including red emission (λem ≈600 nm), enhanced Stokes shift (Δλ ≈80 nm) and sensitivity to solvent polarity. Although the probe 3 exhibited almost no emission in aqueous environments (φfl ≈0.002), its fluorescence could be increased by ≈50 fold in organic solvents (φfl ≈0.10), making it possible for live cell imaging under wash-free conditions. Probe 3 exhibited excellent ability to visualize cellular mitochondria and lysosomes simultaneously, as observed from fluorescence confocal microscopy. In addition, probe 3 also exhibited good biocompatibility (calculated LC50 > 20 µM) and high photostability.

Challenges Encountered Applying Equilibrium and Nonequilibrium Binding Free Energy Calculations.

Binding free energy calculations have become increasingly valuable to drive decision making in drug discovery projects. However, among other issues, inadequate sampling can reduce accuracy, limiting the value of the technique. In this paper, we apply absolute binding free energy calculations to ligands binding to T4 lysozyme L99A and HSP90 using equilibrium and nonequilibrium approaches. We highlight sampling problems encountered in these systems, such as slow side chain rearrangements and slow changes of water placement upon ligand binding. These same types of challenges are also likely to show up in other protein-ligand systems, and we propose some strategies to diagnose and test for such problems in alchemical free energy calculations. We also explore similarities and differences in how the equilibrium and the nonequilibrium approaches handle these problems. Our results show the large amount of work still to be done to make free energy calculations robust and reliable and provide insight for future research in this area.

Softening of the chronic hemi-section spinal cord injury scar parallels dysregulation of cellular and extracellular matrix content.

Regeneration following spinal cord injury (SCI) is challenging in part due to the modified tissue composition and organization of the resulting glial and fibrotic scar regions. Inhibitory cell types and biochemical cues present in the scar have received attention as therapeutic targets to promote regeneration. However, altered Young's modulus of the scar as a readout for potential impeding factors for regeneration are not as well-defined, especially in vivo. Although the decreased Young's modulus of surrounding tissue at acute stages post-injury is known, the causation and outcomes at chronic time points remain largely understudied and controversial, which motivates this work. This study assessed the glial and fibrotic scar tissue's Young's modulus and composition (scar morphometry, cell identity, extracellular matrix (ECM) makeup) that contribute to the tissue's stiffness. The spatial Young's modulus of a chronic (~18-wks, post-injury) hemi-section, including the glial and fibrotic regions, were significantly less than naïve tissue (~200 Pa; p < 0.0001). The chronic scar contained cystic cavities dispersed in areas of dense nuclei packing. Abundant CNS cell types such as astrocytes, oligodendrocytes, and neurons were dysregulated in the scar, while epithelial markers such as vimentin were upregulated. The key ECM components in the CNS, namely sulfated proteoglycans (sPGs), were significantly downregulated following injury with concomitant upregulation of unsulfated glycosaminoglycans (GAGs) and hyaluronic acid (HA), likely altering the foundational ECM network that contributes to tissue stiffness. Our results reveal the Young's modulus of the chronic SCI scar as well as quantification of contributing elastic components that can provide a foundation for future study into their role in tissue repair and regeneration.

Metabolomic and Signaling Programs Induced by Immobilized versus Soluble IFN γ in Neural Stem Cells.

Neural stem cells (NSCs) provide a strategy to replace damaged neurons following traumatic central nervous system injuries. A major hurdle to translation of this therapy is that direct application of NSCs to CNS injury does not support sufficient neurogenesis due to lack of proper cues. To provide prolonged spatial cues to NSCs IFN-γ was immobilized to biomimetic hydrogel substrate to supply physical and biochemical signals to instruct the encapsulated NSCs to be neurogenic. However, the immobilization of factors, including IFN-γ, versus soluble delivery of the same factor, has been incompletely characterized especially with respect to activation of signaling and metabolism in cells over longer time points. In this study, protein and metabolite changes in NSCs induced by immobilized versus soluble IFN-γ at 7 days were evaluated. Soluble IFN-γ, refreshed daily over 7 days, elicited stronger responses in NSCs compared to immobilized IFN-γ, indicating that immobilization may not sustain signaling or has altered ligand/receptor interaction and integrity. However, both IFN-γ delivery types supported increased ßIII tubulin expression in parallel with canonical and noncanonical receptor-signaling compared to no IFN-γ. Global metabolomics and pathway analysis revealed that soluble and immobilized IFN-γ altered metabolic pathway activities including energy, lipid, and amino acid synthesis, with soluble IFN-γ having the greatest metabolic impact overall. Finally, soluble and immobilized IFN-γ support mitochondrial voltage-dependent anion channel (VDAC) expression that correlates to differentiated NSCs. This work utilizes new methods to evaluate cell responses to protein delivery and provides insight into mode of action that can be harnessed to improve regenerative medicine-based strategies.

Large-Scale Assessment of Binding Free Energy Calculations in Active Drug Discovery Projects.

Accurate ranking of compounds with regards to their binding affinity to a protein using computational methods is of great interest to pharmaceutical research. Physics-based free energy calculations are regarded as the most rigorous way to estimate binding affinity. In recent years, many retrospective studies carried out both in academia and industry have demonstrated its potential. Here, we present the results of large-scale prospective application of the FEP+ method in active drug discovery projects in an industry setting at Merck KGaA, Darmstadt, Germany. We compare these prospective data to results obtained on a new diverse, public benchmark of eight pharmaceutically relevant targets. Our results offer insights into the challenges faced when using free energy calculations in real-life drug discovery projects and identify limitations that could be tackled by future method development. The new public data set we provide to the community can support further method development and comparative benchmarking of free energy calculations.

Efficient synthesis of NIR emitting bis[2-(2'-hydroxylphenyl)benzoxazole] derivative and its potential for imaging applications.

Unassymetric bis[2-(2'-hydroxyphenylbenzoxole)] bis(HBO) derivatives with a DPA functionality for zinc binding have been developed with an efficient synthetic route, using the retrosynthetic analysis. Comparison of bis(HBO) derivatives with different substitution patterns allows us to verify and optimize their unique fluorescence properties. Upon binding zinc cation, bis(HBO) derivatives give a large fluorescence turn-on in both visible (λem ≈ 536 nm) and near-infrared (NIR) window (λem ≈ 746 nm). The probes are readily excitable by a 488 nm laser, making this series of compounds a suitable imaging tool for in vitro and in vivo study on a confocal microscope. The application of zinc binding-induced fluorescence turn-on is successfully demonstrated in cellular environments and thrombus imaging.

Synthesis of highly selective lysosomal markers by coupling 2-(2'-hydroxyphenyl)benzothiazole (HBT) with benzothiazolium cyanine (Cy): the impact of substituents on selectivity and optical properties.

HBT-Cy 1 has been previously reported as a highly selective fluorescent probe for lysosome visualization in live cells. To further investigate the role of the structural components of HBT-Cy in lysosome selectivity, cyanine based fluorescent probe series (2-5) have been synthesized in good yields by connecting benzothiazolium cyanine (Cy) with 2-hydroxyphenylbenzothiazole (HBT) via a meta phenylene ring. Probes 2-5 exhibited exceptional photophysical properties including bright red-emission (λem≈ 630-650 nm), a large Stokes shift (Δλ > 130 nm) and high fluorescence quantum yields (φfl≈ 0.1-0.5). Probes 2, 3, and 5 exhibited exceptional selectivity towards cellular lysosomes in NHLF and MO3.13 cells. Our further study revealed that the phenyl benzothiazolium cyanine component (6) was the lysosome directing group in the HBT-Cy probe structure. The attachment of the hydroxyphenyl benzothiazole (HBT) component to the HBT-Cy probe structure has significantly improved its photophysical properties. Lysosome probes 2, 3 and 5 exhibited excellent biocompatibility, quick staining, bright red fluorescence, and wash-free application for live cell imaging. These probes further exhibited excellent characteristics for bioimaging experiments including a non-alkalinizing nature, high biocompatibility, high photostability and long-term imaging ability (>4 hours).

NIR-emitting benzothiazolium cyanines with an enhanced stokes shift for mitochondria imaging in live cells.

A series of benzothiazolium-based hemicyanines (3a-3f) have been synthesized. Evaluation of their photophysical properties shows that they exhibit improved photophysical characteristics. In comparison with the available commercial MitoTrackers, the new probes revealed an enhanced Stokes shift (Δλ ∼ 80 nm) and minimized aggregation for increased sensitivity. The synthesized probes are found to exhibit excellent selectivity for mitochondrial staining in an oligodendrocyte cell line. Probes show almost no fluorescence in aqueous environments, while the fluorescence is increased by ∼10-fold in organic solvents, making it possible for mitochondrial imaging without the need for post-staining washing. Since the absorption peaks of probes are close to the laser wavelengths of 561 and 640 nm on a commercial confocal microscope, e.g.3a exhibits λabs ∼ 620 nm and λem ∼ 702 nm, they could be useful probes for mitochondrial tracking in live cells.

A fluorescent flavonoid for lysosome detection in live cells under "wash free" conditions.

Lysosomes are vital organelles in living cells, which have acidic environments (pH 4.0-5.0) where macrobiomolecules and malfunctioning organelles are broken down into monomers by hydrolase activity. The majority of the currently reported fluorescent probes for detecting lysosomes suffer from small Stokes shifts (Δλ < 20 nm) and higher cytotoxicity due to an "alkalinizing effect". An interesting flavonoid-based lysosome probe is synthesized by introducing a morpholine moiety onto the flavonoid skeleton. This new probe has shown excellent selectivity to detect lysosomes in MO3.13 oligodendrocytes and normal human lung fibroblast cell lines. Probes 1a and 1b have shown excellent fluorescence quantum yield (φfl up to 0.43 in non-aqueous solvents) and large Stokes shifts (120-150 nm). These new fluorescent probes also exhibit a large quantum yield difference from an aqueous to organic environment, making them potentially useful as "wash-free" stains for visualizing lysosomes. Cell viability evaluation of these probes shows excellent biocompatibility with the median lethal concentration being LC50 ≈ 50 µM.

Time-dependence in mixture toxicity prediction.

The value of time-dependent toxicity (TDT) data in predicting mixture toxicity was examined. Single chemical (A and B) and mixture (A+B) toxicity tests using Microtox(®) were conducted with inhibition of bioluminescence (Vibrio fischeri) being quantified after 15, 30 and 45-min of exposure. Single chemical and mixture tests for 25 sham (A1:A2) and 125 true (A:B) combinations had a minimum of seven duplicated concentrations with a duplicated control treatment for each test. Concentration/response (x/y) data were fitted to sigmoid curves using the five-parameter logistic minus one parameter (5PL-1P) function, from which slope, EC25, EC50, EC75, asymmetry, maximum effect, and r(2) values were obtained for each chemical and mixture at each exposure duration. Toxicity data were used to calculate percentage-based TDT values for each individual chemical and mixture of each combination. Predicted TDT values for each mixture were calculated by averaging the TDT values of the individual components and regressed against the observed TDT values obtained in testing, resulting in strong correlations for both sham (r(2)=0.989, n=25) and true mixtures (r(2)=0.944, n=125). Additionally, regression analyses confirmed that observed mixture TDT values calculated for the 50% effect level were somewhat better correlated with predicted mixture TDT values than at the 25 and 75% effect levels. Single chemical and mixture TDT values were classified into five levels in order to discern trends. The results suggested that the ability to predict mixture TDT by averaging the TDT of the single agents was modestly reduced when one agent of the combination had a positive TDT value and the other had a minimal or negative TDT value.