Search

High-throughput assays are a common strategy for the identification of compounds able to modulate a certain cellular activity. Here, we show an automatized analysis platform for the quantification of nuclear foci as inhibitory effect of compounds on a target protein labeled by fluorescent antibodies. Our experience led us to a fast analysis platform that combines cell-based assays, high-content screening, and confocal microscopy, with an automatic and user-friendly statistical analysis of plate-based assays including positive and negative controls, able to identify inhibitory effect of compounds tested together with the Z-prime and Window of individual plate-based assays to assess the reliability of the results. The platform integrates a web-based tool implemented in Pipeline Pilot and R, and allows computing the inhibition values of different parameters obtained from the high-content screening and confocal microscopy analysis. This facilitates the exploration of the results using the different parameters, providing information at different levels as the number of foci observed, the sum of intensity of foci, area of foci, etc, the detection and filtering of outliers over the assay plate, and finally providing a set of statistics of the parameters studied together with a set of plots that we believe significantly helps to the interpretation of the assay results.

16 pags., 4 figs., 3 tabs. -- Open Access funded by Creative Commons Atribution Licence 4.0, Enhancement, control, and tuning of hydrolytic activity and specificity of lipases are major goals for the industry. Thermoalkaliphilic lipases from the I.5 family, with their native advantages such as high thermostability and tolerance to alkaline pHs, are a target for biotechnological applications. Although several strategies have been applied to increase lipases activity, the enhancement through protein engineering without compromising other capabilities is still elusive. Lipases from the I.5 family suffer a unique and delicate double lid restructuration to transition from a closed and inactive state to their open and enzymatically active conformation. In order to increase the activity of the wild type Geobacillus thermocatenulatus lipase 2 (BTL2) we rationally designed, based on its tridimensional structure, a mutant (ccBTL2) capable of forming a disulfide bond to lock the open state. ccBTL2 was generated replacing A191 and F206 to cysteine residues while both wild type C64 and C295 were mutated to serine. A covalently immobilized ccBTL2 showed a 3.5-fold increment in esterase activity with 0.1% Triton X-100 (2336 IU mg) and up to 6.0-fold higher with 0.01% CTAB (778 IU mg), both in the presence of oxidizing sulfhydryl agents, when compared to BTL2. The remarkable and industrially desired features of BTL2 such as optimal alkaliphilic pH and high thermal stability were not affected. The designed disulfide bond also conferred reversibility to the enhancement, as the increment on activity observed for ccBTL2 was controlled by redox pretreatments. MD simulations suggested that the most stable conformation for ccBTL2 (with the disulfide bond formed) was, as we predicted, similar to the open and active conformation of this lipase., Financial and logistic support from Colombian Universidad del Valle and COLCIENCIAS (CI 71083—Grant 745-2016-Project 110671250425), Spanish CICYT project BIO-2005-6018576, BFU2017-90030-P, and BFU2011-25326, B. Di G. In addition, thanks to the Spanish MINECO for a FPU fellowship.

Myeloid differentiation factor 2 (MD-2) is an extracellular protein, associated with the ectodomain of TLR4, that plays a critical role in the recognition of bacterial LPS. Despite high overall structural and functional similarity, human (h) and murine (m) MD-2 exhibit several species-related differences. hMD-2 is capable of binding LPS in the absence of TLR4, whereas mMD-2 supports LPS responsiveness only when mMD-2 and mTLR4 are coexpressed in the same cell. Previously, charged residues at the edge of the LPS binding pocket have been attributed to this difference. In this study, site-directed mutagenesis was used to explore the hydrophobic residues within the MD-2 binding pocket as the source of functional differences between hMD-2 and mMD-2. Whereas decreased hydrophobicity of residues 61 and 63 in the hMD-2 binding pocket retained the characteristics of wild-type hMD-2, a relatively minor change of valine to alanine at position 135 completely abolished the binding of LPS to the hMD-2 mutant. The mutant, however, retained the LPS binding in complex with TLR4 and also cell activation, resulting in a murine-like phenotype. These results were supported by the molecular dynamics simulation. We propose that the residue at position 135 of MD-2 governs the dynamics of the binding pocket and its ability to accommodate lipid A, which is allosterically affected by bound TLR4.

PIM kinase family (PIM-1, PIM-2 and PIM-3) is an appealing target for the discovery and development of selective inhibitors, useful in various disease conditions in which these proteins are highly expressed, such as cancer. The significant effort put, in the recent years, towards the development of small molecules exhibiting inhibitory activity against this protein family has ended up with several molecules entering clinical trials. As part of our ongoing exploration for potential drug candidates that exhibit affinity towards this protein family, we have generated a novel chemical series of triazolo[4,3-b]pyridazine based tricycles by applying a scaffold hopping strategy over our previously reported potent pan-PIM inhibitor ETP-47453 (compound 2). The structure-activity relationship studies presented herein demonstrate a rather selective PIM-1/PIM-3 biochemical profile for this novel series of tricycles, although pan-PIM and PIM-1 inhibitors have also been identified. Selected examples show significant inhibition of the phosphorylation of BAD protein in a cell-based assay. Moreover, optimized and highly selective compounds, such as 42, did not show significant hERG inhibition at 20 μM concentration, and proved its antiproliferative activity and utility in combination with particular antitumoral agents in several tumor cell lines.

Toll-like receptors (TLRs) are a family of proteins with a key role in the innate immune system. They are specialized in the recognition of molecular patterns present in microbial components, through mechanisms not yet unraveled at atomic level. Improvement in the understanding of the molecular mechanisms that drive TLR signaling is of paramount importance to grasp key aspects of immunity, potentially leading to the design of new molecules able to modulate their functions. Toll-like receptor 4 (TLR4), along with its accessory protein myeloid differentiation factor 2 (MD-2), builds a heterodimeric complex that specifically recognizes lipopolysaccharides (LPS), which are present on the cell wall of gramnegative bacteria, activating the immune response. Some TLR4 modulators are undergoing preclinical and clinical evaluation for the treatment of sepsis, inflammatory diseases, cancer, and rheumatoid arthritis. Reported X-ray crystal structures together with molecular modeling studies, not reviewed before in the literature, have recently contributed to the elucidation of key interactions at atomic level of the binding between the TLR4/MD-2 system and different TLR4/MD-2 ligands. The purpose of this review is to summarize these reported studies which may account for the SAR rationalization of natural/ synthetic agonist/antagonist TLR4 binders and may also guide further design of novel TLR4 modulators.

ALFA is a fast computational tool for the conformational analysis of small molecules that uses a custom-made iterative algorithm to provide a set of representative conformers in an attempt to reproduce the diversity of states in which small molecules can exist, either isolated in solution or bound to a target. The results shown in this work prove that ALFA is fast enough to be integrated into massive high-throughput virtual screening protocols with the aim of incorporating ligand flexibility and also that ALFA reproduces crystallographic X-ray structures of bound ligands with great accuracy. Furthermore, the application includes a graphical user interface that allows its use through the popular molecular graphics program PyMOL to make it accessible to nonexpert users. ALFA is distributed free of charge upon request from the authors.

Heterotrimeric G proteins play an essential role in the initiation of G protein-coupled receptor (GPCR) signaling through specific interactions with a variety of cellular effectors. We have recently reported that GPCR activation promotes a direct interaction between Gαq and protein kinase C ζ (PKCζ), leading to the stimulation of the ERK5 pathway independent of the canonical effector PLCβ. We report herein that the activation-dependent Gαq/PKCζ complex involves the basic PB1-type II domain of PKCζ and a novel interaction module in Gαq different from the classical effector-binding site. Point mutations in this Gαq region completely abrogate ERK5 phosphorylation, indicating that Gαq/PKCζ association is required for the activation of the pathway. Indeed, PKCζ was demonstrated to directly bind ERK5 thus acting as a scaffold between Gαq and ERK5 upon GPCR activation. The inhibition of these protein complexes by G proteincoupled receptor kinase 2, a known Gαq modulator, led to a complete abrogation of ERK5 stimulation. Finally, we reveal thatGαq/ PKCζ complexes link Gαq to apoptotic cell death pathways. Our data suggest that the interaction between this novel region in Gαq and the effector PKCζ is a key event in Gαq signaling., Ministerio de Educación y Ciencia (SAF2011-23800, SAF2014-55511-R), Fundación Ramón Areces, The Cardiovascular Diseases Network of Ministerio Sanidad y Consumo-Instituto Carlos III (RD12/0042/0012), Comunidad de Madrid (S-2011/BMD-2332), and Instituto de Salud Carlos III (PI11/00126, PI14/00201) (to F. M. and C. R.). This work was also supported in part by the NIGMS, National Institutes of Health Grant R01-GM088242 (to G. T.), the Canadian Institutes of Health Research (CIHR) (MOP-GMX-231013) (to S. M.), an EMBO Short Fellowship (to G. S. F.), and Fondo Europeo de Desarrollo Regional (FEDER, European Union)

We present gCOMBINE, a Java-written graphical user interface (GUI) for performing comparative binding energy (COMBINE) analysis (Ortiz et al. J Med Chem 1995; 38:2681-2691) on a set of ligand-receptor complexeswith the aim of deriving highly informative quantitative structure-activity relationships. The essence of the method is to decompose the ligand-receptor interaction energies into a series of terms, explore the origins of the variance within the set using Principal Component Analysis, and then assign weights to selected ligandresidue interactions using partial least squares analysis to correlate with the experimental activities or binding affinities. The GUI allows plenty of interactivity and provides multiple plots representing the energy descriptors entering the analysis, scores, loadings, experimental versus predicted regression lines, and the evolution of parameterssuch as r(2) (correlation coefficient), q(2) (cross-validated r(2)), and prediction errors as the number of extracted latent variables increases. Other representative features include the implementation of a sigmoidal dielectric function for electrostatic energy calculations, alternative cross-validation procedures (leave-N-out and random groups), drawing of confidence ellipses, and the possibility to carry out several additional tasks such as optional truncation of positive interaction energy values and generation of ready-to-use PDB files containing information related to the importance for activity of individual protein residues. This information can be displayed and color-coded using a standard molecular graphics program such as PyMOL. It is expected that this user-friendly tool will expand the applicability of the COMBINE analysis method and encourage more groups to use it in their drug design research programs.

The cationic glycolipid IAXO-102, a potent TLR4 antagonist targeting both MD-2 and CD14 co-receptors, has been used as scaffold to design new potential TLR4 modulators and fluorescent labels for the TLR4 receptor complex (membrane TLR4.MD-2 dimer and CD14). The primary amino group of IAXO-102, not involved in direct interaction with MD-2 and CD14 receptors, has been exploited to covalently attach a fluorescein (molecules 1 and 2) or to link two molecules of IAXO-102 through diamine and diammonium spacers, obtaining 'dimeric' molecules 3 and 4. The structure-based rational design of compounds 1-4 was guided by the optimization of MD-2 and CD14 binding. Compounds 1 and 2 inhibited TLR4 activation, in a concentration-dependent manner, and signaling in HEK-Blue TLR4 cells. The fluorescent labeling of murine macrophages by molecule 1 was inhibited by LPS and was also abrogated when cell surface proteins were digested by trypsin, thus suggesting an interaction of fluorescent probe 1 with membrane proteins of the TLR4 receptor system.

18 pags, 8 figs, 1 tab, Streptococcus pneumoniae is not only a commensal of the nasopharyngeal epithelium, but may also cause life-threatening diseases. Immune- electron microscopy studies revealed that the bacterial glycolytic enzyme, phosphoglycerate kinase (PGK), is localised on the pneumococcal surface of both capsulated and non-capsulated strains and colocalises with plasminogen. Since pneumococci may concentrate host plasminogen (PLG) together with its activators on the bacterial cell surface to facilitate the formation of plasmin, the involvement of PGK in this process was studied. Specific binding of human or murine PLG to strain-independent PGK was documented, and surface plasmon resonance analyses indicated a high affinity interaction with the kringle domains 1-4 of PLG. Crystal structure determination of pneumococcal PGK together with peptide array analysis revealed localisation of PLG-binding site in the N-terminal region and provided structural motifs for the interaction with PLG. Based on structural analysis data, a potential interaction of PGK with tissue plasminogen activator (tPA) was proposed and experimentally confirmed by binding studies, plasmin activity assays and thrombus degradation analyses. © Schattauer 2014., The research leading to these results has received funding from the European Community’s Seventh Framework Program under Grant Agreement no. HEALTH-F3–2009–223111. This work was also supported by grants from the Spanish Ministry of Economy and Competitiveness (BFU2011–25326) and Comunidad Autónoma de Madrid (CAM) S2010-BMD-2457 (BIPEDD2). J.K. is funded by a grant from Ministerio de Economía y Competitividad (BFU2011–24595). A.M. also acknowledges CAM for financial support to the Fundación Severo Ochoa through the AMAROUTO program

Ten novel taxanes bearing modifications at the C2 and C13 positions of the baccatin core have been synthesized and their binding affinities for mammalian tubulin have been experimentally measured. The design strategy was guided by (i) calculation of interaction energy maps with carbon, nitrogen and oxygen probes within the taxane-binding site of β-tubulin, and (ii) the prospective use of a structure-based QSAR (COMBINE) model derived from an earlier series comprising 47 congeneric taxanes. The tubulin-binding affinity displayed by one of the new compounds (CTX63) proved to be higher than that of docetaxel, and an updated COMBINE model provided a good correlation between the experimental binding free energies and a set of weighted residue-based ligand-receptor interaction energies for 54 out of the 57 compounds studied. The remaining three outliers from the original training series have in common a large unfavourable entropic contribution to the binding free energy that we attribute to taxane preorganization in aqueous solution in a conformation different from that compatible with tubulin binding. Support for this proposal was obtained from solution NMR experiments and molecular dynamics simulations in explicit water. Our results shed additional light on the determinants of tubulin-binding affinity for this important class of antitumour agents and pave the way for further rational structural modifications. © The Royal Society of Chemistry 2013., National Natural Science Foundation of China (NSFC grant 30930108 to W.S.F.); Spanish Ministerio de Ciencia e Innovación (BIO2008-04384 and BIO-2010-16351); Comisión Interministerial de Ciencia y Tecnología (SAF2009-13914-C02-02), and Comunidad de Madrid (S-BIO/0214/2006 and S2010-BMD-2457)

The relationship between functional conformation changes and thermal dynamics of proteins is investigated with the help of the torsional network model (TNM), an elastic network model in torsion angle space that we recently introduced. We propose and test a null-model of >random> conformation changes that assumes that the contributions of normal modes to conformation changes are proportional to their contributions to thermal fluctuations. Deviations from this null model are generally small. When they are large and significant, they consist in conformation changes that are represented by very few low frequency normal modes and overcome small energy barriers. We interpret these features as the result of natural selection favoring the intrinsic protein dynamics consistent with functional conformation changes. These >selected> conformation changes are more frequently associated to ligand binding, and in particular phosphorylation, than to pairs of conformations with the same ligands. This deep relationship between the thermal dynamics of a protein, represented by its normal modes, and its functional dynamics can reconcile in a unique framework the two models of conformation changes, conformational selection and induced fit. The program TNM that computes torsional normal modes and analyzes conformation changes is available upon request. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly., Comunidad de Madrid; the Spanish CSIC; Spanish Ministry of Science (Consolider grants CSD2006-00023 and BFU2011-24595)

An ultrafast docking and virtual screening program, CRDOCK, is presented that contains (1) a search engine that can use a variety of sampling methods and an initial energy evaluation function, (2) several energy minimization algorithms for fine tuning the binding poses, and (3) different scoring functions. This modularity ensures the easy configuration of custom-made protocols that can be optimized depending on the problem in hand. CRDOCK employs a precomputed library of ligand conformations that are initially generated from one-dimensional SMILES strings. Testing CRDOCK on two widely used benchmarks, the ASTEX diverse set and the Directory of Useful Decoys, yielded a success rate of ∼75% in pose prediction and an average AUC of 0.66. A typical ligand can be docked, on average, in just ∼13 s. Extension to a representative group of pharmacologically relevant G protein-coupled receptors that have been recently cocrystallized with some selective ligands allowed us to demonstrate the utility of this tool and also highlight some current limitations. CRDOCK is now included within VSDMIP, our integrated platform for drug discovery. © 2012 American Chemical Society., CICYT (SAF2009-13914-C02-02); Comunidad Autónoma de Madrid; AMAROUTO; Ministerio de Economia y Competitividad

An efficient chemoselective method for the creation of semisynthetic lipases by site-specific incorporation of tailor-made peptides on the lipase-lid site was developed. These new enzymes showed excellent improved specificity and regio- or enantioselectivity in different biotransformations. © 2012 The Royal Society of Chemistry., Spanish National Research Council (CSIC) and by grants BFU2011-25326, BFU2011-24595

An ultrafast and accurate scoring function for protein-protein docking is presented. It includes (1) a molecular mechanics (MM) part based on a 12-6 Lennard-Jones potential; (2) an electrostatic component based on an implicit solvent model (ISM) with individual desolvation penalties for each partner in the protein-protein complex plus a hydrogen bonding term; and (3) a surface area (SA) contribution to account for the loss of water contacts upon protein-protein complex formation. The accuracy and performance of the scoring function, termed MM-ISMSA, have been assessed by (1) comparing the total binding energies, the electrostatic term, and its components (charge-charge and individual desolvation energies), as well as the per residue contributions, to results obtained with well-established methods such as APBSA or MM-PB(GB)SA for a set of 1242 decoy protein-protein complexes and (2) testing its ability to recognize the docking solution closest to the experimental structure as that providing the most favorable total binding energy. For this purpose, a test set consisting of 15 protein-protein complexes with known 3D structure mixed with 10 decoys for each complex was used. The correlation between the values afforded by MM-ISMSA and those from the other methods is quite remarkable (r 2 ∼ 0.9), and only 0.2-5.0 s (depending on the number of residues) are spent on a single calculation including an all vs all pairwise energy decomposition. On the other hand, MM-ISMSA correctly identifies the best docking solution as that closest to the experimental structure in 80% of the cases. Finally, MM-ISMSA can process molecular dynamics trajectories and reports the results as averaged values with their standard deviations. MM-ISMSA has been implemented as a plugin to the widely used molecular graphics program PyMOL, although it can also be executed in command-line mode. MM-ISMSA is distributed free of charge to nonprofit organizations. © 2012 American Chemical Society., CICYT (SAF2009-13914-C02-02); Comunidad Autónoma de Madrid; AMAROUTO; Ministerio de Economía y Competividad; Ministerio de Educación y Ciencia

The conformational preferences of epothiloneA (EPA) and a 12,13-cyclopropyl C12-epimerized analogue were explored in aqueous solution using molecular dynamics simulations. The simulated conformers that provided an optimal fit in the paclitaxel binding site of mammalian β-tubulin were then selected. The resulting modeled complexes were simulated before and after refinement of the M-loop to improve the fitting and assess ligand stability within the binding pocket. The tubulin-bound conformation of EPA was found to be unlike a previously reported solution obtained through mixed crystallographic/NMR/modeling studies. However, our conformation was in agreement with an NMR-based proposal although the exact binding pose within the site was different. Molecular models were built for the complexes of 14 epothilone derivatives with β-tubulin. A projection to latent structures regression method succeeded in providing a good prediction of the experimentally measured binding enthalpies for the whole set of ligands by assigning weights to a selection of interaction energy terms. These receptor-based, quantitative structure-activity relationships support the proposed binding mode, help confirm and interpret previously acquired experimental data, shed additional light on the effect of several β-tubulin mutations on ligand binding, and can potentially direct further experimental studies. © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim., Comisión Interministerial de Ciencia y Tecnología de España; Comunidad de Madrid; Ministerio de Ciencia e Innovación