Bioinformatics I (BIF401) book HANDOUTS / POWER POINT SLIDES in pdf

<br /><b><span style="font-size: large;">Bioinformatics I (BIF401) book HANDOUTS / POWER POINT SLIDES in pdf</span></b><br />Handouts | Lectures | Contents | Books<div> <br /></div> <div> <a href="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhPpZpS5uOo-fjphx1nsKKD-ET3IT50S_kDfrbQ4rCEsjpnISbH5gFJrSDdsFOZ_-NbYRzRCdguGpw3MLIms4DkH6zwGb4oHbXM0BkAp6GPDibO96NaS52doLMsOgw6XGpI5NlUs51ry6Q/s1600/Capture.PNG" imageanchor="1"><img border="0" src="https://blogger.googleusercontent.com/img/b/R29vZ2xl/AVvXsEhPpZpS5uOo-fjphx1nsKKD-ET3IT50S_kDfrbQ4rCEsjpnISbH5gFJrSDdsFOZ_-NbYRzRCdguGpw3MLIms4DkH6zwGb4oHbXM0BkAp6GPDibO96NaS52doLMsOgw6XGpI5NlUs51ry6Q/s1600/Capture.PNG" /></a></div> <div> <div> Chapter 1- Introduction to Bioinformatics</div> <div> Module 001: INTRODUCTION TO BIOINFORMATICS</div> <div> <div> AN INTRODUCTORY COURSE</div> <div> BIOINFORMATICS-I</div> <div> A STUDENT HANDOUT</div> <div> 2016</div> <div> VER 1.0.0.0</div> <div> LAST REVISED ON May 2, 2016</div> <div> Table of Contents</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 1 of 320</div> <div> Chapter 1- Introduction to Bioinformatics</div> <div> Module 001: INTRODUCTION TO BIOINFORMATICS</div> <div>  BACKGROUND</div> <div> Bioinformatics is an interdisciplinary science at the cross-roads of biology, mathematics, computer science, chemistry and physics. With the digitalization of the biological information, doors have been wide opened towards the analysis of this information using computer algorithms and software.</div> <div> Now we know well that the human genome has over 25,000 genes and these genes code for thousands of different proteins which perform day-to-day functions in the living cell. Furthermore, these proteins may take on various post-translational modifications leading to a very large number of functionally unique molecules. This presents us with a huge challenge in identification of genes and proteins.</div> <div>  EXPERIMENTS IN BIOLOGY</div> <div> With the advancements in experimental protocols, now we have several next generation instruments and techniques available for obtaining digitalized biological information on genes and proteins etc. These instruments include:</div> <div> 1. Next Generation Sequencers (NGS) for whole genome sequencing</div> <div> 2. High Resolution Mass Spectrometry for whole proteome profiling</div> <div> 3. Nuclear Magnetic Resonance Spectroscopy for structural studies</div> <div>  DIGITALIZATION OF BIOLOGY</div> <div> In today’s world, when a biologist performs an experiment in the wet-lab, he or she in fact produces digital data which is continuously being stored on computer disks. The data may include text, numbers, symbols or images.</div> <div>  SPEED OF DATA GROWTH</div> <div> Due to advancement in instrumentation used in biological experiments, data is being accumulated at exponentially increasing rates. For example; genome sequences in genome databases are doubling every few years.</div> <div>  CONCLUSION</div> <div> Human brain is limited in recalling information from memory. First, we have to commit all information to our memory followed by its recall. To overcome our ability to memorize and recall, computers can come to our rescue. This is because computers have an infinite ability to recall this information and process it quickly towards results.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 2 of 320</div> <div> Module 002: INTRODUCTION TO BIOINFORMATICS</div> <div>  MOTIVATION</div> <div> Bioinformatics is a becoming a popular science due to several reasons.</div> <div>  It is an interdisciplinary field as it covers the information of biological digital information including human, plants, animals and microorganisms.</div> <div>  Although it is a new field but it is rapidly developing field.</div> <div>  It demands a very low cost infrastructure and hardly any lab equipment.</div> <div>  As bioinformatics data concerns a wide range of species such as humans, plants and micro-organisms, it presents us with plenty of opportunities in scientific discovery.</div> <div>  SCOPE OF BIOINFORMATICS</div> <div> Bioinformatics primarily deals with digitalized biological information as well as data reported from biology experiments. Computational methods, data processing techniques and algorithms are employed in addressing the following issues:</div> <div>  Storage of data</div> <div>  Organization data</div> <div>  Analysis of many experiments</div> <div>  For representation of biological information</div> <div>  ACTIVITIES</div> <div> In modern biological sciences, bioinformatics is used for activities such as:</div> <div>  Developing algorithms for organizing data collected from experiments</div> <div>  Writing software and tools for data analysis</div> <div>  Data processing to determine the role of underlying biomolecules</div> <div>  Statistical evaluation of data using methods such as t-test and ANOVA</div> <div>  Data visualization for meaningful presentation of biological information</div> <div>  CONCLUSION</div> <div> In Pakistan, the field of biology is undergoing a rapid change due to the onset of bioinformatics. New research and educational programs are being constructed which is opening new door of opportunities for our future generations.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 3 of 320</div> <div> Module 003: INTRODUCTION TO BIOINFORMATICS</div> <div>  NEED FOR BIOINFORMATICS –I</div> <div> If we look at the pace of development in the area of bioinformatics then we can easily observe that from year’s 2000 to 2015, the number of online tools for processing genomics and proteomics information are rapidly increasing. This is just a reflection of the need for bioinformatics in modern day biology. The field of Bioinformatics and Computational Biology is characterized by a highly diverse confluence of traditional academic disciplines. Informatics and Bio-science are the umbrella terms given to a set of allied disciplines which make up the field, but a much larger array of traditional areas contribute to the set of tools needed by individuals training for this new and expanding interdisciplinary field. Biomedical Engineering, Electrical and Computer Engineering, Computer Science, Applied Mathematics, Genetics, Biology, Anatomy and Cell Biology, Micro Biology, and Biostatistics are the principal allied disciplines.</div> <div>  CONCLUSION</div> <div> The need for bioinformatics is on a rapid rise as biological data is rapidly increasing and becoming available online, free of any cost.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 4 of 320</div> <div> Module 004: NEED FOR BIOINFORMATICS –II</div> <div> If we observe the growth of gene bank than from 1982 it comprised of 2 billion base pairs but by year 2002 it had risen to 56 billion base pairs. With the data in our hands, there is an urgent need to interpret this data. For instance, analysis of this data can help us in developing an understanding of the phylogenetic “tree of life” which consist of:</div> <div>  Bacteria</div> <div>  Archaea</div> <div>  Eucarya</div> <div> Towards exploring the possible benefits of using bioinformatics, one needs to answer the following question:</div> <div>  WHAT IS IT THAT BIOINFORMATICS CAN DELEIVER?</div> <div> The simple answer to that bioinformatics is:</div> <div>  Provide us better understanding of life, evolution, molecular mechanisms as well as disease.</div> <div>  Moreover, we can make better drugs with the availability of an enhanced molecular understanding of disease.</div> <div>  POSSIBLE CONTRIBUTIONS</div> <div>  It can help us to organize the large datasets from new experiments instruments</div> <div>  Bioinformatics can help store and process this data as well.</div> <div>  It can provide insights into the meanings of our research results and findings.</div> <div>  Overall, it can help us to better understand paradoxes defining the life forms.</div> <div>  CONCLUSION</div> <div> From gene sequencing to protein sequencing, bioinformatics is providing us with an improved understanding of the genes, proteins, protein interaction and signaling pathways involved in biological functioning and disease.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 5 of 320</div> <div> Module 005: APPLICATIONS OF BIOINFORMATICS – I</div> <div> When we look at bioinformatics, it seems to be a very complex and abstract field. How and where can bioinformatics be applied specifically? How does it improve the fundamental understanding of biological phenomenon? Most importantly, how can its benefits be delivered to the society at large?</div> <div> The answers to these questions are categorized as follows:</div> <div>  GENOMICS</div> <div>  Bioinformatics can help in assembling DNA sequencing data.</div> <div>  It can help in gene finding (markers).</div> <div>  Gene assembly can be performed using bioinformatics tools (nucleotide alignments)</div> <div>  It can help transcribe the gene data to RNA data</div> <div>  Also, databases can be generated from such data.</div> <div>  EVOLUTIONARY STUDIES</div> <div>  Evolutionary relationships between different organisms can be derived from data.</div> <div>  Evolutionary distance among species can be computed by using bioinformatics tools.</div> <div>  Phylogenetic trees can be constructed to find relationships between species.</div> <div>  Ancestry can be better understood between several species and organisms.</div> <div>  PROTEOMICS</div> <div>  Bioinformatics can help us in decoding protein sequences.</div> <div>  It can also help us in understanding protein structure.</div> <div>  We can also understand post translational changes in proteins with the help of bioinformatics.</div> <div>  We can better understand the protein-protein interaction in different biological reactions.</div> <div>  It can also help us in generating databases of these sequences and structures.</div> <div>  SYSTEMS BIOLOGY</div> <div>  Bioinformatics can assist us in modelling regulatory mechanisms in gene and protein networks.</div> <div>  Such models can be analyzed to identify the key regulators in these networks.</div> <div>  Moreover, the models can help evaluate drugs to treat these key regulators.</div> <div>  CONCLUSION</div> <div> Bioinformatics can be applied to life in many ways it helps us to understand the sequence and function of biomolecules and their relationships. Recent trends in bioinformatics involve development of personalized therapeutics for cancer and diabetes.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 6 of 320</div> <div> Module 006: APPLICATIONS OF BIOINFORMATICS - II</div> <div> Bioinformatics is being applied in routine life in many ways like in Genomics, transcriptomics, Proteomics, Metabolomics, Structural Proteomics, Designing Drugs, System Biology and in personalization of medicines for cure.</div> <div> Except these applications Bioinformatics introduced us the techniques which enabled us to generate the large data regarding biology and also its use. And step by step the applications of bioinformatics increased from genomic level to entire system level.</div> <div>  SMALL TO BIG</div> <div>  Bioinformatics helps us to understand the systems from small to big like from gene findings to entire system prediction</div> <div>  In structure findings and modeling of many biological system to understand them in better ways.</div> <div>  Bioinformatics helped the human to understand the protein, protein interaction in many biological systems.</div> <div>  And provide us the concept how these biological process are interconnected with each other and how they affect each other.</div> <div>  Now we are able to understand the modeling of molecules and genome at cell level.</div> <div>  Signaling pathways are easy just because of bioinformatics.</div> <div>  Now morphology of tissue can be understand by creating the models with help of bioinformatics tools.</div> <div>  CONCLUSION</div> <div> Bioinformatics not only just collect, analyze and store the data it process it in very authentic way and validates our hypothesis and very soon in future it will help us to understand that which disease is coming in future and how to tackle it with personalize medicine.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 7 of 320</div> <div> Module 007: FRONTIERS IN BIOINFORMATICS - I</div> <div>  INTROCDUCTION</div> <div> Bioinformatics is new and emerging field of science having vast opportunities and with innovation in tools it is increasing the scale of biological data, but still there are many unsolved challenges which are pending in the field of life science and for which bioinformatics is doing new innovative ideas.</div> <div>  FRONTIER IN GENOMICS</div> <div> Now we are able to sequence the whole genome with the bioinformatics tool of Next generation sequencing (NGS)</div> <div> We are able to save, store and analyze the massive amount of biological data which is in (Terabyte files)</div> <div> We can handle the large number of data easily and can process it as well in easy way.</div> <div> Whole genome can be assemble in sequence and can flaws can be identified easily.</div> <div>  FRONTIER IN TRANSCRIPTOMICS</div> <div> Now in genomics we are able to identify those matters which are unknown yet or under discussion.</div> <div> Role of RNA in making proteins and its dynamics can be understood easily now.</div> <div> Interactions of RNA molecule can be easily understood by simple model.</div> <div>  FRONTIER IN PROTEOMICS</div> <div> Deficiency of low proteins in any patient tissue sample can be identify.</div> <div> Expression and manufacture of protein in large molecular level in any organism can be identified.</div> <div> Pathways before and after any biological reaction are easy to design.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 8 of 320</div> <div>  CONCLUSION</div> <div> Bioinformatics is literally a science full of challenges and opportunities having a revolution in field of biology and routine life.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 9 of 320</div> <div> Module 008: FRONTIER IN BIOINFORMATICS-II</div> <div> Frontier in Bioinformatics includes</div> <div>  Next generation genomics</div> <div>  Transcriptomics</div> <div>  Proteomics</div> <div>  FRONTIER IN PROTEIN STURUCTURE</div> <div> Bioinformatics helps us to understand the layer folding of proteins that how they are processed, and helps to know that how protein interact with each other and how a drug can affect or stimulate a protein.</div> <div>  FRONTIER IN SYSTEM BIOLOGY</div> <div> It helps us to understand the whole system of a single cell, in that cell how organelles, gene, proteins and metabolites are interconnected in a single unified system (cell). And bioinformatics also give us the idea how these models can be applied to real-time.</div> <div>  FRONTIER IN PERSONALIZED MEDICINE</div> <div> This is the important thing for this century and upcoming generation that personalize the medicine for exact cure of a disease. Because all the medicine cannot work exact some effect patient badly therefor with the help of Bioinformatics we are now able to personalize some medicines for some diseases. And bioinformatics helps us to evaluate the medicine.</div> <div>  CONCLUSION</div> <div> If we talk about the 21st century than it’s the century of bioinformatics it will enable the human to cure many disease with one drug by personalizing it.</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 10 of 320</div> <div> Module 009: Overview of Course Contents - I</div> <div> Philosophy behind the Course Outlay</div> <div> 1. Introduce the classical algorithms in bioinformatics</div> <div> 2. Link them to latest developments in the field</div> <div> 3. Evaluate the future applications</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 11 of 320</div> <div> Sequences and operations such as alignment and comparison will be covered along with phylogenetic and RNA structure modelling. Next up we will delve into protein sequences and structures!</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 12 of 320</div> <div> Module 010: Overview of Course Contents - II</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 13 of 320</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 14 of 320</div> <div> Summary</div> <div> • Protein sequence and structure topics will be dealt in these modules</div> <div> • Next set of modules is about the homology modelling and systems biology topics!</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 15 of 320</div> <div> Module 011: Overview of Course Contents – III</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 16 of 320</div> <div> Handout Chapter 01: Introduction to Bioinformatics</div> <div> Page 17 of 320</div> <div> Conclusion</div> <div> • These contents will give you an initial exposure to the variety of topics in bioinformatics</div> <div> • After covering these topics, you should have a basic conceptual foundation for further studies into Bioinformatic</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 18 of 320</div> <div> Chapter 2 - Sequence Analysis</div> <div> Module 001: Gene, mRNA and Protein Sequences</div> <div>  INTRODUTION</div> <div> We all know that all the living things are composed of cells and here a question arise that how cells are made? For composition of cell DNA has blueprints for building cells along with the information of cell’s protein, carbohydrate and vitamins production.</div> <div> And transfer of this information from DNA to these molecules is termed as “Central Dogma” which is</div> <div> DNA RNA Protein.</div> <div> Proteins are than use in constructing the cell.</div> <div>  DNA</div> <div> Figure 0.1 DNA Double helix</div> <div> DNA molecule is double helix structure contain base pairs composed of nucleotides and these nucleotides are composed of sugar phosphate group and are bind with each other with hydrogen bonds.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 19 of 320</div> <div> Normally all the nucleotides are same in both DNA and RNA except one position in RNA which is U (Uracil) and in DNA it is T (Thiamin)</div> <div> DNA sends the information to cell via mRNA and that sequence the amino acids according to coded information and protein structure is formed and that protein form a cell.</div> <div>  CONCLUSION</div> <div> According to the central dogma DNA codes information for RNA and RNA makes the Protein and that protein along with some organelles make cells and its systems.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 20 of 320</div> <div> Module 002: TRANSCRIPTION</div> <div> All cells are made of carbohydrates and proteins and for these cells DNA codes the information which makes the RNA and protein both.</div> <div> Figure 0.2 Flow of information from DNA to Proteins</div> <div> The above mechanism explains the process of transcription in very simple way, DNA codes the information and converted into RNA where mRNA copies the information and it execute the information in cell and amino acids combine with each other according to coded information of DNA and protein formation takes place. Which is known as Translation.</div> <div> Molecule of DNA contains only four base pairs (A, T, C, and G) which are repeated thousands of time and Adenine “A” pairs with Cytosine “C”, While Thymine “T” binds with Guanine “G” and all pairings are with the help of Hydrogen bonding.</div> <div> Same like DNA, the RNA contains four base pairs but Thymine is replaced with Uracil “U” and RNA is single stranded.</div> <div> DNA just codes the information for protein but RNA helps in making protein.</div> <div> DNA</div> <div> RNA</div> <div> Proteins Translation Transcription Information Copy of Information Execution of Information</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 21 of 320</div> <div> Module 003: NUCLEOTIDES</div> <div> If we talk about the composition of DNA and RNA molecule than these are composed of four other molecules which are named as Nucleotides.</div> <div> These molecules are Adenine (A), Cytosine (C), Thymine (T), Uracil (U), and Guanine (G).</div> <div> DNA molecule although is double stranded and RNA is single stranded but there is difference in sugar composition.</div> <div> RNA has Ribose sugar and DNA has de-oxyribose sugar:</div> <div> Figure 0.3 Difference between RNA and DNA sugar</div> <div> RNA DNA</div> <div> Adenine and Guanine collectively called Purines while Cytosine, Uracil, and Thymine are called as Pyrimidine.</div> <div> when phosphate, nitrogen base and sugar come together if there is (OH) than molecule is RNA and if there is (H) in sugar than molecule is DNA. As figure shows.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 22 of 320</div> <div>  CONCLUSION</div> <div> DNA molecule make RNA and RNA make the protein and DNA differ from RNA in nature due to sugar and nucleotide.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 23 of 320</div> <div> Module 004: TRANSLATION</div> <div> Cells are built of proteins and carbohydrates and these proteins are made in results of transformation of RNA molecule and this transformation is called as translation.</div> <div> Translation takes place in ribosome of cell and ribosomes after reading the information of mRNA collects the amino acids from cell cytosol which is the part of the cytoplasm that is not held by any of the organelles in the cell.</div> <div>  MECHANISM</div> <div> At ribosome three nucleotides are read at a time from mRNA, this set of three nucleotide is called as codon and each codon correspond to a specific amino acid.</div> <div> Figure 0.4 sixty four codons combinations</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 24 of 320</div> <div>  CONCLUSION</div> <div> RNA codes for protein and codons of here nucleotide code for specific amino acid on ribosomes and this process is called as translation.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 25 of 320</div> <div> Module 005: AMINO ACIDS</div> <div> RNA decodes the information at ribosomes in form of Codons each codon select a specific amino acid. Because there are 20 different amino acids in nature therefore they fold together and make a protein structure by polymerizing themselves.</div> <div> If we observe the structure of amino acid it contains nitrogen, hydrogen, oxygen and two carbon atoms. And a variable group R.</div> <div> Figure 0.5 structure of amino acid</div> <div> When polymerizations takes place water is formed and if any compound attached with R group than structure of protein is changed.</div> <div> These amino acids are joined with each other with peptide bonds and fold with each other in 3D form they make protein structure.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 26 of 320</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 27 of 320</div> <div> Module 006: STORAGE OF BIOLOGICAL SEQUENCE INFORMATION</div> <div> We know that sequence of DNA contain A,C,T&G nucleotides and sequence of RNA contains A,C,U&G while sequence of protein contain A,R,N,D,C,E,Q,G,H,I,L,K,M,F,P,S,T,W,Y&P these are actually 20 different amino acids in nature which compose a protein.</div> <div> When both DNA and RNA or mRNA are sequenced in lab their sequences contains larger number of nucleotides with variety</div> <div> And when we talk about protein its sequences contain large number of bases as they are complex in nature.</div> <div>  SOLUTIONS DATABASES</div> <div> This large number of sequence or bases cannot be stored in a single computer that’s why solution lies in public sequence data bases for DNA & RNA the public database is GenBank (by NIH).</div> <div> For proteins the public database is UniProt (by Uniprot Consortium)</div> <div> Both GenBank and UniProt are online database and the DNA, RNA and Protein sequences are available here online for public and researchers.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 28 of 320</div> <div> Module 007: USING GENBANK</div> <div> GenBank is online database where researcher can get access to the sequences of DNA, RNA and proteins.</div> <div> To find any sequence we go online to NCBI GenBank website which is Public database site. Which is;</div> <div> www.ncbi.nlm.nih.gov/genbank</div> <div> And for example we want to find the sequence for Immunoglobulin which is responsible for Glycoprotein antibodies in white blood cells plasma and act for immunity.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 29 of 320</div> <div> Sequences can be searched from GenBank by typing;</div> <div> o Sequence name</div> <div> o ID</div> <div> o Name</div> <div> o Species</div> <div> o Locus</div> <div> o Accession Number</div> <div> o Author</div> <div> o Journal</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 30 of 320</div> <div> Module 008: USING UNIPROT</div> <div> UniProt is public database which is being used to search the sequence of proteins.</div> <div> www.Uniprot.org</div> <div> For example we want to search a sequence of a protein which is Ubiquitin which plays an important role in cytosol for recycling the proteins. We have to go online to the website www.Uniprot.org and above page will appear.</div> <div> We have to write the name of protein in search box and press enter. You will get the searched results like this one.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 31 of 320</div> <div> By clicking on any result you can download or Blast the sequence.</div> <div> In home page there is a box named “Swiss Prot” which contains human curated protein information, molecular mass, observed and predicted modifications etc.</div> <div> Uniprot can be searched by typing amino acid, Name, ID or sequence.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 32 of 320</div> <div> Module 009: COMPARING SEQUENCES</div> <div> There are millions sequences on GenBank and UniProt what will happen if we will compare them? By comparing sequences of DNA, RNA and Proteins we can get</div> <div>  Similarity among sequences</div> <div>  There might be some specific difference due to some disease or mutation</div> <div>  There may be some evolutionary relationship.</div> <div> As there nucleotides can be similar or differ from each other</div> <div> Figure 0.6 BLAST is used to compare the nucleotides sequences</div> <div> While UniProt is used in case of amino acids sequence comparison.</div> <div> By comparison of nucleotides and Amino acids of any DNA, RNA and protein sequence we can find many evolutionary facts and relations among species.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 33 of 320</div> <div> Module 010: SIMILARITIES & DIFFERENCES IN SEQUENCES</div> <div> When we compare the sequences of DNA and RNA we can get the similarity and differences or relationship in evolution. And same case is with amino acids of proteins.</div> <div> In compression not only they have the same number of nucleotides but they have same order or arrangements.</div> <div> If some sequence are exactly similar to each other it means;</div> <div>  They might have some regular expression in cell or system.</div> <div>  Or they indicate some specific presence like signature of any protein or gene.</div> <div>  Or they might have similar nucleotide just one or two between them are different from rest.</div> <div>  CONCLUSION</div> <div> If there is exact match in sequences it means their order or arrangement and maximum numbers of nucleotides match to each other not all of those. While the genome of each created kind is unique, many animal kinds share some specific types of genes that are generally similar in DNA sequence. When comparing DNA sequences between animal taxa, evolutionary scientists often hand-select the genes that are commonly shared and more similar (conserved), while giving less attention to categories of DNA sequence that are dissimilar. One result of this approach is that comparing the more conserved sequences allows the scientists to include more animal taxa in their analysis, giving a broader data set so they can propose a larger evolutionary tree. Although these types of genes can be easily aligned and compared, the overall approach is biased towards evolution. It also avoids the majority of genes and sequences that would give a better understanding of DNA similarity concepts. http://www.icr.org/article/common-dna-sequences-evidence-evolution/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 34 of 320</div> <div> Module 011: SIMILARITIES & DIFFERENCES IN SEQUENCES</div> <div> When we compare the sequences of DNA and RNA we can get the similarity and differences or relationship in evolution. And same case is with amino acids of proteins.</div> <div> In compression not only they have the same number of nucleotides but they have same order or arrangements.</div> <div> If some sequence are exactly similar to each other it means;</div> <div>  They might have some regular expression in cell or system.</div> <div>  Or they indicate some specific presence like signature of any protein or gene.</div> <div>  Or they might have similar nucleotide just one or two between them are different from rest.</div> <div>  CONCLUSION</div> <div> If there is exact match in sequences it means their order or arrangement and maximum numbers of nucleotides match to each other not all of those. While the genome of each created kind is unique, many animal kinds share some specific types of genes that are generally similar in DNA sequence. When comparing DNA sequences between animal taxa, evolutionary scientists often hand-select the genes that are commonly shared and more similar (conserved), while giving less attention to categories of DNA sequence that are dissimilar. One result of this approach is that comparing the more conserved sequences allows the scientists to include more animal taxa in their analysis, giving a broader data set so they can propose a larger evolutionary tree. Although these types of genes can be easily aligned and compared, the overall approach is biased towards evolution. It also avoids the majority of genes and sequences that would give a better understanding of DNA similarity concepts. http://www.icr.org/article/common-dna-sequences-evidence-evolution/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 35 of 320</div> <div> Module 012: PAIR WISE ALIGNMENT –II</div> <div> In pair wise alignment of nucleotides the nucleotides comes in pairs and matching are colored while missing amino acids are indicated with “” and this empty space is called as gap.</div> <div> And these Gaps are inserted for deletion or insertion of any nucleotide. Increase in Gaps can increase the chance of plenty in sequencing and less number of Gaps can increase the similarity rate of sequences.</div> <div> There are two types of pair alignments.</div> <div> 1. Global</div> <div> 2. Local</div> <div> In Global ways of sequence pair alignment we introduce the Gaps in all sequence to know over all matching. While in Local type of sequence pair alignment we find those regions where nucleotides are maximum matching with each other it is used to find the similarity or some nutation.</div> <div> Most important the Gaps are introduced so that we may add the missing nucleotides. Pairwise Sequence Alignment is used to identify regions of similarity that may indicate functional, structural and/or evolutionary relationships between two biological sequences (protein or nucleic acid).</div> <div> http://www.ebi.ac.uk/Tools/psa/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 36 of 320</div> <div> Module 013: PAIR WISE SEQUENCE ALIGNMENT –III</div> <div> Pair wise alignment helps us to find the similarity and differences there are three ways according to which sequences can differ from each other.</div> <div> Which are</div> <div>  Substitutions ACGA  AGGA</div> <div>  Insertions ACGA  ACCGA</div> <div>  Deletions ACGA  AGA</div> <div> By applying all above ways to any sequence the matching and mismatching can be increased or decreased between to different comparing sequencing.</div> <div> Both local and Global ways of alignments give us different results.</div> <div> But among above Substitution increases mismatch of sequence.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 37 of 320</div> <div> Module 014: DOT PLOTS</div> <div> To visualize the sequence alignment we have a method called Dot Plots in this method the sequence is written top and left side of dot matrix grid.</div> <div> Where one nucleotide or amino acid match with each other the dot is placed in grid position in each row for one time.</div> <div> Similar dots are match with diagonal pattern and which remain separate differ from similar sequence</div> <div> A C A C G A C C G G</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 38 of 320</div> <div> Dots on diagonal repeats the alignments and separate one give difference to the sequence. Figure 0.7 dot plot diagonal pattern</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 39 of 320</div> <div> Module 015: EXEMPLE OF DOT PLOTS</div> <div> In dot plot the matching nucleotides are connected in diagonal way and represent the sequence alignments.</div> <div> When we compare the human Cytochrome and Tuna Fish Cytochrome than the diagonal alignment of sequence we find is in this below diagram.</div> <div>  BENEFITS</div> <div> Dot plots provides us the Global similarity between the two sequences and helps us to visualize the alignments of sequences and sequence repeats appear as diagonal stacks in plot.</div> <div>  CONCLUSION</div> <div> Dot plot help us to find the threshold difference among two sequences.</div> <div> Figure 0.8 tuna fish vs Human</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 40 of 320</div> <div> Module 016: IDENTY VS SIMILARITY</div> <div> When we talk about the comparison of two sequences than question arise that how we can compare the biological sequences and after comparison what will be the degree of comparison.</div> <div> There are two concepts for sequence analysis</div> <div> 1. Identity</div> <div> 2. Similarity</div> <div> Identity means the counting number of nucleotides or amino acids which exactly match when two biological sequences are matched.</div> <div> For example:</div> <div> Number of match = 5</div> <div> Smaller length = 5</div> <div> Sequence (1) = 7</div> <div> Sequence (2) = 5</div> <div> Formula for Identity:</div> <div> Identity = No. of Matches / smaller length × 100</div> <div> =</div> <div> And Similarity means the comparison between two different sequences calculated by alignment approach.</div> <div> In both identity and similarity the dots are not counted.</div> <div> 1: CATGCTT 2: CATGC</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 41 of 320</div> <div> Module 017: INTRODUCTION TO ALIGNMENT APPROACHES</div> <div> When we align the sequence that may be vary due to insertion and deletion of nucleotides and to calculate the similarity we need to align the sequence first. And there are two different approaches to align the sequence.</div> <div> 1. Global Alignment</div> <div> 2. Local Alignment</div> <div> In Local alignment we compare one whole sequence with the one portion of other like this.</div> <div> While in Global alignment we compare both sequence from end to end completely.</div> <div> Local alignment just focus on highly matching portions of sequence while in Global one whole sequence is compared with other one.</div> <div> Figure 0.9 local alignment Figure 0.10 Global alignment</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 42 of 320</div> <div> Module 018: WHY LOCAL ALIGNMENT</div> <div> When there is Global alignment which compare the whole sequence from end to end than why local alignment is done question arise.</div> <div> Because Local alignment have power to detect the smaller regions with high similarity and such matches are motifs or domains which remain hidden in case of protein function.</div> <div>  DOMAIN SHUFFLING</div> <div> Aligned portions of sequence can be considered in varying orders and this process is called as domain shuffling.</div> <div>  ADVENTAGES</div> <div>  We can compare the different length sequences</div> <div>  Conserved domains can be determined from proteins</div> <div>  Common function features can be identified.</div> <div>  CONCLUSION</div> <div> Local alignment is used to compare the segments for high matching sequencing.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 43 of 320</div> <div> Module 019: ALIGNING, INSERTION & DELETION</div> <div> Insertion means addition of amino acids in protein sequence and addition of nucleotides in DNA sequences.</div> <div> And deletion means removal of amino acids from protein sequence and removal of nucleotides from DNA or RNA sequences.</div> <div>  ALIGNING INSERTION</div> <div> For example we have following two sequences</div> <div> To add the nucleotide in sequence 2 we will add gap first. And same happens with the deletion alignment we add gap where we delete the nucleotide from sequence. And such insertion of gap is called as –ve or plenty.</div> <div> 1: A C T G A C T G 2: A C G A C T G 1: A C T G A C T G 2: A C G A C T G 1: A C T G A C T G 2: A C . G A C T G</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 44 of 320</div> <div> Module 020: ALIGNING MUTATION IN SEQUENCES</div> <div> Removal and addition of amino acids in proteins and nucleotides in DNA, RNA by using Gaps named as Indels.</div> <div> Mutation is totally different from Indels, because in Mutation we replace the amino acid with other amino acids and replace the nucleotides with other and we don’t use Gaps is inserted in template or target for mutation.</div> <div>  CONCLUSION</div> <div> In identity alignment we use Gaps and in mutation we use substitution penalties and penalties depend upon the substitution.</div> <div> 1: A C T G A C T 2: A C G G A C T 1: A C T G A C T 2: A C G G A C T</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 45 of 320</div> <div> Module 021: INTRODUCTION TO DYNAMIC PROGRAMMING</div> <div> To find matching in nucleotides and amino acids of two sequences we use dot plot method. But dot plot cannot capture the insertions, deletions and gaps in the sequences.</div> <div> To deal with this situation we modify the dot plot.</div> <div> We represent the matching nucleotides with +1 while gaps, substitutions, insertions and mutations can be represented as -1 in dot plot. Dynamic programming is an algorithmic technique used commonly in sequence analysis. Dynamic programming is used when recursion could be used but would be inefficient because it would repeatedly solve the same sub problems. http://www.ibm.com/developerworks/library/j-seqalign/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 46 of 320</div> <div> Module 022: DYNAMIC PROGRAMMING ESSENTIALS</div> <div> When we talk about the compression of two sequences one by one it need time and is computationally expensive method. That’s why we need algorithm.</div> <div> In algorithm we calculate the step involve in sequence compression for example if we if we compare two sequences of length “n” than it would be “n2 “</div> <div> And its order is O (n2)</div> <div> Figure 0.11 -1 represent deletion, insertion and gaps while +1 represent matching nucleotides or amino acids</div> <div> One by one sequence compression is costly and time consuming process we minimize the cost with the help of algorithm.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 47 of 320</div> <div> Module 023: DYNAMIC PROGRAMMING METHODOLOGY</div> <div> Dynamic programming helps us to reduce the computational cost in sequence comparisons and it works on the method of “scoring function”.</div> <div> For example</div> <div> Match = +a</div> <div> Mismatch = -b</div> <div> Gap = -c</div> <div> Score = #match + #Mismatch +#Gaps</div> <div> Total =11</div> <div> All the alignments are done in diagonal way in dot plot matrix. For total score we make calculations in diagonal way and after calculation best one is selected.</div> <div> Match Rewards = 10 Mismatch penalty = 2 Gap penalty = 5 C T G T C G – C T G C - T G C – C G – T G - -5 10 10 -2 -5 -2 -5 -5 10 10 -5</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 48 of 320</div> <div> Module 024: NEEDLEMAN WUNSCH ALGORITHM-I</div> <div> In two different sequences alignments are arranged in a diagonal pattern of dot matrix. Total scores are captured for each alignment and at the end the best one is selected.</div> <div> Needleman Wunsch Algorithm is the way for alignment. The method is same like dot plot but it computes the scores in different way. We Start with a zero in the second row, second column. Move through the cells row by row, calculating the score for each cell. The score is calculated as the best possible score (i.e. highest) from existing scores to the left, top or top-left (diagonal). When a score is calculated from the top, or from the left this represents an indel in our alignment. When we calculate scores from the diagonal this represents the alignment of the two letters the resulting cell matches to. Given there is no 'top' or 'top-left' cells for the second row we can only add from the existing cell to the left. Hence we add -1 for each shift to the right as this represents an indel from the previous score. This results in the first row being 0, -1, -2, -3, -4, -5, -6, and -7. The same applies to the second column as we only have existing scores above. Thus we have:</div> <div> https://en.wikipedia.org/wiki/Needleman%E2%80%93Wunsch_algorithm</div> <div> Figure 11 Needleman-Wunsch pairwise sequence alignment</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 49 of 320</div> <div> Figure 12 Needleman wunsch algorithm way of computation of nucleotides</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 50 of 320</div> <div> Module 025: NEEDLEMAN WUNSCH ALGORITHM-II</div> <div> Alignments are represented by unbroken diagonal dot matrix plot. In this way we can create numerous combinations.</div> <div> Figure 13 various combinations of sequences through dot plot</div> <div> If the sequence is too long then there will be many diagonal alignments and at the end we select the best alignment by combinations of all. And for this we use Needleman Algorithm</div> <div> In Needleman Algorithm we use 0, 0 in first row and first column.</div> <div> Figure 14 initial column and row are kept zero (0)</div> <div> Left to right and top to bottom the best element (having high score) is selected.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 51 of 320</div> <div> Figure 15 maximum score element is selected from all three sides comparison</div> <div> The terms for match, mismatch are:</div> <div> The matrix is computed progressively until the bottom right element</div> <div> Alpha = Match reward Beta = Mismatch penalty Gamma = Gap penalty</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 52 of 320</div> <div> Module 026: NEEDLEMAN WUNSCH ALGORITHM-III</div> <div> We follow the diagonal route for scoring in Needleman Wuncsh Algorithm. Left to right and top to bottom in diagonal way we select the highest score.</div> <div> Figure 16 filling order of sequence alignment</div> <div> We select the best score to make best alignment and matrix is computed progressively until we reach to the bottom right.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 53 of 320</div> <div> Module 027: NEEDLEMAN WUNSCH ALGORITHM EXAMPLE</div> <div> Top left and diagonal element are considered to calculate an element in the matrix. Match, mismatch and gap penalty is computed from all there sides (Left to right) (Top to bottom) and (Diagonal).</div> <div> For example:</div> <div> Figure 17 score computation from all sides</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 54 of 320</div> <div> Figure 18 the best score is computed in diagonal way</div> <div> DNA, RNA and Protein sequences can be computed by using Needleman algorithm.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 55 of 320</div> <div> Module 028: BACKTRACKING ALIGNMENTS</div> <div> To find an optimal alignment in Needleman Wunsch Algorithm we use traceback method.</div> <div> Figure 19 traceback method starts from end maximum score.</div> <div> After completely matrix calculations we apply traceback to find the optimal alignment and traceback starts from bottom right (maximum score) to top side.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 56 of 320</div> <div> Module 029: REVISITING LOCAL AND GLOBAL ALIGNMENTS</div> <div> We use Needleman Algorithm to align the sequences in scoring way and traceback method to find the optimal alignment.</div> <div> In Needleman Algorithm we start traceback from bottom right towards top left progressively and this provides us the global alignment.</div> <div> Figure 20 Traceback in Needleman Algorithm</div> <div> We can start traceback from any point in matrix and smith waterman algorithm helps elicit the local alignments.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 57 of 320</div> <div> Module 030: OVERLAP MATCHES</div> <div> Dot plot and Needleman wucsch are algorithm method with little difference. Dot plot help us in finding matching residues of two sequences while Needleman wunsch helps us to find the global alignments.</div> <div> If some sequences have different regions of nucleotides which does not match to any other for that alignment we prefer Global alignment not local, but that does not penalize leading or trailing end.</div> <div> Figure 21 leading and trailing edge mismatches versus global alignment by gap-insertion (stretching) of sequences</div> <div> And “Traceback” is the technique by which we can check the sequences from any end of the matrix box. And such “Tracebacks” helps us to find the overlaps in aligned sequences.</div> <div> Figure 22 Traceback method</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 58 of 320</div> <div> Module 031: EXAMPLE</div> <div> A slight variation in traceback can helps us to find the overlaps in sequences and can apply some interesting strategies in sequences alignments.</div> <div> In following example of amino acids alignment we can understand the ways of tracback.</div> <div> Figure 23 Traceback in amino acid sequence alignment</div> <div> Scoring stagey is:</div> <div> Sequences are:</div> <div> Match = +2. Mismatch = -1, Gap = -2</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 59 of 320</div> <div> Module 032: MOVING FROM GLOBAL TO LOCAL ALIGNMENT</div> <div> DNA has coding and noncoding regions. Coding regions are called “EXON” expressed as protein and they remain more conserved due to their role in making functional proteins.</div> <div> And noncoding regions of DNA are called as “INTRONS” which are more likely involved in mutations than coding ones. It means high degree of alignment can be find among two exons.</div> <div> In local alignment we use small segments of sequences and through which we can find exons. Through this we can find “functional subunits”. However, the term exon is often misused to refer only to coding sequences for the final protein. This is incorrect, since many noncoding exons are known in human genes.</div> <div> (Zhang 1998)</div> <div> Zhang, M. Q. (1998). "Statistical features of human exons and their flanking regions." Human molecular genetics 7(5): 919-932.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 60 of 320</div> <div> Module 033: SMITH WATERMAN ALGORITHM</div> <div> In global alignment we compare the sequence from end to end but in local alignment we compare the sequences in segments.</div> <div> For Global sequences we use Needleman and Wunsch algorithm while for local pairwise alignment we use Smith and waterman.</div> <div> The Smith Freshman algorithm is different from Needle man.</div> <div>  Top row and Colum are set to zero.</div> <div>  Alignment can end anywhere.</div> <div>  Traceback starts from highest score.</div> <div> Local alignments can identify the coding portions of DNA and in this way we can find the functional domains from protein sequences.</div> <div> Figure 24 Global and local sequence alignment comparison.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 61 of 320</div> <div> Module 034: EXAMPLE OF SMITH WATERMAN ALGORITHM</div> <div> The only difference between Needleman and Smith Waterman is that zero “0” is placed in the</div> <div> relationship.</div> <div>  </div> <div>  </div> <div>  </div> <div>  </div> <div> </div> <div> </div> <div> </div> <div>  </div> <div> </div> <div> </div> <div>  </div> <div>  </div> <div>   </div> <div> </div> <div> 0</div> <div> , 1</div> <div> 1,</div> <div> [ 1, 1] ,</div> <div> , max</div> <div> </div> <div> </div> <div> C i j</div> <div> C i j</div> <div> C i j score i j</div> <div> C i j</div> <div> And in the matrix we place top line of zero and first Colum of zero.</div> <div> Figure 25 top line and first Colum are filled with zero in Smith Freshman Algorithm</div> <div> Local alignments can be extracted by starting from a high score till reaching ‘0’</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 62 of 320</div> <div> Module 035: REPEATED ALIGNMENTS</div> <div> We can find the best local alignments by using Smith Waterman algorithm.</div> <div> By making some change in strategy of traceback we can find the repeated sequences.</div> <div> We use threshold “T” score for matching and it avoids low scoring local alignment. And traceback can help us to find multiple aligned regions in multiple ways.</div> <div> Figure 27 (-5) is threshold score in table</div> <div> This threshold scoring method with some modifications in waterman algorithm can help us to find many matching sequence of amino acids or DNA.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 63 of 320</div> <div> Module 036: EXAMPLES OF REPEATED ALIGNMENTS</div> <div> Slight modification in waterman model can help us to find the Exons as well as the functional units in any sequence. Matches should be end at the threshold score or we should keep track of maximum score in sequence.</div> <div> Figure 28 Trackback from different sides to find maximum or Threshold (T) Score.</div> <div> Traceback should start from last element of the row and should reach at the top of row element and then move to the highest score of the Column. And this traceback is done twice and end at the point where score become “0, 0”</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 64 of 320</div> <div> Module 037: INTRODUCTION TO SCORING ALIGNMENTS</div> <div> There are two types of alignments;</div> <div>  Optimal Alignments</div> <div>  Best Alignment</div> <div> Scoring scheme used in sequences matches play crucial role in producing optimal alignment. An optimal alignment should be:</div> <div>  Appropriately rewarded for matches and mismatches.</div> <div>  INTRODUCTION</div> <div> We identify the pairs of symbol which most frequently appear in a sequences it helps us to find the substitution of specific pair of amino acid or nucleotide with other on in a sequence.</div> <div> For example AA nucleotide have a specific pattern of substitution. And same pairs of amino acids does in protein sequences because it can help to preserve the function of protein.</div> <div>  CONCLUSION</div> <div> Statically we can better align any sequence of protein or DNA, optimal gaps, penalties, insertions and deletions can be computed statically better.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 65 of 320</div> <div> Module 038: MEASURING ALIGNMENTS SCORES</div> <div> Score of match and mismatch both are equally observed while sequence alignments.</div> <div> For example:</div> <div> The matrix has positive and negative scores both, matches and mismatches therefore are all considered because it’s a diagonal pattern.</div> <div> If we build such scoring matrixes with matches and mismatches we can we can sequence in according to real life.</div> <div> Figure 29 Needleman wunsch algorithm match, mismatch scoring</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 66 of 320</div> <div> Module 039: SCORING MATRICES</div> <div> Alignments are used to align the biological sequences. Amino acids and nucleotides are more easily substituted because they have similar chemical nature.</div> <div> As amino acids are substituted with many probabilities that’s why we need flexible scoring. And we use Scoring Matrices contain such flexible scoring during alignment.</div> <div> To build the Scoring Matrices we analyze the amino acids and nucleotides which are substituted in single gene and protein sequence.</div> <div> Scoring Matrices have both values +ve and –ve. Positive value for matches and negative value for mismatches.</div> <div> Figure 0.12 Ubiquitin Protein where amino acids matching</div> <div> Different type of scoring matrices can be developed based on underlying strategy.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 67 of 320</div> <div> Module 040: DERIVING SCORING MATRICES</div> <div> Each amino acid have different property.</div> <div> Figure 0.13 properties of amino acids (Image Esquivel et al. (2013)</div> <div> And each amino acid have different frequency.</div> <div> When we compare the sequences they match and mismatch according to their frequency.</div> <div> For example.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 68 of 320</div> <div> Based on frequencies we match and mismatch the sequence alignments for scoring.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 69 of 320</div> <div> Module 041: PAM MATRICES</div> <div> Alignment matrices scoring is very useful method to score the sequences alignment for match and mismatch.</div> <div> There are two types of scoring matrices.</div> <div>  PAM</div> <div>  BLOSUM</div> <div> PAM means “Point Accepted Mutation”</div> <div> Point accepted mutations means the substitution of one amino acid in a sequence with another that protein function remain conserved.</div> <div>  PAM UNIT</div> <div> PAM unit is actually that time during which 1% amino acid undergo for acceptable mutation. If two sequences diverge by 100 PAM units, it does not mean that they will be at totally different positions.</div> <div>  STEP TOCOMPUTE PAM MATRICES</div> <div> 1. Align the protein sequence which are 1-PAM Unit diverge.</div> <div> 2. Let Ai,i be the number of times Ai is substituted by Ai.</div> <div> 3. Compute the frequency fi of amino acid Ai.</div> <div> Then, PAM1=pii=</div> <div> PAM ‘n’= (PAM1)n</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 70 of 320</div> <div> Module 042: BLOSUM MATRICES</div> <div> BLOSUM matrices can be used to align the protein sequences. BLOSUM matrices was first purposed in 1992 by Henikoff et al.</div> <div> BLOSUM matrices is also called the Block substitution matrix without any gap although it has mismatches in sequences.</div> <div> There are three steps to compute the BLOUSM Matrices.</div> <div> Step 1: Eliminate sequences that are identical in x% positions</div> <div> Step 2: Compute observed frequency f i, j of aligned pair Ai to Aj. Hence, f i,j becomes the probability of aligning Ai and Aj in the selected blocks.</div> <div> Step 3: Compute fi which is the frequency of observing Ai in the entire block</div> <div> Typically used matrices: BLOSUM62 or PAM120 in PAMx, larger x detects more divergent sequences.</div> <div> Figure 0.14 sequence of amino acids which have mismatch but no gap. Figure 0.15 formula for computation of BLOSUM MATRICES.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 71 of 320</div> <div> Module 043: MULTIPLE SEQUENCE ALIGNMENT</div> <div> In pair wise sequence alignments we use pairs of sequence to compare them. And scoring matrices were used to score the sequence ranks.</div> <div> In Multiple sequence Alignments we compare multiple number of protein and DNA sequences to identify the matches and mismatches.</div> <div> For pair wise alignment we use Dynamic programming but for multiple alignment it would be very expensive computationally. So solution for this is progressive alignment.</div> <div> M Q V K L F T P L H D K S D H G K Y H M Q V K I F T P L H D K S - H G K S H M Q V H L Y - P L H D K S - T G K S H M Q V H L F - P L H D K S D T G K S H M Q V K L Y T P L H D K S D H G K Y HFigure 0.16 multiple sequence alignment</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 72 of 320</div> <div> Module 044: MORE ON MULTIPLE SEQUENCE ALIGNMENT</div> <div> MSA helps compare several sequences by aligning them. MSA can extract consensus sequences from several aligned sequences. Characterize protein families based on homologous regions.</div> <div> APPLICATION OF MSA</div> <div>  Predict secondary and tertiary structures of new protein sequences</div> <div>  Evaluate evolutionary order of species or “Phylogeny”</div> <div> METHODOLOGY</div> <div>  Pairwise alignment is the alignment of two sequences</div> <div>  MSA can be performed by repeated application of pairwise alignment</div> <div> Figure 0.17 Methodology</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 73 of 320</div> <div> Figure 0.18 sequence alignments</div> <div> CONCLUSION</div> <div> MSA can help align multiple sequences. Progressive alignment can help perform MSA. Need to remove sequences with >80% similarity.</div> <div> Figure 0.19 CLUSTAL – Online tool</div> <div> http://www.ebi.ac.uk/Tools/msa/clustalo/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 74 of 320</div> <div> Module 045: PROGRESSIVE ALIGNMENT FOR MSA</div> <div> MSA involves progressive alignment of sequences. Doing so many progressive alignments can be slow.</div> <div> STEPS:</div> <div>  Step 1 : Pairwise Alignment of all sequences</div> <div> Example: S1, S2, S3, S4, so that is 6 pairwise comparisons.</div> <div>  Step 2: Construct a Guide Tree (Dendogram) using a Distance Matrix.</div> <div>  Step 3: Progressive alignment following branching order in tree.</div> <div> Figure 0.20 Similarity Matrix</div> <div>  SHORTCOMING OF THIS APPROACH</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 75 of 320</div> <div>  Dependence upon initial alignments</div> <div>  If sequences are dissimilar, errors in alignment are propagated</div> <div>  Solution: Begin by using an initial alignment, and refine it repeatedly</div> <div> Progressive alignments are used in aligning multiple sequences. Iterative approaches can help refine results from progressive alignments.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 76 of 320</div> <div> Module 046: MSA-EXAMPLE</div> <div> MSA involves progressive alignment of sequences. Doing so many progressive alignments can be slow.</div> <div> For example:</div> <div> Figure 0.21 MSA on globin sequences</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 77 of 320</div> <div> Figure 0.22 Progressive alignment using sequential branching</div> <div> Figure 0.23 Progressive alignment following a guide tree</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 78 of 320</div> <div> Figure 0.24 Alignment results</div> <div> MSA can be better performed using clustering strategies followed by alignment of the alignments later. CLUSTAL is a free online tool that does all of this for us!</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 79 of 320</div> <div> Module 047: CLUSTALW</div> <div> MSA involves progressive alignment of sequences. Doing so many progressive alignments can be slow. CLUSTALW is an online tool to perform MSA.</div> <div> Developed by European Molecular Biology Laboratory & European Bioinformatics Institute. Performs alignment in:</div> <div>  slow/accurate</div> <div>  fast/approximate</div> <div> SCOPE</div> <div>  create multiple alignments,</div> <div>  optimize existing alignments,</div> <div>  profile analysis &</div> <div>  create phylogenetic trees</div> <div> http://www.genome.jp/tools/clustalw</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 80 of 320</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 81 of 320</div> <div> Module 048: INTRODUCTION TO BLAST-I</div> <div> National Center for the Biotechnology Information (NCBI) – USA. BLAST developed in 1990. “Basic Local Alignment Search Tool”. Searches databases for query protein and nucleotide sequences. Also searches for translational products etc. Online availability</div> <div> www.blast.ncbi.nlm.nih.gov/Blast.cgi</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 82 of 320</div> <div> BLAST can be used to search for local alignment of protein and nucleotide sequences. It is available online. Can perform searches across species and organisms</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 83 of 320</div> <div> Module 049: INTRODUCTION TO BLAST-II</div> <div> National Center for the Biotechnology Information (NCBI) – USA. BLAST developed in 1990. “Basic Local Alignment Search Tool”. Searches databases for query protein and nucleotide sequences. Also searches for translational products etc. Online availability</div> <div> www.blast.ncbi.nlm.nih.gov/Blast.cgi</div> <div> Smith Waterman can align complete sequences. BLAST does it in an approximate way. Hence, BLAST is faster BUT does not ensure optimal alignment. BLAST provides for approximate sequence matching. Input to BLAST is a FASTA formatted sequence and a set of search parameters</div> <div> OUTPUT OF BLAST</div> <div> Results are shown in HTML, plain text, and XML formats. A table lists the sequence hits found along with scores. Users can read this table off and evaluate results</div> <div> Figure 0.25Input to BLAST: Gene IDs</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 84 of 320</div> <div> Figure 0.26 Input to BLAST: Protein IDs</div> <div> Figure 0.27 Results from BLAST</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 85 of 320</div> <div> Module 050: BLAST ALGORITHM</div> <div> BLAST can search sequence databases and identify unknown sequences by comparing them to the known sequences. This can help identify the parent organism, function and evolutionary history.</div> <div> For example:</div> <div> Query sequence: PQGELV</div> <div> Make list of all possible worlds (length 3 for proteins)</div> <div> PQG (score 15)</div> <div> QGE (score 9)</div> <div> GEL (score 12)</div> <div> ELV (score 10)</div> <div> Assign scores from Blosum62, use those with score> 11: PQG & GEL</div> <div> Mutate words such that score still > 11</div> <div> PQG (score 15) similar to PEG (score 13)</div> <div> At the end, we get: PQG, GEL and PEG</div> <div> Find all database sequences that have at least 2 matches among our 3 words: PQG, GEL & PEG. Find database hits and extend alignment (High-scoring Segment Pair):</div> <div> High Scoring Pair: PQGI (score 8+5+5+2)</div> <div> If 2 HSP in query sequence are < 40 positions away</div> <div> Full dynamic alignment on query and hit sequences</div> <div> BLAST performs quick alignments on sequences. The results are tabulated with alignment regions overlapping each other. Statistical evaluation is also provided alongside</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 86 of 320</div> <div> Module 051: TYPES OF BLAST</div> <div> BLAST can search sequence databases and identify unknown sequences by comparing them to the known sequences. This can help identify the parent organism, function and evolutionary history.</div> <div> There are two main types of BLAST.</div> <div> Nucleotides</div> <div> • Blastn: Compares a nucleotide query sequence against a nucleotide database.</div> <div> Proteins</div> <div> • Blastp: Compares an amino acid query sequence against a protein database.</div> <div> There are also many other types of BLAST:</div> <div>  Blastx:</div> <div> Compares a nucleotide query sequence against a protein sequence database.</div> <div> Helps find potential translation products of unknown nucleotide sequences</div> <div>  tblastn:</div> <div> Compares a protein query sequence against a nucleotide sequence database</div> <div> Nucleotide sequence dynamically translated into all reading frames</div> <div>  tblastx:</div> <div> Compares the six-frame translated proteins of a nucleotide query sequence against the six-frame translated proteins of a nucleotide sequence database.</div> <div> • BLAST performs quick alignments on biological sequences</div> <div> • Several types of BLAST exist which can assist in comparing nucleotide sequences with amino acids and vice versa</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 87 of 320</div> <div> Module 052: SUMMERY OF BLAST</div> <div> BLAST can search sequence databases and identify unknown sequences by comparing them to the known sequences. This can help identify the parent organism, function and evolutionary history.</div> <div> Step1: obtain a query of sequence</div> <div> Step2: choose a type of BLAST</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 88 of 320</div> <div> Step3: search parameter</div> <div> Step4: tabulated search results</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 89 of 320</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 90 of 320</div> <div> Figure 0.28 tabulated search results</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 91 of 320</div> <div> Module 053: INTRODUCTION TO FASTA</div> <div> For comparing two sequences we use pair wise sequencing and for the comparison of many sequences we use multiple sequence alignment. To handle the multiple alignments we perform alignment through smith-waterman algorithm for local one. And for global alignment we use Needleman-wunsch algorithm.</div> <div> Both local and global alignments are the dynamic approaches. Many of the sequences are compared, which takes time and we use BLAST which is an approximate local alignment search tool BLAST compares a large number of sequences, quickly. FASTA took a similar approach.</div> <div> Developed in 1988.it does Fast Alignment .Searches databases for query protein and nucleotide sequences. Was later improved upon in BLAST.</div> <div> Figure 0.29 Regions of absolute identity</div> <div> http://www.ebi.ac.uk/Tools/sss/fasta/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 92 of 320</div> <div> Figure 0.30 Nucleotide FASTA Figure 0.31 Protein FASTA</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 93 of 320</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 94 of 320</div> <div> Module 054: INTRODUCTION TO FASTA-II</div> <div> FASTA – Fast Alignment Algorithm. Classical global and local alignment algorithms are time consuming. FASTA achieves alignment by using short lengths of exact matches.</div> <div>  USES OF FASTA</div> <div> FASTA relies on aligning subsequences of absolute identity. Input to FASTA search can be in FASTA, EMBL, GenBank, PIR, NBRF, PHYLIP or UniProt formats</div> <div>  OUTPUT OF BLAST</div> <div> Results are output in visual format along with functional prediction. Makes table lists the sequence hits found along with scores. Users can click on each reported match to look at the details.</div> <div> Figure 0.32 Input to FASTA: Gene IDs</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 95 of 320</div> <div> Figure 0.33 Input to FASTA: Protein Sequence</div> <div> Figure 0.34 Results from FASTA</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 96 of 320</div> <div> Module 055: FASTA ALGORITHM</div> <div> FASTA can search sequence databases and identify unknown sequences by comparing them to the known sequence databases. This can help obtain information on the parent organism, function and evolutionary history.</div> <div> STEP1: Local regions of identity are found</div> <div> STEP2: Rescore the local regions using PAM or BLOSUM matrix</div> <div> STEP3: Eliminate short diagonals below a cutoff score</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 97 of 320</div> <div> STEP4: Create a gapped alignment in a narrow segment and then perform Smith Watermann alignment</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 98 of 320</div> <div> Module 056: TYPES OF FASTA</div> <div> There are six types of FASTS:</div> <div>  fasts35</div> <div> Compare unordered peptides to a protein sequence database</div> <div>  fastm35</div> <div> Compare ordered peptides (or short DNA sequences) to a protein (DNA) sequence database</div> <div>  Fasta35</div> <div> Scan a protein or DNA sequence library for similar sequences</div> <div>  Fastx35</div> <div> Compare a translated DNA sequence (6 ORFs) to a protein sequence database</div> <div>  tfastx35</div> <div> Compare a protein sequence to a DNA sequence database (6 ORFs)</div> <div>  fasty35</div> <div> Compare a DNA sequence (6ORFs) to a protein sequence database</div> <div> FASTA performs quick alignments on biological sequences. Several types of FASTA exist which can assist in comparing DNA/RNA/Protein sequences with each other</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 99 of 320</div> <div> Module 057: SUMMERY OF FASTA</div> <div> FASTA can briskly perform sequence search databases if given a query sequence. Multiple types of FASTA exist which assist in aligning DNA/RNA/Protein sequences</div> <div> Figure 0.35 Step 1: Obtain a query sequence</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 100 of 320</div> <div> Figure 0.36 Step 2: Choose a type of FASTA</div> <div> Figure 0.37 Type of FASTA</div> <div> http://fasta.bioch.virginia.edu/fasta_docs/fasta35.shtml</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 101 of 320</div> <div> Figure 0.38 Step 3: Setup Search Parameters</div> <div> Figure 0.39 Step 4: Tabulated Search Results</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 102 of 320</div> <div> Figure 0.40 Tabulated data</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 103 of 320</div> <div> Module 058: BIOLOGICAL DATABASE AND ONLINE TOOLS</div> <div> All molecular information of RNA, DNA, Proteins have need to be stored and retrieved. Sequences are obtained from genome sequencing and mass spectrometry</div> <div> Structures are obtained from X-Ray Crystallography, Atomic Force Microscopy & Nuclear Magnetic Resonance Spectroscopy.</div> <div> Vast amounts of such data exists. Moreover, this data is rapidly accumulating. Online Databases are formed to store and share this data.</div> <div>  OBJECTIVE</div> <div>  Make biological data available to scientists in computer-readable form</div> <div>  For handling, sharing and analysis of the data</div> <div>  The best way to share is to keep this data on the web</div> <div> Several sequence, structure and molecular interaction databases exist. These are available online on the web. Users can freely access and download such data</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 104 of 320</div> <div> Module 059: EXPASY</div> <div> It is developed by Swiss Bioinformatics Institute (SIB). Website provides access to databases and tools Proteomics, genomics, phylogeny, systems biology, population genetics, transcriptomics etc. can be searched.</div> <div> http://www.expasy.org/</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 105 of 320</div> <div> Figure 0.41 flowchart</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 106 of 320</div> <div> Figure 0.42 prosite scanning section</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 107 of 320</div> <div> Figure 0.43 peptide mass finding</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 108 of 320</div> <div> Figure 0.44 for local use of protein sequencing</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 109 of 320</div> <div> Figure 0.45 potential protein finding tool</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 110 of 320</div> <div> Module 060: UNIPROT AND SWISSPROT</div> <div> Both UniProt and SwissProt are the online database for proteins.</div> <div> Figure 0.46 gene, protein or chemical can be find</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 111 of 320</div> <div> Figure 0.47 online database for proteins</div> <div> Swiss-Prot contains human curated protein information</div> <div>  Accession number, unique identifier</div> <div>  The sequence</div> <div>  Molecular mass</div> <div>  Observed and predicted modifications</div> <div> Protein sequences from various species and organisms can be found in uniprot. SwissProt is the manually annotated version of the UniProt Database.</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 112 of 320</div> <div> Module 061: PROTEIN DATA BANK</div> <div> Protein Data Bank is the premier resource of protein structures. These structures have been determined using experimental techniques. It’s Open & Free</div> <div> Figure 0.48 protein data bank</div> <div> Figure 0.49 P0CG47 - UBB_HUMAN</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 113 of 320</div> <div> Figure 0.50 P0CG47 - UBB_HUMAN</div> <div> Figure 0.51 searched results</div> <div> Protein Data Bank provides Cartesian coordinates of each atom in the protein structure. Over 50,000 protein structures are reported and present in this database</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 114 of 320</div> <div> Module 062: REVIEW OF SEQUENCE ALIGNMENT</div> <div> We use next generation sequencing and whole genome sequencing to obtain the genetic information. For protein sequencing we use Mass Spectrometry and Edman Degradation</div> <div> STORAGE:</div> <div>  Sequence information is stored digitally</div> <div>  Databases are designed to store sequence data</div> <div>  Several databases exist depending on the type of sequence data</div> <div> SHARING AND ACCESS:</div> <div>  Sequence databases are shared via online websites</div> <div>  Access to several such websites is free</div> <div>  Data can be downloaded or searched on these website</div> <div> USAGE OF DATA:</div> <div> Sequence data can be used to obtain:</div> <div>  Similarity of sequences</div> <div>  Evolutionary History</div> <div>  Predict the function of molecules</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 115 of 320</div> <div> Module 063: GENBANK</div> <div>  Developed by Swiss Bioinformatics Institute (SIB)</div> <div>  Website provides access to databases and tools</div> <div>  Proteomics, genomics, phylogeny, systems biology, population genetics, transcriptomics etc.</div> <div> Figure 0.52 http://www.ncbi.nlm.nih.gov/genbank/</div> <div> Several sequence, structure and molecular interaction databases exist. These are available online on the web. Users can freely access and download such data</div> <div> Handout Chapter 02: Sequence Analysis</div> <div> Page 116 of 320</div> <div> Module 064: ENSEMBLE</div> <div> As human brain is limited to remember and store the information for long time that’s why we use online database for the storage of Molecular information.</div> <div> ESEMBLE is genome search engine which is used to search the genome of every recorded species.</div> <div> http://asia.ensembl.org/index.html</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 117 of 320</div> <div> Chapter 3 - Molecular Evolution</div> <div> Module 001: Molecular Evolution & Phylogeny Molecular evolution is the process of change in the sequence composition of cellular molecules such as DNA, RNA, and proteins across generations. The field of molecular evolution uses principles of evolutionary biology and population genetics to explain patterns in these changes. Genes and Proteins are modified in this process.</div> <div> All molecules have an evolutionary history. Phylogenetics is the science of studying evolutionary relationships. Phylogenetics has led to the creation of relationship trees between various species of Bacteria, Archaea, and Eukaryota.</div> <div> (Page and Holmes 2009)</div> <div> Figure 0.1 Phylogenetic tree</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 118 of 320</div> <div> Types of Phylogenetic Trees</div> <div> Scaled Trees</div> <div> • Branch lengths are equal to the magnitude of change in the nodes</div> <div> Unscaled Trees</div> <div> • Only representing the relationship between sequences</div> <div> Figure 0.2 phylogenetic tree interference</div> <div> Conclusion</div> <div> Phylogenetics is the study of extracting evolutionary relationships between species. Sequence information from each species is used to measure the difference between the species.</div> <div> Page, R. D. and E. C. Holmes (2009). Molecular evolution: a phylogenetic approach, John Wiley & Sons.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 119 of 320</div> <div> Module 002: Evolution of Sequences</div> <div> DNA acts as cellular memory unit and protein are the translated product of DNA coded information. And evaluation is very important to survive in different type of environments. There are some methods which brings change or evolution in any organism. (Kluger 2015)</div> <div> Method of Change</div> <div> DNA gets modified by:</div> <div>  Mutation & Substitution</div> <div>  Insertion</div> <div>  Deletion</div> <div> Discussion</div> <div> Over time, species evolve to adapt to their circumstances. Since the environment and circumstances may be different for each species, they evolve uniquely. Unique evolutionary pressures may be encountered by each cell for struggle of life. However, in which sequence they are presented to the cells is also unique. Combinations of evolutionary factors are involve in evolution. The evolutionary events and their combination impart relationships between sequences. These relationships are explored in Phylogenetics .Several algorithms exist for finding such relationships</div> <div> Kluger, M. J. (2015). Fever: its biology, evolution, and function, Princeton University Press.</div> <div> Page, R. D. and E. C. Holmes (2009). Molecular evolution: a phylogenetic approach, John Wiley & Sons.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 120 of 320</div> <div> Module 003: Concepts and Terminologies - I</div> <div> To understand the concept of evolution we follow some rules. Phylogenetics involves processing sequence information from different species to find evolutionary relationships. Output from such studies include Phylogenetic Trees</div> <div> Figure 3 phylogenetic tree from ancestor to evolution</div> <div> In above figure the point A stands for ancestor and with the passage of time the evolution occurred with and the genome sequence of organisms changed.</div> <div> Figure 4 layout of trees</div> <div> All trees have same meanings.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 121 of 320</div> <div> Figure 5 rooted tree</div> <div> Root node is the ancestor of all other nodes. The direction of evolution is from ancestor to the terminal nodes.</div> <div> Conclusion</div> <div> Phylogenetics specifies evolutionary relationship with the help of trees. Trees can be rooted or unrooted. Rooted trees can show temporal evolutionary direction.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 122 of 320</div> <div> Module 004: Concepts and Terminologies - II</div> <div> Rooted and Unrooted trees can be used to show phylogenetic relationships between sequences. Let’s examine the properties of these trees further.</div> <div> Figure 6 rooted tree vs unrooted</div> <div> Rooted trees are computationally expensive. http://everything.explained.today/Computational_phylogenetics/</div> <div> https://github.com/joey711/phyloseq/issues/597</div> <div> Figure 7 computation comparison</div> <div> Conclusion</div> <div> Rooted and Unrooted trees have their own advantages and disadvantages. Depending on our requirement, we can choose between them.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 123 of 320</div> <div> Module 005: Algorithms and Techniques</div> <div> Rooted and Unrooted trees can be used to show phylogenetic relationships between sequences. Several types of algorithms exist which are divided into two classes. There are many methods for constructing evolutionary trees.</div> <div> Figure 9 construction methods UPGMA (Unweighted Pair Group Method with Arithmetic Mean) is a simple agglomerative (bottom-up) hierarchical clustering method. The method is generally attributed to Sokal and Michener.</div> <div> In this method two sequences with with the shortest evolutionary distance between them are considered and these sequences will be the last to diverge, and represented by the most recent internal node.</div> <div> Least Squares Distance Method. Branch lengths, represent the “observed” distances between sequences (i & j).</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 124 of 320</div> <div> Find X, Y and Z such that D (i, j) are conserved?</div> <div> Conclusion</div> <div>  Several methods exist for constructing phylogenetic trees.</div> <div>  Broadly, they belong to objective methods or clustering methods.</div> <div>  We will study UPGMA and Distance Methods.</div> <div> Module 006: Introduction to UPGMA</div> <div> Phylogenetic trees can be used to show phylogenetic relationships between sequences. To construct these trees, several types of algorithms exist which are divided into two classes.</div> <div> UPGMA: Unweighted Pair – Group Method using arithmetic Averages</div> <div>  Calculating distance between two clusters:</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 125 of 320</div> <div> Cluster X + Cluster Y = Cluster Z</div> <div> Calculate the distance of a cluster (e.g. W) to the new cluster Z</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> Nx is the number of sequences in cluster x</div> <div>  Calculating distance between two trees:</div> <div> Assume we have N sequences</div> <div> Cluster X has NX sequences, cluster Y has NY sequences</div> <div> dXY : the evlotionary distance between X and Y</div> <div> </div> <div>  </div> <div> </div> <div> i X j Y</div> <div> ij</div> <div> X Y</div> <div> XY d</div> <div> N N</div> <div> d</div> <div> ,</div> <div> 1</div> <div>  Methods for constructing trees</div> <div> The distance matrix is obtained using pairwise sequence alignment.</div> <div>  Calculating distance between two clusters:</div> <div> Cluster X + Cluster Y = Cluster Z</div> <div> Calculate the distance of a cluster (e.g. W) to the new cluster Z</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> Nx is the number of sequences in cluster x</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 126 of 320</div> <div>  Calculating distance between two trees:</div> <div> Assume we have N sequences</div> <div> Cluster X has NX sequences, cluster Y has NY sequences</div> <div> dXY : the evlotionary distance between X and Y</div> <div> </div> <div>  </div> <div> </div> <div> i X j Y</div> <div> ij</div> <div> X Y</div> <div> XY d</div> <div> N N</div> <div> d</div> <div> ,</div> <div> 1</div> <div>  Methods for constructing trees</div> <div> A – D becomes a new cluster lets say V. We have to modify the distance matrix. What are the</div> <div> distances between:</div> <div>  V and B (Calculate),</div> <div>  V and C,</div> <div>  V and E,</div> <div>  V and F.</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> 6</div> <div> 1 1</div> <div> 1*6 1*6</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> A D</div> <div> A AB D DB</div> <div> VB N N</div> <div> N d N d</div> <div> d</div> <div> Conclusion</div> <div> UPGMA is a clustering algorithm which can help us compute phylogenetic trees. We will see the</div> <div> detailed working of this approach in later modules.</div> <div> Module 007: UPGMA-I</div> <div> UPGMA has two components to it. These include distance calculations between two clusters</div> <div> and between two trees.</div> <div>  Building trees using UPGMA</div> <div> Combining Clusters: Cluster X + Cluster Y = Cluster Z</div> <div> Calculate the distance of each cluster (e.g. W) to the new cluster Z</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 127 of 320</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> Nx is the number of sequences in cluster x</div> <div>  Calculating the distance between two trees</div> <div> Assume we have N sequences</div> <div> Cluster X has NX sequences, cluster Y has NY sequences</div> <div> dXY : the evlotionary distance between X and Y</div> <div> </div> <div>  </div> <div> </div> <div> i X j Y</div> <div> ij</div> <div> X Y</div> <div> XY d</div> <div> N N</div> <div> d</div> <div> ,</div> <div> 1</div> <div> Figure 10 the distance matrix is obtained using pairwise sequence alignment</div> <div>  Methods for constructing trees</div> <div> A – D becomes a new cluster lets say V.</div> <div> We have to modify the distance matrix!</div> <div> What are the distances between:</div> <div>  V and B (Calculate),</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 128 of 320</div> <div>  V and C,</div> <div>  V and E,</div> <div>  V and F.</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> 6</div> <div> 1 1</div> <div> 1*6 1*6</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> A D</div> <div> A AB D DB</div> <div> VB N N</div> <div> N d N d</div> <div> d</div> <div> Conclusion</div> <div> UPGMA starts with creating clusters of sequences which are the closest. Next, distance is</div> <div> computed between the new cluster and the remaining sequences. The process is repeated for</div> <div> all sequences.</div> <div> Module 008: UPGMA-II</div> <div> UPGMA steps include distance calculations between two clusters and between two trees. We</div> <div> formed clusters from sequences which had the shortest distance.</div> <div> Building trees using UPGMA</div> <div> Combining Clusters: Cluster X + Cluster Y = Cluster Z</div> <div> Calculate the distance of each cluster (e.g. W) to the new cluster Z</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 129 of 320</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> Nx is the number of sequences in cluster x</div> <div> Calculating the distance between two trees</div> <div> Assume we have N sequences</div> <div> Cluster X has NX sequences, cluster Y has NY sequences</div> <div> dXY : the evlotionary distance between X and Y</div> <div> </div> <div>  </div> <div> </div> <div> i X j Y</div> <div> ij</div> <div> X Y</div> <div> XY d</div> <div> N N</div> <div> d</div> <div> ,</div> <div> 1</div> <div> Methods for constructing trees</div> <div> The distance matrix is obtained using pairwise sequence alignment.</div> <div> 6</div> <div> 1 1</div> <div> 1*6 1*6</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> A D</div> <div> A AB D DB</div> <div> VB N N</div> <div> N d N d</div> <div> d</div> <div> 8</div> <div> 1 1</div> <div> 1*8 1*8</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> A D</div> <div> A AC D DC</div> <div> VC N N</div> <div> N d N d</div> <div> d</div> <div> 2</div> <div> 1 1</div> <div> 1*2 1*2</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> A D</div> <div> A AE D DE</div> <div> VE N N</div> <div> N d N d</div> <div> d</div> <div> 6</div> <div> 1 1</div> <div> 1*6 1*6</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> A D</div> <div> A AF D DF</div> <div> VF N N</div> <div> N d N d</div> <div> d</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 130 of 320</div> <div> V – E becomes a new cluster lets say W</div> <div> Now we have to modify the distance matrix again.</div> <div> What are the distances between:</div> <div> W and B,</div> <div> W and C,</div> <div> W and F.</div> <div> Conclusion</div> <div> Once a cluster is selected and its distance is computed with all other sequences, we update the distance matrix. Next, we select the shortest distance from the new matrix and repeat the process.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 131 of 320</div> <div> Module 009: UPGMA-III</div> <div> UPGMA has two components to it. These include progressive distance calculations between</div> <div> two clusters or between two trees.</div> <div> Building trees using UPGMA</div> <div> Combining Clusters: Cluster X + Cluster Y = Cluster Z. Calculate the distance of each cluster (e.g.</div> <div> W) to the new cluster Z.</div> <div> X Y</div> <div> X XW Y YW</div> <div> ZW N N</div> <div> N d N d</div> <div> d</div> <div> </div> <div> </div> <div> </div> <div> Nx is the number of sequences in cluster x</div> <div> Calculating the distance between two trees</div> <div>  Assume we have N sequences</div> <div>  Cluster X has NX sequences, cluster Y has NY sequences</div> <div>  dXY : the evlotionary distance between X and Y</div> <div> </div> <div>  </div> <div> </div> <div> i X j Y</div> <div> ij</div> <div> X Y</div> <div> XY d</div> <div> N N</div> <div> d</div> <div> ,</div> <div> 1</div> <div> V – E becomes a new cluster lets say W. Now we have to modify the distance matrix again.</div> <div> What are the distances between:</div> <div> W and B,</div> <div> W and C,</div> <div> W and F.</div> <div> 6</div> <div> 2 1</div> <div> 2*6 1*6</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> V E</div> <div> V VB E EB</div> <div> WB N N</div> <div> N d N d</div> <div> d</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 132 of 320</div> <div> 8</div> <div> 2 1</div> <div> 2*8 1*8</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> V E</div> <div> V VC E EC</div> <div> WC N N</div> <div> N d N d</div> <div> d</div> <div> 6</div> <div> 2 1</div> <div> 2*6 1*6</div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> </div> <div> V E</div> <div> V VF E EF</div> <div> WF N N</div> <div> N d N d</div> <div> d</div> <div> New matrix</div> <div> Cluster according to min distance</div> <div> Conclusion</div> <div> Now we have formed three clusters. Also, two separate trees have been formed. Next, we need</div> <div> to join these trees to create a complete tree.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 133 of 320</div> <div> Module 010: UPGMA-IV</div> <div> Application of UPGMA resulted in formation of two sub-trees. The need now was to join them</div> <div> into a single tree. Let’s see how that is done.</div> <div> F – B becomes a new cluster lets say X. We have to modify the distance matrix yet again. What</div> <div> is the distance between trees: W and X.</div> <div> *(6 6 6 6 6 6) 6</div> <div> 3*2</div> <div> 1</div> <div> ( )</div> <div> 1</div> <div> 1</div> <div> ,</div> <div>      </div> <div>      </div> <div>   </div> <div>  </div> <div> AB AF DB DF EB EF</div> <div> W X</div> <div> i W j X</div> <div> ij</div> <div> W X</div> <div> WX</div> <div> d d d d d d</div> <div> N N</div> <div> d</div> <div> N N</div> <div> d</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 134 of 320</div> <div> X – W becomes a new cluster lets say Y. We have to modify the distance matrix</div> <div> What is the distance between: Y and C.</div> <div> Conclusion</div> <div> We have now seen how trees are generated and connected. Next, we need to finalize the tree by adding the last two clusters.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 135 of 320</div> <div> Module 011: UPGMA-V</div> <div> Application of UPGMA resulted in formation of two sub-trees. The need now was to join them into a single tree. Let’s see how that is done.</div> <div> X – W becomes a new cluster lets say Y. We have to modify the distance matrix. What is the distance between: Y and C.</div> <div> Handout Chapter 03: Molecular Evolution</div> <div> Page 136 of 320</div> <div> Conclusion</div> <div> Un-weighted Pair Group Method using Arithmetic Averages is a clustering method to construct phylogenetic trees. Non-clustering methods such as Maximum Parsimony may be used for making trees as well.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 137 of 320</div> <div> Chapter 4 - RNA Secondary Structure Prediction</div> <div> Module 001: DNA TO RNA SEQUENCES</div> <div>  MOTIVATION</div> <div> In early days RNA was a considered as a structure which was involve between DNA and protein, means takes information from DNA and converts that information into protein synthesis. Now we know that it has multiple types like mRNA, tRNA, miRNA and siRNA. And they perform most of the work in gene expression and proteins. Not All RNA molecules are same, they differ in nucleotide sequences and functions also.</div> <div> Many viruses assemble their genomes from RNAs. They are therefore called RNA viruses. Examples include Human Immunodeficiency Virus and Hepatitis C Virus.</div> <div> There is little difference between RNA and DNA:</div> <div>  Thymine is replaced by Uracil in RNA molecule.</div> <div>  RNA molecule is single stand.</div> <div>  RNA contain ribose sugar.</div> <div> Because RNA has two (OH) groups that’s why it has shot life spam because of both (OH) repulsion.</div> <div> Figure 0.1 Ribose sugar has (OH) and Deoxyribose (H)</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 138 of 320</div> <div> Module 002: TYPES OF RNA &THEIR FUNCTIONS</div> <div> There are two categories of RNA:</div> <div>  Coding RNA</div> <div>  Non-Coding RNA</div> <div> Coding RNA perform their coded function in protein synthesis. And Non-coding RNA helps in translation process.</div> <div>  TYPES OF RNA</div> <div> There are many types of RNA according to their funtions like:</div> <div>  Messenger RNA (mRNA)</div> <div>  Transfer RNA (tRNA)</div> <div>  Ribosomal RNA (rRNA)</div> <div>  Micro RNAs (miRNA)</div> <div>  Small Interfering RNA (siRNA)</div> <div> </div> <div>  MESSENGER RNA</div> <div> Only 5-10% of this RNA type is present in cell. Which has variable sequence, variable size and it carries the genetic information form DNA to Ribosomes where proteins to be assembled. Messenger RNA 5’ end is capped with (7-Methyl Guanosine Triphosphate) which helps the Ribosomes to identify the mRNA. And 3’ end of the mRNA is poly A tail (around 30-200 adenylate residues) which help shield against 3’ exonucleases)</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 139 of 320</div> <div> As RNA has differ in nucleotides sequences therefore differ in functions.</div> <div> Figure 2 RNA sequence is complementary to the DNA sequence and is translated as codons of three nucleotides</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 140 of 320</div> <div> Module 003: SIGNIFICANCES OF RNA STRUCTURE</div> <div> RNA can form 3D structures {Sarver, 2008 #5}, such structural properties helps the RNA molecule to perform different functions.</div> <div> As RNA is composed of sugars, phosphate and nucleotides and these nucleotides have ability to form hydrogen bonds.</div> <div>  A’ can make hydrogen bonds with ‘U’</div> <div>  ‘G’ makes hydrogen bonds with ‘C’</div> <div>  ‘G’ can also make hydrogen bonds with ‘U’ (Wobble Pair)</div> <div> Due to this ability of bonding RNA forms many structures and due to variety of structures RNA performs many functions in cell like:</div> <div>  DNA information transfer Figure 3 In RNA ribose is used in place of deoxyribose 3 In RNA uracil is used in place of thymine</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 141 of 320</div> <div>  Regulatory roles</div> <div>  Catalytic roles</div> <div>  Defense & immune response</div> <div>  Structure-based special roles</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 142 of 320</div> <div> Module 004: RNA FOLDING AND ENERGY FOLDING</div> <div> RNA molecules form many structures for stability and different functions. “Gibs Free Energy” (LANGRIDGE and KOLLMAN 1987) is the free energy available for RNA molecule for reactions and RNA structure formation takes place at this lower energy. Incase if RNA has two structure we can select the one with lowest energy state.</div> <div> http://chemwiki.ucdavis.edu/Core/Physical_Chemistry/Thermodynamics/State_Functions/Free_Energy/Gibbs_Free_Energy</div> <div> We can calculate the overall energy of RNA structures by summing up energies given out during the process of folding. For knowing the positive and negative values of calculations of stabilizing and destabilizing energies we may factor in ways in which RNA can be destabilized.</div> <div> Figure 4 Energy is continuously given out as the RNA molecule folds by pairing complementary bases Figure 5 calculation of stabilizing and destabilizing values</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 143 of 320</div> <div> Module 005: CALCULATING ENERGIES OF FOLDING-AN EXAMPLE</div> <div> RNA is composed of four nucleotides (A, U, C and G) and these nucleotides are attached with ribose sugar in backbone. And these nucleotides have hydrogen bonding between them. G always bond with C and Always bonds with U through hydrogen bonding and energy is released.</div> <div> That’s why RNA molecule become more stable.</div> <div> Figure 0.2nucleotides are held together by strong bonds which are created by the release of energy</div> <div> Figure 7 RNA Sequence energy</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 144 of 320</div> <div> 5 nucleotides formed H-Bonds. This bond formation released energy (-12.0 kcal/mol) RNA molecule took up a 2’ structure. Hence became more stable.</div> <div> Module 006: TYPES OF RNA SECONDARY STRUCTURES-I</div> <div> All the complimentary bases of RNA combine together to form RNA secondary structures. A simple nucleotide sequences of RNA is called as Primary structure and denoted by 1’ while when these nucleotides fold together and form a complex structure that is called secondary structure and denoted by 2’.</div> <div> The preferred structure of RNA is 2’ which has many structural patterns like Helices, Loops, Bulges and Junctions</div> <div> Figure 8 RNA sequence extends from its 5’ end to 3’ end. Upon folding, 3’ end may fold on to the 5’ end</div> <div> The first 2’ RNA structure is called helix. Unlike the DNA helix, the RNA helix is formed when the RNA folds onto itself.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 145 of 320</div> <div> The second 2’ structure is the hairpin loop</div> <div> The loop of the hairpin must at least four bases long to avoid steric hindrance with base-pairing in the stem part of the structure. Note that hairpins reverses the chemical direction of the RNA molecule.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 146 of 320</div> <div> Module 007: TYPES OF RNA SECONDARY STRUCTURES-II</div> <div> RNA 1’ structure fold the (5’-3’) ends and make RNA 2’ structure just like helix and hairpin structure.</div> <div> The third type of 2’ structure is bulge loop.</div> <div> Bulges, are formed when a double-stranded region cannot form base pairs perfectly. Bulges can be asymmetric with varying number of base pairs on one side of the loop. Bulge loops are commonly found in helical segments of cellular RNAs and used to measure the helical twist of RNA in solution. (Tang and Draper 1990)The forth type of 2’ RNA structure is interior loop.</div> <div> Interior loops are formed by an asymmetric number of unpaired bases on each side of the loop.(Turner, Sugimoto et al. 1988)</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 147 of 320</div> <div> Module 008: TYPES OF RNA SECONDARY STRUCTURES-III</div> <div> Another 2’ RNA structure is the Junction or Intersection.</div> <div> Figure 9 2' RNA structure called junction</div> <div> Junctions include two or more double-stranded regions converging to form a closed structure. The unpaired bases appear as a bulge.(Zuker and Sankoff 1984)</div> <div> Figure 10 Unpaired bases in two 2’ structures form hydrogen bonds with each other</div> <div> RNA tertiary structures are formed when RNA unpaired base bond in 2’ region bond.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 148 of 320</div> <div> Module 009: RNA TERTIARY STRUCTURES</div> <div> 2’ RNA structures is formed due to folding of nucleotide with in RNA molecule but after folding some nucleotides remain open for interaction. And they form hydrogen bonds together.</div> <div> Figure 11 Hydrogen bonding formation in open nucleotides.</div> <div> These unpaired nucleotides of 2’ structure interact with other unpaired nucleotides and form a third structure called tertiary 3’ structure. For example 4 nucleotides in hairpin loop structure does.</div> <div> The above figure:</div> <div> 1. Indicate how these 2’ structures come together</div> <div> 2. Indicate the difference between internal loop and multi loop</div> <div> 3. Indicate the yet unpaired bases</div> <div> The unpaired bases in 3’ structure remain paired by abnormal folding called (pseudoknots) but instead of pairing they remain available or pairing.</div> <div> Figure 12 pseudoknots</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 149 of 320</div> <div> Module 010: CIRCULAR REPRESENTATION OF STRUCTURES</div> <div> Tertiary or 3’ structure of RNA may form pseudoknots to detect the pseudoknots in RNA structure we need “circular plot” which is a graphical approach.</div> <div> Intersecting arcs in circular plot are the pseudoknot.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 150 of 320</div> <div> Module 011: EXPERIMENTAL METHODS TO DETERMINE RNA STRUCTURES</div> <div> RNA has 1’, 2’ and 3’ structures. 1’ has simple nucleotide sequence and 2’ has nucleotides folding and 3’ has knots.</div> <div> For measuring the RNA structure we use X-ray crystallography (Smyth and Martin 2000), which works according to the principle of diffraction. Crystallized RNA diffracts X-rays which helps estimate atomic positions</div> <div> All isotopes that contain an odd number of protons and/or of neutrons (see Isotope) have an intrinsic magnetic moment and angular momentum, in other words a nonzero spin, while all nuclides with even numbers of both have a total spin of zero. The most commonly studied nuclei are 1H and 13C, al</div> <div> Figure 13 X-ray Crystallography https://260h.pbworks.com/w/page/30814223/X%20Ray%20Crystallography</div> <div> Another method to measure the RNA structure is called as Atomic Force microscopy in this technique a laser connected to a Si3N4 piezoelectric probe scans an RNA sample. It works well in air and liquid environment.</div> <div> Figure 14 Atomic microscopy</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 151 of 320</div> <div> The third method for measuring the RNA structure is Nuclear Magnetic Resonance Imaging in this method Hydrogen atoms in RNA resonate upon placement in a high magnetic field. It Works well without crystallizing RNA</div> <div> http://www.slideshare.net/Oatsmith/13-nuclear-magnetic-resonance-spectroscopy-wade-7th</div> <div>  STORAGE OF STRUCTURES</div> <div> Reported structures are stored in online databases. Example includes RNA Bricks and RMDB etc. Bioinformaticians can refer to these databases for RNA structure studies</div> <div> RNA Bricks is a database of RNA 3D structure motifs and their contacts, both with themselves and with proteins</div> <div> Stanford University’s RNA Mapping Database is an archive that contains results of diverse structural mapping experiments performed on ribonucleic acids.</div> <div> Module 012: Strategies for RNA Structure Prediction</div> <div> RNA structure 2’ and 3’ can be measured experimentally, but RNA molecule readily degrade due to their short shelf life.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 152 of 320</div> <div> Give 1’ RNA structure creates the 2’ structure because the simple nucleotides folds and form 2’ structure. And on the base of folding we can predict the stability of the RNA molecule.</div> <div> For example.</div> <div> Figure 15 pairs represent the stability of the RNA molecule</div> <div> Maximizing the number of nucleotides can increase the structure and we have to select the structure according to the stability.</div> <div> Module 013: Dot Plots for RNA 2' Structure Prediction</div> <div> Structure measurement through experiments is slow and costly and there is maximum chances of more than one structure existence.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 153 of 320</div> <div> The dot plot method for RNA structure prediction is easy. Draw a square and partition by drawing gridlines. Put RNA sequence on top and left sides of the square. Put a “DOT” on complementary nucleotides</div> <div> For example:</div> <div> Figure 16 dot are placed at complimentary base pair place.</div> <div> Connect regions of paired nucleotides to form 2’ structures in following image.</div> <div> Figure 17 Potato Tuber Spindle Viroid</div> <div> In longest RNA nucleotides the gaps between complementary nucleotides becomes bulges and loops of the structure.</div> <div> Module 014: ENERGY BASES METHODS</div> <div> Experimental prediction of RNA structure is slow and costly that’s why a few 2’ RNA structures are reported experimentally.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 154 of 320</div> <div> While prediction we get many possible 2’ structures of RNA and for optimal structure selection we calculate their overall stability.</div> <div> Figure 18 energy table</div> <div>  STABILIZING ENERGY</div> <div> Energy table helps us to find the optimal prediction of structure because energy is released when complementary nucleotides make bonds.</div> <div>  DESTABILIZING ENERGY</div> <div> Remaining unpaired nucleotides destabilized the RNA structure in form of hairpin or bulge structure.</div> <div> Figure 19 Hairpin+IntLoop+ExtLoop+Bulge+Hairpin</div> <div>  SUM OF ENERGIES</div> <div> Sum of stabilizing and destabilizing energies can help determine the quality of a 2’ RNA structure. 2’ structure with longest coupled sequences vs. one with lowest energy</div> <div> Module 015: Zuker’s Algorithm</div> <div> Energy based methods involve evaluating the free energy structures. To compute the RNA sequence for 1’ or 2’ optimal structure prediction we use Zuker’s Algorithm.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 155 of 320</div> <div> Zuker’s Algorithm helps us to compute the stabilizing energies (-ve) and also destabilizing energies (+ve values). And also compute the sum of +ve and –ve energies.</div> <div> Figure 20 stacking energies</div> <div> Figure 21 destabilizing energies</div> <div> Figure 22 working principle of Zuker’s Algorithm (2003)</div> <div> It Compute energies of all possible 2’ structures. Generate combinations of all computed 2’ structures. Select the one with lowest energy.</div> <div> Module 016: Zuker’s Algorithm EXAMPLE</div> <div> Zuker’s Algorithm involves computing stabilizing and destabilizing energies of a 2’ structure. All possible 2’ structures are generated. The best 2’ structure is selected!</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 156 of 320</div> <div> Figure 23 Calculation of all possible structure combinations</div> <div> We need to construct all the possible combinations of nucleotides for selection of optimal 2’ RNA structure.</div> <div> Module 017: Zuker’s Algorithm – A Flow Chart</div> <div> Zuker’s Algorithm involves computing stabilizing and destabilizing energies of 2’ structure. And it also computes the overall energy by summing up the positive and negative energies.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 157 of 320</div> <div> The flow chart for energies is:</div> <div> Figure 24 flow chart</div> <div> The diagonal combination from all possible is selected with overall lowest energy.</div> <div> The two diagonals (‘D’) given above include: 1.A/U, C/G, G/C, U/G 2.G/U, U/G</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 158 of 320</div> <div> Module 018: Martinez Algorithm</div> <div> Zuker’s Algorithm involves computing stabilizing and destabilizing energies of 2’ structure. And it also computes the overall energy by summing up the positive and negative energies. Martinez Algorithm is improvement on it.</div> <div> Making combination of all possible structures is time consuming, Martinez Algorithm favors those 2’ structure which are energetically more feasible.</div> <div> Figure 25 Martinez Algorithm flow chart</div> <div> In Martinez algorithm all the 2’ structures are weighed by its stability and optimal one is sorted out. Monte Carlo methods (or Monte Carlo experiments) are a broad class of computational algorithms that rely on repeated random sampling to obtain numerical results. They are often used in physical and mathematical problems and are most useful when it is difficult or impossible to use other mathematical methods.</div> <div> And Monto Carlo method do not provide a definitive solution.</div> <div> Module 019: Dynamic Programming Approaches</div> <div> RNA sequence contains 4 type of nucleotides G/C, G/U and A/U and it may contain hundreds of nucleotides it means there is possibility of many combinations.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 159 of 320</div> <div> In 2’ RNA structure there may be large number of nucleotide sequence with large number of combinations hence it is hard to find the optimal one and for this prediction we us Dynamic Programming (DP) which breaks the larger problems into smaller one.</div> <div>  PRINCIPLE OF DYNAMIC PROGRAMMING</div> <div> For optimal structure combination selection we use the Dynamic Programming (DP) and we select the sequence of RNA nucleotides and list all the possible complementary positions for nucleotides in the given complete sequence.</div> <div> For example:</div> <div> Dynamic Programming then recombines such combinations in a process called “Traceback” to ensure that the highest coupled 2’ structure is reported</div> <div> Figure 26 all possible complementary bases combinations.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 160 of 320</div> <div> Module 020: Nussinov-Jacobson Algorithm –</div> <div> An Overview</div> <div> Nussinov-Jacobson (NJ) Algorithm is a Dynamic Programming (DP) strategy to predict optimal</div> <div> RNA 2’ structures, Proposed in 1980. Computes 2’ structures with most nucleotide coupling.</div> <div> http://ultrastudio.org/en/Nussinov_algorithm</div> <div>  HOW IT WORKS</div> <div>  Create a matrix with RNA sequences on top and right</div> <div>  Set diagonal & lower tri-diagonal to zero</div> <div>  Start filling each empty position in matrix by choosing the maximum of 4 scores</div> <div> J 1 2 3 4 5 6 7 8 9</div> <div> I G G C A A A U G C</div> <div> 1 G 0</div> <div> 2 G 0 0</div> <div> 3 C 0 0</div> <div> 4 A 0 0</div> <div> 5 A 0 0</div> <div> 6 A 0 0</div> <div> 7 U 0 0</div> <div> 8 G 0 0</div> <div> 9 C 0 0</div> <div> Figure 27 A.Note the I and J labels, B. Initialize tri-diagonal and lower tri-diagonals to zero.</div> <div> The score S ( i , j ) is the maximum of the following four possibilities</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 161 of 320</div> <div> Module 021: Nussinov-Jacobson Algorithm – The Flowchart</div> <div> NJ algorithm is actually a dynamic programming (DP) approach to predict the 2’ RNA structure. A scoring matrix is initialized to record scores in NJ Algorithm .For filling scoring matrix, the maximum score from 4 matrix positions is chosen.</div> <div> Figure 0.3 for maximum score 4 positions are used in scoring</div> <div> Figure 28 flow chart of NJ Algorithm.</div> <div> Traceback is used to report the coupling of structures in sequences.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 162 of 320</div> <div> Module 022: Nussinov-Jacobson Algorithm – EXAMPLE</div> <div> The main points to be focused in N-J Algorithm are:</div> <div>  Scoring Matrix</div> <div>  Matrix Initialization</div> <div>  Scoring method</div> <div>  The 4 different positions to be considered for calculating matrix</div> <div> Figure 29 N-J Algorithm scoring</div> <div> The matrix is filled by four different positions. Left, Bottom, Diagonal, and Left/Bottom elements. In this way all complementary nucleotides coupling is catered.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 163 of 320</div> <div> Module 023: Score Calculations & Traceback</div> <div> From four positions the score is calculated and from each position we calculate the score contribution. And maximum score is sorted out.</div> <div> Figure 30 scoring and traceback in N-J Algorithm.</div> <div> There can be many traceback. Each traceback is used to make the RNA secondary structure. And traceback with highest number of nucleotide coupling is selected.</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 164 of 320</div> <div> Module 024: Comparison of Algorithms</div> <div> RNA has three different structures 1’, 2’ and 3’. For these structures predictions there are many algorithms. But in all algorithm there are two main strategies:</div> <div> 1. Nucleotides stacking</div> <div> 2. Energy minimization</div> <div>  ENERGY BASED ALGORITHM.</div> <div> Zuker’s Algorithm involves energy minimization. It is updated version and incorporate the phylogenetic information. It is improved. Overcomes the pseudoknots assumes them and accommodate them. And this algorithm helps to predict the structures of RNA based on nucleotides.</div> <div>  NUCLEOTIDES STACKING ALGORITHM.</div> <div> NJ’s Algorithm comes under this category. It involves the maximizing the nucleotides pairing. Traceback helps to find best 2’ structure.</div> <div> It predict the 75% accurate 2’ structure. Because there may be more than two equal scores as it is calculated from four different positions. To get best results we need to combine the stacking and energy minimization methods together.</div> <div> For further improvements in results we take help from:</div> <div>  Sequences</div> <div>  Comparison</div> <div>  Nucleotide</div> <div>  Covariance analysis</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 165 of 320</div> <div> Module 025: WEB RESOURCES-I</div> <div> For prediction of 1’ and 2’ structure of RNA we use different algorithm like Zuker’s, Martinez and N-J. Online tools also. The mfold web server is one of the oldest web servers in computational molecular biology. Mfold is upgraded version of Zuker’s algorithm. MFOLD is computationally expensive and can give results for 1’ and 2’ structures that have sequences less than 8000 nucleotides.</div> <div> Figure 31 http://unafold.rna.albany.edu/?q=mfold</div> <div> Figure 32 http://unafold.rna.albany.edu/?q=mfold/RNA-Folding-Form</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 166 of 320</div> <div> Figure 33 http://unafold.rna.albany.edu/?q=mfold/Structure-display-and-free-energy-determination</div> <div> MFOLD helps fold an RNA nucleotide sequence into its possible 2’ structures. MFOLD gives out several structures along with their energetic stability!</div> <div> Handout Chapter 04: RNA Secondary Structure Prediction</div> <div> Page 167 of 320</div> <div> Module 026: WEB RESOURCES-II</div> <div> RNA nucleotides folds to form 2’ structure from simple portion of 1’ nucleotides. For example CUUCGG occurs a wide variety in RNA and it mostly forms the stable hairpin loop. So we can make the list of all likely 2’ structures arising from 1’.</div> <div> Figure 34 http://www.rnasoft.ca/strand/</div> <div> Figure 35 http://iimcb.genesilico.pl/rnabricks</div> <div> RNA 1’ folds and makes RNA 2’ structure and this online database is established for 2’ RNA structure and it act as dictionary for 2’ RNA structure.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 168 of 320</div> <div> Chapter 5 - Protein Sequences</div> <div> Module 001: FROM DNA/RNA SEQUENCES TO PROTEIN</div> <div> We are aware that DNA has four nucleotides bases (A, C, T & G). RNA contains (A, C, U & G). And protein contain 20 different amino acids. DNA to RNA then Protein is called as central dogma. Which includes translation, transcription and protein modifications.</div> <div> A set of three nucleotides called codon, codes the information for specific amino acids in protein synthesis.</div> <div> Figure 0.1 amino acids letters information according to codons</div> <div> Figure 0.2 codon (set of three nucleotide) codes for specific amino acid.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 169 of 320</div> <div> Codons select the amino acids and ribosomes make the protein by polymerization process and these nucleotides coil together to form 3D structure.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 170 of 320</div> <div> Module 002: CODING OF AMINO ACIDS</div> <div> Nucleotides (A, G, C, and T) make set of three called codons for amino acid selection in protein synthesis. More than one codon can code for same amino acids as there are 20 amino acids involved in protein synthesis.</div> <div> Figure 3 coding of amino acids</div> <div> Figure 4 Start Codon ATG and Stop Codon TAG, TGA or TAA</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 171 of 320</div> <div> Module 003: OPEN READING FRAMES</div> <div> Codons codes information for amino acid and there are three stop codons and one start codon. For the valid open reading frame it must have longest sequence. In molecular genetics, an open reading frame (ORF) is the part of a reading frame that has the potential to code for a protein or peptide. An ORF is a continuous stretch of codons that do not contain a stop codon (usually UAA, UAG or UGA).</div> <div> https://en.wikipedia.org/wiki/Open_reading_frame</div> <div> Figure 5 ORF 1 is valid, as it is the longest</div> <div> There is online tool from which we can find ORF in any sequence.</div> <div> Figure 6 NCBI, ORF Finder</div> <div> Six ORF exist in any DNA sequence and longest one is marked and first stop codon will be marked end of the protein.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 172 of 320</div> <div> Module 004: ORF Extraction – A Flowchart</div> <div> Codons of 3 nucleotides code for each Amino Acid. There are 1 start and 3 stop codons. Selection of ORF is based on its length if it the longest one from others than it would be suitable for protein synthesis reaction.</div> <div> Figure 7 ORF extraction flowchart</div> <div> Both reverse and forward RNA sequences are considered which may have many ORF and selection is based upon longest protein sequences having.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 173 of 320</div> <div> Module 005: SEQUENCING PROTEINS</div> <div> Given the DNA/RNA sequence, ORFs can be extracted and protein sequence can be determined. But there are chances that protein may be unknown, that’s why we use Adam degradation method in protein sequencing. Edman degradation, developed by Pehr Edman, is a method of sequencing amino acids in a peptide. In this method, the amino-terminal residue is labeled and cleaved from the peptide without disrupting the peptide bonds between other amino acid residues.it starting from the N-terminal and removing one amino acid at a time</div> <div> Figure 8 Mechanism</div> <div> https://en.wikipedia.org/wiki/Edman_degradation</div> <div> Cyclic degradation of peptides by Phenyl-iso-thio-cyanate (PhNCS). PhNCS attaches to the free amino group at N-terminal residue. 1 amino acid is removed as a PhNCS derivative.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 174 of 320</div> <div> Figure 9 working of Edam degradation</div> <div>  DRAWBACKS</div> <div>  It is restricted to chain of 60 residues.</div> <div>  It is very time consuming process 40-50 amino acids per day.</div> <div> Modern techniques for this is Tandem mass spectrometry.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 175 of 320</div> <div> Module 006: Application of Mass Spectrometry in Protein Sequencing</div> <div> Edam Degradation methods helps us to sequence the protein which is unknown. But it is restricted to 60 amino acids only.</div> <div> Protein can be charged with electrons or protons and if moving charges are placed in between the magnetic field they get deflected. And their deflection is proportional to their momentum.</div> <div> Figure 10 Moving charged particles in a magnetic field</div> <div> Where:</div> <div>  F is the force applied to the ion</div> <div>  m is the mass of the particle,</div> <div>  a is the acceleration</div> <div>  Q is the electric charge,</div> <div>  E is the electric field</div> <div>  v × B is the cross product of the ion's velocity and the magnetic flux density.</div> <div> Figure 11 equation for MS application in protein sequencing</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 176 of 320</div> <div>  COMPUNENTS</div> <div>  Sample Injection</div> <div>  Ionization Source</div> <div>  Mass Analyzer</div> <div>  Ion Detector</div> <div>  Spectra search using computational tools</div> <div>  CONCLUSION</div> <div> Charged proteins can be set into motion within a magnetic field. Their deflections accurately correspond to their molecular mass. Deflections can be measured (hence protein’s mass)</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 177 of 320</div> <div> Module 007: Techniques for MS Proteomics</div> <div> MS proteomics works on the principle of protein ionization which are placed in very high magnetic field. Each protein deflect to its proportion which is equal to its molecular weight in this way molecular mass is measured. The protein mass of unknown protein is compared with the masses of proteins in database and matching one is selected.</div> <div> Example for protein sequences database is uniProt, swissprot etc.</div> <div> Proteins are measured and sequenced if are unknown than matched with existing database if matched than are shortlisted.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 178 of 320</div> <div> Module 008: Types of MS-based proteomics</div> <div> Proteins can be sequenced by Edam’s degradation and Mass spectrometry. MS based proteomics helps us to sequence the larger and bigger proteins more quickly.</div> <div> Following steps are involved in MS:</div> <div>  Separation</div> <div>  Ionization</div> <div>  Mass analysis</div> <div>  Detection</div> <div> Two methodologies are involved</div> <div> 1. Bottom up proteomics</div> <div> 2. Top down proteomics</div> <div> Bottom up proteomics measures the peptide masses produced after protein enzymatic digestion. And Top down proteomics measures the intact proteins followed by peptides after fragmentation.</div> <div>  BOTTOM UP PROTEOMICS</div> <div> In this methodology the protein complex is treated with site specific enzymes which cleaves them into amino acid residue and resultant peptides are measured for their masses. One peptide is selected at one time for processing and when all are processed than protein search engine is used for matches.</div> <div>  TOP DOWN PROTEOMICS</div> <div> In this methodology proteins are ionized and measured for their masses and one protein is mass selected at a time for fragmentation. And resultant peptide fragments are measured for mass.</div> <div> We can say that bottom up proteomics deals with peptides while top down proteomics can handle the whole protein.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 179 of 320</div> <div> Module 009: BOTTOM UP PROTEOMICS</div> <div> There are two types of proteomics protocol that are usually employed.</div> <div> 1. Bottom up proteomics</div> <div> 2. Top down proteomics</div> <div> PROTOCOL</div> <div> 1. Sample containing the mixture of protein from cells and tissues is obtained.</div> <div> 2. Enzymes such as trypsin is use to cleave the proteins.</div> <div> 3. Enzyme cleaves the amino acids at specific sites of amino acid.</div> <div> 4. Several peptides are formed when protein is cleaved.</div> <div> 5. Number of peptide depends upon the number of sites where enzymes cleaved the protein. For example trypsin cleaves the protein at lysine (k)</div> <div> 6. Mass of each peptide is measured.</div> <div> 7. One peptide is selected at a time.</div> <div> 8. Different enzyme is use to cleave the protein at different site.</div> <div> 9. This process keep going until the possible number of peptides are formed or searched.</div> <div> 10. Peptides are searched in data base and matched.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 180 of 320</div> <div> Module 010: Two Approaches for Bottom up Proteomics (BUP)</div> <div> There are two approaches for bottom up proteomics.</div> <div> 1. Peptide Mass Fingerprinting.</div> <div> 2. Shotgun Proteomics</div> <div> Figure 0.33 Peptide mass fingerprinting</div> <div> Figure 14 Shotgun Proteomics</div> <div> Shotgun Proteomics digest the whole protein and mix first and compared with database. And peptide mass finger printing involves in protein separation followed by single protein’s peptide analysis.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 181 of 320</div> <div> Module 011: TOP DOWN PROTEOMICS</div> <div> Bottom up proteomics identifies the proteins by cleaving them into segments at specific sites and was not suitable to measure the direct protein masse.</div> <div>  PROTOCOL</div> <div> 1. Sample containing the protein mixture from cells and tissues is obtained</div> <div> 2. The entire protein is mixed and analyzed for masses.</div> <div> 3. The list of masses is obtained.</div> <div> 4. TDP Measures all post translational masses of protein.</div> <div> 5. After MS1 one protein is selected at a time and fragmented to obtain its peptides.</div> <div> 6. The process is repeated many times.</div> <div> Comparison is done from protein database uniProt and swissProt.</div> <div> TDP also measure the masses of intact proteins and masses of post transcriptional changes.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 182 of 320</div> <div> Module 012: PROTEIN SEQUENCE IDENTIFICATION</div> <div> Mass spectrometry helps us to measure the molecular weight of proteins and peptides, but</div> <div> several proteins can have same masses to identify them we follow the flow chart of following</div> <div> techniques.</div> <div> Figure 0.45 protein sequence identification flowchart</div> <div> The flowcharts discussed above can help us arrive at the sequence of the protein in question.</div> <div> Scoring schemes are required to quantitatively represent the quality of results</div> <div> P62837...51</div> <div> P62837...51</div> <div> Q6ZPJ3...40</div> <div> P21734...33</div> <div> Q02159...21</div> <div> P62837...51</div> <div> Filter Protein</div> <div> Database</div> <div> In Silico Fragmentation</div> <div> of Candidate Proteins</div> <div> Ionization</div> <div> EST Determination MS Spectra</div> <div> Protein Scores</div> <div> =</div> <div> Matching of Experimental</div> <div> and Insilco Peak List</div> <div> Protein Complex Separation</div> <div> + +</div> <div> +</div> <div> +</div> <div> +</div> <div> Mass Spectrometer</div> <div> Fragmentation</div> <div> MS2</div> <div> Post</div> <div> Translational</div> <div> Modifications</div> <div> Compare Theoretical</div> <div> Masses with</div> <div> Experimental</div> <div> Masses</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 183 of 320</div> <div> Module 013: PROTEIN IONIZATION TECHNIQUES</div> <div> Protein ionization is used in Mass spectrometry based on proteomics protocols. Ionization involves loading of proton in protein or removal of protein. Ionizations can increase or decrease the mass of protein or peptide.</div> <div>  SALIENT IONIZATION</div> <div> Is the technique which include Matrix Assisted Laser Desorption Ionization MALDI) & Electro Spray Ionization (ESI)</div> <div> For example:</div> <div>  MALDI</div> <div> In this technique one proton is added to protein or peptide and the molecular weight is increases by one and Mass spectrometry reports the molecule at +1.</div> <div>  ESI</div> <div> ESI adds many protons to protein or peptides and molecular weight is increased by the number of protons added. But it is difficult in ESI to find the molecule with +1.</div> <div>  EXAMPLE</div> <div> Figure 0.56 resolving multiple charges</div> <div> MS data from MALDI ionization is easier to handle as the product ions masses are mostly at “1+mass”. ESI is difficult to use as it does not easily give away the +1 charged ion</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 184 of 320</div> <div> Module 014: MS1 & INTACT PROTEIN MASS</div> <div> When we ionize the protein, it can be deflected by a magnetic field in proportion to its mass and the mass of protein can be measured by spectrometry.</div> <div> Figure 0.67 MS1 Schematic (Image courtesy Wikipedia)</div> <div> Mass/charge helps us to calculate the mass of protein, “Mass Select” can help to select specific MS1 for further analysis.</div> <div> MS1 results the intact masses of the peptides.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 185 of 320</div> <div> Module 015: SCORING INTACT PROTEIN MASS</div> <div> MS1 helps us to obtain the intact masses of precursor molecules which depend upon the proteomics and protocol applied. Protein masses reported by MS1 are matched with protein database, but before match the masses are converted into +1 of all molecule.</div> <div>  SCORING</div> <div> We can score each protein in the way that it get maximum score and low quality matches should get low scores.</div> <div> After filtering the multiple charges we get the only the peaks having charge 1. And after this filter we compare it with protein data base.</div> <div>  SCORING SCHEME</div> <div> Example: A Protein Sequence from Database “MQLF”</div> <div> http://web.expasy.org/compute_pi/</div> <div> All experimental masses are compared with theoretical masses of database and mass is selected on the base of closeness.</div> <div> 𝑀𝑆𝑐𝑜𝑟𝑒 √ 𝑀𝐸𝑥𝑝−𝑀𝑇 2 MW: 537.67</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 186 of 320</div> <div> Module 016: PROTEIN FRAGMENTATION TECHNIQUES</div> <div> We compare the experimental mass with theoretical data base mass of protein and on base of closeness we rank or score it.</div> <div> If several proteins have same score than selection is done by using another technique protein fragmentation. We fragment the protein or peptide and ionize it, it helps us to measure the fragment masses as the same ways as their precursor.</div> <div> There are different techniques for protein fragmentation.</div> <div>  Electron Capture Dissociation (ECD)</div> <div>  Electron Transfer Dissociation (ETD)</div> <div>  Collision Induced Dissociated (CID)</div> <div> Each fragmentation technique gives result of specific type of fragments.</div> <div> ECD gives out ‘C’ and ‘Z’ ions. CID gives out ‘B’ and ‘Y’ ions, etc.</div> <div> Figure 0.78 natural peptide of four residue</div> <div> If we can measure the mass of fragments using MS, Calculate the theoretical mass of the fragments. Then, we can award score on the basis of the similarity of experimental and theoretical mass.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 187 of 320</div> <div> Module 017: TANDEM MS</div> <div> Intact masses can measure the intact proteins or peptides. And this can be followed up by their fragmentation in MS chamber.</div> <div> Tandem MS can be extended to the fragments of the intact fragment. All you need is the MS instrument capability to, (i) select fragment’s mass range. (ii) Fragment the precursor fragment.</div> <div> Tandem MS helps us to measure masses of fragments. By this scoring and protein identifications so easy.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 188 of 320</div> <div> Module 018: MEASURING EXPERIMENTAL FRAGMENT’S MASS</div> <div> In MS1, the molecular weight of intact sample molecule is measure and then intact molecule is fragmented in two afterward, these two fragments are measured by MS or MS2</div> <div> FRAGMENTATION TECHNIQUES AND MOLECULAR WEIGHT</div> <div> Fragmentation techniques include ECD, CID etc. intact molecule fragmentation splits the molecule into two parts.</div> <div> FRAGMENT MASS</div> <div> Mass of fragment is produced by MS2 deepening upon the technique because each techniques splits the protein or peptide at different location.</div> <div> Figure 0.89 Masses after Fragmentation by ECD, CID & ETD</div> <div> Experimental mass reported from MS2 is matched with theoretical peptides of candidate proteins (from DB). Score is awarded on the basis of the closeness between experimental and theoretical masses.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 189 of 320</div> <div> Module 019: Calculating Theoretical Fragment's Mass</div> <div> After measuring the mass of intact molecule from MS2 we compare that mass with theoretical mass of databased proteins.</div> <div> Figure 20 Masses after Fragmentation by CID</div> <div> Figure 0.91 Masses after Fragmentation by ECD</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 190 of 320</div> <div> Figure 22 Masses after Fragmentation by ECD</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 191 of 320</div> <div> Module 020: PEPTIDE SEQUENCE TAG</div> <div> Peptide sequence tag are the sequence of peptide which are produced after MS2. We can obtain the sequence of peptide through variation in fragmentation site.</div> <div> Figure 23 variation sites</div> <div> Precursor proteins or peptides fragmentation leads to formation of multiple ions of the same fragment type. However, fragments have variation in their molecular weights due to variation in site of fragmentation</div> <div> Fragmentation at consecutive sites leads to a mass difference equal to that of a single amino acid. Such consecutive peaks can reveal partial peptide sequence tags</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 192 of 320</div> <div> Module 021: Extracting Peptide Sequence Tags</div> <div> PSTs are formed due to sequential cleavage of precursor protein/peptide’s backbone.</div> <div> Figure 24 peptide sequencing tagging</div> <div> Peptide sequence tags can be extracted from peak list iteratively. A high quality mass spectrum will produce large number of PSTs. The bigger the peptide sequence tags, the better!</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 193 of 320</div> <div> Module 022: Using Peptide Sequence Tags in Protein Search</div> <div> PSTs provide clues of the precursor protein/peptides sequence. Consider that we extract the following PSTs: M, MQ, QV etc. Search protein sequence database (e.g. Uniprot, Swissprot)</div> <div> Sample sequence in protein DB</div> <div> >>sp|Q6GZ4X|0X1R_FRG3G Putative transcription factor 0X1R OS=Random virus 3 (isolate Goorha) GN=FV3-0X1R PE=4 SV=1</div> <div> MAFSAEDVLKEYDRRRRMEALLLSLYYPNDRKLLDYKEWSPPRVQVECPKAPVEWNNPPSEKGLIVGHFSGIKYKGEKAQASEVDVNKMCCWVSKFKDAMRRYQGIQTCKIPGQVLSDLDAKIKAYNTVEGVEGFVRYSRVTKQHVAAFLKELRHSKQYENVNLIHYILTDKRVDIQHLEKDLVKDFKALVESAHRMRQGHMQNVKYILYQLLKKHGHGPDGPDILTVKTGSMQLYDDSFRKIYTDLGWKFTPL</div> <div> For all the proteins in the database, we find out which PSTs exist in which proteins. The protein reporting the most PSTs is more probable to be the precursor protein.</div> <div> If many PSTs report the same number or protein report the longer PSTs than through scoring we find the greater number. After extracting the PSTs we search the entire database for protein who report it.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 194 of 320</div> <div> Module 023: Scoring Peptide Sequence Tags</div> <div> According to scoring scheme if a candidate protein matches ‘n’ PSTs, then its score can be given by:</div> <div> Additionally, if we include RMSE to the scoring system, then it can highlight better PST matches. And RMSE is the root mean square error.</div> <div> Figure 25 root mean square error</div> <div> RMSE for a sequence tag ‘i’ of length ‘n’?</div> <div> So, the updated relationship is:</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 195 of 320</div> <div> Module 024: In silico Fragment Comparison</div> <div> MS1 reports the intact mass of molecules (proteins or peptides) in the sample. Intact mass can be compared with every protein’s mass in database to identify the molecule in the sample.</div> <div> Incase multiple candidate proteins are reported, MS2 can be performed. MS2 helps measure fragment peptide masses. MS2 data can be used to extract peptide sequence tags</div> <div> If the protein identification is still not conformed than each experimentally reported MS2 fragment is compared with the in silico spectra of proteins from database.</div> <div> Fragmentation techniques determine product ions e.g. ECD -> c/z and CID -> b/y ions etc. With known fragment types, we can compute the MW of all possible protein fragments</div> <div> For obtaining all possible theoretical fragments in a protein, we need to compute the MWs of each fragment individually</div> <div> Consider a random protein sequence from DB:</div> <div> Figure 26 random protein sequence</div> <div> Matching experimental fragments with in silico fragments is the final resort in protein search and identification.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 196 of 320</div> <div> Module 025: In silico Fragment Comparison and scoring</div> <div> Experimental MS2 can be compared with the in silico spectra of protein from database.</div> <div>  Count the matches between in silico and in-vitro peaks.</div> <div>  Give an equivalent score to candidate protein.</div> <div>  Weigh each of the aforementioned match by the mass error</div> <div>  Accumulate the score</div> <div> With “all possible” fragments in in silico spectrum, and “reported” fragments in experimental spectrum, we can match and rank.</div> <div> Scoring scheme should also consider the errors in peak matching</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 197 of 320</div> <div> Module 026: Protein Sequence Database Search Algorithm</div> <div> MS1 and MS2 provide us with a host of data towards enabling us in identifying unknown proteins. A step by step approach combining MW, PSTs and insilico spectral matching is required.</div> <div> Figure 27 protein sequencing flowchart</div> <div> Integrating MW, PST and insilico comparison algorithms in a workflow can help create a composite protein search engine. A composite scoring system is also required for this search engine.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 198 of 320</div> <div> Module 027: Integrative Scoring Schemes</div> <div> Three individual scores can be obtained:</div> <div>  MW Match score</div> <div>  PST Match score</div> <div>  In vitro & In silico Match score</div> <div> For overall cumulative score computation:</div> <div>  We simply sum the scores up (a linear function).</div> <div>  Weigh each scoring component up by respective RMSE before summing them up</div> <div>  Complex non-linear functions integrate the scoring components in Mascot etc.</div> <div>  Highly proprietary for commercial proteomics software are used.</div> <div> Composite scoring schemes are needed to combine scores coming in from multiple criteria. The ability of a scoring scheme to better isolate true positives from false positives is important.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 199 of 320</div> <div> Module 028: LARGE SCALE PROTEOMICS</div> <div> Peptide mixtures in bottom up proteomics are very complex. Tryptic peptides may reach up to an order of 300,000–400,000.In whole proteome samples, protein count may be over 10,000. Experiments have shown that it is difficult or even impossible to analyze all these peptides in a single analysis, as the mass spectrometer is essentially overwhelmed.</div> <div> Over half a million peptides reported in a typical LSP experiment are redundant.</div> <div> If we could find a unique peptide for a protein, that would make sequence coverage suffer and we have to strike a compromise between sequences coverage and sample coverage.</div> <div>  TECHNIQUE</div> <div> One way forward would be to transfer peptides to the MS chamber in a step-by-step manner. However, this imposes a precondition that a peptide is not selected earlier as well (i.e. more than once)</div> <div>  STEP BY STEP TECHNIQUE</div> <div> 1. The instrument alternates between MS and MS/MS modes.</div> <div> 2. Three most intense peaks are chosen for MS/MS analysis.</div> <div> 3. After the initial MS scan, an MS/MS spectrum from peptide A is obtained by selectively fragmenting this mass only.</div> <div> 4. Next, a spectrum for peptide B is produced, followed by a recording of the MS/MS spectrum for peptide C.</div> <div> 5. After these three fragmentation spectra have been obtained, a new MS scan is started.</div> <div> From this scan, three more peptides A B C are selected for fragmentation and the cycle starts over again.</div> <div> The number of MS1/2 scans can be limited by carefully selecting the peptide peaks. Once the intense peptides are identified, next batch of peptides is chosen for MS2.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 200 of 320</div> <div> Module 029: PROTEOMICS DATA FILE FORMATS</div> <div> Mass spectrometer is used to measure mass/charge ratio of ionized proteins and peptides. Data output from the MS comprises of m/z ratios and intensities of each molecule that is measured.</div> <div> Followings are the formats for proteomics data:</div> <div> Figure 28 Formats used for proteomics data</div> <div> OPEN FORMTAS:</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 201 of 320</div> <div> mzXML (tools.proteomecenter.org/mzXMLViewer.php)</div> <div> MGF (proteomicsresource.washington.edu/mascot/help/data_file_help.html)</div> <div> Multiple MS data formats exist. Proprietary formats exist which come implemented as software with hardware. Also, open software standards exist for interoperability etc.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 202 of 320</div> <div> Module 030: RAW FILE FORMAT</div> <div> Mass spectrometer outputs data with ionic mass/charge ratios & respective ion intensities. RAW file is a format in which an instrument outputs data in binary form.</div> <div> Figure 29 Raw file formats</div> <div> Figure 30 list of tools for raw data processing</div> <div> Multiple RAW file formats are prevailing in the industry. Each vendor has its own unique RAW file format. You can convert proprietary formats into open formats</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 203 of 320</div> <div> Module 031: MGF FILE FORMAT</div> <div> MGF – Mascot Generic Format. MGF is a simple human-readable format for MS/MS data developed by Matrix Science. Mascot Search Engine available at this link online. http://www.matrixscience.com/search_form_select.html</div> <div> http://www.matrixscience.com/help/data_file_help.html</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 204 of 320</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 205 of 320</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 206 of 320</div> <div> Module 032: OPEN MS DATA FORMAT</div> <div> Mass spectrometer outputs data with mass/charge ratios & respective ion intensities. RAW file formats are specific to each instrument and each vendor has its own unique file format. Once an instrument is upgraded, data output from the instrument is also changed. Hence the underlying RAW file format needs to be upgraded as well.</div> <div> NEED</div> <div> Proprietary RAW formats are binary formats which are difficult to read and parse. If you have the software from the maker of the MS then you can read the RAW data file as well.</div> <div> SOLUTION</div> <div>  mzData was developed by HUPO-PSI</div> <div>  mzXML was developed at the Institute for Systems Biology</div> <div>  To combine them, a joint venture produced mzML</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 207 of 320</div> <div> Figure 31 Formats used for open use</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 208 of 320</div> <div> Several software exist for converting RAW file formats into open software formats. Each open format has its own unique advantages. mzXML and MGF formats are most frequently used</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 209 of 320</div> <div> Module 033: Online Proteomics Tools - MASCOT</div> <div> Matrix Science developed an online Bottom up Proteomics Search Engine. “MASCOT”. Mascot can search peptide mass fingerprinting and shotgun proteomics dataset</div> <div> Figure 32 mzML http://tools.proteomecenter.org/software.php</div> <div> Figure 33 http://www.matrixscience.com/search_form_select.html</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 210 of 320</div> <div> Mascot is the most widely used online search tool for proteomics data. However, it lacks a batch processing mode. Also, it does not cater for top-down proteomics data.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 211 of 320</div> <div> Module 034: Online Proteomics Tools – ProSight PTM</div> <div> Kelleher et al have developed an online Top down Proteomics Search Engine. “Prosight PTM”. ProsightPTM searches top down proteomics data and reports the precursor protein</div> <div> https://prosightptm.northwestern.edu/about_retriever.html</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 212 of 320</div> <div> https://prosightptm.northwestern.edu/about_retriever.html</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 213 of 320</div> <div> ProSight PTM is the state of the art in top down proteomics search. Using Prosight PTM, post-translational modifications can be accurately identified.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 214 of 320</div> <div> Module 035: Example Case Study - I</div> <div> For case study we follow some steps:</div> <div> Step 1 – Monoisotopic Peak Detection</div> <div> Natural elements occur in multiple isotopes. Isotopes differ in their masses.The abundance of each isotopic variant is unique.</div> <div> Figure 34 Isotopic variants of natural elements</div> <div> TYPES OF MASSES</div> <div>  Nominal Mass</div> <div>  Monoisotopic Mass</div> <div>  Average Mass</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 215 of 320</div> <div> Figure 35 Detecting Monoisotopic Peaks</div> <div> Figure 36 Detecting Monoisotopic Peaks</div> <div> MS1 data reports the isotopic distribution of intact molecule’s mass. Monoisotopic mass value has to be selected from this mass distribution. This value is the highest mass value in the distribution</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 216 of 320</div> <div> Module 036: Example Case Study – II</div> <div> The first step in protein identification and characterization using mass spectrometry involves intact protein/peptide mass measurement. Next, we fragment the protein. A protein or peptide backbone may be fragmented anywhere along the peptide backbone.</div> <div> This results in formation of two fragments i.e. N-term fragment and C-Term fragment.</div> <div> For possible fragments let’s take an example protein with 100 residues. Such a molecule’s backbone can be fragmented at 100 different locations. The total number of possible fragments is then 200</div> <div> TANDEM MS</div> <div> The mass of 200 fragments can then be measured by using an MS again. The necessary condition for this measurement is that all 200 fragments are ionized.</div> <div> To ensure that all fragments of precursor molecule are also charged, we can use Electrospray ionization (ESI).ESI induces multiple charges on the intact molecule</div> <div> Role of Electrospray Ionization</div> <div> Since ESI induces multiple charges on the precursor molecule, there is a good chance that upon precursor’s fragmentation, each fragment will have a portion of the charge. ESI allows for production of multiple charged ions. Tandem MS helps measure molecular weights of ionized fragments</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 217 of 320</div> <div> Module 037: Example Case Study – III</div> <div> Tandem MS helps measure the mass of the fragments Those fragments which differ from each other by one amino acid’s mass can provide clues on the sequence of proteins</div> <div> Figure 37 Example peptide sequence tags</div> <div> Peptide sequence tags help derive clues about the sequence of precursor proteins/peptides. The short peptide sequences help us in shortlisting candidate proteins from the database.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 218 of 320</div> <div> Module 038: Example Case Study – IV</div> <div> MS1 helps measure the intact mass of proteins/peptides. A list of candidate proteins/peptides can be formed by comparing MS1 mass to the mass of proteins/peptides in the database. MS2 or Tandem MS was performed after fragmentation of intact proteins.MS2 helped extract peptide sequence tags from MS2 data. Candidate proteins can be further shortlisted by the PSTs</div> <div> Exhaustive matching of all MS2 peaks with the theoretical fragments of candidate proteins. The set of theoretical fragments contains all possibilities of fragmentation</div> <div> Theoretical vs experimental fragments comparison helps as the third stage for shortlisting candidate protein list. This shortlisting will help you arrive at a small number of proteins</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 219 of 320</div> <div> Module 039: Example Case Study – V</div> <div> MS1 and MS2 provide mass of intact molecules and its fragments. This information helps filter proteins from protein database. For a quantitative measure, scoring scheme is required.</div> <div> Figure 38 Example intact protein mass score</div> <div> Figure 39 Example peptide sequence tags</div> <div> Three scoring schemes can be applied to score the match at each stage of protein search. These scoring elements can be integrated to arrive at an overall candidate protein score.</div> <div> Handout Chapter 05: Protein Sequences</div> <div> Page 220 of 320</div> <div> Module 040: Example Case Study – VI</div> <div> Comparisons can be performed at various levels of information. These include MS1, MS2, PSTs and theoretical fragments comparison. Integrated scoring schemes couple these factors.</div> <div> For comprehensive scoring</div> <div> A comprehensive scoring scheme can combine all the scores. Several optimizations can be undertaken on the scoring scheme to further improve protein identification</div> <div> 𝑺𝒄𝒐𝒓𝒆 𝑺𝒄𝒐𝒓𝒆 𝑴𝑾⬚+𝑺𝒄𝒐𝒓𝒆 𝑷𝑺𝑻⬚+𝑺𝒄𝒐𝒓𝒆 𝑬𝒙𝒑<>𝑻𝒉𝒓⬚ Simply sum the scores up (a linear function) 𝑺𝒄𝒐𝒓𝒆⬚ 𝑺𝒄𝒐𝒓𝒆 𝑴𝑾𝑬𝑴𝑾+ 𝑺𝒄𝒐𝒓𝒆 𝑷𝑺𝑻𝑹𝑴𝑺𝑬𝑷𝑺𝑻𝒎𝒊=𝟎+ 𝑺𝒄𝒐𝒓𝒆 𝑬𝒙𝒑<>𝑻𝒉𝒓𝑬𝑬𝑿𝑷<>𝑻𝒉𝒓𝒏𝒊=𝟎</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 221 of 320</div> <div> Chapter 6 - Protein Structures</div> <div> Module 001: PROPERTIES OF AMINO ACIDS-I</div> <div> Proteins are made by polymerization of amino acids on ribosomes and proteins properties are linked to the properties of amino acids. There are 20 amino acids in nature each has different chemical composition and that’s why each protein is different from other.</div> <div> Figure 0.1 chemical structure of amino acid</div> <div> Amino acid have three groups, hydroxyl group, Amine group and R group. The R group is representing any group.</div> <div> Figure 0.2 periodic chart of amino acid</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 222 of 320</div> <div> During polymerization of amino acids the water is formed and amino acids attached with each other.</div> <div> Figure 0.3 polymerization of amino acids</div> <div> Amino acids have unique properties such as polarity, charge states and interactions with water. Each of these properties describes the overall characteristic of an amino acid.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 223 of 320</div> <div> Module 002: PROPERTIES OF AMINO ACIDS-II</div> <div> Amino acids have characteristics like polarity, hydrophobicity, and charge states. These characteristics are governed by the elemental composition of an amino acid’s side chain (R group).</div> <div> Figure 4 R group in amino acid</div> <div> HYDROPHILIC AMINO ACIDS</div> <div> Since H and C introduce very little dipole moments in hydrophobic amino acids, these amino acids are non-polar. Hydrophobic amino acids are mostly found at the inside of folded proteins. Hydrophilic group contain the chain of C and H group in their R group.</div> <div> Figure 5 hydrophilic group</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 224 of 320</div> <div> POLAR AMINO ACID</div> <div> These amino acids are polar but are not charged i.e. no net charge on the amino acid. Prefer to reside / interact with aqueous environments. Mostly found at the surface of folded proteins.</div> <div> Figure 6 Polar amino acids</div> <div> Amino acids have unique properties such as polarity, charge states and interactions with water. Each of these properties describes the overall characteristic of an amino acid.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 225 of 320</div> <div> Module 003: PROPERTIES OF AMINO ACIDS-III</div> <div> Some amino acids are positively charged and some have negative charge.</div> <div> Figure 7 positively charged amino acids</div> <div> Figure 8 negative charge amino acid</div> <div> Upon polymerization of amino acids into polypeptide chains, charged amino acids get neutralized. At pH=7, five amino acids are charged, 2 negatively and 3 positively.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 226 of 320</div> <div> Module 004: PROPERTIES OF AMINO ACIDS-IV</div> <div> Some amino acids are positively charged and some have negative charge. pK is the values for an amino acid is the pH at which exactly half of the chargeable group is charged.</div> <div> If pH < pK for an amino acid, the amine side chains gain a proton (H+) and become positively charged, hence basic.</div> <div> If pH > pK for an amino acid, the carboxyl side chains loses a proton (H+) and become negatively charged, hence acidic.</div> <div> -ve +ve</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 227 of 320</div> <div> Figure 9 properties of amino acids according to pK and PH.</div> <div> Depending on the pH, an amino acid may become charged. This may be positive or negative depending on the amino acid.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 228 of 320</div> <div> Module 005: PROPERTIES OF AMINO ACIDS-V</div> <div> Amino acids may be charged depending on pH. This depends on the charge acceptance or donation from within an amino acid. Additionally, amino acids have structures as well.</div> <div> Figure 10 Aliphatic Amino Acids (Non polar C and H chains)</div> <div> Figure 11 Aromatic R groups</div> <div> Side chain also impact some properties. Side chains comprising merely of Carbon and Hydrogen are:</div> <div>  Chemically inert,</div> <div>  Poorly soluble in water</div> <div> However, side chains containing organic acids are very different. They are chemically reactive and Soluble in water. Elemental composition plays a very important role in determining properties of amino acids. Solubility and reactivity are key factors participating in protein folding.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 229 of 320</div> <div> Module 006: STRUCTURAL TRAIT OF AMINO ACID-I</div> <div> Amino acids have several properties such as charge state, polarity and hydrophobicity. It is important to note that the physical size of each amino acid also varies.</div> <div> EXAMPLE-1: Glycine</div> <div> Glycine residues increase backbone flexibility because they have no R group (only an H), hence agile.</div> <div> EXAMPLE-2: Proline</div> <div> Proline residues reduce the flexibility of polypeptide chains. Proline cis-trans isomerization is often a rate-limiting step in protein folding.</div> <div> Figure 12 cis and Tran’s form of proline</div> <div> EXAMPLE-3: Cystine</div> <div> Cysteines cement together by making disulfide bonds to stabilize 3-D protein structures. In eukaryotes, disulfide bonds can be found in secreted proteins or extracellular domains.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 230 of 320</div> <div> Figure 13 cystine</div> <div> Amino acids not only have physical and chemical properties, but also structural properties. These structural properties are equally important in giving rise to protein structures.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 231 of 320</div> <div> Module 007: STRUCTURAL TRAIT OF AMINO ACID-II</div> <div> Each amino acid has a unique set of properties such as charge state, polarity and hydrophobicity. Moreover, it may have unique structural traits as well which can help in protein folding. Since some amino acids are hydrophobic, they may be employed in forming a stable core in a protein. Also, chemically inactive amino acids reduce chances of destabilizing reactions in core.</div> <div> There comes a problem in burying hydrophobic amino acids in protein core Backbone is highly polar (hydrophilic) due to polar -NH and C=O in each peptide unit; these polar groups must be neutralized.</div> <div> Form regular secondary structures!</div> <div> Such as:</div> <div> • Alpha Helices</div> <div> • Beta Sheets</div> <div> Which are stabilized by H-bonds!</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 232 of 320</div> <div> Module 008: STRUCTURAL TRAIT OF AMINO ACID-III</div> <div> The size and structure of each amino acid is unique. Coupled with their chemical properties, each amino acid can uniquely contribute in the protein folding process.</div> <div> Hydrophobic core formed by packed secondary structural elements provides compact, stable core. Upon establishment of a stable protein core, unstable or reactive groups can be added.</div> <div> "Functional groups" of protein are attached to the hydrophobic core framework. Surface or a protein or its exterior must have more flexible regions (loops) and polar/charged residues.</div> <div> The very few hydrophobic "patches" on protein surface are involved in protein-protein interactions. The active regions in a protein are almost all present on the surface.</div> <div> Figure 0.44 Organization of core and surface in a protein</div> <div> Each component of the protein structure has a unique and precise role in the construction of proteins. Hydrophobic and hydrophilic components have equally useful roles.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 233 of 320</div> <div> Module 09: STRUCTURAL TRAIT OF AMINO ACID-IV</div> <div> The size and structure of each amino acid is unique. Coupled with their chemical properties, each amino acid can uniquely contribute in the protein folding process.</div> <div> Figure 0.55 Alpha Helix C = black O = red N = blue</div> <div> Alpha Helix is an example of amino acid folding. Stabilized by H-bonds between every ~ 4th residues in backbone. Reactive amino acids are exposed for external interactions.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 234 of 320</div> <div> Module 010: INTRODUCTION TO PROTEIN FOLDING</div> <div> Proteins are made by polymerization of amino acids on ribosomes and proteins properties are linked to the properties of amino acids. There are 20 amino acids in nature each has different chemical composition and that’s why each protein is different from other.</div> <div> But how does a protein actually fold? The answer is still unknown. Scientists have spent decades in trying to find a definite answer to this question, but to no avail. After polymerization of amino acids, linear chains are formed. When these chains of amino acids are put in water, the proteins fold spontaneously. The folded protein molecule should have the lowest possible energy. Anfinsen's dogma (also known as the thermodynamic hypothesis) is a postulate in molecular biology that, at least for small globular proteins, the native structure is determined only by the protein's amino acid sequence. Unique, stable and kinetically accessible minimum free energy</div> <div> Figure 0.66 Overall Energy (stability) of the Protein</div> <div> Proteins fold spontaneously in water. Proteins fold to achieve thermodynamic stability. Proteins fold to organize themselves for performing functions in cells.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 235 of 320</div> <div> Module 011: IMPORTANCE OF PROTEIN FOLDING</div> <div> Proteins are like functional machines in cell, therefore understanding the folding behavior of proteins can helps us in designing the suitable drug. If a protein is misfolded, then it can lead to a lack of function in the protein. To study anomalies in structures and to discover newer structural forms, computational algorithms are used.</div> <div> We can study the folding behavior of protein computationally First, we collect clues & evidences from experimentally reported structures. We utilize these observations to analyze unknown structures. The manner in which a newly synthesized chain of amino acids transforms itself into a perfectly folded protein depends both on the intrinsic properties of the amino-acid sequence (Dobson 2003)</div> <div> Dobson, C. M. (2003). "Protein folding and misfolding." Nature 426(6968): 884-890.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 236 of 320</div> <div> Module 012: COMPUTING PROTEIN FOLDING POSSIBILITIES</div> <div> Computing the protein folding can help us study misfolding, interaction between drugs and proteins etc. However, first, it is important to know the number of the protein folding possibilities.</div> <div> Let’s assume that each amino acid can fold into three different conformations. They are Alpha Helices, Beta Sheets and Loops. We know that proteins comprise of 100s of amino acids</div> <div> If each amino acid can take 3 different conformations, and its parent protein has 100 amino acids, then 1003 = 5 x 1047 will be the combination. If it take 1/10th of a Nano-second (10 -10), then to compute all the folding possibilities will take 1.6 x 1030 years.</div> <div> In fact, it take a protein less than a second to fold. It’s the Amazing speed of folding.</div> <div> Figure 0.77 Overall Energy (stability) of the Protein</div> <div> This is called “Levinthal’s Paradox”. We will try to understand this folding process using experimental datasets and algorithms. Molecular simulations are also helpful for it.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 237 of 320</div> <div> Module 013: PROCESSING OF PROTEIN FOLDING</div> <div> Levinthal’s Paradox- enormous time required to compute all folding possibilities. It’s impossible to consider all the possibilities computationally. So, we are trying to understand the folding process.</div> <div> The forces involved in protein folding include:</div> <div>  Electrostatic interactions</div> <div>  van der Waals interactions</div> <div>  Hydrogen bonds</div> <div>  Hydrophobic interactions</div> <div> Figure 0.88 Protein folding</div> <div> Figure 19 Anfinsen’s Experiment</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 238 of 320</div> <div> Figure 20 Anfinsen’s Experiment</div> <div> All the information required for folding a protein into its native structure is present within the protein’s amino acid sequence. The native folded form of protein is thermodynamically most stable as compared to others</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 239 of 320</div> <div> Module 014: MODELS OF PROTEIN FOLDING</div> <div> Information required for folding a protein into its native structure is present within the protein’s amino acid sequence. The native folded form of protein is thermodynamically most stable as compared to others.</div> <div> FRAME WORK MODEL</div> <div> Figure 21 Step 1: Formation of secondary structures</div> <div> Figure 0.92 Step 2: Arrangement of secondary structures</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 240 of 320</div> <div> NUCLEAR CONDENSATION MODEL</div> <div> Figure 20.10 Step 1: Formation of a Hydrophobic Core</div> <div> Figure 20.11 Step 2: Including remaining amino acids and expanding the nucleus</div> <div> Several models exist for folding a protein given its amino acid sequence. The fundamental requirement is that the folding process remain spontaneous. There is still no definitive folding hypothesis.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 241 of 320</div> <div> Module 015: PROTEIN STRUCTURE</div> <div> Proteins spontaneously fold to take 3D forms. It’s a fast yet specific process which leads to a folded protein. Several forces act together to fold the protein structure.</div> <div> Figure 25 Folding funnel</div> <div> Figure 0.126 Energies of Various Bonds & Interactions</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 242 of 320</div> <div> Figure 27 Hydro peroxide resistance protein OsmC (1vla)</div> <div> Figure 28 Cystatin – 3 (C) http://beautifulproteins.blogspot.com/</div> <div> Protein structures are very complex yet they form spontaneously. We will investigate how to develop algorithms to predict such structures.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 243 of 320</div> <div> Module 016: Primary, Secondary, Tertiary and Quaternary Structures</div> <div> Proteins are made by polymerization of amino acids on ribosomes and proteins properties are linked to the properties of amino acids. There are 20 amino acids in nature each has different chemical composition and that’s why each protein is different from other.</div> <div> Complex protein structures form spontaneously as a protein folds. A huge variety of protein structures exist. Each structure is designed to perform a specific function. Interestingly, each protein mega structure gets built out of only a few sub-structures. Combinations from the SMALL substructure set are used to construct larger protein structures.</div> <div> There are many types of structure Single Alphabet Amino acid tags can be put together linearly to represent a protein sequence. This sequence is also called the primary sequence. Primary sequence can also be referred to as 1’ structure. Sub-structures are formed as a result of 1’ structure’s folding. Folded sub structures are called secondary protein structures .Secondary structures are also referred to as 2’ structures.</div> <div> 2’ sub-structures are packed together to form super structures. These protein super structures are called tertiary structures .Tertiary structures are also referred to as 3’ structures.</div> <div> 3’ structures represent the complete monomeric protein structure.3’ structures can combine with other polypeptide units to form a quaternary structure.</div> <div> Quaternary structures are also called 4’ structures. 4’ structures are exemplified by protein complexes etc.</div> <div> Protein structures are organized into 1’, 2’, 3’ and 4’ modular conformations. We will investigate how to develop algorithms to predict these structures</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 244 of 320</div> <div> Module 017: Primary Structures of Protein</div> <div> Protein structures are organized into 1’, 2’, 3’ and 4’ modular conformations. 1’ structures are essentially the amino acid sequence of the proteins.</div> <div> Figure 29 protein folding funnel</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 245 of 320</div> <div> Figure 30 list of amino acids</div> <div> There are two methods for obtaining 1’ structure.</div> <div>  Edman Degradation</div> <div>  Tandem Mass Spectrometry</div> <div> 1’ structure databases are essentially protein sequence databases. Examples include Uniprot, Swissprot amongst several others.</div> <div> Protein sequences are the primary structures of proteins. The primary or 1’ structure of a protein determines its initial properties.1’ structure lays the foundation for 2’ structures</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 246 of 320</div> <div> Module 018: Secondary Structures of Proteins - I</div> <div> The primary or 1’ structure of a protein determines its basic properties and 1’ structure lays the foundation for 2’ structures. 2’ structures are also referred to as secondary structures.</div> <div> Figure 30.13 Organization of Secondary Structure</div> <div> Formation of 2’ structure</div> <div> C- Terminus is negatively charged .N-terminus is positively charged. C and N termini can therefore make Hydrogen Bonds. Hydrogen Bonds are the reason of 2’ structure formation.</div> <div> Figure 30.14 Forming Secondary Structure</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 247 of 320</div> <div> Figure 30.15 Types of Secondary Structures – Alpha Helix</div> <div> Protein sequences fold onto themselves and make H-Bonds to create 2’ structures. Several types of 2’ structures exist. These include Alpha Helices and Beta Sheets.</div> <div> Figure 30.16 Types of Secondary Structures – Beta Sheets</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 248 of 320</div> <div> Module 019: Secondary Structures of Proteins – II</div> <div> 2’ structures or secondary protein structures are formed as a result of H-Bond formation between N and C termini in a protein backbone. Types of 2’ structures include Alpha helices and Beta sheets.</div> <div> Figure 35 A Special Secondary Structure</div> <div> Properties of Loop</div> <div>  Loops connect helices and sheets</div> <div>  Loops vary in length and 3-D configurations</div> <div>  Loops are mostly located on the surface of proteins</div> <div>  Loops are more “acceptable" of mutations</div> <div>  Loops are flexible and can adopt multiple conformations</div> <div>  Loops tend to have charged and polar amino acids</div> <div>  Loops are frequently components of active sites</div> <div> Coils</div> <div>  Secondary structure that are not helices, sheets, or recognizable turns</div> <div>  Disordered regions, but also appear to play important functional roles</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 249 of 320</div> <div> Loops and Coils are also secondary structure which form the first structures after folding of protein’s amino acids. Loops and Coils are very important 2’ structures in that they form active sites of proteins.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 250 of 320</div> <div> Module 020: Tertiary Structures of Proteins</div> <div> 2’ structures include alpha helices, beta sheets, loops and coils. Upon combination of 2’ structures, a tertiary or 3’ structure is formed.3’ structure is next level of structure organization.</div> <div> Figure 36 Example of Tertiary Structure</div> <div> Formation of 3’ structure</div> <div>  Hydrophobic interactions between nonpolar R-groups</div> <div>  Covalent bonds in the form of Disulphide bridges</div> <div> Combinations of Alpha helices, Beta sheets, coils and loops help form 3’ structures. Covalent bonds, Hydrogen bonds and hydrophobic interactions enforce the 3’ structure.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 251 of 320</div> <div> Module 021: Quaternary Structures of Proteins</div> <div> 4’ structures or quaternary structures are formed by different peptide chains that make up the protein. Multimeric proteins which comprise of multiple peptides form 4’ structures.</div> <div> Monomeric vs. Multimeric Proteins</div> <div> Protein comprised of only a single chain (monomeric) do not have a quaternary structure. Proteins with multiple chains can form 4’ structures.</div> <div> Figure 37 Example of Quaternary Structure See how 2’ and 3’ structures come together</div> <div> 4’ structures are kept in conformation by Hydrogen Bonds, Covalent Disulphide Bonds, Hydrophobic Interactions and ionic bonds. In terms of stability 4’ > 3’ > 2’ > 1’</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 252 of 320</div> <div> Module 022: Introduction to Bond Angles in Proteins</div> <div> Protein folding results in a linear chain of amino acids getting packed into a compact 3D structure. This leads to a reduction in bond angles from an initial of 180 degrees (protein’s linear form)</div> <div> Figure 38 Linear Protein</div> <div> Figure 39 Formation of Planar Peptide Bond</div> <div> The resultant chain gets its own set of attributes and Peptide bond is planar & rigid.</div> <div> Dihedral Angles</div> <div>  Angle between two planes (i.e. 4 points)!</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 253 of 320</div> <div>  Considering the middle two points to be aligned (or overlapped), the angle between the 1st, overlapped and the 4th points forms a dihedral angle.</div> <div> Figure 40 Protein after Folding: Phi and Psi Angles</div> <div> Figure 0.171 Protein after Folding: Phi and Psi Angles</div> <div> Φ (phi, involving C'-N-Cα-C') ψ (psi, involving N-Cα-C'-N)</div> <div> Proteins fold into 3D structures. Phi and psi angles are taken up as a result of folding. These angles can be measured towards understanding the protein structure.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 254 of 320</div> <div> Module 023: Ramachandran Plot</div> <div> Phi and Psi angles can be measured with in the folded structures like:</div> <div>  φ - phi</div> <div>  Involves C'-N-Cα-C‘</div> <div>  ψ – psi</div> <div>  Involves N-Cα-C'-N</div> <div> Figure 42 Phi and Psi Angles</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 255 of 320</div> <div> Figure 43 Allowable Phi and Psi Angles</div> <div> Data as in (Lovell et al. 2003) showing about 100,000 data points for several amino-acids</div> <div> A limited range of Phi and psi angles are taken up as a result of folding. This range of angles constitutes the allowable range of torsion or rotation angles that are taken up by the protein.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 256 of 320</div> <div> Module 024: Structure Visualization - I</div> <div> We know that protein backbone takes up specific rotation angles after folding. A protein consists of multiple amino acids. Each amino acid has a C-terminus and an N-Terminus.</div> <div> Figure 44 Protein Backbone and C atoms</div> <div> Figure 45 Omitting Planar bonds and Tracing C-Alpha atoms in backbone</div> <div> Figure 46 C-Alpha Backbone visualization</div> <div> C-Alpha atoms are traced to recreate a 3D protein structure. The choice is made while keeping planar nature of the peptide bond in view. Later we will see how to insert side chains into the visual models as well.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 257 of 320</div> <div> Module 025: Structure Visualization – II</div> <div> C-Alphas can be used to construct the backbone of a protein towards its visualization. We also need a representation of measurements for assigning the atomic distances. The ångström is used to express the size of atoms, molecules and extremely small biological structures, the lengths of chemical bonds, the arrangement of atoms in crystals.</div> <div> 1 angstrom is a unit of length equal to 10−10 m (one ten-billionth of a meter) or 0.1 nm</div> <div> Atoms of phosphorus, sulfur, and chlorine are ~1 Å in covalent radius, while a hydrogen atom is 0.25 Å</div> <div> Figure 47 Ansedel Anders Ångström (1814–1874)</div> <div> C-Alpha atoms are traced to recreate a 3D protein structure. Each C-Alpha atom is at a distance which can be represented in the unit “Angstrom”.1 A resolution is better than 10 A.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 258 of 320</div> <div> Module 026: Experimental Determination of Protein Structure</div> <div> C-Alpha atoms are traced to recreate a 3D protein structure. Distances between C-Alphas are measured in the unit “Angstrom”.</div> <div> X-Ray Crystallography</div> <div>  Crystallography data gives relative positions of atomic coordinates</div> <div>  The data is obtained from diffractions by the atoms in a protein structure</div> <div>  The coordinates of each atom in x,y and z axis are output</div> <div> Figure 48 x-ray crystallography</div> <div> Crystallized proteins are used to determine protein structures. As X-rays diffract from the atoms in a protein, the atomic distances are noted. These distance in 3D are measured in Angstroms.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 259 of 320</div> <div> Module 027: Protein Databank</div> <div> Position of C-Alpha atoms are used to construct 3D protein structure. X-Ray diffraction data helps measure the atomic positions. X, Y and Z positions of several proteins are available online.</div> <div> Figure 49 PDB File Format</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 260 of 320</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 261 of 320</div> <div> PDB contains protein structure information. It has the coordinates of C-Alphas for over 50,000 proteins. Protein structures can be visualized using this information.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 262 of 320</div> <div> Module 028: Visualization Technique</div> <div> Proteins fold into 3D structures. Phi and psi angles are assumed as a result of folding. These angles can be measured and viewed towards understanding the protein structure. To view a protein, we need to evaluate the physical location of its atoms. Proteins have Carbon and Nitrogen in their backbone.</div> <div> CA atomic coordinates</div> <div>  To trace the backbone of a protein, CA atoms trace can be used</div> <div>  Note that CA atoms have the side chains attached to them</div> <div>  A coordinates can be found in the PDB file</div> <div> Protein structures can be visualized by tracing the CA atoms. Coordinates of CA atoms can be obtained from the PDB. Next, we need a tool to plot these coordinates.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 263 of 320</div> <div> Module 029: Online Resources for Protein Visualization</div> <div> Protein structures can be visualized by tracing the CA atoms. CA Coordinates can be taken from PDB.</div> <div> Online Tools</div> <div>  Rasmol and CHIME are basic tools for visualizing proteins</div> <div>  Swiss PDB Viewer offers several features such as protein surface view, alignment of several proteins & modelling secondary structures</div> <div>  PyMOL is a python-script based tool for visualizing the protein structure</div> <div>  Cn3D is another tool which helps us visualize protein structures</div> <div>  It also provides for annotating protein structures</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 264 of 320</div> <div> Protein structures are visualized using several online tools. These tools include Rasmol, CHIME, Swiss PDB Viewer and Cn3D.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 265 of 320</div> <div> Module 030: Types of Protein Visualizations</div> <div> To visualize proteins, we use CA coordinates or positions. We can use several online tools to view the resulting model.</div> <div> CPK: Corey-Paulin-Koltun Diagrams. In CPK diagrams, each atom is represented by a solid sphere. Spheres are equal to atomic van der Waal radius (the volume of the atom).</div> <div> Figure 50 sphere and surface diagrams of protein</div> <div> http://www.danforthcenter.org/smith/MolView/Over/overview.html</div> <div> Ribbon Diagrams</div> <div> Ribbon diagrams are an easy and frequently used technique for representing protein structures. Structure is represented by the secondary structures (fold) using simple cartoon figures.</div> <div> Figure 51 ribbon diagrams</div> <div> http://www.danforthcenter.org/smith/MolView/Over/overview.html</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 266 of 320</div> <div> Balls & Stick (BS) Models</div> <div> BS model is another popular protein structure representation strategy. BS Models have atoms as colored balls and intermediate bonds as sticks.</div> <div> Figure 52 Balls and sticks model</div> <div> http://www.danforthcenter.org/smith/MolView/Over/overview.html</div> <div> Figure 53 Colored Sticks Models</div> <div> http://www.danforthcenter.org/smith/MolView/Over/overview.html</div> <div> Protein Structure Visualization can be performed using several atomic representations. These include CPK, Ribbon and Balls & Stick Diagrams.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 267 of 320</div> <div> Module 031: Introduction to Energy of Protein Structures</div> <div> Proteins come together as a result of peptide bond formation between various amino acids. The resulting polymer then goes through the step of folding which leads to the formation of a 3D structure.</div> <div> https://folding.stanford.edu/home/the-science/</div> <div> Role of Amino Acids</div> <div> We know that amino acids can be polar, charged and hydrophobic. Role of polar and charged amino acids in folding. Role of hydrophobic amino acids in folding.</div> <div> Overall Goal of Folding</div> <div> Anfinsen’s thermodynamic hypothesis: Proteins fold for a unique, stable and minimum free kinetic energy structure. What other factors may come into play for satisfying Anfinsen hypothesis.</div> <div> Minimizing Energy</div> <div> We know that if bonds can be formed between two atoms, then energy is released. This leads to a situation where there is lesser free energy accessible to each atom for further interactions. So, proteins maximize bonds that can be made between the side chains on each of their constituent amino acids</div> <div> Such atomic interactions include:</div> <div>  Disulphide bonds between Cysteine residues</div> <div>  Hydrogen Bonds</div> <div>  Van der Waals Forces</div> <div>  Electrostatic Interactions between polar/charged amino acids</div> <div> The greater the number of these bonds, the more stable a protein becomes. Hence, the basic idea of thermodynamic stability is to maximize bonding in order to minimize the free energy</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 268 of 320</div> <div> Module 032: Calculating Energy of a Protein Structure</div> <div> As we know the greater the number of bonds between the amino acids, the more stable a protein becomes.</div> <div> Figure 54 Energies of Interactions www.ucdavis.edu</div> <div> Comparison of bond energies</div> <div>  Hydrophobic interactions ></div> <div>  Electrostatic interactions></div> <div>  Hydrogen bond > van der Waals</div> <div> Calculating overall energy of a protein structure</div> <div> Given the number of atomic interactions in a protein, you can simply sum the energy in the protein molecule.</div> <div> Energies of protein structures can be computed by first enumerating the types of interactions between each atom. Then, accumulating the energy of each interaction towards calculating an overall energy of a protein.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 269 of 320</div> <div> Module 033: Structure Determination for Energy Calculations</div> <div> The greater the interactions between the amino acids, the more stable a protein becomes. We can calculate energy of a folded protein based on the number and types of atomic interactions.</div> <div> How to find the number of interactions</div> <div>  To determine the number of each type of interaction within a protein, we need to find its inter-atomic distances.</div> <div>  Based on specific atomic distances, we can guess the type of atomic interaction.</div> <div>  By looking up at the bond/energy table, we can compute the overall energy.</div> <div> Techniques for structure determination</div> <div>  X-Ray Crystallography</div> <div>  Nuclear Magnetic Resonance (NMR) Spectroscopy</div> <div> We need to know the structure of the protein to calculate atomic distances. Atomic distances tell us about atomic interactions with neighboring atoms. To determine the structure, we use X-Ray or NMR.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 270 of 320</div> <div> Module 034: Review of Experimental Structure Determination</div> <div> The greater the interactions between the amino acids, the more stable a protein becomes. We can calculate energy of a folded protein based on the number and types of atomic interactions.</div> <div> The structure also dictates which functions a protein can perform via the positioning of hydrophilic & polar amino acids. For determining stability, structure & function, we need to find the amino acid interactions. Several experimental methods exist for structure determination.</div> <div>  X-Ray Crystallography</div> <div>  Nuclear Magnetic Resonance (NMR) Spectroscopy</div> <div> Figure 55 to measure a bond/interaction, we must first see atoms</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 271 of 320</div> <div> Figure 56 Principle of X-Ray Crystallography</div> <div> Figure 57 from Diffraction Patterns to Atomic Positions</div> <div> Upon establishing the atomic positions and distances, we can then check for possible interaction between the different atoms. Atomic distances can help us classify interaction types e.g. hydrogen bonds, electrostatic & polar.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 272 of 320</div> <div> Module 035: Alpha Helices - I</div> <div> Atomic distances can tell us about their existential interactions. Different types of interactions may occur between atoms. E.g. Hydrogen Bonds, Polar etc. If two atoms are participating in a covalent bond, their distance is ~0.96A. In case of hydrogen bond formation between atoms, the inter-atomic distance is ~1.97A. X-Ray data should have a minimum of 1.97A resolution.</div> <div> Figure 58 Hydrogen Bonds to Fold an Amino Acid Chain</div> <div> X-Ray Crystallography data shows that Hydrogen atoms of N-Term may come together with Oxygen atoms of C-term amino acid at 4th neighboring position. Their atomic distance is ~1.9A and hence are considered to be in a hydrogen bonds.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 273 of 320</div> <div> Module 036: Alpha Helices – II</div> <div> X-Ray Crystallography of protein shows that Hydrogen atoms of N-Term come together with Oxygen atoms of C-term amino acid at 4th neighboring position to make Hydrogen bonds.</div> <div> Figure 59 Forming Alpha Helix</div> <div> Every Oxygen bound to 4th neighboring Amino Group’ Hydrogen.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 274 of 320</div> <div> Figure 60 Carbons (Black) & Nitrogen’s (Blue): 1-5, 2-6, 3-7…</div> <div> Figure 61 Preference of Amino Acids for making Alpha Helices</div> <div> Helix Formers</div> <div> From 20 amino acids, anyone can be present in the backbone. Is there a variable preference in amino acids to form helix? Yes, “Helix Formers” are generally hydrophobic amino acids (M, A, L…). Alpha Helices are formed by hydrogen bonding (O-H) between Ci and Ni+4 atoms in the protein backbone.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 275 of 320</div> <div> Module 037: Beta Sheets - I</div> <div> Alpha Helices are formed by hydrogen bonding (O-H) between Ci and Ni+4 atoms in the protein backbone. Beta Sheets are another common secondary structure. They are constituted by several Beta Strands which come together. 5 to 10 resides are needed to make a Beta Strand, typically.</div> <div> Hydrogen Bonds to make in Beta Strands</div> <div> The Beta Sheet is made up of several Beta Strands</div> <div> C-Alpha atoms and the CO and NH groups are shown in blue, yellow, and green, respectively.</div> <div> This is called a parallel beta sheet.</div> <div> This is called an anti-parallel beta sheet.</div> <div> Beta Sheets are another secondary structure that can be formed as a result of hydrogen bonding between the protein back bones. Some amino acids have a preference for making Beta Sheets.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 276 of 320</div> <div> Module 038: Beta Sheets - II</div> <div> Beta strands can make hydrogen bonds with each other and organize as beta sheets.</div> <div> Beta Sheets have different Properties:</div> <div>  Beta Strand</div> <div>  Beta Sheet</div> <div>  Beta Barrel</div> <div>  Beta Sandwiches</div> <div> Beta Barrels</div> <div> Beta Barrel is made of a single beta sheet that twists and coils upon itself. The first strand in the beta sheet makes a hydrogen bonds with the last strand. A beta barrel is a large beta-sheet that twists and coils to form a closed structure in which the first strand is hydrogen bonded to the last. Beta-strands in beta-barrels are typically arranged in an antiparallel fashion. https://en.wikipedia.org/wiki/Beta_barrel</div> <div> Figure 62 beta barrel</div> <div> Beta Sandwiches</div> <div> Beta Sandwiches are made of two beta sheets which are usually twisted and packed so their strands are aligned.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 277 of 320</div> <div> Figure 63 Illustration of the β-sandwich from Tenascin C (PDB entry: 1TEN).</div> <div> Figure 64 Preference of Amino Acids for making Beta strands</div> <div> Beta Sheets are formed by H bonds between of 5–10 consecutive amino acids in one portion of the backbone with another 5–10 farther down the backbone. Beta strands may be adjacent (with a loop in between) or far with other structures in between.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 278 of 320</div> <div> Module 039: LOOPS-I</div> <div> Alpha Helices and Beta Sheets are secondary structures formed as a result of hydrogen bonding in between protein backbone atoms.</div> <div> Protein Backbone and Secondary Structures</div> <div> Loops are formed by amino acids present in the middle of the Alpha Helices and Beta Sheets in a protein backbone.</div> <div> Figure 65 Joining Alpha Helices and Beta Sheets in a Protein Backbone</div> <div> Variability in length and conformation allows loops to join Alpha Helices and Beta Sheets in a variety of ways. Loops are variable in length and 3-D conformations.</div> <div> Characteristics</div> <div>  Loops are mostly located on the surface of protein structure</div> <div>  Mutate in sequence at a much faster rate than Alpha Helices and Beta Sheets</div> <div>  Loops are flexible and can adopt multiple conformations</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 279 of 320</div> <div> Loops dictate the overall structure of protein as they couple Alpha helices and beta sheets</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 280 of 320</div> <div> Module 040: LOOPS-II</div> <div> Loops dictate the overall structure of protein as they couple Alpha helices and Beta sheets. Loops are flexible and have variable lengths so as to successfully bridge between secondary structures.</div> <div> Figure 66 Loops in 3D Conformation</div> <div> Loop Properties</div> <div>  Loops are mostly comprised of charged and polar amino acids</div> <div>  Loops frequently participate as components of active sites</div> <div> Figure 67 Preference of Amino Acids for making Loops</div> <div> Types of Loops</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 281 of 320</div> <div>  Hairpin loops are two amino acids long and join anti-parallel Beta strands</div> <div>  Other Loops may be 3 to 4 amino acids long</div> <div>  Loops fall into various families</div> <div> Loops are the third type of secondary structure after Alpha helices and Beta sheets. Loops are unique in that they are flexible and variable length. Loops constitute active sites.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 282 of 320</div> <div> Module 041: COILS</div> <div> Alpha helices and beta sheets are the regular secondary structures. Loops are flexible secondary structures &connect alpha helices and beta sheets. Coils are another secondary structure. Coils are unstructured and unlike loops. Essentially, a secondary structure which is not a helix, sheet or loop is a coil.</div> <div> Functional Aspects of Coils</div> <div>  Coils are apparently disordered regions</div> <div>  They are oriented randomly while being bonded to adjacent amino acids</div> <div>  However, coils also appear to play important functional roles</div> <div> Figure 68 Coils in Myoglobin</div> <div> Coils are those secondary structure formed by the protein backbone which are neither helices, sheets nor loops. In fact, coils do not have a consistent classifiable structure. Hence, coils are random structure and random length.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 283 of 320</div> <div> Module 042: Structure Classification - I</div> <div> Proteins have primary, secondary, tertiary & quaternary structures. Each level of protein structure organization is known to impart specific characteristics to the protein.</div> <div> Review of the 4 structure levels</div> <div>  Primary Structures</div> <div>  Secondary Structures</div> <div>  Tertiary Structures</div> <div>  Quaternary Structures</div> <div> Structural artifacts tend to be more conserved as compared to their sequences. Therefore, it may be useful to look at the secondary/tertiary structures for conservation study.</div> <div> Classification</div> <div>  The evolution of protein structures and their hierarchy is not systematized</div> <div>  Hence, we need to classify the function of protein by examining their secondary and tertiary structures</div> <div> Motifs (Non-functional Combinations of 2’ structures)</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 284 of 320</div> <div> Figure 69 Domain (Functionally Complete)</div> <div> Domains are semi-independent functional structures in a protein. Have a stable structure. Over ~40 residues. Protein may contain multiple domains.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 285 of 320</div> <div> Module 043: Structure Classification - II</div> <div> Domains are semi-independent functional structures in a protein. Protein may contain multiple domains. Hence, we can try to classify proteins by their domains. Locally Compact – Domains interact (H-bonds) more internally than externally. Domains have a hydrophobic core. Domains are contiguous (min. chain breaks).</div> <div> Domains have a minimal contact with rest of the peptide. Solvent area in contact with each domain should not vary significantly upon separating two separate domains.</div> <div> Types of Domains</div> <div>  Alpha Domains</div> <div>  Beta Domains</div> <div>  Alpha/Beta Domains</div> <div>  Alpha + Beta Domains</div> <div>  Alpha & Beta Multi-Domains</div> <div>  Membrane & cell-surface proteins</div> <div> So, by looking at proteins, we can list the domains present in each protein. Once domains in each protein are listed, we can classify whole proteins into various types and classes.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 286 of 320</div> <div> Module 044: Examples of Protein Domains</div> <div> There are many domains for protein structure prediction.</div> <div>  Alpha Domains</div> <div>  Beta Domains</div> <div>  Alpha/Beta Domains</div> <div>  Alpha + Beta Domains</div> <div>  Alpha & Beta Multi-Domains</div> <div>  Membrane & cell-surface proteins</div> <div> Figure 70 Alpha Domain: Hemoglobin (1bab)</div> <div> Figure 71 Immunoglobulin (8fab)</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 287 of 320</div> <div> Figure 72 Alpha / Beta: Triosephosphate isomerase (1hti)</div> <div> Figure 73 Alpha + Beta: Lysozyme (1jsf)</div> <div> Various types of domain architectures exist in proteins. Such architectures can be classified into general structural classes. Databases can be made from classes.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 288 of 320</div> <div> Module 045: CATH Classification</div> <div> Domains can be classified into structural classes. Classes can be further classified into Architecture and Topologies. Let’s see how it is done in CATH.</div> <div> Figure 74 Structural Classes</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 289 of 320</div> <div> Class</div> <div>  Similar secondary structure content</div> <div>  All α, all β, alternating α/β etc.</div> <div> Architecture</div> <div>  Also called FOLD</div> <div>  Major structural similarity</div> <div>  SSE’s in similar arrangement</div> <div> Topology</div> <div>  Super Family</div> <div>  Probable common ancestry</div> <div>  Family membership</div> <div> Homology</div> <div>  Same Family</div> <div>  Clear evolutionary relationship</div> <div>  Pairwise sequence similarity > 30%</div> <div> CATH classifies proteins by their structural similarity. It also considers the internal organization of the structural components in proteins.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 290 of 320</div> <div> Module 046: Classification Databases</div> <div> Proteins are classified into various structural classes. CATH is one such system in which proteins are organized into classes, architecture, topology and homology.</div> <div> http://scop.mrc-lmb.cam.ac.uk/scop/</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 291 of 320</div> <div> Figure 75 SCOP Classification Statistics</div> <div> http://scop.mrc-lmb.cam.ac.uk/scop/count.html</div> <div> FSSP - Family of Structurally Similar Proteins, based on the DALI algorithm. Pclass - Protein Classification, based on the LOCK and 3Dsearch algorithms.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 292 of 320</div> <div> Module 047: Algorithms for Structure Classification</div> <div> Several algorithms exist for classifying protein structures.</div> <div> Intra-Molecular Distance Algorithms.</div> <div>  Proteins are considered as rigid bodies.</div> <div>  They are placed in a 3D Cartesian coordinate system.</div> <div>  Structural alignment in 3D.</div> <div>  E.g. VAST, LOCK</div> <div> Inter-Molecular Distance Algorithms</div> <div>  Proteins are considered as rigid bodies.</div> <div>  They are placed in 2D.</div> <div>  Structural alignment using internal distances and angles.</div> <div> The basic idea is to capture internal geometry of protein structures. E.g. DALI, and SSAP.</div> <div> Such algorithms are also very useful to compare whole protein structures. They can help determine evolutionary relationship. Also, functional similarity can be estimated.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 293 of 320</div> <div> Module 048: Protein Structure Comparison</div> <div> Proteins are assembled into primary (1’), secondary (2’), tertiary (3’) and quaternary (4’) structures. Protein sequence is less conserved than its structure. Protein structure determines function. Since protein structure dictates function, comparing two structures can help us evaluate if the proteins do the same or similar function.</div> <div> Comparing Whole Protein Structures</div> <div> Proteins contain multiple structural subunits e.g. secondary structures, motifs and domains. Structures of all such subunits are to be considered as one and compared. We know that domains are functionally independent components of the protein structure. Proteins may have multiple domains. So for two different proteins, sharing the same domain, we may want to compare only a portion of the overall structure i.e. a domain. For comparing the complete or partial protein structures, the position of Alpha Carbon atoms can be used. The (x, y, z) positions of Alpha Carbon atoms can be obtained from the PDB.</div> <div> Figure 74 C-Alpha atoms in backbone</div> <div> Figure 76 Tracing and Visualizing C-Alpha Backbone</div> <div> http://www.danforthcenter.org/smith/MolView/Over/overview.html</div> <div> PDB coordinates of Alpha Carbons in the protein back bone can be used for comparison. In this way, whole protein structure or domains etc. can be compared.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 294 of 320</div> <div> Module 049: Strategies for Structure Comparison - I</div> <div> PDB coordinates of Alpha Carbons in the protein back bone can be used for comparison. Thus, two whole protein structures or domains within each structure can be compared.</div> <div> Figure 77 Tracing and Visualizing C-Alpha Backbone</div> <div> http://www.danforthcenter.org/smith/MolView/Over/overview.html</div> <div> Strategy # 1 – Whole Protein Structure Comparison by Intermolecular distances</div> <div>  Two protein sequences are pair-wise aligned with each other</div> <div>  Corresponding Alpha Carbons are identified</div> <div>  Coordinates of corresponding Alpha Carbons are retrieved from PDB</div> <div>  Their individual differences calculated</div> <div>  Root Mean Square Distance is computed to assess the similarity</div> <div> Whole protein structures can be compared by calculating the root mean squared difference (RMSD) between their Alpha Carbons positions. The lower the RMSD, the similar are the proteins.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 295 of 320</div> <div> Module 050: Strategies for Structure Comparison - II</div> <div> Full protein structures can be compared and ranked by the overall differences in positions</div> <div> between their Alpha Carbons. But proteins are 3D and in various conformations.</div> <div> Full Protein Comparison</div> <div> (Translation -> Rotation -> RMSD)</div> <div> Calculating RMSD – An Example</div> <div> </div> <div> </div> <div> N</div> <div> i</div> <div> i i d a b</div> <div> N</div> <div> RMS A B</div> <div> 1</div> <div> 2 ( , )</div> <div> 1</div> <div> ( , )</div> <div> Domain or Motif Comparison</div> <div> (Region Selection ->Translation -> Rotation -> RMSD)</div> <div> Motifs, Domains and Full Proteins can be compared by using the rigid body super-positioning.</div> <div> Depending on the RMSD, proteins, their motifs and domains can be selectively compared.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 296 of 320</div> <div> Module 051: Online Resources for Structure Comparison</div> <div> Multiple types of comparison can be performed between Proteins, Motifs, and Domains by rigid body super-positioning. RMSD tells us about the quality of the matches.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 297 of 320</div> <div> Protein structures can be compared in multiple ways. Till now, we can compare proteins by their motifs, domains and full structures. There are several advanced techniques for this as well.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 298 of 320</div> <div> Module 052: Protein Structure Prediction</div> <div> Complex protein structures enable proteins to perform complex functions. We know over a million protein sequences but only about 100,000 protein structures. Estimating exact protein structures is very difficult. It’s difficult to crystallize proteins. Even if we manage to get protein’s X-Ray, to reconstruct the structure is extremely complex.</div> <div> Since we know so many sequences, they can be used for predicting protein structures. This indeed is possible and helpful.</div> <div> The Basic Idea</div> <div>  Amino acids determine the protein structure</div> <div>  We have a large protein sequence dataset (uniprot)</div> <div> Hence, we can fold protein sequences and predict their structures</div> <div> Why predict and why not exact solutions?</div> <div> A deterministic solution of protein folding is a major unsolved problem in molecular biology. Proteins fold spontaneously or with the help of enzymes or chaperones. To computationally predict protein structures, we need to copy or mimic the natural folding.</div> <div> To fold we must learn the steps</div> <div> Step 1: "Collapse"- leading to burial of hydrophobic AA’s</div> <div> Step 2: Fluid globule - helices & sheets form, but are unorganized</div> <div> Step 3: Compaction, and rearrangement of 2‘structures</div> <div> Protein structure prediction involves learning how the amino acids in primary sequence fold. Using this information, upon getting a protein sequence, we can try to predict how it folds</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 299 of 320</div> <div> Module 053: Predicting Secondary Structures</div> <div> By looking at the structures in PDB, we know that Alanine mostly found in Alpha Helices. So if we have several Alanines in the sequence, then we can anticipate that a helix may be formed by them. What if we survey the entire PDB and check the presence of each amino in each type of secondary structure. If we know which amino acid is found in which specific secondary structure, then we can use it for prediction.</div> <div> Figure 78 Chou & Fasman (1974 & 1978)</div> <div> Several algorithms have been designed to predict 2’ given an amino acid sequence. The first such algorithm was the Chou-Fasman Algorithm. We will see it in the upcoming modules.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 300 of 320</div> <div> Module 054: Introduction to Chou Fasman Algorithm</div> <div> 3D Structure of proteins is determined by their Amino Acid sequence. Note that we only know 100,000 3D protein structures, but 10 times more sequences. For those proteins whose structure is already known, can we evaluate their amino acid sequence?</div> <div> Figure 79 Propensity Table</div> <div> Predicting the 2’ structures</div> <div> Now, let’s consider that if we are given an amino acid sequence, we can simply look up the propensity table and assign the tentative secondary structure.</div> <div> Given an amino acid sequence, look up the propensity table for each amino acid’s propensity for various 2’ structures. Product of these propensity values will give you the overall propensity for formation of each 2’ structure.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 301 of 320</div> <div> Module 055: 2’ Structures in Chou Fasman Algorithm</div> <div> For a primary sequence, and a tentative 2’ structure, propensity table can help us compute the overall propensity. Product of propensity values is computed for overall propensity for each 2’ structure. An important point to note here is that 2’ structures are formed due to hydrogen bonding between amino acids.</div> <div> So, we need to consider the neighboring amino acids as well.</div> <div> You only need to compute propensities for a small number 2’ structures. The highest net propensity will be the most probably secondary structure that will be formed.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 302 of 320</div> <div> Module 056: Chou Fasman Algorithm - I</div> <div> Only a small number of combinations of secondary structures are possible due to their individual properties. Such as 4 amino acids are needed to start an Alpha Helix and 5 amino acids for Beta Sheet. Note that besides the alpha helix and beta sheets, LOOPS are another secondary structure. Loops are small ~ 3-4 amino acids.</div> <div> 1. Scan through the sequence : E M A V I Y P G</div> <div> 2. Identify sequence regions where:</div> <div> • 4 out of 6 contiguous residues give a P(α ) > 1.0</div> <div> • That region is declared as alpha-helix</div> <div> • Extend helix to both sides until 4 out of 6 contiguous residues give a P(α ) < 1.0</div> <div> That is declared end of the helix. For Alpha Helices, 4 contiguous amino acids are required. Their Alpha-Helix propensity should be more than 1.0. Once this propensity falls below 1.0, Alpha-Helix stops.</div> <div> F</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 303 of 320</div> <div> Module 057: Chou Fasman Algorithm - II</div> <div> Alpha Helices are formed from 4 contiguous amino acids having an Alpha-Helix propensity over 1.0. The Alpha-Helix stops if this propensity falls below 1.0. Once Alpha Helices are constructed, and concluded, the remaining amino acids can be evaluated for Beta sheets and turns etc. Let’s see how Beta sheets are evaluated using Chou Fasman Algorithm.</div> <div> 1. Compute P(β) for contiguous regions of 5 Amino Acids</div> <div> 2. From these regions, identify regions where:</div> <div> 3. 5 contiguous residues have P(α ) > P(β)</div> <div> That region is finalized as alpha-helix.Repeat this step for the full amino acid sequence to finalize all possible alpha helical regions in the sequence.</div> <div> Alpha Helices can be finalized if their propensity is higher than the propensity for Beta Sheets in regions of 5 amino acids. For those regions where that is not the case, further evaluation is required.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 304 of 320</div> <div> Module 058: Chou Fasman Algorithm - III</div> <div> Alpha Helices are formed from 4 contiguous amino acids having an Alpha-Helix propensity over 1.0. The Alpha-Helix stops if this propensity falls below 1.0. Alpha Helices were finalized if their propensity was higher than the propensity for Beta Sheets in regions of 5 amino acids.</div> <div> We can evaluate such regions for Beta Sheets. Let us see step by stop how to find a beta sheet and how to differentiate them from alpha helices.</div> <div> Scan the sequence to identify regions where:</div> <div>  3 out of 5 amino acids have P(β) > 1.0</div> <div>  That region is declared as beta sheet</div> <div>  Extend beta sheet to both sides until 4 contiguous residues average P(β) < 1.0</div> <div>  That is declared end of the beta sheet</div> <div>  Those regions are finalized as beta-sheets which have average P(β) > 1.05 and the average P(β) > P(α) for that region.</div> <div> Regions where overlapping alpha-helices and beta-sheets occur are declared helices if</div> <div>  the average P(a-helix) > P(b-sheet) for that region</div> <div> Else, a beta sheet is declared if</div> <div>  average P(b-sheet) > P(a-helix) for that region</div> <div> Using the strategy of higher propensity, alpha helices and beta sheets can be completely resolved. Assignments for each beta sheet and alpha helix can be finalized.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 305 of 320</div> <div> Module 059: Chou Fasman Algorithm - IV</div> <div> After computing the propensity of alpha helices and beta sheets, we need to settle for loops. Let’s see how we can find out the loops using Chou Fasman Algorithm.</div> <div> For any jth residue in sequence, we calculate f (Total) = f(j) f(j+1) f(j+2) f(j+3) (tetrapeptide)</div> <div> If</div> <div>  f(Total) > 0.000075</div> <div>  the average value for P(turn) > 1.00 in the tetra peptide</div> <div>  the averages for the tetra peptide are such P(a-helix) < P(turn) > P(b-sheet)</div> <div> http://fasta.bioch.virginia.edu/fasta_www/chofas.htm</div> <div> Chou Fasman Algorithm helps predict Alpha Helices, Beta Sheets and Turns. The algorithm is based on statistical occurrence of Amino Acids in known structures.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 306 of 320</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 307 of 320</div> <div> Module 060: Chou Fasman Algorithm – Flowchart I</div> <div> Chou Fasman Algorithm helps predict secondary structures such as Alpha Helices, Beta Sheets and Turns. Step by step flowchart of the entire algorithm.</div> <div> Beta sheets can be predicted from primary amino acid sequences. Next, we will see the flowchart of Alpha Helices and Beta Turns.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 308 of 320</div> <div> Module 061: Chou Fasman Algorithm – Flowchart II</div> <div> Chou Fasman Algorithm helps predict secondary structures such as Alpha Helices, Beta Sheets and Turns. Step by step flowchart of the entire algorithm.</div> <div> Figure 80 Beta sheet flowchart</div> <div> Figure 81 Alpha helices flowchart</div> <div> Now we have reviewed flowcharts for Alpha Helices and Beta Sheets. Next up is the flow chart for Beta Turns.</div> <div> Module 062: Chou Fasman Algorithm – Flowchart III</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 309 of 320</div> <div> Chou Fasman Algorithm helps predict secondary structures such as Alpha Helices, Beta Sheets and Turns. Step by step flowchart of the entire algorithm.</div> <div> Figure 82 Beta sheet flowchart</div> <div> Figure 83 Alpha helices flowchart</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 310 of 320</div> <div> Figure 84 Beta turn</div> <div> Alpha helices, beta sheets and turns can be predicted using Chou Fasman Algorithm. This algorithm is based on statistical analysis of amino acid occurrences in proteins.</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 311 of 320</div> <div> Module 063: Chou Fasman Algorithm – Improvements</div> <div> Alpha helices, beta sheets and turns can be predicted using Chou-Fasman Algorithm. The algorithm is based on statistical analysis of amino acid occurrences in proteins.</div> <div> Secondary structure propensity values of alpha helix, beta sheet and turns should be recalculated with the latest protein data sets.</div> <div> IMPROVEMENTS</div> <div> Special consideration for:</div> <div>  Nucleation regions</div> <div>  Membrane proteins</div> <div>  Hydrophobic domains</div> <div>  Consider variable coil and loop sizes besides the from tetra peptide turns</div> <div>  Consider local protein folding environments</div> <div>  Solvent accessibility of residues</div> <div>  Protein structural class</div> <div>  Protein’s organism</div> <div> Chou Fasman can be improved to better predict secondary structures by incorporating biochemical factors and updated statistics!</div> <div> Handout Chapter 06: Protein Structures</div> <div> Page 312 of 320</div> <div> Module 064: Summary of Visualization, Classification & Prediction</div> <div> Structure Classification</div> <div>  relationship between protein structure and function</div> <div>  There is need to classify proteins</div> <div>  Hierarchy of classification</div> <div> Structure visualization, classification and prediction equip us to perform functional evaluation of proteins. This is important for understanding disease and designing drugs for treating them.</div> </div> </div>

Bioinformatics I (BIF401) book HANDOUTS / POWER POINT SLIDES in pdf

Bioinformatics I (BIF401) book HANDOUTS / POWER POINT SLIDES in pdf Handouts | Lectures | Contents | Books Chapter 1- Introduction to ...

Read more »