Biophysical characterization of tubulin tyrosine ligase-like proteins on microtubule dynamics

Report for the scientific sojourn carried out at the Cell Biology and Biophysics Unit from the National Institutes of Health, from 2010 to 2012.

Distributed all-atom simulations of intrinsically unstructured proteins

The Computational Biophysics Group at the Universitat Pompeu Fabra (GRIB-UPF) hosts two unique computational resources dedicated to the execution of large scale molecular dynamics (MD) simulations: (a) the ACMD molecular-dynamics software, used on standard personal computers with graphical processing units (GPUs); and (b) the GPUGRID. net computing network, supported by users distributed worldwide that volunteer GPUs for biomedical research. We leveraged these resources and developed studies, protocols and open-source software to elucidate energetics and pathways of a number of biomolecular systems, with a special focus on flexible proteins with many degrees of freedom. First, we characterized ion permeation through the bactericidal model protein Gramicidin A conducting one of the largest studies to date with the steered MD biasing methodology. Next, we addressed an open problem in structural biology, the determination of drug-protein association kinetics; we reconstructed the binding free energy, association, and dissaciociation rates of a drug like model system through a spatial decomposition and a Makov-chain analysis. The work was published in the Proceedings of the National Academy of Sciences and become one of the few landmark papers elucidating a ligand-binding pathway. Furthermore, we investigated the unstructured Kinase Inducible Domain (KID), a 28-peptide central to signalling and transcriptional response; the kinetics of this challenging system was modelled with a Markovian approach in collaboration with Frank Noe’s group at the Freie University of Berlin. The impact of the funding includes three peer-reviewed publication on high-impact journals; three more papers under review; four MD analysis components, released as open-source software; MD protocols; didactic material, and code for the hosting group.

Digestion of Streptococcus pneumoniae cell walls with its major peptidoglycan hydrolase releases branched stem peptides carrying proinflammatory activity.

The peptidoglycan of Gram-positive bacteria is known to trigger cytokine release from peripheral blood mononuclear cells (PBMCs). However, it requires 100-1000 times more Gram-positive peptidoglycan than Gram-negative lipopolysaccharide to release the same amounts of cytokines from target cells. Thus, either peptidoglycan is poorly active or only part of it is required for PBMC activation. To test this hypothesis, purified Streptococcus pneumoniae walls were digested with their major autolysin N-acetylmuramoyl-L-alanine amidase, and/or muramidase. Solubilized walls were separated by reverse phase high pressure chromatography. Individual fractions were tested for their PBMC-stimulating activity, and their composition was determined. Soluble components had a Mr between 600 and 1500. These primarily comprised stem peptides cross-linked to various extents. Simple stem peptides (Mr <750) were 10-fold less active than undigested peptidoglycan. In contrast, tripeptides (Mr >1000) were >/=100-fold more potent than the native material. One dipeptide (inactive) and two tripeptides (active) were confirmed by post-source decay analysis. Complex branched peptides represented </=2% of the total material, but their activity (w/w) was almost equal to that of LPS. This is the first observation suggesting that peptidoglycan stem peptides carry high tumor necrosis factor-stimulating activity. These types of structures are conserved among Gram-positive bacteria and will provide new material to help elucidate the mechanism of peptidoglycan-induced inflammation.

Docking to heme proteins.

In silico screening has become a valuable tool in drug design, but some drug targets represent real challenges for docking algorithms. This is especially true for metalloproteins, whose interactions with ligands are difficult to parametrize. Our docking algorithm, EADock, is based on the CHARMM force field, which assures a physically sound scoring function and a good transferability to a wide range of systems, but also exhibits difficulties in case of some metalloproteins. Here, we consider the therapeutically important case of heme proteins featuring an iron core at the active site. Using a standard docking protocol, where the iron-ligand interaction is underestimated, we obtained a success rate of 28% for a test set of 50 heme-containing complexes with iron-ligand contact. By introducing Morse-like metal binding potentials (MMBP), which are fitted to reproduce density functional theory calculations, we are able to increase the success rate to 62%. The remaining failures are mainly due to specific ligand-water interactions in the X-ray structures. Testing of the MMBP on a second data set of non iron binders (14 cases) demonstrates that they do not introduce a spurious bias towards metal binding, which suggests that they may reliably be used also for cross-docking studies.

Splicing and the evolution of proteins in mammals.

It is often supposed that a protein's rate of evolution and its amino acid content are determined by the function and anatomy of the protein. Here we examine an alternative possibility, namely that the requirement to specify in the unprocessed RNA, in the vicinity of intron-exon boundaries, information necessary for removal of introns (e.g., exonic splice enhancers) affects both amino acid usage and rates of protein evolution. We find that the majority of amino acids show skewed usage near intron-exon boundaries, and that differences in the trends for the 2-fold and 4-fold blocks of both arginine and leucine show this to be owing to effects mediated at the nucleotide level. More specifically, there is a robust relationship between the extent to which an amino acid is preferred/avoided near boundaries and its enrichment/paucity in splice enhancers. As might then be expected, the rate of evolution is lowest near intron-exon boundaries, at least in part owing to splice enhancers, such that domains flanking intron-exon junctions evolve on average at under half the rate of exon centres from the same gene. In contrast, the rate of evolution of intronless retrogenes is highest near the domains where intron-exon junctions previously resided. The proportion of sequence near intron-exon boundaries is one of the stronger predictors of a protein's rate of evolution in mammals yet described. We conclude that after intron insertion selection favours modification of amino acid content near intron-exon junctions, so as to enable efficient intron removal, these changes then being subject to strong purifying selection even if nonoptimal for protein function. Thus there exists a strong force operating on protein evolution in mammals that is not explained directly in terms of the biology of the protein.

KnotProt: a database of proteins with knots and slipknots.

The protein topology database KnotProt, http://knotprot.cent.uw.edu.pl/, collects information about protein structures with open polypeptide chains forming knots or slipknots. The knotting complexity of the cataloged proteins is presented in the form of a matrix diagram that shows users the knot type of the entire polypeptide chain and of each of its subchains. The pattern visible in the matrix gives the knotting fingerprint of a given protein and permits users to determine, for example, the minimal length of the knotted regions (knot's core size) or the depth of a knot, i.e. how many amino acids can be removed from either end of the cataloged protein structure before converting it from a knot to a different type of knot. In addition, the database presents extensive information about the biological functions, families and fold types of proteins with non-trivial knotting. As an additional feature, the KnotProt database enables users to submit protein or polymer chains and generate their knotting fingerprints.

Impaired fast-spiking, suppressed cortical inhibition and increased susceptibility to seizures in mice lacking Kv3.2 K+ channel proteins

Voltage-gated K+ channels of the Kv3 subfamily have unusual electrophysiological properties, including activation at very depolarized voltages (positive to −10 mV) and very fast deactivation rates, suggesting special roles in neuronal excitability. In the brain, Kv3 channels are prominently expressed in select neuronal populations, which include fast-spiking (FS) GABAergic interneurons of the neocortex, hippocampus, and caudate, as well as other high-frequency firing neurons. Although evidence points to a key role in high-frequency firing, a definitive understanding of the function of these channels has been hampered by a lack of selective pharmacological tools. We therefore generated mouse lines in which one of the Kv3 genes, Kv3.2, was disrupted by gene-targeting methods. Whole-cell electrophysiological recording showed that the ability to fire spikes at high frequencies was impaired in immunocytochemically identified FS interneurons of deep cortical layers (5-6) in which Kv3.2 proteins are normally prominent. No such impairment was found for FS neurons of superficial layers (2-4) in which Kv3.2 proteins are normally only weakly expressed. These data directly support the hypothesis that Kv3 channels are necessary for high-frequency firing. Moreover, we found that Kv3.2 −/− mice showed specific alterations in their cortical EEG patterns and an increased susceptibility to epileptic seizures consistent with an impairment of cortical inhibitory mechanisms. This implies that, rather than producing hyperexcitability of the inhibitory interneurons, Kv3.2 channel elimination suppresses their activity. These data suggest that normal cortical operations depend on the ability of inhibitory interneurons to generate high-frequency firing.

Hepatitis C virus proteins: from structure to function.

Great progress has been made over the past years in elucidating the structure and function of the hepatitis C virus (HCV) proteins, most of which are now actively being pursued as antiviral targets. The structural proteins, which form the viral particle, include the core protein and the envelope glycoproteins E1 and E2. The nonstructural proteins include the p7 viroporin, the NS2 protease, the NS3-4A complex harboring protease and NTPase/RNA helicase activities, the NS4B and NS5A proteins, and the NS5B RNA-dependent RNA polymerase. NS4B is a master organizer of replication complex formation while NS5A is a zinc-containing phosphoprotein involved in the regulation of HCV RNA replication versus particle production. Core to NS2 make up the assembly module while NS3 to NS5B represent the replication module (replicase). However, HCV proteins exert multiple functions during the viral life cycle, and these may be governed by different structural conformations and/or interactions with viral and/or cellular partners. Remarkably, each viral protein is anchored to intracellular membranes via specific determinants that are essential to protein function in the cell. This review summarizes current knowledge of the structure and function of the HCV proteins and highlights recent advances in the field.

SelenoDB 1.0 : a database of selenoprotein genes, proteins and SECIS elements

Selenoproteins are a diverse group of proteinsusually misidentified and misannotated in sequencedatabases. The presence of an in-frame UGA (stop)codon in the coding sequence of selenoproteingenes precludes their identification and correctannotation. The in-frame UGA codons are recodedto cotranslationally incorporate selenocysteine,a rare selenium-containing amino acid. The developmentof ad hoc experimental and, more recently,computational approaches have allowed the efficientidentification and characterization of theselenoproteomes of a growing number of species.Today, dozens of selenoprotein families have beendescribed and more are being discovered in recentlysequenced species, but the correct genomic annotationis not available for the majority of thesegenes. SelenoDB is a long-term project that aims toprovide, through the collaborative effort of experimentaland computational researchers, automaticand manually curated annotations of selenoproteingenes, proteins and SECIS elements. Version 1.0 ofthe database includes an initial set of eukaryoticgenomic annotations, with special emphasis on thehuman selenoproteome, for immediate inspectionby selenium researchers or incorporation into moregeneral databases. SelenoDB is freely available athttp://www.selenodb.org.

Phosphatidylethanolamine critically supports internalization of cell-penetrating protein C inhibitor.

Although their contribution remains unclear, lipids may facilitate noncanonical routes of protein internalization into cells such as those used by cell-penetrating proteins. We show that protein C inhibitor (PCI), a serine protease inhibitor (serpin), rapidly transverses the plasma membrane, which persists at low temperatures and enables its nuclear targeting in vitro and in vivo. Cell membrane translocation of PCI necessarily requires phosphatidylethanolamine (PE). In parallel, PCI acts as a lipid transferase for PE. The internalized serpin promotes phagocytosis of bacteria, thus suggesting a function in host defense. Membrane insertion of PCI depends on the conical shape of PE and is associated with the formation of restricted aqueous compartments within the membrane. Gain- and loss-of-function mutations indicate that the transmembrane passage of PCI requires a branched cavity between its helices H and D, which, according to docking studies, precisely accommodates PE. Our findings show that its specific shape enables cell surface PE to drive plasma membrane translocation of cell-penetrating PCI.

Deletion of the RNA-binding proteins ZFP36L1 and ZFP36L2 leads to perturbed thymic development and T lymphoblastic leukemia.

ZFP36L1 and ZFP36L2 are RNA-binding proteins (RBPs) that interact with AU-rich elements in the 3' untranslated region of mRNA, which leads to mRNA degradation and translational repression. Here we show that mice that lacked ZFP36L1 and ZFP36L2 during thymopoiesis developed a T cell acute lymphoblastic leukemia (T-ALL) dependent on the oncogenic transcription factor Notch1. Before the onset of T-ALL, thymic development was perturbed, with accumulation of cells that had passed through the beta-selection checkpoint without first expressing the T cell antigen receptor beta-chain (TCRbeta). Notch1 expression was higher in untransformed thymocytes in the absence of ZFP36L1 and ZFP36L2. Both RBPs interacted with evolutionarily conserved AU-rich elements in the 3' untranslated region of Notch1 and suppressed its expression. Our data establish a role for ZFP36L1 and ZFP36L2 during thymocyte development and in the prevention of malignant transformation.

A novel family of dehydrin-like proteins is involved in stress response in the human fungal pathogen Aspergillus fumigatus.

 During a search for genes controlling conidial dormancy in Aspergillus fumigatus, two dehydrin-like genes, DprA and DprB, were identified. The deduced proteins had repeated stretches of 23 amino acids that contained a conserved dehydrin-like protein (DPR) motif. Disrupted DprAΔ mutants were hypersensitive to oxidative stress and to phagocytic killing, whereas DprBΔ mutants were impaired in osmotic and pH stress responses. However, no effect was observed on their pathogenicity in our experimental models of invasive aspergillosis. Molecular dissection of the signaling pathways acting upstream showed that expression of DprA was dependent on the stress-activated kinase SakA and the cyclic AMP-protein kinase A (cAMP-PKA) pathways, which activate the bZIP transcription factor AtfA, while expression of DprB was dependent on the SakA mitogen-activated protein kinase (MAPK) pathway, and the zinc finger transcription factor PacC. Fluorescent protein fusions showed that both proteins were associated with peroxisomes and the cytosol. Accordingly, DprA and DprB were important for peroxisome function. Our findings reveal a novel family of stress-protective proteins in A. fumigatus and, potentially, in filamentous ascomycetes.

Stereospecific recognition of pyochelin and enantio-pyochelin by the PchR proteins in fluorescent pseudomonads.

The siderophore pyochelin of Pseudomonas aeruginosa promotes growth under iron limitation and induces the expression of its biosynthesis genes via the transcriptional AraC/XylS-type regulator PchR. Pseudomonas fluorescens strain CHA0 makes the optical antipode of pyochelin termed enantio-pyochelin, which also promotes growth and induces the expression of its biosynthesis genes when iron is scarce. Growth promotion and signalling by pyochelin and enantio-pyochelin are highly stereospecific and are known to involve the pyochelin and enantio-pyochelin outer-membrane receptors FptA and FetA, respectively. Here we show that stereospecificity in signalling is also based on the stereospecificity of the homologous PchR proteins of P. aeruginosa and P. fluorescens towards their respective siderophore effectors. We found that PchR functioned in the heterologous species only if supplied with its native ligand and that the FptA and FetA receptors enhanced the efficiency of signalling. By constructing and expressing hybrid and truncated PchR regulators we showed that the weakly conserved N-terminal domain of PchR is responsible for siderophore specificity. Thus, both uptake and transcriptional regulation confer stereospecificity to pyochelin and enantio-pyochelin biosynthesis.

Transposase and cointegrase: specialized transposition proteins of the bacterial insertion sequence IS21 and related elements.

The bacterial insertion sequence IS21 shares with many insertion sequences a two-step, reactive junction transposition pathway, for which a model is presented in this review: a reactive junction with abutted inverted repeats is first formed and subsequently integrated into the target DNA. The reactive junction occurs in IS21-IS21 tandems and IS21 minicircles. In addition, IS21 shows a unique specialization of transposition functions. By alternative translation initiation, the transposase gene codes for two products: the transposase, capable of promoting both steps of the reactive junction pathway, and the cointegrase, which only promotes the integration of reactive junctions but with higher efficiency. This review also includes a survey of the IS21 family and speculates on the possibility that other members present a similar transpositional specialization.