Aurélien Grosdidier

ExPASy: SIB bioinformatics resource portal

ExPASy (http://www.expasy.org) has worldwide reputation as one of the main bioinformatics resources for proteomics. It has now evolved, becoming an extensible and integrative portal accessing many scientific resources, databases and software tools in different areas of life sciences. Scientists can henceforth access seamlessly a wide range of resources in many different domains, such as proteomics, genomics, phylogeny/evolution, systems biology, population genetics, transcriptomics, etc. The individual resources (databases, web-based and downloadable software tools) are hosted in a ‘decentralized’ way by different groups of the SIB Swiss Institute of Bioinformatics and partner institutions. Specifically, a single web portal provides a common entry point to a wide range of resources developed and operated by different SIB groups and external institutions. The portal features a search function across ‘selected’ resources. Additionally, the availability and usage of resources are monitored. The portal is aimed for both expert users and people who are not familiar with a specific domain in life sciences. The new web interface provides, in particular, visual guidance for newcomers to ExPASy.

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.

SwissDock, a protein-small molecule docking web service based on EADock DSS

Most life science processes involve, at the atomic scale, recognition between two molecules. The prediction of such interactions at the molecular level, by so-called docking software, is a non-trivial task. Docking programs have a wide range of applications ranging from protein engineering to drug design. This article presents SwissDock, a web server dedicated to the docking of small molecules on target proteins. It is based on the EADock DSS engine, combined with setup scripts for curating common problems and for preparing both the target protein and the ligand input files. An efficient Ajax/HTML interface was designed and implemented so that scientists can easily submit dockings and retrieve the predicted complexes. For automated docking tasks, a programmatic SOAP interface has been set up and template programs can be downloaded in Perl, Python and PHP. The web site also provides an access to a database of manually curated complexes, based on the Ligand Protein Database. A wiki and a forum are available to the community to promote interactions between users. The SwissDock web site is available online at http://www.swissdock.ch. We believe it constitutes a step toward generalizing the use of docking tools beyond the traditional molecular modeling community.

The endocrine disruptor monoethyl-hexyl-phthalate is a selective peroxisome proliferator-activated receptor gamma modulator that promotes adipogenesis

The ability of pollutants to affect human health is a major concern, justified by the wide demonstration that reproductive functions are altered by endocrine disrupting chemicals. The definition of endocrine disruption is today extended to broader endocrine regulations, and includes activation of metabolic sensors, such as the peroxisome proliferator-activated receptors (PPARs). Toxicology approaches have demonstrated that phthalate plasticizers can directly influence PPAR activity. What is now missing is a detailed molecular understanding of the fundamental basis of endocrine disrupting chemical interference with PPAR signaling. We thus performed structural and functional analyses that demonstrate how monoethyl-hexyl-phthalate (MEHP) directly activates PPARγ and promotes adipogenesis, albeit to a lower extent than the full agonist rosiglitazone. Importantly, we demonstrate that MEHP induces a selective activation of different PPARγ target genes. Chromatin immunoprecipitation and fluorescence microscopy in living cells reveal that this selective activity correlates with the recruitment of a specific subset of PPARγ coregulators that includes Med1 and PGC-1α, but not p300 and SRC-1. These results highlight some key mechanisms in metabolic disruption but are also instrumental in the context of selective PPAR modulation, a promising field for new therapeutic development based on PPAR modulation.

Fast docking using the CHARMM force field with EADock DSS

The prediction of binding modes (BMs) occurring between a small molecule and a target protein of biological interest has become of great importance for drug development. The overwhelming diversity of needs leaves room for docking approaches addressing specific problems. Nowadays, the universe of docking software ranges from fast and user friendly programs to algorithmically flexible and accurate approaches. EADock2 is an example of the latter. Its multiobjective scoring function was designed around the CHARMM22 force field and the FACTS solvation model. However, the major drawback of such a software design lies in its computational cost. EADock dihedral space sampling (DSS) is built on the most efficient features of EADock2, namely its hybrid sampling engine and multiobjective scoring function. Its performance is equivalent to that of EADock2 for drug‐like ligands, while the CPU time required has been reduced by several orders of magnitude. This huge improvement was achieved through a combination of several innovative features including an automatic bias of the sampling toward putative binding sites, and a very efficient tree‐based DSS algorithm. When the top‐scoring prediction is considered, 57% of BMs of a test set of 251 complexes were reproduced within 2 Å RMSD to the crystal structure. Up to 70% were reproduced when considering the five top scoring predictions. The success rate is lower in cross‐docking assays but remains comparable with that of the latest version of AutoDock that accounts for the protein flexibility.

Rational design of indoleamine 2,3-dioxygenase inhibitors

Indoleamine 2,3-dioxygenase (IDO) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. We have used the evolutionary docking algorithm EADock to design new inhibitors of this enzyme. First, we investigated the modes of binding of all known IDO inhibitors. On the basis of the observed docked conformations, we developed a pharmacophore model, which was then used to devise new compounds to be tested for IDO inhibition. We also used a fragment-based approach to design and to optimize small organic molecule inhibitors. Both approaches yielded several new low-molecular weight inhibitor scaffolds, the most active being of nanomolar potency in an enzymatic assay. Cellular assays confirmed the potential biological relevance of four different scaffolds.

Docking, virtual high throughput screening and in silico fragment-based drug design.

•  Introduction •  Docking ‐  Protein and ligand flexibility ‐  Sampling algorithm ‐  Scoring functions ‐  Performance •  Virtual high throughput screening ‐  Ligand‐based vHTS ‐  Structure‐based vHTS •  In silico fragment‐based drug design ‐  Theoretical advantages of FBD ‐  Existing FBD methods ‐  Synthetic accessibility of molecules proposed by FBD •  Conclusion

EADock: Docking of small molecules into protein active sites with a multiobjective evolutionary optimization

In recent years, protein–ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure‐based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein–ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 Å around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 Å root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 Å, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the αVβ3 integrin, starting far away from the binding pocket. Proteins 2007.

SwissTargetPrediction: a web server for target prediction of bioactive small molecules

Bioactive small molecules, such as drugs or metabolites, bind to proteins or other macro-molecular targets to modulate their activity, which in turn results in the observed phenotypic effects. For this reason, mapping the targets of bioactive small molecules is a key step toward unraveling the molecular mechanisms underlying their bioactivity and predicting potential side effects or cross-reactivity. Recently, large datasets of protein–small molecule interactions have become available, providing a unique source of information for the development of knowledge-based approaches to computationally identify new targets for uncharacterized molecules or secondary targets for known molecules. Here, we introduce SwissTargetPrediction, a web server to accurately predict the targets of bioactive molecules based on a combination of 2D and 3D similarity measures with known ligands. Predictions can be carried out in five different organisms, and mapping predictions by homology within and between different species is enabled for close paralogs and orthologs. SwissTargetPrediction is accessible free of charge and without login requirement at http://www.swisstargetprediction.ch.

Rational design of 4-aryl-1,2,3-triazoles for indoleamine 2,3-dioxygenase 1 inhibition.

Indoleamine 2,3-dioxygenase 1 (IDO1) is an important therapeutic target for the treatment of diseases such as cancer that involve pathological immune escape. Starting from the scaffold of our previously discovered IDO1 inhibitor 4-phenyl-1,2,3-triazole, we used computational structure-based methods to design more potent ligands. This approach yielded highly efficient low molecular weight inhibitors, the most active being of nanomolar potency both in an enzymatic and in a cellular assay, while showing no cellular toxicity and a high selectivity for IDO1 over tryptophan 2,3-dioxygenase (TDO). A quantitative structure-activity relationship based on the electrostatic ligand-protein interactions in the docked binding modes and on the quantum chemically derived charges of the triazole ring demonstrated a good explanatory power for the observed activities.