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.

Crystal packing modifies ligand binding affinity: the case of aldose reductase.

The relationship between the structures of protein-ligand complexes existing in the crystal and in solution, essential in the case of fragment-based screening by X-ray crystallography (FBS-X), has been often an object of controversy. To address this question, simultaneous co-crystallization and soaking of two inhibitors with different ratios, Fidarestat (FID; K(d) = 6.5 nM) and IDD594 (594; K(d) = 61 nM), which bind to h-aldose reductase (AR), have been performed. The subatomic resolution of the crystal structures allows the differentiation of both inhibitors, even when the structures are almost superposed. We have determined the occupation ratio in solution by mass spectrometry (MS) Occ(FID)/Occ(594) = 2.7 and by X-ray crystallography Occ(FID)/Occ(594) = 0.6. The occupancies in the crystal and in solution differ 4.6 times, implying that ligand binding potency is influenced by crystal contacts. A structural analysis shows that the Loop A (residues 122-130), which is exposed to the solvent, is flexible in solution, and is involved in packing contacts within the crystal. Furthermore, inhibitor 594 contacts the base of Loop A, stabilizing it, while inhibitor FID does not. This is shown by the difference in B-factors of the Loop A between the AR-594 and AR-FID complexes. A stable loop diminishes the entropic energy barrier to binding, favoring 594 versus FID. Therefore, the effect of the crystal environment should be taken into consideration in the X-ray diffraction analysis of ligand binding to proteins. This conclusion highlights the need for additional methodologies in the case of FBS-X to validate this powerful screening technique, which is widely used.

The Oncogene Metadherin Modulates the Apoptotic Pathway Based on the Tumor Necrosis Factor Superfamily Member TRAIL (Tumor Necrosis Factor-related Apoptosis-inducing Ligand) in Breast Cancer.

Metadherin (MTDH), the newly discovered gene, is overexpressed in more than 40% of breast cancers. Recent studies have revealed that MTDH favors an oncogenic course and chemoresistance. With a number of breast cancer cell lines and breast tumor samples, we found that the relative expression of MTDH correlated with tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) sensitivity in breast cancer. In this study, we found that knockdown of endogenous MTDH cells sensitized the MDA-MB-231 cells to TRAIL-induced apoptosis both in vitro and in vivo. Conversely, stable overexpression of MTDH in MCF-7 cells enhanced cell survival with TRAIL treatment. Mechanically, MTDH down-regulated caspase-8, decreased caspase-8 recruitment into the TRAIL death-inducing signaling complex, decreased caspase-3 and poly(ADP-ribose) polymerase-2 processing, increased Bcl-2 expression, and stimulated TRAIL-induced Akt phosphorylation, without altering death receptor status. In MDA-MB-231 breast cancer cells, sensitization to TRAIL upon MTDH down-regulation was inhibited by the caspase inhibitor Z-VAD-fmk (benzyloxycarbonyl-VAD-fluoromethyl ketone), suggesting that MTDH depletion stimulates activation of caspases. In MCF-7 breast cancer cells, resistance to TRAIL upon MTDH overexpression was abrogated by depletion of Bcl-2, suggesting that MTDH-induced Bcl-2 expression contributes to TRAIL resistance. We further confirmed that MTDH may control Bcl-2 expression partly by suppressing miR-16. Collectively, our results point to a protective function of MTDH against TRAIL-induced death, whereby it inhibits the intrinsic apoptosis pathway through miR-16-mediated Bcl-2 up-regulation and the extrinsic apoptosis pathway through caspase-8 down-regulation.

A Modulated Closed Form solution for Quantitative Susceptibility Mapping - A thorough evaluation and comparison to iterative methods based on edge prior knowledge.

The aim of this study is to perform a thorough comparison of quantitative susceptibility mapping (QSM) techniques and their dependence on the assumptions made. The compared methodologies were: two iterative single orientation methodologies minimizing the l2, l1TV norm of the prior knowledge of the edges of the object, one over-determined multiple orientation method (COSMOS) and anewly proposed modulated closed-form solution (MCF). The performance of these methods was compared using a numerical phantom and in-vivo high resolution (0.65mm isotropic) brain data acquired at 7T using a new coil combination method. For all QSM methods, the relevant regularization and prior-knowledge parameters were systematically changed in order to evaluate the optimal reconstruction in the presence and absence of a ground truth. Additionally, the QSM contrast was compared to conventional gradient recalled echo (GRE) magnitude and R2* maps obtained from the same dataset. The QSM reconstruction results of the single orientation methods show comparable performance. The MCF method has the highest correlation (corrMCF=0.95, r(2)MCF =0.97) with the state of the art method (COSMOS) with additional advantage of extreme fast computation time. The l-curve method gave the visually most satisfactory balance between reduction of streaking artifacts and over-regularization with the latter being overemphasized when the using the COSMOS susceptibility maps as ground-truth. R2* and susceptibility maps, when calculated from the same datasets, although based on distinct features of the data, have a comparable ability to distinguish deep gray matter structures.

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

ABSTRACT The drug discovery process has been profoundly changed recently by the adoption of computational methods helping the design of new drug candidates more rapidly and at lower costs. In silico drug design consists of a collection of tools helping to make rational decisions at the different steps of the drug discovery process, such as the identification of a biomolecular target of therapeutical interest, the selection or the design of new lead compounds and their modification to obtain better affinities, as well as pharmacokinetic and pharmacodynamic properties. Among the different tools available, a particular emphasis is placed in this review on molecular docking, virtual high throughput screening and fragment-based ligand design.

Transcriptional Activation and Receptor Dimerization is Affected by Mutations in the NR2E3 Ligand-Binding Domain

Purpose: Mutations in the ligand-binding domain (LBD) of NR2E3 cause recessively inherited enhanced short wavelength sensitive (S-) cone syndrome (ESCS), Goldmann-Favre syndrome (GFS) and clumped pigmentary retinal degeneration (CPRD). In addition to ligand binding, the LBD contains also essential amino acid sequences for the oligomerization of nuclear receptors. The aim of our studies is to characterize the impact of mutations in the LBD on receptor oligomerization and transcriptional activity of NR2E3. Methods: The different NR2E3 mutants were generated by QuickChange mutagenesis and analyzed in 293T-based transactivation studies and BRET2 (bioluminescence resonance electron transfer) assays. In silico homology modeling of mutant proteins was also performed using available crystallographic data of related nuclear receptors. Results: The mutants p.W234S, p.A256V, p.A256E, p.L263P, p.R309G, p.R311Q, p.R334G, p.L336P, p.L353V, p.R385P and p.M407K, all located in the LBD, showed impaired receptor dimerization at various degrees. Impaired repressor dimerization as assessed by BRET2 assays did not always correlate with impaired repressor function of NR2E3 as assessed by cell-based reporter assays. There were minor differences of transcriptional activity of mutant proteins on mouse S-opsin (opn1sw), mouse cone arrestin (arr3) and human cone arrestin, suggesting that the effect of LBD mutations was independent of the promoter context. Conclusions: Mutational analysis and homology modeling allowed the characterization of potential oligomerization interfaces of the NR2E3 LBD. Additionally, mutations in NR2E3 LBD may cause recessive retinal degenerations by different molecular mechanisms.

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.

A soluble hexameric form of CD40 ligand activates human dendritic cells and augments memory T cell response.

A strategy to improve the immunogenicity of candidate vaccines is to trigger the innate immune system. Triggering of CD40 at the surface of dendritic cells (DC) is essential in the induction of an efficient immune response. Although CD40 agonist antibodies have been shown to be potent inducers of immune responses in experimental models, serious safety concerns have been raised for their use in humans. In addition, the production of soluble functional CD40 ligand has been challenging and the soluble form existing so far is not developed anymore. Here, we have evaluated the potency of a new soluble form of hexameric CD40 ligand (sCD40L) to serve as an adjuvant for anti-viral T cell responses. sCD40L was able to activate human DC and to enhance virus-specific memory T cell responses. These results demonstrate that this soluble form of CD40 ligand may serve as an adjuvant for T cell response and thus provide the rationale for its potential use in T cell based vaccine strategies.

SwissBioisostere: a database of molecular replacements for ligand design.

The SwissBioisostere database (http://www.swissbioisostere.ch) contains information on molecular replacements and their performance in biochemical assays. It is meant to provide researchers in drug discovery projects with ideas for bioisosteric modifications of their current lead molecule, as well as to give interested scientists access to the details on particular molecular replacements. As of August 2012, the database contains 21 293 355 datapoints corresponding to 5 586 462 unique replacements that have been measured in 35 039 assays against 1948 molecular targets representing 30 target classes. The accessible data were created through detection of matched molecular pairs and mining bioactivity data in the ChEMBL database. The SwissBioisostere database is hosted by the Swiss Institute of Bioinformatics and available via a web-based interface.

Characterization of the ligand-binding site of the serotonin 5-HT3 receptor: the role of glutamate residues 97, 224, AND 235.

Ligand-gated ion channels of the Cys loop family are receptors for small amine-containing neurotransmitters. Charged amino acids are strongly conserved in the ligand-binding domain of these receptor proteins. To investigate the role of particular residues in ligand binding of the serotonin 5-HT3AS receptor (5-HT3R), glutamate amino acid residues at three different positions, Glu97, Glu224, and Glu235, in the extracellular N-terminal domain were substituted with aspartate and glutamine using site-directed mutagenesis. Wild type and mutant receptor proteins were expressed in HEK293 cells and analyzed by electrophysiology, radioligand binding, fluorescence measurements, and immunochemistry. A structural model of the ligand-binding domain of the 5-HT3R based on the acetylcholine binding protein revealed the position of the mutated amino acids. Our results demonstrate that mutations of Glu97, distant from the ligand-binding site, had little effect on the receptor, whereas mutations Glu224 and Glu235, close to the predicted binding site, are indeed important for ligand binding. Mutations E224Q, E224D, and E235Q decreased EC50 and Kd values 5-20-fold, whereas E235D was functionally expressed at a low level and had a more than 100-fold increased EC50 value. Comparison of the fluorescence properties of a fluorescein-labeled antagonist upon binding to wild type 5-HT3R and E235Q, allowed us to localize Glu235 within a distance of 1 nm around the ligand-binding site, as proposed by our model.

Development of improved soluble inhibitors of FasL and CD40L based on oligomerized receptors.

TNF receptor family members fused to the constant domain of immunoglobulin G have been widely used as immunoadhesins in basic in vitro and in vivo research and in some clinical applications. In this study, we assemble soluble, high avidity chimeric receptors on a pentameric scaffold derived from the coiled-coil domain of cartilage oligomeric matrix protein (COMP). The affinity of Fas and CD40 (but not TNFR-1 and TRAIL-R2) to their ligands is increased by fusion to COMP, when compared to the respective Fc chimeras. In functional assays, Fas:COMP was at least 20-fold more active than Fas:Fc at inhibiting the action of sFasL, and CD40:COMP could block CD40L-mediated proliferation of B cells, whereas CD40:Fc could not. In conclusion, members of the TNF receptor family can display high specificity and excellent avidity for their ligands if they are adequately multimerized.

Tools and techniques to study ligand-receptor interactions and receptor activation by TNF superfamily members.

Ligands and receptors of the TNF superfamily are therapeutically relevant targets in a wide range of human diseases. This chapter describes assays based on ELISA, immunoprecipitation, FACS, and reporter cell lines to monitor interactions of tagged receptors and ligands in both soluble and membrane-bound forms using unified detection techniques. A reporter cell assay that is sensitive to ligand oligomerization can identify ligands with high probability of being active on endogenous receptors. Several assays are also suitable to measure the activity of agonist or antagonist antibodies, or to detect interactions with proteoglycans. Finally, self-interaction of membrane-bound receptors can be evidenced using a FRET-based assay. This panel of methods provides a large degree of flexibility to address questions related to the specificity, activation, or inhibition of TNF-TNF receptor interactions in independent assay systems, but does not substitute for further tests in physiologically relevant conditions.

Purification and ligand binding of a soluble class I major histocompatibility complex molecule consisting of the first three domains of H-2Kd fused to beta 2-microglobulin expressed in the baculovirus-insect cell system.

A recombinant baculovirus encoding a single-chain murine major histocompatibility complex class I molecule in which the first three domains of H-2Kd are fused to beta 2-microglobulin (beta 2-m) via a 15-amino acid linker has been isolated and used to infect lepidopteran cells. A soluble, 391-amino acid single-chain H-2Kd (SC-Kd) molecule of 48 kDa was synthesized and glycosylated in insect cells and could be purified in the absence of detergents by affinity chromatography using the anti-H-2Kd monoclonal antibody SF1.1.1.1. We tested the ability of SC-Kd to bind antigenic peptides using a direct binding assay based on photoaffinity labeling. The photoreactive derivative was prepared from the H-2Kd-restricted Plasmodium berghei circumsporozoite protein (P.b. CS) peptide 253-260 (YIPSAEKI), a probe that we had previously shown to be unable to bind to the H-2Kd heavy chain in infected cells in the absence of co-expressed beta 2-microglobulin. SC-Kd expressed in insect cells did not require additional mouse beta 2-m to bind the photoprobe, indicating that the covalently attached beta 2-m could substitute for the free molecule. Similarly, binding of the P.b. CS photoaffinity probe to the purified SC-Kd molecule was unaffected by the addition of exogenous beta 2-m. This is in contrast to H-2KdQ10, a soluble H-2Kd molecule in which beta 2-m is noncovalently bound to the soluble heavy chain, whose ability to bind the photoaffinity probe is greatly enhanced in the presence of an excess of exogenous beta 2-m. The binding of the probe to SC-Kd was allele-specific, since labeling was selectively inhibited only by antigenic peptides known to be presented by the H-2Kd molecule.

Pi*-pi* bonding interactions generated by halogen oxidation of zirconium(IV) redox-active ligand complexes.

The new complex, [Zr(pda)2]n (1, pda2- = N,N'-bis(neo-pentyl)-ortho-phenylenediamide, n = 1 or 2), prepared by the reaction of 2 equiv of pdaLi2 with ZrCl4, reacts rapidly with halogen oxidants to afford the new product ZrX2(disq)2 (3, X = Cl, Br, I; disq- = N,N'-bis(neo-pentyl)-ortho-diiminosemiquinonate) in which each redox-active ligand has been oxidized by one electron. The oxidation products 3a-c have been structurally characterized and display an unusual parallel stacked arrangement of the disq- ligands in the solid state, with a separation of approximately 3 A. Density functional calculations show a bonding-type interaction between the SOMOs of the disq- ligands to form a unique HOMO while the antibonding linear combination forms a unique LUMO. This orbital configuration leads to a closed-shell-singlet ground-state electron configuration (S = 0). Temperature-dependent magnetism measurements indicate a low-lying triplet excited state at approximately 750 cm-1. In solution, 3a-c show strong disq--based absorption bands that are invariant across the halide series. Taken together these spectroscopic measurements provide experimental values for the one- and two-electron energies that characterize the pi-stacked bonding interaction between the two disq- ligands.

Fatty acids, eicosanoids, and hypolipidemic agents identified as ligands of peroxisome proliferator-activated receptors by coactivator-dependent receptor ligand assay.

Peroxisome proliferator-activated receptors (PPARs) are nuclear hormone receptors controlling the expression of genes involved in lipid homeostasis. PPARs activate gene transcription in response to a variety of compounds including hypolipidemic drugs as well as natural fatty acids. From the plethora of PPAR activators, Scatchard analysis of receptor-ligand interactions has thus far identified only four ligands. These are the chemotactic agent leukotriene B4 and the hypolipidemic drug Wy 14,643 for the alpha-subtype and a prostaglandin J2 metabolite and synthetic antidiabetic thiazolidinediones for the gamma-subtype. Based on the hypothesis that ligand binding to PPAR would induce interactions of the receptor with transcriptional coactivators, we have developed a novel ligand sensor assay, termed coactivator-dependent receptor ligand assay (CARLA). With CARLA we have screened several natural and synthetic candidate ligands and have identified naturally occurring fatty acids and metabolites as well as hypolipidemic drugs as bona fide ligands of the three PPAR subtypes from Xenopus laevis. Our results suggest that PPARs, by their ability to interact with a number of structurally diverse compounds, have acquired unique ligand-binding properties among the superfamily of nuclear receptors that are compatible with their biological activity.