Next Article in Journal
Equine Tenocyte Seeding on Gelatin Hydrogels Improves Elongated Morphology
Next Article in Special Issue
Development of Conductive Gelatine-Methacrylate Inks for Two-Photon Polymerisation
Previous Article in Journal
Carbon Nanomaterials Embedded in Conductive Polymers: A State of the Art
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Triborheological Study under Physiological Conditions of PVA Hydrogel/HA Lubricant as Synthetic System for Soft Tissue Replacement

by
Laura C. Duque-Ossa
,
Gustavo Ruiz-Pulido
and
Dora I. Medina
*
Tecnologico de Monterrey, School of Engineering and Science, Atizapan de Zaragoza, Estado de Mexico 52926, Mexico
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(5), 746; https://doi.org/10.3390/polym13050746
Submission received: 20 January 2021 / Revised: 11 February 2021 / Accepted: 13 February 2021 / Published: 28 February 2021

Abstract

:
In soft tissue replacement, hydrophilic, flexible, and biocompatible materials are used to reduce wear and coefficient of friction. This study aims to develop and evaluate a solid/liquid triborheological system, polyvinyl alcohol (PVA)/hyaluronic acid (HA), to mimic conditions in human synovial joints. Hydrogel specimens prepared via the freeze–thawing technique from a 10% (w/v) PVA aqueous solution were cut into disc shapes (5 ± 0.5 mm thickness). Compression tests of PVA hydrogels presented a Young’s modulus of 2.26 ± 0.52 MPa. Friction tests were performed on a Discovery Hybrid Rheometer DHR-3 under physiological conditions using 4 mg/mL HA solution as lubricant at 37 °C. Contact force was applied between 1 and 20 N, highlighting a coefficient of friction change of 0.11 to 0.31 between lubricated and dry states at 3 N load (angular velocity: 40 rad/s). Thermal behavior was evaluated by differential scanning calorimetry (DSC) in the range of 25–250 °C (5 °C/min rate), showing an endothermic behavior with a melting temperature (Tm) around 231.15 °C. Scanning Electron Microscopy (SEM) tests showed a microporous network that enhanced water content absorption to 82.99 ± 1.5%. Hydrogel achieved solid/liquid lubrication, exhibiting a trapped lubricant pool that supported loads, keeping low coefficient of friction during lubricated tests. In dry tests, interstitial water evaporates continuously without countering sliding movement friction.

1. Impact Statement

Musculoskeletal diseases represent one of the main causes of labor and academic absence in the world, so they represent a public health problem that causes economic, time, and production losses [1]. Osteoarthritis is the most common degenerative disease associated with synovial joints, cartilage and bone wear, bone volume loss and deterioration, osteophyte formation, and inflammation [2,3]. Around the world, it is present in more than 25% of the adult population [2]. However, most treatments are ineffective over the long term, and thus, cartilage replacement surgery is commonly used for treating the disease. The procedure involves the use of synthetic elements which emulate soft tissues and synovial joint behavior. In order to determine components that can be used for treatments, this study proposes the use of biocompatible triborheological systems under physiological conditions, operating with synthetic materials with non-isotropic surfaces in terms of roughness.

2. Introduction

The human body suffers degradation and loss of motor skills over time, affecting quality of life in relation to the mental and physical health of individuals in areas of functional independence [4,5]. The joint cavity is composed by cartilaginous tissue, synovial fluid, subchondral bone, and sometimes menisci [3]. The synovial fluid (SF) is a biological lubricant composed by mucopolysaccharides, proteins, and lipids, which form a protective lubricant film for articular joints by reducing friction and wear [6,7,8]. The articular cartilage (AC) creates a fibrous connection between 2 or more bones [9,10]. It is integrated by chondrocytes and an extracellular matrix (ECM), which give strength and help to generate good mechanical properties against stress loads [11,12].
When an articulation begins its failure stage, it means that its AC has lost its properties, decreasing its volume and deteriorating [10,13,14]. This leads to the appearance of diseases such as osteoarthritis (OA) [15,16]. OA is a degenerative disease associated with wear and tear of joints such as the knee, hip, foot, hand, spine, and elbow, among others (Figure 1), presenting loss of AC, osteophyte formation, inflammation, alteration of SF composition, subchondral bone sclerosis, and varying degrees of synovitis [2,3,17,18]. Worldwide, it is present in more than 25% of the population over 18 years of age and is responsible for an economic cost of around 89.1 billion USD per year worldwide. [2,17]. OA pathology is derived not only from age, but from lifestyle, weight, overuse, injury, and genetic and gender conditions [2,19,20]. For example, cervical facet joint OA affects older adults, increasing its incidence with age [21,22]. Its prevalence is around 54–67% of spinal pain symptoms [23,24]. For hip OA, the population over 85 years old presents a risk of around 25% of suffering from the disease. [25]. In the case of knees, it is estimated that 6% of adults worldwide have signs of OA [26].
Medical and engineering researchers have tried to reduce the progress of OA and correct its effects with pharmacological interventions, physical therapy, cell therapy, and tissue engineering through surgical methods [3,27,28]. AC as an avascular tissue does not regenerate easily, so surgical methods are most commonly used for restoring or replacing cartilage surfaces [10,29]. Tissue engineering has generated advances by repairing cartilage and developing implants for surgery, considering human body mechanics and its response to external agents [30,31]. Cell culture techniques, clinical trials, scaffolds, and synthetic materials are also used [30,32].
The use of biomaterials has focused on polymers, mainly hydrogels [10]. They have hydrophilic structures (high water content), softness, flexibility, and biocompatibility [33,34,35]. Synthetic hydrogels can mimic AC behavior in terms of lubrication and mechanical performance. They can even improve upon AC’s toughness and strength, providing extraordinary mechanical properties [36,37,38]. Within the polymer group, polyvinyl alcohol (PVA) is the most relevant in current studies [39,40,41]. PVA has been used to improve and to regenerate different tissues and organs in the human body, such as arterial phantoms, corneal implants, and cartilage tissue replacements [41]. Specifically, for AC substitution, several sliding studies have demonstrated its performance under diverse temperatures, charges, and velocity conditions. Thus, PVA is considered an excellent choice in tissue engineering for friction reduction and low wear results [42,43,44].
For tests, the use of lubricants to minimize damage in artificial cartilage is imperative. Hyaluronic acid (HA), a component of synovial fluid, is used as an alternative to increase joint lubrication due to its viscoelastic properties. [45,46,47]. This viscous fluid provides protection against articular wear [48,49]. Today, HA is commonly used in tribological tests that involve the study of PVA hydrogels [50].
Due to the above, the goal of this research was to evaluate the potential of PVA hydrogels for soft tissue replacement, such as AC, using HA as lubricant, under physiological conditions of sample temperature (37 °C) and using a saline solution as a solvent for sample preparation to recreate fluid conditions in the human body. Additionally, pressures applied during the tests were equivalent to the impact upon joints on a daily basis.

3. Materials and Methods

3.1. Preparation of PVA Hydrogels

The specimens were produced from a 10% PVA aqueous solution (Sigma Aldrich, St. Louis, MO, USA, 13 × 104 Daltons, 99% hydrolyzed, 10% w/w) in distilled water by the freeze–thawing technique (FT) [51]. The mixture was stirred continuously at 80 °C for 3 h using a magnetic stirrer to achieve complete homogenization. The solution was placed in glass petri dishes and four FT cycles were performed. The freezing stage was carried out for 12 h at −20 °C, while the thawing stage was carried out at room temperature for 12 h.
For triborheological and mechanical tests, FT hydrogel samples were cut in a disc shape with 5 ± 0.5 mm thickness and 45 mm diameter.
HA lubricant solution was prepared using 4 mg of HA powder (Lifecore Biomedical, Chaska, MN, USA, 8 × 105 Daltons) per mL of saline solution (0.15 M NaCl) under stirring conditions of 150 rpm for 6 h.

3.2. Equilibrium Water Content Measurements

For equilibrium water content measurement (EWCM), the procedure proposed by Sardinha et al. in 2013 was followed. PVA hydrogels samples were cut into 25 × 25 mm squares with 5 ± 0.5 mm thickness, weighted, and submerged in distilled water at ambient temperature until they reached a stable weight. Samples were removed from water and weighted again. Weight difference between the hydrated hydrogel (Wh) and the dehydrated hydrogel (Wd) gives the EWCM in percentage, using Equation (1) [44].
E W C M   ( % ) =   W h W d W h × 100

3.3. Mechanical Tests

Compression tests were carried out, emulating joint cavity conditions in a Shidmadzu UH universal testing machine (Shimadzu Corporation, Kyoto, Japan): humidity, preheating at 37 °C, compression speed of 1 mm/min with a maximum load of 4500 N. Considering Equations (2)–(4), and values obtained from stress vs. strain curves, the parameters were calculated as described by Gupta et al. [52,53]:
σ =   F A
ε =   δ L 0 =   L 0 L f L 0
E =   σ ε
where σ is the stress, F is the applied load, A is the area where the load is applied, ε is the strain, δ is the calibrated length, L0 is the initial length, Lf is the final length, and E is the modulus of elasticity or Young’s modulus [52].

3.4. Rheological and Friction Tests

Tests were performed with a Discovery Hybrid Rheometer DHR-3 (TA Instruments, New Castle, DE, USA, maximum charge 50 N) using cone-and-plate (0.9969°) configuration for rheological characterization (Figure 2a) and plate-and-plate configuration for the triborheological part (Figure 2b).
All the experiments were performed at 37 °C to mimic the physiological conditions of the human body. HA in 4 mg/mL concentration in saline solution 0.10 M NaCl was used as lubricant fluid [54]. Samples were submerged at all times in lubricant to ensure their complete hydration. Velocity was fixed between 0.1 and 100 rad/s, and contact force was applied in a range of 1–20 N. The coefficient of friction changes were continuously monitored during tests.
Rheology tests were carried out in 3 steps: conditioning sample (37 °C during 10 s), frequency sweep (0.01–100 Hz, 10% strain), and flow sweep (shear rate 0.001–100 s−1). Triborheological tests were conducted in 2 steps: conditioning options active (37 °C for all charges: 1 N, 3 N, 5 N, 10 N, 15 N, 20 N) and flow sweep (angular velocity 0.1–100 rad/s).

3.5. Surface Analysis

The morphology of the PVA hydrogel and identification of wear mechanism traces was observed by JEOLJSM-6360 scanning electron microscopy (SEM, JEOL USA, Peabody, MA, United States) at 20 kV. Samples were managed using high vacuum at level 3 of the beam diameter, working with secondary electrons at magnifications between 50 and 1000×. Samples were cut into squares of 10 × 10 mm (thickness of 5 ± 0.5 mm) and gold-sputtered to enhance conductivity. Samples were analyzed in surface distribution, porosity, and wear tracks before and after triborheological tests.

3.6. Thermal Behavior

The thermal properties of dried hydrogels were measured by differential scanning calorimetry (DSC, Q10 TA Instruments, United States). Dried hydrogels (2–5 mg) were placed in an aluminum pan and heated at a rate of 5 °C/min from 25–250 °C in a nitrogen atmosphere. The melting enthalpy, ΔHm, was determined by integrating the area under the melting peak over the range 208–242 °C. The crystallinity (Xc) was calculated by
X c = Δ H m Δ H c × 100 %
where ΔHc = 138.6 J/g, which represents the heat required for melting a 100% crystalline PVA sample [55].

4. Results and Discussion

4.1. Surface Morphology and Water Content

Naturally, AC has a high water content to mitigate the stress generated from contact forces applied during the movement of synovial joints. Hydrogels must present similar swelling properties to correctly mimic the natural AC performance in terms of super water-absorbing capability, elasticity, compressive mechanical properties, wear reduction, and low coefficient of friction during sliding [44,56,57]. For that reason, the water content absorption of PVA hydrogels was analyzed by measuring the equilibrium water content in percentage (EWCM %). Results showed an average EWCM of 82.99 ± 1.5% in accordance with Table 1.
The percentage of water stored in the hydrogel is not dependent on the initial weight of the sample. Even thickness variation does not affect the swelling range, which is maintained around the same value. The swelling behavior of hydrogels is a consequence of the diffusion of water molecules inside the polymer network [58], which depends mainly on the sample’s porosity and porous size, which were similar in all the samples as all were prepared under the same conditions [57]. Therefore, the obtained swelling values for the material are adequate for its purpose, because the water content of a natural CA is around 80% [59,60].
Unworn surfaces of PVA hydrogel were visualized in SEM equipment for surface distribution and porosity evaluations. According to Figure 3a, the surface was smooth and largely homogeneous despite the presence of inclusions adhered from the environment (dust particles) that adhered during the process of dehydration. Moreover, SEM micrographs of PVA hydrogels showed the porous network that results from the spaces generated by the ice crystals formed during the FT process [56]. Those pores allowed water absorption, reducing hydrogel wear (Figure 3b). Similar results have been reported by Sardinha et al. and Li et al. in studies of PVA hydrogels under similar conditions [43,44].

4.2. Mechanical Tests

The stress–strain curves obtained from 7 tests of compression after 4 FT cycles are reported in Figure 4a. All samples exhibited a non-linear behavior and were plotted from the data obtained by the instrument. Using a linear regression, the average Young’s modulus obtained was 2.26 ± 0.52 MPa. The displacement observed was caused by the standard deviation present in the initial thickness of the samples.
In accordance with the curves, the behavior of the hydrogel samples accurately fits an exponential relationship. For example, test 1 exhibited a correlation coefficient R2 higher than 0.9, indicating a positive proportionality of the effort against deformation, which explains the viscoelastic characteristics of the material (flexibility and rapid response to external loads) [61].
Figure 4b shows that when the deformation rises from 33.7% to 70%, Young’s modulus increments from 0.05 MPa to 2.7 MPa, which represents an increase of 51 times. This mechanical behavior is expected, and comparable to natural AC. When the cartilage is subjected to a considerable high load activity (such as running or jumping), the compression module increases to resist high stress, reducing deformation and preventing cartilage damage [61].
Measurement of Young’s modulus exhibited similar results to previous studies that involved mechanical characteristics of natural and synthetic cartilages. The viscoelastic behavior of the samples during the tests is shown in Figure 4c, through images that reveal the deformation (sample 5) at certain times.

4.3. Rheological and Friction Tests

Rheological test data are displayed in Figure 5. HA presented a non-Newtonian behavior, where the viscosity decreases (Figure 5a) as the shear stress increases (Figure 5b), in a clear example of shear thinning [62].
The evaluation of complex viscosity (η*), elastic (G′), and viscous (G″) modulus (Figure 5c) showed that at low frequencies, the sample exhibited an elastic behavior (G′ > G″), while at high frequencies, it presented a viscous behavior (G″ > G′). This behavior is typical of substances in solution, known as “weak gels”, and is contrary to the results of pure HA, which initially behaves viscously (G″ > G′) [54,63,64,65]. Additionally, the number of FT cycles should increase the elastic modulus, which is related to the gel strength and porosity, because the cycles stimulate the formation of crystalline junction zones among PVA chains and the appearance of hydrogen bonds resulting in PVA densification. As more couplings and bonds are formed, the sol–gel transition occurs and the sample no longer flows as a polymer solution, and its viscoelastic properties exhibit a solid-like behavior as its macromolecules lose their identity and become part of a large three-dimensional network, expanding throughout the entire sample volume [58].
The crossover frequency from a viscous to an elastic state demonstrates the viscoelastic behavior expected from a lubricant. The changing behavior of the modulus indicates an increased deformation of the cartilage at high frequencies, enhancing its resistant properties [66].
Data obtained of the coefficient of friction of PVA synthetic cartilage against stainless steel plates under different factors are shown in Figure 6. Figure 6a,b indicates that coefficient of friction is affected by angular velocity, increasing between 0 and 4 rad/s, followed by a quick decrease among 4–10 rad/s. Finally, it rises slowly again after 10 rad/s. The speed generates the same effect on the sample, regardless of whether it is lubricated or not. Specifically, coefficient of friction increased from 0.19 to 0.3 when the velocity incremented from 0 to 4 rad/s and decreased from 0.3 to 0.14 when velocity moved from 4–10 rad/s under 5 N load.
As seen in Figure 6c, lubrication causes a decrease in the coefficient of friction. For a 3 N load and an angular velocity of 40 rad/s, it goes from 0.11 to 0.31 when switching from HA lubrication to a dry state. In a dry state, the large amount of water accumulated by the samples contributes to counteracting the sliding movement, while in a lubricated state, the transition of the lubricant from elastic to a viscous state minimizes the friction generated at high angular velocity values.
Pressure also generates an effect upon coefficient of friction. For 0.46 PSI with HA lubricant, the coefficient of friction was 0.14 at an angular velocity of 10 rad/s. For the same conditions, the coefficient of friction increased to 0.24 under a 1.82 PSI load as Figure 6d shows. Meanwhile, in dry samples, the coefficient of friction increased exponentially after 0.91 PSI.
Under an applied load, the coefficient of friction decreases, followed by an increase as the water or HA (depending on the conditions of the experiment) is squeezed out from the PVA hydrogel through its porous network. As the interstitial fluid expands from the porous structure to the conjunction and suffers instantaneous pressurization supporting a fraction of the applied load, it produces a coefficient of friction reduction.
This behavior could be explained by comparing 3 parts of the lubrication mechanisms of natural articular cartilage. Boundary lubrication is an important lubrication mechanism in natural synovial joints, forming a molecular film [45]. This film takes up the boundary lubricants (HA, lubricin, among others), providing a defense against contact [67]. PVA hydrogels do not present a molecular layer that would recreate boundary lubrication similar to synovial joints, but the initial administration of lubricant prior to each test mimics this function to an extent. HA, at the beginning of the test, forms a small film that avoids hard contact between PVA hydrogel and stainless-steel plate, diminishing surface wear. This film is provided by a lubricant pool (trapped lubricant) that gets out from the hydrogel, creating an interface between the steel and the hydrogel that reduces the effective contact area between surfaces, and consequently the coefficient of friction (Figure 7a) [68]. Moreover, the liquid layer could be continuously fed by the trapped and external lubricant through hydrogel swelling to maintain its pressurization and to form new layers during sliding [69]. Interstitial fluid rehydration is limited by hydrogel swelling capacity.
After a while, during contact between both surfaces, only a small portion of the rough parts (nano asperities) of the stainless-steel plate interact directly, producing friction. This is a consequence of the presence of mixed lubrication, where the HA supports a fraction of the load. However, continuous increase in the load causes the displacement of the lubricant from the hydrogel surface, enlarging the contact area, which leads to coefficient of friction increase [70]. This behavior is known as hydrodynamic lubrication. In this case, high pressurization of the interstitial fluids occurs under applied normal loads as the fluid supports most of the loads transmitted through articular surfaces [71]. As such, friction force in the contact surface is reduced considerably, ensuring that coefficient of friction maintains low values. By contrast, under dry sliding conditions, the liquid inside the porous structure expands as lubricant during the initial period. Afterward, the superficial liquid evaporates without replacing, which produces an increase in the coefficient of friction and surface wear. In natural articular cartilage, mixed and hydrodynamic behavior are explained by solid/liquid theory that considers the collagen–proteoglycan network as a solid phase interacting with interstitial liquid and ions as a fluidic phase [72]. In Figure 7b, all mechanisms applied to the solid/liquid system form by PVA hydrogel and HA as a lubricant can be observed.
In general, all results reported are comparable to previous multifactor studies. Pan et al. and Shi et al. showed that an increment in the coefficient of friction is related to the increasing load, as well as with the effect of lubrication or dryness of the sample [72,73]. Likewise, values reported for coefficient of friction in natural AC and synthetic PVA AC are similar. According to Kobayashi et al., natural AC has a coefficient of around 0.1 and synthetic PVA cartilage made by FT cycles between 0.7 and 0.8 [74]. Thus, PVA hydrogel shows a higher coefficient of friction than natural articular cartilage.

4.4. Thermal Behavior

The DSC melting curves of dried PVA hydrogels of different weights are shown in Figure 8. Peaks reported at around 232.31 °C represent the approximate melting temperature, Tm, of the samples. The values obtained are related to sample crystallinity, as has been expressed in previous studies [55]. The melting enthalpy and crystallinity percentages of PVA hydrogel samples were calculated by analyzing the area under the curve from the melting peaks using Universal Analysis 2000 (TA Instruments) software, as shown in Table 2.
In general, PVA hydrogels prepared by the FT technique exhibited similar crystallinity degrees to values reported in literature [75]. Crystallinity can be affected by factors like FT cycles, concentration of aqueous solution, and molecular weight of PVA [55,76].

5. Conclusions

The mechanical and triborheological properties of PVA synthetic cartilage/HA lubricant were investigated, and they demonstrated attributes for the generation of potential prototypes for AC replacement by testing, under physiological conditions, the effects of swelling, load, angular velocity, and lubrication. Surface morphologies of the samples were analyzed with SEM and mechanical properties were determined by a universal machine. PVA showed an average EWCM of 82.99 ± 1.5%, which is comparable to natural AC’s EWCM of 80%. Additionally, the compression Young’s modulus was higher than 2, providing good resistance to stress by reducing deformation and preventing cartilage damage. In rheology studies, the crossover frequency influenced the behavior of the lubricant, changing its nature from elastic state to a more viscous form and indicating a good performance.
It was also confirmed, through triborheological tests, that the coefficient of friction exhibited values between 0.1 and 0.4, which are comparable with those of the coefficient of natural AC that is around 0.1. Therefore, 10% PVA hydrogel lubricated with 4 mg/mL HA could represent a great mechanical alternative for future advances in replacing damaged AC in human joints under moderate load.

Author Contributions

Conceptualization, L.C.D.-O. and D.I.M.; methodology, L.C.D.-O. and D.I.M.; validation, L.C.D.-O. and G.R.-P.; formal analysis, L.C.D.-O., G.R.-P. and D.I.M.; investigation, L.C.D.-O., G.R.-P. and D.I.M.; data curation, L.C.D.-O.; writing—original draft preparation, L.C.D.-O. and G.R.-P.; writing—review and editing, D.I.M.; supervision, D.I.M.; project administration, D.I.M.; funding acquisition, D.I.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Research supported by Tecnologico de Monterrey and Mexican National Council for Science and Technology (CONACYT). The authors would like to make a most cordial acknowledgment to the research group of Dora Medina, Carlos Rojo, Ranieri Durante, and Jose Antonio Garcia for their support during the research.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. OMS. Prevención de Trastornos Musculoesqueléticos en el Lugar de Trabajo: Serie Protección de la Salud de los Trabajadores n 5; 2004. Available online: https://www.who.int/occupational_health/publications/muscdisorders/es/ (accessed on 1 February 2021).
  2. Chen, D.; Shen, J.; Zhao, W.; Wang, T.; Han, L.; Hamilton, J.L.; Im, H.-J. Osteoarthritis: Toward a comprehensive understanding of pathological mechanism. Bone Res. 2017, 5, 16044. [Google Scholar] [CrossRef]
  3. Clouet, J.; Vinatier, C.; Merceron, C.; Potvaucel, M.; Maugars, Y.; Weiss, P.; Grimandi, G.; Ghicheux, J. From osteoarthritis treatments to future regenerative therapies for cartilage. Drug Descov. Today 2009, 14, 913–925. [Google Scholar] [CrossRef]
  4. Mahir, L.; Belhaj, K.; Zahi, S.; Zanmasso, H.; Lmidmani, F.; el Fatimi, A. Impact of knee osteoarthritis on the quality of life. Ann. Phys. Rehabil. Med. 2016, 59, 155–159. [Google Scholar] [CrossRef]
  5. Kumar, H.N.H.; Nagaraj, K.; Luthra, K.; Gupta, P.; Sapar, P.; Gupta, S.; Tyagi, A. Health-related quality of life among osteoarthritis patients attending primary care clinics of mangalore city. Int. J. Med. Public Health 2015, 5, 55–58. [Google Scholar] [CrossRef] [Green Version]
  6. Stevenson, M.; Jaggard, H.A.; Akhbari, U.; Vaghela, P.; Gupte, C.; Cann, P. The role of denatured synovial fluid proteins in the lubrication of artificial joints. Biotribology 2019, 17, 49–63. [Google Scholar] [CrossRef]
  7. Al-Shakarchi, I.; Coakley, G. Synovial fluid tests. Medicine 2018, 46, 166–169. [Google Scholar] [CrossRef]
  8. Horibata, S.; Yarimitsu, S.; Fujie, H. Effect of synovial fluid pressurization on the biphasic lubrication property of articular cartilage. Biotribology 2019, 19, 100098. [Google Scholar] [CrossRef]
  9. Fox, A.J.S.; Bedi, A.; Rodeo, S.A. The basic science of articular cartilage: Structure, composition, and function. Sports Health 2009, 1, 461–468. [Google Scholar] [CrossRef]
  10. Liao, I.-C.; Moutos, F.T.; Estes, B.T.; Zhao, X.; Guilak, F. Composite three-dimensional woven scaffolds with interpenetrating network hydrogels to create functional synthetic articular cartilage. Adv. Funct. Mater. 2013, 23, 5833–5839. [Google Scholar] [CrossRef]
  11. Vinatier, C.; Mrugala, D.; Jorgensen, C.; Ghicheux, J.; Noel, D. Cartilage engineering: A crucial combination of cells, biomaterials and biofactors. Trends Biotechnol. 2009, 27, 307–314. [Google Scholar] [CrossRef]
  12. Fritz, J.; Janssen, P.; Gaissmaier, C.; Schewe, B.; Weise, K. Articular cartilage defects in the knee basics, therapies and results. Injury 2008, 39, 50–57. [Google Scholar] [CrossRef] [PubMed]
  13. Grenier, S.; Bhargava, M.M.; Torzilli, P.A. An in vitro model for the pathological degradation of articular cartilage in osteoarthritis. J. Biomech. 2014, 47, 645–652. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Henao-Murillo, L.; Ito, K.; van Donkelaar, C.C. Collagen damage location in articular cartilage differs if damage is caused by excessive loading magnitude or rate. Ann. Biomed. Eng. 2018, 46, 605–615. [Google Scholar] [CrossRef] [Green Version]
  15. Pap, T.; Korb-pap, A. Cartilage damage in osteoarthritis and rheumatoid arthritis—Two unequal siblings. Nat. Rev. Rheumatol. 2015, 11, 606–615. [Google Scholar] [CrossRef] [PubMed]
  16. Stender, M.E.; Carpenter, R.D.; Regueiro, R.A.; Ferguson, V.L. An evolutionary model of osteoarthritis including articular cartilage damage, and bone remodeling in a computational study. J. Biomech. 2016, 49, 3502–3508. [Google Scholar] [CrossRef] [PubMed]
  17. Carlson, A.; Rawle, R.; Wallace, C.; Brooks, E.; Adams, E.; Greenwood, M.; Olmer, M.; Lotz, M.; Bothner, B.; June, R. Characterization of synovial fluid metabolomic phenotypes of cartilage morphological changes associated with osteoarthritis. Osteoarthr. Cartil. 2019, 27, 1174–1184. [Google Scholar] [CrossRef]
  18. Loeser, R.F.; Goldring, S.R.; Scanzello, C.R.; Goldring, M.B. Osteoarthritis: A disease of the joint as an organ. Arthritis Rheum. 2012, 64, 1697–1707. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, Y.; Jordan, J.M. Epidemiology of osteoarthritis. Clin. Geriatr. Med. 2010, 26, 355–369. [Google Scholar] [CrossRef] [Green Version]
  20. Loeser, R.F. Aging and osteoarthritis. Curr. Opin. Rheumatol. 2011, 23, 492–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Park, W.M.; Kim, K.; Kim, Y.H. Changes in range of motion, intradiscal pressure, and facet joint force after intervertebral disc and facet joint degeneration in the cervical spine. J. Mech. Sci. Technol. 2015, 29, 3031–3038. [Google Scholar] [CrossRef]
  22. Gellhorn, A.C.; Katz, J.N.; Suri, P. Osteoarthritis of the spine: The facet joints. Nat. Rev. Rheumatol. 2013, 9, 216–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. del Seguro Social, I.M. Abordaje Diagnóstico del Dolor de Cuello en la Población Adulta en el Primer Nivel de Atención; IMSS: Mexico City, Mexico, 2013. [Google Scholar]
  24. Manchikanti, L.; Boswell, M.V.; Singh, V.; Pampati, V.; Damron, K.S.; Beyer, C.D. Prevalence of facet joint pain in chronic spinal pain of cervical, thoracic, and lumbar regions. BMC Musculoskelet. Disord. 2004, 5, 15. [Google Scholar] [CrossRef]
  25. Murphy, N.J.; Eyles, J.P.; Hunter, D.J. Hip osteoarthritis: Etiopathogenesis and implications for management. Adv. Ther. 2016, 33, 1921–1946. [Google Scholar] [CrossRef] [Green Version]
  26. Wood, A.M.; Brock, T.M.; Heil, K.; Holmes, R.; Weusten, A. A review on the management of hip and knee osteoarthritis. Int. J. Chronic Dis. 2013, 2013. [Google Scholar] [CrossRef] [Green Version]
  27. Vincent, T.L.; Watt, F.E. Osteoarthritis. Medicine 2017, 46, 187–195. [Google Scholar] [CrossRef]
  28. Conaghan, P.G.; Dickson, J.; Grant, R.L. Care and management of osteoarthritis in adults: Summary of nice guidance. BMJ 2008, 336, 502–503. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Moffat, K.L.; Goon, K.; Moutos, F.T.; Estes, B.T.; Oswald, S.J.; Zhao, X.; Guilak, F. Composite cellularized structures created from an interpenetrating polymer network hydrogel reinforced by a 3d woven scaffold. Macromol. Biosci. 2018, 18, 1800140. [Google Scholar] [CrossRef]
  30. Fisher, M.B.; Mauck, R.L. Tissue engineering and regenerative medicine: Recent innovations and the transition to translation. Tissue Eng. Part B Rev. 2013, 19, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Coburna, J.M.; Gibsonb, M.; Monagleb, S.; Pattersonb, Z.; Elisseeff, J.H. Bioinspired nanofibers support chondrogenesis for articular cartilage repair. Proc. Natl. Acad. Sci. USA 2012, 109, 10012–10017. [Google Scholar] [CrossRef] [Green Version]
  32. Lysaght, M.J.; Jaklenec, A.; Deweerd, E. Great expectations: Private sector activity in tissue engineering, regenerative medicine, and stem cell therapeutics. Tissue Eng. Part A 2008, 14, 305–315. [Google Scholar] [CrossRef]
  33. Ahmed, E.M. Hydrogel: Preparation, characterization, and applications: A review. J. Adv. Res. 2015, 6, 105–121. [Google Scholar] [CrossRef] [Green Version]
  34. Calo, E.; Khutoryanskiy, V.V. Biomedical applications of hydrogels: A review of patents and commercial products. Eur. Polym. J. 2015, 65, 252–267. [Google Scholar] [CrossRef] [Green Version]
  35. Kuoa, C.K.; Lia, W.-J.; Maucka, R.L.; Tuan, R.S. Cartilage tissue engineering: Its potential and uses. Curr. Opin. Rheumatol. 2006, 18, 64–73. [Google Scholar] [CrossRef] [PubMed]
  36. Zhao, X. Designing toughness and strength for soft materials. Proc. Natl. Acad. Sci. USA 2017, 114, 8138–8140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Lin, S.; Yuk, H.; Zhang, T.; Parada, G.A.; Koo, H.; Yu, C.; Zhao, X. Stretchable hydrogel electronics and devices. Adv. Mater. 2016, 28, 4497–4505. [Google Scholar] [CrossRef] [PubMed]
  38. Liu, X.; Tang, T.-C.; Tham, E.; Yuk, H.; Lin, S.; Lu, T.K.; Zhao, X. Stretchable living materials and devices with hydrogel—Elastomer hybrids hosting programmed cells. Proc. Natl. Acad. Sci. USA 2017, 114, 2200–2205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Osada, Y. Polymer gels as artificial soft tissue. Polym. Sci. Ser. C 2017, 59, 3–10. [Google Scholar] [CrossRef]
  40. Murakami, T.; Yarimitsu, S.; Nakashima, K.; Sakai, N.; Yamaguchi, T.; Sawae, Y.; Suzuki, A. Biphasic and boundary lubrication mechanisms in artificial hydrogel cartilage: A review. Proc. IMechE Part H J. Eng. Med. 2015, 229, 864–878. [Google Scholar] [CrossRef]
  41. Jiang, S.; Liu, S.; Feng, W. Pva hydrogel properties for biomedical application. J. Mech. Behav. Biomed. Mater. 2011, 4, 1228–1233. [Google Scholar] [CrossRef]
  42. Li, F.; Su, Y.; Wang, J.; Wu, G.; Wang, C. Influence of dynamic load on friction behavior of human articular cartilage, stainless steel and polyvinyl alcohol hydrogel as artificial cartilage. J. Mater. Sci. Mater. Med. 2010, 21, 147–154. [Google Scholar] [CrossRef]
  43. Li, F.; Zhang, G.; Wang, A.; Guo, F. The effects of surface mechanical deformation and bovine serum albumin on the tribological properties of polyvinyl alcohol hydrogel as an artificial cartilage. Adv. Mater. Sci. Eng. 2017. [Google Scholar] [CrossRef] [Green Version]
  44. Sardinha, V.; Lima, L.; Belangero, W.; Zavaglia, C.; Bavaresco, V.; Gomes, J. Tribological characterization of polyvinyl alcohol hydrogel as substitute of articular cartilage. Wear 2013, 301, 218–225. [Google Scholar] [CrossRef]
  45. McNary, M.; Athanasiou, K.; Reddi, A. Engineering lubrication in articular cartilage. Tissue Eng. Part B Rev. 2012, 18, 88–100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Seror, J.; Zhu, L.; Goldberg, R.; Day, A.J.; Klein, J. Supramolecular synergy in the boundary lubrication of synovial joints. Nat. Commun. 2015, 6, 6497. [Google Scholar] [CrossRef]
  47. Corvelli, M.; Che, B.; Saeui, C.; Singha, A.; Elisseeff, J. Biodynamic performance of hyaluronic acid versus synovial fluid of the knee in osteoarthritis. Methods 2015, 84, 90–98. [Google Scholar] [CrossRef] [Green Version]
  48. Singh, A.; Corvelli, M.; Unterman, A.S.; Wepasnick, K.A.; McDonnell, P.; Elisseeff, J.H. Enhanced lubrication on tissue and biomaterial surfaces through peptide-mediated binding of hyaluronic acid. Nat. Mater. 2014, 13, 988–995. [Google Scholar] [CrossRef] [PubMed]
  49. Swann, D.A.; Radin, E.L.; Nazimiec, M.; Weisser, P.A.; Curran, N.; Lewinnek, G. Role of hyaluronic acid in joint lubrication. Ann. Rheum. Dis. 1974, 33, 318–326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Li, F.; Wang, A.; Wang, C. Analysis of friction between articular cartilage and polyvinyl alcohol hydrogel artificial cartilage. Eng. Nano Eng. Approaches Med. Devices 2016, 27, 87–95. [Google Scholar] [CrossRef]
  51. Wan, W.; Bannerman, A.D.; Yang, L.; Mak, H. Poly (vinyl alcohol) cryogels for biomedical applications. Adv. Polym. Sci. 2014, 263, 283–321. [Google Scholar] [CrossRef]
  52. Hibbeler, R.C. Mecca de Materiales; Educaci, P., Ed.; Pearson: Mexico City, Mexico, 2011; ISBN 139786073205597. [Google Scholar]
  53. Gupta, S.; Sinha, S.; Sinha, A. Composition dependent mechanical response of transparent poly (vinyl alcohol) hydrogels. Colloids Surf. B Biointerfaces 2010, 78, 115–119. [Google Scholar] [CrossRef]
  54. Dodero, A.; Williams, R.; Gagliardi, S.; Vicini, S.; Alloisio, M.; Castellano, M. A micro-rheological and rheological study of biopolymers solutions: Hyaluronic acid. Carbohydr. Polym. 2019, 203, 249–355. [Google Scholar] [CrossRef]
  55. Jiang, H.; Zuo, Y.; Zhang, L.; Li, J.; Zhang, A.; Li, Y.; Yang, X. Property-based design: Optimization and characterizationof polyvinyl alcohol (pva) hydrogel and pva-matrix composite for artificial cornea. J. Mater. Sci. Mater. Med. 2014, 25, 941–952. [Google Scholar] [CrossRef]
  56. Huange, M.; Cai, D.; Liu, Y.; Sun, J.; Wang, J. Investigation of a-pva/s-pva hydrogels prepared by freezing-thawing method. Fibers Polym. 2012, 13, 955–962. [Google Scholar] [CrossRef]
  57. Marrella, A.; Lagazzo, A.; Dellacasa, E.; Pasquini, C.; Finocchio, E. 3d porous gelatin/pva hydrogel as meniscus substitute using alginate microparticles as porogens. Polymers 2018, 10, 380. [Google Scholar] [CrossRef] [Green Version]
  58. Bercea, M.; Morariu, S.; Teodorescu, M. Rheological investigation of poly (vinyl alcohol)/poly(n-vinyl pyrrolidone) mixture in aqueous solution solution and hydrogel state. J. Polym. Res. 2016, 23, 142–151. [Google Scholar] [CrossRef]
  59. Moore, D. Articular Cartilage. 2018. Available online: https://www.orthobullets.com/basic-science/9017/articular-cartilage (accessed on 1 February 2021).
  60. Wang, M.; Peng, Z. Wear in human knees. Biosurface Biotribol. 2015, 1, 98–112. [Google Scholar] [CrossRef] [Green Version]
  61. Li, W.; Wang, D.; Yang, W.; Song, Y. Compressive mechanical properties and microstructure of pva-ha hydrogels for cartilage repair. RSC Adv. 2016, 6, 20166–20172. [Google Scholar] [CrossRef]
  62. Mackley, M. Part IIb Chemical Engineering Tripos. Rheology and Processing. UK. 2011. Available online: http://www.malcolmmackley.com/wp-content/uploads/2012/11/Rheology-lectures-2011-section1_.pdf (accessed on 1 February 2021).
  63. Ambrosio, L.; Borzacchiello, A.; Netti, P.A.; Nicolais, L. Rheological study on hyaluronic acid and its derivative solutions. J. Macromol. Sci. Part A Pure Appl. Chem. 2007, 36, 991–1000. [Google Scholar] [CrossRef]
  64. Falcone, S.J.; Palmeri, D.M.; Berg, R.A. Rheological and cohesive properties of hyaluronic acid. J. Biomed. Mater. Res. Part A 2006, 76, 721–728. [Google Scholar] [CrossRef]
  65. Braithwaite, G.J.C.; Daley, M.J.; Toledo-Velasquez, D. Rheological and molecular weight comparisons of approved hyaluronic acid products—Preliminary standards for establishing class III medical device equivalence. J. Biomater. Sci. 2016, 27, 235–246. [Google Scholar] [CrossRef] [Green Version]
  66. Pan, Y.; Xiong, D.; Gao, F. Viscoelastic behavior of nano-hydroxyapatite reinforced poly (vinyl alcohol) gel biocomposite as an articular cartilage. J. Mater. Sci. Mater. Med. 2007, 19, 1963–1969. [Google Scholar] [CrossRef]
  67. Chan, M.; Neu, C.; Duraine, G.; Komvopoulos, K.; Reddi, A. Atomic force microscope investigation of the boundary-lubricant layer in articular cartilage. Osteoarthr. Cartil. 2010, 18, 956. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Kumar, P.; Oka, M.; Ikeuchi, K.; Shimizu, K.; Yamamuro, T.; Okumura, H.; Kotoura, Y. Low wear rate of uhmwpe against zirconia ceramic (y-psz) incomparison to alumina ceramic and sus 316l alloy. J. Biomed. Mater. Res. Part A 1991, 25, 813–828. [Google Scholar] [CrossRef] [PubMed]
  69. Caligaris, M.; Ateshian, G. Effect of sustained internal fluid pressurization under migrating contact area, and boundary lubrication by synovial fluid, on cartilage friction. Osteoarthr. Cartil. 2008, 16, 1220–1227. [Google Scholar] [CrossRef] [Green Version]
  70. Shi, Y.; Xiong, D.S.; Li, J.L. The friction properties of polyvinyl alcohol/graphene oxide hydrogels as cartilage replacement. Key Eng. Mater. 2017, 739, 152–156. [Google Scholar] [CrossRef]
  71. Park, R.K.; Nicoll, S.; Ateshian, G. Cartilage interstitial fluid load support in unconfined compression. J. Biomech. 2003, 36, 1785–1796. [Google Scholar] [CrossRef]
  72. Shi, Y.; Xiong, D. Microstructure and friction properties of pva/pvp hydrogels for articular cartilage repair as function of polymerization degree and polymer concentration. Wear 2013, 305, 280–285. [Google Scholar] [CrossRef]
  73. Pan, Y.; Xiong, D.; Chen, X. Friction characteristics of poly (vinyl alcohol) hydrogel as an articular cartilage biomaterial. Key Eng. Mater. 2007, 330, 1297–1300. [Google Scholar] [CrossRef]
  74. Kobayashi, M.; Hyu, H.S. Development and evaluation of polyvinyl alcohol-hydrogels as an artificial articular cartilage for orthopedic implants. Materials 2010, 3, 2753–2771. [Google Scholar] [CrossRef] [Green Version]
  75. Ricciardi, R.; Auriemma, F.; Gaillet, C.; de Rosa, C.; Laupretre, F. Investigation of the crystallinity of freeze/thaw poly (vinyl alcohol) hydrogels by different techniques. Macromolecules 2004, 37, 9510–9516. [Google Scholar] [CrossRef]
  76. Peppas, A.N.; Hassan, M.C. Structure and morphology of freeze/thawed pva hydrogels. Macromolecules 2000, 33, 2472–2479. [Google Scholar] [CrossRef]
Figure 1. Damage explanation in the joints: (A) Hip osteoarthritis (OA); (B) knee OA; (C) facet OA; (D) elbow OA.
Figure 1. Damage explanation in the joints: (A) Hip osteoarthritis (OA); (B) knee OA; (C) facet OA; (D) elbow OA.
Polymers 13 00746 g001
Figure 2. Characterization in a Discovery Hybrid Rheometer DHR-3. (a) Rheological study using cone-and-plate stainless steel configuration with an angle of 0.09969°, a diameter of 60 mm, and a truncation gap of 23 um. Samples used were 4 mg/mL hyaluronic acid (HA) solution. (b) Triborheological study using plate-and-plate stainless steel configuration with a diameter of 40 mm. Samples used were 10% polyvinyl alcohol (PVA) hydrogel with 5.0 ± 0.5 mm thickness and 45 mm diameter, hydrated with 4 mg/mL HA solution.
Figure 2. Characterization in a Discovery Hybrid Rheometer DHR-3. (a) Rheological study using cone-and-plate stainless steel configuration with an angle of 0.09969°, a diameter of 60 mm, and a truncation gap of 23 um. Samples used were 4 mg/mL hyaluronic acid (HA) solution. (b) Triborheological study using plate-and-plate stainless steel configuration with a diameter of 40 mm. Samples used were 10% polyvinyl alcohol (PVA) hydrogel with 5.0 ± 0.5 mm thickness and 45 mm diameter, hydrated with 4 mg/mL HA solution.
Polymers 13 00746 g002
Figure 3. SEM of PVA hydrogels. (a) Surface distribution, (b) porosity (High vacuum, 50–1000×, secondary electrons).
Figure 3. SEM of PVA hydrogels. (a) Surface distribution, (b) porosity (High vacuum, 50–1000×, secondary electrons).
Polymers 13 00746 g003
Figure 4. Compression tests (Up to 4500 N, compression rate 1 mm/min). (a) Linear regression of all tests; (b) exponential behavior of compression modulus in test 1; (c) test 5 deformation over time.
Figure 4. Compression tests (Up to 4500 N, compression rate 1 mm/min). (a) Linear regression of all tests; (b) exponential behavior of compression modulus in test 1; (c) test 5 deformation over time.
Polymers 13 00746 g004
Figure 5. Rheological test. (a) Viscosity vs. shear rate; (b) shear stress vs. shear rate; (c) modules and complex viscosity vs. frequency (humidity, 37 °C, 0.01–100 Hz, 0.001–100 s1).
Figure 5. Rheological test. (a) Viscosity vs. shear rate; (b) shear stress vs. shear rate; (c) modules and complex viscosity vs. frequency (humidity, 37 °C, 0.01–100 Hz, 0.001–100 s1).
Polymers 13 00746 g005
Figure 6. Coefficient of friction vs. angular velocity. (a) Dry; (b) HA lubricated; (c) coefficient of friction vs. angular velocity: Lubrication comparative; (d) coefficient of friction vs/pressure (0.1–100 rad/s, 1–20 N).
Figure 6. Coefficient of friction vs. angular velocity. (a) Dry; (b) HA lubricated; (c) coefficient of friction vs. angular velocity: Lubrication comparative; (d) coefficient of friction vs/pressure (0.1–100 rad/s, 1–20 N).
Polymers 13 00746 g006
Figure 7. (a) Lubrication mechanism in PVA hydrogel under applied force in solid–solid contact. (b) Lubrication mechanisms applied to a PVA hydrogel/HA lubricant system.
Figure 7. (a) Lubrication mechanism in PVA hydrogel under applied force in solid–solid contact. (b) Lubrication mechanisms applied to a PVA hydrogel/HA lubricant system.
Polymers 13 00746 g007
Figure 8. Differential scanning calorimetry (DSC) thermal diagram of PVA samples at different weights.
Figure 8. Differential scanning calorimetry (DSC) thermal diagram of PVA samples at different weights.
Polymers 13 00746 g008
Table 1. Water content measurements of dry and hydrated samples at different weights.
Table 1. Water content measurements of dry and hydrated samples at different weights.
TestWh (g)Ws (g)EWCM (%)
11.480.2384.46
21.720.2684.88
31.730.2684.97
41.910.3183.77
52.450.4282.86
62.540.4482.68
72.700.4981.85
83.110.5981.03
93.320.5084.94
103.380.6481.07
113.450.6780.58
123.620.5684.53
133.720.6383.06
144.320.7981.50
154.730.8282.66
Average82.99
Standard Deviation1.50
Table 2. DSC results for melting temperature and crystallinity. Measurements were made on a Q10 TA Instruments. Samples were heated at a rate of 5 °C/min from 25–250 °C in a nitrogen atmosphere.
Table 2. DSC results for melting temperature and crystallinity. Measurements were made on a Q10 TA Instruments. Samples were heated at a rate of 5 °C/min from 25–250 °C in a nitrogen atmosphere.
Sample Weight (mg)Tm (°C)ΔHm (J/g)Xc (%)
2.1231.8313.9210.04
3.4232.799.787.06
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Duque-Ossa, L.C.; Ruiz-Pulido, G.; Medina, D.I. Triborheological Study under Physiological Conditions of PVA Hydrogel/HA Lubricant as Synthetic System for Soft Tissue Replacement. Polymers 2021, 13, 746. https://doi.org/10.3390/polym13050746

AMA Style

Duque-Ossa LC, Ruiz-Pulido G, Medina DI. Triborheological Study under Physiological Conditions of PVA Hydrogel/HA Lubricant as Synthetic System for Soft Tissue Replacement. Polymers. 2021; 13(5):746. https://doi.org/10.3390/polym13050746

Chicago/Turabian Style

Duque-Ossa, Laura C., Gustavo Ruiz-Pulido, and Dora I. Medina. 2021. "Triborheological Study under Physiological Conditions of PVA Hydrogel/HA Lubricant as Synthetic System for Soft Tissue Replacement" Polymers 13, no. 5: 746. https://doi.org/10.3390/polym13050746

APA Style

Duque-Ossa, L. C., Ruiz-Pulido, G., & Medina, D. I. (2021). Triborheological Study under Physiological Conditions of PVA Hydrogel/HA Lubricant as Synthetic System for Soft Tissue Replacement. Polymers, 13(5), 746. https://doi.org/10.3390/polym13050746

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop