Research article

The inhibition of colon cancer development by black rice bran on BALB/C Mice

  • Black rice bran (BRB) is well-known for its high antioxidant activity and its pivotal role in preventing colon cancer. The present study aims to investigate the effects of BRB administration on BALB/C mice induced with azoxymethane (AOM) and dextran sodium sulphate (DSS). The 24 mice were divided into three groups: the group induced by colon cancer (C+), the group induced by cancer and given the BRB diet (C+BRB), and the normal group (C−). Both the C− and C+ groups were given a standard AIN-93 M diet containing cellulose fiber. After 16 weeks, the mice were anesthetized, and the colonic tissue was identified for nodule distribution, histopathological observation, and mRNA expression analysis of proliferating cell nuclear antigen (PCNA), caspase-3, and caspase-8 genes using qRT PCR technique. Nodule distribution in the C+BRB group showed a significant decrease compared to the C+ group, with 1.65 ± 0.71 nodule/cm2 vs. 5.73 ± 2.93 nodule/cm2. Then, the colon weight was significantly decreased in the C+BRB group, at 0.19 ± 0.04 (g) compared to 0.25 ± 0.03 (g) in the C+ group. Also, the BRB diet in the C+BRB group significantly decreased PCNA mRNA expression compared to the C+ group, with values of 0.58 ± 0.09-fold change vs. 5.22 ± 0.80-fold change. Conversely, increased the mRNA expression of caspase-3 (0.91 ± 0.20-fold change vs. 0.36 ± 0.15-fold change) and caspase-8 (0.51 ± 0.18-fold change vs. 0.13 ± 0.31-fold change). In conclusion, administration of BRB inhibited the rate of cancer development by suppressing cancer cell proliferation and inducing apoptosis.

    Citation: Slamet Budijanto, Yeni Kurniati, Lilis Nuraida, Fitriya Nur Annisa Dewi, Bambang Pontjo Priosoeryanto, Nancy Dewi Yuliana, Ardiansyah, Uus Saepuloh, Safrida, Hitoshi Shirakawa. The inhibition of colon cancer development by black rice bran on BALB/C Mice[J]. AIMS Agriculture and Food, 2024, 9(3): 789-800. doi: 10.3934/agrfood.2024042

    Related Papers:

    [1] Urszula Ledzewicz, Heinz Schättler . The Influence of PK/PD on the Structure of Optimal Controls in Cancer Chemotherapy Models. Mathematical Biosciences and Engineering, 2005, 2(3): 561-578. doi: 10.3934/mbe.2005.2.561
    [2] Yuyang Xiao, Juan Shen, Xiufen Zou . Mathematical modeling and dynamical analysis of anti-tumor drug dose-response. Mathematical Biosciences and Engineering, 2022, 19(4): 4120-4144. doi: 10.3934/mbe.2022190
    [3] Bertin Hoffmann, Udo Schumacher, Gero Wedemann . Absence of convection in solid tumors caused by raised interstitial fluid pressure severely limits success of chemotherapy—a numerical study in cancers. Mathematical Biosciences and Engineering, 2020, 17(5): 6128-6148. doi: 10.3934/mbe.2020325
    [4] Colette Calmelet, John Hotchkiss, Philip Crooke . A mathematical model for antibiotic control of bacteria in peritoneal dialysis associated peritonitis. Mathematical Biosciences and Engineering, 2014, 11(6): 1449-1464. doi: 10.3934/mbe.2014.11.1449
    [5] Michele L. Joyner, Cammey C. Manning, Whitney Forbes, Michelle Maiden, Ariel N. Nikas . A physiologically-based pharmacokinetic model for the antibiotic ertapenem. Mathematical Biosciences and Engineering, 2016, 13(1): 119-133. doi: 10.3934/mbe.2016.13.119
    [6] Guochun Lou, Jie Dong, Jing Du, Wanyuan Chen, Xianglei He . Clinical significance of lymph node micrometastasis in T1N0 early gastric cancer. Mathematical Biosciences and Engineering, 2020, 17(4): 3252-3259. doi: 10.3934/mbe.2020185
    [7] Urszula Ledzewicz, Heinz Schättler, Mostafa Reisi Gahrooi, Siamak Mahmoudian Dehkordi . On the MTD paradigm and optimal control for multi-drug cancer chemotherapy. Mathematical Biosciences and Engineering, 2013, 10(3): 803-819. doi: 10.3934/mbe.2013.10.803
    [8] Michele L. Joyner, Cammey Cole Manning, Whitney Forbes, Valerie Bobola, William Frazier . Modeling Ertapenem: the impact of body mass index on distribution of the antibiotic in the body. Mathematical Biosciences and Engineering, 2019, 16(2): 713-726. doi: 10.3934/mbe.2019034
    [9] Peter Hinow, Philip Gerlee, Lisa J. McCawley, Vito Quaranta, Madalina Ciobanu, Shizhen Wang, Jason M. Graham, Bruce P. Ayati, Jonathan Claridge, Kristin R. Swanson, Mary Loveless, Alexander R. A. Anderson . A spatial model of tumor-host interaction: Application of chemotherapy. Mathematical Biosciences and Engineering, 2009, 6(3): 521-546. doi: 10.3934/mbe.2009.6.521
    [10] Rebeccah E. Marsh, Jack A. Tuszyński, Michael Sawyer, Kenneth J. E. Vos . A model of competing saturable kinetic processes with application to the pharmacokinetics of the anticancer drug paclitaxel. Mathematical Biosciences and Engineering, 2011, 8(2): 325-354. doi: 10.3934/mbe.2011.8.325
  • Black rice bran (BRB) is well-known for its high antioxidant activity and its pivotal role in preventing colon cancer. The present study aims to investigate the effects of BRB administration on BALB/C mice induced with azoxymethane (AOM) and dextran sodium sulphate (DSS). The 24 mice were divided into three groups: the group induced by colon cancer (C+), the group induced by cancer and given the BRB diet (C+BRB), and the normal group (C−). Both the C− and C+ groups were given a standard AIN-93 M diet containing cellulose fiber. After 16 weeks, the mice were anesthetized, and the colonic tissue was identified for nodule distribution, histopathological observation, and mRNA expression analysis of proliferating cell nuclear antigen (PCNA), caspase-3, and caspase-8 genes using qRT PCR technique. Nodule distribution in the C+BRB group showed a significant decrease compared to the C+ group, with 1.65 ± 0.71 nodule/cm2 vs. 5.73 ± 2.93 nodule/cm2. Then, the colon weight was significantly decreased in the C+BRB group, at 0.19 ± 0.04 (g) compared to 0.25 ± 0.03 (g) in the C+ group. Also, the BRB diet in the C+BRB group significantly decreased PCNA mRNA expression compared to the C+ group, with values of 0.58 ± 0.09-fold change vs. 5.22 ± 0.80-fold change. Conversely, increased the mRNA expression of caspase-3 (0.91 ± 0.20-fold change vs. 0.36 ± 0.15-fold change) and caspase-8 (0.51 ± 0.18-fold change vs. 0.13 ± 0.31-fold change). In conclusion, administration of BRB inhibited the rate of cancer development by suppressing cancer cell proliferation and inducing apoptosis.



    Therapeutic monoclonal antibody (mAb) therapy has exhibited remarkable success over the past decade: not only did mAb-based therapeutics account for four of the top five best-selling drugs in 2018, but the Nobel Prize in Physiology and Medicine was awarded that same year for the discovery of the mAb-based cancer immunotherapy agent ipilimumab [1,2]. Over 500 mAbs are currently being sponsored for clinical trial studies across several disease types, with cancer as the most common indication: approximately 70% of mAbs in Phase I trials are for oncological indications [3].

    Antibodies are vastly larger than small molecule drugs with a molecular weight of approximately 145 kDa, and in contrast to small molecules, are generally less toxic due to their exceptional specificity to target antigens [4]. Due to their size, antibodies are generally restricted to the extracellular space (i.e., vascular and interstitial) and thus, have low volumes of distribution in humans (5–6 L) [5]. Not only are they much larger, they usually remain in circulation for long periods of time (t1/2~3 weeks), mainly due to the presence of the neonatal Fc receptor (FcRn). FcRn recycling of IgG antibodies, which includes all of the therapeutic antibodies on the market, prevents their degradation in the intracellular lysosomal compartment and can facilitate antibody translocation through cells [5]. In terms of administration, therapeutic antibodies are most commonly delivered either intravenously (IV) or subcutaneously. Limitations of IV dosing, such as discomfort to the patient and an increased number of clinical visits, render the subcutaneous (SC) route an increasingly popular option for mAb delivery While the specific mechanisms governing mAb absorption via SC injection are not completely understood, subcutaneous delivery of mAbs and eventual distribution to the systemic circulation is proposed to be largely dependent on the uptake (i.e., drainage) via the local lymphatic vasculature (i.e., lymphatic capillaries/vessels, lymph nodes) surrounding the injection site [6,7]. Therefore, the lymphatic system plays an indispensable role in the absorption, distribution, target accumulation, and elimination properties of mAbs.

    The lymphatic system is a highly complex, multifaceted organ system whose overall main functions are to drain local tissue waste fluid (i.e., lymph) back into the circulation as well as provide host immune defense against foreign pathogens and/or abnormal (i.e., cancerous) cells [8]. It is comprised of an intricate network of lymphoid organs and tissues, including the lymphatic vasculature and lymph nodes interspersed throughout that drain fluid from surrounding tissues. When antibodies are delivered subcutaneously, they first become exposed to the interstitial space of the hypodermis, which is highly vascularized by both blood and lymphatic capillaries [9]. While FcRn-dependent mAb transcytosis through the blood capillary endothelial cell into circulation is possible, drainage via lymphatic capillaries is the primary mechanism of mAb transport out of the injection site interstitial space [9]. Generally speaking, once drained by the lymphatic capillaries, mAbs will travel through various lymph node networks (joined by larger lymphatic vessels) before converging to a central duct, and ultimately, into the systemic circulation (Figure 1). Unlike the circulatory system, the lymphatic system has no central pump; therefore, it must rely on interstitial fluid buildup (to initiate the drainage cascade) along with local muscle contraction to effectively circulate lymph back into the systemic circulation [10]. As a result, lymph flow is dramatically slower than blood flow; on average, 500x slower (0.2% of blood flow) [11], which contributes to the long peaking time (Tmax = 2–14 days) in the absorption of mAbs after subcutaneous dosing [12].

    Figure 1.  Generalized lymphatic drainage map and transit of therapeutic antibodies after either intravenous or subcutaneous administration. (a) Lymph node clusters (green dots), scattered throughout the body, drain lymph from surrounding tissues and into one of two subclavian veins (yellow and red stars). Lymphatic drainage is asymmetrical; the right lymph duct (yellow star) drains the right side of the head and neck, right arm/upper right core, whereas the thoracic duct (red star) drains the entire remainder of the body. (b) The flow after subcutaneous administration to the upper left arm is dependent on the (c) local system lymphatics to reach the systemic circulation.

    The various lymph node networks scattered throughout the body serve as immune checkpoints to survey local tissue for foreign pathogens or cancerous cells in draining lymph. Naïve lymphocytes (T- and B- cells), which are of bone marrow origin, travel through the bloodstream and reside in lymph nodes. Once presented with antigen from afferent lymph, these lymphocytes can induce an immune response [8]. Therapeutic antibodies are also subject to this same route for distribution; a generalized schematic of the mAb lymphatic drainage pathway is shown in Figure 1 below. The co-localizations of antibodies and lymphocytes are critical for antibodies to trigger the effector functions in the tumor-draining lymph nodes (TDLNs) for tumor eradication, which may be interrupted by surgery and the removal of TDLNs, dampening the activation of the immune system by antibodies and anti-tumor efficacy.

    As shown in Figure 1, lymphatic drainage in the human body is asymmetrical; a disproportionate amount of total lymph (~75%) is drained into the left thoracic duct [13]. However, a study by Xu et al in 2010 showed that the SC injection site did not significantly impact the absorption of golimumab in healthy subjects, suggesting it may not be necessary to stratify by anatomical region from a modeling perspective [14]. Furthermore, in addition to providing local tissue fluid and immune homeostasis, the lymph node clusters also vary greatly in the number, size, and even in their local immune signature (i.e., "immunostats") [15]. A table mapping the lymph node clusters draining various organs is summarized in Table 1; the lymph network structure and cluster density are critical to understanding the distribution of immunotherapies as well as their immune activation efficacy.

    Table 1.  Lymphatic drainage pathway by organ, along with literature reported estimates of number and size (if available). LN = lymph nodes.
    Organ Major LN drainage network Number Size
    Brain Cervical 300 [16] 0.48 mL (volume) [16]
    Lung Mediastinal 225 [17] ~12.6 x 8.3 mm [17]
    Breast Axillary 20–50 [18] < 10 mm (diameter) [18]
    Skin Variable --- ---
    Kidney Aortic 30–50 [19] 5–8 mm (diameter) [19]
    Liver Hepatic 30–50 [20] 3–5 mm (diameter) [20]
    GI Mesenteric 100–150 [21] 3–5 mm (diameter) [20]

     | Show Table
    DownLoad: CSV

    While the lymphatic vasculature and nodes play a vital role in tissue fluid/immune homeostasis, it is also highly exploited in cancers [22]. At its core, cancer is a disease of genetic instability; tumors acquire the ability to evade host immune response via immunosuppressive mechanisms in the tumor microenvironment. Moreover, a significant contributor to cancer lethality is the ability to invade (i.e., metastasize) to distal organ and tissue sites. Metastasis through the lymphatic system (as opposed to systemic blood circulation) is observed in four out of every five melanomas and carcinomas [22]. Not only does the lymphatic vasculature provide a physical basis for metastatic spread, but the lymph nodes that directly drain primary tumors (i.e., tumor-draining lymph nodes, TDLNs) become a highly opportunistic, metastatic site. Tumors generally agammaegate as highly compact, hypoxic clusters of cells exerting high interstitial pressure; this high pressure is proposed to be a contributing mechanism to the dissemination of cancer cells to the TDLN [23]. Furthermore, both the primary tumor and TDLN have been shown to induce lymphangiogenesis (i.e., formation of new lymphatic vessels) to promote both local and distal metastatic spread; alarmingly, the presence of tumor cells in the TDLN was not a requirement for TDLN-mediated lymphangiogenesis. This is instead mediated by vascular endothelial growth factor A (VEGF-A) secreted by the primary tumor, which once drained via the lymphatic capillaries and is in the TDLN, induces lymphangiogenesis distally. However, once metastasized to the TDLN, cancer cells are able to promote further lymphatic vessel formation into an increasingly greater number of distal lymph nodes and tissues [24]. Moreover, rather than elicit an anti-tumor immune response, quite paradoxically, TDLNs promote both local and systemic tolerance to tumor-derived antigens. Increasing evidence has shown that the functional activity of antigen-presenting cells (APCs) in the TDLN, along with enhanced tumor-protective host regulatory T-cell (Treg) activation, promotes a tolerizing rather than immunogenic effect [25]. In light of these findings, mAb immunotherapies aiming to reverse tumor-mediated immunosuppression such as the anti-PD-1 mAb pembrolizumab have been developed and have shown great promise across a broad spectrum of cancers (Figure 2). In fact, both the primary tumor and the TDLN have been heavily implicated as a focal point for mAb immunotherapy via mechanisms that will be explained in a later section.

    Figure 2.  FDA-approved anti-PD-1 pembrolizumab cancer indications (2019) [26].

    A second fundamental principle of cancer describes its disease progression as a microevolutionary process; cancer cells are continually outcompeting their neighbors for resources. Those acquiring mutations conferring a selective advantage will predominate over other sub-clones (i.e., survival of the fittest). Over time, this will create a heterogeneous lesion profile even within the same patient; as tumors metastasize to distal organs, a complex interplay between the tumor microenvironment and the local immune profile (i.e., organ-specific immunostats) further exacerbates this heterogeneity. This is reflected in the variability of response of individual lesions to pembrolizumab [27]. In this study, Osorio et al analyzed individual lesions from patients with metastatic non-small cell lung cancer (NSCLC) and mismatch repair deficiency (MMRD) carcinoma treated with PD-1 monotherapy, including pembrolizumab. CT scans of patients treated with PD-1 therapy, an immunosuppressive protein expressed in T-cells, were analyzed at the lesion level via the Response Evaluation Criteria in Solid Tumors (RECIST) v1.1 criteria to evaluate individual lesion pattern of progression. More specifically, the RECIST v1.1 criteria assesses changes in overall patient tumor burden upon treatment using specific response metrics. Predetermined guidelines and tumor shrinkage cut-offs are used to denote whether a tumor has achieved no, partial, or complete response. These metrics may then be correlated to disease progression endpoints such as progression-free survival (PFS) and overall survival (OS) [28]. Results revealed a site-specific difference in response to treatment; lymph node lesions tended to show the greatest response, while bone and liver metastases were among the least responsive [27]. These results strongly suggest that immunologic response is not homogenous and there may likely be organ-specific factors at play influencing the level of response. Furthermore, these results also argue against the pooling of individual patient lesions to summarize overall disease progression, the current procedure for RECIST v1.1 analyses [28]. Given the established role of TDLNs in influencing both local and systemic tolerance to tumor antigens, the unique TDLN networks draining their anatomically distinct primary tumors become a highly attractive target for further exploration.

    Arguably, the lymphatic vasculature in the context of cancer progression is understudied, and its role in influencing mAb PK/PD remains largely unknown. As previously mentioned, one of the most common routes of mAb administration (i.e., subcutaneous injections) is partly dependent on transit via local lymphatic vessels and lymph nodes for entry into the systemic circulation. Thus, the findings mentioned above substantiate the need to develop a platform that systematically analyzes individual, inherently heterogeneous tumors, and those metastasized to TDLNs via lymphatic vessels, while accounting for mAb disposition and binding at each of these lesion sites. PK/PD modeling provides a quantitative exploration of antibody disposition at the lesion level. It offers the potential to model and simulate several clinically relevant scenarios that would otherwise be extremely costly to evaluate experimentally [29]. The following section will discuss modeling approaches, such as physiologically-based pharmacokinetic (PBPK) and target-mediated drug disposition (TMDD) modeling, to assess therapeutic antibody disposition and binding to both primary tumors and TDLNs. A summary of preclinical PK/PD mAb lymphatic distribution studies and the extrapolation of preclinical mAb PBPK approaches to clinical studies will also be provided.

    PBPK modeling provides a powerful, quantitative opportunity to investigate organ and tissue-specific processes that contribute to mAb disposition and binding in both primary tumors and TDLNs. At its core, PBPK modeling aims to incorporate physiologically-relevant parameters, such as the anatomical volumes of lymph nodes and the tumor-specific lymph flow, to characterize a relatively complete PK footprint of the relevant system [30]. This is in contrast to conventional mammillary models where the lumped compartments used to describe the systemic drug PK, rather than drug disposition in target tissues, do not carry much physiological relevance (Figure 3).

    Figure 3.  PBPK modeling allows for evaluation of mAb disposition and binding in tumor-specific organs and TDLNs. (a) A series of mass-balance differential equations describing the time course of antibody concentrations in both primary tumors (PT) and TDLNs are derived and implemented into the model; every organ-specific primary tumor and TDLN serves as its own physiological compartment. Typically, these compartments are connected via blood flow (black arrows), and lymph flow (green structures) is secondary; however, when considering antibody disposition, the lymphatic vasculature must be emphasized. Fr denotes the flow fraction of blood/lymph to a given tissue (i.e., tumors), which varies across tissues, tumor types/sizes, and tumor vascular structure. (b) Tissue- and organ-specific schematic and model parameters in primary tumors and TDLNs: Organ blood flow (Qorgan), vascular reflection coefficient (σv), lymphatic reflection coefficient (σL), afferent lymph flow (Laff), efferent lymph flow (Leff), organ lymph flow (Lorgan), mAb clearance from central compartment (CLp), interstitial space (IS). Green circular structures denote T-cells in both tumor and TDLN IS compartments interacting with draining antibody. (c) Example simulation illustrating antibody disposition (SUV, standard uptake value) in plasma (blue), primary tumor (red), tumor-positive TDLN (solid green), and tumor-negative TDLN (dashed green). The purple and green arrows denote the concentration gradients between PT-TDLN and the effect of antibody binding in the TDLN, respectively (i.e., metastasis to TDLN). For further exploration on PBPK models, please consult the following references [32,34,35].

    As shown in Figure 3, typical PBPK models emphasize blood flow whereas lymph flow is secondary. While both flows have been accounted for in previously published models [31,32,33], in the context of antibody disposition/binding to both primary tumors and lymph nodes, anatomically distinct tumors are ignored and instead lumped into a "tumor" compartment; lymph nodes are also lumped together and antibody binding to antigens in the lymph nodes has not yet been investigated. From a mass balance perspective, previous PBPK models provided a useful approach to capture whole-body antibody disposition; however, tumor and lymph node specificity is sacrificed. Thus, a novel approach to emphasize mAb binding and disposition to organ-specific primary tumors as well as to their draining lymph nodes must be considered. Our group is currently investigating these questions by utilizing existing positron electron tomography (PET) imaging clinical data (Section 2.4) to build a tumor-specific/TDLN PBPK model. This will be critical in understanding the concentration gradient between these two tissues, while also assessing if antibody can achieve sufficient exposure for full immune activation not only in the primary tumors but also in the TDLNs. Therefore, the following section will compartmentalize the primary tumor-TDLN axis and will describe the differential equations necessary to quantify antibody disposition at these tissue sites.

    The model structure described in Figure 4 summarizes the components involved in modeling mAb disposition and binding across both primary tumors and TDLNs. As described in the schematic, instead of using a full-body PBPK model, a simplified minimal PBPK (mPBPK) model can be applied to emphasize antibody PK in both of these tissues. In fact, mPBPK models describing antibody disposition in similar circuits have previously been described [36,37,38,39].

    Figure 4.  mPBPK model structure of antibody disposition in plasma, primary tumors and TDLNs. (a) Sequential mAb flow as it extravastates from the vascular space into the tumor interstitial space (IS), subsequently into the TDLN IS via afferent lymph flow, and then exits via efferent lymph flow. Membrane-bound antigen present at any site (i.e., tumors) will retain bound antibody, and only free (i.e., unbound) antibody will remain in circulation. (b) Antibody-antigen binding in the primary tumor IS space is described via a quasi-equilibrium target-mediated drug disposition (QE-TMDD) model. Free antibody leaving the tumor IS compartment will then travel via the lymphatic vessels to the TDLN where a second QE-TMDD model can be applied in the TDLN, assuming antigen (i.e., tumor) is present. Efferent flow will then deliver remaining antibody back into the plasma. Symbols are defined below. (Figure prepared in BioRender).

    The differential equations describing the model are as follows (assuming IV injection):

    Vp(dCAb,pdt)=Input(1σv)frLorganCAb,pCLpCAb,p+LeffCf2 (1)
    VisfPT(dCAb,PTdt)=(1σv)frLorganCAb,p(1σL)LaffCf1kintARPTVisfPT (2)
    dRPTdt=ksynkdeg(RT,PTARPT)kintARPT (3)
    Visf,TDLN(dCAb,TDLNdt)=(1σL)LaffCf1LeffCf2kintARTDLNVisfTDLN (4)
    dRTDLNdt=ksynkdeg(RT,TDLNARTDLN)kintARTDLN (5)

    CAb, p, CAb, PT, and CAb, TDLN denote antibody concentrations in the central (i.e., plasma) compartment, primary tumor (PT), and TDLN, respectively. RPT and RTDLN refer to total antigen concentration in primary tumor and in TDLN. Initial conditions for each differential equation are as indicated: IC(Cab, p) = Dose (assuming IV route of administration); IC(CAb, PT) = 0; IC(RPT) = ksyn/kdeg; IC(CAb, TDLN) = 0; IC(RTDLN) = ksyn/kdeg. Vp, Visf, PT, and Visf, TDLN refer to plasma volume, primary tumor interstitial fluid volume, and TDLN interstitial fluid volume. Lorgan, Laff, and Leff are the organ lymph flows quantifying antibody extravasation from the vasculature into the primary tumor interstitial space, as well as TDLN afferent and efferent lymph flows, respectively. CLp is the clearance of antibody from the central compartment, and σv and σL are the vascular and lymphatic reflection coefficients. Fr denotes the flow fraction of blood/lymph to a given tissue (i.e., tumors), which varies across tissues, tumor types/sizes, and tumor vascular structure. ARPT and ARTDLN are bound antibody-antigen concentrations at these sites, defined by the following equations:

    ARPT=RT,PTCf1Kd+Cf1 (6)
    ARTDLN=RT,TDLNCf2Kd+Cf2 (7)

    Therapeutic antibodies have been shown to undergo a phenomenon known as target-mediated drug disposition (TMDD), whereby the exceptionally high affinity for antibody to an antigen, along with limited target capacity, influences its pharmacokinetic properties [40]. For membrane-bound targets, this is reflected by the internalization rate (kint) of the antibody-antigen complex, followed by degradation. The TMDD structure is depicted in Figure 3 within both the primary tumor and TDLN compartments, and due to model complexity, can be simplified according to certain assumptions. A quasi-equilibrium TMDD (QE-TMDD) approximation is used here, wherein many cases, equilibrium binding (Kd) is on a much faster time-scale than the other parameters describing target synthesis (ksyn), degradation (kdeg) and internalization (kint). TMDD models are constructed so that only total antibody and total receptor differential equations are required; free antibody concentrations within either the primary tumor (Cf1) or TDLN (Cf2) interstitial space are expressed as a quadratic function dependent on total antibody and total target concentrations as well as the Kd, and are defined as follows:

    Cf1=0.5{(CAb,PTRT,PTKd)+(CAb,PTRT,PTKd)2+4CAb,PTKd} (8)
    Cf2=0.5{(CAb,TDLNRT,TDLNKd)+(CAb,TDLNRT,TDLNKd)2+4CAb,TDLNKd} (9)

    Sources for parameter values to implement in the model can be consulted in the following references [34,41]; however, it must be emphasized that the volumes and flows are largely organ- and tumor-specific. A study by Brown et al in 2019 assessed tumor-specific perfusion which can be used to calculate tumor-specific processes, such as the fr (fraction of tumor-specific lymph flow) [42]. Once built, the model and its parameters may be evaluated by utilizing PET imaging data that details antibody binding to both PT and TDLNs, as discussed in Section 2.4. There are several critical assumptions and notable observations in this model structure; as previously mentioned, due to their large size, antibodies are mainly restricted to the vascular and target site interstitial spaces and should be reflected accordingly in the differential equations. Moreover, this model structure quantifies only one primary tumor to its sentinel TDLN or the first LN that directly drains the primary tumor. It is assumed that all of the free antibody leaving the TDLN interstitial space flows directly and entirely into the central compartment. Additionally, processes such as lymphangiogenesis and FcRn-mediated salvage within the lymphatic capillaries are not accounted for in this model. Previous modeling analyses and bioavailability studies indicate that roughly 2/3 of antibody is recycled (ranging from ~52–100%) through this pathway [43]; thus, FcRn recycling must be implemented into the model if it is critical for the research question or it may be empirically included. Furthermore, as previously stated, these differential equations capture these concentration gradients after IV injection. If presented with SC data, an additional absorption compartment has to be considered; previous models have captured this process as referenced [34,43] Nonetheless, the model schematic and differential equations illustrated here provide a solid starting point for quantitative investigation of antibody disposition in primary tumors and TDLNs. The following sections will summarize the landscape of previous preclinical PK/PD antibody lymphatic studies as well as provide insights on immuno-PET imaging methods as a non-invasive way to study antibody pharmacokinetics in primary tumors and TDLNs in the clinical setting.

    While there have been several published PBPK models assessing mAb disposition in animal tumors, as previously mentioned, none to-date have analyzed lesion-level organ-specific tumor heterogeneity. Instead, these models have grouped all tumors into a "tumor" compartment, and thus, lose the ability to stratify by organ lesion type [32,33,35]. The same can be said about LNs and TDLNs; these lymphatic structures are usually secondary to blood flow, and while accounted for in the model, are seldom systemically assessed. While the scope of this review emphasizes antibody disposition in primary tumors and TDLNs, it is important to get an understanding of all of the processes that affect mAb PK at these tissue sites. As previously mentioned, subcutaneous injections are an increasingly popular route of administration for mAbs, and as such, factors such as the injection site interstitial space and lymphatic transit, among others, must be evaluated. Charman et al. investigated causes of degradation of human growth hormone (hGH), a different class of biologics, via the lymphatics in sheep; they reported that injection site interstitial space only minimally contributed to hGH degradation. However, it was shown that overall degradation via the lymphatics contributed significantly, as only 62% of hGH was recovered in peripheral lymph from the efferent duct of the popliteal lymph nodes [44]. A subsequent analysis then demonstrated that degradation of therapeutic proteins in freshly collected lymph was also minimal in rats and dogs; however, proteolysis in the LNs was shown to contribute to decreased mAb bioavailability [45]. It is not surprising given that resident protease-secreting immune cells in the LNs can catabolize mAbs. Furthermore, the role of FcRn salvage in the lymphatic vessels is vital in maintaining intact antibody as FcRn-knockout mice are devoid of all IgG within four days; over 50% of IgG remained in wild-type mice after several days [46], consistent with the half-life of IgG antibodies. While the exact processes governing antibody lymphatic absorption remain unknown, these findings highlight various components of the lymphatic vasculature that should be considered when modeling antibody disposition in primary tumors and TDLNs.

    There are no direct clinical studies that have assessed mAb disposition in primary tumors and TDLNs. A study published by Varkhede et al. analyzed published data across several clinical mAb PK studies (in healthy volunteers) to build a minimal PBPK model to quantify mAb absorption after subcutaneous administration. They attributed a loss of antibody bioavailability to interstitial proteolysis, likely in the sentinel LNs, supporting previous findings [34]. Interestingly, they attributed antibody-specific factors as a contributing factor to the reduction in mAb absorption; physiochemical factors such as the mAb isoelectric point was shown to correlate with interstitial LN proteolysis, suggesting that factors inherent to the mAb need to be considered as well. However, a couple of limitations to this study were that actual LN concentration data was unavailable and thus was not directly calibrated or assessed. Furthermore, this model was not optimized in a cancer setting.

    Several studies have used immuno-PET imaging as a non-invasive way to quantify antibody uptake and binding across various tissues, including primary tumors and TDLNs [47,48,49]. Immuno-PET imaging takes advantage of a radiolabeled tracer (89Zr-, 64Cu-, etc.) on the mAb and is able to quantify whole-body antibody biodistribution in vivo and through time. Typically, the standard uptake value (SUV) is used as the quantitative indicator of uptake. It is simply a ratio of the target organ/tissue radioactivity divided by whole-body radioactivity, normalized by body weight. SUV values also take into consideration the half-life of the radioisotope and are thus adjusted for time-decay. A higher SUV value denotes higher accumulation and/or binding in target tissues. Therefore, it provides an opportunity to build a PBPK framework for the investigation of antibody disposition and binding in primary tumors and TDLNs. However, some limitations with this approach include sparse sampling, limited TDLN data, and a lack of specific information on anatomical TDLNs collected. Furthermore, SUV data can be variable even with the same lesion and between scans [50]; nonetheless, these empirical values still provide a valuable resource for model calibration and analysis in the clinical setting.

    Cancer surgery undoubtedly remains the gold standard for treating tumors if "resectable for cure, " and is much more tolerated than either chemotherapy or antibody therapy. However, as previously mentioned, metastases, and especially micrometastases that go unnoticed, are the overwhelming cause of cancer-related deaths. In terms of cancer mAb therapies, especially in the immunotherapy context (i.e., checkpoint blockade mAbs such as pembrolizumab), the currently accepted approach is to administer these therapies in the late-stage, adjuvant (i.e., post-primary tumor resection) setting. However, an emerging body of evidence is supporting immunotherapy administration in much earlier settings, before surgery (i.e., neoadjuvant immunotherapy); these studies are being supported across several different cancer types. A study in early-stage colon cancer revealed that 95% of mismatch repair (MMR)-deficient (dMMR) tumors treated with neoadjuvant immunotherapy exhibited major pathological response (≤10% residual viable tumor), with 60% of all lesions showing complete response (100% reduction) [51]. A study in recurrent glioblastoma showed that neoadjuvant pembrolizumab extended median progression-free survival (PFS) from 72.5 days (adjuvant group) to 99.5 days (neoadjuvant group); enhanced clonal expansion of T-cells and decreased PD-1 expression on peripheral T-cells were also reported in the neoadjuvant group [52]. Furthermore, in locally advanced melanomas, a potent and rapid anti-tumor response was observed in patients receiving neoadjuvant immunotherapy, with 8 of 27 patients experiencing a major or complete pathological response. All eight of those patients remained disease after only a single dose of pembrolizumab and were associated with T-cell reinvigoration seen as early as one-week post-treatment [12]. Collectively, these findings strongly suggest that antibody administration in the neoadjuvant setting is modulating a fundamental principle of possibly all cancers.

    The primary tumor-TDLN signaling axis has been highlighted as a focal point for neoadjuvant immunotherapy success [53]; in the neoadjuvant setting, while the primary tumor is still intact, this allows for more efficient antibody-mediated pharmacodynamic T-cell priming at both the primary tumor (PT) site as well as in the TDLN due to the higher levels of endogenous tumor antigen present. Moreover, an intact PT-TDLN signaling axis will ensure an adequate antibody concentration gradient at both sites (i.e., sufficient antibody delivery at both the PT and TDLN). This enhanced systemic immunity is then better able to scavenge and eliminate micrometastases in tissues (via the intact lymphatic vasculature and the blood vasculature) that would otherwise be the source of post-surgical relapse. In the case of immunotherapy administration such as pembrolizumab, which targets the PD-1 receptor in T-cells, this will allow for an even greater reinvigoration of T-cells in the tumor microenvironment, reversing tolerance induction caused by tumor antigens such as the PD-1 receptor ligand PD-L1 [53]. However, careful attention must be paid to the timing of neoadjuvant immunotherapy and allowing adequate time for systemic immunity, as preclinical evidence shows an optimal time window before tumor burden supersedes therapeutic benefit [54]. In cases where antibody is administered post-surgery, insufficient mAb delivery (i.e., lower/inefficient biodistribution gradient) to the TDLN, due to a disruption in the local lymphatic vasculature, may have deleterious consequences, such as micrometastatic tumor relapse/resistance. Thus, maintaining adequate mAb biodistribution gradients to the TDLN are critical in driving both antibody PK and PD response in the clinical setting.

    Sentinel lymph node (SLN) biopsies is a technique that arose in the early 1990s that is used to assess the metastatic potential of the primary tumor and involves excising surrounding proximal lymph nodes that are believed to be the first set of nodes draining the tumor. However, the controversy over the overall risk-benefit utility of this procedure has been raised throughout the years as lymphedema (swelling of the excised area due to compromised lymphatic drainage) is presented in nearly 30% of SLN patients [55]. Not only does this impede the patient's quality of life, but it also provides a further pharmacokinetic barrier of efficient antibody delivery to more distal lymph nodes surrounding the region. Furthermore, as previously mentioned in the neoadjuvant immunotherapy case, this insufficient delivery will likely result in suboptimal mAb-antigen binding found in micrometastases in the TDLN, leading to an attenuated anti-tumor response for other antibody types, such as trastuzumab.

    The Varkhede et al. subcutaneous lymphatic mAb absorption study that was previously discussed noted a co-formulation effect on affecting antibody PK in the interstitial space; hyaluronidase is an enzyme that is commonly co-administered with therapeutic antibodies as it helps ease discomfort with subcutaneous injections. It allows for higher mAb injection volumes. It is also known to disrupt the hyaluronic acid-dependent subcutaneous extracellular matrix present in the interstitial space, allowing for more efficient, rapid delivery [34]. Their original model only estimated the interstitial LN proteolysis (i.e., LN-mediated mAb clearance parameter, CLLN) and, in doing so, overpredicted the Tmax for those mAbs co-formulated with hyaluronidase. Once the subcutaneous injection site lymph flow rate was allowed to be simultaneously estimated with the CLLN parameter, the model fits were much improved; final model estimates for the injection site lymph flow rate were markedly higher and ranged anywhere from 20- to 330- fold than the original value [34], consistent with the known pharmacology of hyaluronidase. Therefore, a detailed formulation description of the therapeutic antibody should be evaluated when modeling disposition in the lymphatic vasculature and LNs.

    Inflammation is the body's natural response to combat pathogens and is essential for tissue immune homeostasis. However, tumor-associated inflammation, which has been proposed as the seventh hallmark of cancer [56], can make profound changes within the tumor microenvironment and in the surrounding lymphatic vasculature and TDLNs. Both pro- and anti-inflammatory cytokine release have been shown to heighten interstitial flow in the lymphatic capillaries and vessels. Coupled with lymphatic vessel enlargement, it can drastically alter the PK of antibody delivery to both primary tumors and TDLNs [57]. Presumably, this would result in a diluting effect that could give rise to suboptimal mAb concentrations in the TDLN.

    The tumor extracellular matrix (ECM) is markedly different from healthy tissue; it is composed of a complex network of myofibroblasts, collagen and other ECM components such as hyaluronan and proteoglycans (PGs), among other components. ECM content, even within different anatomical tumors, can drastically differ [58] and can have dramatic effects on antibody disposition at these sites. These findings corroborate well with a couple of aforementioned studies; lesion-level heterogeneous responses to pembrolizumab have been previously explained and can, in part, be attributed to these anatomical differences. Furthermore, negatively-charged hyaluronan, which is overexpressed in tumors, aids in tumor invasiveness and the epithelial-mesenchymal transition (EMT) [58]. Given that many mAbs usually carry a positive charge, this might enhance electrostatic interaction and thus retain mAbs within the tumor ECM and create an even higher concentration gradient between the primary tumor and TDLN. For these reasons, antibody disposition in organ-specific primary tumors and TDLNs should be emphasized and characterized accordingly.

    In summary, while largely understudied, the TDLN is an invaluable asset to the tumorigenic and metastatic cascade; it not only acts as a physical means of cancer cell dissemination but is also a dynamic regulator of both local and systemic tolerance to tumor antigens. Cancer cells are able to exploit the TDLNs in such a way that ultimately progresses, rather than prohibits, further spread. In the context of mAb therapies that depend on interstitial and lymphatic uptake to both the primary tumor and TDLN, several clinical as well as pathophysiological factors should be considered when quantitatively modeling this axis. While seemingly identical in shape and universal in function, organ-specific immuno-regulation can vastly influence the local immune signature, which can differentially affect antibody responses to anatomically distinct lesions. The work described here has provided a PBPK/TMDD modeling approach to successfully characterize mAb disposition in organ-specific primary tumors as well as in the TDLNs that drain them. Future clinical studies should take advantage of immuno-PET imaging methods and emphasize distinct TDLN networks throughout the body to predict antibody disposition at these sites more reliably.

    This research is supported by NIH grant GM119661.

    The authors declare no conflicts of interest. E.S. is a postdoctoral fellow sponsored by the UNC/IQVIA fellowship.



    [1] World Health Organization (2017) Cancer. Available from: http://www.who.int/mediacentre/.
    [2] American Cancer Society (2017) Cancer Facts and Figures in 2017. American Cancer Society. Inc. Surveillance Research.
    [3] Datsis A, Tsoga A, Langouretos V (2010) The role of functional foods in the prevention of colorectal cancer. Hellenic J Surg 82: 224–232. https://doi.org/10.1007/s13126-010-0035-5 doi: 10.1007/s13126-010-0035-5
    [4] Moongngarm A, Daomukda N, Khumpika S (2012) Chemical compositions, phytochemicals, and antioxidant capacity of rice bran, rice bran layer, and rice germ. Apcbee Proc 2: 73–79. https://doi.org/10.1016/j.apcbee.2012.06.014 doi: 10.1016/j.apcbee.2012.06.014
    [5] Antczak C, Takagi T, Ramirez CN, et al. (2009) Live cell imaging of caspase activation for high content screening. J Biomol Screen 14: 956–969. https://doi.org/10.1177/1087057109343207 doi: 10.1177/1087057109343207
    [6] Shi Y (2002) Mechanisms of caspase activation and inhibition during apoptosis. Mol Cell 9: 459–470. https://doi.org/10.1016/S1097-2765(02)00482-3 doi: 10.1016/S1097-2765(02)00482-3
    [7] Zheng Z, Shetty K (2000) Solid-state bioconversion of phenolics from cranberry pomace and role of Lentinus edodes beta-glucosidase. J Agric Food Chem 48: 895–900. https://doi.org/10.1021/jf990972u doi: 10.1021/jf990972u
    [8] Bravo R, Macdonald-Bravo H (1987) Existence of two populations of cyclin/proliferating cell nuclear antigen during the cell cycle: association with DNA replication sites. J Cell Biol 105: 1549–1554. https://doi.org/10.1083/jcb.105.4.1549 doi: 10.1083/jcb.105.4.1549
    [9] Zulfafamy KE, Ardiansyah, Budijanto S (2018) Antioxidative properties and cytotoxic activity against colon cancer cell WiDr of Rhizopus oryzae and Rhizopus oligosporus-fermented black rice bran extract. Curr Res Nutr Food Sci 6: 23–24. https://doi.org/10.12944/CRNFSJ.6.1.03 doi: 10.12944/CRNFSJ.6.1.03
    [10] Glasauer A, Chandel NS (2014) Targeting antioxidants for cancer therapy. Biochem Pharmacol 92: 90–101. https://doi.org/10.1016/j.bcp.2014.07.017 doi: 10.1016/j.bcp.2014.07.017
    [11] Suzuki R, Kohno H, Sugie S, et al. (2006) Strain differences in the susceptibility to azoxymethane and dextran sodium sulfate-induced colon carcinogenesis in mice. Carcinogenesis 27: 162–169. https://doi.org/10.1093/carcin/bgi205 doi: 10.1093/carcin/bgi205
    [12] Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25: 402–408. https://doi.org/10.1006/meth.2001.1262 doi: 10.1006/meth.2001.1262
    [13] Phutthaphadoong S, Yamada Y, Hirata A, et al. (2010) Chemopreventive effect of fermented brown rice and rice bran (FBRA) on the inflammation-related colorectal carcinogenesis in ApcMin/+ mice. Oncol Rep 23: 53–59. https://doi.org/10.3892/or_00000605 doi: 10.3892/or_00000605
    [14] Trujillo EB, Bergerson ASL, Graf JC, et al. (2005) Cancer in American society for parenteral and enteral nutrition support practice manual. Nutr Care Org 4: 150–170.
    [15] Norazalina S, Norhaizan ME, Hairuszah I, et al. (2010) Anticarcinogenic efficacy of phytic acid extracted from rice bran on azoxymethane-induced colon carcinogenesis in rats. Exp Toxicol Pathol 62: 259–268. https://doi.org/10.1016/j.etp.2009.04.002 doi: 10.1016/j.etp.2009.04.002
    [16] Tomita H, Kuno T, Yamada Y, et al. (2008) Preventive effect of fermented brown rice and rice bran on N-methyl-N'-nitro-N-nitrosoguanidine-induced gastric carcinogenesis in rats. Oncol Rep 19: 11–15. https://doi.org/10.3892/or.19.1.11 doi: 10.3892/or.19.1.11
    [17] Kampa M, Alexaki VI, Notas G, et al (2004) Antiproliferative and apoptotic effects of selective phenolic acids on T47D human breast cancer cells: potential mechanisms of action. Breast Cancer Res 6: R63–R74. https://doi.org/10.1186/bcr752 doi: 10.1186/bcr752
    [18] Rosenberg DW, Giardina C, Tanaka T (2009) Mouse models for the study of colon carcinogenesis. Carcinogen 30: 183–196. https://doi.org/10.1093/carcin/bgn267 doi: 10.1093/carcin/bgn267
    [19] Lupton JR (2004) Microbial degradation products influence colon cancer risk: the butyrate controversy. J Nutr 134: 479–482. https://doi.org/10.1093/jn/134.2.479 doi: 10.1093/jn/134.2.479
    [20] Andoh A, Fujiyama Y (2004) Anti-inflammatory roles of dietary fiber and short-chain fatty acids as regards inflammatory bowel diseases. Agro Food Ind Hi-Tech 15: 42–43.
    [21] Katyama M, Yoshimi N, Yamda Y, et al. (2002) Preventive effect of fermented black rice and rice bran against colon carcinogenesis in male F344 rats. J Oncol 9: 7–22. https://doi.org/10.3892/or.9.4.817 doi: 10.3892/or.9.4.817
    [22] Yi W, Fischer J, Krewer G, et al. (2005) Phenolic compounds from blueberries can inhibit colon cancer cell proliferation and induce apoptosis. J Agric Food Chem 53: 7320–7329. https://doi.org/10.1021/jf051333o doi: 10.1021/jf051333o
    [23] Yin MC, Lin CC, Wu HC, et al. (2009) Apoptotic effects of protocatechuic acid in human breast, lung, liver, cervix, and prostate cancer cells: potential mechanisms of action. J Agric Food Chem 57: 6468–6473. https://doi.org/10.1021/jf9004466 doi: 10.1021/jf9004466
    [24] Muntana N, Prasong S (2010) Study on total phenolic contents and their antioxidant activities of Thai white, red and black rice bran extracts. Pak J Biol Sci 13: 170–174. https://doi.org/10.3923/pjbs.2010.170.174 doi: 10.3923/pjbs.2010.170.174
    [25] Kim HY, Kim JH, Yang SB, et al. (2007) A polysaccharide extracted from rice bran fermented with Lentinus edodes enhances natural killer cell activity and exhibits anticancer effects. J Med Food 10: 25–31. https://doi.org/10.1089/jmf.2006.116 doi: 10.1089/jmf.2006.116
    [26] Zeng H, Lazarova DL, Bordonaro M (2014) Mechanisms linking dietary fiber, gut microbiota and colon cancer prevention. World J Gastrointest Oncol 6: 41–51. https://doi.org/10.4251/wjgo.v6.i2.41 doi: 10.4251/wjgo.v6.i2.41
    [27] Kuno T, Nagano A, Mori Y et al. (2016) Preventive effect of fermented brown rice and rice bran against Prostate Carcinogenesis in TRAP Rats. Nutrient 87: 421. https://doi.org/10.3390/nu8070421 doi: 10.3390/nu8070421
    [28] Gschwantler M, Kriwanek S, Langner E, et al. (2002) High-grade dysplasia and invasive carcinoma in colorectal adenomas: a multivariate analysis of the impact of adenoma and patient characteristics. Eur J Gastroenterol Hepatol 14: 183–188. https://doi.org/10.1097/00042737-200202000-00013 doi: 10.1097/00042737-200202000-00013
    [29] Rosa LDS, Silva NJA, Soares NCP, et al. (2016) Anticancer properties of phenolic acids in colon cancer—A review. J Nutr Food Sci 6: 1000468.
    [30] Ruemmele FM, Schwartz S, Seidman EG, et al. (2003) Butyrate induced Caco-2 cell apoptosis is mediated via the mitochondrial pathway. Gut 52: 94–100. https://doi.org/10.1136/gut.52.1.94 doi: 10.1136/gut.52.1.94
  • This article has been cited by:

    1. Eric Salgado, Yanguang Cao, A Physiologically Based Pharmacokinetic Framework for Quantifying Antibody Distribution Gradients from Tumors to Tumor-Draining Lymph Nodes, 2022, 11, 2073-4468, 28, 10.3390/antib11020028
    2. Mayu Yamamoto, Taiki Kurino, Reiko Matsuda, Haleigh Sakura Jones, Yoshito Nakamura, Taisei Kanamori, Atushi B. Tsuji, Aya Sugyo, Ryota Tsuda, Yui Matsumoto, Yu Sakurai, Hiroyuki Suzuki, Makoto Sano, Kensuke Osada, Tomoya Uehara, Yukimoto Ishii, Hidetaka Akita, Yasushi Arano, Akihiro Hisaka, Hiroto Hatakeyama, Delivery of aPD-L1 antibody to i.p. tumors via direct penetration by i.p. route: Beyond EPR effect, 2022, 352, 01683659, 328, 10.1016/j.jconrel.2022.10.032
    3. Su-Eon Jin, Jino Kim, Jong-Hyuk Sung, Recent approaches of antibody therapeutics in androgenetic alopecia, 2024, 15, 1663-9812, 10.3389/fphar.2024.1434961
    4. Zahra Khademi, Negar Mottaghi-Dastjerdi, Hamed Morad, Narges Hedayati, Ali Sharafi, 2024, Chapter 22, 978-3-031-64827-4, 453, 10.1007/978-3-031-64828-1_22
    5. Xiaoxu Zhong, Yikai Liu, Arezoo M. Ardekani, A compartment model for subcutaneous injection of monoclonal antibodies, 2024, 650, 03785173, 123687, 10.1016/j.ijpharm.2023.123687
    6. Constantinos Harkos, Triantafyllos Stylianopoulos, Investigating the synergistic effects of immunotherapy and normalization treatment in modulating tumor microenvironment and enhancing treatment efficacy, 2024, 583, 00225193, 111768, 10.1016/j.jtbi.2024.111768
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1695) PDF downloads(96) Cited by(0)

Figures and Tables

Figures(3)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog