dSPRINT: predicting DNA,RNA, ion,peptide and small molecule interaction sites within protein domains

Domains are instrumental in facilitating protein interactions with DNA, RNA, small molecules, ions and peptides. Identifying ligand-binding domains within sequences is a critical step in protein function annotation, and the ligand-binding properties of proteins are frequently analyzed based upon whether they contain one of these domains. To date, however, knowledge of whether and how protein domains interact with ligands has been limited to domains that have been observed in co-crystal structures; this leaves approximately two-thirds of human protein domain families uncharacterized with respect to whether and how they bind DNA, RNA, small molecules, ions and peptides. To fill this gap, we introduce dSPRINT, a novel ensemble machine learning method for predicting whether a domain binds DNA, RNA, small molecules, ions or peptides, along with the positions within it that participate in these types of interactions. In stringent cross-validation testing, we demonstrate that dSPRINT has an excellent performance in uncovering ligand-binding positions and domains. We also apply dSPRINT to newly characterize the molecular functions of domains of unknown function. dSPRINT’s predictions can be transferred from domains to sequences, enabling predictions about the ligand-binding properties of 95% of human genes. The dSPRINT framework and its predictions for 6503 human protein domains are freely available at http://protdomain.princeton.edu/dsprint.     

The low expression of Dmrt7 is associated with spermatogenic arrest in cattle-yak

Dmrt7 is a member of the DM domain family of genes. Dmrt7 deficiency is also a strong candidate as a cause for male cattle-yak infertility, as it is regarded as essential for male spermatogenesis, between the pachynema and diplonema stages. In our study, the coding region sequence of yak and cattle-yak Dmrt7 was cloned by molecular cloning techniques, and the sequence, conserved domains, functional sites, and secondary and tertiary structures of the Dmrt7-encoded protein were predicted and analyzed using bioinformatics methods. The coding region sequences of the Dmrt7 gene, encoding 370 amino acids, were consistent in yak and cattle-yak. The protein encoded by yak and cattle-yak Dmrt7 contains a DM domain. We detected Dmrt7 mRNA expression in testis, but not in any other tissue. Dmrt7 mRNA and protein expression was significantly higher in testis of cattle and yak than that in cattle-yak (p < 0.01). Histological analysis indicated that seminiferous tubules in male cattle-yak were highly vacuolated and contained primarily Sertoli cells and spermatogonia, while those of cattle and yak contained abundant primary spermatocytes. Male cattle-yak testis contained a significantly larger number of apoptotic cells than those in cattle and yak assessed by terminal deoxynucleotidyl transferase dUTP nick end-labeling (TUNEL) analysis. The accumulation of SCP3-positive spermatocytes indicated the arrest of spermatogenesis at the pachynema stage in the cattle-yak. These results suggest low levels of Dmrt7 expression lead to male sterility in cattle-yak. The molecular function of Dmrt7 and the regulation of its expression warrant need to be examined in future studies.     

Membrane anchors effectively traffic recombinant human glucocerebrosidase to the protein storage vacuole of Arabidopsis seeds but do not adequately control N-glycan maturation

Key message

Human glucocerebrosidase with vacuolar anchoring domains was targeted to protein storage vacuoles (PSVs) of Arabidopsis seeds, but unexpectedly via the Golgi complex. PSV-targeting to effectively avoid problematic N-glycans is protein dependent.

Abstract

Plant-specific N-glycosylation patterns elaborated within the Golgi complex are a major limitation of using plants to produce biopharmaceuticals as the presence of β1,2 xylose and/or α1,3 fucose residues on the recombinant glycoprotein can render the product immunogenic if administrated parenterally. A reporter protein fused to a vacuolar membrane targeting motif comprised of the BP-80 transmembrane domain (TMD), and the cytoplasmic tail (CT) of α-tonoplast intrinsic protein (α-TIP) is delivered to protein storage vacuoles (PSVs) of tobacco seeds by ER-derived transport vesicles that bypass the Golgi complex. This prompted us to investigate whether a pharmaceutical glycoprotein is targeted to PSVs using the same targeting sequences, thus avoiding the unwanted plant-Golgi-specific complex N-glycan modifications. The human lysosomal acid β-glucosidase (glucocerebrosidase; GCase) (EC 3.2.1.45) fused to the BP-80 TMD and α-TIP CT was produced in Arabidopsis thaliana wild-type (Col-0) seeds. The chimeric GCase became localized in PSVs but transited through the Golgi complex, as indicated by biochemical analyses of the recombinant protein’s N-glycans. Our findings suggest that use of this PSV-targeting strategy to avoid problematic N-glycan maturation on recombinant therapeutic proteins is not consistently effective, as it is likely protein- and/or species-specific.     

Creation of a HEK293 cell line stably expressing the proteasome subunit PSMD14 fused with fluorescent protein EGFP and HTBT tag

In the present work, a stable cell line based on human embryonic kidney HEK293 cells, which expresses the proteasome subunit PSMD14 fused to the fluorescent protein EGFP and HTBH tag, has been selected. This protein, PSMD14-EGFP-HTBH, was shown to be completely incorporated into proteasomes, and such complexes have a chymotrypsin-like peptidase activity that is fundamental for proteasomes. The constructed and selected cell line can be used for further fluorescent studies of localization of proteasomes in the cell.     

Role of α-helical domains in functioning of ATP-dependent Lon protease of Escherichia coli

Homooligomeric ATP-dependent LonA proteases are bifunctional enzymes belonging to the superfamily of AAA⁺ proteins. Their subunits are formed by five successively connected domains, i.e., N-terminal (N), α-helical (HI(CC)), nucleotide-binding (NB), the second α-helical (H), and proteolytic (P) domains. The presence of the inserted HI(CC) domain determines the uniqueness of LonA proteases among the AAA⁺ proteins. The role of the α-helical domains in the LonA protease functioning was studied with an example of E. coli Lon protease (Ec-Lon). The properties of the intact Ec-Lon and its mutant forms, i.e., Lon-R164A and Lon-R542A bearing the substituted arginine residues at the similar positions in the HI(CC) and H domains, were compared. The H domain was shown to play a crucial role in ATP hydrolysis and enzyme binding to the target protein. The HI(CC) domain is not decisive for the manifestation of the catalytic properties of the enzyme. However, it affects the functioning of Lon ATPase and peptidase sites and is involved in maintaining enzyme stability. The participation of the HI(CC) domain in the formation of three-dimensional structures of LonA proteases and/or their complexes with DNA is suggested.     

Role of the virion host shutoff protein in neurovirulence of monkey B virus (Macacine herpesvirus 1)

Monkey B virus(Macacine herpesvirus 1; BV) is noted for its extreme neurovirulence in humans. Since the vhs protein encoded by the UL41 gene has been shown to be a neurovirulence factor in the related human herpes simplex viruses, the role of the UL41 gene in BV neurovirulence was investigated. BV mutants were constructed that lacked the entire UL41 ORF(Δ41) or had the RNase active site mutated(Δ41A). Neither mutant shut off host protein synthesis, degraded β-actin mRNA, or prevented an IFN-β response, indicating that the vhs protein and its RNase activity are both necessary for these activities. Replication of both mutants in primary mouse cells was impaired and they exhibited a prolonged disease course in mice. Whereas Δ41 infected mice were euthanized for symptoms related to central nervous system(CNS) infection, Δ41A infected mice were euthanized primarily for symptoms of autonomic nervous system dysfunction. While neuroinvasiveness was not affected, lesions in the CNS were more limited in size, anatomical distribution, and severity than for wild-type virus. These results indicate that the vhs protein affects the general replicative efficiency of BV in vivo rather than being a specific neurovirulence factor critical for invasion of or preferential replication in the CNS.     

Conformational effects of the substitution of ARG for GLY 13 in theras oncogene-encoded P21 protein

The effect of the substitution of Arg for Gly 13 on the structure of the transforming region decapeptide (Leu 6-Gly 15) of the ras oncogene encoded P21 protein has been investigated using conformational energy analysis. A human malignancy has been identified that contains a ras gene with a single mutation in the thirteenth codon such that the encoded protein would have Arg substituted for Gly at this position, and transfection of cells in culture with this gene results in malignant transformation. Conformational analysis demonstrates that the Arg 13 decapeptide adopts a conformation identical to that for other peptides with substitutions at position 13 (Asp 13, Val 13) from transforming proteins that is distinctively different from that for peptides (Gly 13, Ser 13) from normal, nontransforming proteins. This is found to be an indirect effect resulting from changes in the conformation of Gly 12 produced by substitutions at position 13. These results are consistent with recent analysis of crystallographic data of proteins on conformational preferences for glycine in tripeptide sequences.     

Loss of EphB6 protein expression in human colorectal cancer correlates with poor prognosis

Erythropoietin-producing hepatocyte (Eph) receptor family constitutes the largest family of tyrosine kinase receptors in the human genome. Loss of EphB6, a kinase-deficient receptor, correlated with a negative outcome in several carcinomas. This study aimed to investigate the expression of EphB6 protein and mRNA levels in colorectal cancers (CRCs) and possible correlations with clinicopathological variables and prognosis. To assess protein expression level, 124 CRCs and 57 colorectal adenomas samples were examined by immunostaining, the mRNA level of 43 paired CRC and the adjacent normal tissues were detected by using SYBR Green real-time PCR method. Decreased expression of EphB6 protein was found in CRC as compared with adenoma and normal tissues (χ² = 10.146, P = 0.001 and χ² = 45.333, P < 0.001, respectively). Low EphB6 mRNA expression was detected in 83.8 % of cancers with negative or low EphB6 protein expression. The loss of EphB6 protein in CRC was positively associated with poorly differentiation (P < 0.001), lymph node metastasis (P = 0.006), Dukes stage (P = 0.002) and depth of invasion (P = 0.016). The patients with lymph node metastasis had a worse prognosis independently of gender, age, tumor site, stage and differentiation (RR = 0.404, CI 0.267–0.213, P < 0.001). Low levels of EphB6 protein expression are associated with a shorter mean duration of survival in colorectal cancer. Our results demonstrated that EphB6 may represent a novel, useful tissue biomarker for the prediction of survival rate in CRC.     

Comparative structural analysis of lipid binding START domains

Background

Steroidogenic acute regulatory (StAR) protein related lipid transfer (START) domains are small globular modules that form a cavity where lipids and lipid hormones bind. These domains can transport ligands to facilitate lipid exchange between biological membranes, and they have been postulated to modulate the activity of other domains of the protein in response to ligand binding. More than a dozen human genes encode START domains, and several of them are implicated in a disease.

Principal Findings

We report crystal structures of the human STARD1, STARD5, STARD13 and STARD14 lipid transfer domains. These represent four of the six functional classes of START domains.

Significance

Sequence alignments based on these and previously reported crystal structures define the structural determinants of human START domains, both those related to structural framework and those involved in ligand specificity.

Enhanced version

This article can also be viewed as an enhanced version in which the text of the article is integrated with interactive 3D representations and animated transitions. Please note that a web plugin is required to access this enhanced functionality. Instructions for the installation and use of the web plugin are available in Text S1.     

Multiple correlations of mRNA expression and protein abundance in human cytokine profile

With the development of genomic study, researchers found that it is insufficient to predict protein expression from quantitative mRNA data in large scale, which is contrary to the traditional opinion that mRNA expression correlates with protein abundance at the single gene level. To try to solve the apparent conflicting views, here we set up a series of research models and chose soluble cytokines as targets. First, human peripheral blood mononuclear cell (PBMC) from one health donor was treated with 16 continuously changing conditions, the protein and mRNA profile were analyzed by multiplex Luminex and genomic microarray, respectively. Among the tested genes, around half mRNA correlated well with their corresponding proteins (ρ > 0.8), however if we put all the genes together, the correlation coefficient for the 16 conditions varied from 0.29 to 0.71. Second, PBMC from 14 healthy donors were stimulated with the same condition and it was found that the correlation coefficient went down (ρ < 0.6). Third, 28 rheumatoid arthritis (RA) patients were tested for their response to the same external stimuli and it turned out different individual displayed different protein expression pattern as expect. Lastly, autoimmune disease cohorts (8 diseases including RA, 103 patients in total) were assayed on the whole view. It was observed that there was still some similarity in the protein profile among patients from the single disease type although completely different patterns were displayed across different disease categories. This study built a good bridge between single gene analysis and the whole genome study and may give a reasonable explanation for the two conflicting views in current biological science.     

Substitution of selenocysteine for cysteine in a reticulocyte lysate protein synthesis system

Selenocysteine occurs in the peptide backbone of several selenoenzymes. The mechanism, of selenocysteine incorporation has not been well characterized. The incorporation of selenocysteine into protein in a rabbit reticulocyte lysate (RRL) was studied at high levels of selenocysteine. [⁷⁵Se]Selenocysteine incorporation was inhibited by cycloheximide and by nuclease treatment. Random RNA copolymers were tested for protein synthesis activity in the messenger RNA-dependent RRL system. Of the active polymers, poly CIU and GU most strongly stimulated the incorporation of selenocysteine. In a series of four polymers with different ratios of U to G, incorporation of selenocysteine and cysteine increased with increasing percentages of U, suggesting that selenocysteine and cysteine responded to the same codon, presumably UGU. Of the 20 protein amino acids, only cysteine and cystine competed with selenocysteine incorporation. Selenocysteine was charged to cysteine-accepting tRNA in RRL. These results show that at supraphysiological concentrations selenocysteine can substitute for cysteine in RRL protein synthesis. Misincorporation of selenocysteine could be important when animal tissues contain high levels of selenium.     

An ice nucleation protein from Fusarium acuminatum: cloning,expression, biochemical characterization and computational modeling

Ice nucleation proteins (INP) are a major cause of frost damage in plants and crops. Here, an INP gene from Fusarium acuminatum was optimized, synthesized, expressed in E.coli and subsequently purified and characterized. The protein belongs to the second class of ice nucleation proteins with an optimum pH 5.5, relative activity and stability between pH 5 and 9.5 and up to 45 °C. The protein was fully active and stable in the presence of dimethyl sulfoxide (DMSO), dioxane, acetone and ethyl acetate. Moreover, it retained over 50 % of its original activity in the presence of polyvinyl alcohol. The 3D structure model of the INP-F indicated the protein had three distinct domains as exist in other ice nucleation proteins with some variations. Considering these promising results, INP-F could be a novel candidate for industrial applications.     

Preparation of cinnamoyl-CoA and analysis of the enzyme activity of recombinant CCR protein from Populus tomentosa

Cinnamoyl-CoA reductase (CCR, EC 1.2.1.44), which catalyzes the reduction of cinnamoyl-CoA esters to their respective cinnamaldehydes, is considered as a key enzyme in lignin formation. The substrates of CCR, cinnamoyl-CoA esters, are products of 4-Coumarate-CoA ligase (4CL, EC 6.2.1.12), which is an enzyme upstream of CCR. The PtCCR and Pt4CL were isolated from Populus tomentosa and expressed in E. coli. Results showed that 4CL can catalyze the conversion of hydroxycinnamic acids to cinnamoyl-CoA esters, with high efficiency. The purification of esters using SPE cartridges suggested that 40 % methanol with 0.1 M of acetic acid was the optimal elution buffer for cinnamoyl-CoA esters. The optimization of prokaryotic expression demonstrated that the best expression conditions for recombinant PtCCR was 6 h of 0.4 mM IPTG induction at 37 °C. PtCCR enzyme assay illustrated that the recombinant protein can catalyze the reduction of cinnamoyl-CoA esters. Kinetics analysis showed that feruloyl-CoA has higher affinity to PtCCR with faster reaction speed (Vmax), indicating that feruloyl-CoA was the most favorable substrate for PtCCR catalysis. The recombinant protein was expressed in E. coli, purified through affinity column chromatography, and characterized by SDS-PAGE. SPE cartridges were used to purify the ester products of the Pt4CL reaction. HPLC-MS was used to analyze the structure of esters and evaluate their purity or quantity. Furthermore, the enzyme activity of recombinant CCR to feruloyl-CoA at different pHs indicated that compartmentalization may be an important factor in lignin monomer formation.     

ProteinArchitect: Protein Evolution above the Sequence Level

Background

While many authors have discussed models and tools for studying protein evolution at the sequence level, molecular function is usually mediated by complex, higher order features such as independently folding domains and linear motifs that are based on or embedded in a particular arrangment of features such as secondary structure elements, transmembrane domains and regions with intrinsic disorder. This ‘protein architecture’ can, in its most simplistic representation, be visualized as domain organization cartoons that can be used to compare proteins in terms of the order of their mostly globular domains.

Methodology

Here, we describe a visual approach and a webserver for protein comparison that extend the domain organization cartoon concept. By developing an information-rich, compact visualization of different protein features above the sequence level, potentially related proteins can be compared at the level of propensities for secondary structure, transmembrane domains and intrinsic disorder, in addition to PFAM domains. A public Web server is available at www.proteinarchitect.net, while the code is provided at protarchitect.sourceforge.net.

Conclusions/Significance

Due to recent advances in sequencing technologies we are now flooded with millions of predicted proteins that await comparative analysis. In many cases, mature tools focused on revealing hits with considerable global or local similarity to well-characterized proteins will not be able to lead us to testable hypotheses about a protein''s function, or the function of a particular region. The visual comparison of different types of protein features with ProteinArchitect will be useful when assessing the relevance of similarity search hits, to discover subgroups in protein families and superfamilies, and to understand protein regions with conserved features outside globular regions. Therefore, this approach is likely to help researchers to develop testable hypotheses about a protein''s function even if is somewhat distant from the more characterized proteins, by facilitating the discovery of features that are conserved above the sequence level for comparison and further experimental investigation.     

Biochemical characterization of a metallothionein-like protein in rat brain

Recent investigations from this laboratory have identified a metallothionein-like protein in the rat brain with an elution volume (ve/vo) of 2.08 and a molecular weight of smaller than 10,000. The synthesis of this protein was stimulated following intracerebroventricular (icv, 0.20 μmol zinc/μL/h, 48 h), but not intraperitoneal (ip) administration of ZnSO₄. Furthermore, chronic ip administration of ZnSO₄ (5.0 mg/kg/d/10 d) did not alter the level of the metallothionein-like protein in the brain. However, the hepatic metallothionein was induced following icv administration of ZnSO₄. The chromatofocusing of metallothionein-like protein isolated by gel permeation chromatography on Sephadex G-75 exhibits three zinc-binding peaks, which focus on pH 6.8, 6.2, and 5.3, respectively. It is expected that the protein peak focusing at 5.3 is a metallothionein-like protein. Purification of the zinc stimulated metallothionein-like protein on ion exchange chromatography on DEAE-Sephadex A-25 columns, using a linear gradient elution procedure produced two isoforms, eluting, respectively, at 75 and 137 mM of Tris-acetate buffer, pH 7.5. The comparative high-performance liquid chromatographic (HPLC) profiles of the zinc-induced hepatic metallothionein isoforms I and II (retention times 17.39 and 18.73 min) and that of the zinc-stimulated metallothionein-like protein isoforms I and II (retention times 17.32 and 18.64 min) are very similar. The function(s) of the metallothionein-like protein isoforms in the brain remains to be elucidated.     

Binding of mitochondrial ATPase from ox heart to its naturally occurring inhibitor protein: Localization by antibody binding

The naturally occurring ATPase inhibitor protein from ox heart mitochondria was cross-linked to its binding site on the mitochondrial ATPase using 1-ethyl-3-(dimethylamino)propyl carbodiimide. The cross-linked product, when transferred electrophoretically to a nitrocellulose sheet, reacted with antibodies directed against the inhibitor protein and the β-subunit of the ATPase. It was concluded that the binding site for the inhibitor protein lies on the β-subunit.