Research Article  Open Access
Sondes Bouabdallah, Med Thaieb Ben Dhia, Med Rida Driss, "Study of a Conformational Equilibrium of Lisinopril by HPLC, NMR, and DFT", International Journal of Analytical Chemistry, vol. 2014, Article ID 494719, 8 pages, 2014. https://doi.org/10.1155/2014/494719
Study of a Conformational Equilibrium of Lisinopril by HPLC, NMR, and DFT
Abstract
The isomerization of lisinopril has been investigated using chromatographic, NMR spectroscopic, and theoretical calculations. The NMR data, particularly the NOEDIFF experiments, show that the major species that was eluted first is the trans form. The proportion was 77% and 23% for the trans and cis, respectively. The thermodynamic parameters (, , and ) were determined by varying the temperature in the NMR experiments. The interpretations of the experimental data were further supported by DFT/B3LYP calculations.
1. Introduction
Lisinopril, N(1carboxy3phenylpropyl)LlysylLproline, belongs to a class of antihypertensive agents which inhibit the angiotensinconverting enzyme (ACE) to control blood pressure [1]. The active parts of ACE inhibitors are peptide derivatives containing Cterminal proline residues. Like other prolinecontaining peptides, lisinopril exists as an equilibrium mixture of cis and trans isomers, with respect to the proline amide bond (Figure 1) [2, 3]. Under unstrained conditions most peptide bonds adopt the trans isomeric form, mainly because of the weaker steric repulsion between hydroxyl and carboxyl group effects in the molecule when compared to the cis. The trans form in lisinopril was shown to be the preferred isomer and biologically active [4–6]. The assignment separation of cis and trans form of lisinopril has been carried out by HPLC [5, 7–12], CZE [13–16], and NMR spectroscopy [2, 17–22]. The latter technique is a powerful tool and has been widely applied for structural and stereochemical characterization of amino acids, oligo and polypeptides [23–26]. The cistrans isomerization of peptide bonds is a slow process on the NMR time scale under normal conditions at ambient temperature due to the high barrier resulting from the C–N partial double bond character. NMR spectroscopy has therefore been successfully used to study the cistrans isomerization process of lisinopril in solution [25, 27].
In this paper, we report on the isomerization of lisinopril using a combination of HPLC, NMR spectroscopy, and theoretical approaches. The effect of temperature on the cistrans isomerization process of lisinopril was investigated in order to determine different thermodynamic parameters (, , and ).
2. Experimental
2.1. Samples
Lisinopril was kindly provided from Solvay Pharmaceuticals.
2.2. Reagents
Potassium dihydrogen phosphate, sodium hydroxide, and phosphoric acid were of RP quality from Prolabo (France). Methanol, acetonitrile, and tetrahydrofuran (THF) were of HPLC grade from LabScan (Dublin, Ireland).
The mobile phase was prepared by first preparing a solution of 0.02 M KH_{2}PO_{4}, adjusting its pH to 2 with phosphoric acid and finally mixing the solution with an organic modifier (acetonitrile, methanol, and THF). The mobile phases were always filtered using 0.45 μm membrane filter (Supelco) and degassed by sonication.
2.3. Chromatography
Liquid chromatographic analyses were performed using a Hewlett Packard 1100 HPLC system equipped with a photo diode array UV detector set at 215 nm. Injection was performed using an autoinjector. A Supelco LC 18 (5 μm) column ( mm I.D) and a guard column ( I.D) both from Supelco (Bellefonte, PA, USA) were used. The pH of mobile phase buffers was adjusted by means of a Schott model CG 825 pH meter (Germany).
2.4. Nuclear Magnetic Resonance (NMR Spectroscopy)
NMR spectra were obtained at 300.13 MHz on a Bruker Avance III spectrometer. The probe temperature was 298 K. chemical shifts were measured relative to tetramethylsilane [TMS, (CH_{3})_{4}Si]. Spectral width was 4201.68 Hz, acquisition time 2.818 s, numbers of scans 120, FID: TD 16384, SI 16384, LB 0.100, and relaxation time was 0.5 s.
All measurements were made on lisinopril in CD_{3}CN/D_{2}O (1/9) solution. The variable temperature NMR spectra were acquired using the same instrument. Probe temperatures (±0.5 K) were measured with a calibrated digital thermocouple. Samples were allowed to equilibrate for 10 min at each temperature before recording the spectrum.
2.5. Computational Details
Density functional theory (DFT) calculations were carried out on the cis and trans isomers of lisinopril with the Gaussian 03W suite of programs [28], with the nonlocal hybrid functional denoted as B3LYP [29]. Then basis sets used were zeta 631+G* [30–34], doubly polarized with diffuse functions on all the atoms. The geometries of both the cis and trans isomers were optimized using an analytical gradient. The harmonic vibration frequencies of the different stationary points of the potential energy surfaces (PES) have been calculated at the same level of theory in order to identify the local minima as well as to estimate a corresponding zeropoint vibrational energy (ZPE).
3. Results and Discussion
3.1. HPLC Study
The study of the cis/trans equilibrium of lisinopril by HPLC demonstrates that chromatographic conditions such as flow rate, temperature, pH, and organic modifier have an important effect on peak shape and retention time of lisinopril.
It appears that the separation of the two isomers of lisinopril can be achieved using a mobile phase consisting of a mixture of 20 mM phosphate buffer [pH 7]acetonitrile (90/10; v/v), a column temperature of 279 K, and flow rate of 2 mL/min with retention time min and min. However, a higher temperature is required for the elution of lisinopril as a single sharp peak at 2.76 min (Figure 2).
(a)
(b)
This is because it was found that an elevated temperature led to deterioration in the separation of the two isomers.
Moreover, at 328 K lisinopril was eluted as a narrow single peak due to the high isomerization rate of the two isomers. On the other hand, at low temperature the two isomers were resolved almost completely indicating that the interconversion rate had slowed down.
From ambient temperature chromatograms, the isomer trans/cis ratio was integrated to be 76/24. This result is similar to those reported earlier demonstrating that high temperature was useful for elution of prolinecontaining substances as a single peak [7, 11, 35]. Conversely, a low temperature is known to have a potential effect on the separation of isomers [5, 9, 36–39].
3.2. NMR Studies
The structure of lisinopril (Figure 3) shows 21 carbon atoms with two sets of two chemically equivalent carbons describing the ortho and metapositions on the aromatic ring. So, we expect to observe 19 signals in NMR spectra. However, the obtained spectra showed the doubling of all signals confirming the existence of the two isomers (Figure 4).
(a)
(b)
(a)
(b)
In addition, the NMR spectra of lisinopril in CD_{3}CN/D_{2}O at 298 K (Figure 5(a)) show two sets of triplets of unequal intensities. The multiplicity of each signal set reflected first from the interaction of H58 with H25 and H26, giving the two signals in the 3.8–4.1 ppm region, and second from the interaction of H43 with H45 in the 4.1–4.4 ppm region. The same spectrum recorded at 333 K (Figure 5(b)) shows a better separation of the two signals at 4.1 ppm.
(a)
(b)
These isomers are assigned to a cistrans equilibrium of the rotation around the amide bond. As described earlier, it is worth noting that this equilibrium appears to be slow on the NMR time scale at ambient temperatures [37, 40]. Using the area of resonance signals of proton 58 (3.8–4.4 ppm), the isomer ratio was integrated to be 77/23 at 298 K. The result obtained in a separate experiment recorded at a probe temperature of 298 K is consistent with that determined by HPLC at the same temperature.
We conclude that the major conformer in the NMR spectrum of lisinopril corresponds to the first eluted peak in the HPLC chromatogram at ambient temperature, which exists in a higher proportion. A similar study demonstrated this correspondence in the case of ramiprilat [4], enalaprilat [5], and perindopril [39] in different solvents.
It is well known that the NOE effect () between two dipoledipole interacting nuclei is inversely proportional to the distance between the irradiation site () and the measured one (), respectively [41], according to the following formula: .
Therefore, a low NOE is observed at H_{58} ( ppm) when H_{43} ( ppm) is irradiated in the major conformer. Accordingly, a stronger NOE () is noted at H_{58} ( ppm) when H_{43} ( ppm) is irradiated in the minor conformer.
In addition, when H_{48} ( ppm) is irradiated, we observe at H_{58} a NOE () in the major conformer and a NOE () in the minor conformer (Figure 6).
(a)
(b)
(c)
(d)
This shows that the NOE in the minor conformer for the H43/H58 is more important than the NOE in the major conformer, which implies that in the major conformer, the distance between the nuclei is higher than the one in the minor conformer
Based on the relationship between NOE and internuclear distances, one can give the expression of distance in each conformer
Consequently, the examination of the molecular structure of each conformer of lisinopril confirmed that the distance in the cis conformer is indeed smaller than the distance in the trans form. The NMR intensities suggest that the major conformer is the trans form and the minor conformer is the cis form. The nuclei of the strans conformer are more deshielded than those of scis conformer, in agreement with literature results [41, 42].
3.3. Thermodynamic Study
At slow chemical exchange, the relative proportion of the two conformers at different temperatures in the range (279–333) K and the cis/trans equilibrium constant of lisinopril have been measured by relative integrals of the two resonance signals of the two states of isomerizations of lisinopril.
These signals are well resolved and allowed to determine accurately the equilibrium constants for the cis to trans interconversion at different temperatures and to measure the thermodynamic parameters: enthalpy (), entropy () of the equilibrium on the basis of the van’t Hoff equation. The Gibbs free enthalpy () is deduced at ambient temperature.
The plot of versus the reciprocal of the absolute temperature is a straight line of equation ln (Figure 7). The correlation coefficient for this straight line is . The enthalpy was obtained via the slope and the entropy via the intercept of plot. The thermodynamic parameters obtained from experiment were kJ/mol, J/kmol, and kJ/mol.
Remarkably, the equilibrium cis to trans isomerization was enthalpically and entropically favored in this condition. Consequently, the decrease in the temperature expected a displacement of the conformational equilibrium to the strans conformer. The latter is stabilized by hydrogen bonding between the carbonyl (40) and the hydrogen of hydroxyl group OH of the acid function (56). This result is in agreement with other studies reported on a similar product such as enalapril [5].
3.4. Theoretical Calculations
In order to confirm the NMR data obtained for the cis/trans isomerization of lisinopril, the geometrics of the two conformers were fully optimized at the DFT/B3LYP level of theory using 631++G* basis set. The structures have been identified as local minima on the singlet potential energy surfaces (PES) (Figure 8). Optimized values of selected geometrical parameters are listed in Table 1. The potential energy difference between the two isomers of lisinopril was 11.397 kJ/mol indicating the stability of trans over the cis isomer. This difference is in good agreement with experimental data that the trans is the majorities form.

As shown in Table 1, in the trans conformer the interaction between H56O40 ( ) would be stronger than that in the cis ( ), thus generating a sharp reduction of the valence angle C53O55H56 ( in the trans form and in the cis form) and a strong variation on dihedral angle O40C39C41C50 ( in the trans form and 2.6° in the cis form), indicating possible existence of some hydrogen bond interactions which would be more favored in the trans than in cis isomers.
This stability of trans form over the cis form is further confirmed by the charge density between the same atoms (, ). The strong interaction between O21H56 atoms ( ) in the scis conformer triggers a modification of the dihedral angle O21H56O55C53 () which may explain the low stability of the cis form. On the other hand, the low interaction between H56O21 ( ) in the strans conformer yields a change in the dihedral angle () so the atoms H56 and O40 were far from each other and the trans was more stable.
It is worth noting that there are two hydrogen bonds, one between atoms H56O40 and another between H55O21, but the stability was determined by the first bond because the distance is shorter ( ). In addition, there are strong intramolecular interactions in the trans conformer ( D) than in the cis form ( D) indicating a higher stability of the conformer trans.
The distance between the nuclei in both conformers (cis and trans) was compared. It is revealed that in the cis conformer, H43 is close to H58 (4.058 ), while H48 is at a much greater distance from H58 (6.612 ). On the other hand, in the trans conformer H43 is further away from H58 (4.471 ), while H48 is at much greater distance from H58 (6.225 ). This is consistent with the results obtained with NOE difference experiments giving an enhancement of 8% and 21% in the first irradiation and 34% and 10% in the second irradiation.
4. Conclusions
Lisinopril exists individually as a mixture of cistrans isomers in solution. The two isomers could be easily distinguished by HPLC, , , and ^{1}H^{1}H NOE spectra. The results of the various investigations by NMR and those of the theoretical study are in favor of assigning trans form to the major isomer present under the conditions used in the HPLC studies. Additionally, a relative to successful combination of molecular modeling studies with experimental spectroscopic assays was used in order to elucidate the molecular bases.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The authors want to express their special thanks to Dr. M. A. M. K. Sanhoury of the Faculty of Sciences of Tunis for his valuable help and discussion of this work. They also thank the editors and anonymous referees for valuable comments and suggestions that greatly improved the paper.
References
 E. Mutschler, H. Derendorf, K. SchhaferKorting, K. Elord, and S. Esters, Drug Actions. Basic Principals and Therapeutic Aspects, pp. 385–390, Medpharm Scientific Publishers, Stuttgart, Germany, 1995.
 F. Thunecke, A. KLalman, F. Kalmen, S. Ma, A. S. Rathore, and C. Horvath, “Kinetic study on the cistrans isomerization of peptidylproline dipeptides,” Journal of Chromatography A, vol. 744, no. 12, pp. 259–272, 1996. View at: Publisher Site  Google Scholar
 J. L. Kofron, P. Kuzmic, E. ColonBonilla, and D. H. Rich, “Determination of kinetic constants for peptidyl prolyl cistrans isomerases by an improved spectrophotometric assay,” Biochemistry, vol. 30, no. 6127, p. 6134, 1991. View at: Google Scholar
 A. Skoglof, I. Nilson, S. Gustaffsson, J. Deinum, and P. Gtithe, “Cistrans isomerization of an angiotensin Iconverting enzyme inhibitor. An enzyme kinetic and nuclear magnetic resonance study,” Biochimica et Biophysica Acta, vol. 1041, no. 1, pp. 22–30, 1990. View at: Publisher Site  Google Scholar
 S. Bouabdallah, H. Trabelsi, T. B. Dhia, S. Sabbah, K. Bouzouita, and R. Khaddar, “RPHPLC and NMR study of cistrans isomerization of enalaprilat,” Journal of Pharmaceutical and Biomedical Analysis, vol. 31, no. 4, pp. 731–741, 2003. View at: Publisher Site  Google Scholar
 R. M. Veronika and G. S. Peter, “HPLCNMR of the conformers of captopril,” in Proceedings of the 23rd International Symposium on High Performance Liquid Phase Separations and Related Techniques (HPLC '99), PB12/46, March 1999. View at: Google Scholar
 W. R. Melander, J. Jacobson, and C. Horváth, “Effect of molecular structure and conformational change of prolinecontaining dipeptides in reversedphase chromatography,” Journal of Chromatography A, vol. 234, no. 2, pp. 269–276, 1982. View at: Google Scholar
 J. Jacobson, W. Melander, G. Vaisnys, and C. Horváth, “Kinetic study on cistrans proline isomerization by highperformance liquid chromatography,” Journal of Physical Chemistry, vol. 88, no. 20, pp. 4536–4542, 1984. View at: Google Scholar
 S. Gustafsson, B.M. Eriksson, and I. Nilsson, “Multiple peak formation in reversedphase liquid chromatography of ramipril and ramiprilate,” Journal of Chromatography A, vol. 506, pp. 75–83, 1990. View at: Publisher Site  Google Scholar
 R. Hana and A. Wada, “Analysis of the cistrans isomerization kinetics of lalanyllproline by the elutionband relaxation method,” Journal of Chromatography A, vol. 394, pp. 273–278, 1987. View at: Publisher Site  Google Scholar
 S. Bouabdallah, H. Trabelsi, K. Bouzouita, and S. Sabbah, “Reversedphase liquid chromatography of lisinopril conformers,” Journal of Biochemical and Biophysical Methods, vol. 54, pp. 391–405, 2002. View at: Publisher Site  Google Scholar
 A. Shoji, A. Yanagida, H. Shindo, Y. Ito, and Y. Shibusawa, “Countercurrent chromatographic estimation of hydrophobicity of Z(cis) and E(trans) enalapril and kinetics of cistrans isomerization,” Journal of Chromatography A, vol. 1157, no. 12, pp. 101–107, 2007. View at: Publisher Site  Google Scholar
 X. Z. Qin, P. Ip, and W. E. Tsai, “Determination and rotamer separation of enalapril maleate by capillary electrophoresis,” Journal of Chromatography A, vol. 6626, no. 2, pp. 251–258, 1992. View at: Publisher Site  Google Scholar
 A. Moore and J. Jorgenson, “Resolution of cis and trans isomers of peptides containing proline using capillary zone electrophoresis,” Journal of Analytical Chemistry, vol. 67, no. 19, pp. 3464–3475, 1995. View at: Publisher Site  Google Scholar
 S. Ma, F. Kalman, A. Kalman, F. Thunecke, and C. Horvath, “Capillary zone electrophoresis at subzero temperatures. I. Separation of the cis and trans conformers of small peptides,” Journal of Chromatography A, vol. 716, no. 12, pp. 167–182, 1995. View at: Publisher Site  Google Scholar
 A. S. Rathore and C. Horváth, “Capillary zone electrophoresis of interconverting cistrans conformers of peptidylproline dipeptides: estimation of the kinetic parameters,” Electrophoresis, vol. 18, no. 15, pp. 2935–2943, 1997. View at: Publisher Site  Google Scholar
 A. Christopher, C. A. Evans, and D. L. Rabenstein, “Nuclear magnetic resonance studies of the acidbase chemistry of amino acids and peptides: II. Dependence of acidity of the Cterminal peptide bond,” Journal of the American Chemical Society, vol. 96, pp. 7311–7315, 1974. View at: Google Scholar
 R. E. London, N. A. Matwiyoff, J. M. Stewart, and J. R. Cann, “^{13}C nuclear magnetic resonance study of the cistrans isomerism in XProPro tripeptides,” The Journal of Biochemistry, vol. 17, no. 12, pp. 2277–2283, 1978. View at: Publisher Site  Google Scholar
 Z. P. Liu and L. M. Gierash, “Combined use of molecular dynamics simulations and NMR to explore peptide bond isomerization and multiple intramolecular hydrogenbonding possibilities in a cyclic pentapeptide, cyclo (GlyProDPheGlyVal),” Journal of Biopolymers, vol. 32, no. 12, pp. 1727–1739, 1992. View at: Publisher Site  Google Scholar
 J. Poznanski, A. Ejchart, K. L. Wierzchowski, and M. Ciurak, “^{1}H and ^{13}Cnmr investigations on cistrans isomerization of proline peptide bonds and conformation of aromatic side chains in HTrp(Pro)_{n}TyrOH peptides,” Journal of Biopolymers, vol. 33, no. 5, pp. 781–759, 1993. View at: Publisher Site  Google Scholar
 Y. Sakamoto and T. Ishi, “Conformational studies by ^{1}H and ^{13}C NMR of lisinopril,” Journal of Molecular Structure, vol. 198, pp. 129–136, 1993. View at: Publisher Site  Google Scholar
 A. Saran, S. Srivastava, and E. Coutinho, “Conformation and macromolecular interaction of ACE inhibitors enalapril and lisinopril by NMR,” European Journal of Pharmaceutical Sciences, vol. 1, no. 4, pp. 181–188, 1994. View at: Google Scholar
 J. N. S. Evans, Biomolecular NMR Spectroscopy, Oxford University Press, 1995.
 Q. Teng, Structural Biology: Practical NMR Applications, Springer, 2005.
 G. Ivanova, B. Yakimova, S. Angelova, I. Stoineva, and V. Enchev, “Influence of pH on the cistrans isomerization of ValineProline dipeptide: an integrated NMR and theoretical investigation,” Journal of Molecular Structure, vol. 975, pp. 330–334, 2010. View at: Publisher Site  Google Scholar
 A. A. Isab and M. I. M. Wazeer, “Solid and solution NMR studies of the complexation of Ag^{+} with the trans isomer of captopril: biological activities of this high blood pressure drug along with its Ag^{+} complex,” Spectrochimica Acta A, vol. 65, no. 1, pp. 191–195, 2006. View at: Publisher Site  Google Scholar
 C. A. Evans and D. L. Rabenstein, “Nuclear magnetic resonance studies of the acidbase chemistry of amino acids and peptides. II. Dependence of the acidity of the Cterminal carboxyl group on the conformation of the Cterminal peptide bond,” Journal of the American Chemical Society, vol. 96, no. 23, pp. 7312–7317, 1974. View at: Publisher Site  Google Scholar
 J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian 03, Revision A. 1, Gaussian, Inc., Pittsburgh, Pa, USA, 2003.
 A. D. Becke, “Densityfunctional thermochemistry. III. The role of exact exchange,” The Journal of Chemical Physics, vol. 98, no. 7, p. 5648, 1993. View at: Publisher Site  Google Scholar
 W. J. Hehere, R. Ditchfield, and J. A. Pople, “Selfconsistent molecularorbital methods XII. Further extension of Gaussiantype basis sets for use in molecular orbital studies of organic molecules,” The Journal of Chemical Physics, vol. 56, no. 5, p. 2257, 1972. View at: Publisher Site  Google Scholar
 G. A. Petersson, A. Bennett, T. G. Tensfeldt, M. A. AlLaham, W. A. Shirley, and J. Mantzaris, “A complete basis set model chemistry. I. The total energies of closedshell atoms and hydrides of the firstrow elements,” The Journal of Chemical Physics, vol. 89, no. 4, p. 2193, 1988. View at: Publisher Site  Google Scholar
 G. A. Petersson and M. A. AlLaham, “A complete basis set model chemistry. II. Openshell systems and the total energies of the firstrow atoms,” The Journal of Chemical Physics, vol. 94, no. 9, p. 6081, 1991. View at: Publisher Site  Google Scholar
 P. C. Hariharan and J. A. Pople, “The influence of polarization functions on molecular orbital hydrogenation energies,” Theoretica Chimica Acta, vol. 28, no. 3, pp. 213–222, 1973. View at: Publisher Site  Google Scholar
 V. A. Rassolov, M. A. Ratner, J. A. Pople, P. C. Redfern, and L. A. Curtiss, “631G* basis set for thirdrow atoms,” Journal of Computational Chemistry, vol. 22, no. 9, pp. 976–984, 2001. View at: Publisher Site  Google Scholar
 D. E. Henderson and C. Horvath, “Low temperature highperformance liquid chromatography of cistrans proline dipeptides,” Journal of Chromatography A, vol. 368, pp. 203–213, 1986. View at: Publisher Site  Google Scholar
 D. E. Henderson and J. A. Mello, “Physicochemical studies of biologically active peptides by lowtemperature reversedphase highperformance liquid chromatography,” Journal of Chromatography A, vol. 499, pp. 79–88, 1990. View at: Publisher Site  Google Scholar
 P. K. Owens, L. A. Svensson, and J. Vessman, “Direct separation of captopril diastereoisomers including their rotational isomers by RPLC using a teicoplanin column,” Journal of Pharmaceutical and Biomedical Analysis, vol. 25, no. 34, pp. 453–464, 2001. View at: Publisher Site  Google Scholar
 H. Trabelsi, S. Bouabdallah, S. Sabbah, F. Raouafi, and K. Bouzouita, “Study of the cistrans isomerization of enalapril by reversedphase liquid chromatography,” Journal of Chromatography A, vol. 871, no. 12, pp. 189–199, 2000. View at: Publisher Site  Google Scholar
 S. Bouabdallah, H. Trabelsi, M. T. B. Dhia, and N. B. Hamida, “Kinetic study on the isomerization of perindopril by HPLC,” Chromatographia, vol. 75, no. 2122, pp. 1247–1255, 2012. View at: Publisher Site  Google Scholar
 F. A. Bovey, J. J. Ryan, and F. P. Hood, “Polymer nuclear magnetic resonance spectroscopy. XV. The conformation of polysarcosine,” Macromolecules, vol. 1, no. 4, pp. 305–307, 1968. View at: Google Scholar
 P. A. Hart and J. P. Davis, “Conformation of nucleosides and nucleotides. Application of the nuclear Overhauser effect,” Journal of the American Chemical Society, vol. 91, no. 2, pp. 512–513, 1969. View at: Publisher Site  Google Scholar
 P. E. Shirmer, J. H. Noggle, J. P. Davis, and P. A. Hart, “Determination of molecular geometry by quantitative application of the nuclear Overhauser effect,” Journal of the American Chemical Society, vol. 92, no. 11, pp. 3266–3273, 1970. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Sondes Bouabdallah et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.