ROS Dependent Wnt/β-Catenin Pathway and Its Regulation on Defined Micro-Pillars—A Combined In Vitro and In Silico Study

The physico-chemical surface design of implants influences the surrounding cells. Osteoblasts on sharp-edged micro-topographies revealed an impaired cell phenotype, function and Ca2+ mobilization. The influence of edges and ridges on the Wnt/β-catenin pathway in combination with the cells’ stress response has not been clear. Therefore, MG-63 osteoblasts were studied on defined titanium-coated micro-pillars (5 × 5 × 5 µm) in vitro and in silico. MG-63s on micro-pillars indicated an activated state of the Wnt/β-catenin pathway. The β-catenin protein accumulated in the cytosol and translocated into the nucleus. Gene profiling indicated an antagonism mechanism of the transcriptional activity of β-catenin due to an increased expression of inhibitors like ICAT (inhibitor of β-catenin and transcription factor-4). Cells on pillars produced a significant reactive oxygen species (ROS) amount after 1 and 24 h. In silico analyses provided a detailed view on how transcriptional activity of Wnt signaling is coordinated in response to the oxidative stress induced by the micro-topography. Based on a coordinated expression of regulatory elements of the Wnt/β-catenin pathway, MG-63s are able to cope with an increased accumulation of β-catenin on micro-pillars and suppress an unintended target gene expression. Further, β-catenin may be diverted into other signaling pathways to support defense mechanisms against ROS.


Introduction
Biofunctional materials which replace damaged tissue and stimulate cell regeneration are currently in great demand and will be in the future. Titanium (Ti) and titanium-based alloys are the most common metallic materials in implant applications due to their good biocompatibility, mechanical compliance, and high corrosion resistance [1][2][3]. The rapid cellular acceptance in the tissue can be optimized by physico-chemical surface functionalization, e.g., topographical cues [4][5][6]. Via adhesion receptors, cells are able to perceive the material surface properties and also the extra cellular matrix
The software ZEN 2011 (ZEISS Efficient Navigation, ZEN 2011 SP4, black or blue edition, Carl Zeiss) was used for image analysis. For MG-63s on micro-pillars, several images were acquired in the "Z-stack" function (3-10 slices with an interval of 1 µm, pinhole 0.7 µm) and then generated a 3D-reconstruction (overlay). Confocal microscopy with 3D Z-stack generation was performed to visualize the protein localization on top of and between the micro-pillars to understand protein localization that was not limited to a single dimension. Using the "Profile" function, intensity profiles of the individual fluorescence's of certain regions in the image can be created.

Bio-Plex Luminex Bead Multiplex Analysis
The secreted proteins osteocalcin (OCN) and osteopontin (OPN) from lysates of the cell culture supernatant after 96 h were analyzed using MILLIPLEX MAP Kit for Human Bone Magnetic Bead Panel (Merck KGAA) according to the manufacturer's recommendations. Briefly, the bead mix was transferred to each well of the 96-well plate and mixed with lysates over night at 4 • C shaking with 300 rpm. Subsequently, incubation was performed with a detection antibody solution (30 min at RT, 300 rpm) and afterwards with a Streptavidin-Phycoerythrin solution (30 min at RT, 300 rpm). The protein quantification, mean fluorescence intensity (MFI), was measured with the Bio-Plex ® 200 System (Bio-Rad) and the Bio-Plex Manager ™ 4.1.1 software (Bio-Rad). The MFI values were normalized to the protein concentration determined by Bradford method.

Enzyme-Linked Immunosorbent Assay (ELISA)
For quantitative protein detection of ICAT in the cell fraction lysates of MG-63s within 24 h, the ICAT ELISA kit (human beta catenin interacting protein 1 (CTNNBIP1); MyBioSource, Inc., San Diego, CA, USA) was used according to the manufacturer's instructions. Briefly, the 96-well pre-coated microtiter plate was loaded with equal amounts of the protein as well as the standards. Lysates were labeled with Cells 2020, 9, 1784 6 of 24 CTNNBIP1-HRP conjugates for 1 h and subsequently incubated with the substrate for HRP enzyme at 37 • C for 20 min. Immediately after stopping the reaction, the optical density (O.D., intensity of the color) could be measured using Anthos Reader 2010 (Anthos Mikrosysteme GmbH, Friesoythe, Germany) at 450 nm. For the calculation of ICAT protein expression level, a standard curve was plotted relating the O.D. to the standard concentration.

Detection of Reactive Oxygen Species (ROS) and Lipid Peroxidation
The production of intracellular ROS were detected using DCF-DA Cellular ROS detection assay kit (Abcam, Cambridge, UK). Therefore, a trypsinated MG-63 cell pellet was resuspended in a dye solution consisting of 2 ,7 -dichlorofluorescein diacetate (DCF-DA, 20 µM) and then shaken (300 rpmi) at 37 • C for 30 min in the dark. The stained MG-63s were resuspended in complete DMEM (without phenol red; Life technologies GmbH). Afterwards, 1 × 10 5 osteoblasts were cultivated on the samples within 24 h. In the presence of an ROS, the DCF-DA is oxidized to the fluorescence compound DCF [13]. The CellROX reagents (Thermo Fisher Scientific) were used to evaluate ROS primarily in the nucleus and mitochondria (CellROX ® Green) or cytosol (CellROX ® Orange) according to the manufacturer's instructions. After 24 h adherence on the samples, cells were washed with PBS and loaded with 5 µM of one fluorogenic probe at 37 • C for 30 min.

Statistics
All in vitro experiments were repeated at least three times with independently passaged cells. To evaluate statistical differences and to prepare the graphs, the GraphPad Prism 7 software for Windows (GraphPad Software Inc., La Jolla, CA, USA) was used. Data were normalized either to unstructured bulk material (Ref) or time point (1 h). After determining the normal distribution by Kolmogorov-Smirnov, subsequent statistical test was done: the Mann-Whitney U test (protein expression by FACS; osteogenic marker analysis), paired t-test (ROS and lipid peroxidation detection) or the multiple t-test (β-catenin, ICAT in subcellular fractions). Herein, a probability value of p < 0.05 was accepted as indicating a significant difference. Data were expressed as mean ± s.e.m. (standard error of the mean).

In Silico Experiments: Modeling the Wnt/β-Catenin Pathway in Osteoblasts
The simulation model is defined in ML-Rules, a multi-level, rule-based modeling language [33]. ML-Rules allows the specification and simulation of complex systems biology models at multiple levels of organization. The semantics of ML-Rules is based on continuous-time Markov chains (CTMC) [34]. For all simulation experiments performed in this study, a stochastic simulator that is based on the Stochastic Simulation Algorithm (SSA) was used [34]. For the specification and execution of the model SESSL (The Simulation Experiment Specification on a Scala Layer) [35] was used. Based on the bindings between SESSL and ML-Rules [36], various simulation experiments have been executed. The specifications of these experiments in SESSL, the model in ML-Rules and the R-Script to replicate the figures shown in the paper are available in the GitHub repository https://github.com/SFB-ELAINE/SI_Staehlke_MDPI_Cells. A table containing all parameter values employed in the (fitted) model as well as the model implementation in ML-Rules can be found as part of the supplementary material (Table S2 and Figure S4).

Topography-Dependent β-Catenin Accumulation and Translocation into the Nucleus
The stabilization, accumulation, and translocation of the β-catenin into the nucleus is a marker for the activation of the Wnt/β-catenin pathway [15,16,23] and was examined by flow cytometry, laser scanning microscopy (LSM) and Western Blot (subcellular fractionation) within 24 h MG-63 cell cultivation on defined sharp-edged micro-pillars (P5).
In flow cytometry a significantly increase in total protein expression of β-catenin in cells on P5 compared with Ref after 24 h was observed ( Figure 1A). Additionally, the total β-catenin protein expression increased over time in cells on P5 from 1 h to 24 h, but not on Ref ( Figure 1B).    On Ref β-catenin appears evenly distributed throughout the cytosol with an increased localization along the cell membrane in cell-cell contacts. On P5 a general higher occurrence was detected and β-catenin was distributed mainly in the cytosol ( Figure 1C). Also, translocation of β-catenin into the nucleus was observed on P5 ( Figure 1D).
Also the quantification of β-catenin protein level in MG-63 subcellular fractions using Western Blot revealed a slightly higher level in the cytosol as well as a significantly increased expression in the nucleus on P5 versus Ref. No topographical influence could be observed in the protein expression of β-catenin in the membrane fraction ( Figure 1E,F).
Taken together, MG-63 cells on P5 indicated an activated state of the Wnt/β-catenin pathway due to the accumulation of the key effector β-catenin in the cytosol and their transduction in the nucleus after 24 h.

Gene Regulation of Wnt Signaling Components Dependent on Micro-Pillars
In order to obtain information about the overall composition of the analyzed complex microarray datasets and to determine the direction of variability, the principal component analysis (PCA) was performed (Figure 2A; PCA1-3 of 73.1%). The approach showed differences between the respective group gene expression profiles ( Figure 2A). This indicates that the group profiles are divergent and separately clustered, which allowed for analysis of the differences in gene regulation ( Figure 2A). In addition, similarities of the individual values within a group were shown, indicating proximity within the profile.
In order to obtain information about the overall composition of the analyzed complex microarray datasets and to determine the direction of variability, the principal component analysis (PCA) was performed (Figure 2A; PCA1-3 of 73.1%). The approach showed differences between the respective group gene expression profiles (Figure 2A). This indicates that the group profiles are divergent and separately clustered, which allowed for analysis of the differences in gene regulation ( Figure 2A). In addition, similarities of the individual values within a group were shown, indicating proximity within the profile.  The signal space transformation robust multi-array average (SST-RMA) algorithm identified expression values for 21,448 annotated gene clusters, which were further filtered with the criteria fold change cutoff, fold change > 1.5 or <−1.5, linear models for microarray data (LIMMA) p < 0.05, and for the identification of differentially regulated genes ( Figure 2B). Of these, 315 (1.47%) genes were detected and marked as differently regulated in osteoblasts on P5 compared to Ref under the criteria fold change cutoff. In summary, 157 genes (49.8%) were up-regulated due to P5, while 158 genes (50.2%) were down-regulated due to P5 compared to Ref ( Figure 2B). Under the filter criteria, 6 genes have been identified that were significantly (p < 0.05) differently regulated regarding the Wnt/β-catenin pathway (Table 1). The datasets were evaluated according to the pathways tool WikiPathways. Figure 3 summarizes these representations of Wnt pathways. Both non-canonical Wnt pathways-Wnt/calcium and Wnt/planar cell polarity (PCP)-indicated no altered mRNA expression profile. Regulated mRNA profiles were apparent solely in the Wnt/β-catenin pathway (Table 1). Interestingly, the mRNA expression of β-catenin remained unchanged (fold change 1.08). Further Wnt/β-catenin pathway components such as ligands Wnt, receptors like Frizzled (Fz) or low-density lipoprotein receptor-related protein (LRP) and Dishevelled (DVL) revealed no significant regulation in gene expression due to the underlying micro-topography. Also, target genes such as AXIN2 or osteoblast-specific ones were not influenced on the mRNA expression level (Table S1). All in all, 5 out of 6 significantly expressed components were upregulated in cells on P5. The activated state of the Wnt/β-catenin pathway can further speculated due to the upregulation of transcription factor 4 (TCF4, fold change 1.55, p < 0.0397). In addition, the increased and abundant expression of WLS (wntless, fold change 1.87, p < 0.0299) was observed on P5. Remarkably, osteoblasts on P5 exhibited an upregulation of antagonists of the Wnt/β-catenin pathway, such as SOST (sclerostin, fold change 1.57, p < 0.0006). Further significantly upregulated inhibitors were ICAT (inhibitor of β-catenin and TCF-4, CTNNBIP1; fold change 1.65, p < 0.0164 by weak intensity signals) as well as SOX17 (SRY-box transcription factor 17, fold change 1.58, p < 0.0066), which finally antagonize the transcriptional activation of β-catenin. Another Wnt antagonist, the SFRP4 (secreted frizzled-related protein 4, fold change −1.95, p < 0.0062), showed a significantly lower expression on P5 versus Ref.

Topography-Dependent Protein Expression of β-Catenin Inhibitors ICAT and SOX17
The expression of β-catenin antagonists were also upregulated at protein level in MG-63s on P5 after 24 h (Figure 4). The fluorescence images revealed that on Ref the ICAT distribution appears evenly throughout the whole cell. In contrast, a higher fluorescence signal could be observed in the nucleus of cells on P5 ( Figure S5) ( Figure 4A left). In agreement with the fluorescence images, the quantification of ICAT protein in MG-63 subcellular fractions indicated a significantly increased nuclear and slightly higher cytosolic expression on P5 versus Ref ( Figure 4A right). Also, the total protein expression of inhibitor SOX17 in MG-63s on P5 indicated a significant higher expression by both-fluorescence imaging and flow cytometry ( Figure 4B).

Topography-Dependent Reduction of the Osteoblast-Specific Function after 96 h
In previous studies a significant reduced expression of osteogenic markers could be observed for osteoblasts on micro-pillars after 24 h [9,12,13]. In this study the protein expression of collagen type I (Col1), fibronectin (FN) as well as protein secretion of osteocalcin (OCN) and osteopontin (OPN) displayed a significant reduction after 96 h on P5 ( Figure 5).

Topography-Dependent Reduction of the Osteoblast-Specific Function after 96 h
In previous studies a significant reduced expression of osteogenic markers could be observed for osteoblasts on micro-pillars after 24 h [9,12,13]. In this study the protein expression of collagen type I (Col1), fibronectin (FN) as well as protein secretion of osteocalcin (OCN) and osteopontin (OPN) displayed a significant reduction after 96 h on P5 ( Figure 5).

Micro-Pillars Induce Oxidative Stress
The micro-pillars (P5) induced the reactive oxygen species (ROS) production and as consequence lipid peroxidation in MG-63s ( Figure 6). To investigate the ROS production after 24 h, first the DCF-DA Cellular ROS detection assay indicated elevated ROS (1.19-fold ± 0.027) in MG-63s on P5 versus Ref ( Figure 6A). In addition, it can be demonstrated that P5 provoked the ROS production 1.12-fold (±0.021, * p < 0.05) just after 1 h cultivation (data not shown). Moreover, intracellular ROS was determined significantly in the nucleus, the mitochondria (both by CellROX ® Green) as well as in the cytosol (CellROX ® Orange) ( Figure 6B,C). Figure 6D shows that the oxidation state of the membrane lipids for MG-63s on P5 is significantly lower on Ref, which implies an increase in lipid peroxidation due to ROS.
Antioxidative enzymes like superoxide dismutase (SOD) eliminate ROS by catalyzing reactions. The protein expression of one of these ROS-scavenging enzymes-SOD3-was significantly increased (microarray technology: fold change 1.32, p < 0.0345; Table S1).

Micro-Pillars Induce Oxidative Stress
The micro-pillars (P5) induced the reactive oxygen species (ROS) production and as consequence lipid peroxidation in MG-63s ( Figure 6). To investigate the ROS production after 24 h, first the DCF-DA Cellular ROS detection assay indicated elevated ROS (1.19-fold ± 0.027) in MG-63s on P5 versus Ref ( Figure 6A). In addition, it can be demonstrated that P5 provoked the ROS production 1.12-fold (±0.021, * p < 0.05) just after 1 h cultivation (data not shown). Moreover, intracellular ROS was determined significantly in the nucleus, the mitochondria (both by CellROX ® Green) as well as in the cytosol (CellROX ® Orange) ( Figure 6B,C). Figure 6D shows that the oxidation state of the membrane lipids for MG-63s on P5 is significantly lower on Ref, which implies an increase in lipid peroxidation due to ROS.
Antioxidative enzymes like superoxide dismutase (SOD) eliminate ROS by catalyzing reactions. The protein expression of one of these ROS-scavenging enzymes-SOD3-was significantly increased (microarray technology: fold change 1.32, p < 0.0345; Table S1).

In Silico
To complement the experimental results obtained in this study, an in silico model of the Wnt/β-catenin pathway was developed. The model allows a quantitative exploration of the complex interplay between ROS, the inhibitors and transcription factors that are differentially expressed on micro-topography (P5) compared with unstructured samples (Ref). The model and the assumptions derived from the experimental results are briefly described below.
According to the in vitro observations, P5 induce the production of ROS in MG-63s (Figure 7). Several studies provide evidence of a cross-talk between redox and Wnt/β-catenin signaling [37][38][39][40]. Accordingly, elevated ROS levels release the redox-sensitive association between Nucleoredoxin (NRX) and DVL leading to a spontaneous increase of unbound DVL molecules, which immediately form a binding platform for AXIN, the main component of the β-catenin degradation complex [38,39]. DVL inhibits thereby β-catenin degradation, resulting in an accumulation of β-catenin in the cytosol and nucleus in a Wnt-independent manner [37]. This pathway is incorporated and validated in the in silico model presented in [28]. It forms the basis for the in silico study and is further extended. micro-topography (P5) compared with unstructured samples (Ref). The model and the assumptions derived from the experimental results are briefly described below.
According to the in vitro observations, P5 induce the production of ROS in MG-63s (Figure 7). Several studies provide evidence of a cross-talk between redox and Wnt/β-catenin signaling [37][38][39][40]. Accordingly, elevated ROS levels release the redox-sensitive association between Nucleoredoxin (NRX) and DVL leading to a spontaneous increase of unbound DVL molecules, which immediately form a binding platform for AXIN, the main component of the β-catenin degradation complex [38,39]. DVL inhibits thereby β-catenin degradation, resulting in an accumulation of β-catenin in the cytosol and nucleus in a Wnt-independent manner [37]. This pathway is incorporated and validated in the in silico model presented in [28]. It forms the basis for the in silico study and is further extended.
It is assumed that there is no extracellular Wnt stimulus; i.e., Wnt proteins are neither part of the medium nor being produced and secreted by the cells. This assumption is supported by the mRNA expression profile, which provides no significant difference between Ref and P5 in the expression of Wnt ligands. Therefore, the model disregards any membrane-associated signaling components, and considers solely the intracellular signaling components.
Next, the intracellular components that are differentially expressed on P5, i.e., ICAT, SOX17 and TCF, were added to the model. The role of TCF was thoroughly described in the original Wnt model by Lee et al. [41] and can be seamlessly incorporated in the model with the corresponding concentration values and reaction rate constants. With regard to ICAT and SOX17 and their regulatory mechanisms, applied fitting have been used to obtain the corresponding dynamics ( Figure  8). Please note, for simplicity and based on their similar function, the homologues AXIN and AXIN2 are both considered as AXIN in the model. . Note that double arrows represent reversible reactions/state changes. Reactions that occur more than once in the figure actually refer to one and the same reaction and have been included to increase readability. This applies in particular for complex formation, where for each molecule a separate arrow points to the final complex, while referring to the one reaction (e.g., R11: β-cat + ICAT ↔ β-cat/ICAT, β-cat (β-catenin), ICAT (inhibitor of β-catenin and TCF-4)). Arrows (or horizontal lines) that point towards another arrow indicate a promoting (or attenuating) effect of the molecule on the given reaction (R8, R14, R16). For example, SOX17 (SRY-box transcription factor 17) promotes the degradation of the β-catenin/TCF (transcription factor) complex (R14). Phosphorylated AXIN is labeled by a white P on a black circle while active, self-aggregated Dishevelled (DVL) complexes are indicated by an orange A. Please note that for some reactions the molecule's state of activity is not relevant, e.g., the phosphorylation state of AXIN has no impact on its binding to activated Dvl (Dvl [A]). Please see the legend for a detailed description of the graphical elements used in the figure.
It is assumed that there is no extracellular Wnt stimulus; i.e., Wnt proteins are neither part of the medium nor being produced and secreted by the cells. This assumption is supported by the mRNA expression profile, which provides no significant difference between Ref and P5 in the expression of Wnt ligands. Therefore, the model disregards any membrane-associated signaling components, and considers solely the intracellular signaling components.
Next, the intracellular components that are differentially expressed on P5, i.e., ICAT, SOX17 and TCF, were added to the model. The role of TCF was thoroughly described in the original Wnt model by Lee et al. [41] and can be seamlessly incorporated in the model with the corresponding concentration values and reaction rate constants. With regard to ICAT and SOX17 and their regulatory mechanisms, applied fitting have been used to obtain the corresponding dynamics ( Figure 8). Please note, for simplicity and based on their similar function, the homologues AXIN and AXIN2 are both considered as AXIN in the model. . Schematic representation of the implemented model of the Wnt-independent β-catenin pathway. Arrows indicate the state change of the affected molecule, such as synthesis or degradation reactions (R1, R6, R9, R14, R15), complex formation or dissociation (R2, R3, R3r, R5, R11, R13), activation or inactivation (R4, R7), and shuttling (R10, R12). Note that double arrows represent reversible reactions/state changes. Reactions that occur more than once in the figure actually refer to one and the same reaction and have been included to increase readability. This applies in particular for complex formation, where for each molecule a separate arrow points to the final complex, while referring to the one reaction (e.g., R11: β-cat + ICAT ↔ β-cat/ICAT, β-cat (β-catenin), ICAT (inhibitor of β-catenin and TCF-4)). Arrows (or horizontal lines) that point towards another arrow indicate a promoting (or attenuating) effect of the molecule on the given reaction (R8, R14, R16). For example, SOX17 (SRY-box transcription factor 17) promotes the degradation of the β-catenin/TCF (transcription factor) complex (R14). Phosphorylated AXIN is labeled by a white P on a black circle while active, self-aggregated Dishevelled (DVL) complexes are indicated by an orange A. Please note that for some reactions the molecule's state of activity is not relevant, e.g., the phosphorylation state of AXIN has no impact on its binding to activated Dvl (Dvl [A]). Please see the legend for a detailed description of the graphical elements used in the figure. Here was observed a small, but not significant accumulation of β-catenin, while the target gene expression of AXIN stays constant after a short initial increase (similar to P5).
Wnt signaling is one of the central signaling pathway that controls cell fate specification [14,16,21,56] and to play a decisive role in bone tissue homeostasis [15,28] and in osteoblast differentiation [57][58][59][60]. Wnt signaling consists of two major pathway subgroups defined by their interaction with β-catenin: canonical or Wnt/β-catenin, and non-canonical pathways-Wnt/planar cell polarity (PCP) pathway and Wnt/calcium pathway [14,53,61,62]. Some studies indicated a specificity of Wnt signaling: the activation of both non-canonical Wnt pathways´ components may antagonize the Wnt/β-catenin pathway and vice versa [11,[17][18][19][20]61,63,64]. Due to the earlier findings that actin is re-organized and Ca 2+ -release is decreased in cells on micro-pillars [9,12,13], it is assumed that both of the non-canonical signaling pathways are suppressed (Scheme 1). The question has been The simulation of the model with higher reactive oxygen species (ROS) production, as observed on P5 (see Figure 6), results in a moderate accumulation of β-catenin, whereas the target gene expression of AXIN stays constant after a short initial increase. (B) The simulation of the model with no ROS production, in contrast to P5. Here was observed a small, but not significant accumulation of β-catenin, while the target gene expression of AXIN stays constant after a short initial increase (similar to P5).
Wnt signaling is one of the central signaling pathway that controls cell fate specification [14,16,21,56] and to play a decisive role in bone tissue homeostasis [15,28] and in osteoblast differentiation [57][58][59][60]. Wnt signaling consists of two major pathway subgroups defined by their interaction with β-catenin: canonical or Wnt/β-catenin, and non-canonical pathways-Wnt/planar cell polarity (PCP) pathway and Wnt/calcium pathway [14,53,61,62]. Some studies indicated a specificity of Wnt signaling: the activation of both non-canonical Wnt pathways' components may antagonize the Wnt/β-catenin pathway and vice versa [11,[17][18][19][20]61,63,64]. Due to the earlier findings that actin is re-organized and Ca 2+ -release is decreased in cells on micro-pillars [9,12,13], it is assumed that both of the non-canonical signaling pathways are suppressed (Scheme 1). The question has been raised if the suppression of non-canonical Wnt pathways then could result in an activated state of Wnt/β-catenin pathway?
In the present study, the key effector of the Wnt/β-catenin pathway, the β-catenin-protein, indicated a very pronounced and effective response of MG-63s to P5 after 24 h. β-catenin protein level is strongly enhanced on P5 in cells at time periods of 1 to 24 h. β-catenin in MG-63s on the Ref appeared evenly distributed throughout the cytosol, with a stronger signal along the cell border. This can be attributed to the presence of cell-to-cell contacts due to the increased spreading area, accompanied by increased concentration of cadherin to mediate cell-cell adherence [10,16,65]. Some studies reported that N-cadherin expression is inversely related to the intracellular β-catenin level [55,57,65,66]. The micro-pillars inhibited the spreading and thus the cells were separated from each other, so that fewer cell-to-cell contacts could be established [9,12,60,67]. In addition, an increased accumulation of the β-catenin protein in the cytosol and a translocation into the nucleus was seen in MG-63s on P5. Several studies found also enhance cytosolic and nuclear β-catenin protein levels on topography [11,28,62,65,68] as marker for the activation of the Wnt/β-catenin signaling [15,16,23].
Many studies associate an activated state of the Wnt/β-catenin pathway with increased osteogenesis due to topography [11,18,45,46,54,61,69,70]. In previous work, however, a significant reduced expression of osteoblast function markers such as alkaline phosphatase (ALP), collagen type I (Col1), fibronectin (FN) and osteocalcin (OCN) could be observed on sharp-edged micro-pillars after 24 h [9,12,13]. This initial decrease in osteogenic markers was also confirmed after 96 h cultivation in the present study.
In one study it was found that the topographical structure has an influence on the Wnt/β-catenin pathway: Schernthaner et al. [71] described in endothelial cells that structured polymers induces the translocation of β-catenin into the nucleus, while the ripple-(wave-like) nanostructure further induces the activation of specific β-catenin target genes, the initiation of β-catenin transcriptional activation was not enhanced on a wall nanostructure (pillar-like, dimension 1 µm) and was thus crucial for cell fate. Another important question in mechanotransduction arose: which effect does a defined topography exert on the β-catenin transcriptional activation?
Transcriptome analyses by microarray technology was used to identify the full array of genes whose expression changes under topographical features. In general, the regulation of mRNA (messenger ribonucleic acid) expression in MG-63s in Wnt signaling was only observed in Wnt/β-catenin pathway components. The most important components such as Wnt, receptors or target genes were topography independent after 24 h (Table S1). Interestingly, an increase of the Frizzled-5 receptor on the protein level could be detected over time due to the micro-pillars, which was not visible on the mRNA level ( Figure S6). Furthermore, the mRNA expression of β-catenin was also not affected. However, a distinct impact was detected at the protein level-β-catenin accumulation on P5-which implies a successful activation of the Wnt/β-catenin pathway. Galli et al. [46] demonstrated that topography does not directly regulate the transcription of β-catenin in C2C12 myoblast cells, but rather reduces the expression of AXIN2, which controls the degradation of β-catenin. The mRNA level exhibited a reduced expression of AXIN2 on P5 which was not significant (fold change −1.82, p < 0.3533).
In the present study, six genes were identified as significantly regulated with respect to the Wnt/β-catenin pathway. The transmembrane protein Wntless displayed the highest gene expression (WLS) and abundancy in MG-63s on P5, which is implicated in the extracellular secretion of Wnt ligands [16,23,60]. Zhong et al. [72] could observe that the Wntless inactivation downregulated the protein level of β-catenin. In addition, the transcription factor 4 (TCF4) was significantly enhanced on P5. The proteins T-cell factor/lymphoid enhancer factor (TcF/LEF) family of DNA-bound transcription factors is the major partner for β-catenin that triggers specific expression of Wnt target genes [15,16,23,26]. The stabilization and nuclear accumulation of β-catenin causes TCF4 to form a complex and recruit other co-activators for gene activation, which ultimately control the cell fate [46]. Despite the upregulated TCF4 mRNA level, the expression of target genes (like AXIN2) was topographically unchanged. Otherwise, antagonists like extracellularly secreted antagonist Sclerostin (SOST) was upregulated due to micro-topography at mRNA level. SOST binds to the LRP5/6 receptor and finally inactivates the Wnt/β-catenin pathway [23,27,56,73].
Two β-catenin signaling antagonists, SRY-box transcription factor 17 (SOX17) and inhibitor of β-catenin and TCF-4 (ICAT), were upregulated on P5 at mRNA level. Both antagonists are expressed and react only when β-catenin accumulates and can translocate into the nucleus, where β-catenin acts as a transcriptional co-activator [26]. ICAT binds to β-catenin in cytosol to promote the nuclear export [16]. In the nucleus, ICAT binds β-catenin to blocking the binding site (β-catenin central armadillo repeat [17,74]) for TCF4 and preventing the β-catenin/TCF transcriptional activities [14,75,76], which ultimately impairs differentiation [16,18,25,74,75]. Chen et al. [77], in their study with isolated primary Col2a1-ICAT transgenic chondrocytes, observed the inhibition of Wnt/β-catenin signaling due to ICAT overexpression, which finally delayed cell growth and maturation. Stow [78] published a review stating that the expression of the inhibitor ICAT correlates conversely with the transcriptional activity of β-catenin. Due to the higher level of the ICAT, the β-catenin is no longer able to act as a co-transcription promoter. The other negative regulator SOX17 repress the β-catenin dependent spatial and temporal LEF-1 activation [79]. Sinner et al. [26] demonstrated in their study that SOX17 interacts with both, TCF and β-catenin, and form a complex, which finally repress their transcriptional activities. Also SOX17 could be important for the regulation of the Wnt/β-catenin pathway and the manipulation of β-catenin-dependent transcriptional activation. An increased expression of both antagonists (ICAT and SOX17) in osteoblasts on P5 could also be confirmed at the protein level. As a result of the negative feedback loop, the gene expression could disturbed and thus the cell fate impaired [18,76,80].
The only significantly downregulated gene was secreted frizzled-related protein 4 (SFRP4), which interacts directly with Wnt, preventing binding to their receptor; in this way they antagonize the Wnt/β-catenin pathway [16,81]. Wang et al. [69] indicated in MG-63s on micro/nano surfaces a reduced expression of sFRP which induced enhanced nuclear amount of β-catenin. Dishevelled segment polarity protein3 (DVL3) prevents β-catenin degradation, which finally results in β-catenin accumulation in the cytosol [18,54,62]. DVL3 showed slightly higher mRNA expression levels in MG-63s on P5 (out of cutoff criteria: fold change 1.36, * p < 0.0152; Table S1). In summary, the microarray gene profile provided evidence of modified Wnt/β-catenin pathway in MG-63s on sharp-edged micro-pillars. It could be demonstrate on the one hand the upregulation of factors for Wnt/β-catenin pathway activation like WLS, TCF4, DVL. However, on the other hand an increased mRNA expression of β-catenin inhibitors such as SOST, CTNNBIP1 (ICAT) and SOX17 suggest that the β-catenin-dependent transcriptional activity (as co-activator) is affected due to micro-topography.
Various studies demonstrated the connection between redox and Wnt/β-catenin signaling [37][38][39][40]. The regulation of β-catenin accumulation can also be achieved by the production of reactive oxygen species (ROS) in cellular differentiation processes [29,37]. In an earlier study, internalization of micro-pillars in MG-63s was demonstrated by caveolae-mediated phagocytosis, a process associated with increased ATP turnover (n-fold to Ref: 0.627 ± 0.053, * p < 0.05) and mitochondrial activity [13]. Furthermore, this was associated with higher ROS production and higher expression of ROS-scavenging enzyme (Catalase, n-fold to Ref: 1.151 ± 0.089, * p < 0.05) [13]. The present study demonstrated the increased production of ROS in cells (cytosol and nucleus) and provides new evidence that the enhanced ROS level was detectable on P5 after only 1 h cell cultivation. The generation of ROS induced finally lipid peroxidation in MG-63s on P5, a biomarker for oxidative stress. In addition, the microarray data indicated an increased expression of SOD3 (fold change 1.32; Table S1). Hoogeboom et al. [58] identified DVL as an entry point for ROS to regulate oxidative stress. In their review they argued that under altered ROS conditions the pool of β-catenin is apparently diverted from TCF to Forkhead box-O (FOXO) transcription factors to regulate cell survival by increasing oxidative stress resistance and to support ROS degradation [58]. These antioxidative defense mechanism ultimately antagonizing the Wnt/β-catenin-dependent transcription [27,56]. In addition, Galli et al. [67] showed in their study how a SLA topography modulates the resistance of C2C12 cells to oxidative stress (stimulation with H 2 O 2 ) by inhibiting TCF-mediated transcription for activation of the FOXO/β-catenin pathway. The mRNA expression level of FOXO displayed only a very slight increase due to the topography (like FOXO1, fold change 1.04; FOXO6, fold change 1.22; Table S1). One further possibility for the β-catenin accumulation in MG-63s on P5 without β-catenin-dependent transcriptional activation could be the diversion of β-catenin for the FOXO/β-catenin pathway as a defense mechanism.

Redox-Dependent β-Catenin Accumulation
While the connection between these experimental measurements is not immediately obvious, it becomes apparent when considering the in silico model. As shown in [28,37], elevated ROS levels activate DVL in a redox-dependent and Wnt-independent manner. In its inactive state, DVL is primarily bound to Nucleoredoxin (NRX), a specific member of the thioredoxin-related redox-regulating protein family. The redox-sensitive association between NRX and DVL, however, is released upon elevated intracellular ROS levels [37]. This leads to a spontaneous increase of unbound DVL molecules. Interestingly, DVL proteins are capable of forming self-aggregates if their local concentration passes a certain threshold. DVL self-aggregates present a high-affinity binding platform for AXIN [29,37,82]. Notably, due to its low copy number, AXIN is a crucial component of the β-catenin degradation complex [16,24]. Its binding to DVL leads to an inhibition of β-catenin degradation, resulting in an accumulation of β-catenin in the cytosol and nucleus ( Figure 9A).  Table S1). One further possibility for the β-catenin accumulation in MG-63s on P5 without β-catenindependent transcriptional activation could be the diversion of β-catenin for the FOXO/β-catenin pathway as a defense mechanism.

Redox-Dependent β-Catenin Accumulation
While the connection between these experimental measurements is not immediately obvious, it becomes apparent when considering the in silico model. As shown in [28,37], elevated ROS levels activate DVL in a redox-dependent and Wnt-independent manner. In its inactive state, DVL is primarily bound to Nucleoredoxin (NRX), a specific member of the thioredoxin-related redoxregulating protein family. The redox-sensitive association between NRX and DVL, however, is released upon elevated intracellular ROS levels [37]. This leads to a spontaneous increase of unbound DVL molecules. Interestingly, DVL proteins are capable of forming self-aggregates if their local concentration passes a certain threshold. DVL self-aggregates present a high-affinity binding platform for AXIN [29,37,82]. Notably, due to its low copy number, AXIN is a crucial component of the β-catenin degradation complex [16,24]. Its binding to DVL leads to an inhibition of β-catenin degradation, resulting in an accumulation of β-catenin in the cytosol and nucleus ( Figure 9A). Figure 9. Simplified representation of (A) the signal pathway from reactive oxygen species (ROS) production (P5 induced) to eventual accumulation of β-catenin. Briefly, MG-63s respond to P5 with increased production of ROS, which inhibit AXIN-mediated degradation of β-catenin through Dishevelled (DVL). The width of the lines is proportional to the reaction kinetics (flux) through each reaction compared to Ref. Note that β-catenin accumulation typically induces increased AXIN synthesis. However, this is attenuated by the specific expression pattern, as shown in (B) the transcriptional regulation of β-catenin target gene expression on P5. A briefly elevated gene expression of T-cell specific transcription factor 4 (TCF4), inhibitor of β-catenin and TCF-4 (ICAT) and SRY-box transcription factor 17 (SOX17) leads to attenuated formation of the TCF/β-catenin complex (by ICAT) and proteosomal degradation of the β-catenin complex (by SOX17). The β-catenin target gene expression is thereby kept level despite its intracellular aggregation. The width of the lines is proportional to the reaction kinetics (flux) through each reaction compared to Ref.

Regulation of the Transcriptional Activity of β-Catenin
Typically, β-catenin accumulation in the nucleus leads to an elevated expression of WNT/βcatenin target genes, such as Myc, c-Jun and AXIN2. Notably, AXIN2 has a similar function as its homolog AXIN and induces the degradation of β-catenin. This means that nuclear β-catenin accumulation usually induces a negative feedback by increasing the expression of its own inhibitor ( Figure 9B). Surprisingly, the expression of β-catenin target genes is not increased on P5, despite the significant aggregation of β-catenin in the nucleus. This leads to two major questions that the model aims to answer: how do the inhibitors of the TCF/β-catenin complex such ICAT and SOX17 interact (quantitatively) with β-catenin to inhibit its transcriptional function; and at the same time, what keeps β-catenin aggregation in balance without the elevated transcription of its inhibitor AXIN2? To answer these questions, the role of the differentially expressed intracellular components in the model was analyzed, i.e., ICAT, SOX17 and TCF. Accordingly, ICAT is located in the cytosol as well as in the nucleus. Its inhibiting function is based on the competitive (reversible) binding of β-catenin [14,16]. SOX17, on the other hand, only acts inside the nucleus, directly binds to the TCF/β-catenin complex Figure 9. Simplified representation of (A) the signal pathway from reactive oxygen species (ROS) production (P5 induced) to eventual accumulation of β-catenin. Briefly, MG-63s respond to P5 with increased production of ROS, which inhibit AXIN-mediated degradation of β-catenin through Dishevelled (DVL). The width of the lines is proportional to the reaction kinetics (flux) through each reaction compared to Ref. Note that β-catenin accumulation typically induces increased AXIN synthesis. However, this is attenuated by the specific expression pattern, as shown in (B) the transcriptional regulation of β-catenin target gene expression on P5. A briefly elevated gene expression of T-cell specific transcription factor 4 (TCF4), inhibitor of β-catenin and TCF-4 (ICAT) and SRY-box transcription factor 17 (SOX17) leads to attenuated formation of the TCF/β-catenin complex (by ICAT) and proteosomal degradation of the β-catenin complex (by SOX17). The β-catenin target gene expression is thereby kept level despite its intracellular aggregation. The width of the lines is proportional to the reaction kinetics (flux) through each reaction compared to Ref.

Regulation of the Transcriptional Activity of β-Catenin
Typically, β-catenin accumulation in the nucleus leads to an elevated expression of WNT/β-catenin target genes, such as Myc, c-Jun and AXIN2. Notably, AXIN2 has a similar function as its homolog AXIN and induces the degradation of β-catenin. This means that nuclear β-catenin accumulation usually induces a negative feedback by increasing the expression of its own inhibitor ( Figure 9B). Surprisingly, the expression of β-catenin target genes is not increased on P5, despite the significant aggregation of β-catenin in the nucleus. This leads to two major questions that the model aims to answer: how do the inhibitors of the TCF/β-catenin complex such ICAT and SOX17 interact (quantitatively) with β-catenin to inhibit its transcriptional function; and at the same time, what keeps β-catenin aggregation in balance without the elevated transcription of its inhibitor AXIN2? To answer these questions, the role of the differentially expressed intracellular components in the model was analyzed, i.e., ICAT, SOX17 and TCF. Accordingly, ICAT is located in the cytosol as well as in the nucleus. Its inhibiting function is based on the competitive (reversible) binding of β-catenin [14,16]. SOX17, on the other hand, only acts inside the nucleus, directly binds to the TCF/β-catenin complex and, in contrast to ICAT, induces the degradation of the TCF/β-catenin complex [26]. Note that this degradation is independent of the cytosolic destruction complex, hence independent of AXIN and GSK-3β. The difference between their localization and the inhibiting function of ICAT and SOX17 has major implications on their impact on β-catenin aggregation and transcriptional activity. The competitive binding of ICAT to β-catenin, for instance, reduces the amount of β-catenin available for its binding to TCF, hence attenuating the transcription; at the same time, it allows an accumulation of β-catenin because it does not induce the degradation of β-catenin. However, for a sustained effect (i.e., attenuating the transcriptional activity of β-catenin) the amount of ICAT needs to match the number of β-catenin. Since β-catenin is highly abundant in the cell, a very high quantity of ICAT would be required to attenuate target gene expression. Here SOX17 acts much more efficiently: since SOX17 targets the entire TCF/β-catenin for proteasomal degradation, both β-catenin and TCF are degraded. While β-catenin is highly abundant, TCF is typically available in a much lower quantity. As a result, even a relatively small increase in the quantity of SOX17 can significantly lower the transcription and expression of β-catenin target genes. However, such an increase of SOX17 needs to be accompanied by a similarly increased expression rate of TCF; otherwise, the system reaches a point where no TCF is available and target gene expression is completely inhibited.

Assembling the Entire View
Relating the insights gained from the model to the in vitro data obtained in MG-63s after cultivation on P5 provides a comprehensive picture of how, in contrast to Ref, β-catenin upregulation and translocation with constant/unchanged transcriptional activity is achieved. At first, the increased and constant production of ROS on P5 ( Figure 9) induces a redox-dependent accumulation of β-catenin. To inhibit β-catenin target gene activation, such as Myc, c-Jun or AXIN2, MG-63s express the specific inhibitors ICAT and SOX17 (with similar fold changes; Table 1) that act concurrently. While ICAT prevents a major fraction of β-catenin from binding to TCF, SOX17 targets the remaining, successfully established TCF/β-catenin complexes for proteasomal degradation. However, the cells have to encounter the SOX17-mediated degradation of TCF to prevent a decrease or total loss of transcriptional activity. Notably, the expression of β-catenin target genes is unchanged between P5 and Ref. Therefore, TCF4 expression is also elevated with a similar fold change as SOX17. Based on this coordinated expression of regulatory elements of the Wnt/β-catenin pathway, MG-63s are able to cope with the ROS-induced accumulation of β-catenin and suppress an unintended target gene expression.

In Summary
Defined sharp-edged micro-pillars caused a change in cell morphology ( Figure S1), reduced the calcium (Ca 2+ ) mobilization and induced a cellular stress reaction with low adenosine triphosphate (ATP) levels and increased ROS production [13]. In the present study it was shown that the Wnt/β-catenin pathway was regulated by the micro-topography-β-catenin accumulation and translocation into the nucleus ( Figure 10). The ultimately impairment of osteogenesis [9,12,13] could be caused by an inhibition of the TCF/β-catenin transcriptional activity due the antagonists ICAT and SOX17 or via other pathways involved in the regulation of β-catenin as co-activator (FOXO/β-catenin signaling pathway) [14,20,68]. However, possible mechanisms how the cell fate is then influenced by the regulation of transcriptional co-activity, was revealed by the in silico study, and confirmed by the increased expression of the inhibitors ICAT and SOX17 as well as TCF4. In addition, the higher cellular oxidative stress on P5 initiates the accumulation of β-catenin and subsequent redirection into other signaling pathways, such as the defense mechanism against ROS.
Other antagonistic effects on Wnt signaling in MG-63 osteoblasts, such as hypoxia [83], treatment with inhibitors of Wnt/β-catenin-mediated transcription PRI-724 [84] and Resveratrol [85] could found in other experiments. Upcoming studies are needed to understand the downstream mechanism of molecular events that control cell behavior induced by micro-topography. In this context, topography properties such as dimension, shape, or smoothness should be modulated, and then correlated with their effects on the Wnt/β-catenin pathway. Cell-stimulatory or inhibitory effects of topographical surface properties of a biomaterial can be studied further in silico to better understand cell behavior. It is important to elucidate possible molecular mechanisms for the design of metal biomaterials in regenerative medicine with improved acceptance and osseointegration.
1 Figure 10. Summary of the complex intracellular responses of MG-63s on micro-pillars which all contribute to reduced cell functions. The sharp-edged micro-pillars (P5) changed the cell morphology (FE-SEM Supra25, Carl Zeiss; Figure S1) [9,13] and hampered the calcium mobilization [9]. Furthermore, a reduced adenosine triphosphate (ATP)-content [13] as well as an increased reactive oxygen species (ROS) production due to cellular stress could be observed ( Figure 6) [13]. The protein-level of β-catenin as well as the translocation into the nucleus was increased (Figure 1), but antagonists such as ICAT (inhibitor of β-catenin and TCF-4) (Figure 4) finally lead to an inhibited transcriptional activation. As a consequence, osteoblasts were reduced in their extracellular matrix-synthesis like collagen type I ( Figure 5) [9,12,13]. (Scale bars 10 µm).