The present invention relates to information processing programs, for example.
Genetic recombination manipulations have been performed using gene vectors due to advancement of gene transfer technologies and deeper understanding related immunological mechanisms. Mediums with various added characteristics are used as gene vectors differently depending on sizes of gene segments inserted and purposes of insertion. Gene vectors derived from host organisms, such as colon bacilli and yeasts, are used for these manipulations.
For example, chimeric antigen receptor (CAR) introduced T cell therapy has attracted attention as immunotherapy of cancer using genetically modified T cells. A CAR is a receptor that: is artificially made by fusion of a part that specifically recognizes an antigen and that is derived from an antibody with a part derived from a T cell receptor (TCR), the part having a cytotoxic function; specifically recognizes a cancer antigen; and is able to attack the cancer antigen.
Developing gene therapy drugs using gene vectors is very promising but synthesizing gene therapy drugs using various gene vectors, as is, is difficult.
Accordingly, one may consider synthesizing target gene therapy drugs by substitution of various gene vectors, but searching for substitutable gene vectors and efficient genetic modification are actually difficult.
According to an aspect of the embodiment of the invention, a non-transitory computer-readable recording medium stores therein an information processing program that causes a computer to execute a process including executing training of a trained model based on training data defining relations between vectors corresponding to genomes and vectors respectively corresponding to pluralities of subgenomes composing the genomes and in a case where a genome to be analyzed has been received, calculating vectors of a plurality of subgenomes corresponding to the genome to be analyzed by inputting the genome to be analyzed into the trained model.
The object and advantages of the invention will be realized and attained by means of the elements and combinations particularly pointed out in the claims.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory and are not restrictive of the invention.
Embodiments of an information processing program, an information processing method, and an information processing apparatus disclosed in the present application will hereinafter be described in detail on the basis of the drawings. The present invention is not limited by these embodiments.
Genomes will be described before description of an embodiment.
DNAs and RNAs each have four types of bases that are each denoted by a symbol, “A”, “G”, “C”, “T”, or “U”. Furthermore, sequences each made up of three bases determine 20 types of amino acids. Each of these amino acids is denoted by symbols, “A” to “Y”.
As illustrated in
A protein is also uniquely determined by a sequence of bases. A primary structure of a protein is a sequence of plural amino acids. A secondary structure includes α-helixes and β-sheets that are symmetrical substructures observed locally. A tertiary structure includes plural secondary structures. Furthermore, a higher-order structure includes plural tertiary structures.
A gene vector used in this embodiment is a DNA or RNA molecule that is used to artificially carry a foreign genetic substance to another cell. Gene vectors include, for example, plasmids, cosmids, lambda phages, and artificial chromosomes.
Furthermore, the gene vector is generated by synthesis of plural subvectors. Subvectors are base sequences of DNAs and RNAs and correspond to, for example, secondary structures of proteins described by reference to
An example of a process by an information processing apparatus according to this embodiment will be described next.
The training data 65 define relations between vectors of target genomes (therapeutic drugs) and vectors of pluralities of subgenomes included in the target genomes. For example, a vector of a target genome corresponds to input data and a plurality of subgenomes serves as correct answer values of output data therefor.
The information processing apparatus executes training by error back propagation so that output upon input of a vector of a target genome into the trained model approaches vectors of its subgenomes. The information processing apparatus adjusts parameters of the trained model 70 (executes machine training) by repeatedly executing the above described process on the basis of the relations included in the training data 65, the relations each being between: a vector of a target genome; and vectors of a plurality of subgenomes.
In response to the information processing apparatus receiving an analysis query 80 that specifies a target genome (a therapeutic drug), the information processing apparatus converts the target genome in the analysis query 80 to a vector Vob80. By inputting the vector Vob80 to the trained model 70, the information processing apparatus calculates a plurality of vectors (Vsb80-1, Vsb80-2, Vsb80-3, . . . , Vsb80-n) corresponding to its subgenomes and stores the calculated plurality of vectors into a subgenome table T1.
The information processing apparatus makes a comparison among degrees of similarity between plural vectors (Vt1, Vt2, Vt3, . . . , Vtn) corresponding respectively to alternative gene vectors stored in an alternative gene vector table T2 and the plurality of vectors (Vsb80-1, Vsb80-2, Vsb80-3, . . . , Vsb80-n) to determine vectors of similar alternative gene vectors. The information processing apparatus registers the vector of the target genome, the vectors of the subgenomes, and the vectors of the similar alternative gene vectors, in association with one another, into an alternative management table 85.
As described above, the information processing apparatus according to the embodiment executes training of the trained model 70 beforehand, on the basis of the training data 65 defining the relations between the vectors of the target genomes and the vectors of their subgenomes. By inputting a vector of an analysis query into the trained model 70 that has been trained, the information processing apparatus calculates vectors of subgenomes corresponding to the target compound in the analysis query. Using the vectors of the subgenomes output from the trained model 70 facilitates detection of substitutable gene vectors that are gene vectors similar to the subgenomes included in the target genome.
An example of a configuration of the information processing apparatus according to the first embodiment will be described next.
The communication unit 110 is connected to, for example, an external device by wire or wirelessly and transmits and receives information to and from, for example, the external device. For example, the communication unit 110 is implemented by a network interface card (NIC). The communication unit 110 may be connected to a network not illustrated in the drawings.
The input unit 120 is an input device that inputs various types of information to the information processing apparatus 100. The input unit 120 corresponds to, for example, a keyboard and a mouse, or a touch panel.
The display unit 130 is a display device that displays information output from the control unit 150. The display unit 130 corresponds to, for example, a liquid crystal display, an organic electro luminescence (EL) display, or a touch panel.
The storage unit 140 has a base file 50, a conversion table 51, a dictionary table 52, a compressed file table 53, a vector table 54, and an inverted index table 55. Furthermore, the storage unit 140 has the subgenome table T1, the alternative gene vector table T2, a genome dictionary D2, the training data 65, the trained model 70, the analysis query 80, and the alternative management table 85. The storage unit 140 is implemented by, for example: a semiconductor memory element, such as a random access memory (RAM) or a flash memory; or a storage device, such as a hard disk or an optical disk.
The base file 50 is a file that holds information including a sequence of plural bases.
The conversion table 51 is a table associating codons and codes of the codons with each other. A sequence of three bases is called a “codon”.
The dictionary table 52 is a table that holds various dictionaries.
The protein primary structure dictionary D1-1 is dictionary data defining relations between compressed codes of proteins and sequences of codons composing the proteins.
For example, a compressed code, “C0008000h”, is assigned to a protein primary structure, “type I collagen”. A codon code sequence corresponding to the compressed code, “C0008000h”, is “02h63h78h . . . 03h”.
The secondary structure dictionary D1-2 is dictionary data defining relations between sequences of compressed codes of protein primary structures and compressed codes of secondary structures.
For example, a compressed code, “D0000000h”, is assigned to a secondary structure, “a secondary structure”. A protein primary structure code sequence corresponding to the compressed code, “D0000000h”, is “C0008001hC00 . . . ”.
The tertiary structure dictionary D1-3 is dictionary data defining relations between sequences of compressed codes of secondary structures and compressed codes of tertiary structures.
For example, a compressed code, “E0000000h”, is assigned to a tertiary structure, “aa tertiary structure”. A secondary structure code sequence corresponding to the compressed code, “E0000000h”, is “D0008031hD00 . . . ”.
The higher-order structure dictionary D1-4 is dictionary data defining relations between sequences of compressed codes of tertiary structures and compressed codes of higher-order structures.
For example, a compressed code, “F0000000h”, is assigned to a higher-order structure, “aaa higher-order structure”. A tertiary structure code sequence corresponding to the compressed code, “F0000000h” is “E0000031hE00 . . . ”.
The description of
The codon compression file 53A is a file having the bases compressed in units of codons, the bases being included in the base file 50.
The protein primary structure compression file 53B is a file having the sequences coded in units of protein primary structures, the sequences being of the compressed codes of the codons included in the codon compression file 53A.
The secondary structure compression file 53C is a file having the sequences coded in units of secondary structures, the sequences being of the compressed codes of the protein primary structures included in the protein primary structure compression file 53B.
The tertiary structure compression file 53D is a file having the sequences coded in units of tertiary structures, the sequences being of the compressed codes of the secondary structures included in the secondary structure compression file 53C.
The higher-order structure compression file 53E is a file having the sequences coded in units of higher-order structures, the sequences being of the compressed codes of the tertiary structures included in the tertiary structure compression file 53D.
The vector table 54 is a table that holds vectors corresponding to protein primary structures, secondary structures, tertiary structures, and higher-order structures.
The protein primary structure vector table VT1-1 is a table that holds vectors corresponding to protein primary structures.
The secondary structure vector table VT1-2 is a table that holds vectors corresponding to secondary structures.
The tertiary structure vector table VT1-3 is a table that holds vectors corresponding to tertiary structures.
The higher-order structure vector table VT1-4 is a table that holds vectors corresponding to higher-order structures.
The description of
For example, the compressed code of the protein primary structure at the head of the protein primary structure compression file 53B has an offset of “0”. In a case where the code, “C0008000h (type I collagen)”, of a protein primary structure is included at the eighth position from the head of the protein primary structure compression file 53B, the bit at a position where the column of the offset of “7” in the protein inverted index In1-1 and the row of the code, “C0008000h (type I collagen)”, of the protein intersect each other is “1”.
For example, the compressed code of the secondary structure at the head of the secondary structure compression file 53C has an offset of “0”. In a case where the code, “D000000h (a secondary structure)”, of a secondary structure is included at the eighth position from the head of the secondary structure compression file 53C, the bit at a position where the column of the offset of “7” in the secondary structure inverted index Int-2 and the row of the compressed code, “D0000000h (a secondary structure)”, of the secondary structure intersect each other is “1”.
For example, the compressed code of the tertiary structure at the head of the tertiary structure compression file 53D has an offset of “0”. In a case where the code, “E0000000h (aa tertiary structure)”, of a tertiary structure is included at the eleventh position from the head of the tertiary structure compression file 53D, the bit at a position where the column of the offset of “10” in the tertiary structure inverted index Int-3 and the row of the compressed code, “E0000000h (aa tertiary structure)”, of the tertiary structure intersect each other is “1”.
For example, the compressed code of the higher-order structure at the head of the higher-order structure compression file 53E has an offset of “0”. In a case where the code, “F0000000h (aaa higher-order structure)”, of a higher-order structure is included at the eleventh position from the head of the higher-order structure compression file 53E, the bit at a position where the column of the offset of “10” in the higher-order structure inverted index Int-4 and the row of the compressed code, “F0000000h (aaa higher-order structure)”, of the higher-order structure intersect each other is “1”.
The description of
The genome dictionary D2 defines relations between names of target genomes and names of subgenomes included in these target genomes.
The training data 65 define relations each between a vector of a target genome and vectors of a plurality of subgenomes included in that target genome. The training data 65 have a data structure corresponding to the data structure of the training data described by reference to
The trained model 70 is a model corresponding to, for example, a CNN or an RNN, and parameters are set for the trained model 70.
The analysis query 80 includes information on a target genome (therapeutic drug) to be analyzed. For example, the information on a target genome includes information on a base sequence corresponding to a higher-order structure.
The alternative management table 85 is a table holding vectors of subgenomes included in target genomes and vectors of gene vectors similar to these subgenomes in association with each other, the gene vectors being able to substitute for the subgenomes.
The control unit 150 has a preprocessing unit 151, a training unit 152, a calculation unit 153, and an analysis unit 154. The control unit 150 is implemented by, for example, a central processing unit (CPU) or a micro processing unit (MPU). Furthermore, the control unit 150 may be implemented by, for example, an integrated circuit, such as an application specific integrated circuit (ASIC) or a field programmable gate array (FPGA).
By executing various processes described below, the preprocessing unit 151 calculates a vector of a higher-order structure or tertiary structure corresponding to a target genome (therapeutic drug) and vectors of secondary structures corresponding to subgenomes, for example.
Firstly, the preprocessing unit 151 executes a process of generating the codon compression file 53A, a process of generating the protein primary structure compression file 53B, and a process of generating the protein primary structure vector table VT1-1 and the protein primary structure inverted index Int-1.
The preprocessing unit 151 compares the base file and the conversion table 51 with each other, assigns compressed codes, in units of codons, to the base sequence in the base file 50, and generates the codon compression file 53A.
The preprocessing unit 151 compares the codon compression file 53A and the protein primary structure dictionary D1-1 with each other, assigns compressed codes, in units of protein primary structures, to sequences of the compressed codes of the codons included in the codon compression file 53A, and generates the protein primary compression file 53B.
In response to the preprocessing unit 151 generating the protein primary structure compression file 53B, the preprocessing unit 151 calculates vectors of the protein primary structures (the compressed codes of the protein primary structures) by embedding the compressed codes of the protein primary structures in Poincare space. A process of calculating a vector by embedding into Poincare space is a technique called Poincare Embeddings. For Poincare Embeddings, for example, a technique described in Non-Patent Literature by Valentin Khrulkovl et al., “Hyperbolic Image Embeddings”, published by Cornell University on Apr. 3, 2019 may be used.
Poincare Embeddings are characterized in that a vector is assigned according to the embedded position in Poincare space and the more similar pieces of information are, the nearer their embedded positions are. Therefore, groups having similar characteristics are embedded at positions that are near one another in Poincare space and similar vectors are thus assigned to these groups. Although illustration will be omitted, the preprocessing unit 151 refers to a protein primary structure similarity table defining protein primary structures that are similar to each other, embeds a compressed code of each protein primary structure into Poincare space, and calculates a vector of the compressed code of each protein primary structure. The preprocessing unit 151 may execute Poincare Embeddings beforehand for a compressed code of each protein primary structure defined in the protein dictionary primary structure D1-1.
The preprocessing unit 151 generates the protein primary structure vector table VT1-1 by associating protein primary structures (compressed codes of protein primary structures) and vectors of the protein primary structures with each other. The preprocessing unit 151 generates the protein primary structure inverted index In1-1 on the basis of relations between the vectors of the protein primary structures and positions of the protein primary structures (the compressed codes of the protein primary structures) in the protein primary structure compression file 53B.
Subsequently, the preprocessing unit 151 executes a process of generating the secondary structure compression file 53C and a process of generating the secondary structure vector table VT1-2 and the secondary structure inverted index Int-2.
The preprocessing unit 151 compares the protein primary structure compression file 53B and the secondary structure dictionary D1-2 with each other, assigns compressed codes, in units of secondary structures, to sequences of the compressed codes of the protein primary structures included in the protein primary structure compression file 53B and generates the secondary structure compression file 53C.
The preprocessing unit 151 refers to the secondary structure dictionary D1-2 and determines protein primary structure code sequences (sequences of compressed codes of protein primary structures) corresponding to a compressed code of a secondary structure. The preprocessing unit 151 obtains vectors of the determined compressed codes of the protein primary structures from the protein primary structure vector table VT1-1, adds up the obtained vectors, and thereby calculates a vector of the compressed code of the secondary structure. By repeatedly executing the above described process, the preprocessing unit 151 calculates vectors of secondary structures.
By associating the secondary structures (the compressed codes of the secondary structures) and the vectors of the secondary structures with each other, the preprocessing unit 151 generates the secondary structure vector table VT1-2. The preprocessing unit 151 generates the secondary structure inverted index Int-2 on the basis of relations between the vectors of the secondary structures and positions of the secondary structures (the compressed codes of the secondary structures) in the secondary structure compression file 53C.
Subsequently, the preprocessing unit 151 executes a process of generating the tertiary structure compression file 53D and a process of generating the tertiary structure vector table VT1-3 and the tertiary structure inverted index Int-3.
The preprocessing unit 151 compares the secondary structure compression file 53C and the tertiary structure dictionary D1-3 with each other, assigns compressed codes, in units of tertiary structures, to sequences of the compressed codes of the secondary structures included in the secondary structure compression file 53C.
The preprocessing unit 151 refers to the tertiary structure dictionary D1-3 and determines secondary structure code sequences (sequences of compressed codes of secondary structures) corresponding to a compressed code of a tertiary structure. The preprocessing unit 151 obtains vectors of the determined compressed codes of the secondary structures from the secondary structure vector table VT1-2, adds up the obtained vectors, and thereby calculates a vector of the compressed code of the tertiary structure. By repeatedly executing the above described process, the preprocessing unit 151 calculates vectors of tertiary structures.
By associating the tertiary structures (the compressed codes of the tertiary structures) and the vectors of the tertiary structures with each other, the preprocessing unit 151 generates the tertiary structure vector table VT1-3. The preprocessing unit 151 generates the tertiary structure inverted index Int-3 on the basis of relations between the vectors of the tertiary structures and positions of the tertiary structures (the compressed codes of the tertiary structures) in the tertiary structure compression file 53D.
Subsequently, the preprocessing unit 151 executes a process of generating the higher-order structure compression file 53E and a process of generating the higher-order structure vector table VT1-4 and the higher-order structure inverted index Int-4.
The preprocessing unit 151 compares the tertiary structure compression file 53D and the higher-order structure dictionary D1-4 with each other, assigns compressed codes, in units of higher-order structures, to sequences of the compressed codes of the tertiary structures included in the tertiary structure compression file 53D, and generates the higher-order structure compression file 53E.
The preprocessing unit 151 refers to the higher-order structure dictionary D1-4 and determines tertiary structure code sequences (sequences of compressed codes of tertiary structures) corresponding to a compressed code of a higher-order structure. The preprocessing unit 151 obtains vectors of the determined compressed codes of the tertiary structures from the tertiary structure vector table VT1-3, adds up the obtained vectors, and thereby calculates a vector of the compressed code of the higher-order structure. By repeatedly executing the above described process, the preprocessing unit 151 calculates vectors of higher-order structures.
The preprocessing unit 151 generates the higher-order structure vector table VT1-4 by associating the higher-order structures (the compressed codes of the higher-order structures) and the vectors of the higher-order structures with each other. The preprocessing unit 151 generates the higher-order structure inverted index Int-4 on the basis of relations between the vectors of the higher-order structures and positions of the higher-order structures (the compressed codes of the higher-order structures) in the higher-order structure compression file 53E.
The following description is on an example of a process in which the preprocessing unit 151 generates the alternative gene vector table T2. For example, the preprocessing unit 151 sets, as is, the vectors of the tertiary structures included in the secondary structure vector table VT1-2, in the alternative gene vector table T2. In a case where the preprocessing unit 151 has received specification of a vector via the input unit 120, the preprocessing unit 151 may set the specified vector in the alternative gene genome table T2.
The following description is on an example of a process in which the preprocessing unit 151 generates the training data 65. The preprocessing unit 151 determines a relation between a name of a target genome and names of subgenomes, on the basis of the genome dictionary D2. The preprocessing unit 151 determines a vector of the target genome on the basis of the higher-order structure dictionary D1-4 and the higher-order structure vector table VT1-4, or the tertiary structure dictionary D1-4 and the tertiary structure vector table VT1-3, and the name of the target genome. The preprocessing unit 151 determines vectors of the subgenomes on the basis of the secondary structure dictionary D1-2 and secondary structure vector table VT1-2 and the names of the subgenomes. The preprocessing unit 151 determines the relation between the target genome and the subgenomes through this process and registers the relation into the training data 65.
The preprocessing unit 151 generates the training data 65 by repeatedly executing the above described process. The information processing apparatus 100 may obtain and use the training data 65 that have been generated beforehand, from, for example, an external device.
The description of
The training unit 152 executes training of the trained model 70 by repeatedly executing the above described process for pairs of vectors of target genomes and vectors of subgenomes in the training data 65.
In a case where the calculation unit 153 has received specification by the analysis query 80, the calculation unit 153 calculates vectors of subgenomes included in the target genome in the analysis query 80, by using the trained model 70 that has been trained. A process by the calculation unit 153 corresponds to the process described by reference to
The calculation unit 153 obtains a base sequence of the target genome included in the analysis query 80. The calculation unit 153 compares the base sequence of the target genome and the conversion table 51 with each other, determines codons included in the base sequence of the target genome, and converts the base sequence of the target genome into compressed codes, in units of codons. Furthermore, the calculation unit 153 compares the codon code sequences compressed in units of codons and the protein primary structure dictionary D1-1 with each other and converts the codon code sequences into compressed codes, in units of protein primary structures.
The calculation unit 153 compares the converted compressed codes of the protein primary structures with the protein primary structure vector table VT1-1 and determines vectors of the compressed codes of the protein primary structures. By adding up the determined vectors of the compressed codes of the protein primary structures, the calculation unit 153 calculates the vector Vob80 corresponding to the target genome included in the analysis query 80.
In a case where a target genome has been specified by secondary structures of plural subgenomes, the calculation unit 153 executes the following process. The calculation unit 153 compares the secondary structures of the subgenomes of the target genome with the secondary structure dictionary D1-2 and secondary structure vector table VT1-2, and determines vectors of the secondary structures of the subgenomes included in the target genome. By adding up the determined vectors of the secondary structures of the subgenomes, the calculation unit 153 calculates a vector of the target genome.
The calculation unit 153 calculates plural vectors corresponding to the subgenomes by inputting the vector Vob80 into the trained model 70. The calculation unit 153 outputs the calculated vectors of the subgenomes, to the analysis unit 154. In the description hereinafter, the vectors of the subgenomes calculated by the calculation unit 153 will be respectively referred to as “analysis vectors”. The calculation unit 153 stores the vectors (analysis vectors) of the subgenomes into the subgenome table T1.
The analysis unit 154 makes a search for information on alternative gene vectors having vectors similar to the analysis vectors, on the basis of the analysis vectors. On the basis of a result of the search, the analysis unit 154 registers the vectors of the subgenomes included in the target genome and the vectors (similar vectors described hereinafter) of the alternative gene vectors similar thereto, in association with each other, into the alternative management table 85.
For example, the analysis unit 154 calculates distances between an analysis vector and the vectors included in the alternative gene vector table T2 to determine any vector having a distance less than a threshold, the distance being from the analysis vector. Any vector included in the alternative gene vector table T2 and having the distance from the analysis vector is a “similar vector”, the distance being less than the threshold. A gene vector corresponding to this similar vector is a substitutable gene vector.
The analysis unit 154 may determine a compressed code of the gene vector corresponding to the similar vector, on the basis of the secondary structure vector table VT1-2 and determine a protein primary structure included in the gene vector, on the basis of the determined compressed code of the gene vector, the secondary structure dictionary D1-2, and the protein primary structure dictionary D1-1. By executing this process, the analysis unit 154 makes a search for characteristics of the substitutable gene vector corresponding to the similar vector and registers the characteristics into the alternative management table 85. The characteristics of the substitutable gene vector are the protein included in the gene vector and the primary structure of the protein.
By repeatedly executing the above described process for the analysis vectors, the analysis unit 154 may make a search, for each of the analysis vectors, for characteristics of the gene vector corresponding to the similar vector and register the characteristics into the alternative management table 85. The analysis unit 154 may output the alternative management table 85 to the display unit 130 to cause the display unit 130 to display the alternative management table 85 or may transmit the alternative management table 85 to an external device connected to a network.
An example of a procedure by the information processing apparatus 100 according to the embodiment will be described next.
On the basis of the base file 50, the conversion table 51, and the dictionary table 52, the preprocessing unit 151 generates the compressed file table 53, the vector table 54, and the inverted index table 55 (Step S102).
The preprocessing unit 151 generates the training data 65 (Step S103). On the basis of the training data 65, the training unit 152 of the information processing apparatus 100 executes training of the trained model 70 (Step S104).
The calculation unit 153 calculates vectors of subgenomes by inputting the calculated vector of the analysis query 80 into the trained model 70 that has been trained (Step S203). The analysis unit 154 of the information processing apparatus 100 compares the vectors of the subgenomes and the alternative gene vector table T2 with each other (Step S204).
The analysis unit 154 makes a search for substitutable gene vectors corresponding to the subgenomes (Step S205). The analysis unit 154 registers a result of the search into the alternative management table 85 (Step S206).
Effects of the information processing apparatus 100 according to the embodiment will be described next. In the training phase, the information processing apparatus 100 executes training of the trained model 70 beforehand, on the basis of the training data 65 defining relations between vectors of target compounds (therapeutic drugs) and vectors of subgenomes. In the analysis phase, the information processing apparatus 100 calculates vectors of subgenomes corresponding to an analysis query (a target genome) by inputting the vector of the analysis query into the trained model 70 that has been trained. Using the vectors of the subgenomes output from the trained model 70 facilitates detection of substitutable gene vectors similar to the subgenomes included in the target genome.
For example, in a case where a subgenome included in a target genome is a rare subgenome, executing the process by the information processing apparatus 100 facilitates the search for an inexpensive gene vector serving as an alternative to that subgenome.
In the above described embodiment, the comparison is made at granularity of subgenomes (secondary structures) for a search to be made for a substitutable gene vector, but the embodiment is not limited to this example. For example, the information processing apparatus 100 may make a comparison at granularity of plural primary structures composing a subgenome for a search to be made for alternative primary structures.
A second embodiment will be described next.
The training data 90 define relations between vectors of pluralities of subgenomes for synthesis of target genomes (therapeutic drugs) and vectors of common structures maintained in genetic modification based on gene vectors. For example, vectors of subgenomes correspond to input data, and vectors of plural common structures serve as correct answer values.
The information processing apparatus executes training by error back propagation, so that output upon input of a vector of a subgenome to the trained model 91 approaches the vector of each common structure. The information processing apparatus adjusts parameters of the trained model 91 (executes machine training) by repeatedly executing the above described process on the basis of the relations between the vectors of the subgenomes included in the training data 90 and the vectors of the common structures.
In response to the information processing apparatus receiving the analysis query 92 specifying the subgenome, the information processing apparatus converts the subgenome in the analysis query 92 into a vector Vsb92-1 by using a subgenome vector table T1. By inputting the vector Vsb92-1 of the subgenome into the trained model 91, the information processing apparatus calculates a vector Vcm92-1 corresponding to its common structure.
The information processing apparatus compares the vector Vsb92-1 of the subgenome and vectors of plural gene vectors included in an alternative gene vector table T2 with each other. The alternative gene vector table T2 corresponds to the alternative gene vector table T2 described with respect to the first embodiment.
The information processing apparatus determines a vector of a gene vector similar to the vector Vsb92-1 of the subgenome. For example, a vector Vt92-1 is determined as the vector of the gene vector similar to the vector Vsb92-1 of the subgenome. A vector of a common structure common to the subgenome having the vector Vsb92-1 and the gene vector having the vector Vt92-1 is then found to be the vector Vcm92-1 output from the trained model 91. Furthermore, a result of subtraction of the vector Vcm92-1 of the common structure from the vector Vt92-1 of the gene vector is a vector of a “genetically modified structure” corresponding to difference between the similar gene vector and the subgenome.
The information processing apparatus registers the relation between the vector of the common structure and the vector of the genetically modified structure into a common structure and genetically modified structure table 93. By repeatedly executing the above described process for vectors of subgenomes, the information processing apparatus generates the common structure and genetically modified structure table 93.
As described above, the information processing apparatus according to the second embodiment inputs the vector of the analysis query 92 into the trained model 91 that has been trained to calculate a vector of each common structure corresponding to the subgenome in the analysis query. Furthermore, by subtraction of a vector of a common structure from each of vectors of gene vectors similar to a subgenome, a vector of a genetically modified structure corresponding to difference between the similar subgenome and the gene vector is calculated. Using the above described vector of the common structure and vector of the genetically modified structure facilitates analysis for a better gene vector usable for synthesis and manufacture of the target genome.
An example of a configuration of the information processing apparatus according to the second embodiment will be described next.
Description related to the communication unit 210, the input unit 220, and the display unit 230 is similar to the description related to the communication unit 110, the input unit 120, and the display unit 130 described with respect to the first embodiment.
The storage unit 240 has a base file 50, a conversion table 51, a dictionary table 52, a compressed file table 53, a vector table 54, and an inverted index table 55. Furthermore, the storage unit 240 has the subgenome table T1, the alternative gene vector table T2, a genome dictionary D2, the training data 90, the trained model 91, the analysis query 92, and the common structure and genetically modified structure table 93. The storage unit 240 is implemented by, for example: a semiconductor memory element, such as a random access memory (RAM) or a flash memory; or a storage device, such as a hard disk or an optical disk.
Description related to the base file 50, the conversion table 51, the dictionary table 52, the compressed file table 53, the vector table 54, the inverted index table 55, the subgenome table T1, the alternative gene vector table T2, and the genome dictionary D2 is similar to what has been described with respect to the first embodiment. The training data 90 are similar to those described by reference to
As described by reference to
The description of
Description related to the preprocessing unit 251 is similar to the description of the process related to the preprocessing unit 151 described with respect to the first embodiment. The base file 50, the conversion table 51, the dictionary table 52, the compressed file table 53, the vector table 54, the inverted index table 55, the subgenome table T1, and the alternative gene vector table T2 are generated by the preprocessing unit 251. The preprocessing unit 251 may obtain the training data 90 from an external device or the preprocessing unit 251 may generate the training data 90.
In a case where the calculation unit 253 has received specification by the analysis query 92, the calculation unit 253 calculates a vector of each common structure to be subjected to a genetic modification via a synthetic pathway for the subgenome in the analysis query 92, by using the trained model 91 that has been trained. The calculation unit 253 outputs the calculated vector of each common structure, to the analysis unit 254.
In the description hereinafter, the vectors of common structures calculated by the calculation unit 253 will each be referred to as the “common structure vector”. The analysis unit 254 generates the common structure and genetically modified mechanism table 93 on the basis of the vector of the subgenome in the analysis query 92, the common structure vectors, and the gene vector table T2. An example of a process by the analysis unit 254 will be described hereinafter.
The analysis unit 254 calculates distances between a vector of a subgenome and vectors included in the alternative gene vector table T2 respectively to determine any vector having a distance less than a threshold, the distance being from the vector of the subgenome. Any vector included in the alternative gene vector table T2 and having a distance less than the threshold will be referred to as a “similar vector”, the distance being from the vector of the subgenome.
By subtracting the common structure vector from the similar vector, the analysis unit 254 calculates a vector of a genetically modified structure, and determines a correspondence relation between the common structure vector and the vector of the genetically modified structure. The analysis unit 254 registers the common structure vector and the vector of the genetically modified structure into the common structure and genetically modified structure table 93. By repeatedly executing the above described process, the analysis unit 254 generates the common structure and genetically modified structure table 93. The analysis unit 245 may output the common structure and genetically modified structure table 93 to the display unit 230 to cause the display unit 230 to display the common structure and genetically modified structure table 93, or may transmit the common structure and genetically modified structure table 93 to an external device connected to a network.
An example of a procedure by the information processing apparatus 200 according to the second embodiment will be described next.
On the basis of the subgenome table T1, the calculation unit 253 converts the subgenome in the analysis query 92 into a vector (Step S302).
By inputting the vector of the subgenome into the trained model 91 that has been trained, the calculation unit 253 calculates a vector of a common structure (Step S303). On the basis of distances between the vector of the common structure and vectors in the alternative gene vector table T2, the analysis unit 254 of the information processing apparatus 200 determines any similar vector (Step S304).
The analysis unit 254 calculates a vector of a genetically modified structure by subtracting the vector of the common structure from the vector of each gene vector similar to the subgenome (Step S305). The analysis unit 254 registers a relation between the vector of the common structure and the vector of the genetically modified structure into the common structure and genetically modified structure table 93 (Step S306). The analysis unit 254 outputs information on the common structure and genetically modified structure table (Step S307).
Effects of the information processing apparatus 200 according to the second embodiment will be described next. The information processing apparatus 200 inputs a vector of the analysis query 92, into the trained model 91 that has been trained to calculate a vector of each common structure corresponding to the subgenome in the analysis query. Furthermore, by subtraction of a vector of a common structure from the vector of each gene vector similar to the subgenome, a vector of a genetically modified structure corresponding to difference between the similar subgenome and the gene vector is calculated. Using the above described vector of the common structure and vector of the genetically modified structure facilitates analysis for a better gene vector usable for genetic modification to the target genome and resynthesis and manufacture of the target genome.
Subgenomes and gene vectors are each a secondary structure composed of plural protein primary structures. Furthermore, using dispersion vectors of the protein primary structures enables estimation of protein primary structures adjacent to a certain protein primary structure and is able to be applied to evaluation of bonding and stability of each protein primary structure. With respect to genetic modification from a gene vector to a proven subgenome, executing machine training on the basis of dispersion vectors of plural protein primary structures composing a secondary structure of the subgenome or gene vector enables improvement of: application of the gene vector; the genetic modification; and accuracy of analysis of resynthesis.
An example of a hardware configuration of a computer that implements functions that are the same as those of the above described information processing apparatus 100 (200) according to the embodiment will be described next.
As illustrated in
The hard disk device 307 has a preprocessing program 307a, a training program 307b, a calculation program 307c, and an analysis program 307d. Furthermore, the CPU 301 reads the programs 307a to 307d and load the read programs 307a to 307d into the RAM 306.
The preprocessing program 307a functions as a preprocessing process 306a. The training program 307b functions as a training process 306b. The calculation program 307c functions as a calculation process 306c. The analysis program 307d functions as an analysis process 306d.
A process by the preprocessing process 306a corresponds to the process by the preprocessing unit 151 or 251. A process by the training process 306b corresponds to the process by the training unit 152 or 252. A process by the calculation process 306c corresponds to the process by the calculation unit 153 or 253. A process by the analysis process 306d corresponds to the process by the analysis unit 154.
The programs 307a to 307d are not necessarily stored in the hard disk device 307 beforehand. For example, each program is stored in a “portable physical medium”, such as a flexible disk (FD), a CD-ROM, a DVD, a magneto-optical disk, or an IC card, which is to be inserted in the computer 300. The computer 300 may then read and execute the programs 307a to 307d therefrom.
Determination of a genome is enabled, the genome serving as a substitute for a subgenome included in a target genome.
All examples and conditional language provided herein are intended for the pedagogical purposes of aiding the reader in understanding the invention and the concepts contributed by the inventors to further the art, and are not to be construed as limitations to such specifically recited examples and conditions, nor does the organization of such examples in the specification relate to a showing of the superiority and inferiority of the invention. Although one or more embodiments of the present invention have been described in detail, it should be understood that the various changes, substitutions, and alterations could be made hereto without departing from the spirit and scope of the invention.
This application is a continuation application of International Application PCT/JP2021/015983 filed on Apr. 20, 2021 and designating U.S., the entire contents of which are incorporated herein by reference.
Number | Date | Country | |
---|---|---|---|
Parent | PCT/JP2021/015983 | Apr 2021 | US |
Child | 18468023 | US | |
Parent | PCT/JP21/15983 | Apr 2021 | US |
Child | PCT/JP2021/015983 | US |