Problem to use coevolution server

Q:
I want to use your co-evolution website to find my co-evolution amino
acid in my interested protein. I found my proteins family in different
species by blast. but when i uploaded it to the server. server sent me an
email that told me your sequences are not enough.i loaded an example in your
server to found out how it worked but it had a problem too and error
message sent to me.please help me to how to use this server.

A:
We set the minimum number of sequences as 50 to ensure statistical significance
of the results. If you really want to run with fewer sequences, you can change
the number in the advanced options.

As for the other error, please send me the error message if you still cannot run
tasks successfully.

Query on SCA coevolution

Q:
I have conducted structure based statistical coupling analyses (SCA) on each
of some mitochondrial proteins using 800 multiple sequences (including one
sequence from our organisms, one 3RKO structure sequence, and 788 protein
sequences from different genera), and we could obtain the coevolutionary
scores and spatial distances between any pair of two residues. The aim of
our study is try to analyze the coevolutionary role of some important given
residues (selected by PAML analyses) on key or important residues
responsible for proton translocation in the proton translocating channel of
respiratory Complex I. The problem is we are not sure how to do it in a more
statistical way. Such as, we could have the data of scores and distances of
a given selected residue on these residues in proton channel or other
residues of the same protein. In order to know possible different
coevolutionary role of a given residue i.e. the selective residue on proton
channel residues or other residues, t-test on scores (s), or distances (d)
or sores/distinces (s/d) were compared by us between those types of
residues, we are not sure if this kind of analyse is ok for us. Such as we
don’t know whether the score obtained by SCA analyses in the platform has
alreadly considered the potential role of distance, or it is just the score
obtained no mattter where both residues are? We know the influencing role
between any two given residues might be correlated with both their
characteristics and spatial distance between them.

Do you have any good idea on this, or do you have more reasonable
statistical way to solve our queries and prolem above?

A:
The scores were calculated based on the MSA alone without
considering the spatial distance between residues.

You may want to plot the global distribution of scores, and look
for scores that are significantly larger than the rest but cannot be
explained by the distance on the primary sequence alone. Indirect
coupling between residues though other residues is also something to be
aware of. There have been a lot of new papers about co-evolutionary
analysis lately (e.g., from Rama Ranganathan’s and Debora Marks’s labs).

Java chromod package request CoassociationAnalyzer.java and GSCCoassociationAnalyzer.java scripts that Kevin Yip wrote (April 14, 2011)

Q:
I’m writing to you to see if you could share with me your java "chromod" package – I’m wanting to use the CoassociationAnalyzer.java and GSCCoassociationAnalyzer.java scripts that Kevin Yip wrote (April 14, 2011), but they rely on the chromod package (package org.gersteinlab.chromod)

If you could share this with me if it’s not a top secret lab package, I would be hugely indebted!

A:
Please download it at http://www.cse.cuhk.edu.hk/~kevinyip/outbox/chromod.jar . Let me know if you encounter any problem when using it.

Submitting MSA to coevolution webpage

Q:
I am trying to submit an MSA to your "coevolution" webpage — and it keeps failing on me.
i keep getting an error email saying "could not be completed" Error message:Not enough sequences.

im sure that my MSA is in fasta format… so not really sure what is going wrong…
if you could help me out, it would be much appreciated.

A:
The tool performs a number of filtering steps, to ensure the
reliability of the results. The error messages states that after
filtering, the number of remaining sequences is too small for a
reliable analysis. You may change the filtering criteria using the
advanced options (hidden by default), but please notice that by doing
so you may get results that are unreliable.

Task failed coevolution

Q:
The seed alignment of PF01036 has been changed recently.
BACR_HALSA is no longer one of the seed sequences. I have just changed
the example to use BACR_HALAR instead, and the program ran fine.
Please let me know if you encounter any problem running the program on
your own data.

A:
I’ve been trying to use your server, but evidently it is not working
correctly; your example data even fails to process.

I am eager to use this server (and likewise cite you), so please let
me know at your earliest convenience if/how I can use this server.

Ask for help about your paper (PLoS ONE, 2010)

Q:
Currently, I’m working on reconstructing gene regulatory network. It’s
really an interesting topic and I would like to estimate which tools is
suitable for our experimental data. I have read your published paper
"Improved Reconstruction of In Silico Gene Regulatory Networks by
Integrating Knockout and Perturbation
Data". In this paper, I can’t understand the section of learning noise from
deletion data.
Step 1: Calculate the probability of regulation Pb->a for each pair of genes
(b,a). I want to know how to calculate this probability, and this value of
probability can decide potential regulation or not?
I want to ask you that how to work in this section, and I’m appreciated if
you can help me to figure out.

A: Basically we used the expression levels currently believed to be
unaffected by a deletion to form a Gaussian background. Then if a gene
has an expression level far away from the mean of this Gaussian
distribution (by calculating the probability that the expression is as
extreme or more extreme than the observed one based on the Gaussian), we
consider the gene to be affected by the deletion.

coevolution webserver

Q:

I have read the article “An integrated system for studying residue coevolution in proteins” and start to use the coevolution webserver http://coevolution.gersteinlab.org/coevolution/. The help page says SCA is one of the coevolution score function, which I can‘t find on the webpage. Could you please tell me what is wrong?

A:
SCA is only available in the version for download, but not at the Web site anymore since it sometimes ran for an extremely long time. If you want to use the SCA method, you can use the download version or the latest SCA method by the Ranganathan group, which I think is better than the version we implemented.