Computational Prediction of Candidate Proteins for S-Nitrosylation in Arabidopsis thaliana

Nitric oxide (NO) is an important signaling molecule that regulates many physiological processes in plants. One of the most important regulatory mechanisms of NO is S-nitrosylation—the covalent attachment of NO to cysteine residues. Although the involvement of cysteine S-nitrosylation in the regulation of protein functions is well established, its substrate specificity remains unknown. Identification of candidates for S-nitrosylation and their target cysteine residues is fundamental for studying the molecular mechanisms and regulatory roles of S-nitrosylation in plants. Several experimental methods that are based on the biotin switch have been developed to identify target proteins for S-nitrosylation. However, these methods have their limits. Thus, computational methods are attracting considerable attention for the identification of modification sites in proteins. Using GPS-SNO version 1.0, a recently developed S-nitrosylation site-prediction program, a set of 16,610 candidate proteins for S-nitrosylation containing 31,900 S-nitrosylation sites was isolated from the entire Arabidopsis proteome using the medium threshold. In the compartments “chloroplast,” “CUL4-RING ubiquitin ligase complex,” and “membrane” more than 70% of the proteins were identified as candidates for S-nitrosylation. The high number of identified candidates in the proteome reflects the importance of redox signaling in these compartments. An analysis of the functional distribution of the predicted candidates showed that proteins involved in signaling processes exhibited the highest prediction rate. In a set of 46 proteins, where 53 putative S-nitrosylation sites were already experimentally determined, the GPS-SNO program predicted 60 S-nitrosylation sites, but only 11 overlap with the results of the experimental approach. In general, a computer-assisted method for the prediction of targets for S-nitrosylation is a very good tool; however, further development, such as including the three dimensional structure of proteins in such analyses, would improve the identification of S-nitrosylation sites.     

Identification of a Novel Proline-Rich Antimicrobial Peptide from Brassica napus

Proline-rich antimicrobial peptides (PR-AMPs) are a group of cationic host defense peptides that are characterized by a high content of proline residues. Up to now, they have been reported in some insects, vertebrate and invertebrate animals, but are not found in plants. In this study, we performed an in silico screening of antimicrobial peptides, which led to discovery of a Brassica napus gene encoding a novel PR-AMP. This gene encodes a 35-amino acid peptide with 13 proline residues, designated BnPRP1. BnPRP1 has 40.5% identity with a known proline-rich antimicrobial peptide SP-B from the pig. BnPRP1 was artificially synthetized and cloned into the prokaryotic expression vector pET30a/His-EDDIE-GFP. Recombinant BnPRP1 was produced in Escherichia coli and has a predicted molecular mass of 3.8 kDa. Analysis of its activity demonstrated that BnPRP1 exhibited strong antimicrobial activity against Gram-positive bacterium, Gram-negative bacterium, yeast and also had strong antifungal activity against several pathogenic fungi, such as Sclerotinia sclerotiorum, Mucor sp., Magnaporthe oryzae and Botrytis cinerea. Circular dichroism (CD) revealed the main secondary structure of BnPRP1 was the random coil. BnPRP1 gene expression detected by qRT-PCR is responsive to pathogen inoculation. At 48 hours after S. sclerotiorum inoculation, the expression of BnPRP1 increased significantly in the susceptible lines while slight decrease occurred in resistant lines. These suggested that BnPRP1 might play a role in the plant defense response against S. sclerotiorum. BnPRP1 isolated from B. napus was the first PR-AMP member that was characterized in plants, and its homology sequences were found in some other Brassicaceae plants by the genome sequences analysis. Compared with the known PR-AMPs, BnPRP1 has the different primary sequences and antimicrobial activity. Above all, this study gives a chance to cast a new light on further understanding about the AMPs’ mechanism and application.     

Genome-Wide Delineation of Natural Variation for Pod Shatter Resistance in Brassica napus

Resistance to pod shattering (shatter resistance) is a target trait for global rapeseed (canola, Brassica napus L.), improvement programs to minimise grain loss in the mature standing crop, and during windrowing and mechanical harvest. We describe the genetic basis of natural variation for shatter resistance in B. napus and show that several quantitative trait loci (QTL) control this trait. To identify loci underlying shatter resistance, we used a novel genotyping-by-sequencing approach DArT-Seq. QTL analysis detected a total of 12 significant QTL on chromosomes A03, A07, A09, C03, C04, C06, and C08; which jointly account for approximately 57% of the genotypic variation in shatter resistance. Through Genome-Wide Association Studies, we show that a large number of loci, including those that are involved in shattering in Arabidopsis, account for variation in shatter resistance in diverse B. napus germplasm. Our results indicate that genetic diversity for shatter resistance genes in B. napus is limited; many of the genes that might control this trait were not included during the natural creation of this species, or were not retained during the domestication and selection process. We speculate that valuable diversity for this trait was lost during the natural creation of B. napus. To improve shatter resistance, breeders will need to target the introduction of useful alleles especially from genotypes of other related species of Brassica, such as those that we have identified.     

Dissecting Quantitative Trait Loci for Boron Efficiency across Multiple Environments in Brassica napus

High yield is the most important goal in crop breeding, and boron (B) is an essential micronutrient for plants. However, B deficiency, leading to yield decreases, is an agricultural problem worldwide. Brassica napus is one of the most sensitive crops to B deficiency, and considerable genotypic variation exists among different cultivars in response to B deficiency. To dissect the genetic basis of tolerance to B deficiency in B. napus, we carried out QTL analysis for seed yield and yield-related traits under low and normal B conditions using the double haploid population (TNDH) by two-year and the BQDH population by three-year field trials. In total, 80 putative QTLs and 42 epistatic interactions for seed yield, plant height, branch number, pod number, seed number, seed weight and B efficiency coefficient (BEC) were identified under low and normal B conditions, singly explaining 4.15–23.16% and 0.53–14.38% of the phenotypic variation. An additive effect of putative QTLs was a more important controlling factor than the additive-additive effect of epistatic interactions. Four QTL-by-environment interactions and 7 interactions between epistatic interactions and the environment contributed to 1.27–4.95% and 1.17–3.68% of the phenotypic variation, respectively. The chromosome region on A2 of SYLB-A2 for seed yield under low B condition and BEC-A2 for BEC in the two populations was equivalent to the region of a reported major QTL, BE1. The B. napus homologous genes of Bra020592 and Bra020595 mapped to the A2 region and were speculated to be candidate genes for B efficiency. These findings reveal the complex genetic basis of B efficiency in B. napus. They provide a basis for the fine mapping and cloning of the B efficiency genes and for breeding B-efficient cultivars by marker-assisted selection (MAS).     

Gene Silencing of BnTT10 Family Genes Causes Retarded Pigmentation and Lignin Reduction in the Seed Coat of Brassica napus

Yellow-seed (i.e., yellow seed coat) is one of the most important agronomic traits of Brassica plants, which is correlated with seed oil and meal qualities. Previous studies on the Brassicaceae, including Arabidopsis and Brassica species, proposed that the seed-color trait is correlative to flavonoid and lignin biosynthesis, at the molecular level. In Arabidopsis thaliana, the oxidative polymerization of flavonoid and biosynthesis of lignin has been demonstrated to be catalyzed by laccase 15, a functional enzyme encoded by the AtTT10 gene. In this study, eight Brassica TT10 genes (three from B. napus, three from B. rapa and two from B. oleracea) were isolated and their roles in flavonoid oxidation/polymerization and lignin biosynthesis were investigated. Based on our phylogenetic analysis, these genes could be divided into two groups with obvious structural and functional differentiation. Expression studies showed that Brassica TT10 genes are active in developing seeds, but with differential expression patterns in yellow- and black-seeded near-isogenic lines. For functional analyses, three black-seeded B. napus cultivars were chosen for transgenic studies. Transgenic B. napus plants expressing antisense TT10 constructs exhibited retarded pigmentation in the seed coat. Chemical composition analysis revealed increased levels of soluble proanthocyanidins, and decreased extractable lignin in the seed coats of these transgenic plants compared with that of the controls. These findings indicate a role for the Brassica TT10 genes in proanthocyanidin polymerization and lignin biosynthesis, as well as seed coat pigmentation in B. napus.     

Dynamic Metabolic Profiles and Tissue-Specific Source Effects on the Metabolome of Developing Seeds of Brassica napus

Canola (Brassica napus) is one of several important oil-producing crops, and the physiological processes, enzymes, and genes involved in oil synthesis in canola seeds have been well characterized. However, relatively little is known about the dynamic metabolic changes that occur during oil accumulation in seeds, as well as the mechanistic origins of metabolic changes. To explore the metabolic changes that occur during oil accumulation, we isolated metabolites from both seed and silique wall and identified and characterized them by using gas chromatography coupled with mass spectrometry (GC-MS). The results showed that a total of 443 metabolites were identified from four developmental stages. Dozens of these metabolites were differentially expressed during seed ripening, including 20 known to be involved in seed development. To investigate the contribution of tissue-specific carbon sources to the biosynthesis of these metabolites, we examined the metabolic changes of silique walls and seeds under three treatments: leaf-detachment (Ld), phloem-peeling (Pe), and selective silique darkening (Sd). Our study demonstrated that the oil content was independent of leaf photosynthesis and phloem transport during oil accumulation, but required the metabolic influx from the silique wall. Notably, Sd treatment resulted in seed senescence, which eventually led to a severe reduction of the oil content. Sd treatment also caused a significant accumulation of fatty acids (FA), organic acids and amino acids. Furthermore, an unexpected accumulation of sugar derivatives and organic acid was observed in the Pe- and Sd-treated seeds. Consistent with this, the expression of a subset of genes involved in FA metabolism, sugar and oil storage was significantly altered in Pe and Sd treated seeds. Taken together, our studies suggest the metabolite profiles of canola seeds dynamically varied during the course of oil accumulation, which may provide a new insight into the mechanisms of the oil accumulation at the metabolite level.     

Overexpression of phyA and appA Genes Improves Soil Organic Phosphorus Utilisation and Seed Phytase Activity in Brassica napus

Phytate is the major storage form of organic phosphorus in soils and plant seeds, and phosphorus (P) in this form is unavailable to plants or monogastric animals. In the present study, the phytase genes phyA and appA were introduced into Brassica napus cv Westar with a signal peptide sequence and CaMV 35S promoter, respectively. Three independent transgenic lines, P3 and P11 from phyA and a18 from appA, were selected. The three transgenic lines exhibited significantly higher exuded phytase activity when compared to wild-type (WT) controls. A quartz sand culture experiment demonstrated that transgenic Brassica napus had significantly improved P uptake and plant biomass. A soil culture experiment revealed that seed yields of transgenic lines P11 and a18 increased by 20.9% and 59.9%, respectively, when compared to WT. When phytate was used as the sole P source, P accumulation in seeds increased by 20.6% and 46.9% with respect to WT in P11 and a18, respectively. The P3 line accumulated markedly more P in seeds than WT, while no significant difference was observed in seed yields when phytate was used as the sole P source. Phytase activities in transgenic canola seeds ranged from 1,138 to 1,605 U kg^–1 seeds, while no phytase activity was detected in WT seeds. Moreover, phytic acid content in P11 and a18 seeds was significantly lower than in WT. These results introduce an opportunity for improvement of soil and seed phytate-P bioavailability through genetic manipulation of oilseed rape, thereby increasing plant production and P nutrition for monogastric animals.     

Prediction of Protein Binding Regions in Disordered Proteins

Many disordered proteins function via binding to a structured partner and undergo a disorder-to-order transition. The coupled folding and binding can confer several functional advantages such as the precise control of binding specificity without increased affinity. Additionally, the inherent flexibility allows the binding site to adopt various conformations and to bind to multiple partners. These features explain the prevalence of such binding elements in signaling and regulatory processes. In this work, we report ANCHOR, a method for the prediction of disordered binding regions. ANCHOR relies on the pairwise energy estimation approach that is the basis of IUPred, a previous general disorder prediction method. In order to predict disordered binding regions, we seek to identify segments that are in disordered regions, cannot form enough favorable intrachain interactions to fold on their own, and are likely to gain stabilizing energy by interacting with a globular protein partner. The performance of ANCHOR was found to be largely independent from the amino acid composition and adopted secondary structure. Longer binding sites generally were predicted to be segmented, in agreement with available experimentally characterized examples. Scanning several hundred proteomes showed that the occurrence of disordered binding sites increased with the complexity of the organisms even compared to disordered regions in general. Furthermore, the length distribution of binding sites was different from disordered protein regions in general and was dominated by shorter segments. These results underline the importance of disordered proteins and protein segments in establishing new binding regions. Due to their specific biophysical properties, disordered binding sites generally carry a robust sequence signal, and this signal is efficiently captured by our method. Through its generality, ANCHOR opens new ways to study the essential functional sites of disordered proteins.     

Computational Design of Binding Proteins to EGFR Domain II

We developed a process to produce novel interactions between two previously unrelated proteins. This process selects protein scaffolds and designs protein interfaces that bind to a surface patch of interest on a target protein. Scaffolds with shapes complementary to the target surface patch were screened using an exhaustive computational search of the human proteome and optimized by directed evolution using phage display. This method was applied to successfully design scaffolds that bind to epidermal growth factor receptor (EGFR) domain II, the interface of EGFR dimerization, with high reactivity toward the target surface patch of EGFR domain II. One potential application of these tailor-made protein interactions is the development of therapeutic agents against specific protein targets.     

BraLTP1, a Lipid Transfer Protein Gene Involved in Epicuticular Wax Deposition,Cell Proliferation and Flower Development in Brassica napus

Plant non-specific lipid transfer proteins (nsLTPs) constitute large multigene families that possess complex physiological functions, many of which remain unclear. This study isolated and characterized the function of a lipid transfer protein gene, BraLTP1 from Brassica rapa, in the important oilseed crops Brassica napus. BraLTP1 encodes a predicted secretory protein, in the little known VI Class of nsLTP families. Overexpression of BnaLTP1 in B. napus caused abnormal green coloration and reduced wax deposition on leaves and detailed wax analysis revealed 17–80% reduction in various major wax components, which resulted in significant water-loss relative to wild type. BnaLTP1 overexpressing leaves exhibited morphological disfiguration and abaxially curled leaf edges, and leaf cross-sections revealed cell overproliferation that was correlated to increased cytokinin levels (tZ, tZR, iP, and iPR) in leaves and high expression of the cytokinin biosynthsis gene IPT3. BnaLTP1-overexpressing plants also displayed morphological disfiguration of flowers, with early-onset and elongated carpel development and outwardly curled stamen. This was consistent with altered expression of a a number of ABC model genes related to flower development. Together, these results suggest that BraLTP1 is a new nsLTP gene involved in wax production or deposition, with additional direct or indirect effects on cell division and flower development.     

转录因子结合位点的计算预测方法研究进展

转录因子结合位点的计算预测是研究基因转录调控的重要环节,但现有算法的预测特异性偏低.在深入分析转录因子结合位点生物特征的基础上,对当前基于保守模体和基于比较基因组学的两类计算预测方法进行了综述,指出了方法各自的优点和不足,并探讨了可能的改进方向.     

Predicting DNA-Binding Proteins and Binding Residues by Complex Structure Prediction and Application to Human Proteome

As more and more protein sequences are uncovered from increasingly inexpensive sequencing techniques, an urgent task is to find their functions. This work presents a highly reliable computational technique for predicting DNA-binding function at the level of protein-DNA complex structures, rather than low-resolution two-state prediction of DNA-binding as most existing techniques do. The method first predicts protein-DNA complex structure by utilizing the template-based structure prediction technique HHblits, followed by binding affinity prediction based on a knowledge-based energy function (Distance-scaled finite ideal-gas reference state for protein-DNA interactions). A leave-one-out cross validation of the method based on 179 DNA-binding and 3797 non-binding protein domains achieves a Matthews correlation coefficient (MCC) of 0.77 with high precision (94%) and high sensitivity (65%). We further found 51% sensitivity for 82 newly determined structures of DNA-binding proteins and 56% sensitivity for the human proteome. In addition, the method provides a reasonably accurate prediction of DNA-binding residues in proteins based on predicted DNA-binding complex structures. Its application to human proteome leads to more than 300 novel DNA-binding proteins; some of these predicted structures were validated by known structures of homologous proteins in APO forms. The method [SPOT-Seq (DNA)] is available as an on-line server at http://sparks-lab.org.     

Pasteurellaceae ComE1 Proteins Combine the Properties of Fibronectin Adhesins and DNA Binding Competence Proteins

A novel fibronectin-binding protein from Pasteurella multocida (PM1665) that binds to the fibronectin type III_9-10 modules via two helix-hairpin-helix motifs has recently been described [1]. This protein shares homology with competence-related DNA-binding and uptake proteins (ComEA and ComE) from Gram-positive and Gram-negative bacteria. Here, we show that recombinant PM1665 (now designated ComE1) also binds to DNA through the same helix-hairpin-helix motifs required for fibronectin-binding. This binding to DNA is non sequence-specific and is confined to double-stranded DNA. We have cloned and expressed ComE1 proteins from five members of the Pasteurellaceae in order to further investigate the function(s) of these proteins. When expressed as recombinant GST-fusion proteins, all of the homologues bound both to fibronectin and to double-stranded DNA. Inactivation of the gene encoding the ComE1 homologue in Actinobacillus pleuropneumoniae indicates major roles for these proteins in at least two processes: natural transformation, and binding of bacteria to fibronectin.     

Comparison of Sugars,Iridoid Glycosides and Amino Acids in Nectar and Phloem Sap of Maurandya barclayana,Lophospermum erubescens,and Brassica napus

Background

Floral nectar contains sugars and amino acids to attract pollinators. In addition, nectar also contains different secondary compounds, but little is understood about their origin or function. Does nectar composition reflect phloem composition, or is nectar synthesized and/or modified in nectaries? Studies where both, the nectar as well as the phloem sap taken from the same plant species were analyzed in parallel are rare. Therefore, phloem sap and nectar from different plant species (Maurandya barclayana, Lophospermum erubescens, and Brassica napus) were compared.

Methodology and Principal Findings

Nectar was collected with microcapillary tubes and phloem sap with the laser-aphid-stylet technique. The nectar of all three plant species contained high amounts of sugars with different percentages of glucose, fructose, and sucrose, whereas phloem sap sugars consisted almost exclusively of sucrose. One possible reason for this could be the activity of invertases in the nectaries. The total concentration of amino acids was much lower in nectars than in phloem sap, indicating selective retention of nitrogenous solutes during nectar formation. Nectar amino acid concentrations were negatively correlated with the nectar volumes per flower of the different plant species. Both members of the tribe Antirrhineae (Plantaginaceae) M. barclayana and L. erubescens synthesized the iridoid glycoside antirrhinoside. High amounts of antirrhinoside were found in the phloem sap and lower amounts in the nectar of both plant species.

Conclusions/Significance

The parallel analyses of nectar and phloem sap have shown that all metabolites which were found in nectar were also detectable in phloem sap with the exception of hexoses. Otherwise, the composition of both aqueous solutions was not the same. The concentration of several metabolites was lower in nectar than in phloem sap indicating selective retention of some metabolites. Furthermore, the existence of antirrhinoside in nectar could be based on passive secretion from the phloem.     

Computational Prediction of Alanine Scanning and Ligand Binding Energetics in G-Protein Coupled Receptors

Site-directed mutagenesis combined with binding affinity measurements is widely used to probe the nature of ligand interactions with GPCRs. Such experiments, as well as structure-activity relationships for series of ligands, are usually interpreted with computationally derived models of ligand binding modes. However, systematic approaches for accurate calculations of the corresponding binding free energies are still lacking. Here, we report a computational strategy to quantitatively predict the effects of alanine scanning and ligand modifications based on molecular dynamics free energy simulations. A smooth stepwise scheme for free energy perturbation calculations is derived and applied to a series of thirteen alanine mutations of the human neuropeptide Y1 receptor and series of eight analogous antagonists. The robustness and accuracy of the method enables univocal interpretation of existing mutagenesis and binding data. We show how these calculations can be used to validate structural models and demonstrate their ability to discriminate against suboptimal ones.     

氧化固醇结合蛋白结构、功能与应用

氧化固醇结合蛋白(oxysterol binding protein,OSBP)是存在于真核细胞内的一类参与脂质代谢的非囊泡运输蛋白质,在哺乳动物中被称为氧化固醇结合蛋白相关蛋白质(oxysterol binding protein-related proteins,ORPs),而在酵母中被称为氧化固醇结合蛋白同源物质(oxysterol-binding protein homologues,OSH）。近年来人们对氧化固醇结合蛋白的研究不断深入,特别是对其同源蛋白质（例如,ORP5/8、Osh3/4、ORP4L等）的结构功能差异和其在信号转导中的作用的相关研究,以及在生物医药方面的应用更成为了本领域的热点。本文综述了关于OSBP及其同源蛋白质结构和功能的相关研究,指出了该领域存在的一些关键问题。与此同时,对OSH和ORPs在细胞内的膜接触位点（membrane contact sites,MCS）进行对比,以及对今后OSBP的研究方向做了展望。