Binding of non-canonical peptidoglycan controls Vibrio cholerae broad spectrum racemase activity.

Broad-spectrum amino acid racemases (Bsrs) enable bacteria to generate non-canonical D-amino acids (NCDAAs), whose roles and impact on microbial physiology, including modulation of cell wall structure and dissolution of biofilms, are just beginning to be appreciated. Here we used a diverse array of structural, biochemical and molecular simulation studies to define and characterize how BsrV is post-translationally regulated. We discovered that contrary to Vibrio cholerae alanine racemase AlrV highly compacted active site, BsrV's is broader and can be occupied by cell wall stem peptides. We found that peptidoglycan peptides modified with NCDAAs are better stabilized by BsrV's catalytic cavity and show better inhibitory capacity than canonical muropeptides. Notably, BsrV binding and inhibition can be recapitulated by undigested peptidoglycan sacculi as it exists in the cell. Docking simulations of BsrV binding the peptidoglycan polymer generate a model where the peptide stems are perfectly accommodated and stabilized within each of the dimers active sites. Taking these biochemical and structural data together, we propose that inhibition of BsrV by peptidoglycan peptides underlies a negative regulatory mechanism to avoid excessive NCDAA production. Our results collectively open the door to use "à la carte" synthetic peptides as a tool to modulate DAAs production of Bsr enzymes.

Screening protocol for the identification of modulators by immunofluorescent cell-based assay.

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.

Disulfide Engineered Lipase to Enhance the Catalytic Activity: A Structure-Based Approach on BTL2.

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-1) and up to 6.0-fold higher with 0.01% CTAB (778 IU mg-1), 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.

Discovery of novel triazolo[4,3-b]pyridazin-3-yl-quinoline derivatives as PIM inhibitors.

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.

Antibiotic Capture by Bacterial Lipocalins Uncovers an Extracellular Mechanism of Intrinsic Antibiotic Resistance.

The potential for microbes to overcome antibiotics of different classes before they reach bacterial cells is largely unexplored. Here we show that a soluble bacterial lipocalin produced by Burkholderia cenocepacia upon exposure to sublethal antibiotic concentrations increases resistance to diverse antibiotics in vitro and in vivo These phenotypes were recapitulated by heterologous expression in B. cenocepacia of lipocalin genes from Pseudomonas aeruginosa, Mycobacterium tuberculosis, and methicillin-resistant Staphylococcus aureus Purified lipocalin bound different classes of bactericidal antibiotics and contributed to bacterial survival in vivo Experimental and X-ray crystal structure-guided computational studies revealed that lipocalins counteract antibiotic action by capturing antibiotics in the extracellular space. We also demonstrated that fat-soluble vitamins prevent antibiotic capture by binding bacterial lipocalin with higher affinity than antibiotics. Therefore, bacterial lipocalins contribute to antimicrobial resistance by capturing diverse antibiotics in the extracellular space at the site of infection, which can be counteracted by known vitamins.IMPORTANCE Current research on antibiotic action and resistance focuses on targeting essential functions within bacterial cells. We discovered a previously unrecognized mode of general bacterial antibiotic resistance operating in the extracellular space, which depends on bacterial protein molecules called lipocalins. These molecules are highly conserved in most bacteria and have the ability to capture different classes of antibiotics outside bacterial cells. We also discovered that liposoluble vitamins, such as vitamin E, overcome in vitro and in vivo antibiotic resistance mediated by bacterial lipocalins, providing an unexpected new alternative to combat resistance by using this vitamin or its derivatives as antibiotic adjuvants.

Glycolipid-based TLR4 Modulators and Fluorescent Probes: Rational Design, Synthesis, and Biological Properties.

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.

Protein Kinase C ζ Interacts with a Novel Binding Region of Gαq to Act as a Functional Effector.

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 protein-coupled receptor kinase 2, a known Gαq modulator, led to a complete abrogation of ERK5 stimulation. Finally, we reveal that Gα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.

Molecular Basis of the Functional Differences between Soluble Human Versus Murine MD-2: Role of Val135 in Transfer of Lipopolysaccharide from CD14 to MD-2.

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.

Modulation of toll-like receptor 4. Insights from x-ray crystallography and molecular modeling.

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: automatic ligand flexibility assignment.

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.

A structure-based design of new C2- and C13-substituted taxanes: tubulin binding affinities and extended quantitative structure-activity relationships using comparative binding energy (COMBINE) analysis.

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.

Characterizing conformation changes in proteins through the torsional elastic response.

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.

Semisynthetic peptide-lipase conjugates for improved biotransformations.

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.

CRDOCK: an ultrafast multipurpose protein-ligand docking tool.

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.

Comparative binding energy (COMBINE) analysis supports a proposal for the binding mode of epothilones to ß-tubulin.

The conformational preferences of epothilone A (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.

MM-ISMSA: An Ultrafast and Accurate Scoring Function for Protein-Protein Docking.

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.

gCOMBINE: A graphical user interface to perform structure-based comparative binding energy (COMBINE) analysis on a set of ligand-receptor complexes.

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.