Next Article in Journal
Additive Manufactured Parts Produced Using Selective Laser Sintering Technology: Comparison between Porosity of Pure and Blended Polymers
Next Article in Special Issue
Recent Approaches to the Plasticization of Poly(lactic Acid) (PLA) (A Review)
Previous Article in Journal
Developing Mixed Matrix Membranes with Good CO2 Separation Performance Based on PEG-Modified UiO-66 MOF and 6FDA-Durene Polyimide
Previous Article in Special Issue
Green, Eco-Friendly, Highly Biocompatible and Bioactive Nanocomposite-Based Biopolymers Loaded with ZnO@Fe3O4 Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magnetic Resonance-Based Analytical Tools to Study Polyvinylpyrrolidone–Hydroxyapatite Composites

1
Institute of Physics, Kazan Federal University, Kremlyovskaya St.18, 420008 Kazan, Russia
2
A.A. Baikov Institute of Metallurgy and Material Science, Russian Academy of Sciences, Leninsky Avenue 49, 119334 Moscow, Russia
3
Faculty of Material Science and Engineering, National University of Science and Technology Politehnica Bucharest, 313 Splaiul Independentei Street, District 6, 060042 Bucharest, Romania
4
Academy of Romanian Scientists, 54 Splaiul Independentei Street, District 5, 050094 Bucharest, Romania
5
Istituto di Struttura della Materia, Consiglio Nazionale delle Ricerche (ISM-CNR), Via del Fosso del Cavaliere, 100, 00133 Rome, Italy
6
Department of Analytical, Physical and Colloid Chemistry, I.M. Sechenov First Moscow State Medical University, Trubetskaya Str., Build. 8/2, 119048 Moscow, Russia
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(22), 4445; https://doi.org/10.3390/polym15224445
Submission received: 5 October 2023 / Revised: 13 November 2023 / Accepted: 13 November 2023 / Published: 17 November 2023
(This article belongs to the Special Issue Advances in Biocompatible and Biodegradable Polymers III)

Abstract

:
The synthesis of biocompatible and bioresorbable composite materials, such as a “polymer matrix-mineral constituent,” stimulating the natural growth of living tissues and the restoration of damaged parts of the body, is one of the challenging problems in regenerative medicine and materials science. Composite films of bioresorbable polymer of polyvinylpyrrolidone (PVP) and hydroxyapatite (HA) were obtained. HA was synthesized in situ in the polymer solution. We applied electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) approaches to study the composite films’ properties. The application of EPR in two frequency ranges allowed us to derive spectroscopic parameters of the nitrogen-based light and radiation-induced paramagnetic centers in HA, PVP and PVP-HA with high accuracy. It was shown that PVP did not significantly affect the EPR spectral and relaxation parameters of the radiation-induced paramagnetic centers in HA, while light-induced centers were detected only in PVP. Magic angle spinning (MAS) 1H NMR showed the presence of two signals at 4.7 ppm and −2.15 ppm, attributed to “free” water and hydroxyl groups, while the single line was attributed to 31P. NMR relaxation measurements for 1H and 31P showed that the relaxation decays were multicomponent processes that can be described by three components of the transverse relaxation times. The obtained results demonstrated that the applied magnetic resonance methods can be used for the quality control of PVP-HA composites and, potentially, for the development of analytical tools to follow the processes of sample treatment, resorption, and degradation.

1. Introduction

Calcium phosphates (CaPs) are the main inorganic component of bone and dental tissues [1,2,3,4,5]. Materials based on CaPs—powders, ceramics, cements, coatings, and composites—are widely used in various applications, especially in medicine for the replacement and restoration of damaged bone tissues. The most widespread representative of CaPs is hydroxyapatite (HA, Ca10(PO4)6(OH)2, Figure 1), which is present in bone tissues in the form of nanocrystals. HA is known to be biocompatible, non-toxic, and the most important mineral that constitutes bones and teeth [6].
It is known that the synthesis of calcium phosphates in the presence of polymers (collagen, gelatin, starch, chitosan, etc.) leads to the formation of CaP nanocrystals with controlled size, morphology, and improved mechanical properties [7,8,9,10]. Currently, 3D-printing technologies are rapidly developing in medicine, in particular for printing bone implants. CaP frames printed using a 3D printer for biomedical applications have good mechanical strength [11]. In the standard 3D technology for creating such objects, composites made of a biodegradable polymer and a calcium-containing filler are applied for layer-by-layer deposition of the material [12]. Polyvinylpyrrolidone (PVP, Figure 2) is widely employed as a synthetic polymer for various biomedical purposes in material engineering due to its diverse properties, including solubility in water and in a broad range of liquid media. PVP also has both hydrophilic and hydrophobic functional groups, due to which it interacts with various solvents (easily soluble in cold water and also soluble in many organic solvents, including alcohols; some chlorinated compounds such as chloroform, methylene chloride and ethylene dichloride; nitroparaffins; and amines [13,14,15,16]). Thanks to its biocompatibility, absence of toxicity and high capacity to form interpolymer complexes, PVP is widely used for designing materials for different applications, such as biomaterials for medical uses.
Bone cements are often used to fix artificial prostheses to the human skeletal system. These materials provide immediate fixation of the implant and ensure a better distribution of body load between the prosthesis and bone. Morejón et al. [17] studied polymer particles with particle size distribution and molecular weight, allowing them to develop bone cement according to the international standards. Needle-shaped HA crystals were obtained in an aqueous solution of PVP via precipitation from solutions of calcium nitrate and orthophosphoric acid at 60 °C. PVP was used to regulate the nucleation and growth of HA crystals. Based on the results obtained from X-ray diffraction and high-resolution transmission electron microscopy measurements, it was concluded that the addition of PVP promotes the growth of HA nanocrystals [18,19]. It was shown that the number of HA crystals formed on the polymer surface depends on the number of negative charges on this surface and the pH of the reaction medium. It was established that to obtain a controlled morphology, the temperature and duration of the reaction are of primary importance [20].
Polymer—CaPs nanocrystal composites can be obtained in situ, without mixing the components, similar to the process happening in the body during bone remodeling. In ref. [21], such a composite was obtained via biomimetic deposition of HA on a PVP matrix. A colorimetric test to assess the metabolic activity of cells (MTT-test) showed that the composites were biocompatible and could be used to fill bone defects. Guesmi Y. et al. [22] investigated the material obtained by grafting PVP onto the surface of HA microcrystals, suitable for bone tissue engineering. Based on the results of the bioactivity study of the composites, it was concluded that the PVP seeded on their surface cells had a higher viability compared to those cultivated on the HA crystals.
Obviously, the properties of polymer-CaP composites depend on the chemical interaction of the components. It is supposed that the inorganic particles interact with the organic component by establishing the inter- and intramolecular hydrogen bonds and via the ion–dipole forces, formed, for example, between the calcium ions of CaPs and the functional groups of polymers [23,24,25,26]. However, the interaction of PVP with calcium phosphate microcrystals has not yet been sufficiently studied.
A lot of work has been conducted on the development and application of various characterization methods that can provide details about calcium phosphate-based biomaterials at different scales, such as X-ray diffraction, infrared and electron microscopy, rheological analysis, water content determination, etc. [8]. In the present work, PVP-HA composites were investigated using electron paramagnetic resonance (EPR) and nuclear magnetic resonance (NMR) techniques. Both of these analytical methods are non-invasive, widely applied in material science, polymer research and for biomedical purposes, but very rarely exploited for polymer-CaPs studies.
Pure HA and PVP are EPR silent. Therefore, due to its inherent high sensitivity, EPR can be used to check for the presence of various paramagnetic impurities, like metals and metal complexes, or defects in the polymer chains both in the starting materials for synthesis and in the resulting products [27,28,29,30,31,32,33,34]. Stable and unstable defects created by light or ionizing irradiation (X-ray irradiation, for example) can be exploited not only for detecting the presence and concentration of impurities but also as paramagnetic probes to study the local environment, material structure and changes in the local lattice [27]. Several stable radiation-induced anions were detected and investigated using EPR in synthetic and biogenic CaPs, like carbonates, nitrogen and oxygen-containing species [30,32,33,34]. It was shown, for example, that the spectra of nitrogen-containing stable radicals depend on the type of CaP and, therefore, EPR may be used to follow the processes of calcium phosphates growth [28,29]. In HA, the mentioned anions can substitute both hydroxyl (OH) and phosphate (PO4) groups (see Figure 1), as well as occupy interstitial sites. The localization of anions depends on the nature of biogenic materials, the method used for the synthesis of artificial compounds, co-doping effects, etc. Advanced EPR techniques, together with theoretical calculations, often allow one to define the location of these anions in the HA structure, and to propose the most probable charge compensation schemes, in the case of non-isovalent substitutions [27,29].
EPR techniques are used to investigate intentional or non-intentional HA cation doping with Mn2+, Cu2+, etc. Two distinct sites, Ca(1) and Ca(2), are distinguished in the HA lattice (Figure 1), and EPR data allow one to confirm or refute assumptions about the introduction of the dopants into the HA structure, preferential substitution or re-distribution of dopants between the calcium sites [27]. Even in the case of non-paramagnetic impurities, such as magnesium and aluminum, for example, it was shown that the EPR of anionic impurities is sensitive to the lattice distortion caused by these cations [29].
Due to the chemical complexity of polymer-CaPs composites, hyperfine EPR and NMR can be engaged, taking into account that different nuclei, such as 31P (with the nuclear spin of I = 1/2), 1H (I = 1/2), 13C (I = 1/2), 14N (I = 1) and even 43Ca, can act as nuclear spectroscopic probes. Since HA is practically insoluble in liquid media, magic angle spinning (MAS) NMR techniques are used to investigate materials containing hydroxyapatite [35,36,37,38,39,40,41,42]. Thanks to the advancements in instrumentation, pulse sequence development and data interpretation (through combined experimental–computational approaches) made by the NMR community since the beginning of 2000s, materials scientists have been able to extend the limits of NMR research in order to obtain a deeper knowledge of CaPs biomaterials [36]. NMR studies of native bone and dental tissues have allowed us to measure the size of the amorphous calcium phosphate layer at the surface of apatite bone crystallites, to detect minor organic species (like nucleic acids) and to measure the distance and binding geometry of small molecular ions at the surface of mineral particles, like citrate. Furthermore, significant information about inorganic and organic components, as well as the exclusion of previously suggested models regarding the bone structure and its formation, has been obtained by comparing the NMR signatures of natural materials with those of synthetically derived models, including in vitro-cultured tissues [36]. The achievements in this field give hope that in the near future, novel approaches can be implemented into clinical magnetic resonance imaging (MRI) for the evaluation of skeleton quality, 3D and 4D synthetic tissue constructs and the diagnosis of bone diseases in order to provide additional biomarkers for the assessment of bone microarchitecture, etc. [37].
The effect of PVP on the morphology and size of nanocomposites was studied using NMR in refs. [17,40]. In several papers from our group, we demonstrated the feasibility of combined EPR/NMR studies to investigate PVP-HA-sodium alginate composites, including the influence of divalent metals on the cross-linking of the polymer and HA [35,41]. However, to the best of our knowledge, no comprehensive analysis of synthesized in situ PVP-HA composites using the EPR and NMR techniques is present in literature.

2. Materials and Methods

2.1. Synthesis of Composite Materials Based on Polyvinylpyrrolidone with Hydroxyapatite

Solutions of Ca(NO3)2 × 4H2O (chemical grade, PanReacAppliChem, Barcelona, Spain), (NH4)2HPO4 (chemical grade, Chimmed, Moscow, Russia) and (PVP) Mw = 12 kDa (Boai NKY Pharmaceuticals Ltd., Tianjin, China) were used for the in situ synthesis of HA powders in the PVP solution. Precipitation of calcium phosphates was performed at room temperature (20–25 °C) at pH of ~11.5, in accordance with Equation (1) [41]:
10Ca2+ + 6HPO42− + 8NH3 H2O → Ca10(PO4)6(OH)2 + 8NH4+ + 6H2O
The in situ synthesis of PVP-HA composites was carried out according to the procedure described in ref. [40]. Briefly, an aqueous PVP solution was prepared by dissolving 14.35 g of PVP in 200 mL of distilled water via stirring, in order to obtain 5.52 wt% concentration. After the formation of a homogeneous water–polymer mixture, 20 mL of 0.1 mol/L calcium nitrate solution was added. To regulate the reaction acidity, 20 mL of 25% NH4OH aqueous solution was added as well. Then, 20 mL of 0.06 mol/L diammonium phosphate solution was added drop by drop while stirring at 500 rpm. To remove the by-products of the HA formation (reaction (1)) (i.e., nitrate ions) from the composite, it was washed using dialysis.
The product of reaction (1) was examined with X-ray diffraction (XRD) and Fourier transform infrared spectroscopy (FTIR). The XRD analysis was carried out by using the DRONE 3M diffractometer (S.-Petersburg, Russia) with Cu Kα radiation (λ = 0.154 nm). The FTIR spectra were obtained with the NicoletAvatar-330 (ThermoScientific, Waltham, MA, USA) spectrometer in the range of 400–4000 cm−1 on the samples mixed with potassium bromide. The annealing of samples before XRD and FTIR was performed for 1 h at 400 and 900 °C. The details of the samples characterization were reported in ref. [41].
From XRD, it was found that the resulting composites are X-ray amorphous. As a result of heat treatment at 400 °C, an apatite structure was formed, and after heat treatment at 900 °C, only crystallized HA without admixtures of other phases was detected.
FTIR measurements showed that in the sample without heat treatment, despite a partial overlap of PVP and phosphate group bands, it is still possible to distinguish the bands related to PVP. After the comparison of FTIR spectra with the data reported in ref. [41], it was found that the most intense bands at 1660 cm−1 correspond to the valence vibrations of the carbonyl group (C=O). In the spectra related to HA, the oscillations for the PO43− group were pronounced (triplets of valence oscillations at 1090, 1053 and 965 cm−1 and at 632, 572 and 472 cm−1). The intense oscillation band of the OH-group at 3570 cm−1, which is present even in the sample without heat treatment, may indicate the beginning of the formation of the HA structure. With an increase in the annealing temperature, the intensity of this band significantly increased, allowing it to be associated with the crystalline hydroxyapatite formation. The intensity and resolution of the bands corresponding to the PO43− group also increased. It should be noted that the PVP bands practically disappear after heat treatment at 400 °C, and completely disappear after annealing at 900 °C.
Scanning electron microscopy (SEM) measurements were performed using a Tescan VEGA3 (Kohoutovice, Czech Republic). Samples were coated with a thin layer of gold for the SEM examinations. The SEM images were acquired using secondary electron (SE) and backscattered electron (BSE) imaging modalities. The results of SEM investigations are presented in Section 3.1.

2.2. EPR

The EPR measurements were performed in two microwave ranges, X-band (with the microwave frequency νMW = 9 GHz) and W-band (νMW = 94 GHz), by exploiting the Bruker Elexsys E580/E680 spectrometer (Karlsruhe, Germany) in conventional (cw) and pulsed modes. Pulse methods are able to provide much more information compared to that obtained using conventional EPR [43]. In the pulsed mode, the Hahn sequence was applied to detect electron spin echo (ESE): π/2–τ–π–τ–; here, the duration of π/2 was equal to 64 ns and τ = 250 ns. Registration of the EPR spectra was performed by detecting the ESE integral intensity depending on the magnetic field value B0. The choice of the W-band was justified by the need to achieve a higher spectroscopic resolution (that allows identifying distinct EPR signals with close g-factors) and high sensitivity (to register weak EPR signals). The spectra were recorded both at room temperature (T = 297 K) and at lower temperature (T = 200 K). Stable paramagnetic centers were formed under X-ray irradiation of the synthesized powders using the URS-55 source (G = 55 kV, I = 16 mA, W-anticathode) at room temperature for 30 min with an estimated dose of 5 kGy and by illuminating the samples in the EPR cavity via laser sources with different wavelengths (λ) from 266 nm up to 355 nm.
The most abundant metal impurities in CaP detected with the EPR method are related to manganese, chromium, iron and copper. Low-temperature measurements and the use of high-frequency methods allowed us to increase the sensitivity of EPR. Details of the EPR measurements, including electronic relaxation times for HA, are given in ref. [43]. Simulations of the EPR spectra were performed by means of the EasySpin package [44].

2.3. Nuclear Magnetic Resonance

For all the samples, the MAS-NMR spectra, the times of spin–lattice and spin–spin relaxation on phosphorus (31P) nuclei and hydrogen (1H) protons were measured. For measurements, an AVANCE400WB NMR spectrometer (Karlsruhe, Germany) with a MAS 4BL CP BB DVT sensor was employed. The resonant frequency on protons was about 400.27 MHz, and on 31P nuclei, it was 162.034 MHz. The duration of a π/2 pulse on protons was 2.5 µs (output power of 94 W), and on 31P nuclei, it was 3.4 µs at a power of 50 W. The standard pulse programs, such as onepulseq, hpdec, cp, inversion-recovery and cpmg, were used to obtain NMR spectra and relaxation on 1H, 31P and 13C. The spectrum width for all measurements was 100 kHz. Powder samples were densely packed in a 4 mm zirconium oxide rotor and spun up to a rotation frequency of 7 kHz. The measurements were carried out at 295 K. The chemical phenomena related to the NMR signals were calibrated by measuring the signal of water (δ = 4.67 ppm) for protons and 85% phosphoric acid for phosphorus (δ = 0 ppm).

3. Results

3.1. SEM

As it follows from the SEM experiments (Figure 3), the microstructure of the PVP-HA composite is represented by rod-shaped particles about 30 µm long, coated with polymer. It resembles the results reported in ref. [18]. The polymer facilitates the formation of the needle-like HA aggregates, whereas, in the absence of PVP, it is difficult to obtain a rod-like structure. This can be attributed to the Ostwald ripening process occurring in the solution and leading to a typical oriented attachment process, especially in the presence of PVP.

3.2. EPR

For HA, PVP and PVP-HA, the EPR signal was not observed due to the absence of paramagnetic centers in the structure of the materials. It additionally proved the purity of the initial materials and the reaction product. After the X-ray irradiation of HA powders in the X-band, the EPR spectrum containing three lines was registered (Figure 4). The spectroscopic parameters of the EPR spectrum allowed us to refer the detected signal to the well-studied NO32 stable species in HA [45,46,47] (see Table 1), due to the hyperfine interaction (A) with the 14N nuclei (I = 1). The PVP samples were also studied after X-ray irradiation (Figure 5, red lines) and photo-induced EPR at low temperatures (Figure 5), since it is known that exposure to the visible light does not lead to the formation of stable radicals in HA [43], while it does in PVP polymer chains [48]. In the W-band, due to the higher spectral resolution, a three-line pattern of PVP was observed (Figure 6). It accounted for the localization of an unpaired electron on the PVP nitrogen atom (Figure 2) with a constant A|| = 106 ± 10 MHz. In the X-band, the three-line pattern for PVP is not resolved.
As can be seen from Figure 5 (blue lines), the EPR spectrum of the PVP-HA contains only the NO32 radicals from HA, and the component from PVP was not observed. It can be supposed that in the PVP-HA mixture, the radiation-induced centers in PVP have a competitive electron trap channel, which is possible only in the presence of a chemical bond between the components. In ref. [49], we also showed that for the PVP-HA samples obtained ex situ, the values of the A-components and their distribution for the electron-14N interaction grew with the increase in the amount of PVP. A change in the crystal structure of HA in the reaction with PVP is unlikely. Therefore, we can assume that during the formation of chemical bonds, the PVP molecules create a positively charged layer around the HA particles, which increases the electron density in the near-surface layer of the HA particle.
Electronic spin–lattice relaxation measurements did not reveal an influence of PVP on the relaxation parameters of the radiation-induced NO32− radicals from HA, as is demonstrated in Table 2 for the X-band experiments. This correlates with the tiny influence of PVP on the EPR spectra of HA, as described above.

3.3. NMR

The 1H NMR spectra for PVP and PVP-HA are shown in Figure 7. The 1H MAS-NMR spectrum of the PVP consists of one unresolved NMR signal at 2.4 ppm, associated with the CH2 groups of the pyridine ring (Figure 2) and the polymer chain. The CH groups of the polymer chain contribute to the signal at about 5 ppm, while the CH3 end groups give a small contribution to the region of 0.5–1.5 ppm. On the 1H NMR spectrum in Figure 7, the contributions of these chemical groups are seen as small shoulders to the right and left of the main NMR signal.
The 1H NMR spectrum of a PVP-HA sample consists of two 1H NMR signals at 4.7 ppm and −2.15 ppm [50,51]. The integrated intensity of the NMR signal at 4.7 ppm is significantly higher compared to the intensity of the NMR signal at −2.15 ppm. The 1H NMR signal at 4.7 ppm for the PVP-HA is associated with «free» water adsorbed on the surface of HA, while the 1H NMR signal at −2.15 ppm is associated with hydroxyls present in the HA structure. The shape of the NMR signal at 4.7 ppm is asymmetric due to the superposition of the proton signals of water and polymer.
The 13C{1H} MAS-NMR spectra of the PVP on 13C during proton decoupling, obtained using the pulse program cp (red) and hpdec (blue), respectively, are shown in Figure 8. The MAS-NMR 13C{1H} spectrum of PVP-HA has the same form as in Figure 7, only the signal-to-noise ratio is several times lower. The NMR signals of the aliphatic region for cross-polarization have lower resolution than those present the 13C NMR spectrum with proton suppression.
The data obtained for the spin–lattice relaxation rates on protons for the two samples (PVP and PVP-HA) are shown in Table 3. The table lists the chemical shifts of the signals for which the spin–lattice relaxation rates were measured.
The relaxation decays obtained via the Carr–Purcell–Meiboom–Gill (cpmg) sequence for the NMR signals of the chemical shift regions, according to Table 3, are presented in Figure 9. As can be observed, the proton relaxation decays for the NMR signals from the PVP and PVP-HA differ from each other and can be described by exponent approximation. The decay curves of the transverse magnetization of protons in the PVP and PVP-HA samples were approximated by the following formula:
A ( t ) = A ( 0 ) i = 1 n p i exp ( t T 2 ( i ) ) ,
where A(0) is the initial amplitude of the NMR signal and T2(i) is the spin–spin relaxation times of components with populations pi. The values of the spin–spin relaxation times and their populations calculated from the relaxation decays are given in Table 4.
It is known that multicomponent relaxation can manifest itself as a result of different types of molecular motion or different surroundings of the nuclei. We assume that the presence of three transverse relaxation components for protons is due to the different positions of the OH groups in the PVP-HA composite. So, hydroxyls, included in the structure of HA (δ = −2.15 ppm), should be characterized by the longest time of transverse (spin–spin) magnetization. In the environment of the HA crystal defects, hydroxyls should be described by high rates of nuclear spin–spin relaxation. The same approach applies to water localized on the HA surface.
The 31P NMR spectrum of the PVP-HA sample is shown in Figure 10. It is a single line with a width of about 290 Hz with an isotropic chemical shift of 3 ppm. The shape of the spectrum indicates that all the phosphorus nuclei in HA are in the same chemical environment. It should be noted that the stationary 31P NMR spectrum of the PVP-HA sample (without rotation) has the same chemical shift and a linewidth of 3000 Hz.
The spin–lattice relaxation time for 31P was measured as T1 = 238 ± 21 s. The relaxation decay of the 31P NMR signal for the PVP-HA and the fitting curve based on the relaxation decay decomposition parameters, according to Equation (2), are shown in Figure 11.
Figure 11. Relaxation decay of the 31P NMR signal for the PVP-HA. White/red circles are the experimental data. Black line is the approximation according to Equation (2) with the parameters listed in Table 5.
Figure 11. Relaxation decay of the 31P NMR signal for the PVP-HA. White/red circles are the experimental data. Black line is the approximation according to Equation (2) with the parameters listed in Table 5.
Polymers 15 04445 g011
The values of spin–spin relaxation times and their populations obtained from the relaxation decay for the line in the 31P NMR spectrum of the PVP-HA are listed in Table 5. The presence of three transverse relaxation components both for phosphorus atoms and protons is probably due to the presence of crystal lattice defects in the HA.
Table 5. The values of 31P spin–spin (transverse) relaxation times of components with populations pi, according to Equation (2).
Table 5. The values of 31P spin–spin (transverse) relaxation times of components with populations pi, according to Equation (2).
Chemical Shifts of NMR Signals (ppm)p1T2(1),
ms
p2T2(2),
ms
p3T2(3),
ms
PVP-HA
30.24455.5 ± 0.50.1894.2 ± 0.30.5670.85 ± 0.05

4. Conclusions

  • Using two microwave frequencies, the EPR parameters for the light- and radiation-induced paramagnetic centers in the PVP, HA and PVP-HA (such as components of g-factors and hyperfine constants A between electrons and 14N nuclei) were defined with high accuracy (Table 1). This shows the possibility of using EPR for the qualitative (synthesis, presence of impurities) and quantitative (concentration of impurities and defects) control of the initial materials (PVP, HA) and the final composite PVP-HA product.
  • In the PVP-HA composite, the EPR spectra and electronic relaxation times of the radiation-induced paramagnetic centers were very close to those in HA. In the PVP-HA, the distribution of A components for the electron–14N interaction for NO32− radicals was larger than for the pure HA (Table 1). The PVP-HA composites did not contain light-induced radicals characteristic of PVP. This can be ascribed to the re-distribution of the electrical charges between PVP and HA. These results can be used for the control of the quality and success of the in situ synthesis of PVP-HA composites using the EPR techniques.
  • In the 1H MAS NMR spectrum of the PVP-HA, the presence of two signals at 4.7 ppm and −2.15 ppm were attributed to “free” water and hydroxyl groups, and a single line attributed to 31P was registered.
  • The NMR relaxation measurements for 1H and 31P showed that the relaxation decays are multicomponent processes that can be described by three components of the transverse relaxation times. Multicomponent relaxation decay can be ascribed to the presence of defects in the HA lattice. The obtained data can serve as a basis for future NMR applications in clinical MRI for the evaluation of skeleton quality.
  • The obtained results demonstrate that the applied magnetic resonance techniques can be used for the quality control of synthesis products and, potentially, to follow the processes of the samples’ treatment, resorption and degradation [36,37].

Author Contributions

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

Funding

This work was funded by the subsidy allocated to Kazan Federal University for the state assignment in the sphere of scientific activities number FZSM-2023-0016. The synthesis of composites was carried out with the support of the Russian science foundation, grant No. 22-23-00278.

Data Availability Statement

The experimental data on the results reported in this manuscript are available upon an official request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vallet-Regi, M. Ceramics for Medical Applications. J. Chem. Soc. Dalton Trans. 2001, 97–108. [Google Scholar] [CrossRef]
  2. Afshar, A.; Ghorbani, M.; Ehsani, N.; Saeri, M.R.; Sorrell, C.C. Some important factors in the wet precipitation process of hydroxyapatite. Mater. Des. 2003, 24, 197–202. [Google Scholar] [CrossRef]
  3. Riman, R.E.; Suchaned, W.L.; Byrappa, K.; Chen, C.W.; Shuk, P. Solution synthesis of hydroxyapatite designer particulates. Solid State Ion. 2002, 151, 393–402. [Google Scholar] [CrossRef]
  4. Chen, J.D.; Wang, Y.J.; Wei, K.; Zhang, S.H.; Shi, X.T. Self-organization of hydroxyapatite nanorods through oriented attachment. Biomaterials 2007, 28, 2275–2280. [Google Scholar] [CrossRef]
  5. Mobasherpour, I.; Heshajin, M.S.; Kazemzadeh, A.; Zakeri, M. Synthesis of nanocrystalline hydroxyapatite by using precipitation method. J. Alloys Compd. 2007, 430, 330–333. [Google Scholar] [CrossRef]
  6. Barinov, S.M. Calcium phosphate-based ceramic and composite materials for medicine. Russ. Chem. Rev. 2010, 79, 13. [Google Scholar] [CrossRef]
  7. Dziadek, M.; Stodolak-Zych, E.; Cholewa-Kowalska, K. Biodegradable ceramic-polymer composites for biomedical applications. A review. Mater. Sci. Eng. C 2017, 71, 1175–1191. [Google Scholar] [CrossRef]
  8. Vranceanu, M.D.; Saban, R.; Antoniac, I.; Albu, M.; Miculescu, F. Development and characterization of novel porous collagen based biocomposite for bone tissue regeneration. UPB Sci. Bull. Ser. B 2012, 74, 145–156. [Google Scholar]
  9. Armentano, I.; Dottori, M.; Fortunati, E.; Mattioli, S.; Kenny, J.M. Biodegradable polymer matrix nanocomposites for tissue engineering: A review. Polym. Degrad. Stab. 2010, 95, 2126–2146. [Google Scholar] [CrossRef]
  10. Nair, L.; Laurencin, C. Biodegradable polymers as biomaterials. Prog. Polym. Sci. 2007, 32, 762–798. [Google Scholar] [CrossRef]
  11. Bandyopadhyay, A.; Mitra, I.; Bose, S. 3D printing for bone regeneration. Curr. Osteoporos. Rep. 2020, 18, 505–514. [Google Scholar] [CrossRef] [PubMed]
  12. Bergmann, C.; Lindner, M.; Zhang, W.; Koczur, K.; Kirsten, A.; Telle, R.; Fischer, H. 3D printing of bone substitute implants using calcium phosphate and bioactive glasses. J. Eur. Ceram. Soc. 2010, 30, 2563–2567. [Google Scholar] [CrossRef]
  13. Teodorescu, M.; Bercea, M. Poly (vinylpyrrolidone)—A versatile polymer for biomedical and beyond medical applications. Polym. Plast. Technol. Eng. 2015, 54, 923–943. [Google Scholar] [CrossRef]
  14. Haaf, F.; Sanner, A.; Straub, F. Polymers of N-vinylpyrrolidone: Synthesis, characterization and uses. Polym. J. 1985, 17, 143–152. [Google Scholar] [CrossRef]
  15. Anderson, C.C.; Rodriguez, F.; Thurston, D.A. Crosslinking aqueous poly(vinyl pyrrolidone) solutions by persulfate. J. Appl. Polym. Sci. 1979, 23, 2453–2462. [Google Scholar] [CrossRef]
  16. Raimi-Abraham, B.T.; Mahalingam, S.; Edirisinghe, M.; Craig, D.Q.M. Generation of poly (N-vinylpyrrolidone) nanofibres using pressurized gyration. Mater. Sci. Eng. C 2014, 39, 168–176. [Google Scholar] [CrossRef] [PubMed]
  17. Morejón, L.; Mendizábal, E.; Delgado, J.A.; Davidenko, N.; López-Dellamary, F.; Manríquez, R.; Ginebra, M.P.; Gil, F.J.; Planell, J.A. Synthesis and characterization of poly (methyl methacrylate-styrene) copolymeric beads for bone cements. Lat. Am. Appl. Res. 2005, 35, 175–182. [Google Scholar]
  18. Zhang, Y.; Lu, J. A Mild and Efficient Biomimetic Synthesis of Rodlike Hydroxyapatite Particles with a High Aspect Ratio Using Polyvinylpyrrolidone as Capping Agent. Cryst. Growth Des. 2008, 8, 2101–2107. [Google Scholar] [CrossRef]
  19. Jemli, Y.E.; Abdelouahdi, K.; Minh, D.P.; Barakat, A.; Solhy, A. Synthesis and Characterization of Hydroxyapatite and Hydroxyapatite-Based Catalysts. In Design and Applications of Hydroxyapatite-Based Catalysts; Wiley: Hoboken, NJ, USA, 2022; pp. 19–72. [Google Scholar]
  20. Schöller, K.; Ethirajan, A.; Zeller, A.; Landfester, K. Biomimetic Route to Calcium Phosphate Coated Polymeric Nanoparticles: Influence of Different Functional Groups and pH. Macromol. Chem. Phys. 2011, 212, 1165–1175. [Google Scholar] [CrossRef]
  21. Langroudi, M.M.; Saravani, M.G.; Nouri, A. Surfactant-Assisted Synthesis of Polyvinylpyrrolidone-Hydroxyapatite Composites as a Bone Filler. J. Appl. Appl. Biomater. Funct. Mater. 2017, 15, 334–340. [Google Scholar] [CrossRef]
  22. Guesmi, Y.; Agougui, H.; Jabli, M.; Alsharabasy, A.M. Bioactive composites of hydroxyapatite/polyvinylpyrrolidone for bone regeneration applications. Chem. Eng. Commun. 2018, 206, 279–288. [Google Scholar] [CrossRef]
  23. Chen, J.; Peng, Q.; Peng, X.; Zhang, H.; Zen, H. Probing and manipulating noncovalent interactions in functional polymeric systems. Chem. Rev. 2022, 122, 14594–14678. [Google Scholar] [CrossRef]
  24. Lehn, J.-M. Supramolecular Chemistry. Science 1993, 260, 1762–1764. [Google Scholar] [CrossRef] [PubMed]
  25. Riley, K.E.; Hobza, P. Noncovalent Interactions in Biochemistry. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2011, 1, 3–17. [Google Scholar] [CrossRef]
  26. Israelachvili, J.N. Intermolecular and Surface Forces; Academic Press: Cambridge, MA, USA, 2011; pp. 223–250. [Google Scholar]
  27. Murzakhanov, F.F.; Grishin, P.O.; Goldberg, M.A.; Yavkin, B.V.; Mamin, G.V.; Orlinskii, S.B.; Komlev, V.S. Radiation-induced stable radicals in calcium phosphates: Results of multifrequency EPR, EDNMR, ESEEM, and ENDOR studies. Appl. Sci. 2021, 11, 7727. [Google Scholar] [CrossRef]
  28. Shurtakova, D.V.; Yavkin, B.V.; Mamin, G.V.; Orlinskii, S.B.; Sirotinkin, V.P.; Fedotov, A.Y.; Shinkarev, A.; Antuzevics, A.; Smirnov, I.V.; Tovtin, V.I.; et al. X-ray diffraction and multifrequency EPR study of radiation-induced room temperature stable radicals in octacalcium phosphate. Radiat. Res. 2021, 195, 200–210. [Google Scholar] [CrossRef]
  29. Vorona, I.P.; Nosenko, V.V.; Baran, N.P.; Ishchenko, S.S.; Lemishko, S.V.; Zatovsky, I.V.; Strutynska, N.Y. EPR study of radiation-induced defects in carbonate-containing hydroxyapatite annealed at high temperature. Radiat. Meas. 2016, 87, 49–55. [Google Scholar] [CrossRef]
  30. Gilinskaya, L.G. Organic radicals in natural apatites according to EPR data: Potential genetic and paleoclimatic indicators. J. Struct. Chem. 2010, 51, 471–481. [Google Scholar] [CrossRef]
  31. Fisher, B.V.; Morgan, R.E.; Phillips, G.O.; Wardale, H.W. Radiation damage in calcium phosphates and collagen: An interpretation of ESR spectra. Radiat. Res. 1971, 46, 229–235. [Google Scholar] [CrossRef]
  32. Cevc, P.; Schara, M.; Ravnik, Č. Electron paramagnetic resonance study of irradiated tooth enamel. In Radiation Research; Radiation Research Society: Lawrence, UK, 1972; Volume 51, pp. 581–589. [Google Scholar]
  33. Vanhaelewyn, G.C.A.M.; Sadlo, J.; Matthys, P.F.A.E.; Callens, F.J. Comparative X-and Q-band EPR study of radiation-induced radicals in tooth enamel. Radiat. Res. 2002, 158, 615–625. [Google Scholar] [CrossRef]
  34. Vanhaelewyn, G.C.A.M.; Morent, R.A.; Callens, F.J.; Matthys, P.F.A.E. X-and Q-band electron paramagnetic resonance of CO2− in hydroxyapatite single crystals. Radiat. Res. 2000, 154, 467–472. [Google Scholar] [CrossRef] [PubMed]
  35. Forysenkova, A.A.; Ivanova, V.A.; Fadeeva, I.V.; Mamin, G.V.; Rau, J. 1H NMR and EPR Spectroscopies Investigation of Alginate Cross-Linking by Divalent Ions. Materials 2023, 16, 2832. [Google Scholar] [CrossRef] [PubMed]
  36. Gervais, C.; Bonhomme, C.; Laurencin, D. Recent directions in the solid-state NMR study of synthetic and natural calcium phosphates. Solid State Nucl. Magn. Reson. 2020, 107, 101663. [Google Scholar] [CrossRef]
  37. Soldati, E.; Rossi, F.; Vicente, J.; Guenoun, D.; Pithioux, M.; Iotti, S.; Malucelli, E.; Bendahan, D. Survey of MRI usefulness for the clinical assessment of bone microstructure. Int. J. Mol. Sci. 2021, 22, 2509. [Google Scholar] [CrossRef] [PubMed]
  38. Jäger, C.; Welzel, T.; Meyer-Zaika, W. EppleM A solid-state NMR investigation of the structure of nanocrystalline hydroxyapatite. Magn. Reson. Chem. 2006, 44, 573–580. [Google Scholar] [CrossRef] [PubMed]
  39. Manatunga, D.C.; de Silva, R.M.; Nalin de Silva, K.M.; de Silva, N.; Premalal EV, A. Metal and polymer-mediated synthesis of porous crystalline hydroxyapatite nanocomposites for environmental remediation. R. Soc. Open Sci. 2018, 5, 171557. [Google Scholar] [CrossRef]
  40. Hu, Y.-Y.; Liu, X.P.; Ma, X.; Rawal, A.; Prozorov, T.; Akinc, M.; Mallapragada, S.K.; Schmidt-Rohr, K. Biomimetic self-assembling copolymer− hydroxyapatite nanocomposites with the nanocrystal size controlled by citrate. Chem. Mater. 2011, 23, 2481–2490. [Google Scholar] [CrossRef]
  41. Fadeeva, I.V.; Trofimchuk, E.S.; Forysenkova, A.A.; Ahmed, A.I.; Gnezdilov, O.I.; Davydova, G.A.; Kozlova, S.G.; Antoniac, A.; Rau, J.V. Composite Polyvinylpyrrolidone–Sodium Alginate—Hydroxyapatite Hydrogel Films for Bone Repair and Wound Dressings Applications. Polymers 2021, 13, 3989. [Google Scholar] [CrossRef]
  42. Rai, R.K.; Sinha, N. Dehydration-induced structural changes in the collagen–hydroxyapatite interface in bone by high-resolution solid-state NMR spectroscopy. J. Phys. Chem. C 2011, 115, 14219–14227. [Google Scholar] [CrossRef]
  43. Gabbasov, B.; Gafurov, M.; Starshova, A.; Shurtakova, D.; Murzakhanov, F.; Mamin, G.; Orlinskii, S. Conventional, Pulsed and High-Field Electron Paramagnetic Reso¬nance for Studying Metal Impurities in Calcium Phosphates of Biogenic and Synthetic Origins. J. Magn. Magn. Mater. 2019, 470, 109–117. [Google Scholar] [CrossRef]
  44. Stoll, S. CW-EPR Spectral Simulations. Solid State Methods Enzymol. 2015, 563, 121–142. [Google Scholar]
  45. Vorona, I.P.; Ishchenko, S.S.; Baran, N.P.; Rudko, V.V.; Zatovskiĭ, I.V.; Gorodilova, N.A.; Povarchuk, V.Y. NO3 2-centers in synthetic hydroxyapatite. Phys. Solid State 2010, 52, 2364–2368. [Google Scholar] [CrossRef]
  46. Nosenko, V.V.; Vorona, I.P.; Ishchenko, S.S.; Baran, N.P.; Zatovsky, I.V.; Gorodilova, N.A.; Povarchuk, V.Y. Effect of pre-annealing on NO32-centers in synthetic hydroxyapatite. Radiat. Meas. 2012, 47, 970–973. [Google Scholar] [CrossRef]
  47. Fattibene, P.; Callens, F. EPR dosimetry with tooth enamel: A review. Appl. Radiat. Isot. 2010, 68, 2033–2116. [Google Scholar] [CrossRef] [PubMed]
  48. Radtsig, V.A. modern achievements in mechanical enzymative processing of plant raw materials. Khim. Fiz. 2004, 23, 70–109. (In Russian) [Google Scholar]
  49. Fadeeva, I.V.; Forysenkova, A.A.; Trofimchuk, E.S.; Gafurov, M.R.; Ahmed, A.I.; Davidova, G.A.; Antonova, O.S.; Barinov, S.M. Porous matrixes based on polyvinylpyrrolidone containing calcium phosphates for medical application. Russ. Chem. Bull. 2022, 71, 543–548. [Google Scholar] [CrossRef]
  50. Laurencin, D.; Almora-Barrios, N.; de Leeuw, N.H.; Gervais, C.; Bonhomme, C.; Mauri, F.; Chrzanowski, W.; Knowles, J.C.; Newport, R.J.; Wong, A.; et al. Magnesium incorporation into hydroxyapatite. Biomaterials 2011, 32, 1826–1837. [Google Scholar] [CrossRef]
  51. James, P.; Yesinowski Hellmut Eckert, J. Hydrogen Environments in Calcium Phosphates: NMR at High Spinning Speeds. Am. Chem. Soc. 1987, 109, 6274–6282. [Google Scholar]
Figure 1. Structure of HA.
Figure 1. Structure of HA.
Polymers 15 04445 g001
Figure 2. Schematic structure of PVP.
Figure 2. Schematic structure of PVP.
Polymers 15 04445 g002
Figure 3. SEM image of the PVP-HA composite.
Figure 3. SEM image of the PVP-HA composite.
Polymers 15 04445 g003
Figure 4. Comparison of EPR spectra of PVP-HA sample detected via ESE (blue upper curve) and cw (black lower curve) after X-ray irradiation in the X-band at T = 297 K. Middle curve was obtained by differentiating the ESE EPR spectra.
Figure 4. Comparison of EPR spectra of PVP-HA sample detected via ESE (blue upper curve) and cw (black lower curve) after X-ray irradiation in the X-band at T = 297 K. Middle curve was obtained by differentiating the ESE EPR spectra.
Polymers 15 04445 g004
Figure 5. EPR spectra obtained by detecting the ESE of samples after X-ray irradiation: X-band in the the left panel, W-band in the right panel. The red lines show the spectra of PVP, the black lines the spectra of HA, and the blue lines the spectra of PVP-HA. The dotted lines show the approximations of the EPR spectra. The approximation parameters are given in Table 1.
Figure 5. EPR spectra obtained by detecting the ESE of samples after X-ray irradiation: X-band in the the left panel, W-band in the right panel. The red lines show the spectra of PVP, the black lines the spectra of HA, and the blue lines the spectra of PVP-HA. The dotted lines show the approximations of the EPR spectra. The approximation parameters are given in Table 1.
Polymers 15 04445 g005
Figure 6. Photo-induced (using lasers with the wavelengths of 355 and 266 nm) and X-ray induced ESE EPR spectra at T = 200 K of PVP in the W-band. The dotted vertical line shows the signal from the quartz tube of the sample holder.
Figure 6. Photo-induced (using lasers with the wavelengths of 355 and 266 nm) and X-ray induced ESE EPR spectra at T = 200 K of PVP in the W-band. The dotted vertical line shows the signal from the quartz tube of the sample holder.
Polymers 15 04445 g006
Figure 7. 1H NMR spectra of PVP (blue curve) and PVP-HA (red curve) samples.
Figure 7. 1H NMR spectra of PVP (blue curve) and PVP-HA (red curve) samples.
Polymers 15 04445 g007
Figure 8. MAS-NMR spectra of 13C{1H} cp (red) and hpdec (blue) of the PVP. The rotation frequency of the sample is 7 kHz.
Figure 8. MAS-NMR spectra of 13C{1H} cp (red) and hpdec (blue) of the PVP. The rotation frequency of the sample is 7 kHz.
Polymers 15 04445 g008
Figure 9. Relaxation decay curves of protons of NMR signals for PVP (black squares) and PVP-HA (red circles for protons with chemical shift 4.7 ppm and blue triangles for the signal at −2.15 ppm).
Figure 9. Relaxation decay curves of protons of NMR signals for PVP (black squares) and PVP-HA (red circles for protons with chemical shift 4.7 ppm and blue triangles for the signal at −2.15 ppm).
Polymers 15 04445 g009
Figure 10. 31P NMR spectrum of PVP-HA. The sample rotation frequency is 7 kHz.
Figure 10. 31P NMR spectrum of PVP-HA. The sample rotation frequency is 7 kHz.
Polymers 15 04445 g010
Table 1. Spectroscopic parameters of the observed EPR spectra obtained from their simulations in the X- and W-bands.
Table 1. Spectroscopic parameters of the observed EPR spectra obtained from their simulations in the X- and W-bands.
Samplegg||A (MHz)A|| (MHz)ΔA (MHz)ΔA|| (MHz)
PVP2.00222.002638 ± 8106 ± 10--
HA2.00112.005292.4 ± 0.5186 ± 17 ± 112 ± 1
PVP-HA2.00112.005193.5 ± 0.5190 ± 213 ± 118 ± 1
Table 2. Spin–lattice (T1e) and spin–spin (T2e) relaxation times of NO32− radicals measured in the X-band (v = 9.6 GHz) at the central peak (transition) of the hyperfine structure of the nitrogen radical at B0 = 341.2 mT (corresponds to g).
Table 2. Spin–lattice (T1e) and spin–spin (T2e) relaxation times of NO32− radicals measured in the X-band (v = 9.6 GHz) at the central peak (transition) of the hyperfine structure of the nitrogen radical at B0 = 341.2 mT (corresponds to g).
T1e (μs)T2 (μs)
HA28.53
PVP-HA27.63.2
Table 3. The values of 1H spin–lattice relaxation rates for PVP and PVP-HA.
Table 3. The values of 1H spin–lattice relaxation rates for PVP and PVP-HA.
Sample2.4 ppm4.7 ppm−2.15 ppm
PVP0.67 ± 0.02 s−1--
PVP-HA-0.49 ± 0.02 s−10.55 ± 0.02 s−1
Table 4. The values of 1H spin–spin (transverse) relaxation times of components with populations pi, according to Equation (2).
Table 4. The values of 1H spin–spin (transverse) relaxation times of components with populations pi, according to Equation (2).
Chemical Shifts of NMR Signals (ppm)p1T2(1),
ms
p2T2(2),
ms
p3T2(3),
ms
PVP
2.40.001529.6 ± 0.50.00673.6 ± 0.10.9920.59 ± 0.02
PVP-HA
4.70.045.42 ± 0.20.181.67 ± 0.050.780.58 ± 0.02
−2.150.34629.2 ± 0.50.336.02 ± 0.050.3340.46 ± 0.03
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Petrova, A.; Mamin, G.; Gnezdilov, O.; Fadeeva, I.; Antonova, O.; Forysenkova, A.; Antoniac, I.V.; Rau, J.V.; Gafurov, M. Magnetic Resonance-Based Analytical Tools to Study Polyvinylpyrrolidone–Hydroxyapatite Composites. Polymers 2023, 15, 4445. https://doi.org/10.3390/polym15224445

AMA Style

Petrova A, Mamin G, Gnezdilov O, Fadeeva I, Antonova O, Forysenkova A, Antoniac IV, Rau JV, Gafurov M. Magnetic Resonance-Based Analytical Tools to Study Polyvinylpyrrolidone–Hydroxyapatite Composites. Polymers. 2023; 15(22):4445. https://doi.org/10.3390/polym15224445

Chicago/Turabian Style

Petrova, Alina, Georgy Mamin, Oleg Gnezdilov, Inna Fadeeva, Olga Antonova, Anna Forysenkova, Iulian V. Antoniac, Julietta V. Rau, and Marat Gafurov. 2023. "Magnetic Resonance-Based Analytical Tools to Study Polyvinylpyrrolidone–Hydroxyapatite Composites" Polymers 15, no. 22: 4445. https://doi.org/10.3390/polym15224445

APA Style

Petrova, A., Mamin, G., Gnezdilov, O., Fadeeva, I., Antonova, O., Forysenkova, A., Antoniac, I. V., Rau, J. V., & Gafurov, M. (2023). Magnetic Resonance-Based Analytical Tools to Study Polyvinylpyrrolidone–Hydroxyapatite Composites. Polymers, 15(22), 4445. https://doi.org/10.3390/polym15224445

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