Hot answers tagged bioinformatics
You saw my last answer of your question? you can do that easily, just use prepare_ligand4.py -l my.mol2 -A "hydrogens"
There is no perfect cut-off. It always depends on what you're doing. The e-value is basically a measure of how many such alignments you would expect to find in a database this size by chance. Therefore, e-values greater than 1 mean that you'd expect at least one alignment similar to what you've found by chance alone. As others have stated, the e-value is ...
The e-value is supposed to be a metric for the chance that an alignment could occur at random, but it is a crude estimate. As pointed out in other answers, this significantly does not include the length of the query sequence. It also does not include the conservation of the gene or the frequency of amino acids (in protein blasts). It does take into ...
E-value refers to the expected number of random hits for a given alignment score. Smaller it is more reliable is your match. There is no hard and fast rule for e-value cutoff. You can keep whatever you want depending on the level of stringency that you require. But you should note that for smaller sequences (< 30nt) there is always a higher likelihood of ...
http://ocw.mit.edu/courses/health-sciences-and-technology/hst-508-genomics-and-computational-biology-fall-2002/audio-lectures/ Download all the slides if you want good resources.
You can use Autodock Vina. It provides an option to calculate local score only. displaying the individual contributions to the intermolecular score, before weighting (these are shown with "--score_only") Autodock is better tool in speed and accuracy than Autodock itself. For using it, you need : 1) Protein.pdbqt 2) Ligand.pdbqt 3) Config.txt (the ...
Your problem will finally boil down to searching your sequence in the Blast databases. Performing Blast seems to be probably the best way to find out if your bacteria has that specific protein expressed or not. If you could not find it in the nearest species using Blast, then try running PSI-BLAST, which would return you distant homologs, by which you can ...
for biomart goto below link http://central.biomart.org/converter/#!/ID_converter/gene_ensembl_config_2 Also there is one more converter which i found pretty useful http://biodbnet.abcc.ncifcrf.gov/db/db2db.php#biodb
Have a look at these ChIP-seq data for HP1 in Drosophila: 1, 2 and 3. From ChIP-seq data you can find the distance between the TFBS peaks and the TSS of the gene. You can also look for nucleosome positioning and DNAse hypersensitvity regions; for the former, I am sure that data is available for Drosophila.
If you're working with BioCyc pathways you can use their REST-API for batch-downloading all genes for their pathways. You can run a wide range of queries including BioVelo-queries for all genes/compounds in a specified pathway, all pathways in organism etc. A query for all the pathways in B. subtilis would look like: ...
A couple of months ago I listened to a plant physiologist who strongly recommended recombination-based mapping over sequencing based mapping. The main reason he gave was the error-rate of 2nd generation sequencing. The error rate on Illumina platforms are about 1 % if I recall. In a small genome like that of Arabidopsis thaliana (157 Mbp) that accounts for ...
Only top voted, non community-wiki answers of a minimum length are eligible