Blind docking of 260 protein-ligand complexes with EADock 2.0.

Molecular docking softwares are one of the important tools of modern drug development pipelines. The promising achievements of the last 10 years emphasize the need for further improvement, as reflected by several recent publications (Leach et al., J Med Chem 2006, 49, 5851; Warren et al., J Med Chem 2006, 49, 5912). Our initial approach, EADock, showed a good performance in reproducing the experimental binding modes for a set of 37 different ligand-protein complexes (Grosdidier et al., Proteins 2007, 67, 1010). This article presents recent improvements regarding the scoring and sampling aspects over the initial implementation, as well as a new seeding procedure based on the detection of cavities, opening the door to blind docking with EADock. These enhancements were validated on 260 complexes taken from the high quality Ligand Protein Database [LPDB, (Roche et al., J Med Chem 2001, 44, 3592)]. Two issues were identified: first, the quality of the initial structures cannot be assumed and a manual inspection and/or a search in the literature are likely to be required to achieve the best performance. Second the description of interactions involving metal ions still has to be improved. Nonetheless, a remarkable success rate of 65% was achieved for a large scale blind docking assay, when considering only the top ranked binding mode and a success threshold of 2 A RMSD to the crystal structure. When looking at the five-top ranked binding modes, the success rate increases up to 76%. In a standard local docking assay, success rates of 75 and 83% were obtained, considering only the top ranked binding mode, or the five top binding modes, respectively.

Different concentrations of Mg++ ions affect nuclear matrix protein distribution during thermal stabilization of isolated nuclei.

The nuclear matrix, a proteinaceous network believed to be a scaffolding structure determining higher-order organization of chromatin, is usually prepared from intact nuclei by a series of extraction steps. In most cell types investigated the nuclear matrix does not spontaneously resist these treatments but must be stabilized before the application of extracting agents. Incubation of isolated nuclei at 37C or 42C in buffers containing Mg++ has been widely employed as stabilizing agent. We have previously demonstrated that heat treatment induces changes in the distribution of three nuclear scaffold proteins in nuclei prepared in the absence of Mg++ ions. We studied whether different concentrations of Mg++ (2.0-5 mM) affect the spatial distribution of nuclear matrix proteins in nuclei isolated from K562 erythroleukemia cells and stabilized by heat at either 37C or 42C. Five proteins were studied, two of which were RNA metabolism-related proteins (a 105-kD component of splicing complexes and an RNP component), one a 126-kD constituent of a class of nuclear bodies, and two were components of the inner matrix network. The localization of proteins was determined by immunofluorescent staining and confocal scanning laser microscope. Mg++ induced significant changes of antigen distribution even at the lowest concentration employed, and these modifications were enhanced in parallel with increase in the concentration of the divalent cation. The different sensitivity to heat stabilization and Mg++ of these nuclear proteins might reflect a different degree of association with the nuclear scaffold and can be closely related to their functional or structural role.

T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension

This article introduces a new interface for T-Coffee, a consistency-based multiple sequence alignment program. This interface provides an easy and intuitive access to the most popular functionality of the package. These include the default T-Coffee mode for protein and nucleic acid sequences, the M-Coffee mode that allows combining the output of any other aligners, and template-based modes of T-Coffee that deliver high accuracy alignments while using structural or homology derived templates. These three available template modes are Expresso for the alignment of protein with a known 3D-Structure, R-Coffee to align RNA sequences with conserved secondary structures and PSI-Coffee to accurately align distantly related sequences using homology extension. The new server benefits from recent improvements of the T-Coffee algorithm and can align up to 150 sequences as long as 10 000 residues and is available from both http://www.tcoffee.org and its main mirror http://tcoffee.crg.cat.

Protein oligomers studied by solid-state NMR the case of the full-length nucleoid-associated protein histone-like nucleoid structuring protein

Members of the histone-like nucleoid structuring protein (H-NS) family play roles both as architectural proteins and as modulators of gene expression in Gram-negative bacteria. The H-NS protein participates in modulatory processes that respond to environmental changes in osmolarity, pH, or temperature. H-NS oligomerization is essential for its activity. Structural models of different truncated forms are available. However, high-resolution structural details of full-length H-NS and its DNA-bound state have largely remained elusive. We report on progress in characterizing the biologically active H-NS oligomers with solid-state NMR. We compared uniformly ((13)C,(15)N)-labeled ssNMR preparations of the isolated N-terminal region (H-NS 1-47) and full-length H-NS (H-NS 1-137). In both cases, we obtained ssNMR spectra of good quality and characteristic of well-folded proteins. Analysis of the results of 2D and 3D (13)C-(13)C and (15)N-(13)C correlation experiments conducted at high magnetic field led to assignments of residues located in different topological regions of the free full-length H-NS. These findings confirm that the structure of the N-terminal dimerization domain is conserved in the oligomeric full-length protein. Small changes in the dimerization interface suggested by localized chemical shift variations between solution and solid-state spectra may be relevant for DNA recoginition.

Vacuolar SNARE Protein Transmembrane Domains Serve as Nonspecific Membrane Anchors with Unequal Roles in Lipid Mixing.

Membrane fusion is induced by SNARE complexes that are anchored in both fusion partners. SNAREs zipper up from the N to C terminus bringing the two membranes into close apposition. Their transmembrane domains (TMDs) might be mere anchoring devices, deforming bilayers by mechanical force. Structural studies suggested that TMDs might also perturb lipid structure by undergoing conformational transitions or by zipping up into the bilayer. Here, we tested this latter hypothesis, which predicts that the activity of SNAREs should depend on the primary sequence of their TMDs. We replaced the TMDs of all vacuolar SNAREs (Nyv1, Vam3, and Vti1) by a lipid anchor, by a TMD from a protein unrelated to the membrane fusion machinery, or by artificial leucine-valine sequences. Individual exchange of the native SNARE TMDs against an unrelated transmembrane anchor or an artificial leucine-valine sequence yielded normal fusion activities. Fusion activity was also preserved upon pairwise exchange of the TMDs against unrelated peptides, which eliminates the possibility for specific TMD-TMD interactions. Thus, a specific primary sequence or zippering beyond the SNARE domains is not a prerequisite for fusion. Lipid-anchored Vti1 was fully active, and lipid-anchored Nyv1 permitted the reaction to proceed up to hemifusion, and lipid-anchored Vam3 interfered already before hemifusion. The unequal contribution of proteinaceous TMDs on Vam3 and Nyv1 suggests that Q- and R-SNAREs might make different contributions to the hemifusion intermediate and the opening of the fusion pore. Furthermore, our data support the view that SNARE TMDs serve as nonspecific membrane anchors in vacuole fusion.

Worldwide Population Structure, Long-Term Demography, and Local Adaptation of Helicobacter pylori.

Helicobacter pylori is an important human pathogen associated with serious gastric diseases. Owing to its medical importance and close relationship with its human host, understanding genomic patterns of global and local adaptation in H. pylori may be of particular significance for both clinical and evolutionary studies. Here we present the first such whole genome analysis of 60 globally distributed strains, from which we inferred worldwide population structure and demographic history and shed light on interesting global and local events of positive selection, with particular emphasis on the evolution of San-associated lineages. Our results indicate a more ancient origin for the association of humans and H. pylori than previously thought. We identify several important perspectives for future clinical research on candidate selected regions that include both previously characterized genes (e.g., transcription elongation factor NusA and tumor necrosis factor alpha-inducing protein Tipα) and hitherto unknown functional genes.

Uncovering the functional constraints underlying the genomic organisation of the Odorant-Binding Protein genes

Animal olfactory systems have a critical role for the survival and reproduction of individuals. In insects, the odorant-binding proteins (OBPs) are encoded by a moderately sized gene family, and mediate the first steps of the olfactory processing. Most OBPs are organized in clusters of a few paralogs, which are conserved over time. Currently, the biological mechanism explaining the close physical proximity among OBPs is not yet established. Here, we conducted a comprehensive study aiming to gain insights into the mechanisms underlying the OBP genomic organization. We found that the OBP clusters are embedded within large conserved arrangements. These organizations also include other non-OBP genes, which often encode proteins integral to plasma membrane. Moreover, the conservation degree of such large clusters is related to the following: 1) the promoter architecture of the confined genes, 2) a characteristic transcriptional environment, and 3) the chromatin conformation of the chromosomal region. Our results suggest that chromatin domains may restrict the location of OBP genes to regions having the appropriate transcriptional environment, leading to the OBP cluster structure. However, the appropriate transcriptional environment for OBP and the other neighbor genes is not dominated by reduced levels of expression noise. Indeed, the stochastic fluctuations in the OBP transcript abundance may have a critical role in the combinatorial nature of the olfactory coding process.

Mechanisms of OGT-mediated HCF-1 protein maturation

Post-translational protein modifications are crucial for many fundamental cellular and extracellular processes and greatly contribute to the complexity of organisms. Human HCF-1 is a transcriptional co-regulator that undergoes complex protein maturation involving reversible and irreversible post-translational modifications. Upon synthesis as a large precursor protein, HCF-1 undergoes extensive reversible glycosylation with β-N-acetylglucosamine giving rise to O-linked-β-N-acetylglucosamine (O-GlcNAc) modified serines and threonines. HCF-1 also undergoes irreversible site-specific proteolysis, which is important for one of HCF-1's major functions - the regulation of the cell-division cycle. HCF-1 O-GlcNAcylation and site-specific proteolysis are both catalyzed by a single enzyme with an unusual dual enzymatic activity, the O-GlcNAc transferase (OGT). HCF-1 is cleaved by OGT at any of six highly conserved 26 amino acid repeated sequences (HCF-1PRO repeats), but the mechanisms and the substrate requirements for OGT-mediated cleavage are not understood. In the present work, I characterized substrate requirements for OGT-mediated cleavage and O-GlcNAcylation of HCF-1. I identified key elements within the HCF-1PRO-repeat sequence that are important for proteolysis. Remarkably, an invariant single amino acid side-chain within the HCF-1PRO-repeat sequence displays particular OGT-binding properties and is essential for proteolysis. Additionally, I characterized substrate requirements for proteolysis outside of the HCF-1PRO repeat and identified a novel, highly O-GlcNAcylated OGT-binding sequence that enhances cleavage of the first HCF-1PRO repeat. These results link OGT association and its O-GlcNAcylation activities to HCF-1PRO-repeat proteolysis.

Genetic variation in the European hake, Merluccius merluccius. Description of protein loci and genetic divergence between Atlantic and Mediterranean populations = Anàlisi de la variabilitat genètica en el lluç europeu, Merluccius merluccius. Descripció dels loci electroforètics i divergència genètica entre les poblacions atlàntica i mediterrània

We examined the genetic population structure of the european hake (Merluccius merluccius) using electrophoretically detectable population markers in 35 protein loci. Samples were collected from 7 locations in the Atlantic Ocean and Mediterranean Sea. Six loci were polymorphic using the 0.05 criterion of polymorphism. Sample heterozigosities ranged from 0.052 to 0.072 and averaged 0.0625. In this study, significant allele frequency differences were detected between Atlantic and Mediterranean populations in three polymorphic loci: GAPDH-1*, GPI-2* and SOD-1*. Two major genetic groups were considered: a North-Atlantic stock and the Mediterranean stock. The Nei genetic distance, D, (based on 33 loci) between samples from these two groups ranged from 0.002 to 0.006. Genetic differenciation between these areas appears to reflect the barrier effect of Strait of Gibraltar. On average over loci, 96.92 % of the total gene diversity was contained within samples, 0.23 % expressed differences among locations within areas, and 2.64 % differences between regions. A review of morphological variation together with the genetic data presented here suggest that the populations of hake from these areas are subdivided into two different stocks: the North-Atlantic stock and the Mediterranean stock. The most conservative approach to the management of these stocks is to consider the Atlantic and Mediterranean stocks independently from oneanother

Structure-Based, Rational Design of T Cell Receptors.

Adoptive cell transfer using engineered T cells is emerging as a promising treatment for metastatic melanoma. Such an approach allows one to introduce T cell receptor (TCR) modifications that, while maintaining the specificity for the targeted antigen, can enhance the binding and kinetic parameters for the interaction with peptides (p) bound to major histocompatibility complexes (MHC). Using the well-characterized 2C TCR/SIYR/H-2K(b) structure as a model system, we demonstrated that a binding free energy decomposition based on the MM-GBSA approach provides a detailed and reliable description of the TCR/pMHC interactions at the structural and thermodynamic levels. Starting from this result, we developed a new structure-based approach, to rationally design new TCR sequences, and applied it to the BC1 TCR targeting the HLA-A2 restricted NY-ESO-1157-165 cancer-testis epitope. Fifty-four percent of the designed sequence replacements exhibited improved pMHC binding as compared to the native TCR, with up to 150-fold increase in affinity, while preserving specificity. Genetically engineered CD8(+) T cells expressing these modified TCRs showed an improved functional activity compared to those expressing BC1 TCR. We measured maximum levels of activities for TCRs within the upper limit of natural affinity, K D = ∼1 - 5 μM. Beyond the affinity threshold at K D < 1 μM we observed an attenuation in cellular function, in line with the "half-life" model of T cell activation. Our computer-aided protein-engineering approach requires the 3D-structure of the TCR-pMHC complex of interest, which can be obtained from X-ray crystallography. We have also developed a homology modeling-based approach, TCRep 3D, to obtain accurate structural models of any TCR-pMHC complexes when experimental data is not available. Since the accuracy of the models depends on the prediction of the TCR orientation over pMHC, we have complemented the approach with a simplified rigid method to predict this orientation and successfully assessed it using all non-redundant TCR-pMHC crystal structures available. These methods potentially extend the use of our TCR engineering method to entire TCR repertoires for which no X-ray structure is available. We have also performed a steered molecular dynamics study of the unbinding of the TCR-pMHC complex to get a better understanding of how TCRs interact with pMHCs. This entire rational TCR design pipeline is now being used to produce rationally optimized TCRs for adoptive cell therapies of stage IV melanoma.

Disassembly of a Medial Transenvelope Structure by Antibiotics during Intracellular Division.

Chlamydiales possess a minimal but functional peptidoglycan precursor biosynthetic and remodeling pathway involved in the assembly of the division septum by an atypical cytokinetic machine and cryptic or modified peptidoglycan-like structure (PGLS). How this reduced cytokinetic machine collectively coordinates the invagination of the envelope has not yet been explored in Chlamydiales. In other Gram-negative bacteria, peptidoglycan provides anchor points that connect the outer membrane to the peptidoglycan during constriction using the Pal-Tol complex. Purifying PGLS and associated proteins from the chlamydial pathogen Waddlia chondrophila, we unearthed the Pal protein as a peptidoglycan-binding protein that localizes to the chlamydial division septum along with other components of the Pal-Tol complex. Together, our PGLS characterization and peptidoglycan-binding assays support the notion that diaminopimelic acid is an important determinant recruiting Pal to the division plane to coordinate the invagination of all envelope layers with the conserved Pal-Tol complex, even during osmotically protected intracellular growth.

Attracting cavities for docking. Replacing the rough energy landscape of the protein by a smooth attracting landscape.

Molecular docking is a computational approach for predicting the most probable position of ligands in the binding sites of macromolecules and constitutes the cornerstone of structure-based computer-aided drug design. Here, we present a new algorithm called Attracting Cavities that allows molecular docking to be performed by simple energy minimizations only. The approach consists in transiently replacing the rough potential energy hypersurface of the protein by a smooth attracting potential driving the ligands into protein cavities. The actual protein energy landscape is reintroduced in a second step to refine the ligand position. The scoring function of Attracting Cavities is based on the CHARMM force field and the FACTS solvation model. The approach was tested on the 85 experimental ligand-protein structures included in the Astex diverse set and achieved a success rate of 80% in reproducing the experimental binding mode starting from a completely randomized ligand conformer. The algorithm thus compares favorably with current state-of-the-art docking programs. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Structure and dynamics of the membrane attaching nitric oxide transporter nitrophorin 7

Nitrophorins represent a unique class of heme proteins that are able to perform the delicate transportation and release of the free-radical gaseous messenger nitric oxide (NO) in a pH-triggered manner. Besides its ability to bind to phospholipid membranes, the N-terminus contains an additional Leu-Pro-Gly stretch, which is a unique sequence trait, and the heme cavity is significantly altered with respect to other nitrophorins. These distinctive features encouraged us to solve the X-ray crystallographic structures of NP7 at low and high pH and bound with different heme ligands (nitric oxide, histamine, imidazole). The overall fold of the lipocalin motif is well preserved in the different X-ray structures and resembles the fold of other nitrophorins. However, a chain-like arrangement in the crystal lattice due to a number of head-to-tail electrostatic stabilizing interactions is found in NP7. Furthermore, the X-ray structures also reveal ligand-dependent changes in the orientation of the heme, as well as in specific interactions between the A-B and G-H loops, which are considered to be relevant for the biological function of nitrophorins. Fast and ultrafast laser triggered ligand rebinding experiments demonstrate the pH-dependent ligand migration within the cavities and the exit route. Finally, the topological distribution of pockets located around the heme as well as from inner cavities present at the rear of the protein provides a distinctive feature in NP7, so that while a loop gated exit mechanism to the solvent has been proposed for most nitrophorins, a more complex mechanism that involves several interconnected gas hosting cavities is proposed for NP7.

Nuclear Matrix in Apoptotic Cell Death and Cell Proliferation. The role of Nuclear Mitotic Apparatus (NuMA) Protein

The nucleus is a membrane enclosed organelle containing most of the genetic information of the cell in the form of chromatin. The nucleus, which can be divided into many sub-organelles such as the nucleoli, the Cajal bodies and the nuclear lamina, is the site for several essential cellular functions such as the DNA replication and its regulation and most of the RNA synthesis and processing. The nucleus is often affected in disease: the size and the shape of the nucleus, the chromatin distribution and the size of the nucleoli have remained the basis for the grading of several cancers. The maintenance of the vertebrate body shape depends on the skeleton. Similarly, in a smaller context, the shape of the cell and the nucleus are mainly regulated by the cytoskeletal and nucleoskeletal elements. The nuclear matrix, which by definition is a detergent, DNase and salt resistant proteinaceous nuclear structure, has been suggested to form the nucleoskeleton responsible for the nuclear integrity. Nuclear mitotic apparatus protein, NuMA, a component of the nuclear matrix, is better known for its mitotic spindle organizing function. NuMA is one of the nuclear matrix proteins suggested to participate in the maintenance of the nuclear integrity during interphase but its interphase function has not been solved to date. This thesis study concentrated on the role of NuMA and the nuclear matrix as structural and functional components of the interphase nucleus. The first two studies clarified the essential role of caspase-3 in the disintegration of the nuclear structures during apoptosis. The second study also showed NuMA and chromatin to co-elute from cells in significant amounts and the apoptotic cleavage of NuMA was clarified to have an important role in the dissociation of NuMA from the chromatin. The third study concentrated on the interphase function of NuMA showing NuMA depletion to result in cell cycle arrest and the cytoplasmic relocalization of NuMA interaction partner GAS41. We suggest that the relocalization of the transcription factor GAS41 may mediate the cell cycle arrest. Thus, this study has given new aspects in the interactions of NuMA, chromatin and the nuclear matrix.

Prediction of 3D structure and ligand-binding properties : ABC transporters and flavodiiron proteins from Synechocystis sp. strain PCC 6803

Cyanobacteria are unicellular, non-nitrogen-fixing prokaryotes, which perform photosynthesis similarly as higher plants. The cyanobacterium Synechocystis sp. strain PCC 6803 is used as a model organism in photosynthesis research. My research described herein aims at understanding the function of the photosynthetic machinery and how it responds to changes in the environment. Detailed knowledge of the regulation of photosynthesis in cyanobacteria can be utilized for biotechnological purposes, for example in the harnessing of solar energy for biofuel production. In photosynthesis, iron participates in electron transfer. Here, we focused on iron transport in Synechocystis sp. strain PCC 6803 and particularly on the environmental regulation of the genes encoding the FutA2BC ferric iron transporter, which belongs to the ABC transporter family. A homology model built for the ATP-binding subunit FutC indicates that it has a functional ATPbinding site as well as conserved interactions with the channel-forming subunit FutB in the transporter complex. Polyamines are important for the cell proliferation, differentiation and apoptosis in prokaryotic and eukaryotic cells. In plants, polyamines have special roles in stress response and in plant survival. The polyamine metabolism in cyanobacteria in response to environmental stress is of interest in research on stress tolerance of higher plants. In this thesis, the potd gene encoding an polyamine transporter subunit from Synechocystis sp. strain PCC 6803 was characterized for the first time. A homology model built for PotD protein indicated that it has capability of binding polyamines, with the preference for spermidine. Furthermore, in order to investigate the structural features of the substrate specificity, polyamines were docked into the binding site. Spermidine was positioned very similarly in Synechocystis PotD as in the template structure and had most favorable interactions of the docked polyamines. Based on the homology model, experimental work was conducted, which confirmed the binding preference. Flavodiiron proteins (Flv) are enzymes, which protect the cell against toxicity of oxygen and/or nitric oxide by reduction. In this thesis, we present a novel type of photoprotection mechanism in cyanobacteria by the heterodimer of Flv2/Flv4. The constructed homology model of Flv2/Flv4 suggests a functional heterodimer capable of rapid electron transfer. The unknown protein sll0218, encoded by the flv2-flv4 operon, is assumed to facilitate the interaction of the Flv2/Flv4 heterodimer and energy transfer between the phycobilisome and PSII. Flv2/Flv4 provides an alternative electron transfer pathway and functions as an electron sink in PSII electron transfer.