On the scaling interpretation of exponents in hyperboloid models of delay and probability discounting

Previously, we (McKerchar et al., 2009) showed that two-parameter hyperboloid models ( [Green and Myerson, 2004] and [Rachlin, 2006]) provide significantly better fits to delay discounting data than simple, one-parameter hyperbolic and exponential models. Here, we extend this effort by comparing fits of the two-parameter hyperboloid models to data from a larger sample of participants (N = 171) who discounted probabilistic as well as delayed rewards. In particular, we examined the effects of amount on the exponents in the two hyperboloid models of delay and probability discounting in order to evaluate key theoretical predictions of the standard psychophysical scaling interpretation of these exponents. Both the Rachlin model and the Green and Myerson model provided very good fits to delay and probability discounting of both small and large amounts at both the group and individual levels (all R²s > .97 at the group level; all median R²s > .92 at the individual level). For delay discounting, the exponent in both models did not vary as a function of delayed amount, consistent with the psychophysical scaling interpretation. For probability discounting, however, the exponent in both models increased as the probabilistic amount increased—a finding inconsistent with the scaling interpretation.     

Tolerance of apexes of coral Pocillopora damicornis L. to cryoprotectant solutions

In this study, we investigated the tolerance of Pocillopora damicornis apexes to treatments with solutions containing penetrating and non-penetrating cryoprotective agents (CPAs). CPAs were employed individually or in binary, tertiary or quaternary solutions. In some experiments apexes were treated successively with two CPA solutions with increasing total concentration. P. damicornis apexes withstood exposure for up to 30 min to solutions containing 0.6–0.8 M sucrose (Suc) or trehalose (Tre). When apexes were treated with binary cryoprotectant solutions containing Suc and ethylene glycol (EG), methanol (Meth), dimethyl sulfoxide (Me₂SO) or glycerol (Gly), the CPAs employed in combination with Suc could be ranked in the following order of decreasing tolerance: EG > Meth > Me₂SO > Gly. P. damicornis apexes tolerated exposure to complex CPA solutions containing Suc, Me₂SO, EG and/or Meth with a total molarity of 2.45 M. In experiments where two successive CPA solutions were employed, apexes withstood treatment with the second, more concentrated solution at 0 °C for up to 10 min. These preliminary results pave the way to the development of a cryopreservation protocol for P. damicornis apexes.     

Systematic multiscale parameterization of heterogeneous elastic network models of proteins

We present a method to parameterize heterogeneous elastic network models (heteroENMs) of proteins to reproduce the fluctuations observed in atomistic simulations. Because it is based on atomistic simulation, our method allows the development of elastic coarse-grained models of proteins under different conditions or in different environments. The method is simple and applicable to models at any level of coarse-graining. We validated the method in three systems. First, we computed the persistence length of ADP-bound F-actin, using a heteroENM model. The value of 6.1 ± 1.6 μm is consistent with the experimentally measured value of 9.0 ± 0.5 μm. We then compared our method to a uniform elastic network model and a realistic extension algorithm via covariance Hessian (REACH) model of carboxy myoglobin, and found that the heteroENM method more accurately predicted mean-square fluctuations of α-carbon atoms. Finally, we showed that the method captures critical differences in effective harmonic interactions for coarse-grained models of the N-terminal Bin/amphiphysin/Rvs (N-BAR) domain of amphiphysin, by building models of N-BAR both bound to a membrane and free in solution.     

Modeling local and global intracellular calcium responses mediated by diffusely distributed inositol 1,4,5-trisphosphate receptors

Considerable insight into intracellular Ca²⁺ responses has been obtained through the development of whole cell models that are based on molecular mechanisms, e.g., single channel kinetics of the inositol 1,4,5-trisphosphate (IP₃) receptor Ca²⁺ channel. However, a limitation of most whole cell models to date is the assumption that IP₃ receptor Ca²⁺ channels (IP₃Rs) are globally coupled by a “continuously stirred” bulk cytosolic [Ca²⁺], when in fact open IP₃Rs experience elevated “domain” Ca²⁺ concentrations. Here we present a 2N+2-compartment whole cell model of local and global Ca²⁺ responses mediated by N=100,000 diffusely distributed IP₃Rs, each represented by a four-state Markov chain. Two of these compartments correspond to bulk cytosolic and luminal Ca²⁺ concentrations, and the remaining 2N compartments represent time-dependent cytosolic and luminal Ca²⁺ domains associated with each IP₃R. Using this Monte Carlo model as a starting point, we present an alternative formulation that solves a system of advection-reaction equations for the probability density of cytosolic and luminal domain [Ca²⁺] jointly distributed with IP₃R state. When these equations are coupled to ordinary differential equations for the bulk cytosolic and luminal [Ca²⁺], a realistic but minimal model of whole cell Ca²⁺ dynamics is produced that accounts for the influence of local Ca²⁺ signaling on channel gating and global Ca²⁺ responses. The probability density approach is benchmarked and validated by comparison to Monte Carlo simulations, and the two methods are shown to agree when the number of Ca²⁺ channels is large (i.e., physiologically realistic). Using the probability density approach, we show that the time scale of Ca²⁺ domain formation and collapse (both cytosolic and luminal) may influence global Ca²⁺ oscillations, and we derive two reduced models of global Ca²⁺ dynamics that account for the influence of local Ca²⁺ signaling on global Ca²⁺ dynamics when there is a separation of time scales between the stochastic gating of IP₃Rs and the dynamics of domain Ca²⁺.     

Differential in vivo response of soft-shell clam hemocytes against two strains of Vibrio splendidus: Changes in cell structure, numbers and adherence

Host-pathogen interaction models in aquatic species are useful tools for understanding the pathogenicity of diseases in cultured and wild populations. In this study we report the differential in vivo response of soft-shell clam (Mya arenaria) hemocytes against two strains of Vibrio splendidus. Responses were measured 24 h after injecting into the posterior adductor muscle either an endemic wild-type strain (7SHRW) or a strain associated with oyster mortalities (LGP32-GFP). Changes in hemocyte structure (percentage of rounded cells) were assessed microscopically. Changes in adherence and hemocyte numbers were analyzed by flow-cytometric cell counting. Increased percentages of rounded cells were found in response to both strains. However, values from the group infected with LGP32-GFP were significantly higher (p < 0.01) than with 7SHRW. The cell adherence was markedly diminished (p < 0.001) by LGP32-GFP whereas 7SHRW did not change it significantly. Increased numbers of hemocytes (p < 0.001) were induced by LGP32-GFP, while no significant changes were found after infection with 7SHRW. These results show the regulatory capacity of soft-shell clams hemocytes to perform specific responses against different strains of V. splendidus.     

Application of screen-printed microband biosensors to end-point measurements of glucose and cell numbers in HepG2 cell culture

Microband glucose biosensors were produced by insulating and sectioning through a screen-printed, water-based carbon electrode containing cobalt phthalocyanine redox mediator and glucose oxidase enzyme. Under quiescent conditions at 37 °C, at an operating potential of +0.4 V, they produced an amperometric response to glucose in buffer solutions with a sensitivity of 26.4 nA/mM and a linear range of 0.45 to 9.0 mM. An optimal pH value of 8.5 was obtained under these conditions, and a value for activation energy of 40.55 kJ mol⁻¹ was calculated. In culture medium (pH 7.3), a sensitivity of 13 nA/mM was obtained and the response was linear up to 5 mM with a detection limit of 0.5 mM. The working concentration was up to 20 mM glucose with a precision of 11.3% for replicate biosensors (n = 4). The microband biosensors were applied to determine end-point glucose concentrations in culture medium by monitoring steady-state current responses 400 s after transfer of the biosensors into different sample solutions. In conjunction with cultures of HepG2 (human Caucasian hepatocyte carcinoma) cells, current responses obtained in 24-h supernatants showed an inverse correlation (R² = 0.98) with cell number, indicating that the biosensors were applicable for monitoring glucose metabolism by cells and of quantifying cell number. Glucose concentrations determined using the biosensor assay were in good agreement, for concentrations up to 20 mM, with those determined spectrophotometrically (R² = 0.99). This method of end-point glucose determination was used to provide an estimated rate of glucose uptake for HepG2 cells of 7.9 nmol/(10⁶ cells min) based on a 24-h period in culture.     

Omiganan interaction with bacterial membranes and cell wall models. Assigning a biological role to saturation

Omiganan pentahydrochloride (ILRWPWWPWRRK-NH₂·5Cl) is an antimicrobial peptide currently in phase III clinical trials. This study aims to unravel the mechanism of action of this drug at the membrane level and address the eventual protective role of peptidoglycan in cell walls. The interaction of omiganan pentahydrochloride with bacterial and mammalian membrane models - large unilamellar vesicles of different POPC:POPG proportions - was characterized by UV-Vis fluorescence spectroscopy. The molar ratio partition constants obtained for the two anionic bacterial membrane models were very high ((18.9 ± 1.3) × 10³ and (43.5 ± 8.7) × 10³) and about one order of magnitude greater than for the neutral mammalian models ((3.7 ± 0.4) × 10³ for 100% POPC bilayers). At low lipid:peptide ratios there were significant deviations from the usual hyperbolic-like partition behavior of peptide vesicle titration curves, especially for the most anionic systems. Membrane saturation can account for such observations and mathematical models were derived to further characterize the peptide-lipid interaction under those conditions; a possible relation between saturation and MIC was deduced; this was supported by differential quenching studies of peptide internalization. Interaction with the bacterial cell wall was assessed using Staphylococcus aureus peptidoglycan extracts as a model. A strong partition towards the peptidoglycan mesh was observed, but not as large as for the membrane models.     

The effects of the Brazilian antDinoponera quadriceps venom on chemically induced seizure models

Arthropod venoms are potential sources of neuroactive substances, providing new tools for the design of drugs. The aim of this study was to evaluate the effects of Dinoponera quadriceps venom (DqV) on seizure models in mice induced by pentylenetetrazole (PTZ), pilocarpine, and strychnine. In the PTZ model, intraperitoneal treatment with DqV (0.5 mg/kg) increased the time until the first seizure and the percentage of survival (155.4 ± 27.7 s/12.5%, p < 0.05) compared to the control group (79.75 ± 3.97 s/0%), whereas endovenous treatment (0.1 and 0.5 mg/kg) decreased the time until the first seizure (0.1 mg/kg: 77.83 ± 5.3 s versus 101.0 ± 3.3 s in the control group; 0.5 mg/kg: 74.43 ± 3.9 s versus 101.0 ± 3.3 s for the control group, p < 0.05). We did not observe significant changes in the pilocarpine- and strychnine-induced seizure models. In assays that measured oxidative parameters in the PTZ model, intraperitoneal treatment with DqV (0.5 and 2.0 mg/kg) only decreased the levels of MDA and nitrite in the cortex. However, endovenous treatment with DqV (0.1 and 0.5 mg/kg) increased the levels of MDA in the cortex and hippocampus and at a dose of 0.5 mg/kg in the striatum. Moreover, increased in nitrite content was observed in all three of the brain regions analyzed. Taken together, the D. quadriceps venom caused both neuroprotective and neurotoxic effects in a PTZ-induced seizure model, and this effect was dependent on the route of administration used.     

The viscoelasticity of membrane tethers and its importance for cell adhesion

Cell adhesion mechanically couples cells to surfaces. The durability of individual bonds between the adhesive receptors and their ligands in the presence of forces determines the cellular adhesion strength. For adhesive receptors such as integrins, it is a common paradigm that the cell regulates its adhesion strength by altering the affinity state of the receptors. However, the probability distribution of rupture forces is dependent not only on the affinity of individual receptor-ligand bonds but also on the mechanical compliance of the cellular anchorage of the receptor. Hence, by altering the anchorage, the cell can regulate its adhesion strength without changing the affinity of the receptor. Here, we analyze the anchorage of the integrin VLA-4 with its ligand VCAM-1. For this purpose, we develop a model based on the Kelvin body, which allows one to quantify the mechanical properties of the adhesive receptor's anchorage using atomic force microscopy on living cells. As we demonstrate, the measured force curves give valuable insight into the mechanics of the cellular anchorage of the receptor, which is described by the tether stiffness, the membrane rigidity, and the membrane viscosity. The measurements relate to a tether stiffness of k_t = 1.6 μN/m, an initial membrane rigidity of k_i = 260 μN/m, and a viscosity of μ = 5.9 μN·s/m. Integrins exist in different activation states. When activating the integrin with Mg²⁺, we observe altered viscoelastic parameters of k_t = 0.9 μN/m, k_i = 190 μN/m, and μ = 6.0 μ N·s/m. Based on our model, we postulate that anchorage-related effects are common regulating mechanisms for cellular adhesion beyond affinity regulation.     

The effects of temperature and diet on age grading and population age structure determination in Drosophila

The age structure of natural population is of interest in physiological, life history and ecological studies but it is often difficult to determine. One methodological problem is that samples may need to be invasively sampled preventing subsequent taxonomic curation. A second problem is that it can be very expensive to accurately determine the age structure of given population because large sample sizes are often necessary. In this study, we test the effects of temperature (17 °C, 23 °C and 26 °C) and diet (standard cornmeal and low calorie diet) on the accuracy of the non-invasive, inexpensive and high throughput near-infrared spectroscopy (NIRS) technique to determine the age of Drosophila flies. Composite and simplified calibration models were developed for each sex. Independent sets for each temperature and diet treatments with flies not involved in calibration model were then used to validate the accuracy of the calibration models. The composite NIRS calibration model was generated by including flies reared under all temperatures and diets. This approach permits rapid age measurement and age structure determination in large population of flies as less than or equal to 9 days, or more than 9 days old with 85–97% and 64–99% accuracy, respectively. The simplified calibration models were generated by including flies reared at 23 °C on standard diet. Low accuracy rates were observed when simplified calibration models were used to identify (a) Drosophila reared at 17 °C and 26 °C and (b) 23 °C with low calorie diet. These results strongly suggest that appropriate calibration models need to be developed in the laboratory before this technique can be reliably used in field. These calibration models should include the major environmental variables that change across space and time in the particular natural population to be studied.     

Qualitative analysis of models with different treatment protocols to prevent antibiotic resistance

This paper is concerned with the qualitative analysis of two models [S. Bonhoeffer, M. Lipsitch, B.R. Levin, Evaluating treatment protocols to prevent antibiotic resistance, Proc. Natl. Acad. Sci. USA 94 (1997) 12106] for different treatment protocols to prevent antibiotic resistance. Detailed qualitative analysis about the local or global stability of the equilibria of both models is carried out in term of the basic reproduction number R₀. For the model with a single antibiotic therapy, we show that if R₀ < 1, then the disease-free equilibrium is globally asymptotically stable; if R₀ > 1, then the disease-endemic equilibrium is globally asymptotically stable. For the model with multiple antibiotic therapies, stabilities of various equilibria are analyzed and combining treatment is shown better than cycling treatment. Numerical simulations are performed to show that the dynamical properties depend intimately upon the parameters.     

Genetic and non-genetic factors associated with milking order in lactating dairy cows

The ability to rapidly identify temporal deviations of an animal from its norm will be important in the management of individual cows in large herds. Furthermore, predictors of genetic merit for especially health traits are useful to augment the accuracy of selection, and thus genetic gain, in breeding programs. The objective of this study was to estimate the repeatability of milking order and to quantify the contribution of differences in additive genetic variation to phenotypic differences (i.e., heritability). The data used in this study included 9813 herd milk recording test-day records with time of milking from 85,532 cows in 1143 herds across an 8-year period. Milking order was available for both morning and evening milking for each cow with, on average, 3.33 milk test-day records (i.e., 6.66 milking events) per lactation, and on average 1.62 lactations per cow. Variance components for milking order were estimated using animal linear mixed models; covariance components between milking order and milk yield, milk composition and somatic cell score (i.e., logarithm₁₀ somatic cell count) were estimated also using animal linear mixed models. The heritability of milking order was 0.20 indicating partial genetic control of milking order. The repeatability of milking order within test-day, within lactation, and across lactations was 0.63, 0.51, and 0.47, respectively. Milking order was positively (P < 0.001), but weakly, phenotypically correlated with milk yield (r = 0.04), and milk fat concentration (r = 0.01) and negatively (P < 0.001), but weakly, correlated with milk protein concentration (r = −0.02) and somatic cell score (r = −0.05). Milking order was positively (P < 0.05), although weakly, genetically correlated with milk yield (r = 0.07) and negatively (P < 0.05), but also weakly, genetically correlated with somatic cell score (r = −0.08). This study is the first to show a contribution of additive genetics to milking order in dairy cattle but the genetic correlation between milking order and somatic cell score was weak.     

The RTK/ERK pathway is associated with prostate cancer risk on the SNP level: A pooled analysis of 41 sets of data from case–control studies

Prostate cancer (PCa) is a malignant disease influencing numerous men worldwide every year. However, the exact pathogenesis and the genes, environment, and other factors involved have not been explained clearly. Some studies have proposed that cell signaling pathways might play a key role in the development and progression of PCa. According to our previous study, the RTK/ERK pathway containing nearly 40 genes was associated with PCa risk. On the basis of these genes, we conducted a meta-analysis with our own Chinese Consortium for Prostate Cancer Genetics (ChinaPCa) study and available studies in the databases to describe the association between the pathway and PCa on the SNP level. The results suggested that rs4764695/IGF1 (recessive model: pooled OR = 0.92, 95%CI = 0.852–0.994, P = 0.034; I² = 0%, P = 0.042; allele analysis: pooled OR = 0.915, 95%CI = 0.874–0.958, P = 0; I² = 0%, P = 0.424; codominant model: OR = 0.835, 95%CI = 0.762–0.916, P = 0; I² = 0%, P = 0.684) and rs1570360/VEGF (recessive model: OR = 0.596, 95%CI = 0.421–0.843, P = 0.003; I² = 23.9%, P = 0.269; codominant model: OR = 0.576, 95%CI = 0.404–0.820, P = 0.002; I² = 49.1%, P = 0.140) were significantly associated with PCa. In subgroup analysis, the relationship was also found in Caucasians for IGF1 (dominant model: OR = 0.834, 95%CI = 0.769–0.904, P = 0; allele analysis: OR = 0.908, 95%CI = 0.863–0.955, P = 0; AA vs CC: OR = 0.829, 95%CI = 0.750–0.916, P = 0; AC vs CC: OR = 0.837, 95%CI = 0.768–0.912, P = 0). In addition, in Asians (allele analysis: OR = 0.21, 95%CI = 0.168–0.262, P = 0) and Caucasians (recessive model: OR = 0.453, 95%CI: 0.240–0.855, P = 0.015; codominant model: OR = 0.464, 95%CI = 0.240–0.898, P = 0.023) for VEGF, the association was significant. The results indicated that rs4764695/IGF1 and rs1570360/VEGF might play a key role in the development and progression of PCa. On the SNP level, we suggest that the study gives us a new view of gene-pathway analysis and targeted therapy for PCa.     

Bio-kinetic analysis on treatment of textile dye wastewater using anaerobic batch reactor

An anaerobic digestion technique was applied to textile dye wastewater aiming at the colour and COD removal. Pet bottles of 5 L capacity were used as reactor which contains methanogenic sludge of half a liter capacity which was used for the treatment of combined synthetic textile dye and starch wastewater at different mixing ratios of 20:80, 30:70, 40:60, 50:50 and 60:40 with initial COD concentrations as 3520, 3440, 3360, 3264 and 3144 mg L⁻¹, respectively. The reactor was maintained at room temperature (30 ± 3 °C) with initial pH of 7. The maximum COD and colour removal were 81.0% and 87.3% at an optimum mixing ratio of 30:70 of textile dye and starch wastewaters. Both Monod’s and Haldane’s models were adopted in this study. The kinetic constants of cell growth under Haldane’s model were satisfactory when compared to Monod’s model. The kinetic constants obtained by Haldane’s model were found to be in the range of μ_max = 0.037-0.146 h⁻¹, K_s = 651.04-1372.88 mg L⁻¹ and K_i = 5681.81-18727.59 mg L⁻¹.     

n − 3, but not n − 6 lipid particle uptake requires cell surface anchoring

Omega-3 (n − 3) fatty acids are emerging as bioactive agents protective against cardiovascular disease. However, their cellular delivery pathways are poorly defined. Here we questioned whether the uptake of n − 3 triglyceride-rich particles (TGRP) is mediated by cell surface proteoglycans (PG) using LDL receptor (LDLR)+/+ and LDLR−/− cell models. LDLR+/+ but not LDLR−/− cells showed higher n − 6 over n − 3 TGRP uptake. Removal of cell surface proteins and receptors by pronase markedly enhanced the uptake of n − 3 but not n − 6 TGRP. Lactoferrin blockage of apoE-mediated pathways decreased the uptake of n − 6 TGRP by up to 85% (p < 0.05) but had insignificant effect on n − 3 TGRP uptake. PG removal by sodium chlorate in LDLR+/+ cells substantially reduced n − 3 TGRP uptake but had little effect on n − 6 TGRP uptake. Thus, while n − 6 TGRP uptake is preferentially mediated by LDLR-dependent pathways, the uptake of n − 3 TGRP depends more on PG and non-LDLR cell surface anchoring.     

Determination of the membrane permeability characteristics of Pacific oyster, Crassostrea gigas, oocytes and development of optimized methods to add and remove ethylene glycol

In order for cryopreservation to become a practical tool for aquaculture, optimized protocols must be developed for each species and cell type. Knowledge of a cell’s osmotic tolerance and membrane permeability characteristics can assist in optimized protocol development. In this study, these characteristics were determined for Pacific oyster oocytes and modified methods for loading and unloading ethylene glycol (EG) were tested. Oocytes were found to behave as ideal osmometers and their osmotically inactive fraction (V_b) was calculated to be 0.48. Oocytes exposed to NaCl solutions of 0.6 to 2.3 Osm fertilized at rates equivalent to oocytes left in seawater. This corresponds to volume changes of +27.3 and −38.1 ± 1.2%. The permeability of the oocytes to water (L_p) was determined to be 3.8 ± 0.4 × 10⁻², 5.7 ± 0.8 × 10⁻², and 13.2 ± 1.3 × 10⁻² μm min⁻¹ atm⁻¹, when measured at temperatures of 5, 10 and 20 °C. The respective EG permeability values (P_s) were 9.5 ± 0.1 × 10⁻⁵, 14.6 ± 1.2 × 10⁻⁵, and 41.7 ± 2.4 × 10⁻⁵ cm min⁻¹. The activation energies for L_p and P_s were determined to be 14.5 and 17.5 kcal mol⁻¹, respectively. Different models for EG loading and unloading from oocytes were developed and tested. Post-thaw fertilization did not differ significantly between a published step addition method and single step addition at 20 °C. This represents a considerable reduction in handling. The results of this study demonstrate that the cryobiological characteristics of a given cell type should be taken into account when developing cryopreservation methods.     

Multi-objective optimization of glycopeptide antibiotic production in batch and fed batch processes

Fermentation optimization involves potentially conflicting multiple objectives such as product concentration and production media cost. Simultaneous optimization of these objectives would result in a multiobjective optimization problem, which is characterized by a set of multiple solutions, knows as pareto optimal solutions. These solutions gives flexibility in evaluating the trade-offs and selecting the most suitable operating policy. Here, ε-constraint approach was used to generate the pareto solutions for two objectives: product concentration and product per unit cost of media, for batch and fed batch operations using process model for Amycolatopsis balhimycina, a glycopeptide antibiotic producer. This resulted in a set of several pareto optimal solutions with the two objectives ranging from (0.75 g l⁻¹, 3.97 g $^-1) to (0.44 g l⁻¹, 5.19 g $^-1) for batch and from (1.5 g l⁻¹, 5.46 g $^-1) to (1.1 g l⁻¹, 6.34 g $^-1) for fed batch operations. One pareto solution each for batch and for fed batch mode was experimentally validated.     

A novel in vivo regulatory role of P-glycoprotein in alloimmunity

P-glycoprotein (P-gp) is required for adaptive immunity through defined functions in T cell activation and antigen presenting cell (APC) maturation. The potential role of P-gp as an in vivo regulator of alloimmunity is currently unknown. Here we show that P-gp blockade prolongs graft survival in a murine heterotopic cardiac allotransplantation model through in vivo inhibition of the T helper 1 (Th1) cytokine IFN-γ and the Th2 product IL-4, and via downregulation of the APC-expressed positive costimulatory molecule CD80. In vitro, the P-gp antagonist PSC833, a non-calcineurin-inhibitory cyclosporine A analogue, specifically inhibited cellular efflux of the P-gp substrate rhodamine-123 in wild-type CD3⁺ T cells and MHC class II⁺ APCs but not their P-gp knockout counterparts that lacked rhodamine-123 efflux capacity. Additionally, P-gp blockade significantly inhibited murine alloimmune T cell activation in a dose-dependent fashion. In vivo, P-gp blockade significantly prolonged graft survival in Balb/c recipients of C57BL/6 cardiac allografts from 8.5 ± 0.5 to 11.7 ± 0.5 days (P < 0.01), similar in magnitude to the effects of monotherapy with cyclosporine A. Moreover, P-gp blockade, compared to controls, attenuated intragraft expression of CD3 and CD80, but not CD86, and inhibited IFN-γ and IL-4 production (P < 0.05). In the setting of systemic CD86 inhibition, P-gp blockade suppressed IFN-γ and IL-4 production significantly further (to 98% and 89% inhibition, respectively) compared to either P-gp or anti-CD86 blockade alone, and markedly prolonged allograft survival compared to anti-CD86 blockade alone (40.5 ± 4.6 versus 22.5 ± 2.6 days, respectively, P < 0.01). Our findings define a novel in vivo regulatory role of P-gp in alloimmunity and identify P-gp as a potential therapeutic target in allotransplantation.