Structure of Lmaj006129AAA, a hypothetical protein from Leishmania major

The N-terminal domain of Staphylothermus marinus McrB shares structural homology with PUA-like RNA binding proteins

McrBC is a conserved modification-dependent restriction system that in Escherichia coli specifically targets foreign DNA containing methylated cytosines. Crystallographic data show that the N-terminal domain of Escherichia coli McrB binds substrates via a base flipping mechanism. This region is poorly conserved among the plethora of McrB homologs, suggesting that other species may use alternative binding strategies and/or recognize different targets. Here we present the crystal structure of the N-terminal domain from Stayphlothermus marinus McrB (Sm3-180) at 1.92 Å, which adopts a PUA-like EVE fold that is closely related to the YTH and ASCH RNA binding domains. Unlike most PUA-like domains, Sm3-180 binds DNA and can associate with different modified substrates. We find the canonical 'aromatic cage' binding pocket that confers specificity for methylated bases in other EVE/YTH domains is degenerate and occluded in Sm3-180, which may contribute to its promiscuity in target recognition. Further structural comparison between different PUA-like domains identifies motifs and conformational variations that correlate with the preference for binding either DNA or RNA. Together these data have important implications for PUA-like domain specificity and suggest a broader biological versatility for the McrBC family than previously described.

Molecular basis for the recognition of methylated adenines in RNA by the eukaryotic YTH domain

Methylation of the N6 position of selected internal adenines (m(6)A) in mRNAs and noncoding RNAs is widespread in eukaryotes, and the YTH domain in a collection of proteins recognizes this modification. We report the crystal structure of the splicing factor YT521-B homology (YTH) domain of Zygosaccharomyces rouxii MRB1 in complex with a heptaribonucleotide with an m(6)A residue in the center. The m(6)A modification is recognized by an aromatic cage, being sandwiched between a Trp and Tyr residue and with the methyl group pointed toward another Trp residue. Mutations of YTH domain residues in the RNA binding site can abolish the formation of the complex, confirming the structural observations. These residues are conserved in the human YTH proteins that also bind m(6)A RNA, suggesting a conserved mode of recognition. Overall, our structural and biochemical studies have defined the molecular basis for how the YTH domain functions as a reader of methylated adenines.

Crystal structures of three protozoan homologs of tryptophanyl-tRNA synthetase

Tryptophanyl-tRNA synthetase (TrpRS) is an essential enzyme that is recognizably conserved across all forms of life. It is responsible for activating and attaching tryptophan to a cognate tRNA(Trp) molecule for use in protein synthesis. In some eukaryotes this original core function has been supplemented or modified through the addition of extra domains or the expression of variant TrpRS isoforms. The three TrpRS structures from pathogenic protozoa described here represent three illustrations of this malleability in eukaryotes. The Cryptosporidium parvum genome contains a single TrpRS gene, which codes for an N-terminal domain of uncertain function in addition to the conserved core TrpRS domains. Sequence analysis indicates that this extra domain, conserved among several apicomplexans, is related to the editing domain of some AlaRS and ThrRS. The C. parvum enzyme remains fully active in charging tRNA(Trp) after truncation of this extra domain. The crystal structure of the active, truncated enzyme is presented here at 2.4Å resolution. The Trypanosoma brucei genome contains separate cytosolic and mitochondrial isoforms of TrpRS that have diverged in their respective tRNA recognition domains. The crystal structure of the T. brucei cytosolic isoform is presented here at 2.8Å resolution. The Entamoeba histolytica genome contains three sequences that appear to be TrpRS homologs. However one of these, whose structure is presented here at 3.0Å resolution, has lost the active site motifs characteristic of the Class I aminoacyl-tRNA synthetase catalytic domain while retaining the conserved features of a fully formed tRNA(Trp) recognition domain. The biological function of this variant E. histolytica TrpRS remains unknown, but, on the basis of a completely conserved tRNA recognition region and evidence for ATP but not tryptophan binding, it is tempting to speculate that it may perform an editing function. Together with a previously reported structure of an unusual TrpRS from Giardia, these protozoan structures broaden our perspective on the extent of structural variation found in eukaryotic TrpRS homologs.

The structure of tryptophanyl-tRNA synthetase from Giardia lamblia reveals divergence from eukaryotic homologs

The 2.1A crystal structure of tryptophanyl-tRNA synthetase (TrpRS) from the diplomonad Giardia lamblia reveals that the N-terminus of this class I aminoacyl-tRNA synthetase forms a 16-residue alpha-helix. This helix replaces a beta-hairpin that is required by human TrpRS for normal activity and has been inferred to play a similar role in all eukaryotic TrpRS. The primary sequences of TrpRS homologs from several basal eukaryotes including Giardia lack a set of three residues observed to stabilize interactions with this beta-hairpin in the human TrpRS. Thus the present structure suggests that the activation reaction mechanism of TrpRS from the basal eukaryote G. lamblia differs from that of higher eukaryotes. Furthermore, the protein as observed in the crystal forms an (alpha(2))(2) homotetramer. The canonical dimer interface observed in all previous structures of tryptophanyl-tRNA synthetases is maintained, but in addition each N-terminal alpha-helix reciprocally interlocks with the equivalent helix from a second dimer to form a dimer of dimers. Although we have no evidence for tetramer formation in vivo, modeling indicates that the crystallographically observed tetrameric structure would be compatible with the tRNA binding mode used by dimeric TrpRS and TyrRS.

Prediction of protein crystallization outcome using a hybrid method

The great power of protein crystallography to reveal biological structure is often limited by the tremendous effort required to produce suitable crystals. A hybrid crystal growth predictive model is presented that combines both experimental and sequence-derived data from target proteins, including novel variables derived from physico-chemical characterization such as R(30), the ratio between a protein's DSF intensity at 30°C and at T(m). This hybrid model is shown to be more powerful than sequence-based prediction alone - and more likely to be useful for prioritizing and directing the efforts of structural genomics and individual structural biology laboratories.

The YTH domain is a novel RNA binding domain

The YTH (YT521-B homology) domain was identified by sequence comparison and is found in 174 different proteins expressed in eukaryotes. It is characterized by 14 invariant residues within an alpha-helix/beta-sheet structure. Here we show that the YTH domain is a novel RNA binding domain that binds to a short, degenerated, single-stranded RNA sequence motif. The presence of the binding motif in alternative exons is necessary for YT521-B to directly influence splice site selection in vivo. Array analyses demonstrate that YT521-B predominantly regulates vertebrate-specific exons. An NMR titration experiment identified the binding surface for single-stranded RNA on the YTH domain. Structural analyses indicate that the YTH domain is related to the pseudouridine synthase and archaeosine transglycosylase (PUA) domain. Our data show that the YTH domain conveys RNA binding ability to a new class of proteins that are found in all eukaryotic organisms.

Membrane protein SMP-1 is required for normal flagellum function in Leishmania

Eukaryotic flagella and cilia are surrounded by a membrane that is continuous with, but distinct from, the rest of the plasma membrane. In Leishmania parasites, the inner leaflet of the flagellar membrane is coated with the acylated membrane protein, SMP-1. Here, we provide evidence that SMP-1 stabilizes the flagellar membrane and is required for flagella elongation and function. The expression and flagella targeting of SMP-1 is tightly associated with flagella elongation during amastigote to promastigote differentiation. Deletion of the genes encoding SMP-1 and the flagellar pocket protein SMP-2, led to the production of short flagella and defects in motility. Alterations in the physical properties of the smp-1/smp-2(-/-) flagellar membrane were suggested by: (1) the accumulation of membrane vesicles in the flagellar matrix, and (2) further retraction of flagella following partial inhibition of sterol and sphingolipid biosynthesis. The flagella phenotype of the smp-1/smp-2(-/-) null mutant was reversed by re-expression of SMP-1, but not SMP-2. SMP-1 contains a jelly-roll beta-sheet structure that is probably conserved in all SMP proteins, and forms stable homo-oligomers in vivo. We propose that the SMP-1 coat generates and/or stabilizes sterol- and sphingolipid-rich domains in the flagellar membrane.

Prediction of enzyme classes from 3D structure: a general model and examples of experimental-theoretic scoring of peptide mass fingerprints of Leishmania proteins

The number of protein and peptide structures included in Protein Data Bank (PDB) and Gen Bank without functional annotation has increased. Consequently, there is a high demand for theoretical models to predict these functions. Here, we trained and validated, with an external set, a Markov Chain Model (MCM) that classifies proteins by their possible mechanism of action according to Enzyme Classification (EC) number. The methodology proposed is essentially new, and enables prediction of all EC classes with a single equation without the need for an equation for each class or nonlinear models with multiple outputs. In addition, the model may be used to predict whether one peptide presents a positive or negative contribution of the activity of the same EC class. The model predicts the first EC number for 106 out of 151 (70.2%) oxidoreductases, 178/178 (100%) transferases, 223/223 (100%) hydrolases, 64/85 (75.3%) lyases, 74/74 (100%) isomerases, and 100/100 (100%) ligases, as well as 745/811 (91.9%) nonenzymes. It is important to underline that this method may help us predict new enzyme proteins or select peptide candidates that improve enzyme activity, which may be of interest for the prediction of new drugs or drug targets. To illustrate the model's application, we report the 2D-Electrophoresis (2DE) isolation from Leishmania infantum as well as MADLI TOF Mass Spectra characterization and theoretical study of the Peptide Mass Fingerprints (PMFs) of a new protein sequence. The theoretical study focused on MASCOT, BLAST alignment, and alignment-free QSAR prediction of the contribution of 29 peptides found in the PMF of the new protein to specific enzyme action. This combined strategy may be used to identify and predict peptides of prokaryote and eukaryote parasites and their hosts as well as other superior organisms, which may be of interest in drug development or target identification.

Structural genomics reveals EVE as a new ASCH/PUA-related domain

We report on several proteins recently solved by structural genomics consortia, in particular by the Northeast Structural Genomics consortium (NESG). The proteins considered in this study differ substantially in their sequences but they share a similar structural core, characterized by a pseudobarrel five-stranded beta sheet. This core corresponds to the PUA domain-like architecture in the SCOP database. By connecting sequence information with structural knowledge, we characterize a new subgroup of these proteins that we propose to be distinctly different from previously described PUA domain-like domains such as PUA proper or ASCH. We refer to these newly defined domains as EVE. Although EVE may have retained the ability of PUA domains to bind RNA, the available experimental and computational data suggests that both the details of its molecular function and its cellular function differ from those of other PUA domain-like domains. This study of EVE and its relatives illustrates how the combination of structure and genomics creates new insights by connecting a cornucopia of structures that map to the same evolutionary potential. Primary sequence information alone would have not been sufficient to reveal these evolutionary links.

Study of peptide fingerprints of parasite proteins and drug-DNA interactions with Markov-Mean-Energy invariants of biopolymer molecular-dynamic lattice networks

Since the advent of Molecular Dynamics (MD) in biopolymers science with the study by Karplus et al. on protein dynamics, MD has become the by foremost well established, computational technique to investigate structure and function of biomolecules and their respective complexes and interactions. The analysis of the MD trajectories (MDTs) remains, however, the greatest challenge and requires a great deal of insight, experience, and effort. Here, we introduce a new class of invariants for MDTs based on the spatial distribution of Mean-Energy values ξ_k (L) on a 2D Euclidean space representation of the MDTs. The procedure forces one MD trajectory to fold into a 2D Cartesian coordinates system using a step-by-step procedure driven by simple rules. The ξ_k (L) values are invariants of a Markov matrix (¹ Π), which describes the probabilities of transition between two states in the new 2D space; which is associated to a graph representation of MDTs similar to the lattice networks (LNs) of DNA and protein sequences. We also introduce a new algorithm to perform phylogenetic analysis of peptides based on MDTs instead of the sequence of the polypeptide. In a first experiment, we illustrate this algorithm for 35 peptides present on the Peptide Mass Fingerprint (PMF) of a new protein of Leishmania infantum studied in this work. We report, by the first time, 2D Electrophoresis isolation, MALDI TOF Mass Spectroscopy characterization, and MASCOT search results for this PMF. In a second experiment, we construct the LNs for 422 MDTs obtained in DNA-Drug Docking simulations of the interaction of 57 anticancer furocoumarins with a DNA oligonucleotide. We calculated the respective ξ_k (L) values for all these LNs and used them as inputs to train a new classifier with Accuracy = 85.44% and 84.91% in training and validation respectively. The new model can be used as scoring function to guide DNA-Drug Docking studies in drug design of new coumarins for PUVA therapy. The new phylogenetics analysis algorithms encode information different from sequence similarity and may be used to analyze MDTs obtained in Docking or modeling experiments for any classes of biopolymers. The work opens new perspective on the analysis and applications of MD in polymer sciences.

QSAR for RNases and theoretic-experimental study of molecular diversity on peptide mass fingerprints of a new Leishmania infantum protein

The toxicity and low success of current treatments for Leishmaniosis determines the search of new peptide drugs and/or molecular targets in Leishmania pathogen species (L. infantum and L. major). For example, Ribonucleases (RNases) are enzymes relevant to several biologic processes; then, theoretical and experimental study of the molecular diversity of Peptide Mass Fingerprints (PMFs) of RNases is useful for drug design. This study introduces a methodology that combines QSAR models, 2D-Electrophoresis (2D-E), MALDI-TOF Mass Spectroscopy (MS), BLAST alignment, and Molecular Dynamics (MD) to explore PMFs of RNases. We illustrate this approach by investigating for the first time the PMFs of a new protein of L. infantum. Here we report and compare new versus old predictive models for RNases based on Topological Indices (TIs) of Markov Pseudo-Folding Lattices. These group of indices called Pseudo-folding Lattice 2D-TIs include: Spectral moments pi ( k )(x,y), Mean Electrostatic potentials xi ( k )(x,y), and Entropy measures theta ( k )(x,y). The accuracy of the models (training/cross-validation) was as follows: xi ( k )(x,y)-model (96.0%/91.7%)>pi ( k )(x,y)-model (84.7/83.3) > theta ( k )(x,y)-model (66.0/66.7). We also carried out a 2D-E analysis of biological samples of L. infantum promastigotes focusing on a 2D-E gel spot of one unknown protein with M<20, 100 and pI <7. MASCOT search identified 20 proteins with Mowse score >30, but not one >52 (threshold value), the higher value of 42 was for a probable DNA-directed RNA polymerase. However, we determined experimentally the sequence of more than 140 peptides. We used QSAR models to predict RNase scores for these peptides and BLAST alignment to confirm some results. We also calculated 3D-folding TIs based on MD experiments and compared 2D versus 3D-TIs on molecular phylogenetic analysis of the molecular diversity of these peptides. This combined strategy may be of interest in drug development or target identification.

HP-Lattice QSAR for dynein proteins: experimental proteomics (2D-electrophoresis, mass spectrometry) and theoretic study of a Leishmania infantum sequence

The toxicity and inefficacy of actual organic drugs against Leishmaniosis justify research projects to find new molecular targets in Leishmania species including Leishmania infantum (L. infantum) and Leishmaniamajor (L. major), both important pathogens. In this sense, quantitative structure-activity relationship (QSAR) methods, which are very useful in Bioorganic and Medicinal Chemistry to discover small-sized drugs, may help to identify not only new drugs but also new drug targets, if we apply them to proteins. Dyneins are important proteins of these parasites governing fundamental processes such as cilia and flagella motion, nuclear migration, organization of the mitotic splinde, and chromosome separation during mitosis. However, despite the interest for them as potential drug targets, so far there has been no report whatsoever on dyneins with QSAR techniques. To the best of our knowledge, we report here the first QSAR for dynein proteins. We used as input the Spectral Moments of a Markov matrix associated to the HP-Lattice Network of the protein sequence. The data contain 411 protein sequences of different species selected by ClustalX to develop a QSAR that correctly discriminates on average between 92.75% and 92.51% of dyneins and other proteins in four different train and cross-validation datasets. We also report a combined experimental and theoretic study of a new dynein sequence in order to illustrate the utility of the model to search for potential drug targets with a practical example. First, we carried out a 2D-electrophoresis analysis of L. infantum biological samples. Next, we excised from 2D-E gels one spot of interest belonging to an unknown protein or protein fragment in the region M<20,200 and pI<4. We used MASCOT search engine to find proteins in the L. major data base with the highest similarity score to the MS of the protein isolated from L. infantum. We used the QSAR model to predict the new sequence as dynein with probability of 99.99% without relying upon alignment. In order to confirm the previous function annotation we predicted the sequences as dynein with BLAST and the omniBLAST tools (96% alignment similarity to dyneins of other species). Using this combined strategy, we have successfully identified L. infantum protein containing dynein heavy chain, and illustrated the potential use of the QSAR model as a complement to alignment tools.

Assessment of predictions submitted for the CASP7 domain prediction category

This paper details the assessment process and evaluation results for the Critical Assessment of Protein Structure Prediction (CASP7) domain prediction category. Domain predictions were assessed using the Normalized Domain Overlap score introduced in CASP6 and the accuracy of prediction of domain break points. The results of the analysis clearly demonstrate that the best methods are able to make consistently reliable predictions when the target has a structural template, although they are less good when the domain break occurs in a region not covered by a template. The conditions of the experiment meant that it was impossible to draw any conclusions about domain prediction for free modeling targets and it was also difficult to draw many distinctions between the best groups. Two thirds of the targets submitted were single domains and hence regarded as easy to predict. Even those targets defined as having multiple domains always had at least one domain with a similar template structure.

The structure of Plasmodium vivax phosphatidylethanolamine-binding protein suggests a functional motif containing a left-handed helix

The structure of a putative Raf kinase inhibitor protein (RKIP) homolog from the eukaryotic parasite Plasmodium vivax has been studied to a resolution of 1.3 A using multiple-wavelength anomalous diffraction at the Se K edge. This protozoan protein is topologically similar to previously studied members of the phosphatidylethanolamine-binding protein (PEBP) sequence family, but exhibits a distinctive left-handed alpha-helical region at one side of the canonical phospholipid-binding site. Re-examination of previously determined PEBP structures suggests that the P. vivax protein and yeast carboxypeptidase Y inhibitor may represent a structurally distinct subfamily of the diverse PEBP-sequence family.