- 1Center of Physiology, Institute of Vegetative Physiology, University of Cologne, Cologne, Germany
- 2Institute of Physiology, Brandenburg Medical School Theodor Fontane, Neuruppin, Germany
- 3Research Cluster, Molecular Mechanisms of Cardiovascular Diseases, Brandenburg Medical School Theodor Fontane, Neuruppin, Germany
- 4Center of Physiology, Institute of Neurophysiology, University of Cologne, Cologne, Germany
- 5Institute of Anatomy, University of Cologne, Cologne, Germany
- 6Institute of Anatomy, Brandenburg Medical School Theodor Fontane, Neuruppin, Germany
- 7Department of Biochemistry and Molecular Biology, Brody School of Medicine at East Carolina University, Greenville, NC, United States
- 8Biological Sciences, Southampton General Hospital, University of Southampton, Southampton, United Kingdom
Stretch-induced vascular tone is an important element of autoregulatory adaptation of cerebral vasculature to maintain cerebral flow constant despite changes in perfusion pressure. Little is known as to the regulation of tone in senescent basilar arteries. We tested the hypothesis, that thin filament mechanisms in addition to smooth muscle myosin-II regulatory-light-chain-(MLC20)-phosphorylation and non-muscle-myosin-II, contribute to regulation of stretch-induced tone. In young BAs (y-BAs) mechanical stretch does not lead to spontaneous tone generation. Stretch-induced tone in y-BAs appeared only after inhibition of NO-release by L-NAME and was fully prevented by treatment with 3 μmol/L RhoA-kinase (ROK) inhibitor Y27632. L-NAME-induced tone was reduced in y-BAs from heterozygous mice carrying a point mutation of the targeting-subunit of the myosin phosphatase, MYPT1 at threonine696 (MYPT1-T696A/+). In y-BAs, MYPT1-T696A-mutation also blunted the ability of L-NAME to increase MLC20-phosphorylation. In contrast, senescent BAs (s-BAs; >24 months) developed stable spontaneous stretch-induced tone and pharmacological inhibition of NO-release by L-NAME led to an additive effect. In s-BAs the MYPT1-T696A mutation also blunted MLC20-phosphorylation, but did not prevent development of stretch-induced tone. In s-BAs from both lines, Y27632 completely abolished stretch- and L-NAME-induced tone. In s-BAs phosphorylation of non-muscle-myosin-S1943 and PAK1-T423, shown to be down-stream effectors of ROK was also reduced by Y27632 treatment. Stretch- and L-NAME tone were inhibited by inhibition of non-muscle myosin (NM-myosin) by blebbistatin. We also tested whether the substrate of PAK1 the thin-filament associated protein, caldesmon is involved in the regulation of stretch-induced tone in advanced age. BAs obtained from heterozygotes Cald1+/− mice generated stretch-induced tone already at an age of 20–21 months old BAs (o-BA). The magnitude of stretch-induced tone in Cald1+/− o-BAs was similar to that in s-BA. In addition, truncation of caldesmon myosin binding Exon2 (CaD-▵Ex2−/−) did not accelerate stretch-induced tone. Our study indicates that in senescent cerebral vessels, mechanisms distinct from MLC20 phosphorylation contribute to regulation of tone in the absence of a contractile agonist. While in y-and o-BA the canonical pathways, i.e., inhibition of MLCP by ROK and increase in pMLC20, predominate, tone regulation in senescence involves ROK regulated mechanisms, involving non-muscle-myosin and thin filament linked mechanisms involving caldesmon.
Introduction
We have recently shown that aging of cerebral arteries is associated with a hypercontractile state of smooth muscle (Lubomirov et al., 2017). The hypercontractile state was blunted by NO released from the endothelium and was associated with increased inhibitory phosphorylation of MYPT1 and increased filamentous F-actin protein content. However, the mechanism leading to this hypercontractile response was undetermined as there are many potential sites at which aging might act.
Vascular smooth muscle tone is primarily the result of the cyclic interaction of thick myosin filaments with filamentous actin (F-actin) in a cycle driven by ATP hydrolysis (Arner and Pfitzer, 1999). Although smooth muscle myosin II is the primary motor protein of smooth muscle contraction, there is increasing evidence that non-muscle myosin II is involved in slowly developing contractions (Morano et al., 2000) and tension maintenance ((Zhang and Gunst, 2017) reviewed in (Brozovich et al., 2016)). Regulation of smooth muscle contraction occurs by altering the level of phosphorylation of myosin, by altering the extent of actin polymerization into thin filaments, and by proteins that bind directly to actin, such as caldesmon (Morgan and Gangopadhyay, 2001). All of these points of regulation are potentially altered by aging leading to the hypercontractile state.
Both myosin types, smooth muscle and non-muscle myosin II, are activated by phosphorylation of regulatory light chains (MLC20) on serine19 (S19) (Zhang and Gunst, 2017). Phosphorylation is catalyzed by Ca2+-calmodulin dependent myosin light chain kinase (MLCK) in response to an increase in cytosolic [Ca2+]. The rise in Ca2+ is triggered by various extracellular stimuli such as neurotransmitters, metabolites and mechanical forces (Brozovich et al., 2016). MLC20 is phosphorylated additionally by Ca2+-independent, non-canonical MLC-kinases, which include RhoA-kinase (ROK) (Amano et al., 1996).
Dephosphorylation of MLC20 is catalyzed by a constitutively active type 1 phosphatase (MLCP) the activity of which is modulated by phosphorylation of its myosin targeting (regulatory) subunit MYPT1 at several serine and threonine residues (Matsumura and Hartshorne, 2008). This allows force to increase or decrease without changes in cytosolic [Ca2+], a phenomenon known as Ca2+-sensitization and desensitization, respectively [reviewed in (Somlyo and Somlyo, 2003)].
The RhoA-ROK kinase pathway is the key mechanism of increasing Ca2+-sensitivity by phosphorylation of MYPT1 at threonine-696 and threonine-853 (Matsumura and Hartshorne, 2008; Puetz et al., 2009). The nitric oxide (NO) protein kinase G (PKG) signaling cascade decreases Ca2+-sensitivity by antagonizing Rho-PKG signaling, which disinhibits MLCP by dephosphorylation of MYPT1 (Puetz et al., 2009) and/or activation of MLCP by phosphorylation of MYPT1 at the PKG site serine-668 (Yuen et al., 2011).
Force is also diminished when actin filaments are disassembled or depolymerized. RhoA and ROK promote adhesome complex formation and polymerization of a small G-actin pool at the cell cortex upon contractile stimulation (Zhang and Gunst, 2019). Large actin filaments are anchored to the cell membrane via adhesion proteins, which contain some cortical actin polymers. They allow force to be transmitted from myofilaments to the extracellular matrix and act independently of pathways that activate cross-bridge cycling (Gunst and Zhang, 2008). Interfering with actin polymerization results in diminished active force (Zhang et al., 2018). Cortical actin polymerization induced by biomechanical stimulation and contractile agonists requires recruitment of inactive adhesome proteins to membrane adhesion junctions, where they are activated (Gunst and Zhang, 2008; Zhang et al., 2015). The underlying complex signaling mechanisms have been delineated in detail in airway smooth muscle [reviewed in (Gunst and Zhang, 2008)]. They involve the small GTPase, RhoA upstream of non-muscle myosin II (Zhang and Gunst, 2017), and ROK upstream of p21-activated protein kinase (PAK) (Zhang et al., 2018), and eventual activation of the N-WASp and Arp2/3 complex.
A growing body of evidence supports the notion that F-actin dynamics is also involved in regulation of vascular smooth muscle tone (Yamin and Morgan, 2012) and specifically in the myogenic responses of rat brain vasculature (Cipolla et al., 2002; Moreno-Dominguez et al., 2013; Walsh and Cole, 2013). Moreover, an increased F-actin content was associated with a hypercontractile state of old murine brain arteries (Lubomirov et al., 2017).
PAK1 is not only an upstream activator of actin polymerization (Zhang et al., 2018), but also a stimulator of force production by phosphorylating caldesmon and reversing its inhibitory effect on actomyosin interaction (Van Eyk et al., 1998). Caldesmon is a myosin-, actin, calmodulin-, and tropomyosin-binding protein, and is expressed in two isoforms generated by alternative splicing [(Ruegg and Pfitzer, 1991; Gusev, 2001) rev. (Pfitzer et al., 1993a)]. The smooth muscle specific isoform, h-CaD, is found in the actomyosin domain of smooth muscle cells, where it presumably tethers actin- and myosin filaments and stabilizes the contractile machinery (Pfitzer et al., 1993b; Smolock et al., 2009; Putz et al., 2021). Both isoforms inhibit actomyosin MgATPase activity and contraction at constant MLC20 phosphorylation. CaD may act as a molecular brake of smooth muscle contraction especially at low MLC20 phosphorylation and cytosolic [Ca2+] (Smolock et al., 2009; Putz et al., 2021). Phosphorylation of caldesmon at several threonine and serine residues by various protein kinases, including PAK1, facilitates activation of MgATPase activity and cross bridge cycling. PAK1 is activated by ROK-dependent phosphorylation (Zhang et al., 2018), giving rise to the interesting hypothesis that ROK may relieve the normal inhibition of tension development caused by caldesmon.
To conclude there is evidence primarily derived from smooth muscle tissue from young animals that contractile activity is regulated in addition to the canonical pathway via phosphorylation of MLC20 by, actin based mechanisms, i.e., actin filament dynamics and/or the thin filament linked protein caldesmon. In these pathways, ROK is center stage and regulates different signaling cascades that target these mechanisms. Little is currently known about the contribution of these pathways in brain arteries from senescent mice.
In the present work we explored whether changes in basal phosphorylation levels of the Rho-kinase effectors MLCP, non-muscle myosin, as well as actin filament dynamics parallel the augmented in stretch-induced tone in old basilar arteries. Further, we investigated whether ablation of the PAK1 effector, caldesmon augments stretch induced tone. Specifically we tested the following hypotheses: 1) stretch induced tone is augmented in senescent basilar arteries in a ROK-dependent manner, 2) ROK is upregulated in senescence and serves as a central signaling hub, which not only increases inhibitory phosphorylation of MYPT1 but also of non-muscle myosin and PAK1 3) the PAK1-target caldesmon acts as a molecular brake to attenuate stretch induced tension rise.
Materials and methods
Ethic statement and detail description of all materials and methods can be found in the online supplement.
Results
Three different types of muscle tone were examined in this work: 1) stretch-induced tone refers to tone that is present after stretching y-BA to IC90 (90% of the internal circumference at wall tension corresponding to a transmural pressure of 100 mmHg (IC90); see expandet “Methods” section in online supplement) and in s-BA it refers to the slow tension rise that develops in 20–25 min after stretching the preparations to IC90 in the absence of treatment. This tone is maintained throughout the experiment. 2) L-NAME-induced tone develops in response to addition of L-NAME, a pan-NOS inhibitor; 3) agonist-induced tone occurs after stimulation with contractile agonist, U46619. Our first goal was to determine the effect of ROK on each of these types of muscle tone.
ROK inhibition blunts stretch-, L-NAME-, and agonist-induced tone in young murine basilar arteries
Tone in the absence of L-NAME amounted to ∼1.7 mN, present after the normalization procedure. To demonstrate that this increase in tone involves ROK, ring preparations from BAs were treated after the normalization procedure (20–25 min equilibration in PSS see “Methods”) with 3 μmol/L Y27632 for 10 min, followed by addition of 100 μmol/L of the pan-NOS inhibitor L-NAME for another 20 min in the presence of Y27632. Pretreatment with Y27632 did not alter stretch-induced tone (1.69 ± 0.1 mN in the presence of Y27632 versus 1.73 ± 0.28 in non-treated time matched controls), but prevented the increase in L-NAME induced tone, which was 1.66 ± 0.1 mN in the presence of Y27632 and 2.7 ± 0.5 mN in the absence of Y27632 (Figures 1A, C).
FIGURE 1. Stretch-, L-NAME-, and agonist-induced force in young murine basilar arteries obtained from WT and MYPT1-T696A/+ -mice. (A–C) Original force tracings and statistical summary from measurement of stretch-induced tone stretch-induced tone in y-BAs from WT-animals treated with 100 μmol/L of pan-NOS-inhibitor L-NAME in the presence of 3 μmol/L Y27632 (gray) or vehicle (H2O; black). Depicted p-values represent results from 5 independent experiments (***p < 0.001; n.s.—not significant, p > 0.05; one way ANOVA n = 5). (D–F) Original force tracings and statistical evaluation of stretch-induced tone, L-NAME- and agonist-induced tone of y-BAs from WT (black) and MYPT1-T696A/+ -mice (red) under control conditions (controls) and after application of 100 μmol/l L-NAME (n = 6–7; *p < 0.05; unpaired t-test; n.s.—not significant). Bars represent mean ± SEM.
We then tested whether ROK inhibition also attenuated the agonist-induced tone by the thromboxaneA2 analogue, U46619. In the presence of Y27632, the cumulative concentration-response relation was shifted to the right; the pEC50 value was significantly lower with than without Y27632 in time-matched controls (6.3 ± 0.1 logarithmic units vs. 7.1 ± 0.18; p < 0.01; n = 5). Maximal force was also lower (Fmax 3.68 ± 0.37 mN with vs. 5.89 ± 1.3 mN without Y27632, p < 0.05; t-test; n = 5; Figure 1B). However, Y27632 at a concentration of 3 μmol/L (concentration able to nearly completely dephosphorylate ROK-site of targeting subunit of MLCP, MYPT1-T853; Figures 5A, D; Figures 8C, D) attenuated U46619 contractile response by only ∼40% compared to complete inhibition of L-NAME induced tone, suggesting that a part of the U46619-induced tone involves one or more additional mechanisms, distinct from ROK activation.
Effect of blocking phosphorylation at threonine 696 of MYPT1 on agonist-induced tone in y-BA
The inhibitory effect of the ROK inhibitor on L-NAME and U46619 induced tone may partially result from reduced phosphorylation of MYPT1-T696 (Puetz et al., 2009), a site that is phosphorylated by several protein kinases including ROK. Phosphorylation of this site inhibits MLCP activity in vitro (Feng et al., 1999). To test the involvement of phosphorylation this site in L-NAME and U46619 induced contractile activity, we used BAs from 2 months old wild type and heterozygous MYPT1-T696A/+ mice, in which the phosphorylation propensity of MYPT1 at T696 was genetically lowered by mutating the threonine-residue into a non-phosphorylatable alanine (MYPT1-T696A/+).
Ring preparations from MYPT1-T696A/+ and WT BA were stimulated after mounting and equilibration in PSS by cumulatively increasing concentrations of U46619. U46619 induced a concentration dependent rise in tone with similar Fmax and pEC50 values in mutant and WT littermates (p > 0.05; n = 7–6; Figures 1D–E). After wash-out of the agonist with PSS, the vessels were pre-incubated for 20 min in L-NAME (100 μmol/L), followed by a second U46619 concentration-response relation (Figures 1D, E). In time-matched controls, the second concentration-response relation was obtained in the absence of L-NAME (n = 4–5; Supplementary Figure S2). L-NAME shifted the U46619-concentration-response curve to the left in both WT and mutant BA and increased Fmax, whereby neither Fmax nor the pEC50 values differed between MYPT1-T696A/+ and WT littermates (p > 0.05; n = 7–6; Figures 1D, E). No significant difference between the first and second U46619 concentration response curve was observed in time matched controls (n = 5–4; Supplementary Figure S2). L-NAME-induced tone was significantly less in mutant BA than in WT (mutant: 1.95 ± 0.1 mN vs. WT 2.82 ± 0.2 mN, n = 7–6; Figures 1D, F). Thus, the effect of the mutation of MYPT1 and Y27632 on L-NAME induced tone were of similar magnitude. All contractile parameters of BAs from MYPT1-T696A/+ and WT mice are given in Table 1.
TABLE 1. Contractile parameters of intact preparations from young basilar arteries from wild type (WT) and heterozygous (MYPT1-T696A/+) mice.
L-NAME increased tone and MLC20 phosphorylation in WT but not in MYPT1-T696A/+ BA from young mice
Next, we tested in separate experiments whether the reduction in tone induced by L-NAME in y-BAs from MYPT1-T696A/+ mice was related to lower phosphorylation levels of the regulatory light chain of myosin (MLC20) under basal conditions. For technical reasons, the vessels were mounted on two wires (ø 25 µm) but not stretched as for mechanical experiments. Basal phosphorylation levels of MLC20-S19 in WT and MYPT1-T696A/+ were similar (Figures 2A, B). In line with the attenuated force, incubation with L-NAME was associated with increased MLC20-S19 phosphorylation in BA of WT, but not in BA of MYPT1-T696A/+.
FIGURE 2. Basal and L-NAME-induced phosphorylation of pMLC20-S19, pMYPT1-T696, pMYPT1-T853, pMYPT1-S695 in basilar arteries from young WT and MYPT1-T696A/+ -mice. Original chemiluminograms (A) and statistical evaluation (B–E) of the immunoreactive signals of lysates of y-BAs from WT and MYPT1-T696A/+ mice transferred on nitrocellulose membranes and probed with antibodies against pMLC20-S19, pMYPT1-T696, pMYPT1-T853, and pMYPT1-S695. Statistical figures represent the ratio of the immunoreactive signal of pMLC20-S19 normalized to β-Actin or pMYPT1-T696, pMYPT1-T853, and pMYPT1-S695 normalized to MYPT1-total (n = 4). Arteries were treated using same experimental protocol as in Figure 1. Statistic comparison on (B) represents evaluation using two-way ANOVA (***p < 0.0001; n = 4). ***p < 0.0001, *p < 0.05 represents evaluation by unpaired t-test (n = 4). n.s.—not significant.
We further tested, whether observed differences in MLC20-S19 phosphorylation were caused by differences in the phosphorylation of MYPT1 at T696 and T853 and thus in activity of MLCP. Under control conditions, phosphorylation of MYPT1 at T696 (PSS, no treatment) was higher in BA from WT-mice than in those from MYPT1-T696A/+-mice (p < 0.05; unpaired t-test; n = 4; Figures 2A, C). Nominal maximal pMYPT1-T696 phosphorylation was assessed by stimulation with 0.1 μmol/L calyculin, a type1 phosphatase inhibitor (Chen et al., 2015). As expected, MYPT1-T696 phosphorylation in calyculin treated mutant BA amounted to only ∼50% of that in WT BA (Supplementary Figure S3). We did not observe a significant change in phosphorylation of pMYPT1-T853 between WT and mutant mice, under control conditions (PSS, no treatment) or after incubation with L-NAME (Figures 2A, D). This result suggests that this site, which is phosphorylated only by ROK, is not affected by silencing phosphorylation in MYPT1-T696.
To test whether altered PKG activity may contribute to the lower L-NAME induced tone, we determined phosphorylation of MYPT1-S695, which is specifically phosphorylated by PKG. L-NAME decreased phosphorylation of this site, indicated that NO-PKG signaling is constitutively active in this preparation. Phosphorylation of pMYPT1-S695 was similar in WT and mutant BA (Figures 2A, E). These results also indicated that the alanine mutation of T696 did not affect the phosphorylation of the neighboring S695, and further, that the NO-PKG pathway was not affected by the mutation.
Ca2+-sensitivity is lower in young, α-toxin permeabilized MYPT1-T696A/+ BA
The level of phosphorylation of MYPT1-T696 was proposed to define the intrinsic Ca2+-sensitivity of the contractile machinery of urinary bladder smooth muscle (Khromov et al., 2009). To test whether this applies to basilar arteries, we obtained Ca2+-force relations in α-toxin permeabilized preparations. Permeabilization with α-toxin generates pores into the plasma membrane which allow diffusion of molecules <1 kDa into the intracellular space. This permits intracellular [Ca2+] to be controlled by EGTA so that the Ca2+-responsiveness of the contractile machinery can be measured without confounding changes in cytosolic [Ca2+]. After permeabilization, the preparations were equilibrated for 20 min in pCa >8 relaxing solution and then challenged with cumulatively increasing concentrations of [Ca2+] (Figure 3). We found that force in MYPT1-T696A/+ BA was lower than in WT BA at intermediate, but not at maximal [Ca2+] (Figure 3; Table 2). The pEC50, the neg. logarithm of concentration of Ca2+ for half maximal contraction, and the Hill coefficient of the Ca2+-force relation were lower than in WT (Table 2). Thus, under Ca2+-clamped conditions, the Ca2+-sensitivity was lower in mutant compared to WT vessels, consistent with a reduced inhibitory phosphorylation of MLCP. This result indicates that a lower Ca2+-sensitivity contributes to the lower L-NAME-induced tone in mutant y-BA.
FIGURE 3. Ca2+-activated contraction of basilar arteries from young WT- and MYPT1-T696A/+-mice. Original force tracing (A) and statistical summary (B) depicting Ca2+-activated contraction of young BA from WT- and MYPT1-T696A/+-mice (n = 4–5). Tone normalized to pCa 4.3 accepted as 100%.
TABLE 2. Ca2+-activated force in α-toxin permeabilized basilar arteries from WT and MYPT1-T696A/+ mice.
Expression of MYPT1-T696A/+ does not alter stretch-induced tone- and L-NAME-induced tone in senescent BA
In contrast to y-BA, in senescent BA (s-BA, >24 months old) a slow spontaneous rise in force was observed (Figure 4A) after stretching them to IC90 and equilibration in PSS for 20 min during the normalization procedure (see Methods). This tone amounted to 3.5 ± 0.9 mN in s-BAs from WT mice, which was not reduced in s-BA MYPT1-T696A/+ mice (3.1 ± 1.3 mN; Figures 4A, C; n = 8). The preparations from both genotypes were then incubated with 100 μmol/L L-NAME for 20 min, followed by stimulation with increasing concentrations of U46619 (0.001–3 μmol/L, Figure 4A). As in the experiments with y-BA, treatment of s-BA with 100 μM L-NAME increased tone in both mouse lines to a similar extent (4.5 ± 0.5 mN in WT vs. 4.2 ± 0.5 mN in mutant; n = 8). However, unlike the case with y-BA, the mutation did not blunt L-NAME induced tone (Figures 4A, C). U46619 concentration-response curves were similar in s-BAs from both mouse lines (Figures 4A, C, D; Table 3).
FIGURE 4. Stretch-, L-NAME-, and U46619-induced tone in arterial rings from senescent basilar arteries from WT and MYPT1-T696A/+ -mice under control conditions and under inhibition of ROK. (A,B) Original force tracings of the experiments performed with s-BAs from WT and MYPT1-T696A/+ -mice pretreated with vehicle [1% H2O; (A)] or 3 μmol/L Y27632 (B) and followed by inhibition of endogenous NO by treatment with 100 μmol/l L-NAME and cumulative application of U46619 (conc. 0.001–3 μmol/L). (C) Statistical evaluation of stretch-induced, L-NAME-induced and maximal tone (n = 8). (D) Statistical evaluation of tone induced by cumulative application of U46619 (n = 8). Data are represented as absolute force ±SEM; n.s. p > 0.05; two-way ANOVA. §p < 0.0001, stretch-induced tone WT controls (no treatment) vs. stretch-induced tone WT after 10 min treatment with 3 μmol/L Y27632, or stretch-induced tone L-NAME WT no treatment vs. stretch-induced tone L-NAME WT after 10 min treatment with 3 μmol/L Y27632, or stretch-induced tone L-NAME MYPT1-T696A/+ no treatment vs. stretch-induced tone L-NAME L-NAME MYPT1-T696A/+ after 10 min treatment with 3 μmol/L Y27632. *p < 0.05 Fmax WT vs. Fmax WT after treatment with Y27632; and **p < 0.001 stretch-induced tone MYPT1-T696A/+ controls (no treatment) vs. stretch-induced tone MYPT1-T696A/+ after 10 min treatment with 3 μmol/L Y27632 or Fmax MYPT1-T696A/+ vs. Fmax MYPT1-T696A/+ after treatment with Y27632. All values were calculated by unpaired t-test. n.s.—not significant; pEC50; ***p < 0.0001; two-way ANOVA.
TABLE 3. Contractile parameters of intact preparations from senescent basilar arteries from wild type (WT) and heterozygous (MYPT1-T696A/+) mice.
In y-BAs, inhibition of ROK blunted completely L-NAME-induced tone. We therefore investigated the effect of Y27632 in WT and MYPT1-T696A/+ s-BAs. U46619 concentration relationships and Fmax after Y27632 treatment were similar between s-BAs from both groups (Figures 4B–D).
Y27632 reduced phosphorylation of pMYPT1-T696/T853 and pMLC20-S19 in s-BAs
In s-BAs from MYPT1-T696A/+ mice, the basal pMLC20-S19 and pMYPT1-T696 were lower than in WT arteries (Figures 5A–C). Incubation with 3 μmol/L of Y27632 reduced pMLC20-S19 only in WT s-BAs (Figures 5A, B). pMYPT1-T853 in WT and mutant s-BAs did not differ under basal conditions or L-NAME-, or Y27631-treatment (Figures 5A, D).
FIGURE 5. Effect of ROK inhibition on MLC20-S19, MYPT1-T696, and MYPT1-T853 in senescent basilar arteries from WT and MYPT1-T696A/+ mice. (A) Original western blots and statistical evaluation of the phosphorylation of (B) MLC20-S19 (pMLC20), and MYPT1 at (C) T696 (pMYPT1-T696) and (D) T853 (pMYPT1-T853) in s-BAs from WT and MYPT1-T696A/+ mice under control conditions (no stimulations), and after preincubation with 3 μmol/L Y27632, 100 μmol/l L-NAME or Y27632 plus L-NAME. Data are represented as ratio of pMYPT1-T696/853 toward MYPT1-total or MLC20-S19 toward GAPDH (n = 6). **p < 0.01, *p < 0.05, n.s.—not significant; two-way ANOVA. n.s.—not significant; ***p < 0.0001 and p = 0.06 have been calculated by using unpaired t-test.
Non-muscle myosin II is involved in stretch-induced tone and L-NAME-induced tone
Stretch-induced tone and L-NAME induced tone development under resting conditions, i.e., in PSS, is much slower than agonist induced contraction. Such slow force development has been observed in neonatal mice, in which smooth muscle myosin II was genetically ablated and force development was ascribed to non-muscle myosin (NM-II) (Morano et al., 2000). Blebbistatin has been used as a small molecule inhibitor of NM-II (Limouze et al., 2004; Rhee et al., 2006), but its specificity has recently been questioned as it was reported to also inhibit smooth muscle myosin II (Eddinger et al., 2007). The active (−) form of blebbistatin reduced stretch-induced tone and L-NAME induced tone in s-BAs to a similar degree as Y27632 (Figures 6A–C; Table 3). Under (+) blebbistatin treatment stretch-induced tone was 2.0 ± 0.3 mN vs. 2.2 ± 0.3 mN in heterozygous s-BAs from MYPT1-T696A/+ animals and amounted to 3.1 ± 0.3 mN and 3.0 ± 0.2 mN after 20 min treatment with L-NAME (n = 5; Figure 6A; and Table 3). In 2 of 5 experiments application of 30 µmol/L (+) Blebbistatin had no additional effect on tone. The concentration responsiveness and Fmax of U46619 in the presence of blebbistatin (−) or (+) were similar to those measured previously in Y27632 treated group (Figures 6A, C, D; Table 3).
FIGURE 6. Stretch-, L-NAME-, and U46619-induced tone in senescent basilar arteries from WT and MYPT1-T696A/+ -mice under inhibition of cross-bridge cycling of non-muscle myosin. (A,B) Original force tracings representing the effect of the inhibitor of cross-bridge cycling of non-muscle myosin, blebbistatin on tone. Vessels were pretreated with the active (−) and inactive (+) enantiomers of blebbistatin and further incubated with 100 μmol/l L-NAME and cumulative application of U46619. (C) Statistic evaluation of stretch-induced, L-NAME-induced and maximal tone in presence of (+) or (−) blebbistatin (n = 5). (D) Statistic evaluation of tone induced by cumulative application of U46619 (n = 5). Data represented as absolute force ±SEM; n.s. p > 0.05; two-way ANOVA. *p < 0.05 and **p < 0.001 calculated by unpaired t-test.
Blebbistatin (−) unlike Y27632 did not reduce phosphorylation of pMYPT1-T696/T853 and pMLC20-S19 in s-BAs
Next we tested whether blebbistatin reduced basal MLC20-S19 phosphorylation or inhibitory phosphorylation of MYPT1. Phosphorylation of pMLC20-S19, pMYPT1-T696, and pMYPT1-T853 in the presence of blebbistatin (−) was determined in s-BAs from WT and from MYPT1-T696A/+ mice using phosphospecific antibodies. Neither pMLC20-S19, nor pMYPT1-T696/T853, were reduced by blebbistatin (−) in BAs from both mouse lines (Figures 7A–D). These experiments support the notion that the hyper-contractile phenotype of senescent BAs involves increase in ROK activity and likely non-muscle myosin-II cross-bridge cycling.
FIGURE 7. Phosphorylation of pMYPT1-T696, pMYPT1-T853, and pMLC20-S19 in senescent basilar arteries from WT and MYPT1-T696A/+ mice under inhibition of cross-bridge cycling of non-muscle myosin. (A–D) Original Western blots (A) and statistic evaluation of the phosphorylation of pMLC20-S19 (B), pMYPT1-T696 (C) and pMYPT1-T853 (D) in s-BAs from WT and MYPT1-T696A/+ mice under control conditions and after inhibition of cross-bridge cycling of non-muscle myosin by blebbistatin (−). Preparations mounted as previously and treated by vehicle (0.3% DMSO; controls), or 10 μmol/L blebbistatin (−), or 100 μmol/l L-NAME, or blebbistatin (−) plus L-NAME. Data represented as ratio of pMYPT1-T696/853 toward MYPT1-total or MLC20-S19 toward GAPDH (n = 5). n.s.—not significant; unpaired t-test.
Phosphorylation of non-muscle myosin-II in young and senescent BAs
Contractile activity requires that NM myosin-II proteins assemble into filaments. A predictor of filament competent NM myosin is the level of phosphorylation of its heavy chains (MHC) at S1943 (Zhang and Gunst, 2017). Here, we tested 1) whether NM-II S1943 phosphorylation was higher in senescent than in young arteries, 2) whether Y27632 decreased phosphorylation of this site and 3) whether the expression pattern of global or IIa and IIb isoforms of NM-II are altered by senescence. Phosphorylation of NM-MHC-S1943 was quantified by Western blot analysis using phosphospecific antibodies. In s-BAs, basal S1943 phosphorylation was significantly higher under control conditions (time matched controls to the Y-27632 treated preparations) than in y-BAs (Figures 8A, B). Incubation with 3 μmol/L Y27632 reduced the immunoreactivity against pNM-II-S1943 in s-BAs, but not in y-BAs (Figures 8A, B). No difference in the expression levels of NM-II-total between y- and s-BAs was observed (Supplementary Figures S5A, B). Interestingly, using the same scan intensity, the intensity of the immunoreactive signal obtained with NM-IIa antibodies was larger than with NM-IIb antibodies, supporting the view that NM-IIa is the predominant isoform in BAs of both ages (Supplementary Figures S5A, C, D).
FIGURE 8. Phosphorylation of NM-II-S1943, MYPT1-T853, and PAK-T423 in basilar arteries from young and senescent mice under control conditions and after inhibition of ROK. (A–C) Original western blots and statistic evaluation of the phosphorylation of MYPT1 at T853 (A,B), non-muscle-myosin II (NM-Myosin) at S1943 (C,D) and PAK-T423 (E,F) in y- and s-BAs from WT animals (n = 10–6). Preparations were treated either by vehicle (0.5% H2O; controls) or by 3 μmol/L Y27632, shock-frozen and subjected to Western blot as described in “Methods.” **p < 0.01 and *p < 0.05; control y-BAs vs. control s-BAs; two-way ANOVA. **p < 0.001 and *p < 0.05; unpaired t-test.
In these experiments, phosphorylation of MYPT1-T853 as a marker for ROK-activity was monitored in parallel. Phosphorylation of T853 was significantly higher in s-BAs than in y-BAs. In both age groups, Y27632 strongly reduced the immunoreactivity of pMYPT1-T853 (Figures 8 C, D). These experiments indicate that ROK regulates NM-II S1943 phosphorylation.
Senescence increased basal phosphorylation of PAK1-T423, a down-stream effector of ROK
PAK1 ist another target of ROK, which phosphorylates PAK1 at threonine 432 and thereby activates the enzyme (Zhang et al., 2018). We assessed the phosphorylation level of PAK1 using phosphospecific antibodies against pPAK1-T432 in Western Blots prepared with the same samples we used for analysis of NM-II-S1943- and MYPT1-T853-phosphorylation. The immunoreactive signal of PAK1-T423 in s-BAs was higher than in y-BAs. In both groups, incubation with Y27632 reduced the PAK1-T423 phospho-signal (Figures 8E, F).
Early onset of stretch-induced tone in Cald1+/− old basilar arteries (o-BAs)
PAK1 has several downstream targets involved in regulation of tone. One of them, caldesmon was reported to elicit a contraction without an increase in MLC20 phosphorylation when phosphorylated by PAK1 (Van Eyk et al., 1998). Thus, phosphorylation of caldesmon by PAK1 and reversal of its inhibitory effect on cross-bridge cycling is a potential mechanism that could account for the stretch-induced tone in s-BAs. Unfortunately, there are no commercial phosphospecific antibodies available to test whether caldesmon was phosphorylated by PAK1. Therefore, we investigated the involvement of caldesmon in the regulation of stretch-induced tone in arteries from heterozygous mice, carrying one allele in which Cald1 gene was deleted by homologous recombination of the Cald1 gene (Putz et al., 2021). As homozygosity is prenatally lethal, only heterozygous mice (Cald1+/−) were employed (Putz et al., 2021). In heterozygous old basilar arteries (20–21 months old; o-BAs) expression of caldesmon was ∼50% lower than in WT mice (n = 4, Figures 9A, B). Interestingly, stretch-induced tone was observed in o-BAs obtained from heterozygotes Cald1+/− mice at an earlier age of 20–21 months, compared to the BAs from controls MYPT1-T696A/+-mice in which such a tone appeared at senescence (age >24 months). At this age, o-BAs from WT mice did not exhibit stretch-induced tone yet (Figures 9C, D). The magnitude of stretch-induced tone of o-BAs from Cald1+/− mice was similar to that in WT s-BAs (3.5 ± 0.9 mN in s-BAs from WTs vs. 3.5 ± 0.6 mN in o-BAs from Cald1+/−; n = 8). There was no difference between both genotypes considering L-NAME- and U46619-induced tone (Figures 9C, E). As in previous experiments, application of 3 μmol/L Y27632 reduced stretch-induced tone in both groups (Figures 9C–F). We did not observe significant differences in stretch-induced or L-NAME-induced tone, nor in maximal -induced tone between y-BAs from WT and Cald1+/− mice (12–14 weeks old) (Supplementary Figure S4).
FIGURE 9. Effect of Caldesmon targeting on tone maintenance in basilar arteries from old mice. (A) Representative Western blot of expression of h-CaD in WT and heterozygous s-BAs from 2 litters (B) statistical evaluation from arteries of four independent measurements (n = 3 litters with WT and Het littermates, n = 1 WT and Het from different litters). M, marker; W, wilde type; H, heterozygotes. Caldesmon expression the Ponceau Red stained actin band. (B) Statistic evaluation. (C,D) Original force tracings representing the effect of caldesmon mutation on tone. (E) Statistic evaluation of stretch-induced, L-NAME-induced and maximal tone in BAs from WT and Cald1+/− BAs (n = 6–4). (F) Statistic evaluation of tone induced by cumulative application of U46619 (n = 6–4). Data represented as absolute force ±SEM; *p < 0.05; n.s. p > 0.05; two-way ANOVA.
Inhibition of cross-bridge cycling by caldesmon is exerted by the C-terminal, actin-binding domain (Pfitzer et al., 1993b). Binding to myosin was proposed to position caldesmon in such a way that the actin-binding domain can inhibit cross-bridge cycling (Lee et al., 2000). We tested this in a second mouse line, CaD-ΔEx2−/−, in which expression of the strong myosin binding domain, which is encoded by Exon 2 of the Cald1 gene, was ablated (Figure 10) (Pfitzer et al., 2005). Expression of the truncated CaD did not differ from WT (Figure 10A). Whereas WT CaD bound to myosin in solution, binding was nearly abolished in truncated caldesmon (CaD-ΔEx2−/−) (Figure 10B). Binding to actin was not affected (Figure 10B). o-BAs from homozygous mice expressing the truncated caldesmon were subjected to the same protocols used in previous experiments. No significant change of vascular tone in o-BAs from CaD-ΔEx2−/− mice compared to their time matched WT controls was observed (Figures 10C–F).
FIGURE 10. Effect of ablation of exon2 (Ex2) of Caldesmon on tone maintenance basilar arteries (A) Original luminogram and statistic evaluation showing caldesmon expression in s-BAs obtained from wild type animals, homozygous animals lacking the strong myosin-binding domain of caldesmon encoded by Exon 2 (−/−) and heterozygous animals (+/−) (n = 3-3-2). (B) Analysis of two independent experiments confirming the lack of myosin binding of the truncated caldesmon isolated from Cald1 −/− mice (left panel), whereas actin-binding was not affected (right panel). (C,D) Original force tracings representing the effect of ablation of Exon2 in caldesmon protein on tone. (E) Left: Statistic evaluation of stretch-induced (n = 7), L-NAME-induced (n = 3–7), and maximal tone (n = 3–7) in BAs from WT (CaD_ΔEx2 +/+) and CaD ΔEx2−/− BAs under control conditions. Right: Statistic evaluation of stretch-induced (n = 4–3), L-NAME-induced (n = 4–3), and maximal tone (n = 4–3) in BAs from same groups in the presence of 3 μmol/L Y27632. (F) Statistic evaluation of tone induced by cumulative application of U46619 (n = 6–4). Data represented as absolute force ±SEM.
Basal G/F-actin ratio and localization of F-actin in senescent BAs after maximal stimulation
The G/F-actin decreased in a ROK-dependent manner by agonist stimulation of airway smooth muscle cells (Zhang and Gunst, 2017) or under basal conditions in old BAs (Lubomirov et al., 2017). To validate that targeting of MYPT1-T696 does not influence F-actin polymerization, we next measured the G/F-actin ratio in s-BAs from WT and MYPT1-T696A/+ animals. No change of G-actin immunoreactivity in supernatants, used as read out for the G-actin fractions, was detected in BAs from WT compared to MYPT1-T696A/+ animals (Figures 11A, B).
FIGURE 11. G/F actin ratio in basilar arteries from senescent WT and MYPT1-T696A/+ animals under basal conditions. Localization of F-actin in basilar arteries from senescent WT mice after Calyculin stimulation (A) Original Western blot representing the global actin-immunoreactivity in supernatant (S) used as read out for globular (G) -actin fraction and pellet (P) used as read out for fibrillar (F) -actin in s-BAs from WT and MYPT1-T696A/+ animals. Supernatant and pellet fractions were obtained by ultracentrifugation as described in methods. (B) Statistic evaluation of n = 6 s-BAs, each group. n.s.—not significant; p > 0.05; two-way ANOVA. Confocal images y-BAs (C) and s-BAs (D) arteries from WT animals stained with phalloidin for F-actin and Hoechst for Nuclei. Transmission light images are denoted as “T-light.” Vessels were isolated, mounted and normalized as described in methods and treated for 30 min with 0.1 μmol/L Calyculin. After stimulation, preparations were fixed and F-actin was stained with Fluor™ 555 conjugated phalloidin, nuclei were stained with Hoechst 33342 (see methods). Representative images from four independent vessels/animals.
In the next series of experiments, we investigated whether maximal stimulation with calyculin resulted in a redistribution of F-actin in favor of subcortical actin. Calyculin treatment induced a maximal increase in force, starting 7–10 min after application and reached its maximum after 30 min (not shown). In the previous experiments, we showed that in BAs, calyculin treatment maximally activates ROK as seen by phosphorylation of MYPT1-T853 (Figure 7A). Alexa Fluor™ 555 phalloidin staining of F-actin filaments showed changes in localization in endothelial and smooth muscle cells. In s-BA endothelial cells F-actin localized to the cortex, while F-actin, stained by Alexa Fluor™ 555 phalloidin, in endothelium of y-BAs was diffuse and could not be attributed to a subcellular niche (Figures 11C, D). Interestingly, in smooth muscle cells from y-BAs F-actin localized predominantly to the cell-poles (Figure 11C).
Discussion
Stretch-induced vascular tone is an important element of autoregulatory adaptation of the cerebral circulatory system. It maintains constant cerebral flow, despite changes in perfusion pressure, and thus ensures adequate brain perfusion. The present work provides novel insights into the mechanisms of regulation of the stretch-induced tone of murine BA in senescence. 1) In young BAs (y-BAs) mechanical stretch did not lead to spontaneous tone generation, while the senescent-BAs (s-BAs) developed stable stretch-induced tone. 2) Stretch-induced tone in y-BAs appeared only after inhibition of NO-release by L-NAME, whereas in s-BAs such an inhibition led to an additive effect. 3) In BAs from both age groups, stretch-induced tone was fully inhibited by treatment with a ROK-inhibitor, Y27632. 4) Lowering of the phosphorylation propensity of MYPT1 in heterozygous mice carrying the point mutation MYPT1-T696A/+ prevented stretch-induced tone in y-BAs but not in aged BAs. Basal MLC20 phosphorylation was lower in MYTP1-T696A/+ than in WT BAs in both age groups. Based on this, we propose that regulation of stretch-induced tone in senescent arteries involves alternative mechanism. 5) The basal phosphorylation level of NM-II and PAK1, which phosphorylates and inhibits caldesmon, were higher in s-BAs than in y-BAs, and were reduced by the ROK inhibitor, Y27632. 6) Ablation of caldesmon augmented stretch-induced tone in old BAs suggesting that it acts as a molecular brake on stretch-induced tone in aged BAs. We propose that in senescent cerebral vessels, mechanisms distinct from MLC20 phosphorylation contribute to regulation of tone in the absence of a contractile agonist.
The onset of senescence is associated with MLC20-independent mechanisms of cerebrovascular tone regulation
There is ample evidence that stretch-induced cerebrovascular tone in young arteries relates to two important components, namely, 1) the plasma membrane depolarization, leading to a rise in intracellular Ca2+-concentration and activation of the specific Ca2+/calmodulin-dependent protein kinase, MLCK, and 2) Ca2+-sensitization via inhibition of MLCP (Walsh and Cole, 2013), both inducing an increase in phosphorylation of MLC20. Here, we show that there is a marked difference in the mechanical behavior of y-BAs and s-BAs and its regulation in the absence of a contractile agonist. In both y-BAs and s-BAs, there is a rapid rise in tone followed by a decline to a lower level upon stepwise stretching the vessels to IC90 during the normalization procedure as expected from the length tension relation. The tone at IC90 is stable in y-BAs during the whole experiment, i.e., there is no slow increase in tone, and basal tone recovers after washout of the contractile agonist U46619. In contrast, in s-BAs tone at IC90 is not stable, but rather slowly increases. This rise in tone we denote as stretch-induced tone. In y-BAs, such a stretch-induced tone was observed only after inhibition of NO-signaling by L-NAME.
At both ages, the developed stretch-induced tone and L-NAME induced tone were inhibited by the ROK inhibitor Y27632, putting ROK center stage in regulating tone in basilar arteries, noteworthy, even in the absence of a contractile agonist. Our experiments indicate that not only NO-PKG signaling (Lubomirov et al., 2017; Lubomirov et al., 2018), but also ROK signaling is constitutively active in BAs. In y-BAs, NO-PKG signaling, known to inhibit RhoA-ROK-MLCP signaling at different levels along the pathway, prevents the spontaneous rise in force, while in s-BAs, the balance between NO-PKG and ROK signaling is shifted in favor of the latter. The mechanisms that underlie the constitutive activity of ROK unmasked by inhibition of NO-release are currently unknown. It is possible that stretching the vessels, either by opening of stretch sensitive cation channels (Mederos y Schnitzler et al., 2011) and/or by activating G-protein receptors (Momotani et al., 2011; Takefuji et al., 2013) leads to a Ca2+-influx, which then via tyrosine kinase Pyk2 might activate Rho-ROK signaling [(Mills et al., 2015); reviewed in (Brozovich et al., 2016)]. It is of interest that in aged renal vessels phosphorylation and activity of 90 kDa ribosomal S6 kinase (RSK2) are higher than in young renal vessels (Lubomirov et al., in press). RSK2 is a non-canonical MLC20-kinase (Suizu et al., 2000; Artamonov et al., 2018) and also activates ROK (Shi et al., 2018), and hence may be responsible for stretch-induced tone in s-BAs. We are currently exploring this possibility.
The downstream targets of ROK appear to differ between y-BA and s-BA. In y-BA, ROK appears to act through the canonical pathway, i.e., ROK-induced inhibitory phosphorylation of MYPT1 at T696 and T853, resulted in an increase in phosphorylation of MLC20 [reviewed in Puetz et al. (2009)]. We do not exclude the possibility that a rise in intracellular Ca2+-concentration also leads to MLC20 phosphorylation and increase in tone. The ROK inhibitor, Y27632, lowers phosphorylation of both sites, thus increasing the activity of MLCP, which in turn lowers MLC20 phosphorylation. Y27632 prevents the increase in L-NAME-induced force. However, which one of the two residues or perhaps both are responsible for MLCP inhibition is under debate (Feng et al., 1999; Khromov et al., 2009). To address this, we used a mouse model (MYPT1-T696A/+ mice), in which the phosphorylation propensity of T696 was genetically lowered. In y-BAs, the mutation prevented L-NAME induced tension and phosphorylation of MLC20, as did Y27632.
Surprisingly, stretch-induced tone in s-BAs, although ROK dependent, appears not to be mediated by MLC20 phosphorylation. This notion is based on the observation that basal MLC20 phosphorylation but not stretch-induced tone was reduced in s-BAs carrying the MYPT1-T696A/+ mutation. In other words, stretch-induced tone develops in s-BAs despite low levels of MLC20 phosphorylation. As the vessels for phosphorylation determination were not stretched because of technical reasons, albeit pharmacological treatment was identical to the mechanical experiment, we cannot exclude the possibility that stretching the vessel would have increased MLC20-phosphorylation. These data are in contrast to our previous observation, where the T696A mutation reduced both, phosphorylation of MLC20 and tone in young and old cerebral arteries (Lubomirov et al., 2017; Lubomirov et al., 2018). It appears from our study that regulation of cerebrovascular tone in young and old age occurs predominantly through MLC20 phosphorylation, the canonical pathway, whereas in senescence these pathways play a minor role.
NM-II and PAK1, important ROK downstream-effector-proteins, regulate stretch-induced tone in senescence
Several non-canonical pathways, in addition to the canonical one, have been shown to regulate tone in tracheal (Zhang and Gunst, 2019), and vascular smooth muscle (Kim et al., 2010; Moreno-Dominguez et al., 2013) by acting through both actin and NM-II. A series of recent publications demonstrated that ROK regulates polymerization of subcortical actin at cell adhesomes in response to contractile stimuli (Zhang et al., 2016; Zhang et al., 2018). Actin polymerization and assembly are multistep processes involving a large number of other regulatory proteins.
Furthermore, according to Gunst and co-workers (Zhang et al., 2016; Zhang et al., 2018), key events are RhoA mediated NM-II filament assembly, as well as ROK mediated phosphorylation and activation of PAK1. Our study suggests that these events are also involved in regulating stretch-induced tone and L-NAME induced contractions in senescence, notably in the absence of a contractile agonist. Phosphorylation of both NM-II and PAK1 was much more prominent in s-BAs than in y-BAs and was reduced by Y27632 in s-BAs along with the attenuation of stretch-induced tone and L-NAME induced tone. We note that global inhibition of phosphatase activity using calyculin resulted in increased subcortical F-actin (Figure 11). Our findings are in keeping with a previous report that myogenic tone was associated with a small increase in filamentous actin (Moreno-Dominguez et al., 2013). Different from airway smooth muscle, where these phosphorylation events play an important role in agonist induced contraction at young age [reviewed in Zhang et al. (2015)], they appear to be silenced in cerebral arteries in young age, but they are markers of cerebrovascular senescence and intrinsically active, i.e., in the absence of a contractile agonist.
Caldesmon acts as a molecular brake on stretch-induced tone
In addition to activation of signaling molecules that eventually lead to actin polymerization, PAK1 also works through the actin binding protein caldesmon. Phosphorylation of caldesmon reverses its inhibitory effect on actomyosin MgATPase activity thereby increasing force in skinned smooth muscle at low levels of MLC20 phosphorylation (Van Eyk et al., 1998). Technically, it was not possible to monitor PAK1 mediated phosphorylation of caldesmon due to lack of commercial antibodies. Nor was it possible to inhibit caldesmon-phosphorylation by PAK1. Therefore, we used BAs from heterozygous mice, in which the Cald1 gene was ablated, thus reducing the caldesmon content to ∼50% of WT mice (Figure 9). The decrease in caldesmon level was expected to reduce its inhibitory activity. Old BAs from Cald1+/− mice exhibited stretch-induced tone at an age in which stretch-induced tone was not yet present in WT. Similar to our findings, downregulation of caldesmon in carotid arteries ex vivo by siRNA resulted in cross-bridge cycling in unstimulated tissue (Smolock et al., 2009). It is noteworthy that caldesmon immunoreactivity dramatically decreased as a long-term consequence of subarachnoid hemorrhage (Oka et al., 1996; Doi et al., 1997) a status, which is typically associated with a hyper-contractile response of the vasculature. Our finding corroborates previous findings (Smolock et al., 2009; Putz et al., 2021) showing that caldesmon acts as a molecular brake on tone in aged cerebral vessels (Figure 9). Inhibition of caldesmon by phosphorylation, e.g., by upregulation of ROK-PAK1 signaling cascade in senescence, might therefore be a mechanism that leads to a hyper-contractile state in the absence of a contractile stimulus in senescence. Interestingly, myogenic constriction in young rat cerebral arteries was not associated with altered caldesmon phosphorylation, and was proposed not to be involved in tone regulation (Moreno-Dominguez et al., 2014).
Caldesmon, which is expressed in smooth muscle cells of the walls of visceral organs as well as in the vasculature including the medial layer of human cerebral vessels (Kacem et al., 2006), has two important functional domains (Morgan and Gangopadhyay, 2001). Inhibition of actomyosin interaction is exerted by the C-terminal actin-binding domain (Katsuyama et al., 1992; Pfitzer et al., 1993a), while the N-terminal domain binds to myosin and was proposed to be important for positioning caldesmon (Hemric and Chalovich, 1990; Lee et al., 2000). Blocking of myosin binding by inhibitory peptides resulted in an increase in resting tone in permeabilized arteries (Pfitzer et al., 2005). Ablation of exon2 of the Cald1 gene in mice (CaD-ΔEx2−/−), which encodes the strong myosin binding domain, resulted in expression of a truncated caldesmon but the expression level did not differ from WT mice. In vitro binding of the truncated caldesmon, isolated from these mice, had lost myosin but not actin binding (Figure 10). However, vascular tone in o-BAs from homozygous CaD-ΔEx2−/− mice was not different from that in WT mice. This indicates that the C-terminal domain is sufficient for inhibition of force.
Taken together, our work corroborated the notion of the central role ROK has as a signaling node that orchestrates several downstream processes, essential for regulating vascular tone both in y- and s-BAs. Our observations expand previous reports that inhibition of ROK reduced vascular tone of pulmonary, cerebral, coronary, and mesenteric arteries from different species including humans, typically performed in young arteries (Fu et al., 1998; Matsumura et al., 2001; Kitazono et al., 2002; Grisk et al., 2012). However, our work suggests that the downstream ROK activated processes differ especially in senescent cerebral arteries from those in young and old ones. We propose that in y-and o-BAs the canonical pathways, i.e., inhibition of MLCP by ROK and increase in phosphorylation of MLC20, predominate tone regulation, while in senescence regulation of tone involves ROK related thin filament linked mechanisms.
Clinical relevance and perspectives
Several in vivo studies in animal models demonstrated that ROK inhibition prevented coronary and cerebral artery vasospasm and neointimal formation after balloon injury of carotid arteries (Kandabashi et al., 2000; Sato et al., 2000; Sawada et al., 2000). However, ROK inhibition may have a number of unwanted side effects. This is because ROK serves as a master regulator in non-muscle cells, regulating such diverse functions as cell-motility, -migration, and -adhesion, thrombocyte aggregation, tumor cell activity as well as the barrier function of endothelial cells (Aslan and McCarty, 2013; Constantin, 2016; Kataoka and Ogawa, 2016; Shimokawa et al., 2016).
1. Therefore, one goal of the present work was to delineate the ROK substrates NM-II and PAK1-phosphorylation as regulators of spontaneous tone generation in the senescent cerebrovascular system. Based on our study, it is reasonable to conduct further studies, investigating whether targeting of these phosphoproteins would be an approach to treat cerebrovascular dysfunction in the elderly.
2. Genetically lowering caldesmon expression leads to early onset of spontaneous contractions in o-BAs, typical for senescence. Therefore, it will be intriguing to investigate, at least in animal models, whether re-expression or overexpression of caldesmon would be able to prevent hypercontractility in cerebral vasculature.
Data availability statement
The authors acknowledge that the data presented in this study must be deposited and made publicly available in an acceptable repository, prior to publication. Frontiers cannot accept a manuscript that does not adhere to our open data policies.
Ethics statement
The animal study was reviewed and approved by State Office for Nature, Environment and Customer Protection North Rhine Westphalia: AZ 84-02.05.50.15.029.
Author contributions
Idea and design of experiments: MS, LL, and GP. Preparation and writing of manuscript: LL, MS, GP, JC, OG, NS, JH, and DM. Generation of Caldesmon Ex2−/− mice: NS, MS, and VH. Biochemical analysis of the truncated caldesmon protein: MS, MF, and JC. Western blot analysis: MB, DM, and LL. Contractile experiments: LL. G-F-actin analysis: MS and LL. Confocal microscopy: SP, GPr, and LL. Biometrical planning: DM, LL, and GP.
Funding
This work was funded by the Brandenburg Medical School Theodor Fontane (Medizinische Hochschule Brandenburg Theodor Fontane, MHB) publication fund supported by the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG). This work was supported by the German Research Association (DFG); the Center of Molecular Medicine (CMMC) of the Medical Faculty of the University of Cologne to Gabriele Pfitzer, and the German Research Research Association (DFG, project number 5400900 and SFB612) to Gabriele Pfitzer.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2023.1099278/full#supplementary-material
References
Amano, M., Ito, M., Kimura, K., Fukata, Y., Chihara, K., Nakano, T., et al. (1996). Phosphorylation and activation of myosin by Rho-associated kinase (Rho-kinase). J. Biol. Chem. 271 (34), 20246–20249. doi:10.1074/jbc.271.34.20246
Arner, A., and Pfitzer, G. (1999). Regulation of cross-bridge cycling by Ca2+ in smooth muscle. Rev. Physiol. Biochem. Pharmacol. 134, 63–146. doi:10.1007/3-540-64753-8_3
Artamonov, M. V., Sonkusare, S. K., Good, M. E., Momotani, K., Eto, M., Isakson, B. E., et al. (2018). RSK2 contributes to myogenic vasoconstriction of resistance arteries by activating smooth muscle myosin and the Na(+)/H(+) exchanger. Sci. Signal 11 (554), eaar3924. doi:10.1126/scisignal.aar3924
Aslan, J. E., and McCarty, O. J. (2013). Rho GTPases in platelet function. J. Thromb. Haemost. 11 (1), 35–46. doi:10.1111/jth.12051
Brozovich, F. V., Nicholson, C. J., Degen, C. V., Gao, Y. Z., Aggarwal, M., and Morgan, K. G. (2016). Mechanisms of vascular smooth muscle contraction and the basis for pharmacologic treatment of smooth muscle disorders. Pharmacol. Rev. 68 (2), 476–532. doi:10.1124/pr.115.010652
Chen, C. P., Chen, X., Qiao, Y. N., Wang, P., He, W. Q., Zhang, C. H., et al. (2015). In vivo roles for myosin phosphatase targeting subunit-1 phosphorylation sites T694 and T852 in bladder smooth muscle contraction. J. Physiol. 593 (3), 681–700. doi:10.1113/jphysiol.2014.283853
Cipolla, M. J., Gokina, N. I., and Osol, G. (2002). Pressure-induced actin polymerization in vascular smooth muscle as a mechanism underlying myogenic behavior. FASEB J. 16 (1), 72–76. doi:10.1096/cj.01-0104hyp
Constantin, B. (2016). Role of scaffolding proteins in the regulation of TRPC-dependent calcium entry. Adv. Exp. Med. Biol. 898, 379–403. doi:10.1007/978-3-319-26974-0_16
Doi, M., Kasuya, H., Weir, B., Cook, D. A., and Ogawa, A. (1997). Reduced expression of calponin in canine basilar artery after subarachnoid haemorrhage. Acta Neurochir. (Wien) 139 (1), 77–81. doi:10.1007/BF01850872
Eddinger, T. J., Meer, D. P., Miner, A. S., Meehl, J., Rovner, A. S., and Ratz, P. H. (2007). Potent inhibition of arterial smooth muscle tonic contractions by the selective myosin II inhibitor, blebbistatin. J. Pharmacol. Exp. Ther. 320 (2), 865–870. doi:10.1124/jpet.106.109363
Feng, J., Ito, M., Ichikawa, K., Isaka, N., Nishikawa, M., Hartshorne, D. J., et al. (1999). Inhibitory phosphorylation site for Rho-associated kinase on smooth muscle myosin phosphatase. J. Biol. Chem. 274 (52), 37385–37390. doi:10.1074/jbc.274.52.37385
Fu, X., Gong, M. C., Jia, T., Somlyo, A. V., and Somlyo, A. P. (1998). The effects of the Rho-kinase inhibitor Y-27632 on arachidonic acid-GTPgammaS-and phorbol ester-induced Ca2+-sensitization of smooth muscle. FEBS Lett. 440 (1-2), 183–187. doi:10.1016/s0014-5793(98)01455-0
Grisk, O., Schluter, T., Reimer, N., Zimmermann, U., Katsari, E., Plettenburg, O., et al. (2012). The Rho kinase inhibitor SAR407899 potently inhibits endothelin-1-induced constriction of renal resistance arteries. J. Hypertens. 30 (5), 980–989. doi:10.1097/HJH.0b013e328351d459
Gunst, S. J., and Zhang, W. (2008). Actin cytoskeletal dynamics in smooth muscle: A new paradigm for the regulation of smooth muscle contraction. Am. J. Physiol. Cell Physiol. 295 (3), C576–C587. doi:10.1152/ajpcell.00253.2008
Gusev, N. B. (2001). Some properties of caldesmon and calponin and the participation of these proteins in regulation of smooth muscle contraction and cytoskeleton formation. Biochem. (Mosc) 66 (10), 1112–1121. doi:10.1023/a:1012480829618
Hemric, M. E., and Chalovich, J. M. (1990). Characterization of caldesmon binding to myosin. J. Biol. Chem. 265 (32), 19672–19678. doi:10.1016/s0021-9258(17)45424-x
Kacem, K., Sercombe, C., Hammami, M., Vicaut, E., and Sercombe, R. (2006). Sympathectomy causes aggravated lesions and dedifferentiation in large rabbit atherosclerotic arteries without involving nitric oxide. J. Vasc. Res. 43 (3), 289–305. doi:10.1159/000093010
Kandabashi, T., Shimokawa, H., Miyata, K., Kunihiro, I., Kawano, Y., Fukata, Y., et al. (2000). Inhibition of myosin phosphatase by upregulated rho-kinase plays a key role for coronary artery spasm in a porcine model with interleukin-1beta. Circulation 101 (11), 1319–1323. doi:10.1161/01.cir.101.11.1319
Kataoka, K., and Ogawa, S. (2016). Variegated RHOA mutations in human cancers. Exp. Hematol. 44 (12), 1123–1129. doi:10.1016/j.exphem.2016.09.002
Katsuyama, H., Wang, C. L., and Morgan, K. G. (1992). Regulation of vascular smooth muscle tone by caldesmon. J. Biol. Chem. 267 (21), 14555–14558. doi:10.1016/s0021-9258(18)42075-3
Khromov, A., Choudhury, N., Stevenson, A. S., Somlyo, A. V., and Eto, M. (2009). Phosphorylation-dependent autoinhibition of myosin light chain phosphatase accounts for Ca2+ sensitization force of smooth muscle contraction. J. Biol. Chem. 284 (32), 21569–21579. doi:10.1074/jbc.M109.019729
Kim, H. R., Graceffa, P., Ferron, F., Gallant, C., Boczkowska, M., Dominguez, R., et al. (2010). Actin polymerization in differentiated vascular smooth muscle cells requires vasodilator-stimulated phosphoprotein. Am. J. Physiol. Cell Physiol. 298 (3), 559–571. doi:10.1152/ajpcell.00431.2009
Kitazono, T., Ago, T., Kamouchi, M., Santa, N., Ooboshi, H., Fujishima, M., et al. (2002). Increased activity of calcium channels and Rho-associated kinase in the basilar artery during chronic hypertension in vivo. J. Hypertens. 20 (5), 879–884. doi:10.1097/00004872-200205000-00022
Lee, Y. H., Gallant, C., Guo, H., Li, Y., Wang, C. A., and Morgan, K. G. (2000). Regulation of vascular smooth muscle tone by N-terminal region of caldesmon. Possible role of tethering actin to myosin. J. Biol. Chem. 275 (5), 3213–3220. doi:10.1074/jbc.275.5.3213
Limouze, J., Straight, A. F., Mitchison, T., and Sellers, J. R. (2004). Specificity of blebbistatin, an inhibitor of myosin II. J. Muscle Res. Cell Motil. 25 (4-5), 337–341. doi:10.1007/s10974-004-6060-7
Lubomirov, L. T., Papadopoulos, S., Filipova, D., Baransi, S., Todorovic, D., Lake, P., et al. (2018). The involvement of phosphorylation of myosin phosphatase targeting subunit 1 (MYPT1) and MYPT1 isoform expression in NO/cGMP mediated differential vasoregulation of cerebral arteries compared to systemic arteries. Acta Physiol. (Oxf). 224 (1), e13079. doi:10.1111/apha.13079
Lubomirov, L. T., Papadopoulos, S., Putz, S., Welter, J., Klockener, T., Weckmuller, K., et al. (2017). Aging-related alterations in eNOS and nNOS responsiveness and smooth muscle reactivity of murine basilar arteries are modulated by apocynin and phosphorylation of myosin phosphatase targeting subunit-1. J. Cereb. Blood Flow. Metab. 37 (3), 1014–1029. doi:10.1177/0271678X16649402
Matsumura, F., and Hartshorne, D. J. (2008). Myosin phosphatase target subunit: Many roles in cell function. Biochem. Biophys. Res. Commun. 369 (1), 149–156. doi:10.1016/j.bbrc.2007.12.090
Matsumura, Y., Kita, S., and Okui, T. (2001). Mechanisms of endothelin-1-induced potentiation of noradrenaline response in rat mesenteric artery. Clin. Exp. Pharmacol. Physiol. 28 (7), 540–544. doi:10.1046/j.1440-1681.2001.03484.x
Mederos y Schnitzler, M., Storch, U., and Gudermann, T. (2011). AT1 receptors as mechanosensors. Curr. Opin. Pharmacol. 11 (2), 112–116. doi:10.1016/j.coph.2010.11.003
Mills, R. D., Mita, M., Nakagawa, J., Shoji, M., Sutherland, C., and Walsh, M. P. (2015). A role for the tyrosine kinase Pyk2 in depolarization-induced contraction of vascular smooth muscle. J. Biol. Chem. 290 (14), 8677–8692. doi:10.1074/jbc.M114.633107
Momotani, K., Artamonov, M. V., Utepbergenov, D., Derewenda, U., Derewenda, Z. S., and Somlyo, A. V. (2011). p63RhoGEF couples Gα(q/11)-mediated signaling to Ca2+ sensitization of vascular smooth muscle contractility. Circ. Res. 109 (9), 993–1002. doi:10.1161/CIRCRESAHA.111.248898
Morano, I., Chai, G. X., Baltas, L. G., Lamounier-Zepter, V., Lutsch, G., Kott, M., et al. (2000). Smooth-muscle contraction without smooth-muscle myosin. Nat. Cell Biol. 2 (6), 371–375. doi:10.1038/35014065
Moreno-Dominguez, A., Colinas, O., El-Yazbi, A., Walsh, E. J., Hill, M. A., Walsh, M. P., et al. (2013). Ca2+ sensitization due to myosin light chain phosphatase inhibition and cytoskeletal reorganization in the myogenic response of skeletal muscle resistance arteries. J. physiology 591 (5), 1235–1250. doi:10.1113/jphysiol.2012.243576
Moreno-Dominguez, A., El-Yazbi, A. F., Zhu, H. L., Colinas, O., Zhong, X. Z., Walsh, E. J., et al. (2014). Cytoskeletal reorganization evoked by Rho-associated kinase- and protein kinase C-catalyzed phosphorylation of cofilin and heat shock protein 27, respectively, contributes to myogenic constriction of rat cerebral arteries. J. Biol. Chem. 289 (30), 20939–20952. doi:10.1074/jbc.M114.553743
Morgan, K. G., and Gangopadhyay, S. S. (2001). Invited review: Cross-bridge regulation by thin filament-associated proteins. J. Appl. Physiol. (1985) 91 (2), 953–962. doi:10.1152/jappl.2001.91.2.953
Oka, Y., Ohta, S., Todo, H., Kohno, K., Kumon, Y., and Sakaki, S. (1996). Protein synthesis and immunoreactivities of contraction-related proteins in smooth muscle cells of canine basilar artery after experimental subarachnoid hemorrhage. J. Cereb. Blood Flow. Metab. 16 (6), 1335–1344. doi:10.1097/00004647-199611000-00031
Pfitzer, G., Fischer, W., and Chalovich, J. M. (1993). Phosphorylation-contraction coupling in smooth muscle: Role of caldesmon. Adv. Exp. Med. Biol. 332, 195–202. Discussion-3. doi:10.1007/978-1-4615-2872-2_18
Pfitzer, G., Schroeter, M., Hasse, V., Ma, J., Rosgen, K. H., Rosgen, S., et al. (2005). Is myosin phosphorylation sufficient to regulate smooth muscle contraction? Adv. Exp. Med. Biol. 565, 319–328. doi:10.1007/0-387-24990-7_24
Pfitzer, G., Zeugner, C., Troschka, M., and Chalovich, J. M. (1993). Caldesmon and a 20-kDa actin-binding fragment of caldesmon inhibit tension development in skinned gizzard muscle fiber bundles. Proc. Natl. Acad. Sci. U. S. A. 90 (13), 5904–5908. doi:10.1073/pnas.90.13.5904
Puetz, S., Lubomirov, L. T., and Pfitzer, G. (2009). Regulation of smooth muscle contraction by small GTPases. Physiol. (Bethesda) 24, 342–356. doi:10.1152/physiol.00023.2009
Putz, S., Barthel, L. S., Frohn, M., Metzler, D., Barham, M., Pryymachuk, G., et al. (2021). Caldesmon ablation in mice causes umbilical herniation and alters contractility of fetal urinary bladder smooth muscle. J. Gen. Physiol. 153 (7), e202012776. doi:10.1085/jgp.202012776
Rhee, A. Y., Ogut, O., and Brozovich, F. V. (2006). Nonmuscle myosin, force maintenance, and the tonic contractile phenotype in smooth muscle. Pflugers Arch. 452 (6), 766–774. doi:10.1007/s00424-006-0091-4
Ruegg, J. C., and Pfitzer, G. (1991). Contractile protein interactions in smooth muscle. Blood Vessels 28 (1-3), 159–163. doi:10.1159/000158856
Sato, M., Tani, E., Fujikawa, H., and Kaibuchi, K. (2000). Involvement of Rho-kinase-mediated phosphorylation of myosin light chain in enhancement of cerebral vasospasm. Circ. Res. 87 (3), 195–200. doi:10.1161/01.res.87.3.195
Sawada, N., Itoh, H., Ueyama, K., Yamashita, J., Doi, K., Chun, T. H., et al. (2000). Inhibition of rho-associated kinase results in suppression of neointimal formation of balloon-injured arteries. Circulation 101 (17), 2030–2033. doi:10.1161/01.cir.101.17.2030
Shi, G. X., Yang, W. S., Jin, L., Matter, M. L., and Ramos, J. W. (2018). RSK2 drives cell motility by serine phosphorylation of LARG and activation of Rho GTPases. Proc. Natl. Acad. Sci. U. S. A. 115 (2), E190–E199. doi:10.1073/pnas.1708584115
Shimokawa, H., Sunamura, S., and Satoh, K. (2016). RhoA/Rho-kinase in the cardiovascular system. Circ. Res. 118 (2), 352–366. doi:10.1161/CIRCRESAHA.115.306532
Smolock, E. M., Trappanese, D. M., Chang, S., Wang, T., Titchenell, P., and Moreland, R. S. (2009). siRNA-mediated knockdown of h-caldesmon in vascular smooth muscle. Am. J. Physiol. Heart Circ. Physiol. 297 (5), H1930–H1939. doi:10.1152/ajpheart.00129.2009
Somlyo, A. P., and Somlyo, A. V. (2003). Ca2+ sensitivity of smooth muscle and nonmuscle myosin II: modulated by G proteins, kinases, and myosin phosphatase. Physiol. Rev. 83 (4), 1325–1358. doi:10.1152/physrev.00023.2003
Suizu, F., Ueda, K., Iwasaki, T., Murata-Hori, M., and Hosoya, H. (2000). Activation of actin-activated MgATPase activity of myosin II by phosphorylation with MAPK-activated protein kinase-1b (RSK-2). J. Biochem. 128 (3), 435–440. doi:10.1093/oxfordjournals.jbchem.a022771
Takefuji, M., Krüger, M., Sivaraj, K. K., Kaibuchi, K., Offermanns, S., and Wettschureck, N. (2013). RhoGEF12 controls cardiac remodeling by integrating G protein- and integrin-dependent signaling cascades. J. Exp. Med. 210 (4), 665–673. doi:10.1084/jem.20122126
Van Eyk, J. E., Arrell, D. K., Foster, D. B., Strauss, J. D., Heinonen, T. Y., Furmaniak-Kazmierczak, E., et al. (1998). Different molecular mechanisms for Rho family GTPase-dependent, Ca2+-independent contraction of smooth muscle. J. Biol. Chem. 273 (36), 23433–23439. doi:10.1074/jbc.273.36.23433
Walsh, M. P., and Cole, W. C. (2013). The role of actin filament dynamics in the myogenic response of cerebral resistance arteries. J. Cereb. Blood Flow. Metab. 33 (1), 1–12. doi:10.1038/jcbfm.2012.144
Yamin, R., and Morgan, K. G. (2012). Deciphering actin cytoskeletal function in the contractile vascular smooth muscle cell. J. Physiol. 590 (17), 4145–4154. doi:10.1113/jphysiol.2012.232306
Yuen, S., Ogut, O., and Brozovich, F. V. (2011). MYPT1 protein isoforms are differentially phosphorylated by protein kinase G. J. Biol. Chem. 286 (43), 37274–37279. doi:10.1074/jbc.M111.282905
Zhang, W., Bhetwal, B. P., and Gunst, S. J. (2018). Rho kinase collaborates with p21-activated kinase to regulate actin polymerization and contraction in airway smooth muscle. J. physiology 596 (16), 3617–3635. doi:10.1113/JP275751
Zhang, W., and Gunst, S. J. (2019). Molecular mechanisms for the mechanical modulation of airway responsiveness. J. Eng. Sci. Med. Diagn Ther. 2 (1), 0108051–0108058. doi:10.1115/1.4042775
Zhang, W., and Gunst, S. J. (2017). Non-muscle (NM) myosin heavy chain phosphorylation regulates the formation of NM myosin filaments, adhesome assembly and smooth muscle contraction. J. physiology 595 (13), 4279–4300. doi:10.1113/JP273906
Zhang, W., Huang, Y., and Gunst, S. J. (2016). p21-Activated kinase (Pak) regulates airway smooth muscle contraction by regulating paxillin complexes that mediate actin polymerization. J. physiology 594 (17), 4879–4900. doi:10.1113/JP272132
Keywords: stretch-induced tone, senescence, basilar artery, non-muscle myosin, caldesmon
Citation: Lubomirov LT, Schroeter MM, Hasse V, Frohn M, Metzler D, Bust M, Pryymachuk G, Hescheler J, Grisk O, Chalovich JM, Smyth NR, Pfitzer G and Papadopoulos S (2023) Dual thick and thin filament linked regulation of stretch- and L-NAME-induced tone in young and senescent murine basilar artery. Front. Physiol. 14:1099278. doi: 10.3389/fphys.2023.1099278
Received: 15 November 2022; Accepted: 09 March 2023;
Published: 28 March 2023.
Edited by:
Steven Fisher, University of Maryland, United StatesReviewed by:
Rudolf Schubert, Augsburg University, GermanyAvril Somlyo, University of Virginia, United States
Copyright © 2023 Lubomirov, Schroeter, Hasse, Frohn, Metzler, Bust, Pryymachuk, Hescheler, Grisk, Chalovich, Smyth, Pfitzer and Papadopoulos. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lubomir T. Lubomirov, Lubomir.Lubomirov@mhb-fontane.de
†These authors have contributed equally to this work