Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations.

We here report on nonequilibrium targeted molecular dynamics simulations as a tool for the estimation of protein-ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90, we show that the mean nonequilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein-ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and/or van der Waals interactions with the target. In the case of ligands with a rigid chemical scaffold that exhibit longer residence times, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.

The use of virtual screening and differential scanning fluorimetry for the rapid identification of fragments active against MEK1.

We report the analysis of an in-house fragment screening campaign for the oncology target MEK1. The application of virtual screening (VS) as a primary fragment screening approach, followed by biophysical validation using differential screening fluorimetry (DSF), with resultant binding mode determination by X-ray crystallography (X-ray), is presented as the most time and cost-effective combination of in silico and in vitro methods to identify fragments. We demonstrate the effectiveness of the VS-DSF workflow for the early identification of fragments to both 'jump-start' the drug discovery project and to complement biochemical screening data.

Predictive models assessing the response to ustekinumab highlight the value of therapeutic drug monitoring in Crohn's disease.

BACKGROUND: Despite the therapeutic efficacy of Ustekinumab (UST) in Crohn's disease (CD), loss of response (LOR) is observed over time. This study aims to evaluate the impact of the UST pharmacokinetics (PK) at induction on clinical and endoscopic outcomes, as well as to find predictive markers of UST response. METHODS: This retrospective study included 80 CD patients. Pharmacokinetics data (trough levels (TLs)) combined with clinical and biological parameters were fed into tailored logistic regression and tree-based ensemble techniques to predict clinical and endoscopic outcomes at one year of follow-up. RESULTS: TLs at week 16 were significantly lower among patients with moderate to severe endoscopic activity during the follow-up (p = 0.04). The best model to predict endoscopic outcome was obtained at week 16 by Random Forest with an area under the receiver operating characteristic curve of 0.92 ± 0.08, sensitivity 91% and specificity 75%, with key inputs such as lymphocyte and monocyte counts at week 8, and UST TLs and CRP at week 16. CONCLUSIONS: This real-world study confirms the relationship between early UST TLs and both clinical and endoscopic outcomes. Models were developed for the task of predicting clinical and endoscopic remission in CD patients treated with UST, highlighting the clinical relevance of UST TLs at week 16.

Collecting New Peak and Intermediate Infliximab Levels to Predict Remission in Inflammatory Bowel Diseases.

BACKGROUND: The loss of response to infliximab is a challenge for clinicians in the management of inflammatory bowel disease (IBD). Mounting evidence suggests that therapeutic drug monitoring at induction may predict remission during maintenance. The aim of the study was to improve predictive models of remission by exploring new peak and intermediate infliximab measurements during induction. METHODS: This was a prospective multicenter study evaluating the pharmacokinetics of infliximab during induction in a pioneer cohort of 63 patients with IBD. Pharmacokinetics data including peak, intermediate, and trough levels were combined with clinical and biological parameters and were subsequently fed into tailored logistic regression and tree-based techniques to predict remission at week 30. RESULTS: Infliximab peak levels at week 2, intermediate levels at week 3, and trough levels at week 6 were correlated with remission at week 30. Predictive models exhibited an increased accuracy over the successive timepoints of the induction with key inputs such as albumin, C-reactive protein, eosinophils, neutrophils, lymphocytes, intermediate level at week 3, trough level at week 6, and age at diagnosis. Our predictive model of remission at week 30 was obtained with an area under the receiver operating characteristic curve of 0.9 ±â€…0.12, a sensitivity of 89%, and a specificity of 75%. CONCLUSIONS: This study showed the clinical relevance of measuring new infliximab levels to predict remission in patients with IBD. These findings lay the foundation for a personalized medicine in which biotherapies could be monitored at an early stage, thereby improving patients' clinical management.

New approach to determine the healthy immune variations by combining clustering methods.

Immune-mediated inflammatory diseases are characterized by variability in disease presentation and severity but studying it is a challenging task. Defining the limits of a healthy immune system is therefore a prior step to capture variability in disease conditions. The goal of this study is to characterize the global immune cell composition along with their influencing factors. Blood samples were collected from 2 independent cohorts of respectively 389 (exploratory) and 208 (replication) healthy subjects. Twelve immune cells were measured in blood together with biological parameters. Three complementary clustering approaches were used to evaluate if variability related to the immune cells could be characterized as clusters or as a continuum. Large coefficients of variation confirmed the inter-individual variability of immune cells. Considering all subset variations in an overall analysis, it appeared that the immune makeup was organized as a continuum through the two cohorts. Some intrinsic and environmental factors affected the inter-individual variability of cells but without unveiling separable groups with similar features. This study provides a framework based on complementary clustering approach for analyzing inter-individual variability of immune cells. Our analyses support the absence of clusters in our two healthy cohorts. Also, our study reports some influence of age, gender, BMI, cortisol, season and CMV infection on immune variability.

Discovery of 6-(2,4-Dichlorophenyl)-5-[4-[(3S)-1-(3-fluoropropyl)pyrrolidin-3-yl]oxyphenyl]-8,9-dihydro-7H-benzo[7]annulene-2-carboxylic acid (SAR439859), a Potent and Selective Estrogen Receptor Degrader (SERD) for the Treatment of Estrogen-Receptor-Positive Breast Cancer.

More than 75% of breast cancers are estrogen receptor alpha (ERα) positive (ER+), and resistance to current hormone therapies occurs in one-third of ER+ patients. Tumor resistance is still ERα-dependent, but mutations usually confer constitutive activation to the hormone receptor, rendering ERα modulator drugs such as tamoxifen and aromatase inhibitors ineffective. Fulvestrant is a potent selective estrogen receptor degrader (SERD), which degrades the ERα receptor in drug-resistant tumors and has been approved for the treatment of hormone-receptor-positive metastatic breast cancer following antiestrogen therapy. However, fulvestrant shows poor pharmacokinetic properties in human, low solubility, weak permeation, and high metabolism, limiting its administration to inconvenient intramuscular injections. This Drug Annotation describes the identification and optimization of a new series of potent orally available SERDs, which led to the discovery of 6-(2,4-dichlorophenyl)-5-[4-[(3S)-1-(3-fluoropropyl)pyrrolidin-3-yl]oxyphenyl]-8,9-dihydro-7H-benzo[7]annulene-2-carboxylic acid (43d), showing promising antitumor activity in breast cancer mice xenograft models and whose properties warranted clinical evaluation.

Fast NMR Methods for Measuring in the Direct and/or Competition Mode the Dissociation Constants of Chemical Fragments Interacting with a Receptor.

Ligand-based NMR screening represents a powerful method in fragment-based drug discovery for the identification of chemical matter interacting with the receptor of interest. The large dynamic range of these methods allows the detection of weakly binding ligands. However, the methodology has not been extensively used for quantifying the strength of these interactions. This knowledge is important for ranking fragments according to their binding strength and for prioritizing structure-based and medicinal chemistry activities. Rapid NMR methods for measuring the dissociation constant in direct and competition modes are presented here. The theory underpinning these methods are presented, along with their application to the measurement of the binding affinities of several ligands of the heat shock protein 90.

Estimation of Drug-Target Residence Times by τ-Random Acceleration Molecular Dynamics Simulations.

Drug-target residence time (τ), one of the main determinants of drug efficacy, remains highly challenging to predict computationally and, therefore, is usually not considered in the early stages of drug design. Here, we present an efficient computational method, τ-random acceleration molecular dynamics (τRAMD), for the ranking of drug candidates by their residence time and obtaining insights into ligand-target dissociation mechanisms. We assessed τRAMD on a data set of 70 diverse drug-like ligands of the N-terminal domain of HSP90α, a pharmaceutically important target with a highly flexible binding site, obtaining computed relative residence times with an accuracy of about 2.3τ for 78% of the compounds and less than 2.0τ within congeneric series. Analysis of dissociation trajectories reveals features that affect ligand unbinding rates, including transient polar interactions and steric hindrance. These results suggest that τRAMD will be widely applicable as a computationally efficient aid to improving drug residence times during lead optimization.

An evaporation-free solid-phase extraction method for rapid and accurate analysis of sumatriptan in human plasma by LC-MS/MS.

An evaporation-free solid-phase extraction (SPE) method was developed and validated for sumatriptan. High organic washing (50% methanol) and low organic elution (20% methanol) were used and the recovery was greater than 92%. The eluate was injected into a C18 column without evaporation and reconstitution. Sumatriptan was monitored in positive ion mode with mass transition of m/z 296.4-58.1 amu. The calibration curve was 1-100 ng/mL (r>or=0.9923). The inter-day and intra-day precisions ranged from 4.53 to 9.12% and 1.72 to 6.93%, respectively. This method features reduced cost and pollution, clean extract, high speed, and most importantly overall method reliability.

Reconstructing Demography and Social Behavior During the Neolithic Expansion from Genomic Diversity Across Island Southeast Asia.

Archaeology, linguistics, and increasingly genetics are clarifying how populations moved from mainland Asia, through Island Southeast Asia, and out into the Pacific during the farming revolution. Yet key features of this process remain poorly understood, particularly how social behaviors intersected with demographic drivers to create the patterns of genomic diversity observed across Island Southeast Asia today. Such questions are ripe for computer modeling. Here, we construct an agent-based model to simulate human mobility across Island Southeast Asia from the Neolithic period to the present, with a special focus on interactions between individuals with Asian, Papuan, and mixed Asian-Papuan ancestry. Incorporating key features of the region, including its complex geography (islands and sea), demographic drivers (fecundity and migration), and social behaviors (marriage preferences), the model simultaneously tracks a full suite of genomic markers (autosomes, X chromosome, mitochondrial DNA, and Y chromosome). Using Bayesian inference, model parameters were determined that produce simulations that closely resemble the admixture profiles of 2299 individuals from 84 populations across Island Southeast Asia. The results highlight that greater propensity to migrate and elevated birth rates are related drivers behind the expansion of individuals with Asian ancestry relative to individuals with Papuan ancestry, that offspring preferentially resulted from marriages between Asian women and Papuan men, and that in contrast to current thinking, individuals with Asian ancestry were likely distributed across large parts of western Island Southeast Asia before the Neolithic expansion.

An Automated Microscale Thermophoresis Screening Approach for Fragment-Based Lead Discovery.

Fragment-based lead discovery has proved to be an effective alternative to high-throughput screenings in identifying chemical matter that can be developed into robust lead compounds. The search for optimal combinations of biophysical techniques that can correctly and efficiently identify and quantify binding can be challenging due to the physicochemical properties of fragments. In order to minimize the time and costs of screening, optimal combinations of biophysical techniques with maximal information content, sensitivity, and robustness are needed. Here we describe an approach utilizing automated microscale thermophoresis (MST) affinity screening to identify fragments active against MEK1 kinase. MST identified multiple hits that were confirmed by X-ray crystallography but not detected by orthogonal methods. Furthermore, MST also provided information about ligand-induced aggregation and protein denaturation. The technique delivered a large number of binders while reducing experimentation time and sample consumption, demonstrating the potential of MST to execute and maximize the efficacy of fragment screening campaigns.

CODV-Ig, a universal bispecific tetravalent and multifunctional immunoglobulin format for medical applications.

Bispecific immunoglobulins (Igs) typically contain at least two distinct variable domains (Fv) that bind to two different target proteins. They are conceived to facilitate clinical development of biotherapeutic agents for diseases where improved clinical outcome is obtained or expected by combination therapy compared to treatment by single agents. Almost all existing formats are linear in their concept and differ widely in drug-like and manufacture-related properties. To overcome their major limitations, we designed cross-over dual variable Ig-like proteins (CODV-Ig). Their design is akin to the design of circularly closed repeat architectures. Indeed, initial results showed that the traditional approach of utilizing (G4S)x linkers for biotherapeutics design does not identify functional CODV-Igs. Therefore, we applied an unprecedented molecular modeling strategy for linker design that consistently results in CODV-Igs with excellent biochemical and biophysical properties. CODV architecture results in a circular self-contained structure functioning as a self-supporting truss that maintains the parental antibody affinities for both antigens without positional effects. The format is universally suitable for therapeutic applications targeting both circulating and membrane-localized proteins. Due to the full functionality of the Fc domains, serum half-life extension as well as antibody- or complement-dependent cytotoxicity may support biological efficiency of CODV-Igs. We show that judicious choice in combination of epitopes and paratope orientations of bispecific biotherapeutics is anticipated to be critical for clinical outcome. Uniting the major advantages of alternative bispecific biotherapeutics, CODV-Igs are applicable in a wide range of disease areas for fast-track multi-parametric drug optimization.

Discovery and Pharmacokinetic and Pharmacological Properties of the Potent and Selective MET Kinase Inhibitor 1-{6-[6-(4-Fluorophenyl)-[1,2,4]triazolo[4,3-b]pyridazin-3-ylsulfanyl]benzothiazol-2-yl}-3-(2-morpholin-4-ylethyl)urea (SAR125844).

The HGF/MET pathway is frequently activated in a variety of cancer types. Several selective small molecule inhibitors of the MET kinase are currently in clinical evaluation, in particular for NSCLC, liver, and gastric cancer patients. We report herein the discovery of a series of triazolopyridazines that are selective inhibitors of wild-type (WT) MET kinase and several clinically relevant mutants. We provide insight into their mode of binding and report unprecedented crystal structures of the Y1230H variant. A multiparametric chemical optimization approach allowed the identification of compound 12 (SAR125844) as a development candidate. In this chemical series, absence of CYP3A4 inhibition was obtained at the expense of satisfactory oral absorption. Compound 12, a promising parenteral agent for the treatment of MET-dependent cancers, promoted sustained target engagement at tolerated doses in a human xenograft tumor model. Preclinical pharmacokinetics conducted in several species were predictive for the observed pharmacokinetic behavior of 12 in cancer patients.

Selective blockade of the hydrolysis of the endocannabinoid 2-arachidonoylglycerol impairs learning and memory performance while producing antinociceptive activity in rodents.

Monoacylglycerol lipase (MAGL) represents a primary degradation enzyme of the endogenous cannabinoid (eCB), 2-arachidonoyglycerol (2-AG). This study reports a potent covalent MAGL inhibitor, SAR127303. The compound behaves as a selective and competitive inhibitor of mouse and human MAGL, which potently elevates hippocampal levels of 2-AG in mice. In vivo, SAR127303 produces antinociceptive effects in assays of inflammatory and visceral pain. In addition, the drug alters learning performance in several assays related to episodic, working and spatial memory. Moreover, long term potentiation (LTP) of CA1 synaptic transmission and acetylcholine release in the hippocampus, two hallmarks of memory function, are both decreased by SAR127303. Although inactive in acute seizure tests, repeated administration of SAR127303 delays the acquisition and decreases kindled seizures in mice, indicating that the drug slows down epileptogenesis, a finding deserving further investigation to evaluate the potential of MAGL inhibitors as antiepileptics. However, the observation that 2-AG hydrolysis blockade alters learning and memory performance, suggests that such drugs may have limited value as therapeutic agents.

Data mining crystallization databases: knowledge-based approaches to optimize protein crystal screens.

Protein crystallization is a major bottleneck in protein X-ray crystallography, the workhorse of most structural proteomics projects. Because the principles that govern protein crystallization are too poorly understood to allow them to be used in a strongly predictive sense, the most common crystallization strategy entails screening a wide variety of solution conditions to identify the small subset that will support crystal nucleation and growth. We tested the hypothesis that more efficient crystallization strategies could be formulated by extracting useful patterns and correlations from the large data sets of crystallization trials created in structural proteomics projects. A database of crystallization conditions was constructed for 755 different proteins purified and crystallized under uniform conditions. Forty-five percent of the proteins formed crystals. Data mining identified the conditions that crystallize the most proteins, revealed that many conditions are highly correlated in their behavior, and showed that the crystallization success rate is markedly dependent on the organism from which proteins derive. Of the proteins that crystallized in a 48-condition experiment, 60% could be crystallized in as few as 6 conditions and 94% in 24 conditions. Consideration of the full range of information coming from crystal screening trials allows one to design screens that are maximally productive while consuming minimal resources, and also suggests further useful conditions for extending existing screens.

A pilot study to evaluate the pharmacokinetics of sibutramine in healthy subjects under fasting and fed conditions.

PURPOSE: The purpose of this study was to characterize the food effect on the pharmacokinetics of sibutramine and its pharmacologically active metabolites. METHODS: This was an open label, single dose, crossover study completed by six healthy males. A single dose of sibutramine 15 mg was administered orally under fasting and fed conditions. Plasma concentrations of sibutramine and its metabolites were determined by LC/MS/MS: Non-compartmental pharmacokinetics and statistical analysis were performed using SAS. RESULTS: The food intake increased significantly AUCs and C(max) of sibutramine and its M1 metabolite, but did not affect M2 metabolite. When sibutramine was administered with food, the T(max) was delayed by 2 to 4 hours for sibutramine, M1 and M2 metabolites as stated in the literature. CONCLUSIONS: The results of this study indicate that sibutramine is measurable using a sensitive bioanalytical method. In contrast with what is reported in the product monograph, this study demonstrated that the bioavailability of sibutramine and M1 metabolite was significantly increased with administration with food. The results confirmed lack of food effect on the pharmacokinetics of M2 metabolite. These relatively large food effect observed for sibutramine and M1 metabolite, could have implication for the efficacy and safety of the drug.

SAR650984, a novel humanized CD38-targeting antibody, demonstrates potent antitumor activity in models of multiple myeloma and other CD38+ hematologic malignancies.

PURPOSE: The CD38 cell surface antigen is expressed in diverse hematologic malignancies including multiple myeloma, B-cell non-Hodgkin lymphoma (NHL), B-cell chronic lymphocytic leukemia, B-cell acute lymphoblastic leukemia (ALL), and T-cell ALL. Here, we assessed the antitumor activity of the anti-CD38 antibody SAR650984. EXPERIMENTAL DESIGN: Activity of SAR650984 was examined on lymphoma, leukemia and multiple myeloma cell lines, primary multiple myeloma samples, and multiple myeloma xenograft models in immunodeficient mice. RESULTS: We identified a humanized anti-CD38 antibody with strong proapoptotic activity independent of cross-linking agents, and potent effector functions including complement-dependent cytotoxicity, antibody-dependent cell-mediated cytotoxicity, and antibody-dependent cellular phagocytosis (ADCP), equivalent in vitro to rituximab in CD20+ and CD38+ models. This unique antibody, termed SAR650984, inhibited the ADP-ribosyl cyclase activity of CD38, likely through an allosteric antagonism as suggested by 3D structure analysis of the complex. In vivo, SAR650984 was active in diverse NHL, ALL, and multiple myeloma CD38+ tumor xenograft models. SAR650984 demonstrated single-agent activity comparable with rituximab or cyclophosphamide in Daudi or SU-DHL-8 lymphoma xenograft models with induction of the proapoptotic marker cleaved capase-7. In addition, SAR650984 had more potent antitumor activity than bortezomib in NCI-H929 and Molp-8 multiple myeloma xenograft studies. Consistent with its mode of action, SAR650984 demonstrated potent proapoptotic activity against CD38+ human primary multiple myeloma cells. CONCLUSION: These results validate CD38 as a therapeutic target and support the current evaluation of this unique CD38-targeting functional antibody in phase I clinical trials in patients with CD38+ B-cell malignancies.

Comparison between elementary flux modes analysis and 13C-metabolic fluxes measured in bacterial and plant cells.

BACKGROUND: (13)C metabolic flux analysis is one of the pertinent ways to compare two or more physiological states. From a more theoretical standpoint, the structural properties of metabolic networks can be analysed to explore feasible metabolic behaviours and to define the boundaries of steady state flux distributions. Elementary flux mode analysis is one of the most efficient methods for performing this analysis. In this context, recent approaches have tended to compare experimental flux measurements with topological network analysis. RESULTS: Metabolic networks describing the main pathways of central carbon metabolism were set up for a bacteria species (Corynebacterium glutamicum) and a plant species (Brassica napus) for which experimental flux maps were available. The structural properties of each network were then studied using the concept of elementary flux modes. To do this, coefficients of flux efficiency were calculated for each reaction within the networks by using selected sets of elementary flux modes. Then the relative differences - reflecting the change of substrate i.e. a sugar source for C. glutamicum and a nitrogen source for B. napus - of both flux efficiency and flux measured experimentally were compared. For both organisms, there is a clear relationship between these parameters, thus indicating that the network structure described by the elementary flux modes had captured a significant part of the metabolic activity in both biological systems. In B. napus, the extension of the elementary flux mode analysis to an enlarged metabolic network still resulted in a clear relationship between the change in the coefficients and that of the measured fluxes. Nevertheless, the limitations of the method to fit some particular fluxes are discussed. CONCLUSION: This consistency between EFM analysis and experimental flux measurements, validated on two metabolic systems allows us to conclude that elementary flux mode analysis could be a useful tool to complement (13)C metabolic flux analysis, by allowing the prediction of changes in internal fluxes before carbon labelling experiments.

Tricyclic series of heat shock protein 90 (Hsp90) inhibitors part I: discovery of tricyclic imidazo[4,5-c]pyridines as potent inhibitors of the Hsp90 molecular chaperone.

A novel class of heat shock protein 90 (Hsp90) inhibitors was developed after a low throughput screen (LTS) of a focused library containing approximately 21K compounds selected by virtual screening. The initial [1-{3-H-imidazo[4-5-c]pyridin-2-yl}-3,4-dihydro-2H-pyrido[2,1-a]isoindole-6-one] (1) compound showed moderate activity (IC(50) = 7.6 µM on Hsp82, the yeast homologue of Hsp90). A high-resolution X-ray structure shows that compound 1 binds into an "induced" hydrophobic pocket, 10-15 Å away from the ATP/resorcinol binding site. Iterative cycles of structure-based drug design (SBDD) and chemical synthesis led to the design and preparation of analogues with improved affinity. These optimized molecules make productive interactions within the ATP binding site as reported by other Hsp90 inhibitors. This resulted in compound 8, which is a highly potent inhibitor in biochemical and cellular assays (K(d) = 0.35 nM on Hsp90; IC(50) = 30 nM on SKBr3 mammary carcinoma cells) and in an in vivo leukemia model.

In vitro disintegration and dissolution and in vivo bioequivalence of two alendronate once-weekly formulations.

Bioequivalence of two tablet formulations of 70 mg alendronate (CAS 121268-17-5) was assessed in a single-dose, open-label, randomised, fasted state crossover trial, with a washout period of 21 days, in 80 healthy subjects. Urine samples were collected up to +36 h post dosing and the concentrations of alendronic acid were assessed using a high-performance liquid chromatographic method with pre-derivatization and fluorescence detection (HPLC/FL) method. The 90 % confidence intervals (90 % CI) obtained for Ae0-36 (cumulative urinary excretion) and Rmax (maximum rate of urinary excretion) were 98.67-118.99 % and 102.22-122.46 %, respectively. The intra-subject coefficient of variation was between 32-35 % for both parameters. No relevant tolerability problems were detected. Both formulations can be considered bioequivalent. In vitro testing was performed to confirm the adequacy of the quality control conditions and no significant differences were detected neither in the disintegration test nor in the dissolution tests conducted in HCl 0.1 N and H2O and thus in these conditions the lack of statistically significant differences in vitro was accompanied by in vivo bioequivalence.