N-terminal domain of ARF-GEF GNOM prevents heterodimerization with functionally divergent GNL1 in Arabidopsis.

Evolutionary change following gene duplication can lead to functionally divergent paralogous proteins. If comprising identical subunits their random assortment would also form potentially detrimental heteromeric proteins. In Arabidopsis, the ARF GTPase guanine-nucleotide exchange factor GNOM is essential for polar recycling of auxin-efflux transporter PIN1 from endosomes to the basal plasma membrane whereas its paralog GNL1 mediates retrograde Golgi-endoplasmic reticulum traffic. Here we show that both GNOM and GNL1 form homodimers but no heterodimers. To assess the biological significance of this, we generated transgenic plants expressing engineered heterodimer-compatible GNOM variants. Those plants showed developmental defects such as the failure to produce lateral roots. To identify mechanisms underlying heterodimer prevention, we analyzed interactions of the N-terminal dimerization and cyclophilin-binding (DCB) domain. Each DCB domain interacted with the complementary fragment (ΔDCB) both of their own and of the paralogous protein. However, only DCBGNOM interacted with itself whereas DCBGNL1 failed to interact with itself and with DCBGNOM . GNOM variants in which the DCB domain was removed or replaced by DCBGNL1 revealed a role for DCB-DCB interaction in the prevention of GNOM-GNL1 heterodimers whereas DCB-ΔDCB interaction was essential for dimer formation and GNOM function. Our data suggest a model of early DCB-DCB interaction that facilitates GNOM homodimer formation, indirectly precluding formation of detrimental heterodimers.

Determinants of Lipid Domain Size.

Lipid domains less than 200 nm in size may form a scaffold, enabling the concerted function of plasma membrane proteins. The size-regulating mechanism is under debate. We tested the hypotheses that large values of spontaneous monolayer curvature are incompatible with micrometer-sized domains. Here, we used the transition of photoswitchable lipids from their cylindrical conformation to a conical conformation to increase the negative curvature of a bilayer-forming lipid mixture. In contrast to the hypothesis, pre-existing micrometer-sized domains did not dissipate in our planar bilayers, as indicated by fluorescence images and domain mobility measurements. Elasticity theory supports the observation by predicting the zero free energy gain for splitting large domains into smaller ones. It also indicates an alternative size-determining mechanism: The cone-shaped photolipids reduce the line tension associated with lipid deformations at the phase boundary and thus slow down the kinetics of domain fusion. The competing influence of two approaching domains on the deformation of the intervening lipids is responsible for the kinetic fusion trap. Our experiments indicate that the resulting local energy barrier may restrict the domain size in a dynamic system.

Effects of conserved Arg20, Glu74 and Asp77 on the structure and function of a tau class glutathione S-transferase in rice.

KEY MESSAGE: The relative position of domains is critical for enzymatic properties of tau class glutathione S-transferases, and altering the position of linker far away from the active center affects catalytic property. Glutathione S-transferases (GSTs) are a family of phase II detoxification enzymes whose main function is to improve plant resistance to stresses. To understand the structural effects of tau class GSTs on their function, using OsGSTU17 as an example, we predicted the residues involved in the interactions between its domains and linker region. We further detected the structural changes in mutants and the corresponding changes in terms of substrate activity and kinetic parameters. Four pairs of residues, including Ala14 and Trp165, Arg20 and Tyr154, Glu74 and Arg98, Asp77 and Met87, forming hydrogen bonds and salt bridges were found to play important roles in maintaining the relative position between the domains and linker region inside the protein. The hydrogen bond between Trp165 and Ala14 affected the structural stability has been demonstrated in our previous study. The mutant R20A lost almost all catalytic activity. Interestingly, the mutant E74A exhibited a significant decrease in activity towards 7-chloro-4-nitrobenzo-2-oxa-1, 3-diazole, 1-chloro-2, 4-dinitrobenzene and 4-nitrobenzyl chloride, while its activity towards substrate cumene hydroperoxide remained unchanged. Compared with other mutants, the mutant D77A exhibited decreased affinity to its substrates and increased activity towards 1-chloro-2, 4-dinitrobenzene and cumene hydroperoxide, but its thermodynamic stability did not change significantly. The relative position of individual domain was critical for enzymatic properties, and the linker which is far away from the active site could change the enzymatic properties of GSTs via altering the relative position of the individual domain. Our results provide insights into the relationship between structure and function of tau class GSTs.

Domain interactions reveal auto-inhibition of the deubiquitinating enzyme USP19 and its activation by HSP90 in the modulation of huntingtin aggregation.

Ubiquitin-specific protease 19 (USP19) is a member of the deubiquitinating (DUB) enzymes that catalyze removing the ubiquitin signals from target proteins. Our previous research has demonstrated that USP19 up-regulates the protein level and aggregation of polyQ-expanded huntingtin through the involvement of heat shock protein 90 (HSP90). Here, we present solution structures of the CS1, CS2 and UbL domains of USP19 and structural insights into their domain interactions. We found that the tandem CS domains fold back to interact with the C-terminal USP domain (USPD) intra-molecularly that leads to inhibition of the catalytic core of USP19, especially CS1 interacts with the embedded UbL domain and CS2 does with the CH2 catalytic core. Moreover, CS2 specifically interacts with the NBD domain of HSP90, which can activate the DUB enzyme. A mechanism of auto-inhibition of USP19 and activation by HSP90 is proposed, on which USP19 modulates the protein level of polyQ-expanded huntingtin in cells. This study provides structural and mechanistic insights into the modulation of protein level and aggregation by USP19 with the assistance of HSP90.

PlaPPISite: a comprehensive resource for plant protein-protein interaction sites.

BACKGROUND: Protein-protein interactions (PPIs) play very important roles in diverse biological processes. Experimentally validated or predicted PPI data have become increasingly available in diverse plant species. To further explore the biological functions of PPIs, understanding the interaction details of plant PPIs (e.g., the 3D structural contexts of interaction sites) is necessary. By integrating bioinformatics algorithms, interaction details can be annotated at different levels and then compiled into user-friendly databases. In our previous study, we developed AraPPISite, which aimed to provide interaction site information for PPIs in the model plant Arabidopsis thaliana. Considering that the application of AraPPISite is limited to one species, it is very natural that AraPPISite should be evolved into a new database that can provide interaction details of PPIs in multiple plants. DESCRIPTION: PlaPPISite (http://zzdlab.com/plappisite/index.php) is a comprehensive, high-coverage and interaction details-oriented database for 13 plant interactomes. In addition to collecting 121 experimentally verified structures of protein complexes, the complex structures of experimental/predicted PPIs in the 13 plants were also constructed, and the corresponding interaction sites were annotated. For the PPIs whose 3D structures could not be modelled, the associated domain-domain interactions (DDIs) and domain-motif interactions (DMIs) were inferred. To facilitate the reliability assessment of predicted PPIs, the source species of interolog templates, GO annotations, subcellular localizations and gene expression similarities are also provided. JavaScript packages were employed to visualize structures of protein complexes, protein interaction sites and protein interaction networks. We also developed an online tool for homology modelling and protein interaction site annotation of protein complexes. All data contained in PlaPPISite are also freely available on the Download page. CONCLUSION: PlaPPISite provides the plant research community with an easy-to-use and comprehensive data resource for the search and analysis of protein interaction details from the 13 important plant species.

Cation Specific Effects on the Domain-Domain Interaction of Heterogeneous Dimeric Protein Revealed by FRET Analysis.

Specific monovalent cation effects on the domain-domain interaction of heterogeneous dimeric protein were investigated using green fluorescent protein (GFP)-glutathione-s-transferase (GST) fusion protein as a model protein. Conjugating N-terminal of GST domain with a fluorescence probe Cyanine3, complementary increase and decrease of fluorescence intensities of Cyanine3 and GFP were recognized on the exclusive excitation of GFP and further the fluorescence decay of GFP was remarkably accelerated to show that an excellent Förster type of resonance excitation energy transfer (FRET) pair was constructed between GFP- and GST-domain. The spectral overlap integral and critical distance of the FRET pair were estimated to be 5.96×1013 M-1cm3 and 62.5 Å, respectively. The FRET rate and efficiency evaluated by fluorescence lifetime of the energy donor, GFP, were influenced by the monovalent cations included in the buffer solution to suggest that the domain-domain interactions of GFP-GST fusion protein would be susceptible to cation species and their concentrations. The order affecting the domain-domain interaction was estimated to be Li+>NH4+ >Na+>K+>Cs+, almost corresponding to the reverse Hofmeister series.

A mechanism for lipid binding to apoE and the role of intrinsically disordered regions coupled to domain-domain interactions.

Relative to the apolipoprotein E (apoE) E3 allele of the APOE gene, apoE4 strongly increases the risk for the development of late-onset Alzheimer's disease. However, apoE4 differs from apoE3 by only a single amino acid at position 112, which is arginine in apoE4 and cysteine in apoE3. It remains unclear why apoE3 and apoE4 are functionally different. Described here is a proposal for understanding the functional differences between these two isoforms with respect to lipid binding. A mechanism is proposed that is based on the full-length monomeric structure of the protein, on hydrogen-deuterium exchange mass spectrometry data, and on the role of intrinsically disordered regions to control protein motions. It is proposed that lipid binds between the N-terminal and C-terminal domains and that separation of the two domains, along with the presence of intrinsically disordered regions, controls this process. The mechanism explains why apoE3 differs from apoE4 with respect to different lipid-binding specificities, why lipid increases the binding of apoE to its receptor, and why specific residues are conserved.

Female mice with apolipoprotein E4 domain interaction demonstrated impairments in spatial learning and memory performance and disruption of hippocampal cyto-architecture.

We have previously reported cognitive impairments in both young and old mice, particularly in female mice expressing mouse Arg-61 apoE, with a point mutation to mimic the domain interaction feature of human apoE4, as compared to the wildtype mouse (C57BL/6J) apoE. In this study, we further evaluated water maze performance in the female Arg-61 mice at an additional time point and then investigated related hippocampal cyto-architecture in these young female Arg-61 apoE mice vs. the wildtype mice. The results of behavioral performance consistently support our previous report that the young female Arg-61 apoE showed cognitive impairment versus C57BL/6J at the same age. The cyto-architectural results showed that volume of the granular cell layer (GCL) was significantly larger in both 5- and 10-month old Arg-61 apoE mice versus C57BL/6J mice. While the number of newborn calretinin-positive neurons was greater in the sub-granular zone (SGZ) in 5-month old Arg-61 mice, this number dropped significantly in 10-month old Arg-61 mice to a lower level than in age-matched C57BL/6J mice. In addition, the amyloid ß species was significantly higher in 5-month old Arg-61 mice versus age-matched C57BL/6J mice. In conclusion, impaired cognitive functions in female Arg-61 apoE mice appear correlated with larger GCL volume and higher calretinin-positive cell number and suggest a compensatory cellular response that may be related to amyloid beta perturbations early in life. Therefore this study suggests a novel cyto-architectural mechanism of apoE4-dependent pathologies and increased susceptibility of APOEε4 subjects to Alzheimer's disease.

Self-binding peptides: Binding-upon-folding versus folding-upon-binding.

Self-binding peptide (SBP) represents a novel biomolecular phenomenon spanning between folding and binding. It is a structurally independent, short peptide segment within a monomeric protein and fulfills biological function by dynamically binding to/unbinding from its target domain in the same monomer. Here, four representative SBP systems, including mouse proto-oncogene Vav, human retinoic acid receptor RARγ, fruit fly scaffold module INAD and crypto 14-3-3 protein Cp14b, are investigated systematically by using atomistic molecular dynamics (MD) simulations and post binding energetics analyses. The native bound structure, artificial unbound state and isolated peptide segment of SBP moieties in the four systems were constructed, analyzed and compared in detail. It is revealed that the SBP interaction with their targets is almost a binding phenomenon at single-molecule level, but presence of a polypeptide linker between the SBP and target can promote the binding efficiency since the linker restriction largely increases the probability of SBP-target encounters in a statistical physics point of view. In this respect, unlike classical peptide-mediated interactions where the intrinsically disordered peptides are folded into an ordered structure upon binding to their protein partners (folding-upon-binding), we herein propose SBPs as a new and reversed biological event that is naturally a folding phenomenon but exhibits a typical binding behavior (binding-upon-folding).

Domain-based Comparative Analysis of Bacterial Proteomes: Uniqueness, Interactions, and the Dark Matter.

BACKGROUND: Proteins may have none, single, double, or multiple domains, while a single domain may appear in multiple proteins. Their distribution patterns may have impacts on bacterial physi-ology and lifestyle.Objective: This study aims to understand how domains are distributed and duplicated in bacterial prote-omes, in order to better understand bacterial physiology and lifestyles. METHODS: In this study, we used 16712 Hidden Markov Models to screen 944 bacterial reference prote-omes versus a threshold E-value<0.001. The number of non-redundant domains and duplication rates of redundant domains for each species were calculated. The unique domains, if any, were also identified for each species. In addition, the properties of no-domain proteins were investigated in terms of physico-chemical properties. RESULTS: The increasing number of non-redundant domains for a bacterial proteome follows the trend of an asymptotic function. The domain duplication rate is positively correlated with proteome size and in-creases more rapidly. The high percentage of single-domain proteins is more associated with small pro-teome size. For each proteome, unique domains were also obtained. Moreover, no-domain proteins show differences with the other three groups for several physicochemical properties analysed in this study. CONCLUSION: The study confirmed that a low domain duplication rate and a high percentage of single-domain proteins are more likely to be associated with bacterial host-dependent or restricted niche-adapted lifestyle. In addition, the unique lifestyle and physiology were revealed based on the analysis of species-specific domains and core domain interactions or co-occurrences.

The polarity protein Scribble positions DLC3 at adherens junctions to regulate Rho signaling.

The spatial regulation of cellular Rho signaling by GAP proteins is still poorly understood. By performing mass spectrometry, we here identify the polarity protein Scribble as a scaffold for the RhoGAP protein DLC3 (also known as StarD8) at cell-cell adhesions. This mutually dependent interaction is mediated by the PDZ domains of Scribble and a PDZ ligand (PDZL) motif in DLC3. Both Scribble depletion and PDZL deletion abrogated DLC3 junctional localization. Using a RhoA biosensor and a targeted GAP domain, we demonstrate that DLC3 activity locally regulates RhoA-ROCK signaling at and Scribble localization to adherens junctions, and is required for their functional integrity. In a 3D model of cyst development, we furthermore show that DLC3 depletion impairs polarized morphogenesis, phenocopying the effects observed upon Scribble knockdown. We thus propose a new function for Scribble in Rho regulation that entails positioning of DLC3 GAP activity at cell junctions in polarized epithelial cells.

Single-molecule spectroscopy reveals how calmodulin activates NO synthase by controlling its conformational fluctuation dynamics.

Mechanisms that regulate the nitric oxide synthase enzymes (NOS) are of interest in biology and medicine. Although NOS catalysis relies on domain motions, and is activated by calmodulin binding, the relationships are unclear. We used single-molecule fluorescence resonance energy transfer (FRET) spectroscopy to elucidate the conformational states distribution and associated conformational fluctuation dynamics of the two electron transfer domains in a FRET dye-labeled neuronal NOS reductase domain, and to understand how calmodulin affects the dynamics to regulate catalysis. We found that calmodulin alters NOS conformational behaviors in several ways: It changes the distance distribution between the NOS domains, shortens the lifetimes of the individual conformational states, and instills conformational discipline by greatly narrowing the distributions of the conformational states and fluctuation rates. This information was specifically obtainable only by single-molecule spectroscopic measurements, and reveals how calmodulin promotes catalysis by shaping the physical and temporal conformational behaviors of NOS.

A multi-network clustering method for detecting protein complexes from multiple heterogeneous networks.

BACKGROUND: The accurate identification of protein complexes is important for the understanding of cellular organization. Up to now, computational methods for protein complex detection are mostly focus on mining clusters from protein-protein interaction (PPI) networks. However, PPI data collected by high-throughput experimental techniques are known to be quite noisy. It is hard to achieve reliable prediction results by simply applying computational methods on PPI data. Behind protein interactions, there are protein domains that interact with each other. Therefore, based on domain-protein associations, the joint analysis of PPIs and domain-domain interactions (DDI) has the potential to obtain better performance in protein complex detection. As traditional computational methods are designed to detect protein complexes from a single PPI network, it is necessary to design a new algorithm that could effectively utilize the information inherent in multiple heterogeneous networks. RESULTS: In this paper, we introduce a novel multi-network clustering algorithm to detect protein complexes from multiple heterogeneous networks. Unlike existing protein complex identification algorithms that focus on the analysis of a single PPI network, our model can jointly exploit the information inherent in PPI and DDI data to achieve more reliable prediction results. Extensive experiment results on real-world data sets demonstrate that our method can predict protein complexes more accurately than other state-of-the-art protein complex identification algorithms. CONCLUSIONS: In this work, we demonstrate that the joint analysis of PPI network and DDI network can help to improve the accuracy of protein complex detection.

Pathogenicity Genes in Ustilaginoidea virens Revealed by a Predicted Protein-Protein Interaction Network.

Rice false smut, caused by Ustilaginoidea virens, produces significant losses in rice yield and grain quality and has recently emerged as one of the most important rice diseases worldwide. Despite its importance in rice production, relatively few studies have been conducted to illustrate the complex interactome and the pathogenicity gene interactions. Here a protein-protein interaction network of U. virens was built through two well-recognized approaches, interolog- and domain-domain interaction-based methods. A total of 20 217 interactions associated with 3305 proteins were predicted after strict filtering. The reliability of the network was assessed computationally and experimentally. The topology of the interactome network revealed highly connected proteins. A pathogenicity-related subnetwork involving up-regulated genes during early U. virens infection was also constructed, and many novel pathogenicity proteins were predicted in the subnetwork. In addition, we built an interspecies PPI network between U. virens and Oryza sativa, providing new insights for molecular interactions of this host-pathogen pathosystem. A web-based publicly available interactive database based on these interaction networks has also been released. In summary, a proteome-scale map of the PPI network was described for U. virens, which will provide new perspectives for finely dissecting interactions of genes related to its pathogenicity.

Channel Gating Regulation by the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) First Cytosolic Loop.

In this study, we present data indicating a robust and specific domain interaction between the cystic fibrosis transmembrane conductance regulator (CFTR) first cytosolic loop (CL1) and nucleotide binding domain 1 (NBD1) that allows ion transport to proceed in a regulated fashion. We used co-precipitation and ELISA to establish the molecular contact and showed that binding kinetics were not altered by the common clinical mutation F508del. Both intrinsic ATPase activity and CFTR channel gating were inhibited severely by CL1 peptide, suggesting that NBD1/CL1 binding is a crucial requirement for ATP hydrolysis and channel function. In addition to cystic fibrosis, CFTR dysregulation has been implicated in the pathogenesis of prevalent diseases such as chronic obstructive pulmonary disease, acquired rhinosinusitis, pancreatitis, and lethal secretory diarrhea (e.g. cholera). On the basis of clinical relevance of the CFTR as a therapeutic target, a cell-free drug screen was established to identify modulators of NBD1/CL1 channel activity independent of F508del CFTR and pharmacologic rescue. Our findings support a targetable mechanism of CFTR regulation in which conformational changes in the NBDs cause reorientation of transmembrane domains via interactions with CL1 and result in channel gating.

A conserved π-cation and an electrostatic bridge are essential for 11R-lipoxygenase catalysis and structural stability.

Lipoxygenases (LOXs) are lipid-peroxidizing enzymes that consist of a regulatory calcium- and membrane-binding PLAT (polycystin-1, lipoxygenase, α-toxin) domain and a catalytic domain. In a previous study, the crystal structure of an 11R-LOX revealed a conserved π-cation bridge connecting these two domains which could mediate the regulatory effect of the PLAT domain to the active site. Here we analyzed the role of residues Trp107 and Lys172 that constitute the π-cation bridge in 11R-LOX along with Arg106 and Asp173-a potential salt bridge, which could also contribute to the inter-domain communication. According to our kinetic assays and protein unfolding experiments conducted using differential scanning fluorimetry and circular dichroism spectroscopy, mutants with a disrupted link display diminished catalytic activity alongside reduced stability of the protein fold. The results demonstrate that both these bridges contribute to the two-domain interface, and are important for proper enzyme activation.

AraPPISite: a database of fine-grained protein-protein interaction site annotations for Arabidopsis thaliana.

Knowledge about protein interaction sites provides detailed information of protein-protein interactions (PPIs). To date, nearly 20,000 of PPIs from Arabidopsis thaliana have been identified. Nevertheless, the interaction site information has been largely missed by previously published PPI databases. Here, AraPPISite, a database that presents fine-grained interaction details for A. thaliana PPIs is established. First, the experimentally determined 3D structures of 27 A. thaliana PPIs are collected from the Protein Data Bank database and the predicted 3D structures of 3023 A. thaliana PPIs are modeled by using two well-established template-based docking methods. For each experimental/predicted complex structure, AraPPISite not only provides an interactive user interface for browsing interaction sites, but also lists detailed evolutionary and physicochemical properties of these sites. Second, AraPPISite assigns domain-domain interactions or domain-motif interactions to 4286 PPIs whose 3D structures cannot be modeled. In this case, users can easily query protein interaction regions at the sequence level. AraPPISite is a free and user-friendly database, which does not require user registration or any configuration on local machines. We anticipate AraPPISite can serve as a helpful database resource for the users with less experience in structural biology or protein bioinformatics to probe the details of PPIs, and thus accelerate the studies of plant genetics and functional genomics. AraPPISite is available at http://systbio.cau.edu.cn/arappisite/index.html .

Association analysis of the perturbation of interactions in biological pathways and anticancer drug activity.

Understanding how different genomic mutational landscapes in patients with cancer lead to different responses to anticancer drugs is an important challenge for realizing precision medicine for cancer. Many studies have analyzed the comprehensive anticancer drug-response profiles and genomic profiles of cancer cell lines to identify the relationship between the anticancer drug response and genomic alternations. However, few studies have focused on interpreting these profiles with a network perspective. In this work, we analyzed genomic alterations in cancer cell lines by considering which interactions in the signaling pathway were perturbed by mutations. With our interaction-centric approach, we identified novel interaction/drug response associations for two drugs (afatinib and ixabepilone) for which no gene-centric association could be found. When we compared the performance of classifiers for predicting the responses to 164 drugs, the classifiers trained with interaction-centric features outperformed the classifiers trained with gene-centric features, despite the smaller number of features (p-value = 2.0 × 10(-3)). By incorporating the interaction information from signaling pathways, we revealed associations between genomic alterations and drug responses that could be missed when using a gene-centric approach.

Inter-domain interactions of TDP-43 as decoded by NMR.

TDP-43 inclusions have been found in ∼97% ALS as well as an increasing spectrum of other neurodegenerative diseases including Alzheimer's. TDP-43 contains an ubiquitin-like fold, two RRMs and a prion-like domain, but whether they interact with each other remains unknown due to being intrinsically aggregation-prone. Nevertheless, this knowledge is pivotal to understanding physiological functions and pathological roles of TDP-43. Here as facilitated by our previous discovery which allowed NMR characterization of TDP-43 and its five dissected fragments, we successfully decoded that TDP-43 does have dynamic inter-domain interactions, which are coordinated by the intrinsically-disordered prion-like domain. Thus, TDP-43 appears to undergo conformational exchanges between "closed" and "open" states which are needed for its functions. Our study thus offers a mechanism by which cellular processes might control TDP-43 physiology and proteinopathy by mediating its inter-domain interactions.

SPAR: a random forest-based predictor for self-interacting proteins with fine-grained domain information.

Protein self-interaction, i.e. the interaction between two or more identical proteins expressed by one gene, plays an important role in the regulation of cellular functions. Considering the limitations of experimental self-interaction identification, it is necessary to design specific bioinformatics tools for self-interacting protein (SIP) prediction from protein sequence information. In this study, we proposed an improved computational approach for SIP prediction, termed SPAR (Self-interacting Protein Analysis serveR). Firstly, we developed an improved encoding scheme named critical residues substitution (CRS), in which the fine-grained domain-domain interaction information was taken into account. Then, by employing the Random Forest algorithm, the performance of CRS was evaluated and compared with several other encoding schemes commonly used for sequence-based protein-protein interaction prediction. Through the tenfold cross-validation tests on a balanced training dataset, CRS performed the best, with the average accuracy up to 72.01 %. We further integrated CRS with other encoding schemes and identified the most important features using the mRMR (the minimum redundancy maximum relevance) feature selection method. Our SPAR model with selected features achieved an average accuracy of 92.09 % on the human-independent test set (the ratio of positives to negatives was about 1:11). Besides, we also evaluated the performance of SPAR on an independent yeast test set (the ratio of positives to negatives was about 1:8) and obtained an average accuracy of 76.96 %. The results demonstrate that SPAR is capable of achieving a reasonable performance in cross-species application. The SPAR server is freely available for academic use at http://systbio.cau.edu.cn/zzdlab/spar/ .