Chemodiversity and molecular plasticity: recognition processes as explored by property spaces

Simultaneously Predicting the Pharmacokinetics of CES1-Metabolized Drugs and Their Metabolites Using Physiologically Based Pharmacokinetic Model in Cirrhosis Subjects

Hepatic carboxylesterase 1 (CES1) metabolizes numerous prodrugs into active ingredients or direct-acting drugs into inactive metabolites. We aimed to develop a semi-physiologically based pharmacokinetic (semi-PBPK) model to simultaneously predict the pharmacokinetics of CES1 substrates and their active metabolites in liver cirrhosis (LC) patients. Six prodrugs (enalapril, benazepril, cilazapril, temocapril, perindopril and oseltamivir) and three direct-acting drugs (flumazenil, pethidine and remimazolam) were selected. Parameters such as organ blood flows, plasma-binding protein concentrations, functional liver volume, hepatic enzymatic activity, glomerular filtration rate (GFR) and gastrointestinal transit rate were integrated into the simulation. The pharmacokinetic profiles of these drugs and their active metabolites were simulated for 1000 virtual individuals. The developed semi-PBPK model, after validation in healthy individuals, was extrapolated to LC patients. Most of the observations fell within the 5th and 95th percentiles of simulations from 1000 virtual patients. The estimated AUC and Cmax were within 0.5-2-fold of the observed values. The sensitivity analysis showed that the decreased plasma exposure of active metabolites due to the decreased CES1 was partly attenuated by the decreased GFR. Conclusion: The developed PBPK model successfully predicted the pharmacokinetics of CES1 substrates and their metabolites in healthy individuals and LC patients, facilitating tailored dosing of CES1 substrates in LC patients.

Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery

Simple comparative correlation analyses and quantitative structure-kinetics relationship (QSKR) models highlight the interplay of kinetic rates and binding affinity as an essential feature in drug design and discovery. The choice of the molecular series, and their structural variations, used in QSKR modeling is fundamental to understanding the mechanistic implications of ligand and/or drug-target binding and/or unbinding processes. Here, we discuss the implications of linear correlations between kinetic rates and binding affinity constants and the relevance of the computational approaches to QSKR modeling.

Approaching Pharmacological Space: Events and Components

With a view to introducing the concept of pharmacological space and its potential applications in investigating and predicting the toxic mechanisms of xenobiotics, this opening chapter describes the logical relations between conformational behavior, physicochemical properties and binding spaces, which are seen as the three key elements composing the pharmacological space. While the concept of conformational space is routinely used to encode molecular flexibility, the concepts of property spaces and, particularly, of binding spaces are more innovative. Indeed, their descriptors can find fruitful applications (a) in describing the dynamic adaptability a given ligand experiences when inserted into a specific environment, and (b) in parameterizing the flexibility a ligand retains when bound to a biological target. Overall, these descriptors can conveniently account for the often disregarded entropic factors and as such they prove successful when inserted in ligand- or structure-based predictive models. Notably, and although binding space parameters can clearly be derived from MD simulations, the chapter will illustrate how docking calculations, despite their static nature, are able to evaluate ligand's flexibility by analyzing several poses for each ligand. Such an approach, which represents the founding core of the binding space concept, can find various applications in which the related descriptors show an impressive enhancing effect on the statistical performances of the resulting predictive models.

Binding Space Concept: A New Approach To Enhance the Reliability of Docking Scores and Its Application to Predicting Butyrylcholinesterase Hydrolytic Activity

Docking simulations are very popular approaches able to assess the capacity of a given ligand to interact with a target. Docking simulations are usually focused on a single best complex even though many studies showed that ligands retain a significant mobility within a binding pocket by assuming different binding modes all of which may contribute to the monitored ligand affinity. The present study describes an innovative concept, the binding space, which allows an exploration of the ligand mobility within the binding pocket by simultaneously considering several ligand poses as generated by docking simulations. The multiple poses and the relative docking scores can then be analyzed by taking advantage of the same concepts already used in the property space analysis; hence the binding space can be parametrized by (a) mean scores, (b) score ranges, and (c) score sensitivity values. The first parameter represents a very simple procedure to account for the contribution of the often neglected alternative binding modes, while the last two descriptors encode the degree of mobility which a given ligand retains within the binding cavity (score range) as well as the ease with which a ligand explores such a mobility (score sensitivity). Here, the binding space concept is applied to the prediction of the hydrolytic activity of BChE by synergistically considering multiple poses and multiple protein structures. The obtained results shed light on the remarkable potential of the binding space concept, whose parameters allow a significant increase of the predictive power of the docking results as revealed by extended correlative analyses. Mean scores are the parameters affording the largest statistical improvement, and all the here proposed docking-based descriptors show enhancing effects in developing predictive models. Finally, the study describes a new score function (Contacts score) simply based on the number of surrounding residues which appears to be particularly productive in the framework of the binding space.

A new approach to measuring protein backbone protection with high spatial resolution using H/D exchange and electron capture dissociation

Inadequate spatial resolution remains one of the most serious limitations of hydrogen/deuterium exchange-mass spectrometry (HDX-MS), especially when applied to larger proteins (over 30 kDa). Supplementing proteolytic fragmentation of the protein in solution with ion dissociation in the gas phase has been used successfully by several groups to obtain near-residue level resolution. However, the restrictions imposed by the LC-MS/MS mode of operation on the data acquisition time frame makes it difficult in many cases to obtain a signal-to-noise ratio adequate for reliable assignment of the backbone amide protection levels at individual residues. This restriction is lifted in the present work by eliminating the LC separation step from the workflow and taking advantage of the high resolving power and dynamic range of a Fourier transform ion cyclotron resonance-mass spectrometer (FTICR-MS). A residue-level resolution is demonstrated for a peptic fragment of a 37 kDa recombinant protein (N-lobe of human serum transferrin), using electron-capture dissociation as an ion fragmentation tool. The absence of hydrogen scrambling in the gas phase prior to ion dissociation is verified using redundant HDX-MS data generated by FTICR-MS. The backbone protection pattern generated by direct HDX-MS/MS is in excellent agreement with the known crystal structure of the protein but also provides information on conformational dynamics, which is not available from the static X-ray structure.

Exploring the space of histidine containing dipeptides in search of novel efficient RCS sequestering agents

The study reports a set of forty proteinogenic histidine-containing dipeptides as potential carbonyl quenchers. The peptides were chosen to cover as exhaustively as possible the accessible chemical space, and their quenching activities toward 4-hydroxy-2-nonenal (HNE) and pyridoxal were evaluated by HPLC analyses. The peptides were capped at the C-terminus as methyl esters or amides to favor their resistance to proteolysis and diastereoisomeric pairs were considered to reveal the influence of configuration on quenching. On average, the examined dipeptides are less active than the parent compound carnosine (βAla + His) thus emphasizing the unfavorable effect of the shortening of the βAla residue as confirmed by the control dipeptide Gly-His. Nevertheless, some peptides show promising activities toward HNE combined with a remarkable selectivity. The results emphasize the beneficial role of aromatic and positively charged residues, while negatively charged and H-bonding side chains show a detrimental effect on quenching. As a trend, ester derivatives are slightly more active than amides while heterochiral peptides are more active than their homochiral diastereoisomer. Overall, the results reveal that quenching activity strongly depends on conformational effects and vicinal residues (as evidenced by the reported QSAR analysis), offering insightful clues for the design of improved carbonyl quenchers and to rationalize the specific reactivity of histidine residues within proteins.