Figure 1.
The motion of the national and regional economies.
Citation: Byron Carpenter. Current applications of mini G proteins to study the structure and function of G protein-coupled receptors[J]. AIMS Bioengineering, 2018, 5(4): 209-225. doi: 10.3934/bioeng.2018.4.209
| [1] | Massimiliano Calvia . Fossil energy use and carbon emissions: An easy-to-implement technical policy experiment. Green Finance, 2024, 6(3): 407-429. doi: 10.3934/GF.2024016 |
| [2] | Yuhang Zheng, Shuanglian Chen, Nan Wang . Does financial agglomeration enhance regional green economy development? Evidence from China. Green Finance, 2020, 2(2): 173-196. doi: 10.3934/GF.2020010 |
| [3] | Hatice Şimşek, Gökhan Öztürk . Evaluation of the relationship between environmental accounting and business performance: the case of Istanbul province. Green Finance, 2021, 3(1): 46-58. doi: 10.3934/GF.2021004 |
| [4] | Vladimir Kurdyukov, Alla Ovcharenko . Assessing the economic efficiency of green investments in the context of sustainable development of the territory. Green Finance, 2023, 5(3): 343-372. doi: 10.3934/GF.2023014 |
| [5] | Sabuj Saha, Ahmed Rizvan Hasan, Kazi Rezwanul Islam, Md Asraful Islam Priom . Sustainable Development Goals (SDGs) practices and firms' financial performance: Moderating role of country governance. Green Finance, 2024, 6(1): 162-198. doi: 10.3934/GF.2024007 |
| [6] | Mengze Wu, Dejun Xie . The impact of ESG performance on the credit risk of listed companies in Shanghai and Shenzhen stock exchanges. Green Finance, 2024, 6(2): 199-218. doi: 10.3934/GF.2024008 |
| [7] | Inna Semenenko, Yana Bilous, Ruslan Halhash . The compliance of the regional development strategies and funding with the sustainable development concept: The case of Ukraine. Green Finance, 2022, 4(2): 159-178. doi: 10.3934/GF.2022008 |
| [8] | Mukul Bhatnagar, Sanjay Taneja, Ercan Özen . A wave of green start-ups in India—The study of green finance as a support system for sustainable entrepreneurship. Green Finance, 2022, 4(2): 253-273. doi: 10.3934/GF.2022012 |
| [9] | Gonzalo Maldonado-Guzmán . Green innovation mediates between financial innovation and business sustainability? Proof in the mexican manufacturing industry. Green Finance, 2024, 6(3): 563-584. doi: 10.3934/GF.2024021 |
| [10] | Yicong Huang, Kaidong Yu, Chao Huang . Green finance engagement: An empirical study of listed companies on Chinese main board. Green Finance, 2023, 5(1): 1-17. doi: 10.3934/GF.2023001 |
Regional agglomeration and environmental changes are main features of contemporary economies. Moreover, economies experience business cycles during development processes. It is a greatly challenging task for economists to develop genuine dynamic models which exhibit business cycles, economic growth, and spatial agglomeration on microeconomic foundation. This paper is to develop a nonlinear multi-regional dynamic model with wealth and pollutant accumulation over time. We generalize the dynamic model by Zhang (2016). The original model examines multi-regional dynamics with wealth and pollutant accumulation over time. This study generalizes the model by making all exogenous constant parameters as exogenous time-dependent parameters. We demonstrate business cycles due to different exogenous periodic perturbations.
This study is concerned with regional agglomeration commonly observed in developed economies as well as in many underdeveloped economies. Kuznets (1966) discussed spatial unequal economic development over space long time ago. He examined interdependence between spatial agglomeration and industrialization. Myrdal (1957) and Hirschman (1958) earlier argued for dynamic interactions between industrial growth and the geographical concentration of industry. There are numerous studies on various aspects of spatial economic developments (e.g., Krugman, 1991; Fujita et al., 1999; Forslid and Ottaviano, 2003). Different from earlier literature on spatial agglomeration, Zhang (2016) extends neoclassical growth theory to multi-regional economies. This approach is different from the so-called new economic geography (e.g., Krugman, 1993; and Tabuchi, 2014). In the new economic geography, the role of physical capital is omitted in spatial dynamics. Moreover, In the literature of the new economic geography, Tabuchi (2014: 50) reviews: "The scopes of most of the theoretical studies published thus far have been limited to two regions in order for researchers to reach meaningful analytical results. The two-region NEG models tend to demonstrate that spatial distribution is dispersed in the early period (high trade costs or low manufacturing share) and agglomerated in one of the two regions in the late period (low trade costs or high manufacturing share). However, it is no doubt that the two-region NEG models are too simple to describe the spatial distribution of economic activities in real-world economies. Since there are only two regions, their geographical locations are necessarily symmetric, and thus diverse spatial distributions cannot occur." This study examines the model with any number of regions in order to properly address issues related to interregional growth and agglomeration.
Dynamic interdependence between economic growth and environmental changes is extensively studied (e.g., Ploude, 1972; John and Pecchenino, 1994; Gassebner et al., 2011; Lin and Liscow, 2012; Apergis and Payne, 2014; Yang and Zhao, 2014; Antonakakis et al., 2017; and Walenta, 2018). It is pointed out that theoretical economics fails to address many issues related economic growth and environmental change with a comprehensive framework with microeconomic foundation. As pointed out by Fullerton and Kim (2008), existing research has proposed a number of different models for analyzing different questions not in an integrated way. This study deals with growth, interregional trade, and environmental change within an integrated framework. It should be noted that following Zhang (1996), we also consider amenity endogenous. Amenity is assumed to be related to environmental quality (and other factors such as historical building and regional population, residential densities) and household location (e.g., Glaeser et al., 2001; Partridge et al., 2008; Chen et al., 2013; Lai, 2019; and Tivadar and Jayet, 2019). Zhang (1993) first introduced spatial amenity into utility in a general equilibrium framework. This study incorporates amenity into the consumer location decision. According to Chen et al. (2013), "Empirical evidence demonstrates a growing link between the presence of high-valued natural amenities—including pleasant climate and proximity to lakes, oceans, forests, and mountains—and higher rates of population and income growth in the U.S." Chen et al. (2013) developed a two-region model in which labor distribution, production externalities and natural resources are endogenous. According to their research, preference for natural amenities tends to lead population dispersion. This study shows how amenity is related to regional agglomeration.
Business cycles occur due to exogenous or endogenous. For instance, real business-cycle theory treats business cycles responses to various exogenous changes in the real economic variables. In the literature of economic theory, only a few formal models can properly deal with exogenous cycles on basis of microeconomics. This study introduces exogenous real shocks to the neoclassical growth model to demonstrate the existence of business cycles due to real shocks. It should be noted that in modern nonlinear dynamic economic theory business cycles occur with or without any exogenous changes (e.g., Zhang, 1991, 2005, 2006; Lorenz, 1993; Flaschel et al. 1997; Chiarella and Flaschel, 2000; Shone, 2002; Gandolfo, 2005; and Puu, 2011). Although there are many models of business cycles for national economies, there are a few nonlinear multi-regional dynamic models with microeconomic foundation for identifying business cycles. Different theories emphasize various origins of economic periodic changes. The rest of the study is organized as follows. Section 2 extends Zhang's model to allow for existence of different time-dependent changeable parameters. In Section 3 we demonstrate that the movement of the economy is described by differential equations. The section also simulates the model for a three-regional economy. Section 4 studies comparative dynamic analysis in some exogenous periodic oscillations. Section 5 concludes the study.
This section is to generalize Zhang's model by making all exogenous constant parameters to be exogenous time-dependent (Zhang, 2016). The economy is composed of J regions, indexed by $j = 1, \ldots, J$ Each region has one production sector and one environment sector. The production sector is the same as in the the Solow growth model. The system produces one (durable) good. Assets are owned by households. Capital and labor as inputs of the production sector. Productivities of the production sector fully utilized. Saving is undertaken only by households. Commodities are traded without any barriers. We neglect transport costs. We measure prices in terms of the commodity and the price of the commodity be unity. We represent wage and interest rates by $w_{j}(t)$ and $r_{j}(t)$ respectively, in the jth region. The interest rate is equal across regions, i.e., $r(t) = r_{j}(t)$ We consider that the national population N(t) is homogenous. Workers are free to choose their work place and residential location. Workers work and reside in the same region. Each region has land $L_{j}(t)$ which is used only for residential purposes. Regions have different productivities and living conditions. We use subscripts, i, e, to denote the production and environment sectors, respectively. Let $F_{q j}(t)$ stand for the output levels of q's sector in region j at time $t, q = i, e$. Each region has one environment sector. The role of the regional government is to collect different taxes on producers and consumers to financially support the sector. The tax income is used up for employing labor and capital to purify environment.
The production sectors use labor and physical capital inputs. Environment affects productivities. The production functions are specified as follows:
|
$
Fij(t)=Aij(t)Γij(Ej(t))Kαij(t)ij(t)Nβij(t)ij(t),Aij(t),αij(t),βij(t)>0,αij(t)+βij=1
$
|
(1) |
where Aij(t), αij(t), and βij(t) are positive parameters. Here, Γij(Ej(t)) is a function of the environmental quality measured by the level of pollution, Ej(t), in region j. It is assumed that the productivity of the capital goods sector is non-positively related to the pollution level. Let δkj stand for the depreciation rate of physical capital in region j. We use τij(t) to stand for the tax rate on the capital goods sector. The marginal conditions of the production sector are:
| $r\left(t \right) + {\delta _{kj}}\left(t \right) = \frac{{{\alpha _{ij}}\left(t \right){{\bar \tau }_{ij}}\left(t \right){F_{ij}}\left(t \right)}}{{{K_{ij}}\left(t \right)}}, {w_j}\left(t \right) = \frac{{{\beta _{ij}}\left(t \right){{\bar \tau }_{ij}}\left(t \right){F_{ij}}\left(t \right)}}{{{N_{ij}}\left(t \right)}}$ | (2) |
where $\overline{\tau}_{i j}(t) \equiv 1-\tau_{i j}(t), 0 < \tau_{i j}(t) < 1$.
Following Zhang (2016) we specify the dynamics of the stock of pollutants as follows:
| ${\dot E_j}\left(t \right) = {\theta _{ij}}\left(t \right){F_{ij}}\left(t \right) + {\theta _{cj}}\left(t \right){C_j}\left(t \right) - {F_{ej}}\left(t \right) - {\bar \theta _j}\left(t \right){E_j}\left(t \right) + {\Omega _j}\left({\left({{E_q}\left(t \right)} \right)} \right){\rm{}}$ | (3) |
in which θijt, θcj(t) and θj(t) are positive parameters and
|
$Fej(t)=Aej(t)Γej(Ej(t))K˜αej(t)ej(t)N˜βejej(t),Aej(t),˜αej(t),˜βej(t)>0 $
|
(4) |
where Aej(t), αej(t), and βej(t) are positive parameters, and Γej(Ej)(≥0) is a function of Ej(t). The term θij(t) Fij(t) means that pollutants that are emitted during production processes are linearly positively proportional to the output level. The parameter, θij(t), means that in consuming one unit of the good the quantity θij(t) is left as waste. The creation of pollutants by consumers is given θcj(t) Cj(t). The parameter θj(t) is called the rate of natural purification.
We take account of effects of transboundary pollution by ${\mathit{\Omega }_j}\left({\left({{E_q}(t)} \right)} \right)$ where is a function of pollutants of all the regions.
All the income from land is evenly redistributed among the population. Region j's land is Lj(t) and the land rent is Rj(t). The land revenue r(t) per worker is as follows:
| $\overline{r}(t) = \frac{1}{N} \sum_{j = 1}^{J} L_{j}(t) R_{j}(t)$ | (5) |
Let kj(t) stand for per capita wealth of region j. We have kj(t) = Kj(t)/Nj(t), where K(t) is the total wealth held by region j. We use τkj(t) and τwj(t) to stand for the tax rates on the wealth income and wage income in region j, respectively. Per capita current income from the land ownership, the interest payment r(t)kj(t), and the wage payment wj(t), is
| $y_{j}(t) = \overline{r}(t)+\overline{\tau}_{k j}(t) r(t) \overline{k}_{j}(t)+\overline{\tau}_{w j}(t) w_{j}(t)$ | (6) |
where τkj(t)≡1-τkj(t) and τwj(t) ≡1-τwj(t). The per sable income an disposable income is the sum of the current disposable income and the value of wealth:
| $\hat{y}_{j}(t) = \overline{k}_{j}(t)+y_{j}(t)$ | (7) |
The disposable income is used for saving and consumption. At each point of time, a consumer would distribute the total available budget between housing, lj(t), saving sj(t) and consumption cj(t). The budget constraint is given by
| $\left(1+\tau_{R j}(t)\right) R_{j}(t) l_{j}(t)+\left(1+\tau_{c j}(t)\right) p_{j}(t) c_{j}(t)+s_{j}(t) = \hat{y}_{j}(t)$ | (8) |
where τRj(t) is the tax rate on housing and τcj(t) is the tax rate on consumption.
Utility level Uj(t) is dependent on the consumption levels of lot size and commodity, and saving. The utility level of the typical consumer in region j is represented by
| ${U_j}\left(t \right) = {\theta _j}\left(t \right)l_j^{{\eta _0}\left(t \right)}\left(t \right)c_j^{{\xi _0}\left(t \right)}s_j^{{\lambda _0}\left(t \right)}\left(t \right), {\eta _0}\left(t \right), {\xi _0}\left(t \right), {\lambda _0}\left(t \right) \gt 0$ | (9) |
in which η0(t), ξ0(t), and λ0(t) are the elasticities of utility with regard to lot size, commodity and saving. We call η0(t), ξ0(t), and λ0(t) propensities to consume lot size, and goods, and to hold wealth, respectively. In (9), θj(t) is region j's amenity level (e.g., Diamond and Tolley, 1981; Landry and Chakraborty, 2009; Lanz and Provins, 2013; Debarsy and Ertur, 2019; and Saracoglu and Roe, 2019). We specify θj(t) as follows:
| $\theta_{j}(t) = \overline{\theta}_{j}(t) E_{j}^{\omega(t)}(t) N_{j}^{d(t)}(t), j = 1, \ldots, J$ | (10) |
where θj(t) (>0), ω(≤0), and d(t) are parameters, and Nj(t) is region j's population. We don't specify signs of d(t) as population may have either positive or negative effects on regional attractiveness.
Maximizing Uj(t) subject to (8) yields
| $l_{j}(t) R_{j}(t) = \eta_{j}(t) \hat{y}_{j}(t), c_{j}(t) = \xi_{j}(t) \hat{y}_{j}(t), s_{j}(t) = \lambda(t) \hat{y}_{j}(t)$ | (11) |
where
| ${\eta _j}\left(t \right) \equiv \frac{{{\eta _0}\left(t \right)\rho \left(t \right)}}{{1 + {\tau _{Rj}}\left(t \right)}}, {\xi _j}\left(t \right) \equiv {\rm{}}\frac{{{\xi _0}\left(t \right)\rho \left(t \right)}}{{1 + {\tau _{cj}}}}, \lambda \left(t \right) \equiv {\lambda _0}\left(t \right)\rho \left(t \right)$ | (12) |
| $\rho \left(t \right) \equiv \frac{1}{{{\eta _0}\left(t \right) + {\rm{}}{\xi _0}\left(t \right) + {\lambda _0}\left(t \right)}}$ | (13) |
From (11) we see
| $s_{j}(t) = \lambda(t)\left[\overline{r}(t)+\overline{k}_{j}(t)+\overline{\tau}_{k j}(t) r(t) \overline{k}_{j}(t)+\overline{\tau}_{w j}(t) w_{j}(t)\right]$ | (14) |
According to the definitions of sj(t), the wealth accumulation of the representative household in region j is:
| $\dot{\overline{k}}_{j}(t) = s_{j}(t)-\overline{k}_{j}(t)$ | (15) |
We assume that households are freely mobile between the regions. The utility levels of homogenous population between regions are assumed to be equal at any point in time. That is
| $U_{j}(t) = U_{q}(t), \quad j, q = 1, \ldots, J$ | (16) |
Each region's governmental tax income consists of the tax incomes on the production sector, consumption, wage income and wealth income as follows:
|
$Yej(t)=τij(t)Fij(t)+τRj(t)lj(t)Rj(t)Nj(t)+τcj(t)cj(t)Nj(t)+τwj(t)wj(t)Nj(t)+τkj(t)r(t)ˉkj(t)Nj(t) $
|
(17) |
The government's income is used up only for the environmental purpose. The government budget is as follows:
| $\left({r\left(t \right) + {\delta _k}\left(t \right)} \right){K_{ej}}\left(t \right) + {w_j}\left(t \right){N_{ej}}\left(t \right) = {Y_{ej}}\left(t \right)$ | (18) |
The government employs labor and capital stocks for purifying environment in such a way that the purification rate is maximized under the given budget constraint. The government's optimal problem is thus formed by:
| $\mathop {{\rm{Max}}}\limits_{\left\{ {{K_{ej}}\left(t \right), {N_{ej}}\left(t \right)} \right\}} {F_{ej}}\left(t \right)\;\;\;s.t.:$ | (19) |
The optimal solution is:
| $\left({r\left(t \right) + {\delta _k}\left(t \right)} \right){K_{ej}}\left(t \right) = {\alpha _{ej}}\left(t \right){Y_{ej}}\left(t \right), {w_j}\left(t \right){N_{ej}}\left(t \right) = {\beta _{ej}}\left(t \right){Y_{ej}}\left(t \right)$ | (20) |
where
| $\alpha_{e j}(t) \equiv \frac{\widetilde{\alpha}_{e j}(t)}{\widetilde{\alpha}_{e j}(t)+\widetilde{\beta}_{e j}(t)}, \quad \beta_{e j}(t) \equiv \frac{\widetilde{\beta}_{e j}(t)}{\widetilde{\alpha}_{e j}(t)+\widetilde{\beta}_{e j}(t)}$ | (21) |
Each region's capital stocks Kj(t) is fully
| $K_{i j}(t)+K_{e j}(t) = K_{j}(t)$ | (22) |
The total capital stock K(t) employed by the production sectors is equal to the total wealth owned by the households of all the regions:
| $K(t) = \sum_{j = 1}^{J} K_{j}(t) = \sum_{j = 1}^{J} \overline{k}_{j}(t) N_{j}(t)$ | (23) |
The assumption that labor force and land are fully employed is represented by:
| $N_{i j}(t)+N_{e j}(t) = N_{j}(t), \Sigma_{j = 1}^{J} N_{j}(t) = N(t)$ | (24) |
| $l_{j}(t) N_{j}(t) = L_{j}(t), j = 1, \ldots, J$ | (25) |
We have thus built the model.
The economic system contains many variables which are nonlinearly interrelated. No general analytical method is available to solve such nonlinear differential equations. For illustration, we simulate the model. In the appendix, we show that the dynamics of the national economy can be expressed as differential equations. First, we introduce a variable z1(t)
| $z_{1}(t) \equiv \frac{r(t)+\delta_{k 1}(t)}{w_{1}(t)}$ | (26) |
The following lemma provides the procedure to simulate the model.
The motion of the national economy is given by the following 2J+2 differential equations with $z_{1}(t), \overline{r}(t), \left(\overline{k}_{j}(t)\right)$ and $\left(E_{j}(t)\right)$ as variables
|
$
˙ˉkj(t)=Φj(z1(t),ˉr(t),(ˉkq(t)),(Eq(t)),t),˙Ej(t)=Λj(z1(t),ˉr(t),(ˉkq(t)),(Eq(t)),t),j=1,…,J,˙z1(t)=ˉΦ(z1(t),ˉr(t),(ˉkq(t)),(Eq(t)),t),˙ˉr(t)=ˉΛ(z1(t),ˉr(t),(ˉkq(t)),(Eq(t)),t)
$
|
(27) |
where ${\mathit{\Phi }_j}(t), \mathit{\bar \Phi }(t), {\mathit{\Lambda }_j}(t)$ and $\mathit{\bar \Lambda }(t)$ are functions of $t, \quad z_{1}(t), \overline{r}(t), \quad\left(\overline{k}_{j}(t)\right)$ and $\left(E_{j}(t)\right)$ defined in the appendix. For any given positive values of $z_{1}(t), \overline{r}(t), \left(\overline{k}_{j}(t)\right)$ and $\left(E_{j}(t)\right)$ at any point in time, the other variables are uniquely determined by the following procedure: r(t) by (A2) → wj(t) by (A4) →$\hat{y}_{j}(t)$ by (A5) → Rj(t) by (A6) → N1(t) by (A11) → Nj(t) by (A15) → Nej(t) and Nij(t) by (A14) → Kej(t) and Kij(t) by (A1) → Fij(t) by (4) → Fej(t) by (4) → cj(t), lj(t), and sj(t) by (11) → θj(t) by (10) →Uj(t) by (9) → Yej(t) by (A12) → Kj(t) = Kij(t)+Kej(t) → K(t) by (23).
We simulate the case that the economy is composed of three regions. We specify the functions dependent on environmental quality as follows
| $ {\mathit{\Gamma }_{mj}}\left({{E_j}(t), t} \right) = E_j^{ - {b_{mj}}(t)}(t), \;\;{\mathit{\Gamma }_j}\left({{E_j}(t), t} \right) = E_j^{ - {b_j}(t)}(t), \;\;j = 1, 2, 3, \;\;m = i, e $ | (28) |
We require bji, bjs, bj ≥ 0 and bje ≤ 0. We specify transboundary pollution functions as:
| $ {\mathit{\Omega }_q}\left({\left({{E_q}(t)} \right), t} \right) = \sum{j, j \ne q}^J {{\theta _{jq}}} (t){E_j}(t) $ | (29) |
The transboundary pollution functions imply that any region is also polluted by other two regions and the speed is linearly related to the pollutant levels of other regions. The rest of this section is concerned with the case that all the parameters are constant. This case is investigated by Zhang (2016) and we simply state the results which are the reference point for following the motion of the system when some parameters are time-dependent. We specify the parameters as follows:
|
$N=20,λ0=0.75,ξ0=0.1,η0=0.1,d=−0.05,ω=−0.08,τij=τcj,θjq=0.01,j≠q,j,q=1,2,3,(L1L2L3)=(343),(Ai1Ai2Ai3)=(1.210.8),(Ae1Ae2Ae3)=(1.110.9),(αi1αi2αi3)=(0.310.310.31),(bi1bi2bi3)=(be1be2be3)=(0.010.020.01),(ˉθ1ˉθ2ˉθ3)=(0.080.10.11),(τk1τk2τk3)=(τw1τw2τw3)=(0.0050.0060.007),(˜αej˜βejτij)=(0.40.20.005),(bjθijθcj)=(0.010.080.03),(δk1δk2δk3)=(0.050.040.04) $
|
(30) |
We plot the motion of the system under (30) with the following initial conditions:
|
$
z1(0)=0.038,ˉr(0)=1.42,ˉk1(0)=12.7,ˉk2(0)=11,ˉk3(0)=9.7,E1(0)=8.3,E2(0)=9.8,E3(0)=6.4.
$
|
The motion of the variables is plotted in Figure 1.
From the plot we see that the system becomes stationary over time. This implies the existence of an equilibrium point. We list the equilibrium values as follows:
|
$
r = 0.001, \;\;\bar r = 1.46, \;\;K = 219.7, \\
(N1N2N3)=(3.458.607.95),(K1K2K3)=(49.98115.7280.71),(E1E2E3)=(8.059.926.52),(Fi1Fi2Fi3)=(9.0318.0912.57),(Ki1Ki2Ki3)=(42.46111.7177.59),(Ni1Ni2Ni3)=(3.428.547.88),(Fe1Fe2Fe3)=(0.420.990.78),(Ke1Ke2Ke3)=(1.524.013.12),(Ne1Ne2Ne3)=(0.0250.0690.071),(R1R2R3)=(1.943.233.49),(w1w2w3)=(1.811.461.1),(θ1θ2θ3)=(0.0640.0750.085),(ˉk1ˉk2ˉk3)=(12.711.39.92),(l1l2l3)=(0.870.470.38),(c1c2c3)=(1.681.501.32)
$
|
(31) |
It is straightforward to calculate the eigenvalues as follows
| $\{-0.18 \pm 0.02 I, -0.2, -0.15, -0.1, \quad-0.07, -0.005, 0\}.$ |
The real parts of the seven eigenvalues are negative and one is zero. The unique equilibrium is locally stable. This implies the validity of conducting comparative dynamic analysis.
We simulated the motion of the national economy under (30). This section identifies business cycles by allowing various exogenous periodic shocks. We use $\overline{\Delta} \mathrm{x}_{\mathrm{j}}(\mathrm{t})$ to stand for the change rate of the variable, $\mathrm{x}_{\mathrm{j}}(\mathrm{t})$ in percentage due to changes in a parameter value.
We now study the effects of the following periodic oscillations in the national population:
| $N(t) = 20+0.2 sin (t).$ |
The simulation results are plotted in Figure 2. The change rates in the regional populations are regionally oscillatory. The per household consumption and wealth levels cross the regions are slightly affected, even though the other variables changes greatly in amplitude. The national economy experiences business cycles over time.
The population may affect amenity in different ways. In this study we consider that the population has negative impact on the amenity. We assume parameter d equal cross regions. We now examine the effects when the negative impact of the population is periodically perturbed as follows:
| $d(t) = -0.05+0.02 sin (t).$ |
The simulation results are plotted in Figure 3. Although the regional variables change periodically, region 1's macro variables vary largely in amplitude.
We now examine the effects when the negative impact of pollutants on amenities is oscillating as follows: $\omega(t) = -0.08+0.04 sin (t).$
The simulation results are plotted in Figure 4. It is worthwhile to mention that the regional macroeconomic variables change, but the amplitudes are much smaller than the other two regions'. No region's microeconomic variables are strongly affected.
We now examine the effects of a rise in the MPR's tax rate on consumption as follows:
| $\tau_{c 1}(t) = 0.005+0.004 sin (t).$ |
The simulation results are plotted in Figure 5. The MPR spends more resource on purifying environment. The two input factors and output level of the MPR's sector are increased. The level of pollutants is reduced and the region's amenity is enhanced. As pollutants are reduced in the MPR, the productivity of the region is improved. The better environment attracts more people to live in the MPR. The lot size falls in the MPR and the land rent rises. The PR's and LPR's populations are enhanced. The lot sizes in the PR and the LPR fall. The two input factors and output levels of the PR's and LPR's two sectors are reduced. The levels of pollutants are enhanced and the regions' amenities are improved. The consumption and wealth levels in all the regions are improved. We see that a rise in the tax rate brings positive effects on all aspects of regional and national economies.
We now analyze the effects of the following oscillations in region 3's consumption on environment: $\theta_{c 3}(t) = 0.03+0.02 sin (t).$
The simulation results are plotted in Figure 6.
In growth theory effects of saving propensity changes are different in different theories. In Keynesian economic theory savings tend to reduce national income, while neoclassical growth theory claims the opposite effect. We study what happen to the economic system when the propensity to save is oscillated as follows: $\lambda_{0}(t) = 0.75+0.03 sin (t).$
The simulation results are plotted in Figure 7. There are business cycles in both macro and economic variables when the propensity to save is oscillating.
This paper made a unique contribution to the literature of business cycles and multi-regional economic growth. As far as I am aware, there is almost mathematical economic model of business cycles with multiple regions and capital accumulation, and endogenous amenities which are built on microeconomic foundation. This paper showed how economic cycles appear across regions in association with spatial agglomeration and economic growth. It generalized the dynamic model of interdependence between regional economic growth and endogenous environment change proposed by Zhang (2016). The original model examined multi-regional dynamics with wealth and pollutant accumulation over time. The dynamics of J-region economy is described by differential equations. This paper generalized the model by making all exogenous constant parameters as exogenous time-dependent parameters. We demonstrate business cycles due to different exogenous periodic perturbations. The unique contribution to the literature of business cycles theory is that it demonstrates existence of business cycles in a multi-regional neoclassical growth model with endogenous environment built on microeconomic foundation. As the model takes account of various sources of tax, it is possible to study the impact of changes in the tax policy on growth and oscillations in the economic system. Although this study does not take account of dynamics of political parties in a democratic society, our results on oscillations in the tax policies which might be due to powershifts tell us that business cycles can also result from political competition. As the model is built on some well-established economic theories and each theory has a vast amount of literature, the model may be extended and generalized in multiple directions.
The author is grateful to the constructive comments of the two anonymous referees.
I, the author, declare that there is no conflict of interests regarding the publication of this paper.
| [1] |
Wise A, Gearing K, Rees S (2002) Target validation of G-protein coupled receptors. Drug Discov Today 7: 235–246. doi: 10.1016/S1359-6446(01)02131-6
|
| [2] |
Cherezov V, Rosenbaum DM, Hanson MA, et al. (2007) High-resolution crystal structure of an engineered human beta2-adrenergic G protein-coupled receptor. Science 318: 1258–1265. doi: 10.1126/science.1150577
|
| [3] |
Serrano-Vega MJ, Magnani F, Shibata Y, et al. (2008) Conformational thermostabilization of the beta1-adrenergic receptor in a detergent-resistant form. Proc Natl Acad Sci USA 105: 877–882. doi: 10.1073/pnas.0711253105
|
| [4] | Caffrey M (2015) A comprehensive review of the lipid cubic phase or in meso method for crystallizing membrane and soluble proteins and complexes. Acta Crystallogr F Struct Biol Commun 71: 3–18. |
| [5] |
Landau EM, Rosenbusch JP (1996) Lipidic cubic phases: A novel concept for the crystallization of membrane proteins. Proc Natl Acad Sci USA 93: 14532–14535. doi: 10.1073/pnas.93.25.14532
|
| [6] | Chan W, Said M, Zhang C, et al. (2018) GPCR-EXP: A semi-manually curated database for experimentally-solved and predicted GPCR structures. Available from: https://zhanglab.ccmb.med.umich.edu/GPCR-EXP/. |
| [7] |
Jazayeri A, Dias JM, Marshall FH (2015) From G Protein-coupled receptor structure resolution to rational drug design. J Biol Chem 290: 19489–19495. doi: 10.1074/jbc.R115.668251
|
| [8] |
Carpenter B, Tate CG (2017) Active state structures of G protein-coupled receptors highlight the similarities and differences in the G protein and arrestin coupling interfaces. Curr Opin Struct Biol 45: 124–132. doi: 10.1016/j.sbi.2017.04.010
|
| [9] |
Rasmussen SG, Devree BT, Zou Y, et al. (2011) Crystal structure of the beta2 adrenergic receptor-Gs protein complex. Nature 477: 549–555. doi: 10.1038/nature10361
|
| [10] |
Kang Y, Zhou XE, Gao X, et al. (2015) Crystal structure of rhodopsin bound to arrestin by femtosecond X-ray laser. Nature 523: 561–567. doi: 10.1038/nature14656
|
| [11] |
Westfield GH, Rasmussen SG, Su M, et al. (2011) Structural flexibility of the G alpha s alpha-helical domain in the beta2-adrenoceptor Gs complex. Proc Natl Acad Sci USA 108: 16086–16091. doi: 10.1073/pnas.1113645108
|
| [12] |
Van EN, Preininger AM, Alexander N, et al. (2011) Interaction of a G protein with an activated receptor opens the interdomain interface in the alpha subunit. Proc Natl Acad Sci USA 108: 9420–9424. doi: 10.1073/pnas.1105810108
|
| [13] |
Scheerer P, Park JH, Hildebrand PW, et al. (2008) Crystal structure of opsin in its G-protein-interacting conformation. Nature 455: 497–502. doi: 10.1038/nature07330
|
| [14] |
Hamm HE, Deretic D, Arendt A, et al. (1988) Site of G protein binding to rhodopsin mapped with synthetic peptides from the alpha subunit. Science 241: 832–835. doi: 10.1126/science.3136547
|
| [15] |
Pardon E, Laeremans T, Triest S, et al. (2014) A general protocol for the generation of Nanobodies for structural biology. Nat Protoc 9: 674–693. doi: 10.1038/nprot.2014.039
|
| [16] | Carpenter B, Tate CG (2016) Engineering a minimal G protein to facilitate crystallisation of G protein-coupled receptors in their active conformation. Protein Eng Des Sel 29: 583–594. |
| [17] |
Blankenship E, Vahedi-Faridi A, Lodowski DT (2015) The high-resolution structure of activated opsin reveals a conserved solvent network in the transmembrane region essential for activation. Structure 23: 2358–2364. doi: 10.1016/j.str.2015.09.015
|
| [18] |
Choe HW, Kim YJ, Park JH, et al. (2011) Crystal structure of metarhodopsin II. Nature 471: 651–655. doi: 10.1038/nature09789
|
| [19] |
Deupi X, Edwards P, Singhal A, et al. (2012) Stabilized G protein binding site in the structure of constitutively active metarhodopsin-II. Proc Natl Acad Sci USA 109: 119–124. doi: 10.1073/pnas.1114089108
|
| [20] |
Park JH, Morizumi T, Li Y, et al. (2013) Opsin, a structural model for olfactory receptors? Angew Chem Int Ed Engl 52: 11021–11024. doi: 10.1002/anie.201302374
|
| [21] |
Singhal A, Ostermaier MK, Vishnivetskiy SA, et al. (2013) Insights into congenital stationary night blindness based on the structure of G90D rhodopsin. EMBO Rep 14: 520–526. doi: 10.1038/embor.2013.44
|
| [22] |
Standfuss J, Edwards PC, D'Antona A, et al. (2011) The structural basis of agonist-induced activation in constitutively active rhodopsin. Nature 471: 656–660. doi: 10.1038/nature09795
|
| [23] |
Huang W, Manglik A, Venkatakrishnan AJ, et al. (2015) Structural insights into micro-opioid receptor activation. Nature 524: 315–321. doi: 10.1038/nature14886
|
| [24] |
Kruse AC, Ring AM, Manglik A, et al. (2013) Activation and allosteric modulation of a muscarinic acetylcholine receptor. Nature 504: 101–106. doi: 10.1038/nature12735
|
| [25] |
Rasmussen SG, Choi HJ, Fung JJ, et al. (2011) Structure of a nanobody-stabilized active state of the beta(2) adrenoceptor. Nature 469: 175–180. doi: 10.1038/nature09648
|
| [26] |
Ring AM, Manglik A, Kruse AC, et al. (2013) Adrenaline-activated structure of beta2-adrenoceptor stabilized by an engineered nanobody. Nature 502: 575–579. doi: 10.1038/nature12572
|
| [27] |
Weichert D, Kruse AC, Manglik A, et al. (2014) Covalent agonists for studying G protein-coupled receptor activation. Proc Natl Acad Sci USA 111: 10744–10748. doi: 10.1073/pnas.1410415111
|
| [28] |
Manglik A, Kobilka BK, Steyaert J (2017) Nanobodies to study G protein-coupled receptor structure and function. Annu Rev Pharmacol Toxicol 57: 19–37. doi: 10.1146/annurev-pharmtox-010716-104710
|
| [29] |
Nehme R, Carpenter B, Singhal A, et al. (2017) Mini-G proteins: Novel tools for studying GPCRs in their active conformation. PLoS One 12: e0175642. doi: 10.1371/journal.pone.0175642
|
| [30] |
Carpenter B, Nehme R, Warne T, et al. (2016) Structure of the adenosine A(2A) receptor bound to an engineered G protein. Nature 536: 104–107. doi: 10.1038/nature18966
|
| [31] |
Hanzal-Bayer M, Renault L, Roversi P, et al. (2002) The complex of Arl2-GTP and PDE delta: From structure to function. EMBO J 21: 2095–2106. doi: 10.1093/emboj/21.9.2095
|
| [32] |
Sunahara RK, Tesmer JJ, Gilman AG, et al. (1997) Crystal structure of the adenylyl cyclase activator Gsalpha. Science 278: 1943–1947. doi: 10.1126/science.278.5345.1943
|
| [33] |
Spiegel AM, Jr BP, Butrynski JE, et al. (1991) The G protein connection: Molecular basis of membrane association. Trends Biochem Sci 16: 338–341. doi: 10.1016/0968-0004(91)90139-M
|
| [34] | Carpenter B, Tate CG (2017) Expression, purification and crystallisation of the adenosine A2A receptor bound to an engineered Mini G protein. Bio Protoc, 7. |
| [35] |
Wan Q, Okashah N, Inoue A, et al. (2018) Mini G protein probes for active G protein-coupled receptors (GPCRs) in live cells. J Biol Chem 293: 7466–7473. doi: 10.1074/jbc.RA118.001975
|
| [36] | Carpenter B, Tate CG (2017) Expression and purification of mini G proteins from escherichia coli. Bio Protoc 7: e2235. |
| [37] |
Garcia-Nafria J, Lee Y, Bai X, et al. (2018) Cryo-EM structure of the adenosine A2A receptor coupled to an engineered heterotrimeric G protein. Elife 7: e35946. doi: 10.7554/eLife.35946
|
| [38] |
Garcia-Nafria J, Nehme R, Edwards PC, et al. (2018) Cryo-EM structure of the serotonin 5-HT1B receptor coupled to heterotrimeric Go. Nature 558: 620–623. doi: 10.1038/s41586-018-0241-9
|
| [39] |
Tsai CJ, Pamula F, Nehmé R, et al. (2018) Crystal structure of rhodopsin in complex with a mini-Go sheds light on the principles of G protein selectivity. Sci Adv 4: eaat7052. doi: 10.1126/sciadv.aat7052
|
| [40] |
Lebon G, Bennett K, Jazayeri A, et al. (2011) Thermostabilisation of an agonist-bound conformation of the human adenosine A(2A) receptor. J Mol Biol 409: 298–310. doi: 10.1016/j.jmb.2011.03.075
|
| [41] |
Lebon G, Warne T, Edwards PC, et al. (2011) Agonist-bound adenosine A2A receptor structures reveal common features of GPCR activation. Nature 474: 521–525. doi: 10.1038/nature10136
|
| [42] |
Warne T, Serrano-Vega MJ, Baker JG, et al. (2008) Structure of a beta1-adrenergic G-protein-coupled receptor. Nature 454: 486–491. doi: 10.1038/nature07101
|
| [43] |
Lebon G, Edwards PC, Leslie AG, et al. (2015) Molecular determinants of CGS21680 binding to the human adenosine A2A receptor. Mol Pharmacol 87: 907–915. doi: 10.1124/mol.114.097360
|
| [44] |
Lebon G, Warne T, Tate CG (2012) Agonist-bound structures of G protein-coupled receptors. Curr Opin Struct Biol 22: 482–490. doi: 10.1016/j.sbi.2012.03.007
|
| [45] |
Xu F, Wu H, Katritch V, et al. (2011) Structure of an agonist-bound human A2A adenosine receptor. Science 332: 322–327. doi: 10.1126/science.1202793
|
| [46] |
Ye L, Van EN, Zimmer M, et al. (2016) Activation of the A2A adenosine G-protein-coupled receptor by conformational selection. Nature 533: 265–268. doi: 10.1038/nature17668
|
| [47] |
Strege A, Carpenter B, Edwards PC, et al. (2017) Strategy for the thermostabilization of an agonist-bound GPCR coupled to a G protein. Method Enzymol 594: 243–264. doi: 10.1016/bs.mie.2017.05.014
|
| [48] |
Carpenter B, Lebon G (2017) Human adenosine A2A receptor: Molecular mechanism of ligand binding and activation. Front Pharmacol 8: 898. doi: 10.3389/fphar.2017.00898
|
| [49] |
Liang YL, Khoshouei M, Radjainia M, et al. (2017) Phase-plate cryo-EM structure of a class B GPCR-G-protein complex. Nature 546: 118–123. doi: 10.1038/nature22327
|
| [50] |
Liang YL, Khoshouei M, Deganutti G, et al. (2018) Cryo-EM structure of the active, Gs-protein complexed, human CGRP receptor. Nature 561: 492–497. doi: 10.1038/s41586-018-0535-y
|
| [51] |
Zhang Y, Sun B, Feng D, et al. (2017) Cryo-EM structure of the activated GLP-1 receptor in complex with a G protein. Nature 546: 248–253. doi: 10.1038/nature22394
|
| [52] |
Koehl A, Hu H, Maeda S, et al. (2018) Structure of the micro-opioid receptor-Gi protein complex. Nature 558: 547–552. doi: 10.1038/s41586-018-0219-7
|
| [53] |
Draper-Joyce CJ, Khoshouei M, Thal DM, et al. (2018) Structure of the adenosine-bound human adenosine A1 receptor-Gi complex. Nature 558: 559–563. doi: 10.1038/s41586-018-0236-6
|
| [54] |
Kang Y, Kuybeda O, de Waal PW, et al. (2018) Cryo-EM structure of human rhodopsin bound to an inhibitory G protein. Nature 558: 553–558. doi: 10.1038/s41586-018-0215-y
|
| [55] | Khoshouei M, Radjainia M, Baumeister W, et al. (2017) Cryo-EM structure of haemoglobin at 3.2 A determined with the Volta phase plate. Nat Commun 8: 16099. |
| [56] |
Renaud JP, Chari A, Ciferri C, et al. (2018) Cryo-EM in drug discovery: Achievements, limitations and prospects. Nat Rev Drug Discov 17: 471–492. doi: 10.1038/nrd.2018.77
|
| [57] | Green SA, Holt BD, Liggett SB (1992) Beta 1- and beta 2-adrenergic receptors display subtype-selective coupling to Gs. Mol Pharmacol 41: 889–893. |
| [58] |
Niesen FH, Berglund H, Vedadi M (2007) The use of differential scanning fluorimetry to detect ligand interactions that promote protein stability. Nat Protoc 2: 2212–2221. doi: 10.1038/nprot.2007.321
|
| [59] |
Yen HY, Hoi KK, Liko I, et al. (2018) PtdIns(4,5)P2 stabilizes active states of GPCRs and enhances selectivity of G-protein coupling. Nature 559: 423–427. doi: 10.1038/s41586-018-0325-6
|
| [60] |
Gales C, Rebois RV, Hogue M, et al. (2005) Real-time monitoring of receptor and G-protein interactions in living cells. Nat Methods 2: 177–184. doi: 10.1038/nmeth743
|
| [61] |
Hein P, Frank M, Hoffmann C, et al. (2005) Dynamics of receptor/G protein coupling in living cells. EMBO J 24: 4106–4114. doi: 10.1038/sj.emboj.7600870
|
| [62] |
Dixon AS, Schwinn MK, Hall MP, et al. (2016) NanoLuc complementation reporter optimized for accurate measurement of protein interactions in cells. ACS Chem Biol 11: 400–408. doi: 10.1021/acschembio.5b00753
|
| [63] |
Irannejad R, Pessino V, Mika D, et al. (2017) Functional selectivity of GPCR-directed drug action through location bias. Nat Chem Biol 13: 799–806. doi: 10.1038/nchembio.2389
|
| [64] |
Irannejad R, Tomshine JC, Tomshine JR, et al. (2013) Conformational biosensors reveal GPCR signalling from endosomes. Nature 495: 534–538. doi: 10.1038/nature12000
|
| [65] | Kenakin T, Christopoulos A (2013) Signalling bias in new drug discovery: Detection, quantification and therapeutic impact. Nat Rev Drug Discov 12: 205–216. |
| [66] |
Gentry PR, Sexton PM, Christopoulos A (2015) Novel allosteric modulators of G protein-coupled receptors. J Biol Chem 290: 19478–19488. doi: 10.1074/jbc.R115.662759
|
| [67] |
Skiba NP, Bae H, Hamm HE (1996) Mapping of effector binding sites of transducin alpha-subunit using G alpha t/G alpha i1 chimeras. J Biol Chem 271: 413–424. doi: 10.1074/jbc.271.1.413
|
| 1. | Bodo Herzog, Optimal sequencing to reform the European economic and Monetary Union: a roadmap, 2021, 5, 2573-0134, 311, 10.3934/QFE.2021014 |
Figures(3) / Tables(1)
Byron Carpenter. Current applications of mini G proteins to study the structure and function of G protein-coupled receptors[J]. AIMS Bioengineering, 2018, 5(4): 209-225. doi: 10.3934/bioeng.2018.4.209
DownLoad:
