Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (.100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

Docking to heme proteins.

In silico screening has become a valuable tool in drug design, but some drug targets represent real challenges for docking algorithms. This is especially true for metalloproteins, whose interactions with ligands are difficult to parametrize. Our docking algorithm, EADock, is based on the CHARMM force field, which assures a physically sound scoring function and a good transferability to a wide range of systems, but also exhibits difficulties in case of some metalloproteins. Here, we consider the therapeutically important case of heme proteins featuring an iron core at the active site. Using a standard docking protocol, where the iron-ligand interaction is underestimated, we obtained a success rate of 28% for a test set of 50 heme-containing complexes with iron-ligand contact. By introducing Morse-like metal binding potentials (MMBP), which are fitted to reproduce density functional theory calculations, we are able to increase the success rate to 62%. The remaining failures are mainly due to specific ligand-water interactions in the X-ray structures. Testing of the MMBP on a second data set of non iron binders (14 cases) demonstrates that they do not introduce a spurious bias towards metal binding, which suggests that they may reliably be used also for cross-docking studies.

EV02: a Phase I trial to compare the safety and immunogenicity of HIV DNA-C prime-NYVAC-C boost to NYVAC-C alone.

The aim of this randomised controlled trial was to see if the addition of 4 mg/ml DNA-C priming given by the intramuscular route at weeks 0 and 4 to NYVAC-C at weeks 20 and 24, safely increased the proportion of participants with HIV-specific T-cell responses measured by the interferon (IFN)-gamma ELISpot assay at weeks 26 and/or 28 compared to NYVAC-C alone. Although 2 individuals discontinued after the first DNA-C due to adverse events (1 vaso-vagal; 1 transient, asymptomatic elevation in alanine transaminase), the vaccines were well tolerated. Three others failed to complete the regimen (1 changed her mind; 2 lost to follow-up). Of the 35 that completed the regimen 90% (18/20) in the DNA-C group had ELISpot responses compared to 33% (5/15) that received NYVAC-C alone (p=0.001). Responses were to envelope in the majority (21/23). Of the 9 individuals with responses to envelope and other peptides, 8 were in the DNA-C group. These promising results suggest that DNA-C was an effective priming agent, that merits further investigation.

Tests pharmacogénétiques: bientôt avant chaque prescription [Pharmacogenetic testing: soon before every prescription?]

Genetic polymorphisms have currently been described in more than 200 systems affecting pharmacological responses (cytochromes P450, conjugation enzymes, transporters, receptors, effectors of response, protection mechanisms, determinants of immunity). Pharmacogenetic testing, i.e. the profiling of individual patients for such variations, is about to become largely available. Recent progress in the pharmacogenetics of tamoxifen, oral anticoagulants and anti-HIV agents is reviewed to discuss critically their potential impact on prescription and contribution/limits for improving rational and safe use of pharmaceuticals. Prospective controlled trials are required to evaluate large-scale pharmacogenetic testing in therapeutics. Ethical, social and psychological issues deserve particular attention.

SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

Molecular basis of selective PPARgamma modulation for the treatment of Type 2 diabetes.

Peroxisome proliferator-activated receptors (PPARs) (alpha, beta/delta and gamma) are lipid sensors capable of adapting gene expression to integrate various lipid signals. As such, PPARs are also very important pharmaceutical targets, and specific synthetic ligands exist for the different isotypes and are either currently used or hold promises in the treatment of major metabolic disorders. In particular, compounds of the class of the thiazolinediones (TZDs) are PPARgamma agonists and potent insulin-sensitizers. The specific but still broad expression patterns of PPARgamma, as well as its implication in numerous pathways, constitutes also a disadvantage regarding drug administration, since this potentially increases the chance to generate side-effects through the activation of the receptor in tissues or cells not affected by the disease. Actually, numerous side effects associated with the administration of TZDs have been reported. Today, a new generation of PPARgamma modulators is being actively developed to activate the receptor more specifically, in a cell and time-dependent manner, in order to induce a specific subset of target genes only and modulate a restricted number of metabolic pathways. We will discuss here why and how the development of such selective PPARgamma modulators is possible, and summarize the results obtained with the published molecules.

Carence en fer sans anémie: où en est-on en 2012 [Iron deficiency without anemia: where are we in 2012?].

Should we treat iron deficiency without anemia? The simple fact that the question can be formulated already leads to controversies. During the past years, the development of a new formulation of intravenous iron has helped fuel the controversy. What is the situation in 2012? This article gives a practical point of view on the actual situation and provides indications on the use of new intravenous medications.

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

In the last few years, a need to account for molecular flexibility in drug-design methodologies has emerged, even if the dynamic behavior of molecular properties is seldom made explicit. For a flexible molecule, it is indeed possible to compute different values for a given conformation-dependent property and the ensemble of such values defines a property space that can be used to describe its molecular variability; a most representative case is the lipophilicity space. In this review, a number of applications of lipophilicity space and other property spaces are presented, showing that this concept can be fruitfully exploited: to investigate the constraints exerted by media of different levels of structural organization, to examine processes of molecular recognition and binding at an atomic level, to derive informative descriptors to be included in quantitative structure--activity relationships and to analyze protein simulations extracting the relevant information. Much molecular information is neglected in the descriptors used by medicinal chemists, while the concept of property space can fill this gap by accounting for the often-disregarded dynamic behavior of both small ligands and biomacromolecules. Property space also introduces some innovative concepts such as molecular sensitivity and plasticity, which appear best suited to explore the ability of a molecule to adapt itself to the environment variously modulating its property and conformational profiles. Globally, such concepts can enhance our understanding of biological phenomena providing fruitful descriptors in drug-design and pharmaceutical sciences.

Averting inflammation by targeting the cytokine environment.

Cytokines are key instigators and regulators of immune responses and therefore hold great potential as targets for new therapeutic strategies. However, the selection of which cytokines to target, and in particular the identification of which cytokines regulate the rate-limiting steps of disease pathways, is crucial to the success of such strategies. Moreover, balancing the need for ablating pathological inflammatory responses and simultaneously maintaining the ability to control infectious agents is a key consideration. Recent advances in our understanding of cytokine networks, as well as technical progress in blocking cytokines in vivo, are likely to be a source for new drugs that can control chronic inflammatory diseases.

Blind docking of 260 protein-ligand complexes with EADock 2.0.

Molecular docking softwares are one of the important tools of modern drug development pipelines. The promising achievements of the last 10 years emphasize the need for further improvement, as reflected by several recent publications (Leach et al., J Med Chem 2006, 49, 5851; Warren et al., J Med Chem 2006, 49, 5912). Our initial approach, EADock, showed a good performance in reproducing the experimental binding modes for a set of 37 different ligand-protein complexes (Grosdidier et al., Proteins 2007, 67, 1010). This article presents recent improvements regarding the scoring and sampling aspects over the initial implementation, as well as a new seeding procedure based on the detection of cavities, opening the door to blind docking with EADock. These enhancements were validated on 260 complexes taken from the high quality Ligand Protein Database [LPDB, (Roche et al., J Med Chem 2001, 44, 3592)]. Two issues were identified: first, the quality of the initial structures cannot be assumed and a manual inspection and/or a search in the literature are likely to be required to achieve the best performance. Second the description of interactions involving metal ions still has to be improved. Nonetheless, a remarkable success rate of 65% was achieved for a large scale blind docking assay, when considering only the top ranked binding mode and a success threshold of 2 A RMSD to the crystal structure. When looking at the five-top ranked binding modes, the success rate increases up to 76%. In a standard local docking assay, success rates of 75 and 83% were obtained, considering only the top ranked binding mode, or the five top binding modes, respectively.

Preparation of alpha-labeled aldehydes by base-catalyzed exchange reactions

A simple, efficient protocol for the preparation of α-labeled aldehydes based on H/D exchange catalyzed by 4-(N,N-dimethylamino)pyridine or Et3N is described. High chemical yields and ratios of isotope incorporation were obtained even when small amounts (1 mmol) of aldehyde were used.

Preparation of alpha-labeled aldehydes by base-catalyzed exchange reactions

A simple, efficient protocol for the preparation of α-labeled aldehydes based on H/D exchange catalyzed by 4-(N,N-dimethylamino)pyridine or Et3N is described. High chemical yields and ratios of isotope incorporation were obtained even when small amounts (1 mmol) of aldehyde were used.

Allosteric conversation in the androgen receptor ligand-binding domain surfaces

Androgen receptor (AR) is a major therapeutic target that plays pivotal roles in prostate cancer (PCa) and androgen insensitivity syndromes. We previously proposed that compounds recruited to ligand-binding domain (LBD) surfaces could regulate AR activity in hormone-refractory PCa and discovered several surface modulators of AR function. Surprisingly, the most effective compounds bound preferentially to a surface of unknown function [binding function 3 (BF-3)] instead of the coactivator-binding site [activation function 2 (AF-2)]. Different BF-3 mutations have been identified in PCa or androgen insensitivity syndrome patients, and they can strongly affect AR activity. Further, comparison of AR x-ray structures with and without bound ligands at BF-3 and AF-2 showed structural coupling between both pockets. Here, we combine experimental evidence and molecular dynamic simulations to investigate whether BF-3 mutations affect AR LBD function and dynamics possibly via allosteric conversation between surface sites. Our data indicate that AF-2 conformation is indeed closely coupled to BF-3 and provide mechanistic proof of their structural interconnection. BF-3 mutations may function as allosteric elicitors, probably shifting the AR LBD conformational ensemble toward conformations that alter AF-2 propensity to reorganize into subpockets that accommodate N-terminal domain and coactivator peptides. The induced conformation may result in either increased or decreased AR activity. Activating BF-3 mutations also favor the formation of another pocket (BF-4) in the vicinity of AF-2 and BF-3, which we also previously identified as a hot spot for a small compound. We discuss the possibility that BF-3 may be a protein-docking site that binds to the N-terminal domain and corepressors. AR surface sites are attractive pharmacological targets to develop allosteric modulators that might be alternative lead compounds for drug design.

Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.