This specification relates to a neural network for processing aptamer data.
Neural networks are machine learning models that employ one or more layers of models to generate an output, e.g., one or more classifications, for a received input. Some neural networks include one or more hidden layers in addition to an output layer. The output of each hidden layer is used as input to the next layer in the network, i.e., the next hidden layer or the output layer of the network. Each layer of the network generates an output from a received input in accordance with current values of a respective set of parameters.
In general, this specification describes a neural network for processing aptamer data.
In general, one innovative aspect of the subject matter described in this specification can be embodied in methods that include the actions of obtaining data defining a sequence for an aptamer, the aptamer comprising a string of nucleobases; encoding the data defining the sequence for the aptamer as a neural network input; and processing the neural network input using a neural network to generate an output that characterizes how strongly the aptamer binds to a particular target molecule, wherein the neural network has been configured through training to receive the data defining the sequence and to process the data to generate predicted outputs that characterize how strongly the aptamer binds to the particular target molecule.
Implementations can include one or more of the following features. The output that characterizes how strongly the aptamer binds comprises a plurality of SELEX round outputs, each SELEX round output corresponding to a respective round in an iterative SELEX process, wherein each SELEX round output characterizes how strongly the aptamer binds to the particular target molecule for the corresponding SELEX round. The SELEX round output corresponding to an earlier SELEX round in the iterative SELEX process is a binary value indicating whether the aptamer exists in a pool for use in a subsequent SELEX round. A SELEX round output for a later SELEX round in the SELEX process comprises a high throughput sequencing count of the aptamer, the high throughput sequencing count indicating how strongly the aptamer binds to the particular target molecule. The output that characterizes how strongly the aptamer binds comprises a vector of fluorescence levels, where a brightness of each fluorescence level indicates how strongly the aptamer binds to the particular target molecule. The neural network is a convolutional neural network. The data defining the sequence comprises data specifying secondary structure of the aptamer. The data specifying the secondary structure comprises data specifying a pattern of hydrogen bonds in the aptamer or an energy dot plot matrix of the aptamer. Encoding the data defining the sequence comprises converting the data defining the sequence to a plurality of one-hot vectors, each one-hot vector corresponding to a respective nucleobase in the aptamer. Iteratively applying mutations to the sequence to generate a mutated sequence that has a maximum strength of binding to the particular target molecule, comprising for each of a plurality of iterations: applying a mutation to a current mutated sequence to generate a new mutated sequence; processing the new mutated sequence using the neural network to generate an output that characterizes how strongly the new mutated sequence binds to a particular target molecule; using the neural network to determine whether the output indicates that the new mutated sequence binds stronger than the current mutated sequence; and selecting another mutation to be applied in a next iteration based on whether the output indicates that the new mutated sequence binds stronger than the current mutated sequence.
Another innovative aspect includes obtaining data defining a plurality of sequences for aptamers, each aptamer comprising a string of nucleobases, each sequence having a respective label comprising characteristics of how strongly the respective aptamer binds to a particular target molecule; encoding the data to generate training data for a neural network; and training the neural network on the training data, the neural network configured to output characteristics of how strongly a particular input sequence binds to the particular target molecule.
Implementations can include one or more of the following features. The neural network is configured to output a plurality of SELEX round outputs, each SELEX round output corresponding to a respective round in an iterative SELEX process, wherein each SELEX round output characterizes how strongly an input aptamer binds to the particular target molecule for the corresponding SELEX round. Each label for a respective aptamer in the training data further comprises, a plurality of binary values for the respective aptamer for a respective plurality of SELEX rounds, each binary value for the respective SELEX round indicating whether the aptamer exists in a pool for use in a subsequent SELEX round, and wherein the neural network output further comprises a plurality of binary values for a respective plurality of SELEX rounds, each binary value for the respective SELEX round indicating whether the aptamer exists in a pool for use in a subsequent SELEX round. Each label for a respective aptamer in the training data further comprises, a plurality of high throughput sequencing counts of the respective aptamer for a respective plurality of SELEX rounds, and wherein the neural network output further comprises a plurality of high throughput sequencing counts of the aptamer, each high throughput sequencing count indicating how strongly the aptamer binds to the particular target molecule. One or more labels for a respective aptamer in the training data further comprises, a vector of fluorescence levels for the respective aptamer, and wherein the neural network output further comprises a vector of fluorescence levels, where a brightness of each fluorescence level indicates how strongly the aptamer binds to the particular target molecule. The data defining the plurality of sequences for aptamers comprises, for each aptamer, data specifying secondary structure of the aptamer. The data specifying the secondary structure comprises data specifying a pattern of hydrogen bonds in the aptamer or an energy dot plot matrix of the aptamer.
Particular embodiments of the subject matter described in this specification can be implemented so as to realize one or more of the following advantages. Conventional methods of identifying binding strength between aptamers and a particular target molecule include applying a SELEX process to randomly selected aptamers, which can be labor-intensive and frequently fails to yield an aptamer with strong specific binding. Advantageously, a neural network can predict a binding strength for an input aptamer to a particular target molecule. In particular, the neural network can predict the SELEX round outputs for novel sequences not found in the original random pool. Therefore, instead of requiring the best aptamer to be present in the initial, relatively small, random pool, scientists can predict the aptamer having the strongest specific binding in the set of all possible sequences. Thus, rather than performing multiple SELEX rounds that may or may not yield an aptamer with strong specific binding, the neural network predicts aptamers with strong specific bindings, which reduces cost and testing time. The neural network can also model binding to multiple targets at once, and therefore identify an aptamer that likely binds to the target molecule and not bind to known background molecules. The details of one or more embodiments of the subject matter of this specification are set forth in the accompanying drawings and the description below. Other features, aspects, and advantages of the subject matter will become apparent from the description, the drawings, and the claims.
Like reference numbers and designations in the various drawings indicate like elements.
The system 100 includes a neural network 104. In some implementations, the neural network 104 is a convolutional neural network. The neural network 104 receives aptamer data 102. The aptamer data 102 can be data defining a sequence for an aptamer. An aptamer is a single-stranded molecule that can bind to particular target molecules including proteins and peptides with a degree of affinity, i.e., strength.
The aptamer data can include a sequence of nucleotide bases for an aptamer. In some implementations, the sequence of nucleotide bases is encoded as a sequence of one-hot vectors. For example, a cytosine molecule in the sequence can be represented by a vector [1 0 0 0] while a thymine molecule can be represented by a vector [0 1 0 0].
The vectors can be appended to each other to form a matrix of values, e.g., represented as an image, which can be provided as a tensor input to the neural network 104. For example, if the sequence of nucleotide bases is CTCT, the matrix of values can be the one-hot vectors of the respective nucleotide bases appended together to form the matrix:
[1 0 1 0 0 1 0 1 0 0 0 0 0 0 0 0].
By way of illustration, in this example the matrix can be represented as a 4×4 image, where the 1's in the matrix are represented as black pixels and the 0's in the matrix are represented as white pixels. The image representing the matrix can be provided as input to the neural network 104. In some implementations, the neural network 104 performs convolutions on the 4×4 image matrix in a horizontal direction and not a vertical direction. This enables predictions to be based on which nucleotide bases are sequentially adjacent to each other and not based on how one-hot encodings of the nucleotides are implemented. In some other implementations, the matrix is provided as a sequence of bits to the neural network.
The aptamer data can also include data characterizing the secondary structure of the aptamer. For example, the data can specify a pattern of hydrogen bonds in the aptamer or an energy dot plot matrix of the aptamer. Characteristics of the secondary structure can be represented as a byte sequence, e.g., as a one-hot vector or a matrix of values, and provided to the neural network 104. In some implementations, the byte sequence of the secondary structure is appended to the matrix of values formed from the sequence of nucleotide bases to generate the aptamer data that is provided as input to the neural network 104.
The neural network 104 can process the aptamer data 102 to generate an output 106. The output 106 characterizes how strongly the aptamer binds to a particular target molecule. In particular, the output 106 can correspond to outputs generated from a SELEX process for the aptamer data 102.
The SELEX process, i.e., systematic evolution of ligands by exponential enrichment (SELEX), identifies binding strengths to a particular target molecule from a pool of aptamers. In particular, the SELEX process starts with an initial pool of aptamers and processes the pool of aptamers through multiple SELEX rounds. Each SELEX round selects a pool of aptamers that bind the target molecule and discards aptamers that do not bind to the target molecule. After each SELEX round, the resulting pool is amplified using polymerase chain reaction (PCR). During PCR, aptamers with stronger binding strengths tend to be replicated more frequently than aptamers with weaker binding strengths.
After numerous SELEX rounds, a composition of the pool of aptamers indicates a convergence to a small number, e.g., on the order of thousands, of aptamers with high counts, which can potentially indicate a strong binding strength to a target molecule.
By way of illustration, for a SELEX process, the initial pool can have a size on the order of 1013 aptamers. The initial pool can be referred to as the library of aptamers. The pool can be composed of nearly entirely unique aptamers. After a first SELEX round, i.e., after discarding aptamers that do not bind to a target molecule and enriching aptamers that do bind, the pool can be 99.3% unique. After a second SELEX round, the pool can be 90% unique and have a size on the order of 108 aptamers. After another SELEX round, the pool can be less than 10% unique and have a size on the order of 106 aptamers. After yet another SELEX round, the pool can be less than 5% unique and have a size on the order of 105 aptamers. After each SELEX round, the pool can undergo high throughput sequencing to identify aptamers in the SELEX round having high counts.
As described above, the output 106 corresponds to the outputs from the SELEX process for the aptamer data 102. In other words, the neural network 104 predicts, after each SELEX round, how strongly the input aptamer in the aptamer data 102 binds to the target molecule.
Thus, the output 106 can include a library output 108. In some implementations, the library output 108 is a binary value, e.g., true or false, of whether the input aptamer exists in the initial pool of aptamers. In some other implementations, the library output 108 is a confidence value between 0 and 1 that represents a likelihood that the input aptamer exists in the initial pool of aptamers.
The output 106 can include SELEX 1 output 110. Selex 1 output 110 is the output after performing a first SELEX round in the SELEX process. Selex 1 output 110 can also be a confidence value or a binary value of whether the input aptamer exists in the pool of aptamers after the first round. In some implementations, the neural network 104 generates binary or confidence values for SELEX outputs that are in pools with a percentage of unique aptamers over a particular threshold. This is because a pool of aptamers with a high unique aptamer percentage generally does not have high aptamer count for one particular aptamer, and thus does not indicate any particular aptamer with a high binding strength to the target molecule.
The output 106 can include a SELEX 2 output and a SELEX 3 output 112, 114. Selex 2 and 3 output 112, 114 are outputs after performing a second and third SELEX round, respectively. The output for a respective one of these SELEX rounds can be a predicted count of a number of input aptamers in the respective SELEX round.
In some implementations, the neural network 104 generates, in the output 106, an array or vector of fluorescent values 116. The fluorescent values represent a brightness, which indicates how strongly the input aptamer binds to the target molecule as measured in a binding array.
To generate the output 106, the system 100 can train the neural network 101 using conventional neural network training techniques, which will be described further below with reference to
The system obtains data defining multiple sequences for aptamers (step 202). Each sequence can be represented as multiple one-hot vectors, as described above with reference to
Each sequence can also have a respective label that includes characteristics of how strongly the respective aptamer binds to a particular target molecule. A label can have the data specified by the output 106 of
In particular, the label can include binary values for the respective aptamer for earlier SELEX rounds. The binary values can indicate whether the aptamer exists in a pool of aptamers for the round, and thus, whether the aptamer exists in a pool for use in a subsequent SELEX round. The label can also include high throughput sequencing counts of the respective aptamer for later SELEX rounds. In some implementations, the labels also include data of vectors of fluorescence levels for the aptamer. These vectors of fluorescence levels can be obtained through binding array measurements. In some other implementations, the labels include vectors of binding affinity values, where each binding affinity value corresponds to a likelihood of binding to a respective target molecule.
The system encodes the data to generate training data for a neural network (step 204). The system can encode the data defining the sequence by converting the data to multiple one-hot vectors, as described above with reference to
The system trains the neural network on the training data (step 206). That is, the system processes the training data, e.g., the encoded sequence data, to generate predicted outputs. The system determines an error between each predicted output and a corresponding portion of the label. The system uses the error to update values of the parameters used in each layer of the neural network using conventional neural network training techniques, e.g., stochastic gradient descent with backpropagation. After training, the neural network can generate outputs, e.g., the output 106 of
The system obtains data defining a sequence for an aptamer (step 302). Each sequence can be represented as multiple one-hot vectors, as described above with reference to
The system encodes the data defining the sequence as a neural network input (step 304). The system can encode the data defining the sequence by converting the data to multiple one-hot vectors, as described above with reference to
The system processes the neural network input using a neural network to generate an output that characterizes how strongly the aptamer binds to a particular target molecule (step 306). The neural network can be trained, as described above with reference to
Each SELEX round output can correspond to a respective round in an iterative SELEX process, as described above with reference to
In some implementations, the output also includes an array output of a vector of fluorescent values.
Although the output, e.g., the output 106 of
For each iteration, the system applies a mutation to a current mutated sequence, i.e., the sequence of the input aptamer in the first iteration, to generate a new mutated sequence (step 402). The mutation can randomly change a number of nucleotide bases in the sequence. For example, the mutation can randomly select two nucleotide bases and change them to two other nucleotide bases.
The system can process the new mutated sequence using the neural network to generate another output that characterizes how strongly the new mutated sequence binds to a particular target molecule (step 404), e.g., as described above with reference to
The system uses the neural network to determine whether the output indicates that the new mutated sequence binds stronger than the current mutated sequence (step 406). That is, the system can compare a high throughput sequencing count for a particular SELEX round of the current mutated sequence with that of the new mutated sequence. If the high throughput sequencing count for the new mutated sequence is higher than that of the current mutated sequence, the system determines the new mutated sequence binds more strongly to the target molecule than the current mutated sequence. Otherwise, the system determines the current mutated sequence binds more strongly to the target molecule than the new mutated sequence. In some implementations, the system factors in high throughput sequencing counts between the current mutated sequence and the new mutated sequence across multiple SELEX rounds. By way of illustration, the system can compute weighted sums of high throughput sequencing counts across multiple SELEX rounds for each sequence to determine which sequence more strongly binds to the target molecule.
The system selects another mutation to be applied in a next iteration based on whether the output indicates that the new mutated sequence binds stronger than the current mutated sequence (step 408). That is, if the new mutated sequence binds more strongly than the current mutated sequence, the system randomly mutates the new mutated sequence and repeats the process described above. Otherwise, the system randomly mutates the current mutated sequence and repeats the process described above. In some implementations, the system mutates the current mutated sequence to a sequence not previously processed. The system can continue mutating sequences until the system generates a threshold number of new mutated sequences or if a higher binding strength is not found for a threshold number of iterations.
The system can repeat the process for a predetermined number of times. Alternatively, the system can repeat the process until the system iterates on the current mutated sequence for a threshold number of times. In other words, the system stops repeating the process when the system is unable to generate, within the threshold number of times, a new mutated sequence that binds more strongly than the current mutated sequence.
Embodiments of the subject matter and the functional operations described in this specification can be implemented in digital electronic circuitry, in tangibly-embodied computer software or firmware, in computer hardware, including the structures disclosed in this specification and their structural equivalents, or in combinations of one or more of them. Embodiments of the subject matter described in this specification can be implemented as one or more computer programs, i.e., one or more modules of computer program instructions encoded on a tangible non-transitory program carrier for execution by, or to control the operation of, data processing apparatus. Alternatively or in addition, the program instructions can be encoded on an artificially-generated propagated signal, e.g., a machine-generated electrical, optical, or electromagnetic signal, that is generated to encode information for transmission to suitable receiver apparatus for execution by a data processing apparatus. The computer storage medium can be a machine-readable storage device, a machine-readable storage substrate, a random or serial access memory device, or a combination of one or more of them. The computer storage medium is not, however, a propagated signal.
The term “data processing apparatus” encompasses all kinds of apparatus, devices, and machines for processing data, including by way of example a programmable processor, a computer, or multiple processors or computers. The apparatus can include special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application-specific integrated circuit). The apparatus can also include, in addition to hardware, code that creates an execution environment for the computer program in question, e.g., code that constitutes processor firmware, a protocol stack, a database management system, an operating system, or a combination of one or more of them.
A computer program (which may also be referred to or described as a program, software, a software application, a module, a software module, a script, or code) can be written in any form of programming language, including compiled or interpreted languages, or declarative or procedural languages, and it can be deployed in any form, including as a stand-alone program or as a module, component, subroutine, or other unit suitable for use in a computing environment. A computer program may, but need not, correspond to a file in a file system. A program can be stored in a portion of a file that holds other programs or data, e.g., one or more scripts stored in a markup language document, in a single file dedicated to the program in question, or in multiple coordinated files, e.g., files that store one or more modules, sub-programs, or portions of code. A computer program can be deployed to be executed on one computer or on multiple computers that are located at one site or distributed across multiple sites and interconnected by a communication network.
As used in this specification, an “engine,” or “software engine,” refers to a software implemented input/output system that provides an output that is different from the input. An engine can be an encoded block of functionality, such as a library, a platform, a software development kit (“SDK”), or an object. Each engine can be implemented on any appropriate type of computing device, e.g., servers, mobile phones, tablet computers, notebook computers, music players, e-book readers, laptop or desktop computers, PDAs, smart phones, or other stationary or portable devices, that includes one or more processors and computer readable media. Additionally, two or more of the engines may be implemented on the same computing device, or on different computing devices.
The processes and logic flows described in this specification can be performed by one or more programmable computers executing one or more computer programs to perform functions by operating on input data and generating output. The processes and logic flows can also be performed by, and apparatus can also be implemented as, special purpose logic circuitry, e.g., an FPGA (field programmable gate array) or an ASIC (application-specific integrated circuit).
Computers suitable for the execution of a computer program include, by way of example, can be based on general or special purpose microprocessors or both, or any other kind of central processing unit. Generally, a central processing unit will receive instructions and data from a read-only memory or a random access memory or both. The essential elements of a computer are a central processing unit for performing or executing instructions and one or more memory devices for storing instructions and data. Generally, a computer will also include, or be operatively coupled to receive data from or transfer data to, or both, one or more mass storage devices for storing data, e.g., magnetic, magneto-optical disks, or optical disks. However, a computer need not have such devices. Moreover, a computer can be embedded in another device, e.g., a mobile telephone, a personal digital assistant (PDA), a mobile audio or video player, a game console, a Global Positioning System (GPS) receiver, or a portable storage device, e.g., a universal serial bus (USB) flash drive, to name just a few.
Computer-readable media suitable for storing computer program instructions and data include all forms of non-volatile memory, media and memory devices, including by way of example semiconductor memory devices, e.g., EPROM, EEPROM, and flash memory devices; magnetic disks, e.g., internal hard disks or removable disks; magneto-optical disks; and CD-ROM and DVD-ROM disks. The processor and the memory can be supplemented by, or incorporated in, special purpose logic circuitry.
To provide for interaction with a user, embodiments of the subject matter described in this specification can be implemented on a computer having a display device, e.g., a CRT (cathode ray tube) monitor, an LCD (liquid crystal display) monitor, or an OLED display, for displaying information to the user, as well as input devices for providing input to the computer, e.g., a keyboard, a mouse, or a presence sensitive display or other surface. Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback, e.g., visual feedback, auditory feedback, or tactile feedback; and input from the user can be received in any form, including acoustic, speech, or tactile input. In addition, a computer can interact with a user by sending resources to and receiving resources from a device that is used by the user; for example, by sending web pages to a web browser on a user's client device in response to requests received from the web browser.
Embodiments of the subject matter described in this specification can be implemented in a computing system that includes a back end component, e.g., as a data server, or that includes a middleware component, e.g., an application server, or that includes a front end component, e.g., a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the subject matter described in this specification, or any combination of one or more such back end, middleware, or front end components. The components of the system can be interconnected by any form or medium of digital data communication, e.g., a communication network. Examples of communication networks include a local area network (“LAN”) and a wide area network (“WAN”), e.g., the Internet.
The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
While this specification contains many specific implementation details, these should not be construed as limitations on the scope of any invention or of what may be claimed, but rather as descriptions of features that may be specific to particular embodiments of particular inventions. Certain features that are described in this specification in the context of separate embodiments can also be implemented in combination in a single embodiment. Conversely, various features that are described in the context of a single embodiment can also be implemented in multiple embodiments separately or in any suitable subcombination.
Moreover, although features may be described above as acting in certain combinations and even initially claimed as such, one or more features from a claimed combination can in some cases be excised from the combination, and the claimed combination may be directed to a subcombination or variation of a subcombination.
Similarly, while operations are depicted in the drawings in a particular order, this should not be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results. In certain circumstances, multitasking and parallel processing may be advantageous. Moreover, the separation of various system modules and components in the embodiments described above should not be understood as requiring such separation in all embodiments, and it should be understood that the described program components and systems can generally be integrated together in a single software product or packaged into multiple software products.
Particular embodiments of the subject matter have been described. Other embodiments are within the scope of the following claims. For example, the actions recited in the claims can be performed in a different order and still achieve desirable results. As one example, the processes depicted in the accompanying figures do not necessarily require the particular order shown, or sequential order, to achieve desirable results. In certain implementations, multitasking and parallel processing may be advantageous.
This is a continuation of U.S. application Ser. No. 16/752,513, filed on Jan. 24, 2020, which is a continuation of U.S. application Ser. No. 14/921,973, filed on Oct. 23, 2015 (now U.S. Pat. No. 10,546,650). The disclosures of the prior applications are considered part of and are incorporated by reference in the disclosure of this application.
Number | Date | Country | |
---|---|---|---|
Parent | 16752513 | Jan 2020 | US |
Child | 18307748 | US | |
Parent | 14921973 | Oct 2015 | US |
Child | 16752513 | US |