1928z-Knock-In

trac-anti-cd19-donor.py

#!/usr/bin/python
import sys
#DNA PARTS
LITR = 'CCTGCAGGCAGCTGCGCGCTCGCTCGCTCACTGAGGCCGCCCGGGCAAAGCCCGGGCGTCGGGCGACCTTTGGTCGCCCGGCCTCAGTGAGCGAGCGAGCGCGCAGAGAGGGAGTGGCCAACTCCATCACTAGGGGTTCCT'

RITR = 'AGGAACCCCTAGTGATGGAGTTGGCCACTCCCTCTCTGCGCGCTCGCTCGCTCACTGAGGCCGGGCGACCAAAGGTCGCCCGACGCCCGGGCTTTGCCCGGGCGGCCTCAGTGAGCGAGCGAGCGCGCAGCTGCCTGCAGG'

POLYA = 'ACGGGTGGCATCCCTGTGACCCCTCCCCAGTGCCTCTCCTGGCCCTGGAAGTTGCCACTCCAGTGCCCACCAGCCTTGTCCTAATAAAATTAAGTTGCATCATTTTGTCTGACTAGGTGTCCTTCTATAATATTATGGGGTGGAGGGGGGTGGTATGGAGCAAGGGGCAAGTTGGGAAGACAACCTGTAGGGCCTGCGGGGTCTATTGGGAACCAAGCTGGAGTGCAGTGGCACAATCTTGGCTCACTGCAATCTCCGCCTCCTGGGTTCAAGCGATTCTCCTGCCTCAGCCTCCCGAGTTGTTGGGATTCCAGGCATGCATGACCAGGCTCAGCTAATTTTTGTTTTTTTGGTAGAGACGGGGTTTCACCATATTGGCCAGGCTGGTCTCCAACTCCTAATCTCAGGTGATCTACCCACCTTGGCCTCCCAAATTGCTGGGATTACAGGCGTGAACCACTGCTCCCTTCCCTGTCCTT'

CAR = 'ccatggctctcccagtgactgccctactgcttcccctagcgcttctcctgcatgcagaggtgaagctgcagcagtctggggctgagctggtgaggcctgggtcctcagtgaagatttcctgcaaggcttctggctatgcattcagtagctactggatgaactgggtgaagcagaggcctggacagggtcttgagtggattggacagatttatcctggagatggtgatactaactacaatggaaagttcaagggtcaagccacactgactgcagacaaatcctccagcacagcctacagcagctcagcggcctaacatctgaggactctgcggtctatttctgtgcaagaaagaccattagttcggtagtagatttctactttgacactggggccaagggaccacggtcaccgtctcctcaggtggaggtggatcaggtggaggtggatctggtggaggtggatctgacattgagctcacccagtctccaaaattcatgtccacatcagtaggagacagggtcagcgtcacctgcaaggccagtcagaatgtgggtactaatgtagcctggtatcaacagaaaccaggacaatctcctaaaccactgatttactcggcaacctaccggaacagtggagtccctgategcttcacaggcagtggatctgggacagatttcactctcaccatcactaacgtgcagtctaaagacttggcagactatttctgtcaacaatataacaggttccgtacacgtccggaggggggaccaagctggagatcaaacgggcggccgcaattgaagttatgtatcctcctccttacctagacaatgagaagagcaatggaaccaatccatgtgaaagggaaacacctttgtccaagtcccctatttcccggaccttctaagcccttttgggtgctggtggtggttggtggagtcctggcttgctatagcttgctagtaacagtggcctttattattttctgggtgaggagtaagaggagcaggctcctgcacagtgactacatgaacatgactccccgccgccccgggcccacccgcaagcattaccagccctatgccccaccacgcgacttcgcagcctatcgctccagagtgaagttcagcaggagcgcagagccccccgcgtaccagcagggccagaaccagctctataacgagctcaatctaggacgaagagaggagtacgatgttttggacaagagacgtggccgggaccctgagatggggggaaagccgagaaggaagaaccctcaggaaggcctgtacaatgaactgcagaaagataagatggcggaggcctacagtgagattgggatgaaaggcgagcgccggaggggcaaggggcacgatggcctttaccagggtctcagtacagccaccaaggacacctacgacgcccttcacatgcaggccctgccccctcgcg'.upper()

# Extracted from the UI, but easier to query DB in practice
TRAC_EXON_1 = 'TCCAGAACCCTGACCCTGCCGTGTACCAGCTGAGAGACTCTAAA'
SA = 'CAGAT'
PA = 'GGCGCCACCAACTTCTCCCTGCTGAAGCAGGCCGGCGACGTGGAGGAGAACCCCGGCCCC'
BGHPA_TERM = 'ATAAAAGATCTTTATTTTCATTAGATCTGTGTGTTGGTTTTTTGTGTG'
# to design
LHA_STOP = 'TAACCCTGATCCTCTTGT'
RHA_START = 'ATCCAGAAC'
PDHG = 'CCCA'

PAM = 'GGG'
PROTOSPACER = 'TCAGGGTTCTGGATATCTGT'
# 1900 bp sequence flanking start of target exon
DONOR = 'TGCTAATCCTCCGGCAAACCTCTGTTTCCTCCTCAAAAGGCAGGAGGTCGGAAAGAATAAACAATGAGAGTCACATTAAAAACACAAAATCCTACGGAAATACTGAAGAATGAGTCTCAGCACTAAGGAAAAGCCTCCAGCAGCTCCTGCTTTCTGAGGGTGAAGGATAGACGCTGTGGCTCTGCATGACTCACTAGCACTCTATCACGGCCATATTCTGGCAGGGTCAGTGGCTCCAACTAACATTTGTTTGGTACTTTACAGTTTATTAAATAGATGTTTATATGGAGAAGCTCTCATTTCTTTCTCAGAAGAGCCTGGCTAGGAAGGTGGATGAGGCACCATATTCATTTTGCAGGTGAAATTCCTGAGATGTAAGGAGCTGCTGTGACTTGCTCAAGGCCTTATATCGAGTAAACGGTAGTGCTGGGGCTTAGACGCAGGTGTTCTGATTTATAGTTCAAAACCTCTATCAATGAGAGAGCAATCTCCTGGTAATGTGATAGATTTCCCAACTTAATGCCAACATACCATAAACCTCCCATTCTGCTAATGCCCAGCCTAAGTTGGGGAGACCACTCCAGATTCCAAGATGTACAGTTTGCTTTGCTGGGCCTTTTTCCCATGCCTGCCTTTACTCTGCCAGAGTTATATTGCTGGGGTTTTGAAGAAGATCCTATTAAATAAAAGAATAAGCAGTATTATTAAGTAGCCCTGCATTTCAGGTTTCCTTGAGTGGCAGGCCAGGCCTGGCCGTGAACGTTCACTGAAATCATGGCCTCTTGGCCAAGATTGATAGCTTGTGCCTGTCCCTGAGTCCCAGTCCATCACGAGCAGCTGGTTTCTAAGATGCTATTTCCCGTATAAAGCATGAGACCGTGACTTGCCAGCCCCACAGAGCCCCGCCCTTGTCCATCACTGGCATCTGGACTCCAGCCTGGGTTGGGGCAAAGAGGGAAATGAGATCATGTCCTAACCCTGATCCTCTTGTCCCACAGATATCCAGAACCCTGACCCTGCCGTGTACCAGCTGAGAGACTCTAAATCCAGTGACAAGTCTGTCTGCCTATTCACCGATTTTGATTCTCAAACAAATGTGTCACAAAGTAAGGATTCTGATGTGTATATCACAGACAAAACTGTGCTAGACATGAGGTCTATGGACTTCAAGAGCAACAGTGCTGTGGCCTGGAGCAACAAATCTGACTTTGCATGTGCAAACGCCTTCAACAACAGCATTATTCCAGAAGACACCTTCTTCCCCAGCCCAGGTAAGGGCAGCTTTGGTGCCTTCGCAGGCTGTTTCCTTGCTTCAGGAATGGCCAGGTTCTGCCCAGAGCTCTGGTCAATGATGTCTAAAACTCCTCTGATTGGTGGTCTCGGCCTTATCCATTGCCACCAAAACCCTCTTTTTACTAAGAAACAGTGAGCCTTGTTCTGGCAGTCCAGAGAATGACACGGGAAAAAAGCAGATGAAGAGAAGGTGGCAGGAGAGGGCACGTGGCCCAGCCTCAGTCTCTCCAACTGAGTTCCTGCCTGCCTGCCTTTGCTCAGACTGTTTGCCCCTTACTGCTCTTCTAGGCCTCATTCTAAGCCCCTTCTCCAAGTTGCCTCTCCTTATTTCTCCCTGTCTGCCAAAAAATCTTTCCCAGCTCACTAAGTCAGTCTCACGCAGTCACTCATTAACCCACCAATCACTGATTGTGCCGGCACATGAATGCACCAGGTGTTGAAGTGGAGGAATTAAAAAGTCAGATGAGGGGTGTGCCCAGAGGAAGCACCATTCTAGTTGGGGGAGCCCATCTGTCAGCTGGGAAAAGTCCAAATAACTTCAGATTGGAATGTGTTTTAACTCAGGGTTGAGAAAACAGCTACCTTCAGGACAAAAGTCAGGGAAGGGCTCTCTGAAGAAATGCTACTTGAAGATACCAGCCCTACCAAGGGCAGGGAGAGGACCCTATAGAGGCCTGGGACAGGAGCTCAATGAGAAAGGAGAAGAGCAGCAGGCA'

if __name__ == '__main__':

    index = DONOR.find(LHA_STOP)
    LHA = DONOR[0:index+len(LHA_STOP)]
    RHA = DONOR[DONOR.find(RHA_START):]
    payload = ''.join([SA,PDHG,PA,CAR,BGHPA_TERM])
    design = ''.join([LITR,LHA,PDHG,payload,RHA,RITR])
    # Generate record
    sys.stdout.write(design)
    sys.exit()