{"id":40106,"date":"2021-09-30T10:30:09","date_gmt":"2021-09-30T10:30:09","guid":{"rendered":"https:\/\/biomedpharmajournal.org\/?p=40106"},"modified":"2021-10-11T10:21:29","modified_gmt":"2021-10-11T10:21:29","slug":"genetic-algorithm-approach-to-find-the-estimated-value-of-hmm-parameters-for-ns5-methyltransferase-protein","status":"publish","type":"post","link":"https:\/\/biomedpharmajournal.org\/staging\/vol14no3\/genetic-algorithm-approach-to-find-the-estimated-value-of-hmm-parameters-for-ns5-methyltransferase-protein\/","title":{"rendered":"Genetic Algorithm Approach to Find the Estimated Value of HMM parameters for NS5 Methyltransferase Protein"},"content":{"rendered":"

Introduction<\/strong><\/p>\n

Dengue virus (DV), the causative agent of dengue, resides in the family Flaviviridae and is transmitted to humans by biting Aedesaegypti mosquitoes. Four serotypes (Dengue Virus serotype 1, Dengue virus serotype 2, Dengue virus serotype 3, and\u00a0Dengue virus serotype 4) are recognized (El Sahili, Lescar2017).The range of dengue disease spans from a flu-like disease termed dengue fever to Dengue hemorrhagic fever. In chronic cases, it causes dengue shock syndrome and sometimes terminates in death. The most prevalent clinical symptoms of acute dengue disease are hemorrhagic\u00a0diathesis, liver involvement, and plasma leakage.The DV genome is prepared into a single open reading frame (ORF)of single-stranded (positive -sense) RNA of 900 kDa and flanked at 5’end by type I cap and at 3’end by untranslated regions and encodes a\u00a0precursor polyprotein. Post-translational modification of precursor protein gives rise to three structural (C, prM, and E) proteins and seven non-structural(NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5)proteins(Anasiret al.,<\/em>2020).<\/p>\n

The NS5, non-structural methyltransferase, is the largest and conserved region of the genome among all serotypes. It is a bifunctional enzyme with a Methyltransferase domain (MTase; residues 1\u2013296) at its N-terminal end and anRNA-dependent RNA\u00a0polymerase (RdRp; residues 320\u2013900) at its C-terminal end (Vannice2016).This protein has three active sites, namely S-adenosylmethionine(SAM)S-adenosyl-homocysteine (SAH), guanosine triphosphate (GTP), and RNA-binding site (Herreroet al<\/em>., 2013).<\/p>\n

In this work, we considerthe NS5 Methyltransferase protein of the Dengue Virus to find out the probability of observation sequence. Firstly we convert the protein sequence into nucleotide sequence. After that, we implement the G.A. (selection, crossover,\u00a0and mutation operator).The forward algorithm of HMM is already well discussed in my previous paper (Katiyaret al<\/em>., 2020) to calculate the probability of the observation sequence. In their work, we compare both (HMM and GA) results in the bases of\u00a0crystallography data the resolution (in Angstroms).<\/p>\n

Here, we define the Hidden Markov Model as probabilistic models, in which sequences are generated from two simultaneous stochastic processes.This model captures the hidden information from observable sequential symbols (e.g., a nucleotide sequence\u00a0AGCT).This model is defined by states (n), state probabilities (m), transition probabilities (a), emission probabilities (b), and initial probabilities (i). Therefore in a hidden model, there are two stochastic processes: moving between states and the\u00a0process of emitting an output sequence. The sequence of state transitions is a hidden process and is observed through the sequence of emitted symbols(Alghamdi R 2016).<\/p>\n

HMM is characterized by the following:<\/p>\n

N – The number of states in the model.<\/p>\n

M- The number of district observation symbols per state.<\/p>\n

The state transition probability distribution A = (aij<\/sub>)<\/p>\n

The observation symbol probability distribution B = {bj<\/sub>(k)}<\/p>\n

The initial state distribution \u03c0 =\u00a0 (\u03c0i<\/sub>)<\/p>\n

Given applicable values of N, M, A, B, and \u03c0, the HMM can be used as a maker to give an observation sequence.<\/p>\n

Three Basic Problems for HMMs<\/sup><\/strong><\/p>\n

Given the form of HMM of the previous section, there are three basic problems of interest that must be solved for the model to be useful in real-world applications (Mor Bet al.,<\/em> 2020).<\/p>\n

These problems are the following:<\/p>\n

Problem 1: Given the observation sequence O = O1<\/sup> O2<\/sup>\u2026\u2026 O.T. <\/sup>and a model \u00a0\u03bb = (A, B, \u03c0), how do we efficiently compute P (O|\u00a0\u03bb), the probability of the observation sequence, given the model?<\/p>\n

Problem 2: Given the observation sequence O = O1<\/sup> O2<\/sup>\u2026\u2026OT <\/sup>and the model \u00a0\u03bb, how do we choose a corresponding state sequence Q = q1<\/sup>, q2<\/sup>\u2026\u2026\u2026. qT, <\/sup>which is optimal in some meaningful sense (i.e., best “explains” the observations)?<\/p>\n

Problem 3: How do we adjust the model parameters \u00a0\u03bb = (A, B, \u03c0) to maximize P (O|\u00a0\u03bb)?<\/p>\n

Genetic Algorithm (GA)<\/strong><\/p>\n

The genetic algorithm is a method of natural selection that belongs to the class of Evolutionary Algorithms (EA).Genetic Algorithms type of optimization algorithm, meaning they are used to find the optimal solution(s) for a given problem that\u00a0maximizes or minimizes solution of problem(Silvaet al., <\/em>2019).A Genetic Algorithm is the biological process of reproduction and natural random selection to find for the best `fittest’ solution. Like evolution, several of a genetic algorithm’s works randomly\u00a0permit us to set the level of randomization and control (Jenningset al.<\/em>\u00a02019). These are expected to be more powerful algorithms providing random and in-depth search. Such features prove GA to be better than other optimization methods, which have\u00a0drawbacks like lack of stability, derivatives, linearity, or other features. GAareoften designed to simulate a biological process.However, the entities that this terminology refers to in genetic algorithms much simpler than their biological complements (L Halduraiet al., <\/em>2016).<\/p>\n

The basic components of Genetic Algorithms are:<\/p>\n

Function of optimization<\/p>\n

Population of chromosomes<\/p>\n

Random selection of chromosomes<\/p>\n

Crossover to next generation of chromosomes<\/p>\n

Mutation of chromosomes in the new generation<\/p>\n

Some programming languages are also used to implement the GA. They are also well suited for modeling occurrences in economics, ecology, the human immune system, population genetics, and social systems optimization, machine learning optimization (Harsh Bhasinet al.<\/em> 2011) etc.<\/p>\n

The necessary steps involved in the genetic algorithm are (i) generate a random population (suitable and possible solutions for the problem) of chromosomes, (ii) calculate the fitness f(x) of each chromosome (x), (iii) repeat the process until a totallynew population is created, (iv) select two parent chromosomes of better\u00a0fitness, (v) taking into consideration of crossover probability, a crossover is performed between the selected parents (without crossover offspring would be same as a parent), (iv) similarly, with some mutation probability, the new offspring is mutated at each locus, (vii) resulting new offspring is placed in a new population, (viii) the\u00a0newly generated population is used for running with the genetic algorithm, (ix) if end condition is satisfied, run is stopped, the resultant best solution is placed in the current population and (x) go back to step number two for its fitness evaluation.<\/p>\n

A list of some genetic algorithm applications is shown in (Table 1). GA methods play a significant role and provide a useful set of tools in different fields and bioinformatics analysis.<\/p>\n

Table 1: Some common applications of Genetic Algorithm<\/strong><\/p>\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
Domain<\/strong><\/td>\n\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 Application Type<\/strong><\/td>\n<\/tr>\n
Business<\/td>\nControl the gas pipeline, missile evasion<\/td>\n<\/tr>\n
Optimization<\/td>\nMolecular structure optimization, Data compression system, \u00a0stochastic optimization<\/td>\n<\/tr>\n
Management and Engineering<\/td>\nComputer-automated design, vehicle routing problem, sorting network, quality control<\/td>\n<\/tr>\n
Design<\/td>\nPower, electronics aircraft design, keyboard configuration,communication network<\/td>\n<\/tr>\n
Robotics<\/td>\nTrajectory planning<\/td>\n<\/tr>\n
Machine Learning<\/td>\nLearning fuzzy rule base using genetic algorithms, neural networks<\/td>\n<\/tr>\n
Game Playing<\/td>\nPoker, checkers<\/td>\n<\/tr>\n
Security<\/td>\nEncryption, Decryption, and code-breaking<\/td>\n<\/tr>\n
Image Processing<\/td>\nDense pixel matching<\/td>\n<\/tr>\n
Bioinformatics<\/td>\nMultiple sequence alignment, motif discovery, building phylogenetic trees, protein folding, and protein\/ligand docking<\/td>\n<\/tr>\n
Other Applications<\/td>\nMultidimensional system, mutation testing, wireless sensor networks, and clustering using a genetic algorithm.<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

The GA-based methodology provides accuracy, efficiency, and potential for growth to data analysis in bioinformatics. The basic principle behind GAs is that itcreates and maintains a population of individuals represented by chromosomes. Chromosomes are\u00a0essentially a character string analogous to the chromosomes appearing in DNA. These chromosomes are typically encoded solutions to a problem, which then undergo a process of evolution according to rules of reproduction and mutation.<\/p>\n

Material and methods<\/strong><\/p>\n

In this work, we use the NS5 Methyltransferase protein sequence. The protein structure was downloaded from RCSB PDB (http:\/\/www .rcsb.org), and accession codes are used in this study. The coordinates and structure factors for DENV-2 NS5\u00a0 Methyltransferase protein complexes have been deposited in the Protein Data Bank under accession code.<\/p>\n\n\n\n
\"Vol14No3_Gen_Nid_fig1\"<\/a><\/td>\nFigure 1: <\/strong>Flow-chart of Genetic Algorithm.<\/strong><\/p>\n

Click here to view figure\u00a0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Find the optimization value of the function. In the process of the crossover method, first, we select any one population and select any two random points of some length shown in (Figure 1). The same process can be done in the second selected population\u00a0Then we exchange that portion from one population to the other. In mutation, we make minor changes in anyone’s population. Now find the estimated values of both population and compare these estimated values with previously estimated values and\u00a0check which will be better estimate value. At last, select the population contains a better-estimated value.<\/p>\n

Table 2: <\/strong>List of the dengue NS5 methyltransferase proteins with resolution power.<\/strong><\/p>\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
S. No<\/strong><\/td>\nPdb code<\/strong><\/td>\nResolution<\/strong><\/p>\n

Power<\/strong><\/td>\n

Residues<\/strong><\/p>\n

Count<\/strong><\/td>\n

References<\/strong><\/td>\n<\/tr>\n
1.<\/td>\n2P1D<\/td>\n2.90 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Benarroch, D.\u00a0\u00a0 (2002)<\/td>\n<\/tr>\n
2.<\/td>\n2P3O<\/td>\n2.76 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Decroly, E.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2007)<\/td>\n<\/tr>\n
3.<\/td>\n2P3Q<\/td>\n2.75 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Decroly, E.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2007)<\/td>\n<\/tr>\n
4.<\/td>\n2P40<\/td>\n2.70 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Decroly, E.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2007)<\/td>\n<\/tr>\n
5.<\/td>\n1R6A<\/td>\n2.60 \u00c5<\/td>\n295<\/td>\nBenarroch, D., Egloff\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2004)<\/td>\n<\/tr>\n
6.<\/td>\n1L9K<\/td>\n2.40 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Benarroch, D.\u00a0\u00a0 (2002)<\/td>\n<\/tr>\n
7.<\/td>\n4V0Q<\/td>\n2.30 \u00c5<\/td>\n892<\/td>\nZhao, Y., Soh, S., Zheng, \u00a0\u00a0\u00a0\u00a0\u00a0\u00a0(2015)<\/td>\n<\/tr>\n
8.<\/td>\n3EVG<\/td>\n2.20 \u00c5<\/td>\n275<\/td>\nGeiss, B.J.,\u00a0Thompson\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2009)<\/td>\n<\/tr>\n
9.<\/td>\n2P3L<\/td>\n2.20 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Decroly, E.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2007)<\/td>\n<\/tr>\n
10.<\/td>\n5EIF<\/td>\n2.00 \u00c5<\/td>\n552<\/td>\nBenmansour, F.,\u00a0Trist, I.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2016)<\/td>\n<\/tr>\n
11.<\/td>\n4R05<\/td>\n2.00 \u00c5<\/td>\n534<\/td>\nBrecher, M.B.,\u00a0Li, Z.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2015)<\/td>\n<\/tr>\n
12.<\/td>\n3MTE<\/td>\n1.80 \u00c5<\/td>\n444<\/td>\nMacmaster, R.,\u00a0Zelinskaya, N. (2010)<\/td>\n<\/tr>\n
13.<\/td>\n2P41<\/td>\n1.80 \u00c5<\/td>\n305<\/td>\nEgloff, M.P.,\u00a0Decroly, E.\u00a0\u00a0\u00a0\u00a0 (2007)<\/td>\n<\/tr>\n
14.<\/td>\n5E9Q<\/td>\n1.79 \u00c5<\/td>\n512<\/td>\nBenmansour, F.,\u00a0Trist, I.\u00a0\u00a0\u00a0\u00a0 (2016)<\/td>\n<\/tr>\n
15.<\/td>\n3P8Z<\/td>\n1.79 \u00c5<\/td>\n534<\/td>\nLim, S.P.,\u00a0Sonntag, L.S.\u00a0\u00a0\u00a0\u00a0\u00a0 (2011)<\/td>\n<\/tr>\n
16.<\/td>\n5EC8<\/td>\n1.76 \u00c5<\/td>\n552<\/td>\nBenmansour, F.,\u00a0Trist, I.\u00a0\u00a0\u00a0\u00a0 (2016)<\/td>\n<\/tr>\n
17.<\/td>\n5EKX<\/td>\n1.76 \u00c5<\/td>\n552<\/td>\nBenmansour, F.,\u00a0Trist, I.\u00a0\u00a0\u00a0\u00a0 (2016)<\/td>\n<\/tr>\n
18.<\/td>\n3P97<\/td>\n1.74 \u00c5<\/td>\n534<\/td>\nLim, S.P.,\u00a0Sonntag, L.S.\u00a0\u00a0\u00a0\u00a0\u00a0 (2011)<\/td>\n<\/tr>\n
19.<\/td>\n5CUQ<\/td>\n1.74 \u00c5<\/td>\n534<\/td>\nBrecher, M.,\u00a0Chen, H.\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0 (2015)<\/td>\n<\/tr>\n
20.<\/td>\n3MQ2<\/td>\n1.69 \u00c5<\/td>\n436<\/td>\nMacmaster, R.,\u00a0Zelinskaya, N. (2010)<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Implementation of HMM using GA<\/strong><\/p>\n

In the implementation of the Genetic algorithm, we have taken the state transition matrix used in the forward algorithm of the HMM (as described inmy previous paper published) (Katiyaret al<\/em>., 2020). (Table 2) shows the data of 20 gene sequences of\u00a0non-structural methyltransferase protein with resolution power, their number of residues count and relevant references.<\/p>\n

Selection of Population<\/strong><\/p>\n

From (Table 2), we select any two sequences, i.e. (PDB code1 and PDB code2), as the initial population and apply the crossover and mutation function and generate the estimated values of both populations using the forward algorithm of HMM. Then select\u00a0the population PDB code of better-estimated value. Again we select the one sequence of PDB code from (Table 2) and other PDB code, which contains a better-estimated value and applies the same crossover and mutation function now. Once again, find\u00a0which PDB code contains a better estimate value; this process will continue until all the PDB code finished from the (Table 2).<\/p>\n

Crossover<\/strong><\/p>\n

To generate a random population, we applied the crossover operator on the given gene sequences. In this method, we have selected the two gene sequences (among 1\u00a0to 20 protein codes in table 2) of NS5 Methyltransferase protein of dengue virus and interchanged their position.<\/p>\n

Mutation<\/strong><\/p>\n

In this process, we slightly changed a small part of the gene sequence. In a simple sentence, we can say that mutation would change one or more genes, also called as interchanging mutation. After applying the mutation operator, we get the mutated\u00a0gene sequence, and after this, we calculate the optimum value of the gene sequence of using the forward method of the HMM approach.<\/p>\n

In this method, we used the two strings of gene sequences of NS5 methyltransferase protein then applied the crossover and mutation operator. Further,the process of selection, crossover, and mutation operator was performedshown in (Figure 2).<\/p>\n\n\n\n
\"Vol14No3_Gen_Nid_fig2\"<\/a><\/td>\nFigure 2:<\/strong> Selection, Crossover, and Mutation operator Applied.<\/strong><\/p>\n

Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Selected two points randomly of some length in seq1, and the same length at seq2 then applied the crossover operation. In crossover operation, exchange the selected part from both gene sequence seq1 and seq2, and at the result, we have got seq1* and seq2*; and after this, we did mutation operator seq2* and got the seq2**. (In mutation operation here, we have selected a single char ‘G’ and char ‘C’ at the different places of seq2* and then interchange both characters, and get the\u00a0seq2**).The same process is applied in the (Figure 3). In which we have got results after crossover and mutation operation.<\/p>\n\n\n\n
\"Vol14No3_Gen_Nid_fig3\"<\/a><\/td>\nFigure 3:<\/strong> Result of after mutation operation.<\/strong><\/p>\n

Click here to view figure\u00a0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Following matrix for the forward algorithm of HMMs, where each sequence have individual tables ofA = (aij<\/sub>), B = {bj<\/sub>(k)} and \u03c0 =\u00a0 (\u03c0i<\/sub>) matrix values.<\/p>\n\n\n\n
\"Vol14No3_Gen_Nid_tab3\"<\/a><\/td>\nTable 3<\/strong><\/p>\n

Click here to view table<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Using matrices of A, B, and Pi, we apply the GA and forward algorithm of HMM to find the global maxima value of P(O|\u00a0\u03bb) and store in (Table 4).<\/p>\n

Table\u00a0 4: <\/strong>Showing the value of P(O|\u00a0\u03bb<\/strong>)using the Forward Algorithm and Genetic Algorithm.<\/strong><\/p>\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n\n
S. No<\/strong><\/td>\nGene Sequence after mutation<\/strong><\/td>\nResidues Count<\/strong><\/td>\nP(O|\u00a0\u03bb<\/strong>) value after
\n<\/strong>Forward Algorithm<\/strong><\/td>\n		P(O|\u00a0\u03bb<\/strong>) value after
\n<\/strong>Genetic Algorithm<\/strong><\/td>\n<\/tr>\n
1.<\/td>\n5EIF<\/td>\n980<\/td>\n3.96E-278<\/sup><\/td>\n1.63E-266<\/sup><\/td>\n<\/tr>\n
2.<\/td>\n5EC8<\/td>\n980<\/td>\n3.50E-282<\/sup><\/td>\n1.78E-276<\/sup><\/td>\n<\/tr>\n
3.<\/td>\n5E9Q<\/td>\n902<\/td>\n8.93E-269<\/sup><\/td>\n1.08E-244<\/sup><\/td>\n<\/tr>\n
4.<\/td>\n5CUQ<\/td>\n936<\/td>\n1.85E-259<\/sup><\/td>\n9.86E-271<\/sup><\/td>\n<\/tr>\n
5.<\/td>\n4VOQ<\/td>\n1625<\/td>\n0.00E+00<\/sup><\/td>\n0.00E+00<\/sup><\/td>\n<\/tr>\n
6.<\/td>\n4R05<\/td>\n936<\/td>\n8.06E-270<\/sup><\/td>\n3.31E-266<\/sup><\/td>\n<\/tr>\n
7.<\/td>\n3P97<\/td>\n936<\/td>\n8.59E-280<\/sup><\/td>\n1.76E-272<\/sup><\/td>\n<\/tr>\n
8.<\/td>\n3P8Z<\/td>\n936<\/td>\n1.16E-279<\/sup><\/td>\n8.39E-274<\/sup><\/td>\n<\/tr>\n
9.<\/td>\n3MTE<\/td>\n755<\/td>\n1.40E-210<\/sup><\/td>\n5.51E-206<\/sup><\/td>\n<\/tr>\n
10.<\/td>\n3MQ2<\/td>\n740<\/td>\n3.48E-213<\/sup><\/td>\n1.81E-206<\/sup><\/td>\n<\/tr>\n
11.<\/td>\n1R6A<\/td>\n514<\/td>\n6.86E-300<\/sup><\/td>\n2.86E-284<\/sup><\/td>\n<\/tr>\n
12.<\/td>\n1L9K<\/td>\n530<\/td>\n1.96E-304<\/sup><\/td>\n1.87E-306<\/sup><\/td>\n<\/tr>\n
13.<\/td>\n2P1D<\/td>\n530<\/td>\n5.29E-321<\/sup><\/td>\n9.35E-296<\/sup><\/td>\n<\/tr>\n
14.<\/td>\n3EVG<\/td>\n481<\/td>\n3.44E-281<\/sup><\/td>\n3.34E-278<\/sup><\/td>\n<\/tr>\n
15.<\/td>\n2P3L<\/td>\n530<\/td>\n1.96E-304<\/sup><\/td>\n1.99E-313<\/sup><\/td>\n<\/tr>\n
16.<\/td>\n2P3O<\/td>\n530<\/td>\n1.96E-304<\/sup><\/td>\n1.20E-304<\/sup><\/td>\n<\/tr>\n
17.<\/td>\n2P3Q<\/td>\n530<\/td>\n1.88E-304<\/sup><\/td>\n7.08E-311<\/sup><\/td>\n<\/tr>\n
18.<\/td>\n2P40<\/td>\n481<\/td>\n1.96E-304<\/sup><\/td>\n8.83E-295<\/sup><\/td>\n<\/tr>\n
19.<\/td>\n2P41<\/td>\n530<\/td>\n1.96E-304<\/sup><\/td>\n1.66E-307<\/sup><\/td>\n<\/tr>\n
20.<\/td>\n5EKX<\/td>\n980<\/td>\n1.15E-282<\/sup><\/td>\n2.74E-274<\/sup><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

Results and Discussion<\/strong><\/p>\n

This paper presents a comparative result of HMM and GA with the resolution power (crystallographic value of protein sequence given in Table 2).Table3 shows the value of P(O|\u00a0\u03bb) for the gene sequences taken in this research work after the crossover and mutation operation of the GA and forward algorithm of HMM. The result also depends\u00a0on the length of the gene sequence. We have taken 20 gene sequences of NS5methyltransferase protein with residues count (i.e.,the total number of A, G, C, and T characters in sequence) and P (O|\u00a0\u03bb) value of the forward algorithm of HMM and\u00a0P (O|\u00a0\u03bb) of the value of forwarding algorithm after applying GA method. We observed that sequences that have the minimum values of protein 3MTE are1.40E-210<\/sup>,and the maximum value of protein 1L9K is 1.96E-304<\/sup> in the forwarding algorithm. But after GA protein code 2P3Q has the minimum value 7.08E-311 <\/sup>while the maximum values of\u00a0protein 3MQ2 is1.81E-206<\/sup>. Using Table 3, we are comparing and showing the value after the forward algorithm and after GA in (Figure4).<\/p>\n\n\n\n
\"Vol14No3_Gen_Nid_fig4\"<\/a><\/td>\nFigure 4: <\/strong>Graphicimage of the comparative analysisof P(O|<\/strong>\ud835\udecc<\/strong>) values between HMM and Genetic Algorithm for 20 gene sequences of proteins.<\/strong><\/p>\n

Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n

(Figure4) shown thecomparative graph of between HMM and GA to find which method shows the global maxima. Here we observed that in GA, global maxima values lie between (9.85 to 9.35) while in HMM, global maxima values lie between (8.93 to\u00a08.59).So we can say for this experiment, the values of GA for global maxima always give a better result as compared to the HMM forward algorithm. So, here we can state that we can use the global maximum values of both the algorithms so as to get the knowledge of which protein to use among the given four protein as a drug target (in\u00a0case of drug design). Thus, we obtained the better protein in all four proteins which have maxima values, can take any protein which has the highest value of, i.e., (9.86).<\/p>\n

In the case study, we discussed the crystallographic data of the resolution (in Angstroms) of a protein structure resolved by X-ray diffraction or Nuclear Magnetic\u00a0Resonance (NMR). This field can be queried for a value or a range of values. But here, we apply the computational methods to find the resolution power of a protein sequence to take the large data of protein sequence.Most of the structure data are\u00a0obtained from X-ray crystallography and NMR-spectroscopy.\u00a0X-ray crystallography\u00a0determines the preparation of atoms within a protein by passing X-rays through a crystallized form of the protein and analyzing the resulting X-ray\u00a0diffraction pattern. This technique provides the highest resolution and usually yields only one model of a structure.\u00a0Nuclear magnetic resonance (NMR)\u00a0determines the structure of a protein in solution and generally yields multiple models, which allow\u00a0for a description of the biomolecule’s signaling solution(Madej T, Lanczycki CJ\u00a0\u00a0et al<\/em>., 2014).In addition to these experimental methods, some researchers use\u00a0computational modeling\u00a0to predict the structure of a protein by simulating the\u00a0forces that act on each atom in a molecule of known structure. However, this method produces non-experimental models and the least dependable results.<\/p>\n

Conclusion<\/strong><\/p>\n

In this work, we are studying and applying the GA and forward algorithm of HMM on the protein sequence of NS5 Methyltransferase protein of the dengue virus. At last, we compared the estimated value of the forward algorithm of HMM with the GA\u00a0approach. In our experiment, we observed that the genetic algorithm provides a better result as compared to the forward algorithm. In the future, we will repeat this experiment with a large number of gene sequences with higher lengths and will\u00a0compare the results of different algorithms like GA, Metropolis, GWW, and Ant Colony algorithms to evaluate their capabilities to find the global maxima. In a further study, we will improve the algorithm and make it more effective for long protein sequence\u00a0prediction using a multi-core computing platform used other machine learning approaches as an above mention\u00a0for the biological data generation in the analysis and discovery for the drug design.<\/p>\n

Acknowledgement<\/strong><\/p>\n

This study is part of Doctoral research work, and both authors read and approved the final manuscript.<\/p>\n

Conflict of Interest<\/strong><\/p>\n

All the authors declare that there are no conflicts of interest regarding the publication of this research paper.<\/p>\n

Funding Source<\/strong><\/p>\n

This work did not provide any source of financial support.<\/p>\n

References<\/strong><\/p>\n

    \n
  1. Anasir M IRamanathan B and Poh C L. Structure-Based Design of Antivirals against Envelope Glycoprotein of Dengue Virus\u00a0Viruses<\/em>., 2020; 12:4: 367.
    \n    CrossRef<\/a><\/li>\n
  2. Amjad M K Butt S I Kousar R Ahmad Agha M H Faping Z Anjum N Asgher U. Recent Research Trends in Genetic Algorithm Based Flexible Job Shop Schedulling Problems Mathematical\u00a0 Problems in Engineering Article ID 9270802., 2018; 32.
    \n    CrossRef<\/a><\/li>\n
  3. Alghamdi R. Hidden Markov Models (HMMs) and security applications Int J AdvComputSci Appl., 2016; 7: 2: 39-47.
    \n    CrossRef<\/a><\/li>\n
  4. Benmansour F Trist I Coutard B Decroly EQuerat G Brancale A and Barral K. Discovery of novel dengue virus NS5 methyltransferase non-nucleoside inhibitors by fragment-based drug design European journal of medicinal chemistry., 2017; 125: 865-880.
    \n    CrossRef<\/a><\/li>\n
  5. Brecher M B Li Z Zhang J ChenH Lin Q Liu B and Li H. Refolding of a fully functional flavivirus methyltransferase revealed that S\u2010adenosyl methionine but not S\u2010adenosylhomocysteine is copurified with flavivirus methyltransferase\u00a0Protein Science<\/em>., 2015; 24: 1: 117-128.
    \n    CrossRef<\/a><\/li>\n
  6. Chuang C H Chiou S J Cheng T L and Wang Y T. A molecular dynamics simulation study decodes the Zika virus NS5 methyltransferase bound to SAH and RNA analogue Scientific reports<\/em>., 2018; 8: <\/em>1: 6336.
    \n    CrossRef<\/a><\/li>\n
  7. El Sahili A and Lescar J. Dengue Virus Non-Structural Protein 5\u00a0Viruses<\/em>., 2017;\u00a09: 4: 91.
    \n    CrossRef<\/a><\/li>\n
  8. Herrero L J et al. , 2013; <\/em>137: 266-82.<\/li>\n
  9. Harsh Bhasinet al.<\/em>Application of Genetic Algorithms in Machine learning International Journal of Computer Science and Information Technologies (IJCSIT)<\/em>. , 2011; 2: 5: 2412-2415.<\/li>\n
  10. Jennings C Lysgaard S Hummelsh\u00f8j JS\u00a0et al.<\/em>Genetic algorithms for computational materials discovery accelerated by machine learning\u00a0npjComput Mater<\/em>., 2019;\u00a05: 46.
    \n    CrossRef<\/a><\/li>\n
  11. Katiyar, Nidhi, Nath, Ravindra and Katiyar, Shashwat. Estimated Value of Hidden Markov Model Parameters for NS5 methyltransferase Protein of Dengue Virus.International Journal on Emerging Technologies<\/em>., 2020; 11:<\/em> 1: 01\u201311.<\/li>\n
  12. L Haldurai T Madhubala and R Rajalakshmi. A Study on Genetic Algorithms and its Applications In ternational Journal of Computer sciences and Engineering., 2016; 4: 10: 139-143.<\/li>\n
  13. Lim S P Sonntag L S Noble C Nilar S H Ng R H Zou G and Bodenreider C. Small molecule inhibitors that selectively block dengue virus methyltransferase Journal of Biological Chemistry., 2011; 286: 8: 6233-6240.
    \n    CrossRef<\/a><\/li>\n
  14. Macmaster R Zelinskaya N Savic M Rankin C R and Conn G L. Structural insights into the function of aminoglycoside-resistance A1408 16S rRNA methyltransferases from antibiotic-producing and human pathogenic bacteria Nucleic acids research., 2010; 38: 21:7791-7799.
    \n    CrossRef<\/a><\/li>\n
  15. Madej T Lanczycki CJ Zhang D Thiessen PA Geer RC Marchler-Bauer A Bryant SH. MMDB and VAST+:tracking structural similarities between macromolecular complexesNucleic Acids Res., 2014;42: 1:297-303.
    \n    CrossRef<\/a><\/li>\n
  16. Mor B Garhwal S and Kumar A. A Systematic Review of Hidden Markov Models and Their Applications Arch Computational Methods Eng<\/em>., 2020; https:\/\/doi.org\/10.1007\/s11831-020-09422-4.
    \n    CrossRef<\/a><\/li>\n
  17. V. Chandrika, K. Visalakshmi and K. Sakthi Srinivasan Application of Hidden Markov Models in Stock Trading 2020 6th International Conference on Advanced Computing and Communication Systems (ICACCS)<\/em> Coimbatore India., 1144-1147.
    \n    CrossRef<\/a><\/li>\n
  18. Qiong C Huang M Xu Q Wang H Wang J. Reinforcement learning-Based Genetic Algorithm in Optimizing Multidimensional Data Discretizatio Scheme Mathematical Problems in Engineering., 2020; 13.
    \n    CrossRef<\/a><\/li>\n
  19. Silva FT Silva MX and Belchior JC (2019). A New Genetic Algorithm Approach Applied to Atomic and Molecular Cluster Studies.
    \n    CrossRef<\/a><\/li>\n
  20. S. F. Tambunan et al. Drug Target Insights., 2017; 1-11.
    \n    CrossRef<\/a><\/li>\n
  21. VanniceKS Durbin AHombach J. Status of vaccine research and development of vaccines for dengueVaccine<\/em>., 2016; 34: 2934\u20138.
    \n    CrossRef<\/a><\/li>\n
  22. Yang Q M Yang Jack Y. Lecture notes and beyond the decade of high-performance computing for the next-generation sequence analysis I J Computational Biology and Drug Design., 2010; 2: 2: 204-206.
    \n    CrossRef<\/a><\/li>\n
  23. Zhao Y Soh T S Zheng J Chan K W K Phoo W W Lee C C and Shi P Y. A crystal structure of the dengue virus NS5 protein reveals a novel inter-domain interface essential for protein flexibility and virus replication\u00a0PLoS pathogens., 2015; <\/em>11: 3: 1-27.
    \n    CrossRef<\/a><\/li>\n
  24. Zhou YH Brooks P Wang X. A two-stage Hidden Markov Model design for biomarker detection, with application to microbiome research Stat Biosci., 2018; 10: 1-18.
    \n    CrossRef<\/a><\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"

    Introduction Dengue virus (DV), the causative agent of dengue, resides […]<\/p>\n","protected":false},"author":14,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[93],"tags":[],"class_list":["post-40106","post","type-post","status-publish","format-standard","hentry","category-vol14no3"],"_links":{"self":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/40106","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/users\/14"}],"replies":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/comments?post=40106"}],"version-history":[{"count":5,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/40106\/revisions"}],"predecessor-version":[{"id":41227,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/40106\/revisions\/41227"}],"wp:attachment":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/media?parent=40106"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/categories?post=40106"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/tags?post=40106"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}