Gneg-mPLoc: A top-down strategy to enhance the quality of predicting subcellular localization of Gram-negative bacterial proteins

By incorporating the information of gene ontology, functional domain, and sequential evolution, a new predictor called Gneg-mPLoc was developed. It can be used to identify Gram-negative bacterial proteins among the following eight locations: (1) cytoplasm, (2) extracellular, (3) fimbrium, (4) flagellum, (5) inner membrane, (6) nucleoid, (7) outer membrane, and (8) periplasm. It can also be used to deal with the case when a query protein may simultaneously exist in more than one location. Compared with the original predictor called Gneg-PLoc, the new predictor is much more powerful and flexible. For a newly constructed stringent benchmark dataset in which none of proteins included has ≥25% pairwise sequence identity to any other in a same subset (location), the overall jackknife success rate achieved by Gneg-mPLoc was 85.5%, which was more than 14% higher than the corresponding rate by the Gneg-PLoc. As a user friendly web-server, Gneg-mPLoc is freely accessible at http://www.csbio.sjtu.edu.cn/bioinf/Gneg-multi/.     

A top-down approach to enhance the power of predicting human protein subcellular localization: Hum-mPLoc 2.0

Predicting subcellular localization of human proteins is a challenging problem, particularly when query proteins may have a multiplex character, i.e., simultaneously exist at, or move between, two or more different subcellular location sites. In a previous study, we developed a predictor called “Hum-mPLoc” to deal with the multiplex problem for the human protein system. However, Hum-mPLoc has the following shortcomings. (1) The input of accession number for a query protein is required in order to obtain a higher expected success rate by selecting to use the higher-level prediction pathway; but many proteins, such as synthetic and hypothetical proteins as well as those newly discovered proteins without being deposited into databanks yet, do not have accession numbers. (2) Neither functional domain nor sequential evolution information were taken into account in Hum-mPLoc, and hence its power may be reduced accordingly. In view of this, a top-down strategy to address these shortcomings has been implemented. The new predictor thus obtained is called Hum-mPLoc 2.0, where the accession number for input is no longer needed whatsoever. Moreover, both the functional domain information and the sequential evolution information have been fused into the predictor by an ensemble classifier. As a consequence, the prediction power has been significantly enhanced. The web server of Hum-mPLoc2.0 is freely accessible at http://www.csbio.sjtu.edu.cn/bioinf/hum-multi-2/.     

Some remarks on protein attribute prediction and pseudo amino acid composition

With the accomplishment of human genome sequencing, the number of sequence-known proteins has increased explosively. In contrast, the pace is much slower in determining their biological attributes. As a consequence, the gap between sequence-known proteins and attribute-known proteins has become increasingly large. The unbalanced situation, which has critically limited our ability to timely utilize the newly discovered proteins for basic research and drug development, has called for developing computational methods or high-throughput automated tools for fast and reliably identifying various attributes of uncharacterized proteins based on their sequence information alone. Actually, during the last two decades or so, many methods in this regard have been established in hope to bridge such a gap. In the course of developing these methods, the following things were often needed to consider: (1) benchmark dataset construction, (2) protein sample formulation, (3) operating algorithm (or engine), (4) anticipated accuracy, and (5) web-server establishment. In this review, we are to discuss each of the five procedures, with a special focus on the introduction of pseudo amino acid composition (PseAAC), its different modes and applications as well as its recent development, particularly in how to use the general formulation of PseAAC to reflect the core and essential features that are deeply hidden in complicated protein sequences.     

Predicting protein fold pattern with functional domain and sequential evolution information

The fold pattern of a protein is one level deeper than its structural classification, and hence is more challenging and complicated for prediction. Many efforts have been made in this regard, but so far all the reported success rates are still under 70%, indicating that it is extremely difficult to enhance the success rate even by 1% or 2%. To address this problem, here a novel approach is proposed that is featured by combining the functional domain information and the sequential evolution information through a fusion ensemble classifier. The predictor thus developed is called PFP-FunDSeqE. Tests were performed for identifying proteins among their 27 fold patterns. Compared with the existing predictors tested by a same stringent benchmark dataset, the new predictor can, for the first time, achieve over 70% success rate. The PFP-FunDSeqE predictor is freely available to the public as a web server at http://www.csbio.sjtu.edu.cn/bioinf/PFP-FunDSeqE/.     

Computational Investigation of Locked Nucleic Acid (LNA) Nucleotides in the Active Sites of DNA Polymerases by Molecular Docking Simulations

Aptamers constitute a potential class of therapeutic molecules typically selected from a large pool of oligonucleotides against a specific target. With a scope of developing unique shorter aptamers with very high biostability and affinity, locked nucleic acid (LNA) nucleotides have been investigated as a substrate for various polymerases. Various reports showed that some thermophilic B-family DNA polymerases, particularly KOD and Phusion DNA polymerases, accepted LNA-nucleoside 5′-triphosphates as substrates. In this study, we investigated the docking of LNA nucleotides in the active sites of RB69 and KOD DNA polymerases by molecular docking simulations. The study revealed that the incoming LNA-TTP is bound in the active site of the RB69 and KOD DNA polymerases in a manner similar to that seen in the case of dTTP, and with LNA structure, there is no other option than the locked C3′-endo conformation which in fact helps better orienting within the active site.     

The solution structure of a locked nucleic acid (LNA) hybridized to DNA

LNA (Locked Nucleic Acids) is a novel oligonucleotide analogue containing a conformationally restricted nucleotide with a 2'-O, 4'-C-methylene bridge that induces unprecedented thermal affinities when mixed with complementary single stranded DNA and RNA. We have used two-dimensional 1H NMR spectroscopy obtained at 750 and 500 MHz to determine a high resolution solution structure of an LNA oligonucleotide hybridized to the complementary DNA strand. The determination of the structure was based on a complete relaxation matrix analysis of the NOESY cross peaks followed by restrained molecular dynamics calculations. Forty final structures were generated for the duplex from A-type and B-type dsDNA starting structures. The root-mean-square deviation (RMSD) of the coordinates for the forty structures of the complex was 0.32A. The structures were analysed by use of calculated helix parameters. This showed that the values for rise and buckle in the LNA duplex is markedly different from canonical B-DNA at the modification site. A value of twist similar to A-DNA is also observed at the modification site. The overall length of the helix which is 27.3 A. The average twist over the sequence are 35.9 degrees +/- 0.3 degrees. Consequently, the modification does not cause the helix to unwind. The bis-intercalation of the thiazole orange dye TOTO to the LNA duplex was also investigated by 1H NMR spectroscopy to sense the structural change from the unmodified oligonucleotide. We observed that the bis-intercalation of TOTO is much less favourable in the 5'-CT(L)AG-3' site than in the unmodified 5'-CTAG-3' site. This was related to the change in the base stacking of the LNA duplex compared to the unmodified duplex.     

The Structure and Synthesis of Aeruginoic Acid (2-o-Hydroxy-phenylthiazole-4-carboxylic Acid)

The structure of aeruginoic acid which was isolated from the culture medium of P. aeruginosa was detremined and confirmed synthetically.     

Toxicity of platinum(II) amino acid (N,O) complexes parallels their binding to DNA as measured in a new solid phase assay involving a fluorescent HMG1 protein construct readout

 The compound [Pt(lysine)Cl₂] (Kplatin) was previously identified in a study of platinum amino acid complexes as a potential antitumor drug candidate. The DNA binding properties, high mobility group (HMG)-domain protein affinity for the platinated DNA, and cytotoxicity against HeLa cells of Kplatin and three related (N,O) chelated platinum(II) amino acid complexes, [Pt(arginine)Cl₂] (Rplatin), K[Pt(N_e-acetyllysine)Cl₂] (NacKplatin), and K[Pt(norleucine)Cl₂] (Norplatin), are reported. The four complexes have identical PtCl₂(N,O) coordination environments. A new solid phase screening methodology was devised in which platinated DNA probes are covalently attached to a nylon support and tested for their ability to bind a fluorescently labeled HMG-domain protein. The fluorescent HMG-domain protein was generated by expressing a fusion of the green fluorescent protein (GFP) with recombinant rat HMG1. Binding revealed by the solid phase method correlated well with the results of gel mobility shift and HeLa cytotoxicity assays. These results suggest that the net charge on the complex, rather than the nature of the side chain, is the most important factor underlying the DNA binding properties and toxicity of amino acid (N,O) chelated platinum complexes. This property explains why Kplatin was previously selected from the pool of platinum amino acid complexes based on the ability of its DNA adducts to bind HMG1. Received: 3 February 1999 / Accepted: 7 April 1999     

Solution structure of switch Arc,a mutant with 3(10) helices replacing a wild-type beta-ribbon

Adjacent N11L and L12N mutations in the antiparallel beta-ribbon of Arc repressor result in dramatic changes in local structure in which each beta-strand is replaced by a right-handed helix. The full solution structure of this "switch" Arc mutant shows that irregular 3(10) helices compose the new secondary structure. This structural metamorphosis conserves the number of main-chain and side-chain to main-chain hydrogen bonds and the number of fully buried core residues. Apart from a slight widening of the interhelical angle between alpha-helices A and B and changes in side-chain conformation of a few core residues in Arc, no large-scale structural adjustments in the remainder of the protein are necessary to accommodate the ribbon-to-helix change. Nevertheless, some changes in hydrogen-exchange rates are observed, even in regions that have very similar structures in the two proteins. The surface of switch Arc is packed poorly compared to wild-type, leading to approximately 1000A(2) of additional solvent-accessible surface area, and the N termini of the 3(10) helices make unfavorable head-to-head electrostatic interactions. These structural features account for the positive m value and salt dependence of the ribbon-to-helix transition in Arc-N11L, a variant that can adopt either the mutant or wild-type structures. The tertiary fold is capped in different ways in switch and wild-type Arc, showing how stepwise evolutionary transformations can arise through small changes in amino acid sequence.     

Origin of fibronectin type II (FN2) modules: structural analyses of distantly-related members of the kringle family idey the kringle domain of neurotrypsin as a potential link between FN2 domains and kringles

Analysis of complete genome sequences has made it clear that fibronectin type II (FN2) modules are present only in the vertebrate lineage, raising intriguing questions about the origin of this module type. Kringle domains display many similarities to FN2 domains; therefore it was suggested previously that they are highly divergent descendants of the same ancestral protein-fold. Since kringles are present in arthropodes, nematodes, and invertebrate chordates as well as in vertebrates, it is suggested that the FN2 domain arose in the vertebrate lineage through major structural modification of the more ancestral kringle fold. To explore this structural transition, in the present work we compare key structural features of two highly divergent kringle domains (the kringle of Caenorhabditis elegans Ror receptor tyrosine kinase and the kringle of rat neurotrypsin) with those of plasminogen kringles and FN2 domains. Our NMR conformation fingerprinting analysis indicates that characteristic (1)H-NMR markers of kringle or FN2 native folding, such as the dispersion of Trp aromatic connectivities and shifts of the Leu(46)/Thr(16) methyl signals, both decrease in the order kringles > neurotrypsin kringle > FN2 domains. These results suggest that the neurotrypsin kringle may represent an intermediate form between typical kringles and FN2 domains.     

The DNA Methyltransferase Dnmt1 Directly Interacts with the SET and RING Finger-associated (SRA) Domain of the Multifunctional Protein Uhrf1 to Facilitate Accession of the Catalytic Center to Hemi-methylated DNA

Dnmt1 is responsible for the maintenance DNA methylation during replication to propagate methylation patterns to the next generation. The replication foci targeting sequence (RFTS), which plugs the catalytic pocket, is necessary for recruitment of Dnmt1 to the replication site. In the present study we found that the DNA methylation activity of Dnmt1 was DNA length-dependent and scarcely methylated 12-bp short hemi-methylated DNA. Contrarily, the RFTS-deleted Dnmt1 and Dnmt1 mutants that destroyed the hydrogen bonds between the RFTS and catalytic domain showed significant DNA methylation activity even toward 12-bp hemi-methylated DNA. The DNA methylation activity of the RFTS-deleted Dnmt1 toward 12-bp hemi-methylated DNA was strongly inhibited on the addition of RFTS, but to a lesser extent by Dnmt1 harboring the mutations that impair the hydrogen bond formation. The SRA domain of Uhrf1, which is a prerequisite factor for maintenance methylation and selectively binds to hemi-methylated DNA, stimulated the DNA methylation activity of Dnmt1. The SRA to Dnmt1 concentration ratio was the determinant for the maximum stimulation. In addition, a mutant SRA, which had lost the DNA binding activity but was able to bind to Dnmt1, stimulated the DNA methylation activity of Dnmt1. The results indicate that the DNA methylation activity of Dnmt1 was stimulated on the direct interaction of the SRA and Dnmt1. The SRA facilitated acceptance of the 12-bp fluorocytosine-containing DNA by the catalytic center. We propose that the SRA removes the RFTS plug from the catalytic pocket to facilitate DNA acceptance by the catalytic center.     

The crystal structure of death receptor 6 (DR6): a potential receptor of the amyloid precursor protein (APP)

Death receptors belong to the tumor necrosis factor receptor (TNFR) super family and are intimately involved in the signal transduction during apoptosis, stress response and cellular survival. Here we present the crystal structure of recombinantly expressed death receptor six (DR6), one family member that was recently shown to bind to the amyloid precursor protein (APP) and hence to be probably involved in the development of Alzheimer's disease. The extracellular cysteine rich region of DR6, the typical ligand binding region of all TNFRs, was refined to 2.2 Å resolution and shows that its four constituting cysteine rich domains (CRDs) are arranged in a rod-like overall structure, which presents DR6-specific surface patches responsible for the exclusive recognition of its ligand(s). Based on the structural data, the general ligand binding modes of TNFRs and molecular modeling experiments we were able to elucidate structural features of the potential DR6-APP signaling complex.     

Solution Structure of S100A1 Bound to the CapZ Peptide (TRTK12)

As is typical for S100-target protein interactions, a Ca²⁺-dependent conformational change in S100A1 is required to bind to a 12-residue peptide (TRTK12) derived from the actin-capping protein CapZ. In addition, the Ca²⁺-binding affinity of S100A1 is found to be tightened (greater than threefold) when TRTK12 is bound. To examine the biophysical basis for these observations, we determined the solution NMR structure of TRTK12 in a complex with Ca²⁺-loaded S100A1. When bound to S100A1, TRTK12 forms an amphipathic helix (residues N6 to S12) with several favorable hydrophobic interactions observed between W7, I10, and L11 of the peptide and a well-defined hydrophobic binding pocket in S100A1 that is only present in the Ca²⁺-bound state. Next, the structure of S100A1-TRTK12 was compared to that of another S100A1-target complex (i.e., S100A1-RyRP12), which illustrated how the binding pocket in Ca²⁺-S100A1 can accommodate peptide targets with varying amino acid sequences. Similarities and differences were observed when the structures of S100A1-TRTK12 and S100B-TRTK12 were compared, providing insights regarding how more than one S100 protein can interact with the same peptide target. Such comparisons, including those with other S100-target and S100-drug complexes, provide the basis for designing novel small-molecule inhibitors that could be specific for blocking one or more S100-target protein interactions.     

Use of a fusion protein to obtain crystals suitable for X-ray analysis: crystallization of a GST-fused protein containing the DNA-binding domain of DNA replication-related element-binding factor, DREF.

Crystals of glutathione-S-transferase (GST)-fused protein containing the DNA-binding domain of DNA replication-related element-binding factor, DREF, were obtained under crystallization conditions similar to those for GST. Preliminary X-ray crystallographic analysis revealed that crystals of the GST-fused protein belong to space group P6(1)22 or P6(5)22 with unit cell dimensions a = b = 140.4 A, c = 93.5 A and gamma = 120 degrees, having one molecule in the crystallographic asymmetric unit. The crystals diffract to 2.5 A resolution. The cell dimensions are related to those of GST crystals thus far reported. Crystallization of the DNA-binding domain that was cleaved from the fused protein by thrombin was also carried out using several methods under numerous conditions, but efforts to produce well-ordered large crystals were unsuccessful. A possible application of GST-fusion proteins for small target proteins or domains to obtain crystals suitable for X-ray structure determination is proposed.     

Use of a renal tubule cell line (LLC-PK1) to study the nephrotoxic potential of a kappa-type bence-jones protein

Summary The cytotoxicity of a Bence-Jones protein was assessed using a porcine renal tubule cell line (LLC-PK1), with the aim of developing a model for studying the potential nephrotoxicity of these proteins. The effects of a kappa Bence-Jones protein on cell viability were studied by means of biochemical methods (supravital dye uptake and measurement of cellular enzyme activities) and morphological electron microscopy. After a 24-h-treatment with Bence-Jones protein, a moderate cytotoxicity (about 15%) was noted but only a minor difference compared to treatment with bovine albumin in the same conditions. The morphological study showed a few cells in the process of lysis, but their numbers were insufficient for the demonstration of a clear cytotoxic effect. Immunocytochemical studies showed Bence-Jones protein fixation on some cells, especially on the outer membrane. Labeling of the hyaloplasm and basal pole of a few cells pointed to internalization of protein by LLC-PK1 cells. Although the cytotoxicity of the Bence-Jones protein tested here was only moderate, the use of this model enabled its cytotoxic effect to be distinguished from that ofβ-lactoglobulin. This isolate could serve as a “moderate control” for a later study with a BJP having caused acute renal failure.     

Solution structure of ASPP2 N-terminal domain (N-ASPP2) reveals a ubiquitin-like fold

Proteins of the ASPP family bind to p53 and regulate p53-mediated apoptosis. Two family members, ASPP1 and ASPP2, have pro-apoptotic functions while iASPP shows anti-apoptotic responses. However, both the mechanism of enhancement/repression of apoptosis and the molecular basis for their different responses remain unknown. To address the role of the N-termini of pro-apoptotic ASPP proteins, we solved the solution structure of N-ASPP2 (1-83) by NMR spectroscopy. The structure of this domain reveals a beta-Grasp ubiquitin-like fold. Our findings suggest a possible role for the N-termini of ASPP proteins in binding to other proteins in the apoptotic response network and thus mediating their selective pro-apoptotic function.     

Use of ribosomal DNA sequence data to characterize and detect a neogregarine pathogen of Solenopsis invicta (Hymenoptera: Formicidae)

A neogregarine parasite of the red imported fire ant, Solenopsis invicta, was discovered recently in Florida and tentatively placed in the Mattesia genus based on morphological characterization. S. invicta infected with this Mattesia species exhibited a characteristic yellowing of the cuticle which was designated Mattesia "yellow-head disease" (YHD). The 18S rRNA gene sequence from Mattesia YHD was elucidated and compared with the neogregarine pathogens, Mattesia geminata and Ophriocystis elektroscirrha. The sequence data support the previous conclusion that Mattesia YHD is a new species that infects S. invicta. Furthermore, high sequence identity between Mattesia YHD, M. geminata (95.7%), and O. elektroscirrha (86.2%) correctly place the YHD organism in the Mattesia genus and Neogregarinorida order. Oligonucleotide primer pairs were designed to unique areas of the 18S rRNA genes of Mattesia YHD and S. invicta. Multiplex PCR resulted in sensitive and specific detection of Mattesia YHD infection of S. invicta.