Korean J Physiol Pharmacol 2023; 27(2): 187-196
Published online March 1, 2023 https://doi.org/10.4196/kjpp.2023.27.2.187
Copyright © Korean J Physiol Pharmacol.
Juyeon Ko#, Jinhyeong Kim#, Jongyun Myeong, Misun Kwak, and Insuk So*
Department of Physiology, Seoul National University College of Medicine, Seoul 03080, Korea
Correspondence to:Insuk So
#These authors contributed equally to this work.
Author contributions: J. Ko designed the study, performed experiments, generated figures, analyzed data, and wrote the manuscript. J. Kim. generated figures and wrote the manuscript. J.M. performed confocal imaging, and M.K. performed Co-IP experiment. I.S. provided the overall experimental advice and coordinated the study. All authors reviewed the manuscript.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License, which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Transient receptor potential canonical (TRPC) channels are non-selective calcium-permeable cation channels. It is suggested that TRPC4β is regulated by phospholipase C (PLC) signaling and is especially maintained by phosphatidylinositol 4,5-bisphosphate (PIP2). In this study, we present the regulation mechanism of the TRPC4 channel with PIP2 hydrolysis which is mediated by a channel-bound PLCδ1 but not by the GqPCR signaling pathway. Our electrophysiological recordings demonstrate that the Ca2+ via an open TRPC4 channel activates PLCδ1 in the physiological range, and it causes the decrease of current amplitude. The existence of PLCδ1 accelerated PIP2 depletion when the channel was activated by an agonist. Interestingly, PLCδ1 mutants which have lost the ability to regulate PIP2 level failed to reduce the TRPC4 current amplitude. Our results demonstrate that TRPC4 self-regulates its activity by allowing Ca2+ ions into the cell and promoting the PIP2 hydrolyzing activity of PLCδ1.
Keywords: Calcium, Fluorescence resonance energy transfer, Phosphatidylinositols, Phospholipase C delta, Transient receptor potential channels
Transient receptor potential (TRP) channels are nonselective cation channels which are permeable to Ca2+. Among TRP channels, this canonical family is composed of 7 types, and they are subdivided into two groups based on their amino acid sequence and structure; TRPC1, 4, 5 and TRPC3, 6, 7. Interestingly, TRPC4 is functionally analogous to TRPC5, but only TRPC5 shows constitutive activity even without stimulation. However, the mechanism that makes this difference still remains obscure. The activity of TRPC4 and TRPC5 are closely related to Gαq-phospholipase C (PLC) signaling, and many studies also demonstrated the effect of downstream molecules on TRPC channels, which include PLC, phosphatidylinositol 4,5-biphosphate (PIP2) [1-4], diacylglycerol (DAG) [5,6], protein kinase C (PKC)  and Ca2+. We previously emphasized the self-limiting mechanism of Gαq pathway on channels where Gαq protein solely activates TRPC1/4 and TRPC1/5 channels, and PIP2 depletion, PKC or Ca2+ inhibit them [8,9]. These data implicate that TRPCs are strongly activated by Gαq, hence fast and tight regulation is necessary afterwards.
PLC is an enzyme that cleaves PIP2 into inositol triphosphate (IP3) and DAG. There are ten kinds of PLC isozymes that are largely classified into three families; β, γ, δ [10-13]. Among them, primordial PLCδ family has basic structure with the shortest length. The δ subtypes have about 10-fold higher sensitivity to calcium than the other isozyme [14,15]. The precise activation mechanism of PLCδ1 is not known, but it is expected that the Ca2+ concentration mainly controls their activity , and IP3, the product of PIP2 hydrolysis, is involved in its negative regulation .
PIP2 is the most important component in TRPC channel’s regulation as they are concerned with channel maintenance. Other TRP channels are also known to be regulated by PIP2 and the effect of PLCδ subtypes on this action has been proposed. PLCδ3 co-expression facilitated the desensitization of TRPV1 in oocytes , and subsequent PIP2 reductions by PLCδ isoform promoted reduction in TRPM8 currents . Thus, understanding the regulatory mechanism that alters TRPC activity related to the change in PIP2 levels will provide insights on how channels are tightly regulated by PLC signaling, and how PIP2 functions differently on TRPC channel functions.
In the present study, we set out to determine an effect of PLCδ on TRPC channel activity, and their working condition in this regulation. We verified that PLCδ1 directly interacts with TRPC4 but not with TRPC5, and it works as a key component of channel’s negative feedback regulation. PLCδ1 is activated by the calcium
Human embryonic kidney (HEK)293 cells (ATCC) were maintained in Dulbecco’s Modified Eagle’s Medium (HyClone) supplemented with 10% heat-inactivated fetal bovine serum (Gibco), 100 U/ml penicillin, and 100 μg/ml streptomycin (HyClone) at 37˚C in a 5% CO2 humidified incubator. Cells were seeded in glass bottom dish for imaging, 12-well plate for whole-cell patch clamp recordings, and 6-well plate for Western blot. The following day, transfection was carried out by using the FuGENE 6 Transfection Reagent (Promega) for patch clamp and imaging experiments, and Lipofectamine 2000 (Invitrogen) for western blotting according to the manufacturer’s instructions. All experiments were performed 20–30 h after transfection. cDNA construct for mouse TRPC4 (GenBank ID: U50921.1, UniProt ID: Q9QUQ5-2) was kindly donated by Dr. V. Flockerzi and Dr. M. Schaefer, and that for human PLCδ1 (GenBank ID: U09117.1, UniProt ID: P51178-1) and human PLCδ3 (GenBank ID: BC072384.1, UniProt ID: Q8N3E9) by Dr. M. Zhu. The mutants of PLCδ1 were generated using a QuickChange site-directed mutagenesis kit (Agilent Technologies) following the manufacturer’s protocol. The sequence was verified by sequencing.
The recording pipette containing standard intracellular solution; 140 mM CsCl, 10 mM HEPES, 10 mM BAPTA, variable CaCl2, 3 mM Mg-ATP, 0.2 mM Tris-GTP, was balanced to pH 7.3 with CsOH. The free Ca2+ concentration was calculated using CaBuf software (G.Droogmans). External solution was perfused constantly as follows; 135 mM NaCl, 5 mM KCl, 10 mM HEPES, 1 mM MgCl2, 2 mM CaCl2, 10 mM glucose, pH 7.4 with NaOH. Rapamycin was purchased from Sigma-Aldrich, and (-)-Englerin A (EA) was purchased from PhytoLab.
The cells were transferred onto a small chamber on the stage of an inverted microscope (IX70; Olympus), and attached to coverslip in the small chamber for 10 min prior to the recording. Transiently transfected cells were identified by their fluorescence tagging. Recording pipettes were pulled from glass capillaries (Harvard Apparatus) using puller (PC-10; Narishige). Whole-cell currents were recorded using Axopatch 200B amplifier (Molecular Devices) and Digidata 1550B interface (Molecular Devices). Experiments were performed at room temperature (18˚C-22˚C). The recording chamber was continuously perfused at a flow rate of 1–2 ml/min. Glass microelectrodes with 2–2.5 megaohms resistance were used to obtain gigaohm seals. The whole cell configuration was used to measure the TRPC4β channel currents in HEK293 cells. Voltage ramps ranging from +100 to –100 mV over a period of 500 ms were imposed every 10 sec with a holding membrane potential of –60 mV. pCLAMP v.10.2 (OriginLab) was used for data acquisition and the data were analyzed using the OriginPro 8 (OriginLab).
HEK293 cells were cultured in a 35-mm glass bottom dishes for imaging. For the confocal images, we used confocal laser scanning microscopy (LSM 710; Zeiss) equipped 63x oil objective lens. To obtain the FRET images, we used an inverted microscope (IX70; Olympus) equipped with 60x oil objective lens (UPlanSApo; Olympus). Each image was captured on an EMCCD camera (iXon3; Andor) and the light at 440 nm and 500 nm wavelengths was illuminated with LED light source (pE-2; CoolLED) under the control of MetaMorph7.6 software (Molecular Devices). Based on this imaging system, FRET measurements were made by the three-cube FRET method  (excitation, dichroic mirror, filter)
We did FRET efficiency computation following methods in Ko
Here, SCUBE(SPECIMEN) denotes an intensity measurement, where CUBE indicates the filter cube (CFP, YFP, or FRET), and SPECIMEN indicates whether the cell is expressing the donor (D; ECFP), acceptor (A; EYFP), or both (DA). RD1 = SFRET(D)⁄(SCFP(D). RD2 = (SYFP (D)⁄(SCFP (D), and RA1 = SFRET (A)/SYFP (A) are predetermined constants from measurements applied to single cells expressing only ECFP or EYFP-tagged molecules. Although three-cube FRET does not require that ECFP and EYFP fusion constructs preserve the spectral features of the unattached fluorophores, similar ratios and recorded spectra furnished two indications that the spectral features of the fluorophores were largely unperturbed by fusion. Since the FR relies on EYFP emission, EYFP should be attached to the presumed limiting moiety in a given interaction. Subsequent quantitative calculations based on FR relied on a presumed 1:1 interaction stoichiometry. The effective FRET efficiency (EEFF) was determined as:
where E is the intrinsic FRET efficiency when fluorophore-tagged molecules are associated with each other, Ab is the fraction of EYFP-tagged molecules that are assoc0iated with ECFP-tagged molecules, and the bracketed term is the ratio of EYFP and ECFP molar extinction coefficients scaled for the FRET cube excitation filter . We determined this ratio to be 0.094 based on maximal extinction coefficients for ECFP and EYFP and excitation spectra .
We performed Western blotting following the methods in Kwak
All statistical analysis and graph generation were done with OriginPro8 (OriginLab). Results were compared using Student’s t-test. A probability value (p) less than 0.05 was considered statistically significant. Data are presented as means ± SEM; *p < 0.05, **p < 0.01, ***p < 0.001.
TRPC4β is functionally very comparable to TRPC5, but only TRPC5 shows constitutive activity even without the stimulation. Previously, we demonstrated that PIP2 binds to TRPC4 and TRPC5 with different affinities . Thus, we hypothesized that TRPC4 and TRPC5 channels have different mechanisms controlling the PIP2 level. To support this hypothesis, we began studying the interaction of PLCδ enzyme with TRPC4 and TRPC5 channels. To identify the expression patterns, we performed fluorescence imaging experiments in a HEK293 cell line transiently co-transfected with ECFP tagged channel (TRPC4β or TRPC5) and EYFP tagged PLCδ (PLCδ1 or PLCδ3) enzyme. In the cases of expressing EYFP-PLCδ alone, PLCδ1 was expressed in the plasma membrane and was also distributed in the cytoplasm, and PLCδ3 was detected in the plasma membrane and the nuclear fraction (Supplementary Fig. 1). However, their expression pattern changed when they were co-expressed with TRPC4β. As shown in Fig. 1A, EYFP-PLCδ1 showed high fluorescence on the same location where TRPC4β-ECFP showed up as puncta on the membrane. The overlay image and line scan in this location showed an analogous fluorescence pattern (Fig. 1A
We used FRET imaging technique to determine whether colocalization of TRPC4β and PLCδ1 resulted from direct interaction or just close distances between these proteins. We used ECFP (donor)-TRPC ion channel and EYFP (acceptor)-PLCδ enzyme constructs to measure FRET efficiencies between them. As it was with the fluorescence imaging experiments, TRPC4β with PLCδ1 showed high FRET efficiency (Fig. 1B), but TRPC4β with PLCδ3, TRPC5 with PLCδ1 and TRPC5 with PLCδ3 showed low FRET efficiency (Fig. 1D).
As a final approach to identify the binding between TRPC and PLCδ, we expressed Flag tagged TRPC ion channel and EYFP tagged PLCδ enzyme and analyzed them
To determine the function of PLCδ1 on TRPC4β, we measured the changes in channel current. We used rapamycin-inducible system to translocate PLCδ1 to the membrane. HEK293 cells were transfected with CFP-FKBP-PLCδ1 and Lyn-FRB, and 50 nM of rapamycin was used, which is sufficient to translocate the PLCδ1 to the membrane (Fig. 2A). To investigate the changes in the TRPC4β current, we performed whole-cell patch clamp experiments with ramp pulse protocol from +100 to –100 mV every 10 sec in cells expressing TRPC4β, CFP-FKBP-PLCδ1, and Lyn-FRB. First, cells were treated with 50 nM rapamycin to translocate PLCδ1 with an efficient amount, then TRPC4β was stimulated by a channel agonist, EA (100 nM). In the cells not expressing Lyn-FRB, rapamycin did not show any current changes compared to basal current, and EA-evoked TRPC4β currents showed the characteristic double-rectifying current-voltage relationship (Fig. 2B). The presence of Lyn-FRB showed the characteristics of TRPC4β, however, the currents with EA stimulation were significantly decreased (Fig. 2C, D). Thus, these results suggest that PLCδ1 functions as a negative regulator of TRPC4β, reducing the current amplitudes.
We also investigated the Ca2+-dependent effect of PLCδ1 on TRPC4β channel currents. TRPC channels are potentiated by intracellular calcium , and especially, Ca2+ is the main regulator of PLCδ1 activity . Thus, we performed whole-cell patch clamp experiments by using different internal Ca2+ buffering solutions. We used 10 mM BAPTA in recording pipette to tightly buffer the concentration of free calcium to be 50, 100 and 500 nM. HEK293 cells were expressed with TRPC4β alone or together with PLCδ1, and channels were stimulated with 100 nM EA. Channel currents were recorded with ramp pulse protocol from +100 to –100 mV every 10 sec. With 50 nM [Ca2+]i buffered internal solution, TRPC4β alone or TRPC4β and PLCδ1 expressing cells produced currents slowly after EA stimulation (Fig. 3A). In this condition, the presence of PLCδ1 did not significantly decrease the TRPC4β channel’s peak current density (Fig. 3B). By increasing [Ca2+]i to 100 nM, TRPC4β expressing cells showed relatively faster time course to reach the peak current (Fig. 3C). Notably, the EA-evoked current was significantly decreased in the cells expressing PLCδ1 together with TRPC4β (Fig. 3D). These aspects of decreased current in the cells expressing PLCδ1 were also observed in the experiments with 500 nM [Ca2+]i pipette solutions (Fig. 3F). Interestingly, they showed robust current activation and desensitized quickly after reaching the peak current (Fig. 3E). Exceptionally, the time course of TRPC4 with or without PLCδ1 was different only in 100 nM calcium buffering condition. In the cells only expressing TRPC4β, higher internal calcium concentration potentiates channel to show a faster time course, and 500 nM [Ca2+]i is sufficient to activate endogenous PLC to cause desensitization. PLCδ1 reduced the overall TRPC4β currents in 100 nM and 500 nM free Ca2+. Therefore, we suggest that physiological calcium concentration is sufficient to activate PLCδ1 after channel stimulation, and the activity of PLCδ1 is tightly regulated in a calcium dependent manner.
PLC is an enzyme that cleaves phospholipids in a Ca2+-dependent manner. TRPC4 channel is a nonselective cation channel which is permeable to calcium as well as monovalent cations. Therefore, Ca2+ influx-induced by TRPC4β activation may act as a stimulator of PLCδ1. To identify the activity of PLC, we used PH-domain of PLCδ1 as a PIP2 sensor to monitor the PIP2 changes. We co-transfected TRPC4β, PLCδ1, and CFP-tagged PH domain in HEK293 cells, and TRPC4β was stimulated by 100 nM EA. At first, CFP-PH was observed on the membrane, and line scan showed dense fluorescence intensity in the membrane before channel stimulation. In the cells expressing TRPC4β without PLCδ1, CFP-PH still showed similar intensity in the membrane even after 30 sec of EA stimulation and PIP2 started to deplete after 60 sec (Fig. 4A). On the other hand, TRPC4β and PLCδ1 co-expressed cells showed relatively faster PIP2 depletion from membrane. In these cells, 30 sec of EA stimulation was enough to deplete the membrane PIP2 (Fig. 4B). These results suggest that calcium influx
As a final approach, we evaluated the effect of PLCδ1 on TRPC4β. We used two mutant forms of PLCδ1 (Supplementary Fig. 2). The first mutant is PLCδ1 (K30A/K32A) which interrupted the PIP2 binding, as Lys30 and Lys32 residues are located in the PH-domain . And the second one is PLCδ1 (H311A) which shows abolished phospholipase activity since His311 residue is located in the X-region of catalytic domain (Fig. 5C) . We measured TRPC4β currents in the cells expressing TRPC4β and PLCδ1 (K30A/K32A) or PLCδ1 (H311A) mutant, respectively. Whole-cell patch clamp experiments were performed using the same protocol as mentioned before and the intracellular free calcium concentration was buffered to be 100 nM. In the cells expressing both mutants, the TRPC4β currents showed a similar time course as it would with cells expressing TRPC4β alone (Fig. 5A, B). Furthermore, EA-evoked currents in the cells expressing PLCδ1 (K30A/K32A) and PLCδ1 (H311A) with TRPC4β were significantly higher than in the cells expressing wild type PLCδ1, similar to the case when TRPC4β was expressed alone (Fig. 5D). In other words, PLCδ1 mutants which cannot regulate PIP2 levels did not have any effect on the TRPC4β activity. Therefore, these results suggest that TRPC4β-bound PLCδ1 regulates PIP2 level to decrease TRPC4β currents.
Activity of TRPC4 ion channel is closely controlled by PIP2. In this study, we found that PLCδ1 binds to TRPC4 and contribute to channel regulation. The conclusions of this study are as follows: 1. TRPC4 directly interacts with PLCδ1, but TRPC5 has no interaction with any PLCδ subtypes, 2. PLCδ1 causes inhibition of TRPC4 currents in a calcium dependent manner, 3. TRPC4-bound PLCδ1 responds to the calcium influx by the channel, which in turn depletes PIP2. Altogether, we propose a negative feedback regulation of TRPC4 by PLCδ1 in a Ca2+ and PIP2 dependent manner.
Previously, we suggested that Gαq-PLC pathway has a self-limiting activation on TRPC channels. The Gαq protein strongly activates TRPC heteromeric channels , and sequential events such as PIP2 depletion, PKC activation, and Ca2+ increase result in channel inhibition . We emphasized several times that regulation of PIP2 level is the most important factor as they are a requisite for channel activity. Accordingly, this important role of PLC and PIP2 on TRPC4 channel is not surprising. Several TRP channels are also known to be regulated
The direct interaction between TRPC4β and PLCδ1 also appeared clearly in expression patterns. PLCδ1 is expressed in both membrane and cytoplasm, and PLCδ3 is expressed in membrane. In our hands, PLCδ1 showed denser fluorescence intensity where TRPC4β is present, and they showed co-localization. Interestingly, PLCδ3 was observed to be empty in the channel puncta region when it was expressed with TRPC4β or TRPC5 (Fig. 1A, C). As PLC is a membrane-associated enzyme and TRPC channels are transmembrane channels, this expression pattern implies that the membrane resident proteins are present on the membrane competitively. Also, this difference presumably reflects that TRPC4β channel has specificity in PLC interaction, and membrane molecules cannot co-exist on the membrane unless they are interacting together.
In this study, we highlight that PLCδ1 binds to TRPC4β but not TRPC5, and it affects PIP2 level to regulate channel activity. TRPC4 and TRPC5 have a high similarity in amino acid sequence, and thus their structures are almost identical. However, they have some differences in their characteristics. First, TRPC5 shows constitutive activity while TRPC4 does not show channel activity in the absence of stimulation. Second, TRPC4 has a higher affinity for PIP2 than TRPC5. In a preceding study, TRPC5 showed robust current inhibition with weak voltage stimulation which activates voltage-sensitive phosphatase. Based on our data, we suggest that PIP2 and PLCδ1 are possible candidates underlying these features . TRPC4β-bound PLCδ1 would continuously regulate PIP2 level in the vicinity of the channel even without the stimulation, and it would attribute to TRPC4β showing no basal activity. In our hands, TRPC5 interacts with neither PLCδ1 nor PLCδ3, thus relatively abundant PIP2 pool would exist around TRPC5. The detailed knowledge of mechanism causing differences between TRPC4β and TRPC5 is important in understanding how they differently function physiologically in human body, but it should be elucidated in more detail.
Electrophysiological experiments suggest that PLCδ1 activity produces decreased TRPC4β currents in a calcium dependent manner. In the presence of PLCδ1, EA-evoked TRPC4β currents were significantly reduced (Fig. 3). Whether calcium increase alone causes the optimal activation of PLCδ1 is still controversial. It was suggested that high extracellular potassium, or stimulations like thapsigargin or ionomycin also causes an increase in IP3, but much less than the amount that is increased by the fully activated PLCδ subtypes [14,15,33]. Through experiments with differing intracellular calcium concentrations, our data agree with the notion that the channel-induced calcium influx within the physiological calcium range is sufficient to activate PLCδ1. In our hands, 50 nM of [Ca2+]i condition was not enough to fully activate PLCδ1, but PLCδ1 significantly reduced channel currents in higher calcium conditions (100 and 500 nM). Thus, we mainly used 100 nM [Ca2+]i buffering solution in electrophysiological experiments which showed prominent differences including reduced current amplitude and altered time course of currents in the presence of PLCδ1. Interestingly, the time course of channel currents was changed depending on the presence of PLCδ1 (Fig. 3C) or intracellular calcium concentration (Fig. 3E). It is well known that the TRPC4 and TRPC5 channels show desensitization after they are activated by GqPCR pathway [34-36]. Therefore, such desensitization with PLCδ1 and calcium might also explain the inclusion of PIP2 in this negative regulation mechanism. Our data showed that the kinetics of the process was much slower with low calcium which has a latency to respond to the channel agonist (Fig. 3). This observation is coincident with the fact that intracellular calcium strongly potentiates TRPC channels . However, this reaction rate was not influenced by the presence of PLCδ1 in the same calcium condition, hence we concluded that when PIP2 levels are diminished with PLCδ1, equivalent stimulus intensity will be less effective.
By monitoring the PIP2 level changes after TRPC4β stimulation, we suggest that PLCδ1 accelerates PIP2 hydrolysis in response to calcium influx (Fig. 4). As a second approach to demonstrate the involvement of PIP2 in this mechanism, we used PLCδ1 mutant which cannot hydrolyze PIP2 (Fig. 5). We expected that the PLCδ1-reduced TRPC4β currents might not be observed when PLCδ1 is unable to degrade PIP2. In electrophysiological experiments, these PLCδ1 mutants did not affect the channel currents. Therefore, it is likely that the PIP2 is the main source of negative feedback regulation by PLCδ1, and we emphasize the PIP2 as a requisite component for TRPC4β activation in this negative feedback mechanism of PLC.
Until now, we only focused on PIP2 as an indispensable component for TRPC4 channel activity, and the channel seems to fall into desensitization when PIP2 level is decreased. It should be noted that other groups have proposed the inhibitory effect of PIP2 on TRPC4 [6,37]. Recent evidence suggests that PIP2 has a dual regulatory role, and yet it remains an open question as to how PIP2 decides its role. TRPV1 has also been proposed to be regulated by PIP2 with a duality in function. Initially, TRPV1 channel was suggested to be potentiated after PIP2 hydrolysis by releasing inhibitory PIP2 from the channel [18,38-40].
Previously, it is suggested that the activity of PLCδ1 does not exhibit receptor-specific activation . In the electrophysiological experiments, we used channel-specific and strong agonist, hence we suggest that PLCδ1 is activated by channel-mediated calcium. However, it is still obscure how EA is involved in the signaling pathway that activates TRPC4β. The other group specified that PLCδ1 is involved in the Gαi/o-mediated TRPC4 activation. It should be elucidated whether the activation of PLCδ1 is confined to the specific pathway or the calcium through the channel itself is enough.
Collectively, our data reveal that TRPC4β-interacting PLCδ1 is sensitively activated by the calcium influx through TRPC4β and in turn, hydrolyzes TRPC4β-PI(4,5)P2 which maintains TRPC4β activity. Consequently, PLCδ1 is a negative regulator of TRPC4β.
Supplementary data including two figures can be found with this article online at https://doi.org/10.4196/kjpp.2023.27.2.187.
We thank Jung Eun Lee (Seoul National University, Seoul, Korea) and Christine Haewon Park (Seoul National University, Seoul, Korea) for helpful discussion and proofreading of manuscript. Careful proofreading was done by native speakers of English.
This work was supported by the National Research Foundation of Korea (NRF) grants funded by the Korean government (2020R1A2C1012670, 2021R1A4A2001857) (I.S.) and the Education and Research Encouragement Fund of Seoul National University Hospital (I.S.). J. Kim. was supported by the BK21 program from the Ministry of Education.
The authors declare no conflicts of interest.
|View Full Text||Article as PDF|
|Abstract||Figure & Table|
|Print this Page||Export to Citation|
ⓒ 2019. The Korean Journal of Physiology & Pharmacology. Powered by INFOrang Co., Ltd