Naúrar 2 / 11

Binciken Jeri na DNA/RNA: Daidaitawa, Motif, Buɗe Tsarin Karatu da Basic Bioinformatics

Riba:

  • Fahimtar ra'ayoyin jeri-jeri, bincike dalili da buɗaɗɗen firam ɗin karantawa (ORF) kuma sami ƙwarewar wucin gadi ta samar da abin aiwatarwa, lambar tantancewar Biopython/nazari.
  • Ikon duba firam, strand da genome zato a cikin jeri da lambar da aka samar ta hanyar basirar wucin gadi da kwatanta sakamakon tare da sanannen tunani.
  • Ability don amfani da horo na rashin amfani da kowane jerin da aka bayar ta hanyar hankali na wucin gadi daga kai da kuma tabbatar da kowane jeri daga tushen farko kamar NCBI / Ensembl

Binciken jeri shine mafi mahimmancin aiki na ilimin halitta: karanta madaidaicin DNA (ko U a cikin RNA) wanda ya ƙunshi haruffa A, T, G, C, kwatanta shi, da gano yankuna masu ma'ana (genes, motifs, jerin tsari) a cikinsa. A cikin wannan rukunin, zaku koyi yadda ake amfani da hankali na wucin gadi (AI) azaman lambar rubutawa da abokin fassara a cikin waɗannan ayyukan; amma zaku koyi dalilin da yasa yakamata ku tabbatar da sakamakon koyaushe tare da lambar aiwatarwa da tushen farko. Babban kayan aikin mu zai zama Biopython (laburaren Python da aka rubuta don aiki tare da jerin halittu) da kayan aikin daidaitawa na hukuma.

Farko gargaɗi: LLM na iya haifar da kurakurai daga ƙwaƙwalwar ajiya, ko da ɗan gajeren jerin. Yana iya haɗa harafi lokacin da aka nemi yin lissafin juzu'i na gaba-gaba. Don haka, kar a taɓa yin ayyukan kirtani ta hanyar dogaro da martanin rubutu na AI, amma tare da lambar da AI ke rubutawa kuma kuna gudana.

Basic Concepts: abin da muke aiki da?

  • Biyu na tushe (bp): Naúrar harafin DNA. Halin halittar ɗan adam ya kai kusan biliyan 3.2 bp.
  • Strand: DNA helix biyu ne; Zaɓuɓɓukan biyu sune anticomplement na juna. Yana da mahimmanci a cikin wane zaren zaren da aka bayyana bambancin.
  • Codon: Rukunin tushe guda uku; kowane codon yayi daidai da amino acid (tushe na gina jiki). Misali, ATG yawanci shine farkon codon (methionine).
  • Buɗe firam ɗin karantawa (ORF): Yankin jeri wanda zai iya yin lamba don furotin, wanda ya tashi daga farkon codon zuwa codon tasha (TAA, TAG, TGA).
  • Motif: Maimaita gajeriyar tsarin tsari wanda ke da takamaiman aiki; misali, yankin da ma'aunin rubutu ya ɗaure.
  • Daidaitawa: Shirya jeri biyu ko fiye ɗaya ƙarƙashin ɗayan don ganin kamanninsu.

Mataki-mataki: aikin bincike na jeri

1. Samo jerin daga tushen abin dogara. Kada ka sanya AI ta ce "tunatarwa" jerin; Zazzage shi azaman FASTA (tsarin rubutu na yau da kullun wanda ke adana jerin abubuwa) daga tushe kamar NCBI, Ensembl, da sauransu kuma ba da wannan jeri ga AI.

2. Yi ciniki tare da code. Yi ayyuka kamar juzu'i, rubutun (DNA → RNA), fassarar (RNA→protein), rabon GC wanda lambar Biopython yayi kuma gudanar da lambar da kanku.

3. Duba tsarin da zaren zato. Tambaye shi/ta ya bayyana a sarari a layin sharhi wanda zaren da kuma wanne tsarin karatun lambar ke gudana.

4. Kwatanta sakamakon tare da sanannun tunani. Daidaita furotin ko ORF da kuka samar tare da sanannen rikodin a cikin bayanan. Tsawo da rashin daidaituwa na farko suna kama mafi yawan kurakurai.

5. Tabbatar da daidaitawa tare da kayan aiki na hukuma. Kada ka bari AI "kwallon ido" kamanceceniya na jerin biyu; Sami maki na lamba tare da BLAST (kayan aikin neman kamanni na jere) ko ɗakin karatu na daidaitawa.

Tukwici: Koyaushe bari tsawon kirtani ya zama abin dubawa na farko. Adadin amino acid na furotin kusan kashi ɗaya bisa uku na adadin sansanonin jeri na coding (ban da codon tasha). Idan tsayin bai dace ba, firam ko zaren ba daidai ba ne.

uku mini lokuta

Case 1 - Kuskuren daidaitawa mai juzu'i. Wani ɗalibi ya tambayi AI game da juyar da madaidaicin jerin 5'-GATTACA-3'; AI ya ba "TGTAATC" (daidai). Koyaya, a cikin jerin mafi tsayi na sansanonin 20, AI ta tsallake tushe kuma sakamakon shine tushe 19. Lokacin da ɗalibin ya gudu tare da Seq ("...") .reverse_complement () a cikin Biopython, ya ɗauki tushe 20 kuma ya kama kuskuren. Lokacin asarar: mintuna 2.

Case 2 - Canjin Frame. Wani mai bincike yana da jerin lambar tushe 900 da aka fassara zuwa furotin; AI "karanta" sunadaran amino acid 280 ta rubutu. Abin da ake tsammani shine amino acid 299 (900/3 - 1 tasha). Bambanci shine cewa AI ya fara ne daga nucleotide na biyu. An sami daidai tsayi lokacin da aka fara lambar daga firam na farko.

Case 3 - Tabbatar da samu. Wani masanin dakin gwaje-gwaje ya bincika jerin 16S rRNA na nau'ikan kwayoyin cuta guda biyu ta hanyar tambaya "sun kasance iri ɗaya?" sai ya tambayi AI; "Mafi yiwuwa iri ɗaya," in ji AI. Lokacin da mai fasaha ya gudanar da BLAST, ya ga kamanni 97.8% da bambance-bambancen tushe guda 12 - wani muhimmin bambanci ga nuna bambancin jinsi. Idan babu maki na lambobi, za a shigar da sakamakon "daidai" kuskure a cikin rahoton.

Misali: rafin Biopython mai tabbatarwa

daga Bio.Seq shigo da Seq # Shigo da jerin daga FASTA da kuka zazzage daga NCBI; Kar a sa AI ta ce "tunani da ni". dna = Seq("ATGGCCATTGTAATGGGCCGCTGAAAGGGTGCCCGATAG") bugu ("Length (bp):", len (dna)) buga ("GC rate (%):", zagaye (100 * (dna.count("G") + dna.count("C")) / len(dna), 1)) buga (#computer) Fassara daga firam 1; har zuwa dakatar da codonprotein = dna.translate(to_stop=Gaskiya) bugu ("Protein:", protein, "| Tsawon (mm):", len (protein))

Ko da yake AI yana rubuta wannan lambar, kuna ganin daidaiton fitarwa ta hanyar gudanar da shi. Tsawon, rabon GC, da furotin suna kwatankwacin abin da aka sani.

Samfura huɗu masu kwafi

1) Tabbatar da tsararrun aiki:

Rubuta lambar Biopython mai aiwatarwa don jerin FASTA mai zuwa: [jeri/aiki]. Yi ƙididdige tsayin, rabon GC, jujjuya ma'amala, da fassarar daga firam 1. Yi sharhi akan wane zaren da firam ɗin aka ɗauka. Kada ku samar da jerin; Kawai yi amfani da jerin da na ba ku.

2) Binciken ORF:

Rubuta lambar Python wanda ke nemo duk buɗaɗɗen firam ɗin karantawa a cikin jerin da aka bayar (da duka ukun akan zaɓin juyawa na zaɓi). Ba da rahoton matsayin farawa, tsayi, da furotin da aka fassara don kowane ORF. Hakanan alamar ORF mafi tsayi.

3) Tabbatar da daidaitawa:

Ina so in kwatanta tsararraki biyu. "Kamar haka?" Kada ku yi hukunci da ido; rubuta lambar daidaitawa ta hanyar hanya biyu kuma bayar da rahoton adadin kamanni da lambar bambanci a lamba. Tushen: [jeri na 1], [jeri na 2].

4) Binciken Motif:

Bincika abin da ke gaba (kuma azaman magana ta yau da kullun) a cikin kirtani da aka bayar: [motif]. Lissafa wurin (mai tushe 1) na duk matches. Rubuta kuma saka matches masu juna biyu kuma.

Rauni mai ƙarfi / Ƙarfi mai ƙarfi

Rarrauna: "Rubuta furotin na wannan jerin: ATGGCC..."

Matsala: AI tana fassara da rubutu, na iya rikitar da firam/zaren, ba zai iya tabbatar da tsayi ba.

Ƙarfafa: "Rubuta lambar Biopython mai aiwatarwa wanda ke fassara jeri mai zuwa daga firam 1, ya ba da rahoton tsayi da dakatar da codon; kar a canza jerin, kawai yi amfani da wanda na ba: ATGGCC..."

Me ya sa yake da ƙarfi: Ana yin aiki a lamba, tsarin a bayyane yake, ana iya tantance fitarwa ta lambobi.

nema

hanya mara kyau

daidai tsarin

baya cikawa

Bari AI ta rubuta da rubutu

Biopython reverse_complement()

fassarar

Bari AI ta fassara daga ƙwaƙwalwar ajiya

Code, ƙayyadaddun tsarin

kamanceceniya

"Kamar haka?" shawarar ido

BLAST/makin daidaitawa

dalili

Bari AI ta ƙidaya da hannu

Lambar, tare da lissafin wuri

tushen tsararru

Bari AI tuna

FASTA daga NCBI/Ensembl

Kuskuren gama gari

  • Ba ƙayyade tsarin ba. Fassara daga firam ɗin da ba daidai ba yana haifar da gajeriyar furotin ko kuskure.
  • Tangling da yarn. Bambancin ko motsi na iya kasancewa a cikin zaren baya; zato zaren ya kamata a rubuta.
  • Yin AI ya haddace jerin. LLM ba zai iya samar da dogon kirtani ba tare da kurakurai ba; Kullum kuna ba da jerin abubuwa.
  • Yin hukunci da ido don kamanni. Kar a ce "iri daya/kama" ba tare da maki na lamba ba.
  • RNA/DNA cakuda. Hada U tare da T yana rushe fassarar; bayyana nau'in shigarwar.
Tsanaki: Ko da babban kaso na kamanni a cikin BLAST da makamantan kayan aikin ba lallai ba ne yana nufin “kamance” a ilimin halitta; Ya kamata a kimanta ƙimar e-daraja (yiwuwar damar) da tsawon yankin da aka haɗa tare.

Zurfi: karanta fitowar BLAST daidai

Samun AI fassara sakamakon BLAST yana adana lokaci; Amma kada ku yanke shawara har sai kun karanta batutuwa uku da kanku. Na farko shine e-darajar (darajar da ake tsammani): adadin da ake tsammanin wannan makin zai iya faruwa kwatsam; Ƙananan ƙima kamar 1e-50 yana nufin ƙarfi, ƙimar kamar 0.1 kusan amo. Na biyu shine kewayon tambaya: kashi nawa ne na jerin tambayoyin wasan ya rufe; 98% kamance amma 20% kawai ɗaukar hoto yana nufin cewa ƙaramin sashi na jerin yayi kama da yaudara. Na uku shine ainihin kashi. Idan ba tare da waɗannan ukun sun karanta tare ba, yawancin kashi ɗaya kawai ba ya tabbatar da komai.

Misalin tabbatacciyar misali: wani mai bincike ya WUTA wani guntun kwayar halittar da ya tsara; "99% daidai da ɗan adam BRCA2, gene iri ɗaya," in ji AI. Lokacin da mai binciken ya kalli fitowar, ya ga cewa ɗaukar hoto shine kawai 15% - ɓangaren madaidaicin ɗan gajeren yanki ne kawai, maimaita yanki daga dubban sansanonin BRCA2. Madaidaicin fassarar ba "jinin halitta iri ɗaya ba" amma "yana raba maƙasudin maimaitawa gama gari". Karatun iyakar ya hana cikakkiyar ganewa.

Rukunin BLAST

me yake cewa

tarko

E-daraja

yuwuwar daidaituwa

Idan yana da girma, wasan na iya zama mara ma'ana

Tambayar ɗaukar hoto

Adadin tambayar da aka rufe

Idan yana da ƙasa, kashi yana ɓatarwa

kashi na ainihi

Madaidaicin ƙimar tushe

kadai bai isa ba

bitscore

Ƙarfin daidaitawa daidai

Fassara bisa ga tsayi

5) Samfurin fassarar fitarwa na BLAST:

Fassara tebur BLAST mai zuwa, amma kar ku yanke shawara: taƙaita ƙimar e-daraja, ɗaukar hoto da ainihin kashi na kowane jere daban kuma nuna abubuwan da ake buƙatar cikawa kafin cimma matsaya kamar "jinin halitta iri ɗaya". Table: [manna].

A takaice

  • Ayyukan jeri (juya baya, fassarar, ORF, rabon GC) yakamata a yi tare da lambar da AI ta rubuta kuma kuna gudana, ba tare da martanin rubutu na AI ba.
  • Tsarin tsari da zato ya kamata koyaushe a bayyana a sarari; duba tsawon shine kayan aikin kama kuskure mafi sauri.
  • Koyaushe samun jerin daga tushen abin dogara (NCBI, Ensembl); Kar a sanya AI ya haddace shi.
  • Ana kimanta kamanni da jeri ta kayan aikin hukuma da makin lambobi, ba ta ido ba.

Aikin aikace-aikace

Zazzage ɗan gajeren jerin lambobin ƙididdigewa daga ingantaccen tushe (misali NCBI). Tare da samfuran 1 da 2 a sama, tambayi AI don lambar Biopython, gudanar da lambar; Kwatanta tsayi da jerin furotin da kuka samar tare da sanannen rikodin a cikin bayanan. Idan kun sami rashin daidaituwa, gwada gyara ta ta canza tsarin tsarin / zato kuma ku lura da tsarin.

jerin abubuwan dubawa

  • [ ] Na sami jerin daga madaidaicin tushe, ba ni da AI ta haddace shi.
  • [ ] Na yi ayyukan tsararru tare da lambar aiwatarwa.
  • [ ] Na fayyace ƙayyadaddun tsari da zato.
  • [] Na kwatanta tsayin furotin/ORF tare da tunani.
  • [ ] Na kimanta kamancen tare da kayan aiki na hukuma da maki na lambobi.
  • [ ] Na duba RNA/DNA da U/T rabuwa.