Constructed Method for and Application of Nucleic Acid Multimerization-Mediated Multivalent Protein Drug and Vaccine

Information

  • Patent Application
  • 20240294674
  • Publication Number
    20240294674
  • Date Filed
    November 25, 2021
    3 years ago
  • Date Published
    September 05, 2024
    5 months ago
Abstract
A construction method for and an application of a nucleic acid multimerization-mediated multivalent protein drug and vaccine. Specifically provided is a multimeric complex based on complementary nucleic acid backbones. The complex is a multimer formed by complexing of 3-6 monomers having complementary nucleic acid backbones, wherein each monomer is a polypeptide having a nucleic acid single strand. In the multimer, the nucleic acid single strand of each monomer and the nucleic acid single strands of the other two monomers form double strands by means of base complementation, so as to form complementary nucleic acid backbone structures. Also provided are a pharmaceutical composition containing the multimeric complex, a nucleic acid sequence library used for constructing the multimeric complex, and a method for optimizing complementary nucleic acid backbones. By means of the method, off-the-shelf short-acting protein drugs or antigens can be used to complete multivalent formation of protein drugs or antigens without the need of reconstruction of fusion proteins or chemical modification and cross-linking, thereby improving their half-life and activity, and/or immunogenicity.
Description
TECHNICAL FIELD

The present invention relates to the field of biotechnology drugs, and specifically relates to a construction method for and an application of a nucleic acid multimerization-mediated multivalent protein drug and vaccine.


BACKGROUND

For many biological macromolecules, their aggregation or multivalent state directly affects their activities and half-lives in vivo. For example, the activation of most immune receptors involves the aggregation of receptors on the cell membrane, thereby activating downstream signaling pathways within the cell.


Therefore, the ability of natural ligands or antibodies of these receptors to activate receptors can often be significantly improved when they form multivalent or high-valent forms. In addition, some protein drugs with lower molecular weights (MW 40 kDa), such as cytokines, growth hormones, synthetic peptides, etc., have high renal clearance efficiencies and short half-lives in vivo; the molecular weights and half-lives of these protein drugs can also be increased by forming high-valent forms.


Therefore, multivalent formation of proteins is a highly concerned process in the field of biomedicine, and there are many existing methods. However, most chemical crosslinking methods have poor connection specificity and uneven connection of multimers. The most widely used method currently is to express and produce proteins in the form of multivalent fusion proteins by cells, which involves fusing drug functional protein regions with proteins that can form oligomers to form chimeras, such as Fc divalent fusion proteins and GCN4 trivalent fusion proteins, etc. These fusion proteins can form even oligomers, but the cellular expression and activity of fusion proteins are often worse than those of original protein drugs, and the presence of Fc regions brings about the activation of the immune system, induction of cytokine release, and risk of cytotoxicity. Therefore, there is an urgent need to develop a simple, flexible, and efficient method in this field that can form even and highly specific multivalent proteins from validated protein drugs in a non-fusion protein manner.


In addition, multivalent formation of proteins is also of great significance for vaccine development. Firstly, in the design of B cell-based vaccines, activating B cell receptors (BCRs) using viral or bacterial proteins as antigens is a crucial step. Like the immune receptors mentioned above, the effective activation of BCR requires the aggregation of receptors on the cell membrane, so high-valent antigens have an absolute advantage over monomeric antigens in activating B cells. Secondly, high-valent antigens may not necessarily form oligomers from a single antigen; high-valent antigens can contain different proteins in a certain virus or mutations and subtypes of the same protein in different virus strains; in this way, diversified antigen presentation can theoretically induce polyclonal response of B cells in the host immune system, producing a wider range of neutralizing antibodies.


SUMMARY OF THE INVENTION

One purpose of the present invention is to provide an efficient and stable assembly backbone design for n-order nucleic acid oligomers, suitable for the efficient and stable assembly of nucleic acid coupled protein drugs to form multivalent drugs or vaccines.


The second purpose of the present invention is to provide a simple and efficient method for forming multivalent macromolecules from protein drugs for extending the half-life of drugs and increasing drug activity.


The third purpose of the present invention is to provide a simple, flexible, efficient, and modular method for forming multivalent macromolecular complexes of the same or different protein antigens for activating immune cells and improving the immunogenicity of vaccines.


In the first aspect of the present invention, it provides a multimeric complex based on a complementary nucleic acid backbone, wherein the complex is a multimer formed by complexing n monomers having the complementary nucleic acid backbone, wherein each monomer is a polypeptide having a nucleic acid single strand, and n is a positive integer of 3-6; in the multimer, the nucleic acid single strand of each monomer and the nucleic acid single strands of the other two monomers form complementary double strands by means of base complementation, so as to form complementary nucleic acid backbone structures.


In another preferred embodiment, n is 3, 4, 5, or 6.


In another preferred embodiment, the complex is a trimer, tetramer, or pentamer, preferably with a structure as shown in FIG. 1.


In another preferred embodiment, the monomer has a structure of formula I:





Z1-W   (I)

    • wherein,
    • Z1 is a polypeptide moiety;
    • W is a nucleic acid single strand sequence; and
    • “-” is a linker or bond.


In another preferred embodiment, “-” is a covalent bond.


In another preferred embodiment, the nucleic acid sequence is selected from the group consisting of: left-handed nucleic acid, peptide nucleic acid, locked nucleic acid, thio-modified nucleic acid, 2′-fluoro modified nucleic acid, 5-hydroxymethylcytosine nucleic acid, phosphorodiamidate morpholino nucleic acid, and combinations thereof;


In another preferred embodiment, in the multimer, the Z1 of each monomer is the same or different.


In another preferred embodiment, in the multimer, the W of each monomer is different.


In another preferred embodiment, the monomer has a structure of formula II:





D-[L-W]m   (II)

    • wherein,
    • D is a protein drug element moiety;
    • W is a nucleic acid sequence;
    • L is not present or a linker;
    • “-” is a covalent bond; and
    • m is 1, 2, or 3.


In another preferred embodiment, m is 1.


In another preferred embodiment, the monomer has a structure of formula III:





A-[L-W]m   (III)

    • wherein,
    • A is a peptide antigen element moiety;
    • W is a nucleic acid sequence;
    • L is not present or a linker;
    • “-” is a covalent bond; and
    • m is 1, 2, or 3.


In another preferred embodiment, m is 1.


In another preferred embodiment, the nucleic acid sequence W has the structure shown in formula 1:





X1-R1-X2-R2-X3   (1)

    • wherein,
    • R1 is a complementary base pairing region 1;
    • R2 is a complementary base pairing region 2;
    • Each of X1, X2, and X3 is independently not present or redundant nucleic acids; and
    • “-” is a bond.


In another preferred embodiment, each of R1 and R2 is independently 10-20 bases, preferably 14-16 bases in length.


In another preferred embodiment, X1 is 0-5 bases in length.


In another preferred embodiment, X3 is 0-5 bases in length.


In another preferred embodiment, X2 is 0-3 bases in length.


In another preferred embodiment, the sequence of X2 is selected from the group consisting of: A, AA, AGA and AAA.


In another preferred embodiment, the R1 of each monomer forms a complementary base pairing structure with the R2 of the left neighbor (or left side) monomer; while the R2 forms a complementary base pairing structure with the R1 of the right neighbor (or right side) monomer.


In another preferred embodiment, the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 1-60 (see Table 9-1) that form a trimer complex based on the complementary nucleic acid backbone.


In another preferred embodiment, the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 61-140 (see Table 9-2) that form a tetramer complex based on the complementary nucleic acid backbone.


In another preferred embodiment, the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 141-240 (see Table 9-3) that forms a pentamer complex based on the complementary nucleic acid backbone.


In another preferred embodiment, the monomer sequence is a phosphorodiamidate morpholino nucleic acid.


In another preferred embodiment, the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 275-278 that forms a tetramer complex based on the complementary nucleic acid backbone.


In the second aspect of the present invention, it provides a pharmaceutical composition comprising:

    • (a) the multimeric complex based on the complementary nucleic acid backbone according to the first aspect; and
    • (b) a pharmaceutically acceptable carrier.


In another preferred embodiment, the pharmaceutical composition comprises a vaccine composition.


In another preferred embodiment, the pharmaceutical composition comprises a therapeutic and/or prophylatic pharmaceutical composition.


In another preferred embodiment, the multimeric complex comprises a trimer complex, a tetramer complex, and a pentamer complex.


In the third aspect of the present invention, it provides a nucleic acid sequence library, which comprises a nucleic acid sequence for forming the multimeric complex based on the complementary nucleic acid backbone according to the first aspect.


In another preferred embodiment, the nucleic acid sequence comprises:

    • (a) a nucleic acid sequence for forming a trimer complex based on the complementary nucleic acid backbone;
    • (b) a nucleic acid sequence for forming a tetramer complex based on the complementary nucleic acid backbone; and/or
    • (c) a nucleic acid sequence for forming a pentamer complex based on the complementary nucleic acid backbone.


In another preferred embodiment, the nucleic acid sequence W has the structure shown in formula 1:





X1-R1-X2-R2-X3   (1)

    • wherein,
    • R1 is the complementary base pairing region 1;
    • R2 is the complementary base pairing region 2;
    • Each of X1, X2, and X3 is independently not present or redundant nucleic acids; and
    • “-” is a bond.


In the fourth aspect of the present invention, it provides use of the nucleic acid sequence library according to the third aspect in the manufacture of the multimeric complex according to the first aspect or a pharmaceutical composition comprising the multimeric complex.


In the fifth aspect of the present invention, it provides a method of determining a nucleic acid single strand sequence for forming a multimeric complex based on a complementary nucleic acid backbone, comprising steps of:

    • (a) setting annealing algorithm parameters:
    • setting the initial annealing temperature, annealing termination temperature, and annealing temperature attenuation coefficient ΔT;
    • setting optimized constraint parameters:
    • {circle around (1)} the number n of the nucleic acid single strand, preferably a positive integer of 3-6;
    • {circle around (2)} the length L of the pairing sequence, preferably the L is of 12-16 bases;
    • {circle around (3)} the dissociation temperature threshold Tm of the pairing region;
    • {circle around (4)} the free energy threshold ΔG°S of the specific pairing region sequence;
    • {circle around (5)} the free energy threshold ΔG°NS of the non-specific pairing;
    • {circle around (6)} the connecting element X2, preferably A, AA, and AAA;
    • {circle around (7)} the dissociation temperature threshold Tm-H of the secondary structure (hairpin);
    • {circle around (8)} the CG proportion PCG in the pairing sequence, preferably the range of PCG is [0.4,0.6);


{circle around (9)} optionally, for n=4, using a symmetric sequence to initialize a sequence set S={S1, S2, . . . , Sn} according to the above parameters;

    • (b) calculating the objective function value E0 of the set S of the previous step, that is, calculating the sum of the non-specific pairing free energies (ΔG°NS) between sequences and of the sequence itself, while obtaining the non-specific pairing free energy matrix Cn×n, searching the Si and Sj(1≤i≤n, 1≤j≤n) corresponding to the minimum value in the upper triangular matrix thereof, randomly selecting Si or Sj for an updated operation according to the non-specific pairing free energy of the Si and Sj ΔG°NS(Si, Sj), and then obtaining a new nucleic acid sequence, thereby obtaining a updated sequence set S′;
    • (c) determining whether the sequences in the set S′ of the previous step meet the optimized constraint parameter conditions set in step (a), verifying the following parameters, including the dissociation temperature Tm of the specific pairing region, the free energy ΔG°S of the specific pairing region sequence, the dissociation temperature Tm-H of the secondary structure and the CG proportion PCG. If the above parameters meet the constraint conditions, the step (d) is proceeded; otherwise, the step (c) is repeated. If the step (b) is performed 15 times continuously at a certain annealing temperature without obtaining the S′ that meets the conditions, then the set S becomes the set S′ and the next step is proceeded to prevent a dead cycle;
    • (d) calculating the objective function value E1 of the set S′ of the previous step, and comparing E0 with E1. If E1≥E0, it indicates that the non-specific pairing free energy has been optimized, and the sequence set S′ becomes the sequence set S. If E1<E0, it indicates that the non-specific pairing free energy has not been optimized, and in this case, it is necessary to determine whether to accept the set S′ as S according to the Metropolis criterion; and
    • (e) the annealing temperature is attenuated according to the attenuation coefficient ΔT set in the step (a), and the steps (b), (c), and (d) are repeated for the S of the previous step, which is the Monte Carlo-based annealing algorithm, until the annealing temperature reaches the annealing termination temperature. The S={S1, S2, . . . , Sn} of the previous step becomes the nucleic acid single strand sequence for forming the multimeric complex based on the complementary nucleic acid backbone.


In another preferred embodiment, in step (a) setting annealing algorithm parameters, it comprises:

    • for example, setting the initial annealing temperature T0=50° C.±2° C., the annealing termination temperature Tf=0.12° C.±0.02° C., and the annealing temperature attenuation coefficient ΔT depends on the situation, usually 0.98±0.01;
    • setting optimized constraint parameters:
    • {circle around (1)} the number n of the nucleic acid single strand is a positive integer of 3-6,
    • {circle around (2)} the length L of the pairing sequence depends on the situation (preferably the L is of 12-16 bases),
    • {circle around (3)} the dissociation temperature threshold Tm of the pairing region is determined by the length of the pairing sequence (e.g., when L=14 bases, Tm>50° C.; when L=16 bases, Tm>52° C.),
    • {circle around (4)} the free energy threshold ΔG°S of the specific pairing region sequence is determined by the length of the pairing sequence (preferably, when L=14 bases, ΔG°S<−27 kcal/mol; when L=16 bases, ΔG°S<−29 kcal/mol),
    • {circle around (5)} the free energy threshold ΔG°NS of the non-specific pairing is determined by the sequence length (preferably, ΔG°NS>−7 kcal/mol),
    • {circle around (6)} the connecting element X2 depends on the situation (can be A, AA, and AAA, etc.),
    • {circle around (7)} the dissociation temperature threshold Tm-H of the secondary structure (hairpin) depends on the situation (preferably, Tm-H<40° C.±2° C.),
    • {circle around (8)} the range of the CG proportion PCG in the pairing sequence is [0.4,0.6),
    • {circle around (9)} specifically, for n=4, a symmetric sequence can be used to initialize the sequence S={S1, S2, . . . , Sn} according to the above parameters.


In another preferred embodiment, each nucleic acid single strand sequence W has the structure shown in formula 1:





X1-R1-X2-R2-X3   (1)

    • wherein,
    • R1 is the complementary base pairing region 1;
    • R2 is the complementary base pairing region 2;
    • Each of X1, X2, and X3 independently is not present or redundant nucleic acids; and
    • “-” is a bond.


In another preferred embodiment, in the step (d), the optimized set is a set that satisfies following conditions:

    • (C1) the free energy (ΔG°S) of the DNA double strand structure formed by the target pairing is smaller or smallest in the complementary nucleic acid backbone structure; and
    • (C2) the ΔG°S of the non-target pairing is larger or the largest in the complementary nucleic acid backbone structure.


In another preferred embodiment, in the step (d), the optimized set also satisfies following conditions:

    • (C3) the pairing dissociation temperature of the R1 and R2 regions Tm>50° C. (when L=14 bases).


In another preferred embodiment, in the step (c), the free energy (ΔG°S) of the DNA oligomer (i.e., the complementary nucleic acid backbone structure) is calculated using the nearest neighbor method.


In another preferred embodiment, in the step (c), the DNA oligomer (i.e., the complementary nucleic acid backbone structure) is decomposed into 10 different nearest neighbor pairing interactions, which are: AA/TT; AT/TA; TA/AT; CA/GT; GT/CA; CT/GA; GA/CT; CG/GC; GC/CG; and GG/CC; and the corresponding ΔG° value is calculated and obtained respectively based on the enthalpy (ΔH°) and entropy (ΔS°) of these pairing interactions; then, the free energies of the pairing interactions included in the complementary nucleic acid backbone structure are merged (or summed) to obtain the free energy of the complementary nucleic acid backbone structure.


In another preferred embodiment, the method comprises repeating the steps (b), (c), and (d) for multiple times (i.e., performing n1 iterations) to obtain the global optimal solution during the iteration process.


In another preferred embodiment, during the iteration process, a poor solution is limitedly accepted according to the Metropolis criterion, and the probability of accepting the poor solution is gradually approaching zero, so as to find the global optimal solution at all possible when the algorithm terminates.


In another preferred embodiment, the following optimized objective function is used for the iteration of the simulated annealing algorithm to optimize the free energy of the non-target pairing region:







E
=






i
=
1




n








j
=
i




n




-
Δ




G


(


S
i

,

S
j


)





,

n

1







    • ΔG°(Si, Sj) is the free energy of the non-target pairing between Si sequence and Sj sequence, and Σi=1nΣj=inΔG°(Si, Sj) is the sum of the free energies of non-target pairing between all sequences, with a negative value; wherein the larger the negative value, the more beneficial it is to reduce the non-target pairing.





In the sixth aspect of the present invention, it provides a nucleic acid single strand sequence set for forming a multimeric complex based on a complementary nucleic acid backbone, which is determined using the method of the fifth aspect.


In another preferred embodiment, the set is selected from the group consisting of:

    • (S1) a nucleic acid single strand sequence for forming a trimer complex based on the complementary nucleic acid backbone:











TABLE 9-1







SEQ




ID




NO:

















Sequence set 3-1




numbering
Optimized sequence





S1

ACACCTGGTTGTTGGATAAATCGTTGAAG

1



GCTAGGA






S2

ATCCTAGCCTTCAACGAAAAAACTAGAGT

2



CCGCCGA






S3

ATCGGCGGACTCTAGTTAAAATCCAACAA

3



CCAGGTG





Sequence set 3-2




numbering
Optimized sequence





S1

ATGCGTTGAGTTCCAGTAAAGGCAACATC

4



ACCACAT






S2

AATGTGGTGATGTTGCCAAATCTGAATCC

5



TCGTGCT






S3

AAGCACGAGGATTCAGAAAAACTGGAAC

6



TCAACGCA





Sequence set 3-3




numbering
Optimized sequence





S1

ATTCCAATCGTCCTGTGAAAAGTTCCGCT

7



CTGAGTT






S2

AAACTCAGAGCGGAACTAAACTGGCAGA

8



TGGATGAA






S3

ATTCATCCATCTGCCAGAAACACAGGACG

9



ATTGGAA





Sequence set 3-4




numbering
Optimized sequence





S1

ACGAGGCAAGTTCTGTGAAAATGACTACC

10



AGGTCCG






S2

ACGGACCTGGTAGTCATAAAATCCACTGA

11



CGCTGAA






S3

ATTCAGCGTCAGTGGATAAACACAGAACT

12



TGCCTCG





Sequence set 3-5




numbering
Optimized sequence





S1

ATAGTTCGTTGCTCGGAAAAGGCATTGAG

13



AGGACCT






S2

AAGGTCCTCTCAATGCCAAAATGGTGATG

14



TCGCTTG






S3

ACAAGCGACATCACCATAAATCCGAGCAA

15



CGAACTA





Sequence set 3-6




numbering
Optimized sequence





S1

AGTCGTGTGCTTCCAAGAAATAGCCAGGT

16



GAGGACT






S2

AAGTCCTCACCTGGCTAAAAAACAGCGGA

17



GTGTCAT






S3

AATGACACTCCGCTGTTAAACTTGGAAGC

18



ACACGAC





Sequence set 3-7




numbering
Optimized sequence





S1

AACGCATCGCTTGATAGAAAAGAGGAGC

19



ACGGTTAT






S2

AATAACCGTGCTCCTCTAAAGTAGGCAAT

20



CCACCAT






S3

AATGGTGGATTGCCTACAAACTATCAAGC

21



GATGCGT





Sequence set 3-8




numbering
Optimized sequence





S1

AGTCGTTCCACCGAACAAAATGGCTCTGG

22



TCATTGA






S2

ATCAATGACCAGAGCCAAAAAATCGCAC

23



ATCTCAGG






S3

ACCTGAGATGTGCGATTAAATGTTCGGTG

24



GAACGAC





Sequence set 3-9




numbering
Optimized sequence





S1

AGCGGAGTGACCATAGTAAAAGGCAGGA

25



CATTGTTC






S2

AGAACAATGTCCTGCCTAAAGTGCTCGTC

26



GTGAAGA






S3

ATCTTCACGACGAGCACAAAACTATGGTC

27



ACTCCGC





Sequence set 3-




10 numbering
Optimized sequence





S1

AATTGGACCGCTCTACTAAAATGGCACCA

28



CAGTCAA






S2

ATTGACTGTGGTGCCATAAACAGGCTATC

29



AGCATCC






S3

AGGATGCTGATAGCCTGAAAAGTAGAGC

30



GGTCCAAT





Sequence set 3-




11 numbering
Optimized sequence





S1

ACCATTGAGCCAGTGATAAAAACCGTTGT

31



GAGTTGC






S2

AGCAACTCACAACGGTTAAATCGCACACC

32



TGTCGTA






S3

ATACGACAGGTGTGCGAAAAATCACTGGC

33



TCAATGG





Sequence set 3-




12 numbering
Optimized sequence





S1

AAGTGAAGAAGCAGCCTAAAGTTGTCATC

34



GCACACC






S2

AGGTGTGCGATGACAACAAAATGTCGTAA

35



CCGTGGA






S3

ATCCACGGTTACGACATAAAAGGCTGCTT

36



CTTCACT





Sequence set 3-




13 numbering
Optimized sequence





S1

AATAGCGTCTTGAGCCTAAATGGAGGACA

37



TACCGAC






S2

AGTCGGTATGTCCTCCAAAAGGTCACAGT

38



TGCTGCT






S3

AAGCAGCAACTGTGACCAAAAGGCTCAA

39



GACGCTAT





Sequence set 3-




14 numbering
Optimized sequence





S1

ATGCCGTGTTCAGATTCAAATGTGCGTCT

40



GGATTGA






S2

ATCAATCCAGACGCACAAAAAGACAGGT

41



GGTCCGAT






S3

AATCGGACCACCTGTCTAAAGAATCTGAA

42



CACGGCA





Sequence set 3-




15 numbering
Optimized sequence





S1

ATTCAGGACAGCGTCATAAAACCGACTGG

43



AGCAACT






S2

AAGTTGCTCCAGTCGGTAAAGATGCCTTC

44



GTGTGAG






S3

ACTCACACGAAGGCATCAAAATGACGCTG

45



TCCTGAA





Sequence set 3-




16 numbering
Optimized sequence





S1

AGCAGCCAAGGTTATCTAAACAATGACAC

46



GGAGGAT






S2

AATCCTCCGTGTCATTGAAAGTGATTCGC

47



ACCAGAC






S3

AGTCTGGTGCGAATCACAAAAGATAACCT

48



TGGCTGC





Sequence set 3-




17 numbering
Optimized sequence





S1

ACCACCGTGTATGACCTAAAAGTGACAGC

49



ACATCGC






S2

AGCGATGTGCTGTCACTAAAACAGGCTCT

50



ACGAGGA






S3

ATCCTCGTAGAGCCTGTAAAAGGTCATAC

51



ACGGTGG





Sequence set 3-




18 numbering
Optimized sequence





S1

AACTACGGAGCGAAGATAAATCCTGACCA

52



ACTTGCT






S2

AAGCAAGTTGGTCAGGAAAAGACTGGCT

53



GAACACGA






S3

ATCGTGTTCAGCCAGTCAAAATCTTCGCT

54



CCGTAGT





Sequence set 3-




19 numbering
Optimized sequence





S1

AGTTCCTGATCCAGCCTAAACATCCTTGTC

55



TTGCCA






S2

ATGGCAAGACAAGGATGAAACACGACCG

56



CTTAGAAG






S3

ACTTCTAAGCGGTCGTGAAAAGGCTGGAT

57



CAGGAAC





Sequence set 3-




20 numbering
Optimized sequence





S1

ATATCGCACTCCAGCATAAACCGTGTGAA

58



CATCAGG






S2

ACCTGATGTTCACACGGAAAAGCCTACGA

59



GACTTGG






S3

ACCAAGTCTCGTAGGCTAAAATGCTGGAG

60



TGCGATA











    • (S2) a nucleic acid single strand sequence for forming a tetramer complex based on the complementary nucleic acid backbone:














TABLE 9-2







SEQ




ID




NO:

















Sequence set 4-1




numbering
Optimized sequence





S1

AAGCGTCGTGAATCCAAATGAGCCTGC

61



CAATG






S2

ACATTGGCAGGCTCAAAACCGAAGTCA

62



ACGCT






S3

AAGCGTTGACTTCGGAAAACTATGGAC

63



GGCGA






S4

ATCGCCGTCCATAGTAAAGGATTCACG

64



ACGCT





Sequence set 4-2




numbering
Optimized sequence





S1

AATGGCGAGCAATCCAAATGAGCCTGG

65



ACCAA






S2

ATTGGTCCAGGCTCAAAACCGAACGCT

66



GTGAT






S3

AATCACAGCGTTCGGAAAACTATCGTG

67



CGGCA






S4

ATGCCGCACGATAGTAAAGGATTGCTC

68



GCCAT





Sequence set 4-3




numbering
Optimized sequence





S1

ATGACCACGCAATCCAAATGAGCCAAC

69



CTCCA






S2

ATGGAGGTTGGCTCAAAACCGAACAGC

70



AGCTT






S3

AAAGCTGCTGTTCGGAAAACTATCTGC

71



CGCCT






S4

AAGGCGGCAGATAGTAAAGGATTGCGT

72



GGTCA





Sequence set 4-4




numbering
Optimized sequence





S1

ATGTCGCACCAATCCAAATGAGCAAGC

73



CTCGT






S2

AACGAGGCTTGCTCAAAACCGAACGCT

74



GTCAT






S3

AATGACAGCGTTCGGAAAACTATGTGG

75



CGGCA






S4

ATGCCGCCACATAGTAAAGGATTGGTG

76



CGACA





Sequence set 4-5




numbering
Optimized sequence





S1

ATGCTGGCACAATCCAAATGAGCGACG

77



AGGTT






S2

AAACCTCGTCGCTCAAAACCGAAGTGC

78



CAGTT






S3

AAACTGGCACTTCGGAAAACTATGAGG

79



CGGCT






S4

AAGCCGCCTCATAGTAAAGGATTGTGC

80



CAGCA





Sequence set 4-6




numbering
Optimized sequence





S1

ATGTCGCACCAATCCAAATGAGCAGGT

81



TGGCA






S2

ATGCCAACCTGCTCAAAACCGAACGCT

82



GTCAA






S3

ATTGACAGCGTTCGGAAAACTATCAGC

83



CGCCT






S4

AAGGCGGCTGATAGTAAAGGATTGGTG

84



CGACA





Sequence set 4-7




numbering
Optimized sequence





S1

ATGTGGTCGCAATCCAAATGAGCACCT

85



GCCAA






S2

ATTGGCAGGTGCTCAAAACCGAACGTG

86



ACGAT






S3

AATCGTCACGTTCGGAAAACTATCAAC

87



GCCGC






S4

AGCGGCGTTGATAGTAAAGGATTGCGA

88



CCACA





Sequence set 4-8




numbering
Optimized sequence





S1

AAGCGTCGTCAATCCAAATGAGCACGG

89



CAATG






S2

ACATTGCCGTGCTCAAAACCGAAGTGA

90



ACGCT






S3

AAGCGTTCACTTCGGAAAACTATGGCT

91



CGCCT






S4

AAGGCGAGCCATAGTAAAGGATTGACG

92



ACGCT





Sequence set 4-9




numbering
Optimized sequence





S1

ATGTGGCGACAATCCAAATGAGCAAGC

93



CTCCA






S2

ATGGAGGCTTGCTCAAAACCGAAGACG

94



CTGTT






S3

AAACAGCGTCTTCGGAAAACTATCGTG

95



CGGCA






S4

ATGCCGCACGATAGTAAAGGATTGTCG

96



CCACA





Sequence set 4-




10 numbering
Optimized sequence





S1

ATGCTGCCACAATCCAAATGAGCCTGG

97



AACCA






S2

ATGGTTCCAGGCTCAAAACCGAACGCA

98



GTCAT






S3

AATGACTGCGTTCGGAAAACTATCGCC

99



GCTCT






S4

AAGAGCGGCGATAGTAAAGGATTGTGG

100



CAGCA





Sequence set 4-




11 numbering
Optimized sequence





S1

ATGCGTCGTCAATCCAAATGAGCTTGG

101



CAAGG






S2

ACCTTGCCAAGCTCAAAACCGAACGTG

102



CTGTT






S3

AAACAGCACGTTCGGAAAACTATGGAG

103



CGGCT






S4

AAGCCGCTCCATAGTAAAGGATTGACG

104



ACGCA





Sequence set 4-




12 numbering
Optimized sequence





S1

AACTGCCAGCAATCCAAATGAGCCTCG

105



TTCCA






S2

ATGGAACGAGGCTCAAAACCGAAGTTG

106



GCAGT






S3

AACTGCCAACTTCGGAAAACTATCGCC

107



GCTTG






S4

ACAAGCGGCGATAGTAAAGGATTGCTG

108



GCAGT





Sequence set 4-




13 numbering
Optimized sequence





S1

ATGCGTCGTCAATCCAAATGAGCCTCC

109



AGGTT






S2

AAACCTGGAGGCTCAAAACCGAATGAC

110



ACGCT






S3

AAGCGTGTCATTCGGAAAACTATGGCG

111



GCAGT






S4

AACTGCCGCCATAGTAAAGGATTGACG

112



ACGCA





Sequence set 4-




14 numbering
Optimized sequence





S1

AAGCGTCGTGAATCCAAATGAGCCATC

113



GTCCA






S2

ATGGACGATGGCTCAAAACCGAATGTG

114



CTGGT






S3

AACCAGCACATTCGGAAAACTATGCGG

115



CAACC






S4

AGGTTGCCGCATAGTAAAGGATTCACG

116



ACGCT





Sequence set 4-




15 numbering
Optimized sequence





S1

ATTGCCAGGATGCTGAATCACGGTCGG

117



ACA






S2

ATGTCCGACCGTGATAGTCGCAGAAGG

118



CAT






S3

AATGCCTTCTGCGACATAGTACAACGC

119



CGC






S4

AGCGGCGTTGTACTAACAGCATCCTGG

120



CAA





Sequence set 4-




16 numbering
Optimized sequence





S1

AGGCGATCACAATCCAAATGAGCGTGT

121



TACGG






S2

ACCGTAACACGCTCAAAACCGAAGTGC

122



CAATT






S3

AAATTGGCACTTCGGAAAACTATGCGG

123



CTGCT






S4

AAGCAGCCGCATAGTAAAGGATTGTGA

124



TCGCC





Sequence set 4-




17 numbering
Optimized sequence





S1

ATGGTCCAACACGCTAAGCCTCACCGT

125



CTT






S2

AAAGACGGTGAGGCTATCGCACAACCT

126



GGT






S3

AACCAGGTTGTGCGAATCGGAGTGGCA

127



GAA






S4

ATTCTGCCACTCCGAAAGCGTGTTGGAC

128



CA





Sequence set 4-




18 numbering
Optimized sequence





S1

AACCTTGGTGTGCGAAACTCCTGGCAG

129



CAA






S2

ATTGCTGCCAGGAGTAAGCGTGTGGTT

130



CCA






S3

ATGGAACCACACGCTATGAGGACCGTC

131



GTT






S4

AAACGACGGTCCTCAATCGCACACCAA

132



GGT





Sequence set 4-




19 numbering
Optimized sequence





S1

ATGCCAAGTCCGAGAATGCTGCGAACT

133



GGT






S2

AACCAGTTCGCAGCAAAGAGCCTGAAC

134



CGT






S3

AACGGTTCAGGCTCTAACGACGCTTGA

135



CCA






S4

ATGGTCAAGCGTCGTATCTCGGACTTGG

136



CA





Sequence set 4-




20 numbering
Optimized sequence





S1

AAGCAGCCTCGTTGAATCGCCAAGACA

137



CCT






S2

AAGGTGTCTTGGCGAAAGTTGCTCCGA

138



CGA






S3

ATCGTCGGAGCAACTAAGCGGTTCTGT

139



GGA






S4

ATCCACAGAACCGCTATCAACGAGGCT

140



GCT











    • (S3) a nucleic acid single strand sequence for forming a pentamer complex based on the complementary nucleic acid backbone:














TABLE 9-3







Sequence set 5-1

SEQ ID


numbering
Optimized sequence
NO:





S1

ATCAGGCGACCTCTTAAAACCACCATCGT

141



TGC






S2

AGCAACGATGGTGGTAAAAATCCAAATGA

142



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

143



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

144



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAA

145



GAGGTCGCCTGA





Sequence set 5-2




numbering
Optimized sequence





S1

AGGCGACGATGTCTTAAAACCTGGTTGCT

146



GGA






S2

ATCCAGCAACCAGGTAAAAATCCAAATGA

147



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

148



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

149



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAA

150



GACATCGTCGCC





Sequence set 5-3




numbering
Optimized sequence





S1

ATGGAACCTGGTGCTAAATGCTCGCCTGT

151



CAA






S2

ATTGACAGGCGAGCAAAAAATCCAAATG

152



AGCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

153



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

154



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAG

155



CACCAGGTTCCA





Sequence set 5-4




numbering
Optimized sequence





S1

ATGGTCAGGCGACTTAAAAGGACGAGGTT

156



GCT






S2

AAGCAACCTCGTCCTAAAAATCCAAATGA

157



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

158



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

159



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAA

160



GTCGCCTGACCA





Sequence set 5-5




numbering
Optimized sequence





S1

ATGCTGGACCACCTTAAATCAGATGGAGG

161



CGA






S2

ATCGCCTCCATCTGAAAAAATCCAAATGA

162



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

163



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

164



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAA

165



GGTGGTCCAGCA





Sequence set 5-6




numbering
Optimized sequence





S1

AAACGTCCAGGAGCTAAATCTCGTCGCCT

166



GAA






S2

ATTCAGGCGACGAGAAAAAATCCAAATG

167



AGCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

168



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

169



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAG

170



CTCCTGGACGTT





Sequence set 5-7




numbering
Optimized sequence





S1

ACCACGACCATTGCTAAAAACTTCAGGCG

171



ACG






S2

ACGTCGCCTGAAGTTAAAAATCCAAATGA

172



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

173



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

174



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAG

175



CAATGGTCGTGG





Sequence set 5-8




numbering
Optimized sequence





S1

AAGGCGAGGTCTTCAAAATGGTTGCTGGA

176



CGA






S2

ATCGTCCAGCAACCAAAAAATCCAAATGA

177



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

178



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

179



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAATGA

180



AGACCTCGCCT





Sequence set 5-9




numbering
Optimized sequence





S1

ATCAAGGCGACCAGTAAAAAGCTCCTCGA

181



CGA






S2

ATCGTCGAGGAGCTTAAAAATCCAAATGA

182



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

183



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

184



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAACT

185



GGTCGCCTTGA





Sequence set 5-




10 numbering
Optimized sequence





S1

ATTCAGGCGACTCCTAAAAGCACGACGAT

186



GGT






S2

AACCATCGTCGTGCTAAAAATCCAAATGA

187



GCGTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCA

188



ATT






S4

AAATTGGCACTTCGGAAAACTATGCGGCT

189



GCT






S5

AAGCAGCCGCATAGTAAAGGATTAAAAG

190



GAGTCGCCTGAA





Sequence set 5-




11 numbering
Optimized sequence





S1

AAGCACCTGCAATCCAAATCGCCAGGACA

191



AGT






S2

AACTTGTCCTGGCGAAAATGAGCAACCAT

192



GCC






S3

AGGCATGGTTGCTCAAAACCGAACGTCGT

193



GAT






S4

AATCACGACGTTCGGAAAACTATGGAGCG

194



GCT






S5

AAGCCGCTCCATAGTAAAGGATTGCAGGT

195



GCT





Sequence set 5-




12 numbering
Optimized sequence





S1

AACCTGCTGCAATCCAAATCGCCACCTCA

196



AGA






S2

ATCTTGAGGTGGCGAAAATGAGCCTGGAC

197



GTT






S3

AAACGTCCAGGCTCAAAACCGAACTGGTG

198



CTT






S4

AAAGCACCAGTTCGGAAAACTATGCCGCT

199



CCT






S5

AAGGAGCGGCATAGTAAAGGATTGCAGC

200



AGGT





Sequence set 5-




13 numbering
Optimized sequence





S1

AAGCTGGTGCAATCCAAATCGCCTCCTGA

201



CAA






S2

ATTGTCAGGAGGCGAAAATGAGCAAGGTT

202



GGC






S3

AGCCAACCTTGCTCAAAACCGAACGCAGA

203



TGT






S4

AACATCTGCGTTCGGAAAACTATGGAGCG

204



GCA






S5

ATGCCGCTCCATAGTAAAGGATTGCACCA

205



GCT





Sequence set 5-




14 numbering
Optimized sequence





S1

ATGCACGCACAATCCAAATCGCCATCAGA

206



GGT






S2

AACCTCTGATGGCGAAAATGAGCTGCCTC

207



CAT






S3

AATGGAGGCAGCTCAAAACCGAACGTCGT

208



CAT






S4

AATGACGACGTTCGGAAAACTATCGAGCG

209



GCT






S5

AAGCCGCTCGATAGTAAAGGATTGTGCGT

210



GCA





Sequence set 5-




15 numbering
Optimized sequence





S1

AAGCGTCGTGAATCCAAATCGCCATCAGA

211



CCA






S2

ATGGTCTGATGGCGAAAATGAGCAAGGCT

212



CGT






S3

AACGAGCCTTGCTCAAAACCGAACCAGCT

213



TGT






S4

AACAAGCTGGTTCGGAAAACTATGCGGCA

214



GGT






S5

AACCTGCCGCATAGTAAAGGATTCACGAC

215



GCT





Sequence set 5-




16 numbering
Optimized sequence





S1

ATCAGCACGCAATCCAAATCGCCAGTTCA

216



ACC






S2

AGGTTGAACTGGCGAAAATGAGCAAGCA

217



GGCT






S3

AAGCCTGCTTGCTCAAAACCGAACGTGGT

218



GTT






S4

AAACACCACGTTCGGAAAACTATGGAGCG

219



GCA






S5

ATGCCGCTCCATAGTAAAGGATTGCGTGC

220



TGA





Sequence set 5-




17 numbering
Optimized sequence





S1

AAGCTGCACCAATCCAAATCGCCAGAAGG

221



TCA








ATGACCTTCTGGCGAAAATGAGCACGACG

222


S2
CAT






S3

AATGCGTCGTGCTCAAAACCGAACAACCT

223



GCT






S4

AAGCAGGTTGTTCGGAAAACTATGGAGCG

224



GCA






S5

ATGCCGCTCCATAGTAAAGGATTGGTGCA

225



GCT





Sequence set 5-




18 numbering
Optimized sequence





S1

AACGCTCGTCAATCCAAATCGCCTCAGGA

226



CAA






S2

ATTGTCCTGAGGCGAAAATGAGCCAACGA

227



CCT






S3

AAGGTCGTTGGCTCAAAACCGAAGCTGGT

228



GTT






S4

AAACACCAGCTTCGGAAAACTATGCCGCA

229



CCT






S5

AAGGTGCGGCATAGTAAAGGATTGACGA

230



GCGT





Sequence set 5-




19 numbering
Optimized sequence





S1

AAGTGCGTCGAATCCAAATCGCCAAGACC

231



TCA






S2

ATGAGGTCTTGGCGAAAATGAGCAGGCTG

232



GAA






S3

ATTCCAGCCTGCTCAAAACCGAAGCAACG

233



TGT






S4

AACACGTTGCTTCGGAAAACTATGCCGCT

234



CCT






S5

AAGGAGCGGCATAGTAAAGGATTCGACG

235



CACT





Sequence set 5-




20 numbering
Optimized sequence





S1

ATCACGCAGCAATCCAAATCGCCATCACA

236



ACG






S2

ACGTTGTGATGGCGAAAATGAGCACGAGC

237



CTT






S3

AAAGGCTCGTGCTCAAAACCGAAGGTTGC

238



ACT






S4

AAGTGCAACCTTCGGAAAACTATGCCGCT

239



CCA






S5

ATGGAGCGGCATAGTAAAGGATTGCTGCG

240



TGA









In the seventh aspect of the present invention, it provides a device for determining the nucleic acid single strand sequence for forming the multimeric complex based on the complementary nucleic acid backbone, which comprises:

    • (M1) an input module, which is used to input annealing algorithm parameters, optimized constraint parameters, and optionally nucleic acid sequences to be optimized;
    • wherein the setting annealing algorithm parameters including: initial annealing temperature, annealing termination temperature, and annealing temperature attenuation coefficient ΔT;
    • the optimized constraint parameters including:
    • {circle around (1)} the number n of the nucleic acid single strand, preferably a positive integer of 3-6;
    • {circle around (2)} the length L of the pairing sequence, preferably the L is of 12-16 bases;
    • {circle around (3)} the dissociation temperature threshold Tm of the pairing region;
    • {circle around (4)} the free energy threshold ΔG°S of the specific pairing region sequence;
    • {circle around (5)} the free energy threshold ΔG°NS of the non-specific pairing;
    • {circle around (6)} the connecting element X2, preferably A, AA, and AAA;
    • {circle around (7)} the dissociation temperature threshold Tm-H of the secondary structure (hairpin);
    • {circle around (8)} the CG proportion PCG in the pairing sequence, preferably the range of PCG is [0.4,0.6),
    • (M2) an optimization operation module, which is configured to perform the following sub steps to obtain optimized nucleic acid single strand sequences or sets thereof:
    • (z1) calculating the objective function value E0 of the initial set S, that is, calculating the sum of the non-specific pairing free energies (ΔG°Ns) between sequences and of the sequence itself, while obtaining the non-specific pairing free energy matrix Cn×n, searching the Si and Sj(1≤i≤n, 1≤j≤n) corresponding to the minimum value in the upper triangular matrix thereof, randomly selecting Si or Sj for an updated operation according to the non-specific pairing free energy of the Si and Sj ΔG°NS(Si, Sj), and then obtaining a new nucleic acid sequence, thereby obtaining a updated sequence set S′;
    • (z2) determining whether the sequences in the set S′ of the previous step meet the set optimized constraint parameter conditions, verifying the following parameters, including the dissociation temperature Tm of the specific pairing region, the free energy ΔG°S of the specific pairing region sequence, the dissociation temperature Tm-H of the secondary structure and the CG proportion PCG. If the above parameters meet the constraint conditions, the step (z3) is proceeded; otherwise, the step (z2) is repeated. If the step is performed 15 times continuously at a certain annealing temperature without obtaining the S′ that meets the conditions, then the set S becomes the set S′ and the next step is proceeded to prevent a dead cycle;
    • (z3) calculating the objective function value E1 of the set S′ of the previous step, and comparing E0 with E1. If E1≥E0, it indicates that the non-specific pairing free energy can be optimized, and the sequence set S′ becomes the sequence set S. If E1<E0, it indicates that the non-specific pairing free energy has not been optimized, and in this case, it is necessary to determine whether to accept the set S′ as S according to the Metropolis criterion; and
    • (z4) the annealing temperature is attenuated according to the set attenuation coefficient ΔT, and the steps (z1), (z2), and (z3) are repeated for the S of the previous step, which is the Monte Carlo-based annealing algorithm, until the annealing temperature reaches the annealing termination temperature. The S={S1, S2, . . . , Sn} of the previous step becomes the nucleic acid single strand sequence for forming the multimeric complex based on the complementary nucleic acid backbone; and
    • (M3) an output module, which is used to output optimized nucleic acid single strand sequences or sets thereof.


In another preferred embodiment, the optimized constraint parameters further comprises: for a tetramer (n=4), using a symmetric sequence to initialize a sequence set S={S1, S2, . . . , Sn} according to the above parameters;


It should be understood that within the scope of the present invention, the above-mentioned technical features of the present invention and the technical features specifically described in the following (such as the embodiments) can be combined with each other to form a new or preferred technical solution, which are not redundantly repeated one by one due to space limitation.





DESCRIPTION OF DRAWINGS


FIG. 1 shows a schematic diagram of the multimer.



FIG. 2 shows a flowchart of the annealing algorithm involved in the present patent.



FIG. 3 shows a schematic diagram of specific sequence pairing of a trimer.



FIG. 4 shows a gel electrophoresis diagram of the nucleic acid backbone assembly of the trimer optimized sequences.



FIG. 5 shows a schematic diagram of specific sequence pairing of the tetramer optimized sequences.



FIG. 6 shows a statistical diagram of the sum of free energies of non-specific pairing regions between the tetramer optimized sequences.



FIG. 7 shows a gel electrophoresis diagram of the nucleic acid backbone assembly of the tetramer optimized sequences.



FIG. 8 shows a schematic diagram of the conversion of a tetramer into a pentamer.



FIG. 9 shows a statistical diagram of the sum of free energies in the unpaired regions of the optimized sequences for the first conversion scheme of pentamers.



FIG. 10 shows a gel electrophoresis diagram of the nucleic acid backbone assembly of the optimized sequences for the first conversion scheme of pentamers.



FIG. 11 shows a schematic diagram of specific sequence pairing of the optimized sequences for the second conversion scheme of pentamers.



FIG. 12 shows a statistical diagram of the sum of free energies in the unpaired regions of the optimized sequences for the second conversion scheme of pentamers.



FIG. 13 shows a gel electrophoresis diagram of the nucleic acid backbone assembly of the optimized sequences for the second conversion scheme of pentamers.



FIG. 14 shows the coupling between G-CSF and L-DNA.



FIG. 15 shows the purification effect of the (L-DNA)-G-CSF conjugate.



FIG. 16 shows the assembly effect of the monovalent, divalent, and trivalent G-CSF complex.



FIG. 17 shows the effect of the L-DNA tetramer framework on the in vitro activity of G-CSF.



FIG. 18 shows the in vitro activity evaluation of the divalent and trivalent G-CSF assembled by L-DNA tetramer.



FIG. 19 shows a. purification of SM(PEG)2-PMO1 using HiTrap Capto MMC; b. gel electrophoresis diagram of coupling efficiency between SM(PEG)2-PMO1 and anti-HSA nanoantibody.



FIG. 20 shows the identification of PMO1(a), SM(PEG)2-PMO1(b), anti-HSA Nb(c), and anti-HSA Nb-PMO1(d) by positive ion mode of liquid chromatography-mass spectrometry.



FIG. 21 shows the separation of nanoantibodies and PMO nanoantibody conjugates using Superdex™ 75 Increase 10/300 GL.



FIG. 22 shows a gel electrophoresis diagram and schematic diagram of the NAPPA-PMO assembly sample. Left: pmo-NAPPA4-HSA(1,2,3); Right: pmo-NAPPA4-HSA(1).



FIG. 23 shows the binding activity of the anti-HSA Nb, anti-HSA Nb-PMO1, and pmo-NAPPA4-HSA(1) with human serum albumin protein detected by ELISA.



FIG. 24 shows the resistance experiment of the pmo-NAPPA4-HSA(1) to nuclease degradation. Left: SDS-PAGE gel electrophoresis diagram of the pmo-NAPPA4-HSA(1) treated with three kinds of nuclease; Right: gel electrophoresis diagram of DDNA-NAPPA4 treated with three kinds of nuclease.





MODES FOR CARRYING OUT THE INVENTION

After extensive and intensive research, the inventors have developed a multivalent protein drug and its library, as well as a preparation method and application thereof for the first time. By using the drug library and preparation method of the present invention, short-acting protein drugs can be quickly, efficiently formed into multivalent complexes at low-cost, and with high yield, to increase the drug half-life, or form high valent antigens with monomer antigens to enhance their immunogenicity according to needs. On this basis, the inventors have completed the present invention.


Specifically, the present invention provides a multivalent protein drug, comprising n protein drug units, wherein each drug unit comprises a drug element moiety of the same kind and different nucleic acid element moieties connected to the drug element moiety; n is a positive integer≥2; n of different nucleic acid element moieties form n-multimers through nucleic acid base complementation, thereby forming the multivalent protein drug; a stable pairing structure with nucleic acid base complementation (rather than complex peptide bonds or other chemical modifications, etc.) of the multivalent protein drug of the present invention can be formed through rapid assembly (such as 1 minute). Experiments have shown that the molecular weight of the drug of the present invention can be increased through high-valent formation, thereby extending its half-life in animals.


In addition, based on the same implementation method, the drug element can also be an antigen used for vaccine development; the difference is that each antigen unit comprises the same or different antigen element moiety, as well as different nucleic acid element moieties connected to the antigen element moiety; n of different nucleic acid element moieties can form a n-multimer through nucleic acid base complementation, thereby forming the multivalent antigen.


Finally, the present invention provides a highly optimized nucleic acid sequence library, including nucleic acid sequence groups that can be efficiently and accurately assembled into 2-5 aggregates, for rapid and accurate self-assembly of the aforementioned drugs or antigen units into multivalent macromolecular complexes.


Term

Unless otherwise defined, all technical and scientific terms used herein have the same meanings as those commonly understood by common technicians in the art to which the present invention belongs. As used herein, when referring to specific listed values, the term “about” means that the values can vary by no more than 1% from the listed values. For example, as used herein, the expression of “about 100” includes all values between 99 and 101 (such as 99.1, 99.2, 99.3, 99.4, etc.).


Drug D

In the present invention, the drug element moiety comprises protein drugs and polypeptide drugs.


Typically, the protein drugs include but are not limited to cytokines, hormones (such as insulin, growth hormone, etc.), antibody drugs, and polypeptides.


In a preferred embodiment of the present invention, the protein drug is G-CSF for treating leukopenia.


Antigen Library A

The present invention provides an antigen library, which comprises N of antigen units;

    • wherein, the antigen units comprise an antigen element moiety and a nucleic acid element moiety connected to the antigen element moiety; different nucleic acid element moieties can form a n-multimer through nucleic acid base complementation, thereby forming the multivalent antigen;
    • wherein, the antigen element moiety is a protein antigen or a polypeptide antigen;
    • typically, the protein antigen or polypeptide antigen includes but is not limited to a virus, bacterial protein, structural region thereof, and fragment;


In the present invention, the antigen element moiety is partially selected from M of different antigen proteins in the library, M≤N; M of different antigen proteins contain different proteins in a certain virus or mutants of the same protein from different virus strains;


In a preferred embodiment of the present invention, the protein antigen is derived from novel coronavirus SARS-COV-2; specifically, it is a high-valent antigen formed by the receptor binding domain (RBD) of the viral spike protein.


Left-Handed Nucleic Acid

Left-handed nucleic acid refers to the mirror image of the natural right-handed nucleic acid (D-nucleic acid), which can be divided into left-handed DNA (L-DNA) and left-handed RNA (L-RNA). Left-handed (chiral center) mainly exists in the deoxyribose or ribose portion of nucleic acids, exhibiting mirror flipping. Therefore, Left-handed nucleic acid cannot be degraded by ubiquitous nuclease (such as exonuclease and endonuclease) in plasma.


Preparation Method
1. Design and Preparation of L-Nucleic Acid Strand Framework

According to the present invention, an L-nucleic acid strand framework is formed by base complementation of two or more L-nucleic acid single strands. The 5′ end or 3′ end of each L-nucleic acid single strand is activated into groups (such as NH2, etc.) that can be subsequently modified, and then one end of the linker (such as SMCC, SBAP, etc.) is coupled with the activate group on the L-nucleic acid single strand. L-nucleic acids with the linkers can be assembled into the desired L-nucleic acid strand framework. In another preferred embodiment, the activated functional groups at the 5′ or 3′ end of the L-nucleic acid single strand (such as aldehide, maleimide, etc.) are already included in nucleic acid synthesis. After confirming that the L-nucleic acid with the linker can successfully self-assemble into a framework, the L-nucleic acid single strand with the linker can be separately coupled with antibodies for subsequent assembly. The L-nucleic acid framework of the present invention can be prepared by the following steps.


1.1 Design of L-Nucleic Acid Single Strand for Rapid Self-Assembly

The required multivalent number n (such as a trimer, tetramer) is determined; the required number n of L-nucleic acid single strands based on the multivalent number n is determined; a corresponding number of L-nucleic acid single strand sequences is designed, and the stability of the target nucleic acid framework is regulated by optimizing base pairing, and the possibility of non-specific pairing between nucleic acid strands is reduced. The details of nucleic acid sequence design are specifically described in the summary of the invention and embodiments.


1.2 Activation of L-DNA or L-RNA

The activation of L-nucleic acid involves modification of its active group at the 5′ end (X1) or 3′ end (X3) and subsequent conjugation with the linkers. The modification of active groups can be customized by nucleic acid synthesis companies; The linkers generally have a bifunctional group, that is, one end can be coupled with an active nucleic acid group, and the other end can be connected to specific sites on the protein (such as NH3, SH).


According to a preferred embodiment of the present invention, all L-nucleic acids that make up the framework are modified with aldehydes at the 5′ end, thereby completing the activation of L-nucleic acids and subsequently coupling to the N-terminal a-amine of the protein.


2. Preparation Method of Protein-L-Nucleic Acid Complex

First, the 5′ or 3′ end of L-nucleic acids is modified with aldehydes, and then, the aldehyde groups of L-nucleic acids are specifically connected to the N-terminal NH3 of the protein through a reductive amination reaction under low pH (5-6) conditions.


Algorithm and Algorithm Optimized Nucleic Acid Sequences

The present invention also provides a method and device of determining a nucleic acid single strand sequence for forming a multimeric complex based on a complementary nucleic acid backbone. Preferably, the method comprises a preferred algorithm of the present invention.


Typically, the nucleic acid sequence library optimized by computer algorithms of the present invention can include: (a) a nucleic acid sequence that can self-assemble into a trimer through base complementation; (b) a nucleic acid sequence that can self-assemble into a tetramer through base complementation; and (c) a nucleic acid sequence that can self-assemble into a pentamer through base complementation.


The complexes formed by representative nucleic acid sequences are trimer molecules, tetramer molecules, and pentamer molecules as shown in FIG. 1.


Preferably, the nucleic acid sequence W of the present invention has the structure shown in formula 1:









W
=


X

1

-

R

2

-

R

2

-

X

3






(
1
)









    • wherein,

    • R1 is the complementary base pairing region 1;

    • R2 is the complementary base pairing region 2;

    • X1, X2, and X3 are none or redundant nucleic acids, independently;

    • the length of RI and R2 is 14-16 bases;

    • the length of X1 and X3 is 0-5 bases;

    • X2 has a length of 0-3 bases and a sequence of A, AA, AGA or AAA;

    • wherein, the R1 of the nucleic acid sequence forms a target pairing with the R2 of different nucleic acid sequences, while the R2 of the nucleic acid sequence forms a target pairing with the R1 of another nucleic acid sequence.





In the present invention, the self-pairing of any region of the nucleic acid sequence belongs to non-target pairing, which needs to be avoided in design.


Preferably, the nucleic acid sequence of the present invention can be designed or optimized using a computer algorithm of Simulated Annealing (SA).


The computer algorithm minimizes the free energy (ΔG°) of the DNA double strand structure formed by target pairing, while maximizing the ΔG° of non-target pairing;

    • specifically, the stability of DNA double strand structure depends on the base pairs of each nearest neighbor in the sequence; there may be 10 different nearest neighbor interactions in any Watson-Crick DNA double strand structure, and these pairing interactions are: AA/TT; AT/TA; TA/AT; CA/GT; GT/CA; CT/GA; GA/CT; CG/GC; GC/CG; and GG/CC;
    • more specifically, the ΔG° values of the 10 base pairs mentioned above can be calculated by enthalpy (ΔH°) and entropy (ΔS°) at any temperature; enthalpy, entropy, and free energy data for 10 sequences summarized in Table 1.









TABLE 1







Thermodynamic data for 10 sequences














ΔH°
ΔS°
ΔG1°
ΔG2°




kcal/
cal K−1
kcal/
kcal/



Sequence
mol
mol−1
mol
mol

















AA/TT
−9.1
−24.0
−1.9
−1.94



AT/TA
−8.6
−23.9
−1.5
−1.47



TA/AT
−6.0
−16.9
−0.9
−0.96



CA/GT
−5.8
−12.9
−1.9
−1.95



GT/CA
−6.5
−17.3
−1.3
−1.34



CT/GA
−7.8
−20.8
−1.6
−1.6



GA/CT
−5.6
−13.5
−1.6
−1.57



CG/GC
−11.9
−27.8
−3.6
−3.61



GC/CG
−11.1
−26.7
−3.1
−3.14



GG/CC
−11.0
−26.6
−3.1
−3.07







Note:



ΔH°, ΔS° and ΔG1° are measured under conditions of 1 M NaCl, 25° C., and pH 7. ΔG2° is measured from the IDT website (https://sg.idtdna.com/calc/analyzer).






Using the thermodynamic values in Table 1, the enthalpy ΔH° and free energy ΔG° values of DNA oligomers can be effectively predicted by the nearest neighbor method. Taking the complementary pairing of GGAATTCC/CCTTAAGG as an example, using the nearest neighbor method, it is calculated to be ΔG°=−14.63 kcal/mol. The annealing algorithm not only optimizes the ΔG°NS values, but also implements constraints on the dissociation temperature of nucleic acid sequence pairing (Tm) to ensure Tm>50° C. (when L=14 bases) of R1 and R2 regions. The nearest neighbor model is based on thermodynamic calculations and accurately predicts the stability of DNA double strands. The prediction of a given base sequence is provided by the model based on the nearest neighbor base pairs. The calculation of the unwinding temperature involves enthalpy (ΔH°) and entropy (ΔS°), and the calculation method is as follows:










T
m

=



Δ


H





Δ



S


[

Na
+

]


+

R
×

ln
(

CT
)




-

2

7


3
.
1


5






(
1
)









    • wherein, the R is a constant (1.987 cal K−1 mol−1), the CT is the strand concentration given as 0.1 μM, the ΔS°[Na+] is the entropy value of the DNA double strand at a given sodium ion concentration, and the ΔH° is the enthalpy value under given conditions.





Simulated annealing is a universal probability algorithm, and it is a method for approximate optimal solution of the problem, which is designed according to Monte Carlo's ideas, aiming to find the approximate optimal solution in a large searching space within a certain period of time. The idea of simulated annealing algorithm originates from the annealing process of solid materials in physics: first, the solid is fully heated, and then slowly cooled. When heated, the internal energy of free motion of particles inside the solid increases. Later, as the temperature gradually decreases, the particles tend to become orderly and reach equilibrium at each temperature. If the temperature drops slowly enough near the condensation point, the ground state can be reached, and the internal energy is minimized. According to the Metropolis criterion, the probability of a particle reaching equilibrium at temperature T is exp(−ΔE/(KT)), wherein the E is the internal energy at temperature T, the ΔE is its variation, and the K is the Boltzmann constant. For the combinatorial optimization problem, there is a similar process. The solid micro state i is simulated as a solution X, the objective function is equivalent to the internal energy Ei of state i, and the solid temperature is simulated with the control parameter T. The iteration process of “generating a new solution→calculating the objective function difference→judging whether to accept→accepting/discarding” for each value of T is repeated, and the T value is gradually attenuated. During the iteration process, a poor solution is limitedly accepted according to the Metropolis criterion, and the probability of accepting the poor solution is gradually approaching zero, so as to find the global optimal solution at all possible when the algorithm terminates.


In the present invention, the free energy of the non-target pairing region between the single strands of a multimeric nucleic acid is simulated as internal energy. While ensuring the free energy and dissociation temperature of the target pairing region between the nucleic acid single strands, the free energy of the non-target pairing region is optimized through iteration of a simulated annealing algorithm. Finally, the optimized single strand is more conducive to the assembly of the multimer. The optimized objective function used herein is as follows:










E
=






i
=
1




n








j
=
i




n




-
Δ




G


(


S
i

,

S
j


)





,

n

1





(
2
)









    • ΔG°(Si, Sj) is the free energy of the non-target pairing between Si sequence and Sj sequence, and Σi=1nΣj=inΔG°(Si, Sj) is the sum of the free energies of non-target pairing between all sequences, with a negative value;





wherein, the larger the negative value, the more beneficial it is to reduce the non-target pairing. Therefore, an objective function is constructed based on the idea of minimizing energy using a degradation algorithm, and a minus sign is added before ΔG°(Si, Sj), converting it into a positive number. At this point, the smaller the value of Σi=1nΣj=i−ΔG°(Si, Sj), the more beneficial it is to reduce non-target pairing.


The flowchart of a representative algorithm of the present invention is shown in FIG. 2.


In the present invention, the simulated annealing algorithm introduces random factors, and in each iteration update process, it will accept a solution that is worse than the current one with a certain probability, so it is possible to jump out of the local optimal solution and reach the global optimal solution.


Multivalent Macromolecular Complexes

The present invention also provides a multivalent macromolecular complex with improved drug half-life and activity, which is formed by using the nucleic acid multimer to mediate protein drugs, wherein the nucleic acid multimer is designed using the above algorithm.


Preferably, the nucleic acid sequence is a nucleic acid sequence set which can be specifically assembled into n-multimers in a nucleic acid sequence library;

    • preferably, in the present invention, the protein drug is a protein drug that needs multivalent formation to increase its half-life or activity.


Typically, each nucleic acid strand of the nucleic acid sequence set is connected to the protein drug, forming a protein drug-nucleic acid strand unit, with the structure shown in formula 2:





D-[L-Wi], i=1 to n   (2)

    • wherein,
    • D is the protein drug element moiety;
    • each Wi is independently a nucleic acid sequence; the nucleic acid sequence is selected from the group consisting of: left-handed nucleic acid, peptide nucleic acid, locked nucleic acid, thio-modified nucleic acid, 2′-fluoro modified nucleic acid, 5-hydroxymethylcytosine nucleic acid, and combinations thereof; the nucleic acid sequence has the structure shown in formula 1, and is selected from the nucleic acid sequence set that can form n-multimers mentioned above;
    • L is a linker; the linker moiety has already been included in the synthesis or preparation of Wi, connected to the X1 or X3 of Wi (see formula 1);
    • “-” is a covalent bond;


In another preferred embodiment, the drug element moiety is selected from the group consisting of: protein drugs and polypeptide drugs that need to increase their molecular weights, thereby increasing their half-lives, and protein drugs and peptide drugs that need multivalent formation to increase their activities;


In another preferred embodiment, the L has aldehyde, NHS ester, or similar functional groups near the D-end, for connecting the N-terminal α-amine or lysine ε-amine on D;


In another preferred embodiment, the L has a maleimide functional group or haloacetyl (such as bromoacetyl, iodoacetyl, etc.) functional group near the D end, for connecting the free thiol (—SH) functional group on D;


In another preferred embodiment, the D is selected from the group consisting of: natural proteins, recombinant proteins, chemically modified proteins, and synthesized polypeptides;


In another preferred embodiment, the D can have a site-directed modification or site-directed addition of non-natural amino acids for connecting the L-Wi of formula 1;

    • the protein drugs are respectively connected to different L-Wi that can be assembled into n-multimers, forming protein drug self-assembly units, D-[L-W1], D-[L-W2], . . . , D-[L-Wn];
    • the protein drug self-assembly units, D-[L-W1], D-[L-W2], . . . , D-[L-Wn], are mixed in equimolar solution to assemble multivalent molecular complexes of protein drugs.


The present invention also provides a multivalent macromolecular complex which enhances the effectiveness of vaccines in inducing neutralizing antibody production in vivo, and it is formed by using the nucleic acid multimer to mediate one or more antigens, wherein the nucleic acid multimer is designed using the above algorithm;

    • wherein, the nucleic acid sequence is a nucleic acid sequence set which can be specifically assembled into n-multimers in a nucleic acid sequence library;
    • wherein, the antigen is an antigen or antigen library; the antigen library comprises M of different antigen proteins, 1≤M≤n;
    • each nucleic acid strand of the nucleic acid sequence set is connected to an antigen in the antigen library, forming an antigen-nucleic acid strand unit with the structure shown in formula 3:











A
k

-

[

L
-

W
i


]


,

i
=

1

n


,

k
=

1

M






(
3
)









    • wherein,

    • Ak is the antigen k in the antigen library; one Ak corresponds to one or more L-Wi (such as A1-[L-W1], A1-[L-W2], A2-[L-W3], A3-[L-W4]);

    • the other aspects of formula 3 are the same as those of the above formula 3;

    • the antigen proteins are respectively connected to different L-Wi that can be assembled into n-multimers, forming antigen self-assembly units, such as A1-[L-W1], A2-[L-W2], . . . , A3-[L-WN];

    • the antigen self-assembly units are mixed in equimolar solution to assemble multivalent antigen complexes.





Main Advantages of the Present Invention Are





    • (1) the present invention can achieve multivalent formation of ready-made short-acting protein drugs without the need of the reconstruction of fusion proteins or complex chemical modification and cross-linking conditions, thereby improving their half-lives and activities; aldehyde modification of L-nucleic acids can specifically connect the N-terminal amine of proteins, forming protein drug units that can self-assemble into oligomers;

    • (2) the protein drug units (protein-nucleic acid connecting products) of the present invention can achieve multivalent formation of protein drugs within one minute through the mediation of a left-handed nucleic acid strand;

    • (3) in terms of vaccine development, the present invention can achieve multivalent formation of monomeric protein antigens, improving their immunogenicity;

    • (4) in terms of vaccine development, the present invention can also achieve the assembly of antigen mutations and subtypes of different virus or bacterial strains into diverse high-valent antigens, inducing a wider range of neutralizing antibodies.





The present invention will be further illustrated below with reference to the specific examples. It should be understood that these examples are only to illustrate the invention, not to limit the scope of the invention. The experimental methods with no specific conditions described in the following examples are generally performed under the conventional conditions (e.g., the conditions described by Sambrook et al., Molecular Cloning: A Laboratory Manual (New York: Cold Spring Harbor Laboratory Press, 1989), or according to the manufacturer's instructions. Unless indicated otherwise, percentages and portions are weight percentages and weight portions.


EXAMPLE 1: ASSEMBLY DESIGN, SYNTHESIS, AND VALIDATION OF A TRIMER NUCLEIC ACID BACKBONE

Three nucleic acids that can be paired according to the shape shown in FIG. 1(A) are designed. Specifically, the specific complementary pairing of nucleic acid single strands R1 with nucleic acid single strands R6, instead of other nucleic acid single strands, is performed. Similarly, the specific complementary pairing of R3, R5 with R2, R4 respectively, instead of other nucleic acid single strands, is performed. And the free energy of specific complementary pairing (ΔG°S) is much smaller than that of non-specific pairing (ΔG°NS). The free energy of specific complementary pairing (ΔG°S) is less than −29 kcal/mol, while the free energy of non-specific pairing (ΔG°NS) is greater than −7 kcal/mol. In this way, the form of trimers is the most stable in the reaction system. The specific implementation steps of trimer optimization are described in FIG. 2; the annealing parameters are: annealing initial temperature To=50° C., annealing termination temperature Tf=0.12° C., annealing temperature decay coefficient ΔT=0.98 (each attenuation of annealing initial temperature is 0.98 of the current value); the optimized constraint parameters are: pairing sequence length L=16 bases, dissociation temperature threshold Tm>54° C., pairing sequence free energy threshold: ΔG°S<−29 kcal/mol, and non-specific pairing free energy threshold: ΔG°NS>−7 kcal/mol. The specific implementation steps for optimization are as follows:


Sequence initialization: pairing sequences R={R1, R3, R5} are initialized based on parameters, and R6, R2 and R4 are obtained according to base complementation. After splicing the six sequences as shown in FIG. 1(A), the initial sequence set S={S1, S2, S3} is finally obtained.


Generation of new solution: new solutions S′ are obtained by updating the sequence set S. First, ∥−αG°(Si, Sj)∥is calculated to identify the two sequences Si and Sj of non-specific pairings that have the greatest impact on the target pairing in the set S. Then, Si or Sj in non-target pairing regions is randomly selected for update according to the ΔG°NS(Si, Sj), thereby obtaining a new nucleic acid sequence. The dissociation temperature of the pairing regions (Tm) of this nucleic acid sequence is tested to see if it is greater than 54° C., and the free energy of the pairing regions (ΔG°S) is tested at the same time to see if it is less than −29 kcal/mol. If the constraint requirements of dissociation temperature and pairing region free energy (ΔG°S) are not met, the update will be repeated. If the constraint requirements of dissociation temperature and pairing region free energy (ΔG°S) are met, the S is updated according to the principle of base complementation, and finally obtaining a new sequence set S′. If the new nucleic acid sequence obtained after fifteen updates still does not meet the constraints of dissociation temperature and pairing region free energy (ΔG°S), in order to prevent a dead cycle, the set S becomes the new solution S′.


Optimization judgment: the objective function values Ei and Ei+1 of the set S and S′ are calculated respectively according to formula 2. If Ei+1−Et≥0, it indicates that the update has optimized the non-target pairing free energy (ΔG°NS), then S←S′, the new solution S′ becomes S. If Ei+1−Ei<0, it indicates that the update has obtained a deteriorating solution. According to the Metropolis criterion, the probability p=exp(−ΔE/To) is calculated and r is randomly generated, rε[0,1). If p>r, then accept the deteriorating solution, otherwise reject the deteriorating solution. Finally, the annealing initial temperature decays to 0.98 of the current value, generating the next new solution until it decays to the annealing termination temperature, thereby obtaining an optimized sequence set S.


The optimized sequences in Table 2_1 are obtained through the above algorithm optimization. The purpose of the 5′ end A of each sequence is to modify the active group for the subsequent coupling with linkers. In the non-target pairing free energy (ΔG°NS) matrix table and the target pairing region parameter index table, some main parameter values are counted. The schematic diagram of specific sequence pairing of the trimer optimized sequences is shown in FIG. 3. From the corresponding gel electrophoresis diagram of the nucleic acid backbone (FIG. 4), it can be seen that Lane9 is an artificially designed trimer with an unclear and trailing band. Lane10 is the main band formed by the optimized sequences S1, S2 and S3, indicating the formation of a trimer and exhibiting high stability.









TABLE 2_1





Trimer initialized sequences and optimized sequences

















Sequence

SEQ ID


numbering
Initialized sequences
NO:





S1

AGTGATCCGAAGTCGACAAACGTATTAGCGC

241



TCGAT






S2

AATCGAGCGCTAATACGAAAGTGCAATGCGT

242



CGATG






S3

ACATCGACGCATTGCACAAAGTCGACTTCGG

243



ATCAC





Sequence




numbering
Optimized sequences





S1

ACCACCGTGTATGACCTAAAAGTGACAGCAC

244



ATCGC






S2

AGCGATGTGCTGTCACTAAAACAGGCTCTAC

245



GAGGA






S3

ATCCTCGTAGAGCCTGTAAAAGGTCATACAC

246



GGTGG
















TABLE 2_2







Trimer initialized sequence and optimized sequence


free energy (ΔGNS°) matrix








Trimer initialized
Trimer optimized sequence


free energy matrix
free energy matrix














Sequence
S1
S2
S3
Sequence
S1
S2
S3

















S1
−13.09
−8.1
−8.24
S1
−3.61
−5.19
−4.89


S2

−13.09
−8.1
S2

−4.89
−5.19


S3


−16.53
S3


−4.95
















TABLE 2_3







Parameter indicators for target pairing


regions of trimer optimized sequences
















Tm
Tm



Pairing


° C.
° C.



sequences
CG %
ΔGS°
(TF)
(IDT)







R1
56.3%
−29.66
55.3
50.9



R3
56.3%
−29.55
57.1
52.5



R5
56.3%
−29.61
54.4
51.6










EXAMPLE 2: ASSEMBLY DESIGN, SYNTHESIS, AND VALIDATION OF A TETRAMER NUCLEIC ACID BACKBONE

Four nucleic acids that can be paired according to the shape shown in FIG. 1(B) are designed, wherein, the specific complementary pairing of nucleic acid single strands R1, R3, R5 and R7 with nucleic acid single strands R8, R2, R4, and R6 respectively, instead of other nucleic acid single strands, is performed. And the free energy of specific complementary pairing (ΔG°S) is much smaller than that of non-specific pairing (ΔG°NS). The free energy of specific complementary pairing (ΔG°S) is less than −27.4 kcal/mol, while the free energy of non-specific pairing (ΔG°NS) is greater than −7 kcal/mol. In this way, the form of tetramers is the most stable in the reaction system. The specific implementation steps of tetramer optimization are described in FIG. 2; the annealing parameters are: annealing initial temperature To=50° C., annealing termination temperature Tf=0.12° C., annealing temperature decay coefficient ΔT=0.98 (each attenuation of annealing initial temperature is 0.98 of the current value); the optimized constraint parameters are: pairing sequence length L=14 bases, dissociation temperature threshold Tm>52° C., pairing sequence free energy threshold: ΔG°S<−27.4 kcal/mol, and non-specific pairing free energy threshold: ΔG°NS>−7 kcal/mol.


The specific implementation steps for optimization are as follows:


Sequence initialization: pairing sequences R={R1, R3, R5, R7} are initialized based on parameters, and R8, R2, R4, and R6 are obtained according to base complementation. After splicing the eight sequences as shown in FIG. 1(B), the initial sequence set S={S1, S2, S3, S4} is finally obtained. During the tetramer experiment, it is found that a fixed core structure can effectively improve the assembly efficiency of nucleic acid sequences. Based on the nucleic acid sequence of formula 1, the nucleic acid sequence W2 is developed, which has the structure of formula 4:










W

2

=


X

1

-

Q

1

-

C

1

-

X

2

-

C

2

-

Q

2

-

X

3






(
4
)









    • wherein, C1 and C2 are fixed core structure parts, and Q1 and Q2 are sequences other than fixed core structure.












TABLE 3





Nucleic acid sequences S1, S2, S3 and


S4 containing fixed core structures


















S1
X1-Q1-AATCC-X2-TGAGC-Q2-X3



S2
X1-Q3-GCTCA-X2-CCGAA-Q4-X3



S3
X1-Q5-TTCGG-X2-ACTAT-Q6-X3



S4
X1-Q7-ATAGT-X2-GGATT-Q8-X3










Generation of new solution: the same as the generation of new solution in Example 1. The difference is that if a fixed core structure is used, the update will not include the fixed core structure part; constraint parameters need to be strictly followed, the dissociation temperature (Tm) of the pairing regions of the new nucleic acid sequence should be greater than 52° C., and the free energy of the pairing regions (ΔG°S) should be less than −27.4 kcal/mol.


Optimization judgment: the same as the optimization judgment in Example 1.


The optimized sequences in Table 4_1 are obtained through the above algorithm optimization, and the specific pairing diagram thereof is shown in FIG. 5. The statistical line graph of the non-target pairing free energy during the optimization process of the tetramer optimized sequences is shown in FIG. 6. In the case of adding connecting elements, the present invention seeks to avoid the significant impact on the free energy (ΔG°NS) of the optimized sequences which are not added with the connecting elements. Regardless of whether the connecting elements are added, the sum of the free energy (ΔG°NS) matrices of the optimized sequences should be as large as possible. Therefore, during the optimization process, the objective function values of the optimized sequences without the addition of the connecting elements are also counted, and the final detection is only performed on the optimized sequences with the addition of the connecting elements. From the corresponding gel electrophoresis diagram of the nucleic acid backbone (FIG. 7), it can be seen that Lane15 is an artificially designed tetramer with a trailing band. Lane16 is the main band formed by the optimized sequences S1, S2, S3 and S4, which is around 100 bp, indicating the formation of a tetramer and exhibiting high stability.


The above implementing steps for tetramers mainly optimize the non-specific pairing free energy between sequences, and the secondary structure formed by the self-folding of sequences also has a significant impact on the assembly of tetramers. If the dissociation temperature of the secondary structure formed by the sequence itself is too high, once this stable secondary structure is formed, it is difficult to break this state, making it difficult for the tetramer to assemble. Therefore, it is necessary to control the dissociation temperature of the secondary structures corresponding to the four sequences of the assembled tetramer to be not too high. For tetramers, it is necessary to control the dissociation temperature of the secondary structures of the four nucleic acid sequences. If a symmetric structure is used (as shown in FIG. 1, R1, R3, R5, R7 maintain symmetry with R4, R6, R8, R2, respectively), similar secondary structures will appear between the two sequences, and the dissociation temperature of the secondary structure between the two sequences is not significantly different. At this point, only the secondary structure of the two nucleic acid sequences needs to be controlled to achieve the previous effect. Therefore, symmetry is beneficial for controlling the secondary structure in tetramer optimization.









TABLE 4_1





Tetramer initialized sequences and optimized


sequences



















SEQ


Sequence

ID


numbering
Initialized sequences
NO:





S1

AACCTGGTACAATCCAAATGAGCTACACTA

247



GC






S2

AGCTAGTGTAGCTCAAAACCGAAGTATCGA

248



TT






S3

AAATCGATACTTCGGAAAACTATAGTGAGT

249



TG






S4

ACAACTCACTATAGTAAAGGATTGTACCAG

250



GT





Sequence




numbering
Optimized sequences





S1

AGGCGATCACAATCCAAATGAGCGTGTTAC

251



GG






S2

ACCGTAACACGCTCAAAACCGAAGTGCCA

252



ATT






S3

AAATTGGCACTTCGGAAAACTATGCGGCTG

253



CT






S4

AAGCAGCCGCATAGTAAAGGATTGTGATC

254



GCC
















TABLE 4_2







Tetramer initialized sequence and optimized sequence free energy (ΔGNS°) matrix








Tetramer initialized
Tetramer optimized sequence


free energy matrix
free energy matrix
















Sequence
S1
S2
S3
S4
Sequence
S1
S2
S3
S4



















S1
−6.34
−6.34
−5.85
−5.13
S1
−4.95
−6.75
−6.97
−6.75


S2

−9.69
−5.19
−5.85
S2

−5.36
−6.75
−5.37


S3


−9.69
−4.99
S3


−5.36
−6.75


S4



−9.27
S4



−4.62
















TABLE 4_3







Parameter indicators for pairing


regions of tetramer optimized sequences














Tm
Tm


Pairing
CG

° C.
° C.


sequences
%
ΔGS°
(TF)
(IDT)














R1
57.1%
−27.76
53.8
46.2


R3
57.1%
−27.44
53.0
48.0


R5
50.0%
−28.62
53.7
45.4


R7
57.1%
−28.58
53.6
50.2









EXAMPLE 3: ASSEMBLY DESIGN, SYNTHESIS, AND VALIDATION OF A PENTAMER NUCLEIC ACID BACKBONE

The tetramer optimized sequences exhibit excellent assembly performance, largely due to the lack of pairing in the central region of the tetramer, which means that the fixed core structure provides sufficient freedom for the tetramer and does not form complicated complexes in the central region. Therefore, in order to integrate the tetramer sequences and the fixed core structure in the optimization of the pentamer sequences, two schemes of utilizing the tetramer sequences and the fixed core structure in Example 2 are explored and designed.


The first conversion scheme: five nucleic acids that can be paired according to the shape shown in FIG. 8(B) are designed. This scheme preserves the partial sequence and complete fixed core structure of the tetramer in Example 2, and the core structure is not opened. Open up the tetramer at R1 and R8 of the original tetramer, except for the core structure, and add two nucleic acid sequences R9 and R10 with a length of 14. Among them, nucleic acid single strands R1, R3, R5, R7 and R9 are specifically complementary paired with nucleic acid single strands R10, R2, R4, R6 and R8, respectively, without pairing with other nucleic acid single strands. The free energy of specific complementary pairing (ΔG°S) is less than −27.4 kcal/mol, while the free energy of non-specific pairing (ΔG°NS) is greater than −7 kcal/mol. In this way, the form of pentamers is the most stable in the reaction system. The specific implementation steps for optimizing the first conversion scheme of pentamers are described in FIG. 2; the annealing parameters are: annealing initial temperature T0=50° C., annealing termination temperature Tf=0.12° C., and annealing temperature decay coefficient ΔT=0.9 (each attenuation of annealing initial temperature is 0.9 of the current value. Due to the use of partial sequences of tetramers, there are fewer updated regions and faster annealing temperature decay); the optimized constraint parameters are: pairing sequence length L=14 bases, dissociation temperature threshold Tm>52° C., pairing sequence free energy threshold: ΔG°S<−27.4 kcal/mol, and non-specific pairing free energy threshold: ΔG°NS>−7 kcal/mol.


In this scheme, a complete tetramer fixed core structure is used, and the nucleic acid sequence W5 is developed based on the nucleic acid sequence of formula 4, while the nucleic acid sequence W6 is developed based on the nucleic acid sequence of formula 4. W5 has the structure of formula 7, and W6 has the structure of formula 8:










W

3

=


X

1

-

R

1

-

X

2

-

C

1

-

X

2

-

C

2

-

Q

1

-

X

3






(
5
)













W

4

=


X

1

-

Q

1

-

C

1

-

X

2

-

C

2

-

X

2

-

R

1

-

X

3







(
6
)








This scheme includes sequences of four structures: formula 1, formula 4, formula 5, and formula 6.









TABLE 5





Nucleic acid sequences S1, S2, S3 , S4


and S5 containing partial core structures


















S1
X1-R9-X2-R10-X3



S2
X1-R1-X2-AATCC-X2-TGAGC-Q1-X3



S3
X1-Q2-GCTCA-X2-CCGAA-Q3-X3



S4
X1-Q4-TTCGG-X2-ACTAT-Q5-X3



S5
X1-Q6-ATAGT-X2-GGATT-X2-R8-X3










The specific implementation steps for optimization are as follows:

    • pairing sequences R={R9, R10} are initialized based on parameters, and R8 and R1 are obtained according to base complementation, R2, R3, R4, R5, R6 and R7 are from the homonymous optimized tetramer sequences in Example 2. S1 is obtained by concatenating R9 with R10, S2 is obtained by concatenating R1 with R2 according to formula 5, S3 is obtained by concatenating R3 with R4, S4 is obtained by concatenating R5 with R6, S5 is obtained by concatenating R7 with R8 according to formula 6, and the initialized sequence set S={S1, S2, S3, S4, S5} is finally obtained.


Generation of new solution: Randomly select a sequence from R1, R8, R9 and R10 to update, and obtain a new nucleic acid sequence. Check whether the dissociation temperature of this nucleic acid sequence is greater than 52° C., and whether the free energy of the pairing region (ΔG°S) is less than −27.4 kcal/mol. If the constraint requirements of dissociation temperature and free energy of the pairing regions (ΔG°S) are not met, repeat the update. If the constraint requirements of dissociation temperature and free energy of the pairing region (ΔG°S) are met, update S according to the principle of base complementation to obtain a new sequence set S′. If the new nucleic acid sequence obtained after fifteen updates does not meet the constraint requirements of dissociation temperature and pairing region free energy (ΔG°S), in order to prevent dead circulation, S becomes a new solution S′. Because this scheme uses a complete fixed core structure and part of tetramer sequences, during the process of generating a new solution, it is necessary to ensure that the fixed core structure and the retained tetramer sequences remain unchanged.


Optimization judgment: The same as the optimization judgment in Example 1. The difference is that the annealing initial temperature decays to 0.9 of the current value.


The optimized sequences of the first pentamer conversion scheme in Table 6_1 are obtained through this optimized algorithm, and FIG. 9 shows the statistical line graph of the sum of the free energy values of the non-target pairing regions between the sequences during this optimization process (the free energy of the non-target pairing regions between S3 and S3, S3 and S4, S4 and S4 are not included in the statistics). Lane9 in FIG. 10 shows the main band formed by the optimized sequence assembly, indicating the formation of pentamers and exhibiting high stability.









TABLE 6_1





Initialized sequences and optimized sequences for the first


conversion scheme of pentamers

















Sequence

SEQ ID


numbering
Initialized sequences
NO:





S1

ATCACAGAGCGCGTAAAAACGCCACTCATGG

255



A






S2

ATCCATGAGTGGCGTAAAAATCCAAATGAGC

256



GTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCAAT

257



T






S4

AAATTGGCACTTCGGAAAACTATGCGGCTGC

258



T






S5

AAGCAGCCGCATAGTAAAGGATTAAATACGC

259



GCTCTGTGA





Sequence




numbering
Optimized sequences





S1

ATTCAGGCGACTCCTAAAAGCACGACGATGG

260



T






S2

AACCATCGTCGTGCTAAAAATCCAAATGAGC

261



GTGTTACGG






S3

ACCGTAACACGCTCAAAACCGAAGTGCCAAT

262



T






S4

AAATTGGCACTTCGGAAAACTATGCGGCTGC

263



T






S5

AAGCAGCCGCATAGTAAAGGATTAAAAGGAG

264



TCGCCTGAA
















TABLE 6_2







Initialized sequence and optimized sequence free energy


(ΔGNS°) matrix for the first conversion scheme of pentamers








Petramer initialized
Pentamer optimized sequence


free energy matrix
free energy matrix


















Sequence
S1
S2
S3
S4
S5
Sequence
S1
S2
S3
S4
S5





















S1
−10.36
−8.09
−9.92
−8.16
−6.75
S1
−4.67
−6.91
−6.75
−6.69
−6.75


S2

−7.85
−8.16
−6.97
−9.92
S2

−4.95
−6.91
−6.97
−6.75


S3


−5.36
−6.75
−6.75
S3


−5.36
−6.75
−5.19


S4



−5.36
−6.75
S4



−5.36
−6.75


S5




−10.36
S5




−4.85
















TABLE 6_3







Parameter indicators for pairing


regions of pentamer optimized sequences 1














Tm
Tm


Pairing
CG

° C.
° C.


sequences
%
ΔGS°
(TF)
(IDT)





R9
57.1%
−27.63
53.2
48.8


R1
57.1%
−27.56
54.1
50.0


R3
53.6%
−27.42
53.0
48.0


R5
53.6%
−28.59
53.7
45.4


R7
57.1%
−28.57
53.6
50.2









The second conversion scheme: five nucleic acids that can be paired according to the shape shown in FIG. 8(C) are designed. This scheme only retains the fixed core structure sequence of the tetramer. To adapt to the pentamer, the core structure is opened up at R1 and R2, the other parts need to be randomly generated and then optimized for the pentamer. Among them, nucleic acid single strands R10, R2, R4, R6 and R8 are specifically complementary paired with nucleic acid single strands R1, R3, R5, R7 and R9, respectively, without pairing with other nucleic acid single strands. The free energy of specific complementary pairing (ΔG°S) is less than −27.4 kcal/mol, while the free energy of non-specific pairing (ΔG°NS) is greater than −7.2 kcal/mol. In this way, the form of pentamers is the most stable in the reaction system. The specific implementation steps for optimizing the second conversion scheme of pentamer are described in FIG. 2; the annealing parameters are: annealing initial temperature To=50° C., annealing termination temperature Tf=0.12° C., and annealing temperature decay coefficient ΔT=0.98 (each attenuation of annealing initial temperature is 0.98 of the current value); the optimized constraint parameters are: pairing sequence length L=14 bases, dissociation temperature threshold Tm>52° C., pairing sequence free energy threshold: ΔG°S<−27.4 kcal/mol, and non-specific pairing free energy threshold: ΔG°NS>−7.2 kcal/mol.









TABLE 7





Nucleic acid sequences S1, S2, S3, S4 and S5


containing partial core structures


















S1
X1-Q1-AATCC-X2-R10-X3



S2
X1-R2-X2-TGAGC-Q2-X3



S3
X1-Q3-GCTCA-X2-CCGAA-Q4-X3



S4
X1-Q5-TTCGG-X2-ACTAT-Q6-X3



S5
X1-Q7-ATAGT-X2-GGATT-Q8-X3










The specific implementation steps for optimization are as follows:


Sequence initialization: according to the optimized constraint parameters, initialize the pairing sequence R9 and the set Q={Q1, Q3, Q5, Q7} of sequences with a length of 9 except for the core structure. Based on the principle of base complementation, R8, Q2, Q4, Q6 and Q8 are obtained and concatenated according to Table 7, and finally the initialized sequence set S={S1, S2, S3, S4, S5} is obtained.


Generation of new solution: the same as the generation of new solution in Example 1. The difference is that because this scheme uses a fixed core structure of tetramers, it is necessary to ensure that the fixed core structure remains partially unchanged during the generation of new solution. At the same time, constraint parameters need to be strictly followed, the dissociation temperature of the updated nucleic acid sequence pairing regions is greater than 52° C., and the free energy of the pairing regions (ΔG°S) is less than −27.4 kcal/mol.


Optimization judgment: the same as the optimization judgment in Example 1.


The optimized sequences of the second pentamer conversion scheme in Table 8_1 are obtained through this optimized algorithm, the specific pairing diagram thereof is shown in FIG. 11, and FIG. 12 shows the statistical line graph of the sum of free energy values of non-target pairing regions between sequences during this optimization process. In FIG. 13, Lane35 is the main band formed by sequences S1, S2, S3, S4 and S5. Although the pentamer shows somewhat trailing, the assembly effect is good.









TABLE 8_1





Initialized sequences and optimized sequences for the second


conversion scheme of pentamers

















Sequence




numbering
Initialized sequences
SEQ ID NO:





S1

ATGAGTGCGCAATCCAAATCGCCAGTCATG

265



CA






S2

ATGCATGACTGGCGAAAATGAGCGCTCGTT

266



GA






S3

ATCAACGAGCGCTCAAAACCGAAGTGCCA

267



ACT






S4

AAGTTGGCACTTCGGAAAACTATCGCGCG

268



ACT






S5

AAGTCGCGCGATAGTAAAGGATTGCGCAC

269



TCA





Sequence




numbering
Optimized sequences





S1

ATCACGCAGCAATCCAAATCGCCATCACAA

270



CG






S2

ACGTTGTGATGGCGAAAATGAGCACGAGC

271



CTT






S3

AAAGGCTCGTGCTCAAAACCGAAGGTTGC

272



ACT






S4

AAGTGCAACCTTCGGAAAACTATGCCGCT

273



CCA






S5

ATGGAGCGGCATAGTAAAGGATTGCTGCG

274



TGA
















TABLE 8_2







Initialized sequence and optimized sequence free energy


(ΔGNS°) matrix for the second conversion scheme of pentamers








Petramer initialized
Pentamer optimized sequence


free energy matrix
free energy matrix


















Sequence
S1
S2
S3
S4
S5
Sequence
S1
S2
S3
S4
S5





















S1
−13.79
−9.89
−9.89
−9.89
−9.89
S1
−3.61
−6.75
−7.04
−5.09
−6.75


S2

−16.23
−8.16
−9.89
−9.89
S2

−6.3
−6.61
−7.13
−6.91


S3


−16.23
−9.89
−9.89
S3


−7.05
−6.61
−6.68


S4



−20.25
−10.36
S4



−7.05
−7.04


S5




−20.25
S5




−5.09
















TABLE 8_3







Parameter indicators for pairing regions of


pentamer optimized sequences 2














Tm
Tm


Pairing
CG

° C.
° C.


sequences
%
ΔGS°
(TF)
(IDT)





R1
57.1
−28.33
56.3
48.8


R10
57.1
−28.55
57.4
48.8


R3
57.1
−28.1 
54.7
50  


R5
57.1
−28.14
53.8
48.3


R7
57.1
−28.49
54.5
49.4









In addition, Examples 1, 2, and 3 are repeated to obtain the nucleic acid single strand sequences and sets thereof shown in Tables 9-1, 9-2, and 9-3 (see above) for forming the trimer, tetramer, and pentamer complexes based on the complementary nucleic acid backbone.


EXAMPLE 4: COUPLING OF G-CSF AND L-DNA

The coupling of G-CSF and L-DNA selectively couples L-DNA with aldehyde group modification at the 5′ end to the N-terminal of G-CSF through a reductive amination reaction. Use dilution method or gel filtration chromatography and other methods to replace the buffer solution of G-CSF with acetate buffer solution (20 mM acetate, 150 mM NaCl, pH 5.0) and concentrate the sample to 30 mg/mL. Dissolve 100 OD (1 OD=33 μg) of L-DNA dry powder with aldehyde group modification at the 5′ end in 60 μL acetate buffer. Take 30 mg of sodium cyanide borohydride, dissolve in acetate buffer, and adjust the concentration to 800 mM. Take 50 μL, 60 μL, and 20 μL of G-CSF, L-DNA, and sodium cyanide borohydride at the above concentrations, mix them evenly, and incubate them at room temperature in dark for 48 hours. The samples before and after the reaction are verified by polyacrylamide gel electrophoresis for the effect of coupling reaction. The (L-DNA)-(G-CSF) coupling compound deviated significantly compared to the uncoupled G-CSF on the electrophoresis gel figure, and the coupling efficiency can reach 70%˜80% (FIG. 14).


EXAMPLE 5: PURIFICATION OF (L-DNA)-(G-CSF) COUPLING COMPOUND

The purification of the (L-DNA)-(G-CSF) conjugate is divided into two steps. The first step is to remove the unreacted G-CSF and the (L-DNA)2-(G-CSF) conjugate connected to two L-DNA strands using Hitrap Q HP (FIG. 15a). The reaction mixture obtained in Example 5 is diluted 10 times with the loading buffer of Q column, and then loaded, eluted 10 times the column volume with the loading buffer to remove the unreacted G-CSF, and then eluted 50 times the column volume with a 0-100% linear gradient to separate the (L-DNA)-(G-CSF) coupling compound and (L-DNA)2-(G-CSF) coupling compound, and identified the component type of each A280 absorption peak by polyacrylamide gel electrophoresis (FIG. 15b). Collect the (L-DNA)-(G-CSF) conjugate (containing unreacted nucleic acids). The second step is to use Hiscreen Capto MMC to remove unreacted nucleic acids, and ultimately isolate the highly purified (L-DNA)-(G-CSF) conjugate (FIG. 15c). Directly load the sample collected in step 1 onto the Hiscreen Capto MMC column, elute 10 times the column volume with the loading buffer to remove unreacted nucleic acids, and then elute the (L-DNA)-(G-CSF) conjugate with 100% elution buffer.


The purification conditions are as follows:


















Loading
Elution

Elution
Flow


Column
buffer
buffer
Gradient
volume
velocity







Hitrap
25 mM
25 mM acetate,
0-100%
50  custom-character
1


Q HP
acetate,
1 M NaCl,

custom-character
mL/min



pH 5.0
pH 5.0





Hiscreen
25 mM
150 mM
  100%
5  custom-character
2


Capto
acetate,
phosphate,

custom-character
mL/min


MMC
pH 5.0
150 mM NaCl,







pH 7.5









The purity of (L-DNA)-(G-CSF) conjugate sample obtained by two-step purification method is identified by 2% agarose gel electrophoresis (FIG. 15d). The gel diagram shows that there is only one nucleic acid band in the sample, and the (L-DNA)-(G-CSF) conjugate has a significant deviation compared to the uncoupled L-DNA on the gel diagram, indicating that the unreacted nucleic acid and (L-DNA)2-(G-CSF) conjugate has been cleaned out.


EXAMPLE 6: ASSEMBLY OF MONOVALENT, DIVALENT, AND TRIVALENT G-CSF COMPLEXES

Measure the nucleic acid concentrations of S1-G-CSF, S3-G-CSF, S4-G-CSF, S2, S3, and S4 using Nanodrop. Take an appropriate amount of the above components according to the structural design of the monovalent, divalent, and trivalent protein complexes, and mix them in a 1:1:1:1 molar ratio. After mixing, each assembly unit automatically completes assembly according to the principle of base complementation. Before and after assembly, samples are identified by polyacrylamide gel electrophoresis for assembly effect and sample purity (FIG. 16).


EXAMPLE 7: IN VITRO ACTIVITY EVALUATION OF G-CSF

M-NFS-60 cells (mouse leukemia lymphocytes/G-CSF dependent cells) are inoculated into the resuscitation culture medium (RPMI1640+10% FBS+15 ng/ml G-CSF+1× penicillin-streptomycin), and the cells are resuscitated at 37° C. and 5% CO2 conditions. When the cell density reaches 80%-90%, the cells are subcultured. After two or three times of subculture, the cells are inoculated into 96 well plates for cell plating. The cell plating experiment uses corning 3599 #96 well plates, with a cell plating density of 6000 cells per well. Gradient dilution is performed on different samples (GCSF, NAPPA4-GCSF, NAPPA4-GCSF2, NAPPA4-GCSF3) at working concentrations of (0.001, 0.01, 0.1, 1, 10, 100 ng/mL), with a final volume of 100 μL. Use PBS as a control. After incubating in a constant temperature incubator for 48 hours, add 10 μL of CCK8 solution to each well, and incubate the culture plate in the incubator for 1-4 hours, measure the absorbance at 450 nm using an enzyme-linked immunosorbent assay, and calculate the cell proliferation rate of different samples.





Cell proliferation rate (%)=[A(dosing)−A(0 dosing)]/[A(0 dosing)−A(blank)]×100

    • A (dosing): absorbance of wells with cells, CCK solution, and drug solution
    • A (blank): absorbance of wells with culture medium and CCK8 solution without cells
    • A (0 dosing): absorbance of wells with cells, CCK8 solution, without drug solution


The activity of a G-CSF linked L-DNA tetramer framework is evaluated using the above activity testing method, and it is found that the L-DNA tetramer framework has no effect on the activity of G-CSF (FIG. 17). Using the same activity testing method, it is found that bivalent and trivalent G-CSFs assembled with L-DNA tetramers do not have a negative impact on the activity of G-CSF (FIG. 18).


EXAMPLE 8: COUPLING AND PURIFICATION OF SM(PEG)2-PMO COUPLING COMPOUND

In this embodiment, phosphorodiamidate morpholino nucleic acid is used for the experiment. Specifically, the following four PMO single strand sequences are selected (from 5′ to 3′)











Strand 1 (PMO1):



SEQ ID NO: 275



5′-AGCAGCCTCGTTGAATCGCCAAGACACC-3′







Strand 2 (PMO2):



SEQ ID NO: 276



5′-AGGTGTCTTGGCGAAAGTTGCTCCGACG-3′







Strand 3 (PMO3):



SEQ ID NO: 277



5′-ACGTCGGAGCAACTAAGCGGTTCTGTGG-3′







Strand 4 (PMO3):



SEQ ID NO: 278



5′-ACCACAGAACCGCTATCAACGAGGCTGC-3′






The 5′ end is modified with an NH2 group, which is used for coupling the NHS active group of SM(PEG)2.


Dissolve PMO single strand containing 5′-terminal NH2 modification with phosphate buffer (50 mM NaH2PO4, 150 mM NaCl, pH 7.4) to prepare a mother solution with a final concentration of 1 mM. Dissolve SM(PEG)2 (linker molecule) powder with dimethyl sulfoxide (DMSO) and freshly prepare 250 mM of SM(PEG)2 mother solution. Add 10 to 50 times molar amounts of SM(PEG)2 mother solution to the PMO single strand mother solution, quickly mix them, and react at room temperature for 30 minutes to 2 hours. After the reaction is completed, add 10% volume of 1M Tris HCl (pH 7.0) to the reaction solution, mix them and incubate at room temperature for 20 minutes to quench the excessive SM(PEG)2 reaction. After incubation, the SM(PEG)2-PMO is purified using Hitrap Capto MMC. Unreacted SM(PEG)2 flows through the column without binding, while SM(PEG)2-PMO bound to the upper column is eluted with buffer (25 mM BICINE, 200 mM NH4Cl, 1M Arginine monohydrochloride, pH 8.5). The elution results are shown in FIG. 19a.


Analyze PMO samples before and after coupling using a positive ion mode of liquid chromatography-mass spectrometry. As shown in FIG. 20a and FIG. 20b, the results show that the final SM(PEG)2-PMO molecular weight is consistent with the theoretical value, and the coupling reaction efficiency is high.


EXAMPLE 9: PREPARATION OF NANO ANTIBODY MUTANTS

Cysteine mutations are introduced into the carboxyl end of nano antibodies for nucleic acid coupling. Optimize the gene sequence of the anti-HSA nano antibody into the yeast preferred codons, and then subclone it into the pPICZ alpha A plasmid. The amino acid sequence of the anti-HSA nano antibody is shown in SEQ ID NO: 279. To facilitate purification, the N-terminal of the nano antibody is labeled with His.


SEQ ID NO: 279, the amino acid sequence of the anti-HSA nano antibody:









HHHHHHAVQLVESGGGLVQPGNSLRLSCAASGFTFRSFGMSWVR


QAPGKEPEWVSSISGSGSDTLYADSVKGRFTISRDNAKTTLYLQMNSL


KPEDTAVYYCTIGGSLSRSSQGTQVTVSSGSC






Linearize the plasmid and transfer it to Pichia pastoris strain X33, and screen the high copy strain of the target gene using Zeocin concentration gradient YPD agar plate. Cultivate monoclonal strains using GMGY at 30° C. and 250 rpm conditions to obtain sufficient strains. Then, induce the expression and secretion of target nano antibodies using GMMY at 20° C. and 250 rpm conditions, and supplement with 1% methanol every 24 hours.


The expression yield of nano antibodies in laboratory grade glass flasks of the high copy-screened strain can reach 40-80 mg/L.


After SDS-PAGE identification and analysis, the culture supernatant after 72 hours of induction contains a large number of target nano antibody monomers and nano antibody dimers. Purify the nano antibodies in the culture supernatant using His labeled affinity column.


EXAMPLE 10: COUPLING AND PURIFICATION OF NANO ANTIBODY-PMO CONJUGATES

Dialyze the nano antibody sample eluted by His label affinity chromatography (Example 9) with a dialysis buffer containing a reducing agent (20 mM Tris, 15 mM NaCl, pH 7.4). During the dialysis process, the C-terminal thiol group is reduced while removing the small impurities such as free-SH groups. Mix the reduced nano antibody with SM(PEG)2-PMO single strand (prepared in Example 8) in a molar ratio of 1:1 to 2, and react at room temperature for 2 hours after mixing evenly.


As shown in FIG. 19b, the coupling efficiency can reach over 90% through SDS-PAGE identification.


Using His labeled affinity column to remove unreacted SM(PEG)2-PMO single strands, nano antibodies and nano antibody-PMO mixtures are collected.


Using Superdex™ 75 Increase 10/300 GL to separate nano antibodies and nano antibody-PMO, and the nano antibodies and nano antibody-PMO are effectively separated as shown in FIG. 21.


The nano antibodies before and after nucleic acid coupling are analyzed using a positive ion mode of liquid chromatography-mass spectrometry. As shown in FIGS. 20c and d, the final nano antibody-PMO molecular weight obtained is consistent with the theoretical value.


EXAMPLE 11: SELF-ASSEMBLY OF NAPPA-PMO DRUGS

Taking pmo NAPPA4-HSA (1) as an example, the self-assembly process of NAPPA-PMO drugs is introduced below.


Measure the concentrations of anti-HSA Nb-PMO1, PMO2, PMO3, and PMO4 respectively. Take an appropriate amount of the above components and preheat at 37° C. for 5 minutes, then mix them in a 1:1 molar ratio at 37° C. condition and incubate for 1 minute. Complete the assembly of pmo-NAPPA4-HSA (1).


Similarly, the pmo-NAPPA4-HSA (1,2,3) is assembled, and the required modules are anti-HSA Nb-PMO1, anti-HSA Nb-PMO2, anti-HSA Nb-PMO3, and PMO4.


Under low temperature conditions, the assembly of the samples is identified using SDS-PAGE. As shown in FIG. 22, the results show that the assembled samples have uniform bands.


EXAMPLE 12: VERIFICATION OF BINDING ACTIVITY OF NANO ANTIBODY-PMO MONOMERS AND ASSEMBLED NAPPA-PMO DRUGS

Taking pmo-NAPPA4-HSA (1) as an example, the binding ability of anti-HSA nano antibody, anti-HSA Nb-PMO1 monomer and assembled pmo-NAPPA4-HSA (1) to HSA protein (ACRO Biosystems, HSA-H5220) is tested by ELISA.


Each well of the 96-well ELISA plate is coated with 100 ng of HSA protein overnight at 4° C. Wash the plate with washing solution (PBS containing 0.05% Tween-20) and block it with blocking solution (PBS containing 3% BSA and 0.05% Tween-20), then add gradient diluted Anti-HSA nano antibodies, nano antibodies coupled with PMO, Anti-HSA Nb-PMO1, and assembled pmo-NAPPA4-HSA (1), and incubate at room temperature for 1 hour. After washing three times, add 1:5000 diluted horseradish peroxidase coupled rabbit-anti-camel VHH antibody (Genscript, A02016) and incubate at room temperature for 1 hour. After washing 3 times, add a tetramethylbenzidine substrate solution (Biyuntian, P0209) for development. Use a termination solution (Biyuntian, P0215) to quench the development. Use an enzyme-linked immunosorbent assay (MolecularDevices, SpectraMax i3x) to read the absorbance at 450 nm for each well and calculate the corresponding EC50.


The calculation results in FIG. 23 show that the EC50 of the binding of Anti-HSA nano antibodies, anti-HSA Nb-PMO1, and assembled pmo-NAPPA4-HSA (1) to HSA protein are 0.577 nM, 0.391 nM, and 0.529 nM, respectively. This indicates that the PMO coupling method of nano antibodies does not affect the binding activity of nano antibodies to corresponding antigens, and the PMO assembly method also does not affect the binding activity of nano antibodies to corresponding antigens.


EXAMPLE 13: RESISTANCE EXPERIMENT OF NAPPA-PMO DRUGS TO NUCLEASE DEGRADATION

PMO, as a nucleic acid derivative, can withstand the degradation of various nuclease. In order to verify whether the assembled NAPPA-PMO drugs can also tolerate nuclease or depolymerization, the following experiment is designed. Three common nuclease DNase I (Thermo Scientific, EN0523), T7 Endonuclease I (NEB, M0302S) and S1 Nuclease (Thermo Scientific, EN0321) are selected to incubate the PMO assembly sample pmo-NAPPA4-HSA (1) and the D-DNA assembly sample DDNA-NAPPA4 (control) for one hour at 37° C. The incubated pmo NAPPA4-HSA (1) and DDNA-NAPPA4 are analyzed using SDS-PAGE and 2% agarose electrophoresis, respectively.


As shown in FIG. 24, the pmo-NAPPA4-HSA(1) is not degraded by the three nuclease (left figure), while DDNA-NAPPA4 is completely degraded by DNase I and S1 Nuclease, and cut into shorter fragments by T7 Endonuclease I. Therefore, the experiment shows that NAPPA-PMO drugs can tolerate the degradation of common nuclease.


All documents mentioned in the present invention are incorporated by reference herein as if each document was incorporated separately by reference. Furthermore, it should be understood that after reading the foregoing teachings of the present invention, various changes or modifications can be made to the present invention by those skilled in the art and that these equivalents also fall in the scope of the claims appended to this application.

Claims
  • 1. A multimeric complex based on a complementary nucleic acid backbone, wherein the complex is a multimer formed by complexing n monomers having the complementary nucleic acid backbone, wherein each monomer is a polypeptide having a nucleic acid single strand, and n is a positive integer of 3-6; in the multimer, the nucleic acid single strand of each monomer and the nucleic acid single strands of the other two monomers form complementary double strands by means of base complementation, so as to form complementary nucleic acid backbone structures.
  • 2. The multimeric complex of claim 1, wherein the monomer has a structure of formula I: Z1-W   (I)wherein,Z1 is a polypeptide moiety;W is a nucleic acid single strand sequence; and“-” is a linker or bond.
  • 3. The multimeric complex of claim 2, wherein the nucleic acid sequence W has the structure shown in formula 1: X1-R1-X2-R2-X3   (1)wherein,R1 is a complementary base pairing region 1;R2 is a complementary base pairing region 2;Each of X1, X2, and X3 is independently not present or redundant nucleic acids; and“-” is a bond.
  • 4. The multimeric complex of claim 2, wherein the sequence of X2 is selected from the group consisting of: A, AA, AGA and AAA.
  • 5. The multimeric complex of claim 1, wherein the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 1-60 that form a trimer complex based on the complementary nucleic acid backbone.
  • 6. The multimeric complex of claim 1, wherein the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 61-140 that form a tetramer complex based on the complementary nucleic acid backbone.
  • 7. The multimeric complex of claim 1, wherein the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 141-240 that forms a pentamer complex based on the complementary nucleic acid backbone.
  • 8. The multimeric complex of claim 1, wherein the monomer sequence is any sequence or a sequence set thereof selected from the nucleic acid single strand sequences as shown in SEQ ID Nos: 275-278 that forms a tetramer complex based on the complementary nucleic acid backbone.
  • 9. A pharmaceutical composition comprising: (a) the multimeric complex based on the complementary nucleic acid backbone of claim 1; and(b) a pharmaceutically acceptable carrier.
  • 10. A nucleic acid sequence library, which comprises a nucleic acid sequence for forming the multimeric complex based on the complementary nucleic acid backbone of claim 1.
  • 11. The nucleic acid sequence library of claim 10, wherein the nucleic acid sequence W has the structure shown in formula 1: X1-R1-X2-R2-X3   (1)wherein,R1 is the complementary base pairing region 1;R2 is the complementary base pairing region 2;Each of X1, X2, and X3 is independently not present or redundant nucleic acids; and“-” is a bond.
  • 12. A method of preparing the multimeric complex of claim 1 or a pharmaceutical composition comprising the multimeric complex of claim 1, comprising the step of using the nucleic acid sequence of a nucleic acid sequence library, which comprises a nucleic acid sequence for forming the multimeric complex based on the complementary nucleic acid backbone of claim 1, to prepare the multimeric complex of claim 1 or a pharmaceutical composition comprising the multimeric complex of claim 1.
  • 13. A method of determining a nucleic acid single strand sequence for forming a multimeric complex based on a complementary nucleic acid backbone, comprising steps of: (a) setting annealing algorithm parameters:setting the initial annealing temperature, annealing termination temperature, and annealing temperature attenuation coefficient ΔT;setting optimized constraint parameters:(i) the number n of the nucleic acid single strand, preferably a positive integer of 3-6;(ii) the length L of the pairing sequence, preferably the L is of 12-16 bases;(iii) the dissociation temperature threshold Tm of the pairing region;(iv) the free energy threshold ΔG°S of the specific pairing region sequence;(v) the free energy threshold ΔG°NS of the non-specific pairing;(vi) the connecting element X2, preferably A, AA, and AAA;(vii) the dissociation temperature threshold Tm-H of the secondary structure (hairpin);(viii) the CG proportion PCG in the pairing sequence, preferably the range of PCG is [0.4,0.6);(ix) optionally, for n=4, using a symmetric sequence to initialize a sequence set S={S1, S2, . . . , Sn} according to the above parameters;(b) calculating the objective function value E0 of the set S of the previous step, that is, calculating the sum of the non-specific pairing free energies (ΔG°NS) between sequences and of the sequence itself, while obtaining the non-specific pairing free energy matrix Cn×n. searching the Si and Sj(1≤i≤n, 1≤j≤n) corresponding to the minimum value in the upper triangular matrix thereof, randomly selecting Si or Sj for an updated operation according to the non-specific pairing free energy of the Si and Sj ΔG°NS(Si, Sj), and then obtaining a new nucleic acid sequence, thereby obtaining a updated sequence set S′;(c) determining whether the sequences in the set S′ of the previous step meet the optimized constraint parameter conditions set in step (a), verifying the following parameters, including the dissociation temperature Tm of the specific pairing region, the free energy ΔG°S of the specific pairing region sequence, the dissociation temperature Tm-H of the secondary structure and the CG proportion PCG. If the above parameters meet the constraint conditions, the step (d) is proceeded; otherwise, the step (c) is repeated. If the step (b) is performed 15 times continuously at a certain annealing temperature without obtaining the S′ that meets the conditions, then the set S becomes the set S′ and the next step is proceeded to prevent a dead cycle;(d) calculating the objective function value E1 of the set S′ of the previous step, and comparing E0 with E1. If E1≥E0, it indicates that the non-specific pairing free energy has been optimized, and the sequence set S′ becomes the sequence set S. If E1<E0, it indicates that the non-specific pairing free energy has not been optimized, and in this case, it is necessary to determine whether to accept the set S′ as S according to the Metropolis criterion; and(e) the annealing temperature is attenuated according to the attenuation coefficient ΔT set in the step (a), and the steps (b), (c), and (d) are repeated for the S of the previous step, which is the Monte Carlo-based annealing algorithm, until the annealing temperature reaches the annealing termination temperature. The S={S1, S2, . . . , Sn} of the previous step becomes the nucleic acid single strand sequence for forming the multimeric complex based on the complementary nucleic acid backbone.
  • 14. A nucleic acid single strand sequence set for forming a multimeric complex based on a complementary nucleic acid backbone, which is determined using the method of claim 13.
  • 15. The nucleic acid single strand sequence set of claim 14, wherein the set is selected from the group consisting of: (S1) a nucleic acid single strand sequence for forming a trimer complex based on the complementary nucleic acid backbone:
Priority Claims (1)
Number Date Country Kind
202011340377.3 Nov 2020 CN national
PCT Information
Filing Document Filing Date Country Kind
PCT/CN2021/133254 11/25/2021 WO