WO2018138727A1 - Séquences d'acide nucléique synthétiques virales et leur utilisation - Google Patents

Séquences d'acide nucléique synthétiques virales et leur utilisation Download PDF

Info

Publication number
WO2018138727A1
WO2018138727A1 PCT/IL2018/050093 IL2018050093W WO2018138727A1 WO 2018138727 A1 WO2018138727 A1 WO 2018138727A1 IL 2018050093 W IL2018050093 W IL 2018050093W WO 2018138727 A1 WO2018138727 A1 WO 2018138727A1
Authority
WO
WIPO (PCT)
Prior art keywords
codon
genes
viral
virus
nucleic acid
Prior art date
Application number
PCT/IL2018/050093
Other languages
English (en)
Inventor
Tamir Tuller
Eli Goz
Oriah MIODUSER
Alon DIAMENT
Original Assignee
Synvaccine Ltd.
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Synvaccine Ltd. filed Critical Synvaccine Ltd.
Priority to US16/480,671 priority Critical patent/US20190359990A1/en
Priority to EP18745104.2A priority patent/EP3574097A4/fr
Publication of WO2018138727A1 publication Critical patent/WO2018138727A1/fr

Links

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09—Recombinant DNA-technology
    • C12N15/63—Introduction of foreign genetic material using vectors; Vectors; Use of hosts therefor; Regulation of expression
    • C12N15/67—General methods for enhancing the expression
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00—Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09—Recombinant DNA-technology
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N7/00—Viruses; Bacteriophages; Compositions thereof; Preparation or purification thereof
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
    • G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2710/00—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA dsDNA viruses
    • C12N2710/00011—Details
    • C12N2710/16011—Herpesviridae
    • C12N2710/16611—Simplexvirus, e.g. human herpesvirus 1, 2
    • C12N2710/16621—Viruses as such, e.g. new isolates, mutants or their genomic sequences
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2710/00—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA dsDNA viruses
    • C12N2710/00011—Details
    • C12N2710/16011—Herpesviridae
    • C12N2710/16611—Simplexvirus, e.g. human herpesvirus 1, 2
    • C12N2710/16634—Use of virus or viral component as vaccine, e.g. live-attenuated or inactivated virus, VLP, viral protein
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2710/00—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA dsDNA viruses
    • C12N2710/00011—Details
    • C12N2710/16011—Herpesviridae
    • C12N2710/16611—Simplexvirus, e.g. human herpesvirus 1, 2
    • C12N2710/16661—Methods of inactivation or attenuation
    • C12N2710/16662—Methods of inactivation or attenuation by genetic engineering
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2770/00—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA ssRNA viruses positive-sense
    • C12N2770/00011—Details
    • C12N2770/24011—Flaviviridae
    • C12N2770/24111—Flavivirus, e.g. yellow fever virus, dengue, JEV
    • C12N2770/24121—Viruses as such, e.g. new isolates, mutants or their genomic sequences
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2770/00—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA ssRNA viruses positive-sense
    • C12N2770/00011—Details
    • C12N2770/24011—Flaviviridae
    • C12N2770/24111—Flavivirus, e.g. yellow fever virus, dengue, JEV
    • C12N2770/24134—Use of virus or viral component as vaccine, e.g. live-attenuated or inactivated virus, VLP, viral protein
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2770/00—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA ssRNA viruses positive-sense
    • C12N2770/00011—Details
    • C12N2770/24011—Flaviviridae
    • C12N2770/24111—Flavivirus, e.g. yellow fever virus, dengue, JEV
    • C12N2770/24161—Methods of inactivation or attenuation
    • C12N2770/24162—Methods of inactivation or attenuation by genetic engineering
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2795/00—Bacteriophages
    • C12N2795/00011—Details
    • C12N2795/10011—Details dsDNA Bacteriophages
    • C12N2795/10311—Siphoviridae
    • C12N2795/10321—Viruses as such, e.g. new isolates, mutants or their genomic sequences
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2795/00—Bacteriophages
    • C12N2795/00011—Details
    • C12N2795/10011—Details dsDNA Bacteriophages
    • C12N2795/10311—Siphoviridae
    • C12N2795/10361—Methods of inactivation or attenuation
    • C12N2795/10362—Methods of inactivation or attenuation by genetic engineering
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2800/00—Nucleic acids vectors
    • C12N2800/22—Vectors comprising a coding region that has been codon optimised for expression in a respective host
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N2840/00—Vectors comprising a special translation-regulating system
    • C12N2840/10—Vectors comprising a special translation-regulating system regulates levels of translation

Definitions

  • the present invention is directed to the field of viral genome optimization.
  • ribosome profiling which enables in vivo genome wide monitoring of ribosome states at a resolution of a single nucleotide, is a useful tool for deciphering the coding complexity of viral (and other organisms) genomes. Specifically, it was shown that ribosome profiling enables detecting novel (possibly very short) coding regions and estimating the translation status of various open reading frames.
  • Gene expression within different DNA viruses is regulated in a temporal fashion and can be divided into early and late stages with respect to the viral replication cycle.
  • the early genes are expressed following the entry into the host cell and code typically for non-structural proteins that are responsible for different regulatory functions in processes such as: viral DNA replication, activation of late genes expression, trans-nuclear transport, interaction with the host cell, induction of the cell's DNA replication machinery necessary for viral replication, etc.
  • Late genes largely code for structural proteins required for virion assembly; they are generally highly expressed and their expression is usually induced or regulated by early genes.
  • Vaccines and their discovery, are topics of singular importance in present-day biomedical science; however, the discovery of vaccines has hitherto been primarily empirical in nature requiring considerable investments of time, efforts and resourced. To overcome the numerous pitfalls attributed to the classical vaccine design strategies, more efficient and robust rational approaches based on computer-based methods are highly desirable.
  • One direction in designing in- silico vaccine candidates may be based on exploiting the synonymous information encoded in the genomes for attenuating the viral replication cycle while retaining the wild type proteins.
  • the present invention provides a nucleic acid molecule comprising a codon substituted to a synonymous codon.
  • the invention also provides modified viruses and vaccine compositions comprising the nucleic acid molecules of the invention, as well as methods for producing and using these viruses and vaccine compositions.
  • a nucleic acid molecule comprising a coding sequence, the coding sequence comprises at least one codon substituted to a synonymous codon, the synonymous codon has a parameter that effects translation rate which differs from the parameter of the at least one codon in a set of coordinately expressed viral genes.
  • the parameter than effects translation rate is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation, codon typical decoding rate (TDR), GC content, average repetitive substring index, codon pair bias, dinucleotide bias, nucleotide bias and amino acid bias.
  • the parameter than effects translation rate is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation and codon typical decoding rate (TDR).
  • the set of coordinately expressed viral genes are a subset of a virus's total genome.
  • the set of viral genes is selected from the group consisting of: genes with a common structure, genes with common temporal expression during a virus's life cycle, genes encoding for proteins having common function, genes translated at a common cellular location and genes encoding for proteins having common cellular localization.
  • At least ten codons of the coding sequence have been substituted.
  • at least 10% of the codons of the coding sequence have been substituted.
  • the synonymous codon has a parameter that results in faster translation in a particular cellular context relative to the at least one codon that has been substituted.
  • a plurality of codons having synonymous codons with faster translation rates in a particular cellular context have been substituted.
  • the at least one codon is substituted to the fastest translating synonymous codon in a particular cellular context.
  • the coding sequence is optimized for expression that is the same as the set of viral genes.
  • the synonymous codon has a parameter that results in slower translation in a particular cellular context relative to the at least one codon that has been substituted.
  • a plurality of codons having synonymous codons with slower translation rates in a particular cellular context have been substituted.
  • the at least one codon is substituted to the slowest translating synonymous codon in a particular cellular context.
  • the coding sequence is de-optimized for expression that is the same as said set of viral genes.
  • the coding sequence is selected from the group consisting of: a viral gene, a mammalian gene, a heterologous transgene or a fragment thereof.
  • the mammalian gene is a human gene.
  • the temporal expression during a virus's life cycle is selected from the group consisting of: early expression in the virus life cycle, intermediate expressing in the virus life cycle, and late expression in the virus life cycle.
  • the optimization is selected from optimizing translation efficiency and optimizing protein yield.
  • a modified virus comprising the nucleic acid molecule of the present invention.
  • an attenuated form of a virus comprising the nucleic acid molecule of the present invention, wherein:
  • the coding sequence is an essential viral sequence endogenously expressed during a defined period of the virus's life cycle
  • the set of genes is a plurality of viral genes expressed during the defined period of the virus's life cycle
  • virus is devoid of the endogenous viral sequence.
  • a vaccine composition for inducing a protective immune response in a subject comprising:
  • a method for vaccinating a subject at risk of viral infection comprising, administering to the subject the vaccine composition of the present invention.
  • nucleic acid molecule optimized for expression in a particular cellular context comprising:
  • a method for producing a nucleic acid molecule deoptimized for expression in a particular cellular context comprising: a. selecting a coding sequence
  • nucleic acid molecule optimized for expression in a particular cellular context thereby producing a nucleic acid molecule optimized for expression in a particular cellular context.
  • the parameter than effects translation is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation, codon typical decoding rate (TOR), GC content, average repetitive substring index, codon pair bias, dinucleotide bias, nucleotide bias and amino acid bias.
  • the parameter than effects translation is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation and codon typical decoding rate (TOR).
  • the particular cellular context is selected from the group consisting of: a specific period of a cell's life cycle, during cellular division, during cellular stress, during apoptosis, during viral infection, or during viral lysogeny of the cell.
  • a method for producing a modified virus comprising:
  • the viral sequence is an essential viral sequence
  • FIG. 1A A flow diagram and illustration of some embodiments of the study disclosed herein.
  • IB A bar chart depicting relative expression levels of each of the gene groups (early late) in read count per nucleotide.
  • FIGS 2A-F Figures 2A-F.
  • (2A) A plot of two principal components of codon relative frequency vectors for all bacteriophage genes and the top 50 E. Coli genes, indicating that early and late genes tend to be clustered into two distinct groups.
  • (2B) Histograms of 500 random average MTDR values of early and late viral genes at different time points; the wild type average MTDR appears as a line.
  • (2C) The wild type average MTDR appears as a line.
  • RTEC (mean MTDR E - mean MTDR L )/(mean MTDR E + mean MTDR L ) as a function of time from the beginning of the lytic stage (0-20 minutes), where MTDR E and MTDR L signify the MTDR of early and late genes respectively.
  • 2D A line graph showing correlation between codon typical decoding rates (TDRs) and codon frequencies at different time points for E. coli genes, just early viral genes, just late viral genes and all viral genes. Time points with significant correlations (Spearman p- values lower than 0.05) are marked by asterisk.
  • TDRs codon typical decoding rates
  • 2E The wild type mean tRNA adaptation index (tAI) (line) vs.
  • FIG. 3A Histograms of 500 random average MTDR values of all viral genes in the grey; the wild type average MTDR appears as a line.
  • FIG. 4A-D Figures 4A-D.
  • A-B Line graph profiles of folding energy (average MFE in all windows of 39-nt length) across the bacteriophage genome (gray) vs. an averaged profile corresponding to 100 randomized variants (black) based on dinucleotide preserving randomization; the window index denotes the distance (in nucleotides) from the beginning of the ORF to the beginning of the window. Regions where the folding energy of the wild type genome is significantly higher (star) or lower (dot) than in randomized variants are marked at the bottom of the figure.
  • the profiles include the 50-UTR near the beginning of the ORF (negative window indexes).
  • the profiles include the 30-UTR near the ending of the ORF (positive window indexes).
  • C Histograms of mean local folding energies (folding energies averaged over all the windows of each gene) compared with randomized mean local folding energies obtained from two models: (i) proteindinucleotides preserving and (ii) proteincodon usage bias preserving.
  • D Histograms of log[ARS index]. Eight analyses were performed: two types of reference genomes; bacterial and viral, two type of randomizations; dinucleotide and codons, two groups of genes; early and late. In each histogram, the wild type distribution is compared with the mean random distribution ( 1 ,000 random genomes).
  • FIG. 5A-H A phylogenetic tree built from complete phage proteomes using ARS distance (see Materials and Methods). Phages with significant differences in temporary codon usage are marked by blue.
  • B Table of viruses with significant (p- value ⁇ 0.05) separation between early and late genes with respect to synonymous codons or AA are marked by stars.
  • C-D Scatter plots of the CUB principal component analysis (PGA) for (C) significantly separated bacteriophages and (D) not significantly separated bacteriophages.
  • E Bar charts of PCA of codon usage variances distribution for all 14 bacteriophages.
  • Figures 6A-B Bar chart comparisons of a mean pairwise distance between the group of phages with a differential codon usage in temporal genes (light gray line) to 100 groups of randomly sampled phages of the same size (dark grey) for (A) all 14 phages and (B) all phages excluding one streptococcus.
  • FIG. Scatter plots of the CUB principal component analysis (PCA) for human viruses.
  • FIG. 8A-H Figures 8A-H.
  • A Table showing significance of separation between early and late genes with respect to additional genomic features estimated by Wilcoxon ranksum p-value. Features/viruses with significant (p-value ⁇ 0.05) separation between the two temporal groups are marked by yellow stars; green is related to higher mean in the case of the early genes and red is related to higher mean in the case of the late genes.
  • B-H Bar chart comparisons of (B) genome length, (C) ENC, (D) CPB, (E) DNTB, (F) NTB, (G) GC content and (H) AAB in viruses with significant temporal separation in codon usage bias (blue) vs. viruses with no such separation (grey). Y axis - probability; X - genomic mean feature values (each virus is represented by an average of feature values over all its genes). Wilcoxon ranksum p-values are specified.
  • D(x,y) Euclidian distance between synonymous codons relative frequency vectors of sequences x and y; CAI(x,R)/ RCDI(x,R) - CA RCDI of sequence x with respect to the reference set R; s/t - wUd-type/deoptimized UL2 coding sequences.
  • the distribution of D/CAI/RCDI values for are given for variants randomized with respect to the average viral CUB and the average human CUB. P-value and Z-score estimations with respect to each random model is specified.
  • FIGS 11A-B Zika virus modification.
  • A A table summarizing the 6 Zika variants, including the algorithm used for codon substitution, the number of codons changed, the predicted virulence, the foci relative size and representative micrographs of the virally-infected foci.
  • B A bar chart showing viral titers at various time points after infection.
  • the present invention provides, in some embodiments, nucleic acid molecules comprising a coding sequence with at least one codon substituted to a synonymous codon, a modified form of a virus comprising the nucleic acid molecules of the invention, and methods for producing these nucleic acid molecules, and viruses.
  • the present invention concerns a nucleic acid molecule comprising a coding sequence, the coding sequence comprises at least one codon substituted to a synonymous codon, and the synonymous codon has a parameter that effects translation rate which differs from the parameter of the at least one codon in a set of viral genes.
  • nucleic acid is well known in the art.
  • a “nucleic acid” as used herein will generally refer to a molecule (i.e., a strand) of DNA, RNA or a derivative or analog thereof, comprising a nucleobase.
  • a nucleobase includes, for example, a naturally occurring purine or pyrimidine base found in DNA (e.g., an adenine "A,” a guanine “G,” a thymine “T” or a cytosine “C”) or RNA (e.g., an A, a G, an uracil “U” or a C).
  • DNA e.g., an adenine "A,” a guanine "G,” a thymine “T” or a cytosine "C”
  • RNA e.g., an A, a G, an uracil "U” or a C.
  • nucleic acid molecule include but not limited to single- stranded RNA (ssRNA), double-stranded RNA (dsRNA), single-stranded DNA (ssDNA), double-stranded DNA (dsDNA), small RNA such as miRNA, siRNA and other short interfering nucleic acids, snoRNAs, snRNAs, tRNA, piRNA, tnRNA, small rRNA, hnRNA, circulating nucleic acids, fragments of genomic DNA or RNA, degraded nucleic acids, ribozymes, viral RNA or DNA, nucleic acids of infectios origin, amplification products, modified nucleic acids, plasmidical or organellar nucleic acids and artificial nucleic acids such as oligonucleotides.
  • ssRNA single- stranded RNA
  • dsRNA double-stranded RNA
  • ssDNA single-stranded DNA
  • dsDNA double-stranded DNA
  • Coding sequence refers to a nucleic acid sequence that when translated results in an expressed protein.
  • the coding sequence is to be used as a basis for making codon alterations.
  • the coding sequence is a viral gene.
  • the coding sequence is a mammalian gene.
  • the coding sequence is a human gene.
  • the coding sequence is a portion of one of the above listed genes.
  • the coding sequence is a heterologous transgene.
  • the above listed genes are wild type, endogenously expressed genes.
  • the above listed genes have been genetically modified or in some way altered from their endogenous formulation. These alterations may be changes to the coding region such that the protein the gene codes for is altered.
  • heterologous transgene refers to a gene that originated in one species and is being expressed in another. In some embodiments, the transgene is a part of a gene originating in another organism.
  • the coding sequence is selected from the group consisting of: a viral gene, a mammalian gene, a heterologous transgene or a fragment thereof.
  • codon refers to a sequence of three DNA or RNA nucleotides that correspond to a specific amino acid or stop signal during protein synthesis.
  • the codon code is degenerate, in that more than one codon can code for the same amino acid.
  • Such codons that code for the same amino acid are known as "synonymous" codons.
  • CUU, CUC, CUA, CUG, UUA, and UUG are synonymous codons that code for Leucine.
  • Synonymous codons are not used with equal frequency. In general, the most frequently used codons in a particular cell are those for which the cognate tRNA is abundant, and the use of these codons enhances the rate and/or accuracy of protein translation.
  • Codon bias refers generally to the non-equal usage of the various synonymous codons, and specifically to the relative frequency at which a given synonymous codon is used in a defined sequence or set of sequences.
  • At least 5, at least 10, at least IS, at least 20, at least 25, at least 30, at least 35, at least 40, at least 45, or at least 50 codons of the coding sequence have been substituted.
  • Each possibility represents a separate embodiment of the present invention.
  • One skilled in the art will be able to determine based on the virus and the host cells the minimum number of codons to be substituted.
  • protein expression after substitution can be measured and compared to protein expression without substitutions to determine if a sufficient number of codons have been substituted.
  • greater than 5%, greater than 10%, greater than 15%, greater than 20%, greater than 25%, greater than 30%, greater than 35%, greater than 40%, greater than 45%, greater than 50%, greater than 55%, greater than 60%, greater than 65%, greater than 70%, greater than 75%, greater than 80%, greater than 85%, greater than 90%, greater than 95%, or 100% of all codons in the target sequence have been substituted.
  • Each possibility represents a separate embodiment of the present invention.
  • Each possibility represents a separate embodiment of the present invention.
  • OSS greater than 5%, greater than 10%, greater than 15%, greater than 20%, greater than 25%, greater than 30%, greater than 35%, greater than 40%, greater than 45%, greater than 50%, greater than 55%, greater than 60%, greater than 65%, greater than 70%, greater than 75%, greater than 80%, greater than 85%, greater than 90%, greater than 95%, or 100% of codons that have synonymous codons with lower frequencies have been substituted.
  • a plurality of codons having synonymous codons with different frequencies have been substituted.
  • a plurality of codons having synonymous codons with higher frequencies have been substituted.
  • a plurality of codons having synonymous codons with lower frequencies have been substituted.
  • the term "parameter that effects translation rate” as used herein refers to any characteristic of a codon that impacts the rate at which a codon is translated in a particular and/or pre-determined cellular milieu. Not all codons translate at the same rate, and the rate of codon translation of an entire coding region is proportional to the amount of protein that can be produced by that coding region. Further a single codon may translate at a different rate depending upon the cellular context. As one non-limiting example, the tRNA pool found in a cell can change depending on the condition of the cell. Greater concentration of a particular tRNA may result in faster translation of the codon to which that tRNA corresponds.
  • the parameter than effects translation rate is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation, codon typical decoding rate (TDR), GC content, average repetitive substring index, codon pair bias, dinucleotide bias, nucleotide bias and amino acid bias.
  • the parameter than effects translation rate is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation, codon typical decoding rate (TDR), GC content, average repetitive substring index, codon pair bias, dinucleotide bias, nucleotide bias, amino acid bias and any measure that correlates with any of these parameters.
  • the parameter than effects translation rate is selected from the group consisting of: relative synonymous codons frequencies (RSCF), relative codon-tRNA adaptation and codon typical decoding rate (TDR). Measuring these parameters are well known in the art and examples of such are described herein. Parameters which are equivalent or correlate with these parameters may also be used.
  • RSCF relative synonymous codons frequencies
  • TDR codon typical decoding rate
  • Relative synonymous codons frequencies refers to the frequency at which a codon is used relative to other synonymous codons within a specific reference set.
  • Relative synonymous codons frequencies can be represented as a vector which entries correspond to each one of 61 coding codons (stop codons are excluded):
  • synfi is a subset of indexes in RSCF pointing at codons synonymous to codon i.
  • the term "relative codon-tRNA adaptation" as used herein refers to how well a codon is adapted to the tRNA pool relative to other synonymous codons within a specific reference set.
  • the tRNA pool in a cell can change over time depending on the cellular context. In some embodiments, the tRNA pool changes due to viral infection of the cell. In such a case, the tRNA pool will differ during the early phase of a virus's life cycle and the later stage of the virus's life cycle. In this way, a codon can be well adapted to the tRNA pool in one cellular context, but poorly adapted in another.
  • Relative codon-tRNA adaptation and the tRNA adaptation index (tAI) quantify the adaptation of one codon, or a coding region, respectively, to the tRNA pool.
  • the S vector [sI:U, sG:C, sU:A, sC:G, sG:U, sI:C, si: A, sU:G, sL:A] was defined for E.coli as [0, 0, 0, 0, 1, 0.25, 0.81, 1, 0.71] according to optimization performed previously (Sabi R, et al., DNA Research, 2014, 21:511-525).
  • the absolute adaptiveness value of a codon of type i (1 ⁇ i ⁇ 61 ; stop codons are excluded) to the tRNA pool is defined by:
  • the weight of each of its codons is computed as the ratio between the absolute adaptiveness value of the codon and the maximal absolute adaptiveness value of the synonymous codons for that amino acid:
  • W t is the absolute adaptiveness of codon ⁇ in a sequence
  • syn[i] is a subset of indexes in pointing at codons synonymous to codon i.
  • w t takes values from 0 (not adapted) to 1 (maximally adapted). If the weight value is zero a value of 0.5 is used
  • tAI is the geometric mean of W; (relative codon-tRNA adaptation) over codons of a coding sequence.
  • typical decoding rate refers to the estimated time it takes to decode a specific codon and place the appropriate amino acid in the chain of the protein.
  • the methodology of calculating the typical decoding rate can be found in the Examples section of the present application as well as in the following reference (Dana and Tuller, 2015, G3 [Genes, Genomes, Genetics], 5(1): 73-80), which is hereby incorporated into the current application by reference.
  • GC content refers to the percentage of nitrogenous bases on a nucleic acid molecule that are either guanine or cytosine.
  • the binding of guanine and cytosine is a stronger bond that thymine and adenosine or uracil and adenosine. Calculation of GC content is well known in the art and is the sum of all G and C bases divided by the total number of bases.
  • the term "average repetitive substring index” as used herein refers to the observed overexpression of a substring (a string of nucleotides of at least a length of 3 bases) with in a set of genes that is different than would be expected by chance. As binding sites are coded into gene sequences optimized binding sites are hypothesized to be more frequently used. If that binding site is for the transcription/translation machinery the use of optimized or deoptimized sites for a particular cellular context will effect translation in that cellular context. [066]
  • the term “codon pair bias” as used herein refers to use of a particular codon pair at a rate that is different than would be expected by chance. Codon pair bias is independent from amino acid and codon bias. It will be understood that the codon pair bias is thus an average of codon pair scores over all codon pairs of all viral coding sequences from a set of viral genes.
  • dinucleotide bias refers to use of a particular pair of nucleotides at a rate that is different than would be expected by chance. Dinucleotide bias can be determined for any set of genes, be that a whole virus genome or a subset of viral genes.
  • nucleotide bias and “amino acid bias” as used herein refer to the use of a nucleotide/AA in a set of sequences that is not what would be expected by chance. These biases can be determined for any set of genes, be that a whole virus genome or a subset of viral genes.
  • viruses and “viral” refer to a small infectious agent that replicates only inside the living cells of other organisms.
  • the viruses described herein are capable of infecting all types of life, including animals, plants and microorganisms.
  • the virus is capable of infecting an animal.
  • the virus is capable of infecting a human.
  • the virus is capable of infecting a plant.
  • Non-limiting examples of the virus used under the method of the invention include dengue virus, poliovirus, rhinovirus, influenza virus, severe acute respiratory syndrome (SARS) coronavirus, Human Immunodeficiency Virus (HIV), Hepatitis C Virus (HCV), infectious bronchitis virus, Ebolavirus, Marburg virus, West Nile disease virus, Zika, Epstein-Barr virus (EB V) and yellow fever virus.
  • dengue virus poliovirus, rhinovirus, influenza virus, severe acute respiratory syndrome (SARS) coronavirus
  • HIV Human Immunodeficiency Virus
  • HCV Hepatitis C Virus
  • infectious bronchitis virus Ebolavirus
  • Marburg virus Marburg virus
  • West Nile disease virus West Nile disease virus
  • Zika Zika
  • Epstein-Barr virus EB V
  • yellow fever virus yellow fever virus
  • the virus is human Simplex virus. In some embodiments, the virus is HIV-1. In some embodiments, the virus is Zika virus.
  • the virus is a bacteriophage. In some embodiments, the virus infects mammalian cells. In some embodiments, the virus infects human cells. In some embodiments, the virus is bacteriophage Lambda. In some embodiments, the virus is bacteriophage T4. In some embodiments, the virus is bacteriophage is Pak P3. In some embodiments, the virus is bacteriophage phi29. In some embodiments, the virus is bacteriophage phiYs40. In some embodiments, the virus is bacteriophage Streptococcus DTI. In some embodiments, the virus is bacteriophage Mu. In some embodiments, the virus is bacteriophage phiEco32. In some embodiments, the virus is bacteriophage p23-45. In some embodiments, the virus is bacteriophage phiRl-37.
  • the set of viral genes is more than 1 gene, at least 2 genes, at least 3 genes, at least 5 genes, at least 10 genes, at least 20 genes, at least 30 genes, at least 40 genes, or at least 50 genes.
  • Each possibility represents a separate embodiment of the present invention.
  • the set of viral genes is coordinately expressed. In some embodiments, the set of viral genes is the entire viral genome. In some embodiments, the set of viral genes is the entire viral genome or transcriptome wherein the entire viral genome or transcriptome is coordinately expressed. In some embodiments, the set of viral genes is not the entire viral genome or viral transcriptome. In some embodiments, the set of viral genes is a subset of the viral genome. In some embodiments, the set of viral genes is a subset of the viral transcriptome. In some embodiments, the set of viral genes are selected randomly. In some embodiments, the set of viral genes are specifically selected. In some embodiments, the set of viral genes share a specific parameter. In some embodiments, the set of viral genes share a common trait.
  • the common trait shared by the set of viral genes is selected from the group consisting of: common gene structure, common temporal expression during a virus's life cycle, common function of the encoded proteins, common cellular localization of translation and common cellular localization of the encoded proteins.
  • coordinated expression comprises at least one of coordinated temporal expression, coordinated location of translation, coordinated location of expressed proteins, and coordinated rate of expression.
  • coordinated expression is coordinated temporal expression. In some embodiments, the temporal expression is during a virus's life cycle.
  • gene structure refers to the number, size, and spacing of exons in the gene body.
  • gene structure refers to base pair motifs and structures within the gene body, such as GC repeats, CpG islands, AT rich region, or DNA/RN A/protein binding sites.
  • gene structure refers to elements within the gene that code for specific protein motifs, or functional regions.
  • genes with common gene structure are member of a particular gene or protein family.
  • the term “temporal expression during a virus's life cycle”, refers to a gene's expression occurring during a defined period in the life of a virus.
  • the terms “express” or “expression” as used herein refers to the biosynthesis of a product, including the transcription and/or translation of said gene product or a non-coding RNA.
  • expression of a nucleic acid molecule may refer to transcription of the nucleic acid fragment (e.g., transcription resulting in mRNA or other functional RNA) and/or translation of RNA into a precursor or mature protein (polypeptide).
  • the defined period is selected from the list consisting of: early in the virus life cycle, the intermediate period of the virus life cycle, and late in the virus life cycle.
  • early in the virus's life cycle is during the first 20 minutes after infection.
  • early in the virus's life cycle is during the first IS minutes after infection.
  • early in the virus's life cycle is during the lysogenic portion of an infection.
  • early in the virus's life cycle is following entry into the host cell, but prior to replication.
  • early in the virus's life cycle is before nuclear export of viral mRNAs.
  • late in the virus's life cycle is during the last IS minutes before lysis of the host cell.
  • late in the virus's life cycle is during the last 10 minutes before lysis of the host cell. In some embodiments, late in the virus's life cycle is during the lysis portion of an infection. In some embodiments, late in the virus's life cycle is after replication. In some embodiments, late in the virus's life cycle is after nuclear export of viral mRNAs.
  • One skilled in the art will be able to determine the different stages of a virus's life cycle. More than one way of dividing a specific virus's life cycle may be useful for determining the set of genes to be used for determining codon frequency.
  • viral life cycle and “virus's life cycle” are synonymous and interchangeable, and as used herein refer to the period spanning entry of a virus into a host cell, stably residence within the host (lysogeny), and a lytic phase during which the virus produces progeny viral particles, and lyses and kills the host cell.
  • the virus is bacteriophage Lambda
  • the viral life cycle is approximately 25 minutes long
  • the early period of the viral life cycle is the first IS minutes
  • the late period of the viral life cycle is the last 10 minutes.
  • the term "host cell” refers to any cell into which a virus has been introduced.
  • the host cell is a bacterial cell.
  • the host cell is an E. coli cell.
  • the host cell is a mammalian cell.
  • the host cell is a human cell.
  • common cellular localization of translation refers to the mRNA from a group of genes all being translated at the same location. Translation, generally occurs in the cytoplasm, but can also occur in the nucleus or at specific locations within the cytoplasm. Non- limiting examples of common cellular localization of translation are translation at the ER, perinuclear translation, and translation near the plasma membrane.
  • the term "common function of the encoded proteins” refers to the proteins for which a group of genes code having a common cellular function. Methods of discerning protein function are well known to one skilled in the art and protein function can be found using several commercially available websites and softwares. Examples of such are not limited to but include: the gene ontology consortium (www.geneontology.org), the Ingenuity Pathway Analysis (www.ingenuity.com), Uniprot (www.uniprot.org), and PredictProtein (www.predictprotein.org).
  • the common function is viral structure. In some embodiments, the common function is transcriptional regulation. In some embodiments, the common function is lysongeny.
  • the term "common cellular localization of the encoded proteins” refers to the proteins for which a group of genes code having a common localization in the cell.
  • the common localization is in the plasma membrane, in the nucleus, in the cytosol, in the mitochondria, in the endoplasmic reticulum, in the Golgi apparatus, or in the nuclear membrane. Each possibility represents a separate embodiment of the present invention.
  • the synonymous codon has a parameter that results in faster translation in a particular cellular context relative to the at least one codon that has been substituted. In some embodiments, a plurality of codons having synonymous codons with faster translation rates in a particular cellular context have been substituted. In some embodiments, the at least one codon is substituted to the fastest translating synonymous codon in a particular cellular context.
  • faster translation in a particular cellular context is at least %, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 100%, 150%, 200%, 250%, 300%, 350%, 400%, 450%, 500%, 1000%, or 10000% faster as compared to the at least one codon that has been replaced.
  • optimizing a codon increases translational efficiency of the nucleic acid sequence in a particular cellular context.
  • optimizing a codon increases expression of a polypeptide encoded by the nucleic acid sequence in a particular cellular context.
  • deoptimizing a codon decreases translational efficiency of the nucleic acid sequence in a particular cellular context. In some embodiments, deoptimizing a codon decreases expression of a polypeptide encoded by the nucleic acid sequence in a particular cellular context.
  • optimizing the codons increases expression by at least 1%, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 100%, 150%, 200%, 250%, 300%, 350%, 400%, 450%, 500%, 1000%, or 10000% relative to expression of the target sequence.
  • Each possibility represents a separate embodiment of the present invention.
  • the synonymous codon has a parameter that results in slower translation in a particular cellular context relative to the at least one codon that has been substituted.
  • a plurality of codons having synonymous codons with slower translation rates in a particular cellular context have been substituted.
  • the at least one codon is substituted to the slowest translating synonymous codon in a particular cellular context.
  • slower translation in a particular cellular context is at least %, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 100%, 150%, 200%, 250%, 300%, 350%, 400%, 450%, 500%, 1000%, or 10000% slower as compared to the at least one codon that has been replaced.
  • deoptimizing the codons decreases expression by at least 1%, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 100%, 150%, 200%, 250%, 300%, 350%, 400%, 450%, 500%, 1000% or 10000% relative to expression of the target sequence.
  • deoptimizing the codons decreases expression to zero.
  • deoptimizing the codons abolishes, or completely removes expression of the protein.
  • a modified virus comprising anyone of the above described nucleic acid molecules.
  • the virus is a virulent virus.
  • the virus infects mammals. In some embodiments, the virus infects humans.
  • the modified virus is for use in creating a vaccine. In some embodiments, the modified virus is for use as an oncolytic virus. In some embodiments, the modified virus is for use as a carrier for a therapeutic agent.
  • the modified virus has a coding sequence that is an essential viral sequence expressed during a defined period of a virus's life cycle and the set of genes used for assessing the frequency of codons is a plurality of viral genes expressed during the defined period of a virus's life cycle. In some embodiments, the modified virus has a coding sequence that is part of an essential viral sequence.
  • the essential viral gene is required for viral replication.
  • the essential viral gene is a structural protein.
  • the essential viral gene is required for lysis of a host cell.
  • the essential viral gene is a capsid.
  • the essential viral gene is an enzyme.
  • the essential viral gene is reverse transcriptase.
  • the essential viral gene is RNA dependent RNA polymerase (RDRP).
  • an attenuated form of a virus comprising a nucleic acid molecule with at least one codon substituted for a less frequent codon, and wherein the coding sequence is an essential viral sequence endogenous ly expressed during a defined period of a virus's life cycle, said set of genes used for calculating codon frequency is a plurality of viral genes expressed during the same defined period of the virus's life cycle and wherein the virus is devoid of the endogenous viral sequence.
  • the viral genome is altered such that the endogenous viral sequence is completely absent. In some embodiments, the viral genome is altered such that the endogenous viral sequence is not expressed.
  • General methods for disruption of the genome as well as general methods in molecular and cellular biochemistry can be found in such standard textbooks as Molecular Cloning: A Laboratory Manual, 3rd Ed. (Sambrook et al., HaRBor Laboratory Press 2001); Short Protocols in Molecular Biology, 4th Ed. (Ausubel et al. eds., John Wiley & Sons 1999); Protein Methods (Bollag et al., John Wiley & Sons 1996); Nonviral Vectors for Gene Therapy (Wagner et al.
  • the viral genome When the viral genome is an RNA genome, they may be isolated from virions or from infected cells, converted to DNA ("cDNA") by the enzyme reverse transcriptase, possibly modified as desired, and reverted, usually via the RNA intermediate, back into infectious viral particles. Most commonly, the entire cDNA copy of the genome is cloned immediately downstream of a phage 17 RNA polymerase promoter that allows the in vitro synthesis of genome RNA, which is then transfected into cells for generation of virus (van der Wert, et al., 1986). Alternatively, the same DNA plasmid may be transfected into cells expressing the 17 RNA polymerase in the cytoplasm.
  • the modifying is achieved by de novo synthesis of DNA containing the synonymous codons and substitution of the corresponding region of the genome with the synthesized DNA.
  • the entire genome is substituted with the synthesized DNA.
  • a portion of the genome is substituted with the synthesized DNA.
  • a vaccine composition for inducing a protective immune response in a subject, the vaccine composition comprising one of the above described attenuated viruses and a pharmaceutically acceptable carrier or adjuvant.
  • the composition is an immunogenic composition.
  • vaccine composition and “vaccine” as used herein are interchangeable and refers to a product, the administration of which is intended to elicit an immune response that is capable of preventing and/or lessening the severity of one or more viral infections.
  • an attenuated virus of the invention where used to elicit a protective immune response (i.e. immunize) in a subject or to prevent a subject from becoming afflicted with a virus-associated disease, is administered to the subject in the form of a composition additionally comprising a pharmaceutically acceptable carrier.
  • a pharmaceutically acceptable carrier refers to any component of a pharmaceutical composition that is not the active agent.
  • pharmaceutically acceptable carrier refers to non-toxic, inert solid, semi-solid liquid filler, diluent, encapsulating material, formulation auxiliary of any type, or simply a sterile aqueous medium, such as saline.
  • sugars such as lactose, glucose and sucrose, starches such as corn starch and potato starch, cellulose and its derivatives such as sodium carboxymethyl cellulose, ethyl cellulose and cellulose acetate; powdered tragacanth; malt, gelatin, talc; excipients such as cocoa butter and suppository waxes; oils such as peanut oil, cottonseed oil, safflower oil, sesame oil, olive oil, corn oil and soybean oil; glycols, such as propylene glycol, polyols such as glycerin, sorbitol, mannitol and polyethylene glycol; esters such as ethyl oleate and ethyl laurate, agar; buffering agents such as magnesium hydroxide and aluminum hydroxide; alginic acid; pyrogen-free water; isotonic saline, Ringer's solution; ethyl
  • substances which can serve as a carrier herein include sugar, starch, cellulose and its derivatives, powered tragacanth, malt, gelatin, talc, stearic acid, magnesium stearate, calcium sulfate, vegetable oils, polyols, alginic acid, pyrogen-free water, isotonic saline, phosphate buffer solutions, cocoa butter (suppository base), emulsifier as well as other non-toxic pharmaceutically compatible substances used in other pharmaceutical formulations.
  • Wetting agents and lubricants such as sodium lauryl sulfate, as well as coloring agents, flavoring agents, excipients, stabilizers, antioxidants, and preservatives may also be present.
  • any non-toxic, inert, and effective carrier may be used to formulate the compositions contemplated herein.
  • Suitable pharmaceutically acceptable carriers, excipients, and diluents in this regard are well known to those of skill in the art, such as those described in The Merck Index, Thirteenth Edition, Budavari et al., Eds., Merck & Co., Inc., Rahway, NJ. (2001); the CTFA (Cosmetic, Toiletry, and Fragrance Association) International Cosmetic Ingredient Dictionary and Handbook, Tenth Edition (2004); and the "Inactive Ingredient Guide," U.S. Food and Drug Administration (FDA) Center for Drug Evaluation and Research (CDER) Office of Management, the contents of all of which are hereby incorporated by reference in their entirety.
  • CTFA Cosmetic, Toiletry, and Fragrance Association
  • Examples of pharmaceutically acceptable excipients, carriers and diluents useful in the present compositions include distilled water, physiological saline, Ringer's solution, dextrose solution, Hank's solution, and DMSO. These additional inactive components, as well as effective formulations and administration procedures, are well known in the art and are described in standard textbooks, such as Goodman and Gillman's: The Pharmacological Bases of Therapeutics, 8th Ed., Gilman et al. Eds. Pergamon Press (1990); Remington's Pharmaceutical Sciences, 18th Ed., Mack Publishing Co., Easton, Pa.
  • compositions may also be contained in artificially created structures such as liposomes, ISCOMS, slow-releasing particles, and other vehicles which increase the half -life of the peptides or polypeptides in serum.
  • liposomes include emulsions, foams, micelies, insoluble monolayers, liquid crystals, phospholipid dispersions, lamellar layers and the like.
  • Liposomes for use with the presently described peptides are formed from standard vesicle-forming lipids which generally include neutral and negatively charged phospholipids and a sterol, such as cholesterol.
  • the selection of lipids is generally determined by considerations such as liposome size and stability in the blood.
  • a variety of methods are available for preparing liposomes as reviewed, for example, by Coligan, J. E. et al, Current Protocols in Protein Science, 1999, John Wiley & Sons, Inc., New York, and see also U.S. Pat. Nos.4,235,871, 4,501,728, 4,837,028, and 5,019,369.
  • the carrier may comprise, in total, from about 0.1 % to about 99.99999% by weight of the pharmaceutical compositions presented herein.
  • a method for vaccinating a subject at risk of viral infection comprising, administering to said subject the vaccine composition described above.
  • subject at risk of viral infection includes but is not limited to a subject that due a likelihood of future exposure to a virus, future exposure to an individual or animal infected with the virus, or future exposure to biological mater infected with the virus, is at a higher risk than the general population of contracting the virus.
  • nucleic acid molecule optimized for expression in a particular cellular context comprising:
  • nucleic acid molecule optimized for expression in a particular cellular context thereby producing a nucleic acid molecule optimized for expression in a particular cellular context.
  • the computing of average relative frequency comprises, selecting a reference set (R) of viral genes expressed during a defined period, for each amino acid coded for by said coding sequence (A), computing a relative frequency of synonymous codons in each reference sequence k:
  • substituting a codon with a faster rate of translation comprises, a. for each codon of said coding sequence (C A T ) replacing C A T with a synonymous codon C Ai j according to
  • optimizing the codons comprises the steps of:
  • optimizing said codons comprises optimizing the expression levels of the wildtype sequence s with respect to the codons Typical Decoding Rate (TDR) basing on available ribosomal profiling data.
  • TDR Typical Decoding Rate
  • a statistical mode which takes into consideration the skewed nature of the ribose read count distribution can be used.
  • This model describes the readcount histogram of each codon as an output of a random variable which is a sum of two random variables: a normal and an exponential variable.
  • the distribution of this new random variable includes three parameters and is called EMG distribution.
  • the typical codon decoding time was described by the normal distribution with two parameters: mean ( ⁇ ) and standard deviation 6; the ⁇ parameter represents the location of the mean of the theoretical Gaussian component that should be obtained if there are no phenomena such as pauses/ biases/ ribosomal traffic jams; ⁇ represents the width of the Gaussian component.
  • the exponential distribution has one parameter ⁇ which represents the skewness of the readcount distribution due to reasons such as ribosomal jamming caused by codons with different decoding times, extreme pauses, incomplete halting of the ribosomes, biases in the experiment, etc.
  • the EMG is defined as follows:
  • These three parameters may be estimated for each codon at different replication stages based on time dependent ribosome profiling data by fitting the suggested model to the given read count distribution (e.g. using the maximal likelihood estimation or any other algorithm). is defined to
  • TDR Typical Decoding Rate
  • a method for producing a nucleic acid molecule deoptimized for expression in a particular cellular context comprising:
  • nucleic acid molecule optimized for expression in a particular cellular context thereby producing a nucleic acid molecule optimized for expression in a particular cellular context.
  • substituting a codon with a slower rate of translation comprises, b. for each codon of said target sequence (C Ai t ) replacing C Ai t with a synonymous codon C A according to
  • deoptimizing said codons comprises the steps of:
  • the deoptimization comprises creating a gradient library of deoptimized variants of the coding sequence with different levels of deoptimization where codon of said target sequence (C Ai i ) is replaced with a synonymous codon
  • the gradient library may also include Cartesian products of codons with different ranks of relative frequencies:
  • deoptimizing said codons comprises replacing the codons of the coding sequence with synonymous codons chosen according to the "inverse rank" rule: if the wildtype codon C A i is the most frequent codon in the reference set it is replaced by the less frequent synonymous codon C AiJ -; if it is the second most frequent it is replaced with the second less frequent; if it is the third most frequent it is replaced with the third less frequent. ... , the second less frequent is replaced by the second most frequent; the less frequent is replaced by the most frequent.
  • the particular cellular context is selected from the group consisting of: a specific period of a cell's life cycle, during cellular division, during cellular stress, during apoptosis, during viral infection, or during viral lysogeny of the cell. It should be understood, that any change in cellular context may result in a change in the tRNA pool. In some embodiments, the particular cellular context is a pre-determined cellular context.
  • the above recited methods further comprise confirming that the nucleic acid or virus has been optimized, deoptimized or modified.
  • confirming optimization, deoptimization or modification comprises expressing the selected sequence and the substituted sequence in the virus, measuring protein expression during the defined time period of the virus's life cycle and comparing the protein levels produced by the selected sequence and substituted sequence. Expression within the virus and measuring protein expression are routine laboratory practices and would be well known to one skilled in the art. Detailed descriptions of these methods can be found in the publications previously enumerated in this application.
  • said viral sequence is an essential viral sequence
  • said synonymous codon with a different translation rate has a slower translation rate, thereby producing an attenuated virus.
  • the present invention may be a system, a method, and/or a computer program product.
  • the computer program product may include a computer readable storage medium (or media) having computer readable program instructions thereon for causing a processor to carry out aspects of the present invention.
  • the computer readable storage medium can be a tangible device that can retain and store instructions for use by an instruction execution device.
  • the computer readable storage medium may be, for example, but is not limited to, an electronic storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the foregoing.
  • a non-exhaustive list of more specific examples of the computer readable storage medium includes the following: a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a static random access memory (SRAM), a portable compact disc read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanically encoded device having instructions recorded thereon, and any suitable combination of the foregoing.
  • a computer readable storage medium is not to be construed as being transitory signals per se, such as radio waves or other freely propagating electromagnetic waves, electromagnetic waves propagating through a waveguide or other transmission media (e.g., light pulses passing through a fiber-optic cable), or electrical signals transmitted through a wire. Rather, the computer readable storage medium is a non- transient (i.e., not-volatile) medium.
  • Computer readable program instructions described herein can be downloaded to respective computing processing devices from a computer readable storage medium or to an external computer or external storage device via a network, for example, the Internet, a local area network, a wide area network and/or a wireless network.
  • the network may comprise copper transmission cables, optical transmission fibers, wireless transmission, routers, firewalls, switches, gateway computers and/or edge servers.
  • a network adapter card or network interface in each computing/processing device receives computer readable program instructions from the network and forwards the computer readable program instructions for storage in a computer readable storage medium within the respective computing/processing device.
  • Computer readable program instructions for carrying out operations of the present invention may be assembler instructions, instruction-set-architecture (ISA) instructions, machine instructions, machine dependent instructions, microcode, firmware instructions, state-setting data, or either source code or object code written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or the like, and conventional procedural programming languages, such as the "C" programming language or similar programming languages.
  • the computer readable program instructions may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server.
  • the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet Service Provider).
  • electronic circuitry including, for example, programmable logic circuitry, field-programmable gate arrays (FPGA), or programmable logic arrays (PLA) may execute the computer readable program instructions by utilizing state information of the computer readable program instructions to personalize the electronic circuitry, in order to perform aspects of the present invention.
  • These computer readable program instructions may be provided to a processor of a general-purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
  • These computer readable program instructions may also be stored in a computer readable storage medium that can direct a computer, a programmable data processing apparatus, and/or other devices to function in a particular manner, such that the computer readable storage medium having instructions stored therein comprises an article of manufacture including instructions which implement aspects of the function act specified in the flowchart and/or block diagram block or blocks.
  • the computer readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable apparatus or other device to produce a computer implemented process, such that the instructions which execute on the computer, other programmable apparatus, or other device implement the functions/acts specified in the flowchart and/or block diagram block or blocks.
  • each block in the flowchart or block diagrams may represent a module, segment, or portion of instructions, which comprises one or more executable instructions for implementing the specified logical function(s).
  • the functions noted in the block may occur out of the order noted in the figures.
  • two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved.
  • Transcript sequences were obtained from EnsEMBL for E. coli (K-12 MG16SS release 121, accessed 28/07/15) and from NCBI for the lambda phage (J02459, accessed 07/12/15).
  • Ribo-seq reads mapping Ribosome footprint sequences were obtained from previously published results (Liu XQ, et al., Proceedings of the National Academy of Sciences of the United States of America, 2013, 110:11928-11933) (GSE47509, induction 0-20min).
  • the poly-A adaptors were trimmed from the reads using the Cutadapt program (version 1.8.3), and the Bowtie program (version 1.1.1) was used to map them to the E. coli-lambda transcriptome.
  • the Cutadapt program version 1.8.3
  • Bowtie program version 1.1.1
  • the remaining reads were mapped to the transcriptome with Bowtie parameters '-v 2 -a—strata—best—norc -m 200'. Alignments were extended to their maximal length by comparing the poly A adaptor with the aligned transcript until reaching the maximal allowed error (2 mismatches across the read, with 3' -end mismatches avoided). Reads were either longer than 31 nucleotides (nt) or shorter than 21 nt. Unique alignments were first assigned to the ribosome occupancy profiles. For multiple alignments, the best alignments in terms of number of mismatches were kept. Then, multiple aligned reads were distributed between locations according to the distribution of unique ribosomal reads in the respective surrounding regions.
  • a 100-nt window was used to compute the read count density RCDi (total read counts in the window divided by length, based on unique reads) in vicinity of the M multiple aligned positions in the transcriptome, and the fraction of a read assigned to each position was RCDj/ ⁇ j t RCD j .
  • the location of the A-site was approximated by an 11-nt shift from the 5' end of the aligned read. This shift maximized the correlation between MTDR (described below) and the observed read densities per E. coli gene.
  • Synonymous codon composition of a coding sequence was represented by a 61- dimensional vector of relative synonymous codons frequencies (RSCF) of each one of 61 coding codons (stop codons are excluded):
  • syn[i] is a subset of indexes in RSCF pointing at codons synonymous to codon i.
  • Each viral sequence was assigned a group label corresponding to its temporal expression stage (Early/Late) (according to the classification known in the literature).
  • the tendency of sequences to cluster according to the codons usage in two different clusters corresponding to their temporal expression stages (early late) was measured using the Davies-Bouldin score (DBS). This score is based on a ratio of within-cluster and between-cluster distances and is defined as:
  • k - is the number of evaluated clusters, is the within-to-between cluster distance ratio for the i-th and j-th clusters;
  • d L is the average Euclidian distance between each point in the i-th cluster and the centroid of the i-th cluster;
  • d j is the average Euclidian distance between each point in the j-th cluster and the centroid of the j-th cluster;
  • dy is the Euclidean distance between the centroids of the ith and jth clusters.
  • j represents the worst-case within- to-between cluster ratio for cluster i.
  • the optimal clustering solution has the smallest Davies- Bouldin score value.
  • PC A principal component analysis
  • Codons composition of a coding sequence was represented by codons relative frequency vector.
  • each entry i corresponds to one specific codon and contains its relative frequency with respect to codon usage in some reference set of coding sequences (e.g. E. coli coding regions, bacteriophage early genes, bacteriophage late genes, etc.):
  • the clustering analysis was performed on the codons relative frequency vectors of each viral coding sequence (with all viral genes as a reference set). A principal component analysis was applied to project the r vectors to a plane spanned by their first two principal components. Two different clusters separating between early and late genes were distinguished (Fig. 2A) using the Davies-Bouldin score. The significance of clusters separation was assessed by comparing this score to the randomized scores obtained from 100 permutations of gene group labels (early or late).
  • codons usage frequency distance CUFD measures similarity between genes in terms of codon and amino acid usage. Codon usage frequency vectors were computed by counting all appearances ni of a codon i in the ORF, and dividing by the total codon count:
  • dKL is the Kullback-Leibler divergence. Larger values related to less similarity / larger distance.
  • tAI quantifies the adaptation of a coding region to the tRNA pool.
  • tCGNi j be the copy number of the j-th anti-codon that recognizes the i-th codon and let Si j be the selective constraint of the codon-anti-codon coupling efficiency.
  • the S vector [sI:U, sG:C, sU:A, sC:G, sG:U, sI:C, sI:A, sU:G, sL:A] was defined for E.coli as [0, 0, 0, 0, 1, 0.25, 0.81, 1, 0.71] according to optimization performed previously (Sabi R, et al., DNA Research, 2014, 21:511-525).
  • the absolute adaptiveness value of a codon of type i (1 ⁇ i ⁇ 61; stop codons are excluded) to the tRNA pool is defined by:
  • the weight of each of its codons is computed as the ratio between the absolute adaptiveness value of the codon and the maximal absolute adaptiveness value of the synonymous codons for that amino acid:
  • W t is the absolute adaptiveness of codon i in a sequence
  • syn[i] is a subset of indexes in pointing at codons synonymous to codon i.
  • W takes values from 0 (not adapted) to 1 (maximally adapted). If the weight value is zero a value of 0.5 is used
  • tAI is the geometric mean of w t (relative codon-tRNA adaptation) over codons of a coding sequence.
  • Ribosome profiles for E. coli and Bacteriophage Lambda expressed genes were reconstructed.
  • the ribosome profiling method produces ribosome footprint counts that are proportional to the time spent in decoding each codon of all translated transcripts in a genome, at single nucleotide resolution.
  • genes profiles with fewer than 30 percent non-zero read counts were further filtered. The first and last 20 codons were excluded when determining these thresholds or when calculating the average RCs per ORF.
  • J is the number of codons in the gene and j is the index of a codon.
  • Each NFC distribution describes the probability (y-axis) of observing each of the codon' s NFC values (x-axis) in the ORFs of the analyzed organism.
  • the typical codon decoding time was described by the normal distribution with two parameters: mean ( ⁇ ) and standard deviation ( ⁇ ); the ⁇ parameter represents the location of the mean of the theoretical Gaussian component that should be obtained if there are no phenomena such as pauses/ biases/ ribosomal traffic jams; ⁇ represents the width of the Gaussian component.
  • the exponential distribution has one parameter ⁇ , which represents the skewedness of the NFC distribution due to reasons such as ribosomal jamming caused by codons with different decoding times, extreme pauses, incomplete halting of the ribosomes, biases in the experiment, etc.
  • the EMG is defined as follows:
  • MTDR Mean typical decoding rate
  • i is an index of a codon and L is the gene length in codon unit.
  • RTEC Relative translation elongation efficiency coefficient
  • Minimum free folding energy is a thermodynamic energy involved in maintaining a secondary structure available to perform physical work whilst being released, and thus is characterized by nonpositive values. mRNA secondary structure is believed to be in the most stable conformation when minimum amount of free energy is exerted (the MFE obtains the most negative value).
  • the local MFE-profiles were constructed by applying a 39 nt length sliding window to a genomic sequence: in each step the MFE of a local subsequence enclosed by the corresponding window was calculated by Vienna (v. 2.1.9) package RNAfold function with default parameters. This function predicts the MFE and the associated secondary structure for the input RNA sequence using a dynamic programming based on the thermodynamic nearest-neighbor approach (the Zucker algorithm).
  • the mean MFE-profile was obtained by averaging the MFE-profiles of all genes (in a position wise manner).
  • 100 randomized mean MFE-profiles were computed by grouping the randomized MFE-profiles of all genes in 100 groups, each group contains a different variant for each gene, and then averaging the profiles in each group in a position-wise manner.
  • mean MFE values were computed for each gene over all windows (by averaging the values in the corresponding MFE-profiles) and compared to the mean MFE values obtained from the corresponding 100 randomized profiles. For each gene its mean MFE value was calculated and as were an average of 100 mean MFE values from its randomized variants; the distributions of the wildtype and randomized mean MFE values of different genes were compared using Wilcoxon signed-rank test. Early and late genes were analyzed separately.
  • ARS Average Repetitive Substring
  • This measure is based on the assumption that evolution shapes the organismal coding sequences (and other part of the gene) to improve their interaction with the intra-cellular gene expression machinery. Since these interactions are mediated via binding of the gene expression machinery (e.g. translation/transcription factors, RNAP, ribosomes, RNA binding proteins, etc.) to the genetic material (DNA or RNA), the genetic material tend to have optimized binding sites (which are sub-sequences of nucleotides). It was expected that binding sites will appear in many coding regions and that more optimal binding sites will tend to appear more times in the genome.
  • the gene expression machinery e.g. translation/transcription factors, RNAP, ribosomes, RNA binding proteins, etc.
  • ARS index was computed for each bacteriophage gene in comparison to the host (E. coli) and in comparison to the rest of the viral genes.
  • the algorithm of ARS index is based on the following steps: (i) For each position i in the coding sequence S find the longest substring Sji that starts in that position, and also appears in at least one of the coding sequences of the reference genome (E. coli/viral). (ii) Let jSj denote the length of a sequence S; the ARS index of S is the mean length of all the substrings Sl'. ARS—j- . Rare codons analysis
  • Rare codons in a reference set of coding sequences were defined as codons with the relative synonymous frequency ⁇ 0.2.
  • Three reference sets were used: E. coli coding sequences, bacteriophage early coding sequences, bacteriophage late coding sequences.
  • a rare codons score (RCS) for a specific early late/E. coli coding sequence with respect to a reference set of all early/all late/all E. coli coding sequences is defined as a percentage of amino acids in that sequence encoded by a rare codon out of all amino acids that are encoded by at least one rare codon in the corresponding reference set (if an AA is not encoded by codons that are rare in the reference set we exclude it from the analysis):
  • N is the total number of such codons.
  • Human Viruses analyzed in this study include Herpes viruses, papilloma viruses, Polyomavirus and HIV.
  • the analyzed bacteriophages include: bacteriophage Lambda, bacteriophage T4, bacteriophage Pak P3, bacteriophage phi29, bacteriophage T7, bacteriophage phiYs40, bacteriophage Fah, bacteriophage xplO, bacteriophage Streptococcus DTI, bacteriophage Streptococcus 2972, bacteriophage Mu, bacteriophage phiC31, bacteriophage phiEco32, bacteriophage p23-45 and bacteriophage phiRl-37.
  • Effective number of codons is a measure that quantifies how far the synonymous codon usage of a gene departs from what is expected under the assumption of uniformity.
  • ENC is a measure that quantifies how far the codon usage of a gene departs from equal usage of synonymous codons. It can be calculated from codon usage data alone and is independent of gene length and amino acid (AA) composition. ENC can take values from 20, in the case of extreme bias where one codon is exclusively used for each aa, to 61 when the use of alternative synonymous codons is equally likely and is defined by:
  • xi is the number of synonymous codons of each type in the sequence
  • n is the number of times the AA appears in the sequence
  • p is the frequency (/probability) of each codon.
  • ENC for the sequence (e.g., gene):
  • GC-content is the percentage of nitrogenous bases on a DNA or RNA molecule that are either guanine or cytosine.
  • GC-content (or guanine-cytosine content) is the percentage of nitrogenous bases on a DNA or RNA molecule that are either guanine or cytosine (from a possibility of four different nucleotides) and is defined as: Where F() is a number of occurrences. We calculated the GC content in gene levels for all viruses, or for sets of viral genes.
  • Codon pair bias (CPB). To quantify the CPB, one follows and defines a codon pair score (CPS) as the log ratio of the observed over the expected number of occurrences of this codon pair in the coding sequence. The CPB of a virus is then defined as an average CPSs over all codon pairs comprising all viral coding sequences. To achieve independence from amino acid and codon bias, the expected frequency is calculated based on the relative proportion of the number of times an amino acid is encoded by a specific codon:
  • codon pair bias (CPB) of a virus, or set of viral genes, is than defined as an avenge of codon pair scores over all codon pairs comprising the viral coding sequences:
  • DNTB Dinucleotide bias
  • DNTS dinucleotide score
  • the DNTB of a virus is defined as an average of DNTSs over all dinucleotides comprising all viral coding sequences. This can also be calculated for a set of viral genes.
  • the dinucleotide pair bias (DNTB) of a virus is defined as an average of dinucleotide scores over all dinucleotides comprising all viral sequences: [0187] Nucleotide (NTB) and amino acid (AAB) biases are defined as a normalized Shannon entropy over the frequencies of the nucleotides / AA in a genomic sequence and is defined as:
  • F is the frequency of nucleotides or amino acids and the number of unique symbols equals to 4 in the case of nucleotides or 20 in the case of amino acids. This measure takes values between 0 and 1, and describes how dispersed the distribution of the alphabet elements is: higher values correspond to more uniform nucleotide usage; lower values correspond to more biased nucleotide, indicating that some nucleotides/amino acids are preferred (positions are conserved)
  • a and B be two proteomes of two different viruses consisting of n and m proteins correspondingly.
  • the total length of each proteomes is:
  • proteome B as:
  • l t (j) is the length of the longest starting at position j in c3 ⁇ 4 that exactly matches a starting at
  • the ARS score of the entire proteome A with respect to the proteome B is defined as:
  • the first equality is a weighted average of ARS scores ⁇ ,( ⁇ ⁇ , B) over all proteins in A; for each protein, c3 ⁇ 4 E A its weight is defined as its relative lengths
  • the second equality means
  • L(A,A) should be computed by the formula of L(A,B).
  • ds is the ARS induced distance for tree building.
  • the tree was built using neighbor joining algorithm.
  • Bacteriophage lambda is a well-known and studied member of the Siphoviridae family of double-stranded DNA viruses in the Caudovirales order (also known as "tailed bacteriophages" due to their characteristic form). During its lifecycle, this phage either stably resides within the genome of its E. coli host through lysogeny, or enters into a lytic phase (which lasts about 25 minutes) during which it produces progeny viral particles, and lyses and kills the host cell.
  • the genome size of bacteriophage lambda is about SO kilobases (kb) and includes 66 known genes that were analyzed in this study. These genes were divided into two groups, “early” and “late” according to the stages in the lytic phase when their expression is dominant.
  • EXAMPLE 3 Differential codon usage in early and late genes can be partially explained by adaptation of translation elongation efficiency to different bacteriophage developmental stages
  • MTDR Mean Typical Decoding Rate
  • RTEC relative translation elongation efficiency coefficient
  • Figure 2D describes the per codon correlation between the TOR and relative synonymous codon frequencies (RSCF) for the two bacteriophage gene groups (early and late) and the E. coli genes at different time points. As can be seen, the correlation is higher and significant for the early points in the case of the E. coli and early bacteriophage genes. For the late genes, the correlation is significant only at the initial point.
  • RSCF relative synonymous codon frequencies
  • EXAMPLE 4 Selection for translation efficiency in bacteriophage genes may be partially explained by adaptation to E. coli tRNA pool and the fact that it changes during bacteriophage development
  • Figure 3B describes the per condition partial correlation (controlling for gene length) between rare codons score and mean read counts for the two viral gene groups (early and late) and the E. coli genes. As can be seen, the correlation decreases in the case of the early genes and E. coli genes and increases for late genes Our analysis demonstrates that early late genes with rare early late genes codons tend to be lowly expressed at the early late stages respectively.
  • EXAMPLE 5 Codon bias in other bacteriophages and human viruses
  • viruses undergo an extensive evolutionary selection for adaptation to their host environment, and thus it can be assumed that their codon composition reflects an efficient adaptation of the viral machinery at specific replication stages and therefore may be used as a reference set for codons substitution.
  • 14 additional phages with known temporal gene classification were analyzed.
  • each coding sequence was represented by its relative synonymous codons frequencies (RSCF) - a 61 dimensional vector expressing each sense codon by its frequency in that sequence normalized relative to the frequencies of other synonymous codons coding for the same AA.
  • RSCF relative synonymous codons frequencies
  • the mean distance between all pairs of phages was compared with the signal to the distribution of mean pairwise distances in 100 randomly sampled groups of 7 viruses (the number of viruses in each sample is equal to the number of viruses in the test group). No significant differences in the test group as compared to the randomly sampled groups with respect to genomic similarity was found (empiric p-value - 0.55).
  • Viruses undergo an extensive evolutionary selection for adaptation to their host's cell environment, and thus it can be assumed that their codon composition reflects an efficient adaptation of the viral genes to specific intracellular conditions (e.g. in terms of gene expression factors such as tRNA molecules, AA concentration, etc) that are prevalent in different gene expression stages, in accordance with the reported results.
  • EXAMPLE 6 Weaker separation between synonymous codon usage in early and late genes in human viruses
  • HIV-1 Human Immunodeficiency Virus 1
  • Figure 7 shows the principal component analysis (PCA) for the 11 human viruses.
  • RNA Ribonucleic Acid
  • RNA Ribonucleic Acid
  • This phenomenon can be related, for example, to the contribution of the CpG stacking basepairs to RNA folding and/or to the enhanced innate immune responses to viruses with elevated CpG.
  • the stability of the RNA secondary structures can also be affected by the genomic composition of nucleotides and in particular by GC content.
  • nucleotide compositions and AA usage bias may affect, among others, the synthesis of viral molecules, and the function and structure of the encoded proteins.
  • codon bias at defined periods of a virus's life cycle could enhance translational elongation and thus expression
  • genes viral genes, human genes, heterologous transgenes
  • codons would be selected with frequencies more similar to the set of early or late genes than a random set of genes. This similarity could be based on any metric of codon bias e.g. L distance, Euclidean distance, LI, etc. or on any codon bias indexes (e.g. CAI).
  • Target sequence t a wild-type endogenous viral coding sequences or a coding sequence originating in some heterologous transgene.
  • a set of reference sequences R viral coding sequences grouped by expression period (early/late/intermediate)
  • R [UL2, UL5, UL8, UL12, UL23, UL29, UL30, UL39, UL40, UL42, UL50, UL52, US3] - a reference set of Simplex virus early genes.
  • the optimization evaluation function L was defined as: (1) Euclidean distance; (2) Codon Adaptation Index (CAI)
  • Gene can also be designed that have deoptimized codons with respect to viral genes expressed during a defined period and this would result in minimal protein expression during a defined time period. To achieve this, codons would be selected with frequencies maximally dissimilar to the set of early or late genes than a random set of genes. This dissimilarity could be based on any metric of codon bias e.g. KL distance, Euclidean distance, LI, etc. or on any codon bias indexes (e.g. CA1).
  • RCDI Relative Codon Deoptimization Index
  • F s (Ci) is the observed relative frequency in the tested sequence of each codon i out of all synonymous codons for the same amino acid (0 to 1)
  • F R (Q) is the relative frequency observed in the Reference set of each codon i * out of all synonymous codons for that amino acid (0 to 1)
  • NQ is the number of occurrences of that codon / in the sequence
  • N is the total number of codons (amino acids) in the sequence. The higher the index, the more extensive deoptimization is.
  • R [UL2, UL5, UL8, UL12, UL23, UL29, UL30, UL39, UL40, UL42, UL50, UL52, US3] - a reference set of Simplex virus early genes (Roizman, 1996)
  • the deoptimization evaluation function L was defined as: (1) Euclidean distance; (2) Codon Adaptation Index (CAI); (3) Relative Codon Deoptimization Index (RCDI)
  • Codon alteration was now applied to the Zika virus. 6 modified variants of the virus Zika were designed; 4 variants comprised codon alterations that changed the folding energy of the Zika mRNA and 2 variants comprised codons altered to synonymous codons with a different relative frequency. As all Zika proteins are translated together, the codon substitutions were made based on a comparison to the entire Zika genome of various Zika strains. To test virus virulence a Focus Forming Assay (FFA) was performed. The monkey Vero cell line, which is infectible by Zika, was grown in a monolayer and then infected with the modified virus as well as wild-type virus. After 3 days, infected cells were visualized by immunostaining for the virus (Fig.
  • FFA Focus Forming Assay
  • variants 32 (RSCF) and 104, 103 and 105 (folding) were all predicted to have reduced virulence and indeed all showed reduced size of virally infected foci.
  • variant 29 (folding) resulted in only a negligible increase in virally infected foci, while variant 97 which was predicted to have increase virulence did indeed show larger viral foci.
  • Viral titers from the infected cells were quantified and the 4 attenuated variants did indeed have reduced viral titers (Fig. 11B).
  • Variant 29 showed a very slight increase in viral titers, and variant 97 showed a marked increase at early time points that appeared similar to the wild-type by day 3 (Fig. 11B).

Landscapes

  • Health & Medical Sciences (AREA)
  • Genetics & Genomics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Chemical & Material Sciences (AREA)
  • Zoology (AREA)
  • Organic Chemistry (AREA)
  • Biotechnology (AREA)
  • Wood Science & Technology (AREA)
  • Biomedical Technology (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Microbiology (AREA)
  • Biophysics (AREA)
  • Biochemistry (AREA)
  • Molecular Biology (AREA)
  • Plant Pathology (AREA)
  • Virology (AREA)
  • Medicinal Chemistry (AREA)
  • Immunology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Medical Informatics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Micro-Organisms Or Cultivation Processes Thereof (AREA)

Abstract

L'invention concerne des molécules d'acide nucléique comprenant une séquence de codage à au moins un codon substitué par un codon synonyme, une forme modifiée d'un virus comprenant les molécules d'acide nucléique de l'invention, des méthodes de production de ces molécules d'acide nucléique, et des virus.
PCT/IL2018/050093 2017-01-25 2018-01-25 Séquences d'acide nucléique synthétiques virales et leur utilisation WO2018138727A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US16/480,671 US20190359990A1 (en) 2017-01-25 2018-01-25 Viral synthetic nucleic acid sequences and use thereof
EP18745104.2A EP3574097A4 (fr) 2017-01-25 2018-01-25 Séquences d'acide nucléique synthétiques virales et leur utilisation

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201762450115P 2017-01-25 2017-01-25
US62/450,115 2017-01-25

Publications (1)

Publication Number Publication Date
WO2018138727A1 true WO2018138727A1 (fr) 2018-08-02

Family

ID=62979086

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/IL2018/050093 WO2018138727A1 (fr) 2017-01-25 2018-01-25 Séquences d'acide nucléique synthétiques virales et leur utilisation

Country Status (3)

Country Link
US (1) US20190359990A1 (fr)
EP (1) EP3574097A4 (fr)
WO (1) WO2018138727A1 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020077395A1 (fr) * 2018-10-16 2020-04-23 Griffith University Vaccin antiviral
WO2022011428A1 (fr) * 2020-07-16 2022-01-20 Griffith University Vaccin contre un virus atténué vivant

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11278617B2 (en) 2020-03-31 2022-03-22 Engimata, Inc Immunogenic composition forming a SARS-CoV-2 vaccine, and a method for its manufacture
US11559577B2 (en) 2020-03-31 2023-01-24 Engimata, Inc. Immunogenic composition forming a vaccine, and a method for its manufacture

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7561973B1 (en) * 2008-07-31 2009-07-14 Dna Twopointo, Inc. Methods for determining properties that affect an expression property value of polynucleotides in an expression system
WO2016077123A1 (fr) * 2014-11-10 2016-05-19 Moderna Therapeutics, Inc. Optimisation d'acides nucléiques à plusieurs paramètres

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
COLEMAN JR ET AL.: "Virus attenuation by genome-scale changes in codon pair bias", SCIENCE, vol. 320, no. 5884, 2008, pages 1784 - 1787, XP002633528, Retrieved from the Internet <URL:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2754401> [retrieved on 20080627] *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2020077395A1 (fr) * 2018-10-16 2020-04-23 Griffith University Vaccin antiviral
US20220088169A1 (en) * 2018-10-16 2022-03-24 Griffith University Virus vaccine
WO2022011428A1 (fr) * 2020-07-16 2022-01-20 Griffith University Vaccin contre un virus atténué vivant

Also Published As

Publication number Publication date
US20190359990A1 (en) 2019-11-28
EP3574097A4 (fr) 2020-11-18
EP3574097A1 (fr) 2019-12-04

Similar Documents

Publication Publication Date Title
WO2018138727A1 (fr) Séquences d&#39;acide nucléique synthétiques virales et leur utilisation
Gribble et al. The coronavirus proofreading exoribonuclease mediates extensive viral recombination
Trobaugh et al. RNA viruses can hijack vertebrate microRNAs to suppress innate immunity
Stern et al. The evolutionary pathway to virulence of an RNA virus
Miller et al. Evolutionary dynamics of Newcastle disease virus
JP5409354B2 (ja) 配列モチーフを同定するための方法、およびその応用
Neverov et al. Coordinated evolution of influenza A surface proteins
Kashiwagi et al. Ongoing phenotypic and genomic changes in experimental coevolution of RNA bacteriophage Qβ and Escherichia coli
Stern et al. Costs and benefits of mutational robustness in RNA viruses
Eschke et al. Attenuation of a very virulent Marek's disease herpesvirus (MDV) by codon pair bias deoptimization
Sellis et al. Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes
ES2749733T3 (es) Atenuación de virus de ARN mediante alteración de la robustez mutacional y el espacio de secuencia
Petterson et al. Natural infection of Atlantic salmon (Salmo salar L.) with salmonid alphavirus 3 generates numerous viral deletion mutants
Kerr et al. Genomic and phenotypic characterization of myxoma virus from Great Britain reveals multiple evolutionary pathways distinct from those in Australia
Warmbrod et al. Viral RNA-dependent RNA polymerase mutants display an altered mutation spectrum resulting in attenuation in both mosquito and vertebrate hosts
US20190060262A1 (en) Enhanced expression of rna vectors
Wu et al. Biased Pol II fidelity contributes to conservation of functional domains in the Potato spindle tuber viroid genome
Yang et al. The mitochondrial genome of the quiet-calling katydids, Xizicus fascipes (Orthoptera: Tettigoniidae: Meconematinae)
Mauroy et al. High-throughput sequencing analysis reveals the genetic diversity of different regions of the murine norovirus genome during in vitro replication
Zhou et al. Characteristics of codon usage bias in two regions downstream of the initiation codons of foot-and-mouth disease virus
Tyagi et al. Genome dynamics, codon usage patterns and influencing factors in Aeromonas hydrophila phages
Kumar et al. Evidence of reassortment between eastern and western topotype strains of bluetongue virus serotype 16 (BTV-16) from India
Kunec et al. Synthetically recoded virus sCPD9–A tool to accelerate SARS-CoV-2 research under biosafety level 2 conditions
US20220098598A1 (en) Methods for modifying the growth rate of a cell
Ma et al. Overall codon usage pattern of enterovirus 71

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18745104

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2018745104

Country of ref document: EP

Effective date: 20190826