Search

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.

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)

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