Demonstrate
1 Abstract
Chronic PKA phosphorylation of RyR2 has been shown to increased diastolic SR Ca2+ leak and lead to cardiac dysfunction. Since the change of phosphorylation level of RyR2 is a biomarker of failing heart, we attempted to verify the hypothesis that intracellular gene delivery of a RyR2 targeting phosphorylation site-specific nanobody could preserve contractility of failing myocardium. In present study, we acquired the RyR2-specific nanobodies from a phage display library which are variable domains of camellidae heavy chain-only antibodies (VHH). One of the monoclonal nanobodies, AR185, inhibiting RyR2 phosphorylation in an in vitro assay was then chosen for further investigation. We investigated the potential of adeno-associated virus (AAV)-9-mediated cardiac expression of AR185 against post-ischemic heart failure. Adeno-associated virus gene delivery elevated the intracellular expression AR185 protein in the ischemic heart failure model of rats, and this treatment normalized the systolic and diastolic dysfunction of the failing myocardium in vivo and in vitro by reversing myocardial Ca2+ handling. Furthermore, AR185 gene transfer to failing cardiomyocytes reduced the frequency of sarcoplasmic reticulum (SR) calcium leak, thereby restoring the attenuated intracellular calcium transients and SR calcium load. Moreover, AR185 gene transfer inhibited PKA phosphorylation of RyR2 in failing cardiomyocytes. Our results provided strong pre-clinical experimental evidence of the cardiac expression of RyR2 nanobody with AAV9 vectors as a promising therapeutic strategy for ischemic heart failure.
2 Introduction
In most countries, heart problems are the leading cause of morbidity and mortality(1). Heart failure is the end stage of most heart diseases, with approximately 38 million heart failure patients worldwide(2). Despite the use of drugs, implanted cardiac assist devices and surgical treatments, many patients' conditions will still irreversibly deteriorate, eventually being difficult to control and rescue.(3) Therefore, studying the molecular mechanism of cardiovascular pathology, developing new therapeutic strategies and preparations have the great significance for the prevention and treatment of heart failure. No matter what the cause, heart failure patients in the late stage of the disease have a common feature that calcium circulation is abnormal in the cardiomyocytes. A hallmark of failing cardiomyocytes is the change of excitatory contraction coupling, including the reduced amplitude of Ca2+ transients, delayed onset and decay kinetics of Ca2+ transients. These changes eventually lead to decreased contractility, delayed contraction, and reduced diastole. In addition, a rare spontaneous Ca2+ release events (calcium leaks) in the resting period of healthy cardiomyocytes occurs more frequently. All the change of Ca2+ handling is attributed to impaired function of RyR2, SERCA2a, Na+-Ca2+ exchanger (NCX).
Changes in expression and function of these ion channel proteins and their negative effect after MI have been reported in many studies(4, 5), which directly attenuate the cardiac contractility. Therefore, preventing these detrimental effects of myocardial infarction, as well as the effects of ischemic injury due to the changes in expression or function of calcium-channel proteins, is an optimal treatment strategy for cardio-protection.
The activity of RyR2 function is regulated by various mechanisms, and many factors will affect it. The post-translational regulation abnormality of RyR2 is the main cause of the dysfunction of RyR2 protein. Phosphorylation is a crucial post-translational modification of RyR2 protein.(6, 7) The "calcium leakage" caused by RyR2 hyperphosphorylation is considered to be an important pathological mechanism for myocardial injury and heart failure development.(8, 9) The study found that the expression of catecholamines in the blood of chronic heart failure patients was up-regulated, and intracellular protein kinase A (PKA) levels and activities continued to increase. After PKA is continuously activated, hyperphosphorylation of RyR2 protein causes an increase of the dissociation rate of FKBP12.6 and RyR2, resulting in a change of RyR2 structure, which increases the sensitivity of RyR2 to calcium ions, resulting in a small amount of calcium ions can stimulate the release of RyR2 channels, causing the calcium leakage in the resting stage of cardiomyocytes. The earliest article demonstrating changes in ryr2 phosphorylation in patients with heart failure was published by Marx et al. in 2000(10). Marx et al.(10), Wehrens et al(11). and Shan et al.(12) proposed a hypothesis: sympathetic excitation in heart failure, activation of PKA-mediated RyR2 S2808 Hyperphosphorylation, while inhibiting the binding of FKBP12.6 to RyR2, increases the probability of RyR2 opening.
But this assumption is currently controversial. The biggest controversy is whether RyR2 S2808 hyperphosphorylation plays a decisive role in heart failure. Wehrens et al.(11) believe that RyR2 S2808 phosphorylation plays a critical role in the experimental myocardial ischemia-induced heart impaired process. They established a heart failure model which induced by myocardial infarction in RyR2 S2808 knockout mice and wild-type mice. Heart function tests were performed 4 weeks after ischemia. The results showed that RyR2 S2808 knockout mice showed significant improvement in cardiac function such as ejection fraction, shortening fraction, maximal rate of left ventricular pressure. At the same time, they also found that RyR2 of RyR2 S2808 knockout mice could not be phosphorylated by PKA, and PKA lost regulation of FKBP12.6, resulting in FKBP12.6 not being able to dissociate from RyR2. However, Zhang et al.(13) and Houser et al.(14) believe that PKA-mediated hyperphosphorylation of RyR2 S2808 does not alter myocardial contractility and does not improve symptoms of heart failure and arrhythmia. They used a similar experimental protocol and found no improvement in the heart function of the mouse knocked out by RyR2 S2808. However, the biggest problem with Hauser et al. is that the number of samples is too small. Inhibition of RyR2 phosphorylation may be one of the effective treatments of ischemic heart failure.
Here, we hypothesis that targeting RyR2 using anti-phosphorylation agents may improve treatment efficacy. We identified a camel single-domain antibody to RyR2 that have the ability to inhibit PKA dependent S2808 phosphorylation in vitro. To evaluate its potential effect in the treatment of heart failure, an adeno-associated virus (AAV) based intracellular antibody delivery strategy were adopt to achieve cardiac-specific gene-therapy and demonstrated therapeutic effect both in cell-based assays and in vivo models.
3 Results
Generation of Anti-RyR2 nanobodies that specifically inhibits the Phosphorylation of RyR2 S2808
We first obtained and purified RyR2 from rat heart by using GST-fused FKBP12 as the published strategies described. To construct the camel VHH library, blood samples of 30 non-immunized, four-year-old male Bactrian camel were collected. B lymphocyte cDNA encoding VHHs was used to construct a phage display VHH library that consisted of approximately 3×108 individual colonies. VHH gene corresponded to the size of insert of over 98% colonies. For confirming the heterogeneity of the individual clones from the library, we sequenced fifty randomly selected clones, and each clone showed a distinct VHH sequence.
In order to select nanobodies with specific ability to bind RyR2, bio-panning was performed with immobilized RyR2 protein. After the third round of panning, the result showed an obvious enrichment of phage particles that carried RyR2-specific VHH (Fig. 1A). Phage clones exhibited increased binding to RyR2 after the second round of panning. During four rounds of panning there was no phage clone that was found binding to BSA (Fig. 1B). VHH fragments of 300 individual colonies that were randomly chosen were expressed in an ELISA for screening colonies which bound to RyR2. Among these clones, 276 antibody fragments specifically bound to RyR2. One antibody fragment which did not bind to RyR2 was choose as a negative control, termed as VHH-AR117. To obtain antibodies that functionally inhibit of RyR2 phosphorylation, each of the antibody fragments was tested for its effect in an ELSA based RyR2 phosphorylation assay. 4 antibody fragments were potent inhibitors of RyR2 phosphorylation. The complementary determining regions (CDRs) were confirmed by sequence analysis and the result revealed that there was only one unique clone in this panel of antibody fragments, termed as VHH-AR185. To investigating the basis of dephosphorylation of RyR2 by VHH-AR185, the binding affinity of VHH-AR185 to RyR2 was measured by surface plasmon resonance. VHHs were purified for these experiments by expressing and secreting from the E. coli cytosol. As shown in Fig. 1D, the affinity (KD) of VHH-AR185 to RyR2 was estimated to be 1.93 nM. The result of affinity studies likely explained the inhibition of RyR2 phosphorylation due to the extremely slow dissociation rate of VHH-AR185 from RyR2.
Fig.1 balablablabla
The interaction of VHH-AR185 to RyR2 in the cytoplasm of eukaryotic cells was examined by co-immunoprecipitation experiments. VHH-AR185 and RyR2 were expressed in neonatal cardiomyocytes cells and the lysates of transfected cells were detected. As the result in fig. S1, anti-his antibody was able to efficiently co-precipitate RyR2 from the cells that expressed VHH-AR185-HIS, but could not co-precipitate RyR2 from cells expressing VHH-AR117-HIS. Conversely, anti-RyR2 antibody was able to co-precipitate VHH-AR185-HIS with RyR2, but not VHH-AR185-HIS. This result indicated that the VHH-AR185 could maintain its antigen binding ability in the cytoplasm and fold as a soluble protein.
Fig.s1 balablablabla
To identify the epitopes recognized by AR185, phage clones were isolated by panning the PhD.-7 phage display peptide library with AR185. Three rounds of selection were performed, and, at each round, the library was pre-cleared with a control AR177 nanobody. After the third round of panning, the binding of the isolated phage clones to AR185 was determined by ELISA. Sequence analysis of AR185-positive phage clones identified five and six distinct amino acid sequences, respectively (fig. S2A). Alignment of these sequences revealed the consensus motifs DKLAC, which could be aligned with the (2725) DKLAN (2729) sequence located at P2 Domain of RyR2 (fig. S2B).
Fig.s1 balablablabla
Intrabody AR185 rescues cardiac function and reverses remodeling in failing rat myocardium in vivo
We constructed an AAV9 vector containing a VHH-AR185 expression between the two AAV2 inverted terminal repeats and the vector was pseudotyped with a capsid of AAV serotype 9, termed as AAV9.AR185. VHH-AR117 were also constructed as a negative control, termed as AAV9.AR117. To access cardiac expression of VHH, we used the method of adding “self-cleaving” T2A peptide to co-expressed a GFP reporter downstream of VHHs (Fig. 2A). We used the HEK-293 cells expression of different AAV9 particles in vitro and Transmission electron microscope was used to access the AAV9 particles (Fig. 2B). Next, we evaluated the efficiency of gene expression delivered by AAV in vivo. A dosage of AAV9.AR185 particles was delivered to each rat at 1×1012 genome containing particles (gcp), whereas AAV9.AR117 vector was given to the control treated group (n = 5) at the same dosage. After four weeks, we removed all organs from the sacrificed rat, weighed, and assayed treated tissue for fluorescence intensity. Efficiency of gene expression and ability of targeting were evaluated by the ratio of fluorescence intensity to mass of tissue under fluorescence microscope (Fig. 2C and D).
Fig.2 balablablabla
To explore the therapeutic potential of VHH, we chose the mode of ischemic heart failure induced by coronary artery ligation for this study. Following the ligation operation, rats were divided into different groups as described in Methods and received control virus (AAV9.AR117), AAV9.AR185 treatment or saline (HF) (n=7-8). The sham-operated animals (Sham) were used as healthy controls. Nine weeks after ligation operation and injection of AAV particles, LV dimensions in the short-axis view was measured by cardiac echo and we also calculated and analyzed the value of ejection fraction and fractional shortening. Our data shows that rats of HF group and AAV9.AR117 group exhibited progressive cardiac dysfunction and LV enlargement, while AAV9.AR185-treated animals showed significant improvement. Moreover, Ejection Fraction and fractional shortening was markedly improved in AAV9.AR185 group compared with HF group and AAV9.AR117 group (Fig. 3A). To determine whether AAV9.AR185 treatment prevented adverse remodeling of the heart after MI, Masson trichrome staining of cardiac sections was performed to measure cardiac fibrosis (Fig. 3B). Whereas there was a significant increase in the development of cardiac fibrosis in Rats of HF group and AAV9.AR117 group after HF, whereas the amount of fibrosis was significant reduced in AAV9.AR185-treated animals. Additionally, HF rat and AAV9.AR117 treated rat had development of a significant increase of heart weight to body weight ratios (HW/BW) after MI compared with sham-operated rat, which is indicative of cardiac remodeling in the context of congestive HF(Fig. 3C). In contrast, there was no significant increase in HW/BW ratio after MI in AAV9.AR185-treated rat compared with sham-operated rat. Sarcomeres and mitochondria were the most important index for analysis of ultra-structures of cardiomyocytes from left ventricle that were observed by transmission electron microscopy (Fig. 3D). In the AAV9.AR185 treated and Sham groups, myofilaments were neatly arranged, sarcomeres were intact and Z lines were clear. Conversely, in the HF and AAV9.AR117 groups, MI leaded to disordered arrangement of sarcomeres, dissolution of myofilaments, and frequent vacuoles. In both HF and AAV9.AR117 groups, a lot of mitochondria were swollen and even ruptured, and the separated mitochondrial cristae frequently appeared. The mitochondria in Sham group were well shaped, and the cristae of the mitochondria were obvious and tightly packed. The observations of mitochondria were improved in the AAV9.AR185 treated group compared with AAV9.AR117 treated group. Comprehensively considering the alteration of cardiac function and changes in structure of different groups, the TEM images further support that VHH-AR185 had therapeutic effect in treating heart failure.
Fig.3 balablablabla
We next accessed the contractile kinetics of isolated LV cardiomyocytes(Table1). When cardiomyocytes were field-stimulated at a frequency of 1 Hz, HF and AAV9.AR117 treated myocytes had significantly slower velocities of shortening and relengthening in than AAV9.AR185 treated myocytes. Fractional shortening of myocytes that were isolated from HF and AAV9.AR117 treated animals also decreased, and time to 50% peak shortening (TPS50%) and time to 50% relengthening (TR50%) became longer. AAV9.AR185 treatment protected cardiomyocytes contractility reserve from the impairment induced by MI. However, only the index of TR50% in myocytes from AAV9.AR185 treated animals returned to a level similar to those of sham operated animals.
Table 1 Contractile properties of cardiomyocytes from groups
Sham | HF | AAV-AR117 | AAV-AR185 | |
Fractional shortening (%) | 10.35±0.53** | 7.03±1.00 | 7.12±1.13 | 9.49±0.49* |
+dl/dt (μm/s) | 3.37±0.13** | 2.30±0.16 | 2.26±0.18 | 2.96±0.18* |
–dl/dt (μm/s) | 3.02±0.23** | 2.01±0.15 | 2.04±0.21 | 2.78±0.15* |
TPS50% | 64.75±2.36** | 75.06±5.22 | 74.25±5.58 | 69.90±3.88* |
TR50% | 155.70±9.56# | 196.37±10.47 | 186.67±15.22 | 162.18±8.91* |
AR185 gene therapy restores cardiomyocyte and SR calcium handling in failing myocardium
We used laser scanning confocal microscopy recorded the fluorescence intensity to measure the sarcoplasmic reticulum Ca2+ content of cardiomyocytes from different groups by incubation in the fluorescent dye Fluo-5N/AM. As shown in Fig. 4A, basal sarcoplasmic reticulum Ca2+ contents in HF and AAV9.AR117 treated animals were lower than in AAV9.AR185 treated and sham-operated animals.
Additionally, we measured amplitude of calcium transient by incubation in Fluo-4/AM and caffeine perfusion. The representative colorful images in Fig. 4B and Table 2 show line-scan results of evoked Ca2+ transients from Shams, HFs, AAV9.AR117s and AAV9.AR185s. When challenged with 20 mM caffeine, less Ca2+ was released from the SR of myocytes from AAV9.AR117 group compared with myocytes from AAV9.AR185 treated rats. The results also showed there was significantly reduction of the amplitude of Ca2+ transients in the HFs and AR117s compared to that in AR185s. Therefore, the decrease in Ca2+ transient amplitude may be the causative factor of the impairment in SR Ca2+ load. Rate of Ca2+ rise also was significantly slower in HF and AAV9.AR117 myocytes than in AAV9.AR185 treated myocytes (Fig. 4C). AAV9.AR185 treatment increased peak of amplitude of evoked Ca2+ release and rate of Ca2+ rise during Ca2+ release.
Fig.4 balablablabla
Table 3 shows representative line-scan images of Ca2+ release during the resting stage of cardiomyocytes from Sham (A), HF (B), AAV9.AR117 (C), and AAV9.AR185 animals (D). The data showed that the frequency of Ca2+ release was significantly higher and Ca2+ sparks occurred frequently in the HF and AAV9.AR117 group compared with the AAV9.AR185 group. The duration of Ca2+ sparks in HF and AR117 myocytes were similar to those in Sham and AR185 myocytes, but the Ca2+ rise rate of sparks was slower, fluorescence intensity of Ca2+ sparks was decreased and T50 decay was longer.
Table 2 Ca2+ transient in cardiomyocytes from groups
Sham | HF | AAV-AR117 | AAV-AR185 | |
Peak Ca2+ amplitude | 45.73±3.32 | 13.06±1.39 | 12.25±1.28 | 23.11±2.12 |
Rate of Ca2+ rise (/sec) | 687.12±73.95 | 87.19±9.35 | 83.90±6.83 | 334.28±59.11 |
Table 3 Spontaneous Ca2+ release in¬¬ cardiomyocytes from groups
Sham | HF | AAV-AR117 | AAV-AR185 | |
Frequency of Ca2+ sparks (per 50um/sec) | 0.7±0.34 | 12.6±1.38 | 12.3±1.25 | 1.1±0.34 |
Duration of Ca2+ sparks(msec) | 44.7±3.1 | 48.1±4.2 | 49.4±3.7 | 46.9±4.0 |
Peak Ca2+ amplitude | 37.5±4.3 | 18.4±2.9 | 18.6±2.1 | 32.3±3.6 |
Rate of Ca2+ rise (/msec) | 31.32±4.05 | 4.67±0.54 | 4.58±0.36 | 27.15±3.59 |
T50 decay (msec) | 7.75±0.15 | 10.85±0.57 | 10.57±0.70 | 8.46±0.19 |
VHH-AR185 inhibits phosphorylation of RyR2 S2808 in failing hearts
To examine whether AAV9.AR185 treatment results in dephosphorylation of RyR2 and in vivo, cardiomyocyte lysates were further subjected to ELISA analysis, our data shows treatment with AAV9.AR185 significantly reduced the level of pRyR2 (S2808) in the cardiomyocytes compared with HF group and AAV9.AR117 treatment (p=0.0003, Dunnett’s test). Moreover, immunohistochemical analysis of the heart tissues in different treatment group also revealed that an increased accumulation of RyR2 phosphorylation was also observed in the AAV9.AR117 treated group, AAV9.AR185 treatment decreased the level of pRyR2 stain of cells in the myocardium, which indicated that VHH185 has blockage effect of RyR2 phosphorylation. Together, these data demonstrate that AAV9.AR185 treatment leads to inhibition of the RyR2 phosphorylation in vivo. (Fig. 5).
Fig.4 balablablabla
4 Discussion
In this experiment, we blocked the PKA regulated phosphorylation site(S2808) by exogenously expressing a intracellular antibody specifically binding to RyR2 in cardiomyocytes. After AAV.AR185 treatment, the Ca2+ handing properties of RyR2 in rat hearts isolated form treatment group were very similar to that from Sham group. When the β-adrenergic signaling pathway was activated in vivo, the PKA-regulated site of RyR2 could not bind to phosphate residues, and the phosphorylation of RyR2-S2808 maintained at a low level. At the organic level, AR185 antibody can improve the contractile function of the failing heart, significantly increase the cardiac ejection fraction of heart failure rats, and protect the heart structure from the effects of poor remodeling (Fig 3). At the cellular level, AR185 antibody can reduce the calcium leakage of cardiomyocytes in heart failure rats, maintain the calcium capacity in the sarcoplasmic reticulum of cardiomyocytes, restore calcium homeostasis, and protect the contractility of cardiomyocytes (Fig 4 and Table 1-3). These data is consistent with the previous transgenic animal experiment that down-regulated phosphorylation level of RyR2 S2808 is associated with improved contractility of cardiomyocytes, indicating that S2808 is one of the most important targets in the pathological state of ischemic heart failure.(11)
To our knowledge, the evidence for efficiency and safety of using antibodies to treat heart disease was limited. Nanobodies are single domain antibodies consisting of the heavy chain variable domain (VHH) in the camelid family which lacks the light chain. Currently, a variety of nanobodies have entered the clinical research stage.(15-18) Compared with traditional antibodies, nanobodies have the advantages of low molecular weight, high affinity, high stability, low immunogenicity and strong penetrability.(19) Based on the characteristics of nanobodies and VHH, the use of adeno-associated virus vectors to mediate nanobody treatment of heart failure has great potential. In this study, we successfully expressed AR185 nanobodies which specifically bind to RyR2 in rat cardiomyocytes. Our experimental data demonstrate that intracellular antibody treatment is effective in heart disease rats and does not present a significant safety risk.
There are several limitations in our study. We provide evidence that targeting RyR2 with AR185 can be achieved, but rats with heart failure may not fully recapitulate RyR2 in human disease, the experiment data are acquired from a small number of animals, and limited by the shorter observation time, no rats died of heart failure before sacrifice. Furthermore, the mechanisms responsible for therapeutic improvement of AR185 against RyR2 have not been well characterized yet. Hence, these findings will need further validation.
The therapy strategy of AAV9-mediated intracellular antibody which takes advantage of the high specificity and affinity of biomacromolecules can achieve the purpose of specificity regulation of RyR2 and treatment of heart failure. RyR2 is a very promising therapeutic target protein of treatment for heart failure, and intracellular antibody technology in gene therapy is considered as a promising approach. The targets and techniques are worthy of further validation and exploration by researchers.
References
- W. Chen et al., Outline of the report on cardiovascular diseases in China, 2014. 18, F2 (2016).
- B. Ziaeian, G. C. J. N. R. C. Fonarow, Epidemiology and aetiology of heart failure. 13, 368 (2016).
- C. W. Yancy et al., 2016 ACC/AHA/HFSA Focused Update on New Pharmacological Therapy for Heart Failure: An Update of the 2013 ACCF/AHA Guideline for the Management of Heart Failure. 68, 1476-1488 (2016).
- J. S. Hulot, K. Ishikawa, R. J. J. E. H. J. Hajjar, Gene therapy for the treatment of heart failure: promise postponed. 37, 1651-1658 (2016).
- M. Luo, M. E. J. C. R. Anderson, Mechanisms of Altered Ca^2^+ Handling in Heart Failure. 113, 690-708 (2013).
- F. V. J. J. o. B. C. Petegem, Ryanodine Receptors: Structure and Function. 287, 31624-31632 (2012).
- G. W. Cho, F. Altamirano, J. A. J. B.-M. B. o. D. Hill, Chronic heart failure: Ca2+, catabolism, and catastrophic cell death. 1862, 763-777 (2016).
- J. Roussel et al., Palmitoyl-carnitine increases RyR2 oxidation and sarcoplasmic reticulum Ca 2 + leak in cardiomyocytes: Role of adenine nucleotide translocase. 1852, 749-758 (2015).
- M. Munkvik et al., Training effects on skeletal muscle calcium handling in human chronic heart failure. 42, 847-855 (2010).
- S. O. Marx et al., PKA Phosphorylation Dissociates FKBP12.6 from the Calcium Release Channel (Ryanodine Receptor). 101, 365-376 (2000).
- X. H. Wehrens et al., Ryanodine receptor/calcium release channel PKA phosphorylation: a critical mediator of heart failure progression. Proc Natl Acad Sci U S A 103, 511-518 (2006).
- J. Shan et al., Role of chronic ryanodine receptor phosphorylation in heart failure and β-adrenergic receptor blockade in mice. 120, 4375-4387 (2010).
- H. Zhang et al., Hyperphosphorylation of the cardiac ryanodine receptor at serine 2808 is not involved in cardiac dysfunction after myocardial infarction. 110, 831 (2012).
- S. R. J. C. R. Houser, Role of RyR2 phosphorylation in heart failure and arrhythmias: protein kinase A-mediated hyperphosphorylation of the ryanodine receptor at serine 2808 does not alter cardiac contractility or cause heart failure and arrhythmias. 114, 1320-1327 (2014).
- S. Steeland, R. E. Vandenbroucke, C. J. D. D. T. Libert, Nanobodies as therapeutics: big opportunities for small antibodies. 21, 1076-1113 (2016).
- F. Peyvandi et al., Caplacizumab reduces the frequency of major thromboembolic events, exacerbations, and death in patients with acquired thrombotic thrombocytopenic purpura. 15, 1448-1452 (2017).
- T. Dörner et al., in European Congress of Rheumatology, 14–17 June. (2017), pp. 575.572-575.
- Safety and efficacy of multiple ascending doses of subcutaneous M1095, an anti-interleukin-17A/F bispecific nanobody, in patients with moderate-to-severe psoriasis %J Retour Au Numéro. (2017).
- L. D. Wang, D. X. Feng, S. M. Zhang, B. M. J. L. i. B. University, Research Progress of Nanobody. (2016).