kacst-bioinfo-lab / labwork

3 stars 1 forks source link

Arabian Came alpha-1-microglobulin/bikunin precursor (AMBP), mRNA #16

Open SULTAN-ALHARBI opened 6 years ago

SULTAN-ALHARBI commented 6 years ago

All the best

MManee commented 6 years ago

For more information about how you can perform and generate your results, please check this issue #14

1sar2003 commented 6 years ago

Dear all
I used ncbi ORF finder to determined AMBP ORF with Search Parameters below, and I end up having 14 ORFs, should I change minimal length or codon to have less result or keep it stander as it is ?

screen shot 2017-08-14 at 1 28 09 pm
SULTAN-ALHARBI commented 6 years ago

Hi @1sar2003
Open Reading Frame (ORF) is the region in your mRNA which will be translated to protein, almost all proteins are started with Met(methionine ) amino acid this amino acid has ATG(AUG) codon. so we should use "ORF start codon to use : "ATG" only this is first. second: if you search for AMBP protein length for any well-annotated species, you will obtain that the length of AMBP protein should be around 350 amino acids. For example Human AMBP protein length is 352 amino acids and Caw AMBP protein is 334 amino acids. Our Arabian Camel AMBP protein length will be around those lengths. You should end up with only one ORF

MManee commented 6 years ago
malsaeed commented 6 years ago

@MManee @1sar2003

We already have small call about the work today

And I'm glad to help any person with geneious

We will see if we can meet tomorrow

Regards

1sar2003 commented 6 years ago

thank you all for the respond the license has been activated the maximum number of times so I will use the trial version
@MManee @malsaeed @SULTAN-ALHARBI

SULTAN-ALHARBI commented 6 years ago

That is weird :( Because we have 5 license users and we just used one so far. I'll check with Geneious team what the problem.

1sar2003 commented 6 years ago

Hello
I used geneious to perform my task, the length of Arabian Camel AMBP protein is 352 aa end up with 1 ORF(length 1059) and I followed the next :

1/Translate AMBP mRNA to Protein (BLAST):

screen shot 2017-08-15 at 10 56 19 pm

2/multiple sequence alignment :

screen shot 2017-08-15 at 11 36 58 pm

3/ phylogenetic relationships (Tree)

screen shot 2017-08-16 at 12 39 32 am

I will continue to update the result soon, just need to double check

SULTAN-ALHARBI commented 6 years ago

Great to see a quick movement let's take your work one by one 1/Translate AMBP mRNA to Protein: yes, Arabian Camel AMBP protein is 352 aa. 2/multiple sequence alignment : you can do this task in two ways A) using Arabian Camel AMBP mRNA (1059 bp) and this is what you already did. Here you are comparing Arabian Camel AMBP mRNA with mRNA from other species. B) using Arabian Camel AMBP protein 352 aa and comparing it with AMBP proteins from other species. and this what we actually want all published articles use this method. so make a multiple protein sequence alignment. 3/ phylogenetic relationships (Tree): it is a good tree, we can see our Arabian Camel grouped with Bactrian camel (Camelus bactrianus) and alpaca (Vicugna pacos). but we need to construct our phylogenetic tree using protein sequence alignment

the the best

1sar2003 commented 6 years ago

Dear @SULTAN-ALHARBI how should I obtain the Arabian Camel AMBP protein seq alignment in geneious? is it by select ORF and blast it in the protein data bank?

1sar2003 commented 6 years ago

like when I blast the trans protein under the Nucleotide collection I end up have a list of organisms does not include my dromedaries Camelus, and the same thing applied in most of Blast database

screen shot 2017-08-16 at 2 11 58 pm

Tree

screen shot 2017-08-16 at 2 12 51 pm

By the way, I obtained the translated protein by click Sequence / Translate

SULTAN-ALHARBI commented 6 years ago

Hi @1sar2003 First/ Select the region in the sequence that you want to translate screen shot 2017-08-16 at 4 26 23 pm Second/ go to icon called sequence ---> translate . You will get your translated sequence as separate file. you this seq in multiple alignment then phylogenetic tree screen shot 2017-08-16 at 4 27 38 pm

1sar2003 commented 6 years ago

@SULTAN-ALHARBI this is exactly what I did before, the only thing that I change now is running the blast program with (tblastn) instant of (blastp) and I get what I need

1/Multiple protein sequences alignment:

screen shot 2017-08-17 at 10 31 50 am

2/phylogenetic relationships (Tree):

screen shot 2017-08-17 at 10 37 02 am

3/prediction of the secondary AMBP protein structure: 8ba26956-8328-11e7-9477-00163e110593.AMBP.psipred.pdf

4/prediction of the 3D AMBP protein structure:

screen shot 2017-08-17 at 12 18 57 pm
1sar2003 commented 6 years ago

5/ Isoelectric point for (AMBP):

isoelectric point for ambp

6/ physical&chemical properties of (AMBP) by using Kolaskar & Tongaonkar method only:

physical chemical properties of ambp by kolaskar tongaonkar

7/Prediction ordered and disordered regions of (AMBP) protein : Prediction ordered and disordered regions (AMBP) alpha-1-microglobulin:bikunin precursor.pdf

8/Conserved domains in (AMBP) protien seq:

conserved domains on ambp protien seq

Have a nice weekend.

MManee commented 6 years ago
1sar2003 commented 6 years ago

Dear all @MManee @SULTAN-ALHARBI
I'm trying to reconstruct multiple alignments and phylogenetic tree and I facing some of the species didn't start with M !! is anyone can explain why? and most importantly how to resolve it? *Most of the feature on geneious is restricted

finel

finl tree

14 go to go seq

14 good to go

SULTAN-ALHARBI commented 6 years ago

Dear @1sar2003
Great what you did, it seems to me that you do not have any problem with making multiple alignment and constructing a phylogenetic tree. If you find any species that doesn't fit with your exceptions just remove it. ( see the one in LaTex )
I think you should move in to 2nd and 3D structure and domains of AMBP protein . we should write and read more about AMBP structure

SULTAN-ALHARBI commented 6 years ago

Dear @1sar2003 @malsaeed @MManee This is the new License Key for Genious

Licensee: SULTAN ALHARBI License Key: 1BDA-15EE-1B52-92E1-C6A1 Valid from 03 October 2017 to 03 October 2018

1sar2003 commented 6 years ago

Dear @SULTAN-ALHARBI thank you for your comment, I think I already write about the domains region and regard 3D structure I admit I having a difficulty of understanding the basic of it and will continue trying my best with it.

screen shot 2017-10-05 at 11 42 15 am

best regard

SULTAN-ALHARBI commented 6 years ago

Dear @1sar2003 The information of domins'A1M is very good I will work on it these days. we need one more step to confirm them as 3D structures. I know, it's easy to annotate novel 3D. we will be side to side helping each other

1sar2003 commented 6 years ago

Dear @MManee @SULTAN-ALHARBI In the writing process, I have some thoughts I would like your consultant in it, should I mention the biosynthesis of the two proteins in the cell and the path they follow to express.Second which category the protein inhibition of the protienes falls? is it in discussion or introduction section?

SULTAN-ALHARBI commented 6 years ago

@1sar2003 Hi Sara As most literatures mention those two Protein (A1M and Bikunin are co-synthesis = coming from same mRNA) our scope in this paper is to annotate the structural properties of those proteins as well as their biophysical features. we can briefly mention their biosynthesis pathways just to support any structural or biophysical characterizations. For example\ the section where we want to explain the trypsin interaction site or glycosylation site on the protein. about where to explain such information either on introduction or discussion, it depends on the coordination of the paper.

All the Best

MManee commented 5 years ago