Establishment of Correlation between in Vitro Enzyme Binding Potency and in Vivo Pharmacological Activity: Application to Liver Glycogen Phosphorylase a Inhibitors

Insights into the Competitive Mechanisms of Per- and Polyfluoroalkyl Substances Partition in Liver and Blood

Some per- and polyfluoroalkyl substances (PFASs) tend to be accumulated in liver and cause hepatotoxicity. However, the difficulty to directly measure liver concentrations of PFASs in humans hampers our understanding of their hepatotoxicity and mechanisms of action. We investigated the partitioning of 11 PFASs between liver and blood in male CD-1 mice. Although accumulation of the perfluoroalkanesulfonic acids (PFSAs) in mice serum was higher than their carboxylic acids (PFCAs) counterparts as expected, the liver-blood partition coefficients (R_L/S) of PFSAs were lower than the PFCAs R_L/S, implying a competition between liver and blood. The in vitro experiments further indicated that the partitioning was dominantly determined by their competitive binding between human liver fatty acid binding protein (hL-FABP) and serum albumin (HSA). The binding affinities (K_d) of PFASs to both proteins were measured. The correlations between the R_L/S and log K_d (hL-FABP)/log K_d (HSA) were stronger than those with log K_d (hL-FABP) alone, magnifying that the partitioning was dominantly controlled by competitive binding between hL-FABP and HSA. Therefore, the liver concentrations of the selected PFASs in humans could be predicted from the available serum concentrations, which is important for assessing their hepatotoxicity.

Does In Vitro Potency Predict Clinically Efficacious Concentrations?

The in vitro affinity of a compound for its target is an important feature in drug discovery, but what remains is how predictive in vitro properties are of in vivo therapeutic drug exposure. We assessed the relationship between in vitro potency and clinically efficacious concentrations for marketed small molecule drugs (n = 164) and how they may differ depending on therapeutic indication, mode of action, receptor type, target localization, and function. Approximately 70% of compounds had a therapeutic unbound plasma exposure lower than in vitro potency; the median ratio of exposure in relation to in vitro potency was 0.32, and 80% had ratios within the range of 0.007 to 8.7. We identified differences in the in vivo–to–in vitro potency ratio between indications, mode of action, target type, and matrix localization, and whether or not the drugs had active metabolites. The in vitro–assay variability contributions appeared to be the smallest; within the same drug target and mode of action the within‐variability was slightly broader; but both were substantially less compared with the overall distribution of ratios. These data suggest that in vitro potency conditions, estimated in vivo potency, required level of receptor occupancy, and target turnover are key components for further understanding the link between clinical drug exposure and in vitro potency.

Deprotection of N-Boc Groups under Continuous-Flow High-Temperature Conditions

The scope of thermolytic, N-Boc deprotection was studied on 26 compounds from the Pfizer compound library, representing a diverse set of structural moieties. Among these compounds, 12 substrates resulted in clean (≥95% product) deprotection, and an additional three compounds gave ≥90% product. The thermal de-Boc conditions were found to be compatible with a large number of functional groups. A combination of computational modeling, statistical analysis, and kinetic model fitting was used to support an initial, slow, and concerted proton transfer with release of isobutylene, followed by a rapid decarboxylation. A strong correlation was found to exist between the electrophilicity of the N-Boc carbonyl group and the reaction rate.

Overexpressing cell systems are a competitive option to primary adipocytes when predicting in vivo potency of dual GPR81/GPR109A agonists

Mathematical models predicting in vivo pharmacodynamic effects from in vitro data can accelerate drug discovery, and reduce costs and animal use. However, data integration and modeling is non-trivial when more than one drug-target receptor is involved in the biological response. We modeled the inhibition of non-esterified fatty acid release by dual G-protein-coupled receptor 81/109A (GPR81/GPR109A) agonists in vivo in the rat, to estimate the in vivo EC₅₀ values for 12 different compounds. We subsequently predicted those potency estimates using EC₅₀ values obtained from concentration-response data in isolated primary adipocytes and cell systems overexpressing GPR81 or GPR109A in vitro. A simple linear regression model based on data from primary adipocytes predicted the in vivo EC₅₀ better than simple linear regression models based on in vitro data from either of the cell systems. Three models combining the data from the overexpressing cell systems were also evaluated: two piecewise linear models defining logical OR- and AND-circuits, and a multivariate linear regression model. All three models performed better than the simple linear regression model based on data from primary adipocytes. The OR-model was favored since it is likely that activation of either GPR81 or GPR109A is sufficient to deactivate the cAMP pathway, and thereby inhibit non-esterified fatty acid release. The OR-model was also able to predict the in vivo selectivity between the two receptors. Finally, the OR-model was used to predict the in vivo potency of 1651 new compounds. This work suggests that data from the overexpressing cell systems are sufficient to predict in vivo potency of GPR81/GPR109A agonists, an approach contributing to faster and leaner drug discovery.

The hepatic PP1 glycogen-targeting subunit interaction with phosphorylaseacan be blocked by C-terminal tyrosine deletion or an indole drug

The inhibition of hepatic glycogen-associated protein phosphatase-1 (PP1-G(L)) by glycogen phosphorylase a prevents the dephosphorylation and activation of glycogen synthase, suppressing glycogen synthesis when glycogenolysis is activated. Here, we show that a peptide ((280)LGPYY(284)) comprising the last five amino acids of G(L) retains high-affinity interaction with phosphorylase a and that the two tyrosines play crucial roles. Tyr284 deletion abolishes binding of phosphorylase a to G(L) and replacement by phenylalanine is insufficient to restore high-affinity binding. We show that a phosphorylase inhibitor blocks the interaction of phosphorylase a with the G(L) C-terminus, suggesting that the latter interaction could be targeted to develop an anti-diabetic drug.

Astrocyte glycogen sustains neuronal activity during hypoglycemia: studies with the glycogen phosphorylase inhibitor CP-316,819 ([R-R*,S*]-5-chloro-N-[2-hydroxy-3-(methoxymethylamino)-3-oxo-1-(phenylmethyl)propyl]-1H-indole-2-carboxamide)

Glycogen in the brain is localized almost exclusively to astrocytes. The physiological function of this energy store has been difficult to establish because of the difficulty in manipulating brain glycogen concentrations in vivo. Here, we used a novel glycogen phosphorylase inhibitor, CP-316,819 ([R-R*,S*]-5-chloro-N-[2-hydroxy-3-(methoxymethylamino)-3-oxo-1-(phenylmethyl)propyl]-1H-indole-2-carboxamide), that causes glycogen accumulation under normoglycemic conditions but permits glycogen utilization when glucose concentrations are low. Rats treated with CP-316,819 had an 88 +/- 3% increase in brain glycogen content. When subjected to hypoglycemia, these rats maintained brain electrical activity 91 +/- 14 min longer than rats with normal brain glycogen levels and showed markedly reduced neuronal death. These studies establish a novel approach for manipulating brain glycogen concentration in normal, awake animals and provide in vivo confirmation that astrocyte glycogen supports neuronal function and survival during glucose deprivation. These findings also suggest an approach for forestalling hypoglycemic coma and brain injury in diabetic patients.

Three-dimensional quantitative structure-activity relationship analysis of human CYP51 inhibitors

CYP51 fulfills an essential requirement for all cells, by catalyzing three sequential mono-oxidations within the cholesterol biosynthesis cascade. Inhibition of fungal CYP51 is used as a therapy for treating fungal infections, whereas inhibition of human CYP51 has been considered as a pharmacological approach to treat dyslipidemia and some forms of cancer. To predict the interaction of inhibitors with the active site of human CYP51, a three-dimensional quantitative structure-activity relationship model was constructed. This pharmacophore model of the common structural features of CYP51 inhibitors was built using the program Catalyst from multiple inhibitors (n = 26) of recombinant human CYP51-mediated lanosterol 14alpha-demethylation. The pharmacophore, which consisted of one hydrophobe, one hydrogen bond acceptor, and two ring aromatic features, demonstrated a high correlation between observed and predicted IC(50) values (r = 0.92). Validation of this pharmacophore was performed by predicting the IC(50) of a test set of commercially available (n = 19) and CP-320626-related (n = 48) CYP51 inhibitors. Using predictions below 10 microM as a cutoff indicative of active inhibitors, 16 of 19 commercially available inhibitors (84%) and 38 of 48 CP-320626-related inhibitors (79.2%) were predicted correctly. To better understand how inhibitors fit into the enzyme, potent CYP51 inhibitors were used to build a Cerius(2) receptor surface model representing the volume of the active site. This study has demonstrated the potential for ligand-based computational pharmacophore modeling of human CYP51 and enables a high-throughput screening system for drug discovery and data base mining.