Next Article in Journal
The Role of Platelets in Infective Endocarditis
Next Article in Special Issue
The Effect of Non-Solvent Nature on the Rheological Properties of Cellulose Solution in Diluted Ionic Liquid and Performance of Nanofiltration Membranes
Previous Article in Journal
Characterization of Porcine Ocular Surface Epithelial Microenvironment
Previous Article in Special Issue
A Novel D-Psicose 3-Epimerase from Halophilic, Anaerobic Iocasia fonsfrigidae and Its Application in Coconut Water
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Substitution Degree and Homogeneity on Cyclodextrin-Ligand Complex Stability: Comparison of Fenbufen and Fenoprofen Using CD and NMR Spectroscopy

1
Department of Pharmaceutical Chemistry, Semmelweis University, Hőgyes Endre utca 9, 1092 Budapest, Hungary
2
ELKH-ELTE Protein Modelling Research Group, Eötvös Loránd University, Pázmány Péter sétány 1A, 1117 Budapest, Hungary
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(8), 7544; https://doi.org/10.3390/ijms24087544
Submission received: 31 March 2023 / Revised: 13 April 2023 / Accepted: 18 April 2023 / Published: 19 April 2023

Abstract

:
The stability of host–guest complexes of two NSAID drugs with similar physicochemical properties, fenbufen and fenoprofen, was investigated by comparing induced circular dichroism and 1H nuclear magnetic resonance methods using eight cyclodextrins of different degrees of substitution and isomeric purity as guest compounds. These cyclodextrins include native β-cyclodextrin (BCyD), 2,6-dimethyl-β-cyclodextrin 50 (DIMEB50), 80 (DIMEB80) and 95% (DIMEB95) isomerically pure versions, low-methylated CRYSMEB, randomly methylated β-cyclodextrin (RAMEB) and 4.5 and 6.3 average substitution grade hydroxypropyl-β-cyclodextrin (HPBCyD). The stability constants obtained by the two methods show good agreement in most cases. For fenbufen complexes, there is a clear trend that the stability constant increases with the degree of substitution while isomer purity has a smaller effect on the magnitude of stability constants. A significant difference was found in the case of DIMEB50 when compared to DIMEB80/DIMEB95, while the latter two are similar. In the fenbufen–fenoprofen comparison, fenbufen, with its linear axis, gives a more stable complex, while fenoprofen shows lower constants and poorly defined trends.

1. Introduction

Cyclodextrins (CyDs) are used in many areas of the pharmaceutical and food industry [1,2,3,4,5,6], some of their most important properties being that they can enhance the dissolution of poorly water-soluble substances, reduce the volatility of volatile molecules, and increase the chemical stability of degradable substances by inclusion complexation. As an indication of the importance of CyDs, the number of publications on CyDs in the Scopus database is increasing year by year. In pharmaceutical applications, they can also influence the pharmacokinetics of active substances and can be used as processing excipients and orphan drugs [7]. Recently, a review of NSAID active substances was published [8] in which the literature of about 60 different inclusion complexes of 24 active substances and 12 types of CyDs was reviewed between 2010 and 2020. CyDs are also often used as chiral selectors in the analysis of chiral molecules (enantiomers) in separation processes [9]. Chiral differentiation and structure elucidation of cyclodextrins have been extensively studied in the literature [10,11,12,13,14,15,16]. In addition to the natural (α, β, γ) CyDs, a number of substituted derivatives are known that affect the properties of the original molecule.
The substituents may be neutral (methyl, hydroxypropyl), positively charged (e.g., amino-, guanidino- [16]), or negatively charged (carboxymethyl, sulfobutyl-ether, sulfated [17]) derivatives, depending on the intended purpose. In the case of substituted derivatives with acid–base active functional groups, the average protonation state is pH-dependent and the determination of the pKa values for each substituent is not straightforward [18]. Substitution can be made on the OH groups 2, 3, and 6 of α-D-glucopyranose units, resulting in mono-, di-, or trisubstituted derivatives. The wider rim of the CyD ring contains the secondary OH groups at position 2 and 3, while the narrower rim of the truncated cone-shaped ring system contains the primary OH groups at position 6. Thus, a β-cyclodextrin (BCyD) with 1–4 glycosidic OH bonds from seven sugar subunits can theoretically have 7 substituents on the molecule in a monosubstituted, 14 in a disubstituted, and 21 in a trisubstituted state. The selectivity of the substitution reaction strongly influences the actual substitution site and the number of substituents [19,20]. It is therefore important that the average number of substituents (DS—Degree of Substitution) and the isomeric purity of the resulting product are determined by appropriate analytical methods [21]. In the period of 2018–2022, Scopus returned only 136 publications for “single isomer” and “cyclodextrin” out of nearly 100,000 hits related to cyclodextrins.
Several methods are available for the qualitative and quantitative characterization of CyD complexes in both solution and solid phases. The most important features are the stoichiometry of the complex (C), i.e., the host (H) to guest (G) ratio, and the stability constant (logK). In our previous work, the stability constants of CyD complexes of antifungal azoles were determined by circular dichroism (CD) spectroscopy and the results were compared with literature data [22] in Table 1. A relatively good agreement was obtained for the native BCyD, but quite significant differences were found for the substituted derivatives.
Differences between stability constants determined by different methods may be due to: (1) in the comparison, the degree of substitution and isomeric composition of the substituted CyD used are not precisely defined; (2) in solution phase experiments, different experimental conditions may cause discrepancies, e.g., temperature, pH, ionic strength, or the presence of other components.
Recognizing these differences, we set out to determine the stability constants of CyD derivatives with different degrees of substitution and isomeric purity for two NSAIDs (fenbufen and fenoprofen depicted in Figure 1) with similar physicochemical properties by CD and 1H NMR spectroscopy. The model compounds were selected according to the following criteria. On the one hand, most NSAIDs are poorly soluble in water. This can lead to problems with bioavailability. An important role of cyclodextrins is to improve the solubility of poorly soluble compounds. Several formulations containing NSAID-cyclodextrin complexes are currently on the market, e.g., Meloxicam/Mobitil; Nimesulide/Nimedex; Piroxicam/Cycladol or Brexin; Diclofenac/Voltaren Ophtalmic, and in this way, the study may have practical relevance in the future. On the other hand, the ICD signal requires the guest molecule to have a chromophore group, a condition that is met by both compounds. The main consideration in the design of the experiments was to measure the same solutions with the two orthogonal methods. Bratu et al., investigated the BCyD-fenbufen system by NMR, and determined the stoichiometry, structure, and stability constant. They reported a 1:1 composition and a value of “several 1000 s (M−1)” [30]. Another literature value could be found (4.72 × 103) in the publication of Zhang et al. [31] where mass spectrometry was applied to determine the stability of the complex. Furthermore, UV-, induced CD-, and fluorescence spectroscopic methods were used by Sortino et al. [32] to determine the BCyD-fenbufen stability constant. They determined K values of 3200 M−1 by UV, 1690 and 2200 M−1 at two different wavelengths by induced CD and 2690 M−1 by fluorescence. The stability of a pure S-isomer of the BCyD-fenoprofen complex was determined by 1H NMR spectroscopy [33]. Uccello-Barretta et al., found a 1:1 stoichiometry and a K value of 2010 M−1. They calculated the result from the chemical shift of the methyl proton.
1H NMR spectroscopy can be used to monitor complexation as the signals of the guest and CyD are usually visible. The NMR spectra of the investigated guest molecules and BCyD together with their assignments are shown in Figure 2.
The change in the chemical shift of the guest signals as complexation occurs (CIS—Complexation Induced Shift) can be used to determine the stability constant.
δobs = δG∙χG + δC∙χC
As shown in Equation (1), the observed chemical shift (δobs) is the average of the chemical shift of the free guest (δG) and the complexed (δC) species weighted by their mole fraction (χ). In addition to the stability constants, the stoichiometry of the complex can also be estimated [34,35] and information on the structure of the complex can be obtained from the NOESY or ROESY spectra, if cross peaks are observed between the guest signals and the CyD H3 and H5 are found inside the cavity. The limitation of NMR is a low signal-to-noise ratio and an inability to discern between different CyD isomers. In case of native cyclodextrins—due to their symmetry—the interpretation of the NMR spectra can be simple, but in the case of modified cyclodextrins, the presence of structural isomers (different substitution degree and position) complicates the structural analysis. Although CD spectroscopy yields less information in terms of exact atomic interactions, it has superior selectivity towards the complex formation, and, as such, is one of the most effective methods for CyD complex analysis [36]. The basis of the CD selectivity is the fact that achiral guests have no CD spectra while CD spectra of CyD is insignificant above 220 nm. Therefore, the requirement for an induced CD spectrum (ICD), which solely originates from the complex, is a chromophore (of the guest) that is under chiral perturbation by the CyD [37]. The appearance of this ICD signal in the region of the absorption bands is a definite proof of the complex formation. A further requirement is that the dipole moment of the electron transition of the chromophore and the axis of the cyclodextrin have a favorable geometric alignment since the rotator strength is derived from these two factors. It may be that the geometry of the complex is so unfavorable that no ICD signal appears. The sign and intensity of the ICD signal can provide information regarding the structure of the complex [38,39,40]. The value of the stability constant can be obtained from logistic non-linear regression on the ICD ~ host concentration data (with constant guest concentration), as described in our previous study [22]. A limitation of CD spectroscopy is the requirement for the presence of a chromophore on the guest and the low intensity of ICD signals when absorption is relatively high. This latter effect is a result of the absorption difference (Δε) of the two circularly polarized light components being proportional to the CD signal, while the sum of the component absorptions is equal to the total absorption. Another disadvantage is that the intensity of the ICD signal is usually an order of magnitude lower than the normal CD signal, and for racemic compounds, only an apparent stability constant can be calculated for the two enantiomers. The intensity of the ICD signal can be influenced by the dynamics of the complex and by possible hydrogen bond interactions.
The use of cyclodextrins in pharmaceutical research and development is becoming increasingly frequent. Accurate characterization of cyclodextrin complex stability constants is essential for correct drug formulations. Of the analytical methods commonly used to characterize cyclodextrin complexation, two well-known methods are CD and NMR, but in some cases (see Table 1), a considerable difference has been observed between the stability constants obtained by the two methods. The purpose of this study was to investigate whether this difference is relevant and systematic. We chose to perform this investigation on two selected guest compounds and a variety of CyD hosts.

2. Results

2.1. Fenbufen Complex Stability Constants

Fenbufen complexation with the various CyD hosts was analyzed using CD and 1H NMR spectroscopy. The ICD spectra of HPBCyD(4.5)-fenbufen are shown in Figure 3, along with the 1H NMR spectra series. Figure 3 also depicts the non-linear regression plots on the CD and NMR data. The CD and NMR data with the other CyDs can be found in Supplementary Tables S1–S11. The determined complex stability constant values in log10 units are summarized in Table 2, together with the difference between the values obtained with the two methods. Furthermore, in Table 3, the non-linear regression results from each individual NMR signals are compiled, where the logK* is handled as a random effects model parameter as opposed to the logK fitted globally to all NMR signals as a fixed effects model parameter presented in Table 2.

2.2. Fenoprofen Complex Stability Constants

Fenoprofen complexation was studied with the same methods as for fenbufen. The ICD and 1H NMR spectra of CyD-fenoprofen complexes are shown in Figure 4 along with select non-linear regression plots. The full set of CD and NMR data are found in Supplementary Tables S12–S23. The determined complex stability constant values (and the difference between the methods) in log10 units are compiled in Table 4. In Table 5, the non-linear regression results from each individual NMR signals are tabulated.
The BCyD-fenoprofen complex was also studied with ROESY NMR in order to elucidate the binding points in the complex. Figure 5 shows the ROESY cross peaks between the H3, H5, and H6 hydrogens of BCyD and the aromatic hydrogens of fenoprofen.

3. Discussion

The combined CD and NMR analysis of CyD-fenbufen and CyD-fenoprofen complexes revealed that in the case of both complexes, the aromatic rings of the guest species entered the hydrophobic cavity of the cyclodextrin hosts. Fenbufen showed similar ICD spectra for all CyD hosts, as seen in Figure 3, with a characteristic spectra series of HPBCyD(4.5) complex. For the fenbufen-cyclodextrin complex, a relatively high intensity ICD signal with a positive sign was obtained for all cyclodextrins. Thus, based on the Kodaka–Harata rules [38,40] it can be assumed that the temporary dipole vector of the excitation transition and the axis of the cyclodextrin cavity coincide, and the chromophore of the guest molecule is predominantly located in the cyclodextrin cavity.
On the other hand, fenoprofen had different ICD profiles for the various cyclodextrins (Figure 4). Only the BCyD complex had higher intensity positive bands, whereas CRYSMEB had only about half the intensity. These ICD bands still suggest that the chromophore is predominantly located in the cyclodextrin cavity, but in the case of CRYSMEB the methyl groups at the rim already inhibit the complete immersion. For the other cyclodextrins, the lower intensity bands make it more difficult to predict the structure. The substituents on the cyclodextrin rim may be spatially inhibited, resulting in less immersion. In the absence of quantum chemical calculations, the temporary dipole vector of the excitation transition is also difficult to estimate but is likely to be non-parallel to the axis of the cyclodextrin cavity. It is also difficult to describe only one structure based on 2D NMR (ROESY) spectra, since the H2, H3, and H6 hydrogens of BCyD gave rise to several cross-peaks with the aromatic protons of fenoprofen (Figure 5). NMR and CD results suggest that the ICD signal is lower in intensity due to the complex dynamics and that several structures may exist simultaneously. The spectral differences between fenbufen- and fenoprofen-CyD complexes are probably due to fenbufen having a rather linear shape and thus not having many degrees of freedom when the biphenyl part of the molecule enters the CyD cavity, while fenoprofen can enter cyclodextrins in various ways as a result of the rotating etheric oxygen bridge.
When comparing the complex stability constants obtained from NMR measurements, one can perform a non-linear regression analysis on the data of the NMR signals with the logK as a random effects model parameter, i.e., allowing different logK values for the regression fit of individual NMR peak data (Table 3 and Table 5). This sometimes leads to considerable differences (or even the fit did not converge for a given NMR signal) in the obtained logK values for the same complex. This can be interpreted mainly as an error effect of certain NMR signals since the deviance of logK values obtained from individual NMR peaks from a globally fitted logK showed a direct relationship with the response effect (maximal chemical shift response) of that individual NMR peak (Figure 6). That is, the greater the response of an NMR signal to complexation, the greater the reliability of logK obtained from that signal. It could also be argued that different NMR signals are sensitive to different aspects of the complexation; therefore they might convey different information and thus different apparent logK values are obtained. It is therefore much more reliable to perform the regression analysis of the NMR data globally, i.e., holding the logK as a fixed effects model parameter only and as a result the same logK value is obtained from the fit of all NMR signals. This globally fitted logK will have a lower standard error and a better agreement with the regression result of CD data. It may be argued therefore, that CD measurements represent the complexation, as a whole since the chromophore of the guest is central to the binding processes (H-bond formation, hydrophobic interactions) and is therefore perturbed by all aspects of complexation.
The complex stability constants obtained from the two methods (CD and NMR global fit) can be compared with a Bland–Altman plot (Figure 6), where the mean of logKCD and logKNMR values obtained for a particular complex are plotted against the differences (logKCD − logKNMR). From this scatter plot one can evaluate the deviation of mean differences (i.e., bias) from the null value and also the 95% confidence interval of the differences gives us an idea on whether the methods differ considerably in variability inherent to the measurement.
It can be seen from the Bland–Altman plot that the 95% confidence interval of the differences falls within 0.2 logK units; given the uncertainty of the obtained data from CD measurements and global fits, it can be stated that this difference between the two methods is not considerable, i.e., methodologically not relevant. There is also no tendency in the scatter of the points in the Bland–Altman plot and no considerable systematic error (bias). The small differences between the two methods can thus be explained with measurement error (differences in the sensitivity of the two methods). The larger difference found in the literature review (Table 1) can therefore be assigned to varying experimental procedures or conditions, such as solvent composition, concentrations achieved during the measurements, or poorly defined substitution degree of CyD.
It is also revealed by the data in Table 2 and Table 4 that the complex stability constants for both guests show a tendency towards lower values for BCyD. In the case of fenbufen, among the methyl substituted CyD hosts, the stability increases with the substitution grade, and also with homogeneity. Among the HPBCyD hosts, the stability increases with the substitution degree. The trends are much less pronounced for fenoprofen.

4. Materials and Methods

4.1. Materials

All cyclodextrins used were products of Cyclolab Ltd. (Budapest, Hungary). Cyclolab Ltd. have provided us with quality certificates for the batches in question: BCyD—beta-cyclodextrin; CRYSMEB randomly methylated methyl-beta-cyclodextrin with low average substitution degree (/CYL4429/DS = 4.9; Purity > 99%); DIMEB—2,6-Di-O-methyl-beta-cyclodextrin as 3 differently derivatized substance: (1). DIMEB50 (/CYL4477/DS = 15.8, Purity > 95%, Isomer purity 35,7%); (2). DIMEB80 (/CYL2305/DS = 14, Purity > 95%, Isomer purity 80%); (3). DIMEB95 (/CYL4622/DS = 14.5, Purity > 99%, Isomer purity 93,4%); RAMEB—Random methyl-beta-cyclodextrin (/CYL4537/DS = 11.2, Purity > 99%); HPBCYD (2-Hydroxypropyl)-beta-cyclodextrin: (1). HPBCyD(4.5) (/CY-2005/.2.27 DS = 4.5); (2). HPBCyD(6.3) (/CY-2005/DS = 6.3).
Fenbufen (γ-Oxo-(1,1′-biphenyl)-4-butanoic acid) and fenoprofen calcium salt (calcium (±)-2-(3-phenoxyphenyl)propanoate dihydrate (Figure 2) as well as methanol-d4 (CD3OD), D2O and sodium-deuteroxide (NaOD) were purchased from Sigma-Aldrich/Merck Group (Budapest, Hungary).

4.2. Preparation of Solutions

To ensure that the conditions were identical, both the CD and NMR measurements were carried out on the same solutions. The solutions had to be prepared so that the concentration was high enough for the less sensitive NMR measurements, but not so high to produce absorbance values of higher than 2 during the CD measurements in order to avoid increasing spectral noise. As a first step, stock solutions of the two NSAIDs and of the CyDs were prepared.
Fenbufen stock solution of 10 mM was prepared in D2O by adding a small amount of NaOD solution to aid dissolution. Since fenoprofen calcium salt cannot be dissolved in pure water, methanol-d4 was used to prepare a 4 mM stock solution. Cyclodextrin stock solutions were prepared in H2O (22–25 mM for fenbufen and 12 mM for fenoprofen).
The difference in concentration between the solutions of the two guest molecules is due to the differing chromophores (Figure 7). Whereas in fenbufen, the conjugation of the biphenyl group and the carbonyl oxygen gives a broad, high intensity absorption band at 225–325 nm, in fenoprofen, the two aromatic rings are separated by an etheric oxygen. This results in a characteristic triple band and low intensity spectrum for the aromatic Lb (benzoid) band.
For the measurements, a series of solutions was prepared, each containing the same amount of fenbufen or fenoprofen, but different volumes of CyD stock solution. The final volume of each solution was 1.00–1.50 mL for fenbufen and 5.00 mL for fenoprofen. The ratio of guest molecule to CyD varied between 1:0.25 and 1:20 and all the solutions contained approximately 10% deuterated solvent (D2O in fenbufen solutions and methanol-d4 in fenoprofen solutions) which provided the lock signal in NMR measurements. As methanol does not have significant interaction with CyDs [41,42] small amount of methanol was also added to the fenbufen solution as an NMR chemical shift reference standard. In the measured samples, the concentration of fenbufen was about 1 mM and that of fenoprofen 0.4 mM. For all cases, the spectrum of an aqueous solution containing the same concentration of NSAID without CyD was also recorded. The exact concentrations of the CyDs and guest compounds in the solution series can be found in Supplementary Tables S1–S23.

4.3. Circular Dichroism (CD) Measurements

CD and UV experiments were performed on a Jasco J-815 spectrometer (Jasco LTD, Tokyo, Japan) in cylindrical Hellma cuvettes with pathlengths of 0.1 cm for fenbufen and 1.0 cm for fenoprofen. The slit was set to 1 nm, the registration speed was set to 50 nm/min, and the accumulation was set to 5. In order to obtain accurate readings of the measured ICD values, noise filtering was performed using the Fast Fourier Transform menu in Spectra Analysis program. For fenbufen, the ellipticity data were read at the peak maximum, for fenoprofen, at the 2 or 3 peak maxima (the number of maxima depends on the type of the CyD).

4.4. Nuclear Magnetic Resonance (NMR) Measurements

1H NMR spectra were recorded at 25 °C on a Varian (Varian, Palo Alto, California, United States) DDR spectrometer (599.9 MHz for 1H) equipped with a dual 5 mm inverse detection gradient (IDPFG) probe-head with z-gradient. 1H NMR spectra were collected in 32,000 data points using 32 scans with a spectral window of 6000 Hz and referenced to the methanol (fenbufen solutions) or the methanol-d4 (fenoprofen solutions) signals. For the suppression of the water signal, the double pulse field gradient spin echo (DPFGSE) sequence was used.
The 2D ROESY spectrum of the fenoprofen-BCyD 1:5 complex was recorded at 25 °C using Bruker (Bruker, Billerica, Massachusetts, United States) Avance III spectrometer (700.1 MHz for 1H) equipped with a 5 mm inverse TCI probehead with z-gradient. The Bruker rosygpph19 pulse sequence was used with a mixing time of 300 ms and a spectral width of 10 ppm in both dimensions. The spectrum was acquired into 1024 complex points in t2 with 128 scans coadded at each of 256 t1 increments. The spectrum was processed to 2048 × 1024 data points.

4.5. Mathematical Equations

For the regression analysis, the following relationships were taken into account. The equilibrium constant (stability constant, K) that describes the formation of the complex (C) from 1:1 stoichiometric association of the host (H) and guest (G) can be written as follows:
K = C H · G
The equilibrium concentrations [ ] in the above equation can be expressed in terms of total concentrations [ ]T as follows:
H = H T C
G = G T C
Substituting Equations (3) and (4) into Equation (2) and rearranging the equation one affords a quadratic solution for the complex concentration:
C = G T + H T + 1 K G T + H T + 1 K 2 4 · G T · H T 2
In CD measurements, the ellipticity depends only on the amount of complex. Since the concentration of guest is kept constant during the experiments, at excess host concentrations, the concentration of the complex is essentially equal to the total guest concentration; as a result, the maximal ellipticity response (θmax) can be written as:
θ m a x G T
For ellipticity values below the maximum Equation (7) holds:
θ = θ m a x · C G T
For the regression analysis of NMR data the following analogous Equation (8) was applied, that was obtained from Equation (1), where Δδ is the Complexation Induced Shift (i.e., δobsδG), and Δδmax is the maximal chemical shift response (i.e., δCδG).
δ = δ m a x · C G T
The non-linear regression analyses were performed using Origin Pro 8 (OriginLab Corp., Northampton, MA, USA). The standard errors of the fitted parameters were used to calculate the Gaussian propagation of uncertainty to the other parameters derived in the Results chapter. The Bland–Altman plot was constructed based on guidelines of the authors [43] in R version 4.1.3 (R Foundation for Statistical Computing, Vienna, Austria) [44].

5. Conclusions

The main purpose of this study was to investigate two methods used for characterizing cyclodextrin complexes: CD and NMR spectroscopy. The difference between the two methods in the determination of CyD-complex stability constants is not significant, as expected. The ICD data provide selective information on the complex formation and the concomitant binding perturbations. For ICD spectra with one maximum, it is advised to obtain the stability constant from the ellipticity data at the maximum. However, for structured ICD spectra, as is the case for fenoprofen, ellipticity data readings at multiple wavelength maxima are more appropriate followed by a global fit with mixed effects. Regression analysis of NMR data requires more precaution; the appropriate signals must be chosen based on the CIS magnitude. The substitution degree of CyD and their isomeric purity influences the stability constant. For both guests, native BCyD affords one of the lowest stability constants, while the substitution degree and isomeric purity increase the stability constant. For fenbufen, the order of stability constants is as follows:
BCyD ≈ CRYSMEB < HPBCyD4.5 < RAMEB ≈ HPBCyD6.3 ≈ DIMEB50 < DIMEB80 ≈ DIMEB95.
No similar trend was found for fenoprofen, although slightly higher constants were obtained for CyD derivatives with higher degrees of substitution. In the comparison between the guests (fenoprofen–fenbufen), it is clear that due to its linear shape, fenbufen forms more stable complexes compared to fenoprofen, since in the latter, the two aromatic rings are connected via an etheric oxygen, granting the molecular structure higher degree of freedom.
The results of the present study show that stability constants should only be compared if the average degree of substitution, the isomeric purity, and the measurement conditions are well defined. Our measured stability constant for the fenbufen-BCyD complex is in excellent agreement with a literature value measured by ICD by Sortino et al. [32] and in good agreement with constants measured by other methods [31,32]. For the fenoprofen-BCyD complex, our measured stability constant also shows good agreement with the literature value determined by 1H NMR [33].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24087544/s1.

Author Contributions

Conceptualization, P.H.; methodology, M.K., D.H. and P.H.; formal analysis, F.Á.; investigation, M.K., F.Á. and D.H.; resources, P.H.; data curation, M.K., A.M. and P.H.; writing—original draft preparation, A.M.; writing—review and editing, M.K., D.H., A.M. and P.H.; visualization, M.K., A.M. and P.H.; supervision, P.H.; project administration, P.H.; funding acquisition, A.M. and P.H. All authors have read and agreed to the published version of the manuscript.

Funding

The research was financed by the Higher Education Institutional Excellence Programme within the framework of the molecular biology thematic programme of Semmelweis University (TKP2021-EGA-24), and the ÚNKP-22-5-SE-3 New National Excellence Program of the Ministry for Culture and Innovation from the source of the National Research, Development and Innovation Fund. Arash Mirzahosseini is grateful for the János Bolyai Research Scholarship of the Hungarian Academy of Sciences. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article and supplementary material. The recorded spectra presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Szejtli, J. Introduction and general overview of cyclodextrin chemistry. Chem. Rev. 1998, 98, 1743–1754. [Google Scholar] [CrossRef] [PubMed]
  2. Del Valle, E.M.M. Cyclodextrins and their uses: A review. Process Biochem. 2004, 39, 1033–1046. [Google Scholar] [CrossRef]
  3. Hedges, A.R. Industrial applications of cyclodextrins. Chem. Rev. 1998, 98, 2035–2044. [Google Scholar] [CrossRef]
  4. Szejtli, J. Past, present and future of cyclodextrin research. Pure Appl. Chem. 2004, 76, 1825–1845. [Google Scholar] [CrossRef]
  5. Crini, G. Review: A history of cyclodextrins. Chem. Rev. 2014, 114, 10940–10975. [Google Scholar] [CrossRef] [PubMed]
  6. Dodziuk, H. Cyclodextines and Their Complexes: Chemistry, Analytical Methods, Applications; Wiley Verlag GmbH& Co. KGaA: Weinheim, Germany, 2016; ISBN 13:978-3-527-31280-1. [Google Scholar]
  7. Marugan, J.; McKew, C.J.; Austin, P.C.; Porter, D.F.; Silber, A.S.; Ory, S.D.; Walkley, U.S.; H.Vite, C.; Heiss, J.; Pavan, J.W.; et al. Collaborative development of 2-hydroxypropyl-β-cyclodextrin for the treatment of Niemann-Pick type C1 disease. Curr. Top. Med. Chem. 2014, 14, 330–339. [Google Scholar] [CrossRef]
  8. Miranda, G.M.; Santos, V.O.; Bessa, J.R.; Teles, Y.C.F.; Yahouédéhou, S.C.; Goncalves, M.S.; Ribeiro-Filho, J. Inclusion complexes of non-steroidal anti-inflammatory drugs with cyclodextrins: A systematic review. Biomolecules 2021, 11, 361. [Google Scholar] [CrossRef]
  9. Bui, C.V.; Rosenau, T.; Hettegger, H. Polysaccharide- and β-cyclodextrin-based chiral selectors for enantiomer resolution: Recent developments and applications. Molecules 2021, 26, 4322. [Google Scholar] [CrossRef]
  10. Choi, H.; Oh, Y.H.; Park, S.; Lee, S.S.; Oh, H.B.; Lee, S. Unveiling host–guest–solvent interactions in solution by identifying highly unstable host–guest configurations in thermal non-equilibrium gas phase. Sci. Rep. 2022, 12, 8169. [Google Scholar] [CrossRef] [PubMed]
  11. Chankvetadze, B.; Scriba, G.K.E. Cyclodextrins as chiral selectors in capillary electrophoresis: Recent trends in mechanistic studies. TrAC Trends Anal. Chem. 2023, 160, 116987. [Google Scholar] [CrossRef]
  12. Li, Z.; Couzijn, E.P.A.; Zhang, X. Intrinsic properties of α-cyclodextrin complexes with benzoate derivatives in the gas phase: An experimental and theoretical study. J. Phys. Chem. B 2012, 116, 943–950. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, S.-S.; Lee, J.-u.; Oh, J.H.; Park, S.; Hong, Y.; Min, B.K.; Lee, H.H.L.; Kim, H.I.; Kong, X.; Lee, S.; et al. Chiral differentiation of d- and l-isoleucine using permethylated β-cyclodextrin: Infrared multiple photon dissociation spectroscopy, ion-mobility mass spectrometry, and DFT calculations. Phys. Chem. Chem. Phys. 2018, 20, 30428–30436. [Google Scholar] [CrossRef] [PubMed]
  14. Szente, L.; Szemán, J. Cyclodextrins in analytical chemistry: Host-guest type molecular recognition. Anal. Chem. 2013, 85, 8024–8030. [Google Scholar] [CrossRef] [PubMed]
  15. Balzano, F.; Uccello-Barretta, G.; Sicoli, G.; Vanni, L.; Recchimurzo, A.; Aiello, F. Chiral Discrimination Mechanisms by Silylated-Acetylated Cyclodextrins: Superficial Interactions vs. Inclusion. Int. J. Mol. Sci. 2022, 23, 13169. [Google Scholar] [CrossRef] [PubMed]
  16. Poulson, B.G.; Alsulami, Q.A.; Sharfalddin, A.; El Agammy, E.F.; Mouffouk, F.; Emwas, A.-H.; Jaremko, L.; Jaremko, M. Cyclodextrins: Structural, chemical, and physical properties, and applications. Polysaccharides 2022, 3, 1–31. [Google Scholar] [CrossRef]
  17. Guo, C.; Xiao, Y. Negatively charged cyclodextrins: Synthesis and applications in chiral analysis-A review. Carbohydr. Polym. 2021, 256, 117517. [Google Scholar] [CrossRef] [PubMed]
  18. Kalydi, E.; Malanga, M.; Ujj, D.; Benkovics, G.; Szakács, Z.; Béni, S. Fully symmetric cyclodextrin polycarboxylates: How to determine reliable protonation constants from NMR titration data. Int. J. Mol. Sci. 2022, 23, 14448. [Google Scholar] [CrossRef] [PubMed]
  19. Liu, J.Y.; Zhang, X.; Tian, B.R. Selective modifications at the different positions of cyclodextrins: A review of strategies. Turk. J. Chem. 2020, 44, 261–278. [Google Scholar] [CrossRef]
  20. Řezanka, M. Synthesis of substituted cyclodextrins. Environ. Chem. Lett. 2019, 17, 49–63. [Google Scholar] [CrossRef]
  21. Szente, L.; Szemán, J.; Sohajda, T. Analytical characterization of cyclodextrins: History, official methods and recommended new techniques. J. Pharm. Biomed. Anal. 2016, 130, 347–365. [Google Scholar] [CrossRef]
  22. Kiss, E.; Szabó, V.A.; Horváth, P. Simple circular dichroism method for selection of the optimal cyclodextrin for drug complexation. J. Incl. Phenom. Macrocycl. Chem. 2019, 95, 223–233. [Google Scholar] [CrossRef]
  23. Castro-Puyana, M.; Crego, A.L.; Marina, M.L.; García-Ruiz, C. Enantioselective separation of azole compounds by EKC. Reversal of migration order of enantiomers with CD concentration. Electrophoresis 2007, 28, 2667–2674. [Google Scholar] [CrossRef]
  24. Van Doorne, H.; Bosch, E.H.; Lerk, C.F. Formation arid antimicrobial activity of complexes of β-cyclodextrin and some antimycotic imidazole derivatives. Pharm. Weekbl. 1988, 10, 80–85. [Google Scholar] [CrossRef] [PubMed]
  25. Penn, S.G.; Bergstrom, E.T.; Goodall, D.M.; Loran, J.S. Capillary electrophoresis with chiral selectors: Optimization of separation and determination of thermodynamic parameters for binding of tioconazole enantiomers to cyclodextrins. Anal. Chem. 1994, 66, 2866–2873. [Google Scholar] [CrossRef]
  26. Kelemen, H.; Csillag, A.; Hancu, G.; Székely-Szentmiklósi, B.; Fülöp, I.; Varga, E.; Grama, L.; Orgován, G. Characterisation of inclusion complexes between bifonazole and different cyclodextrins in solid and solution state. Maced. J. Chem. Chem. Eng. 2017, 36, 81–91. [Google Scholar] [CrossRef]
  27. Mohammed, N.N.; Pandey, P.; Khan, N.S.; Elokely, K.M.; Liu, H.; Doerksen, R.J.; Repka, M.A. Clotrimazole–cyclodextrin based approach for the management and treatment of candidiasis—A formulation and chemistry-based evaluation. Pharm. Dev. Technol. 2016, 21, 619–629. [Google Scholar] [CrossRef]
  28. Tenjarla, S.; Puranajoti, P.; Kasina, R.; Mandal, T. Preparation, characterization, and evaluation of miconazole-cyclodextrin complexes for improved oral and topical delivery. J. Pharm. Sci. 1998, 87, 425–429. [Google Scholar] [CrossRef] [PubMed]
  29. Loftsson, T.; Brewster, M.E. Cyclodextrins as functional excipients: Methods to enhance complexation efficiency. J. Pharm. Sci. 2012, 101, 3019–3032. [Google Scholar] [CrossRef] [PubMed]
  30. Bratu, I.; Gavira-Vallejo, J.M.; Hernanz, A.; Bogdan, M.; Bora, G. Inclusion complex of fenbufen with β-cyclodextrin. Biopolymers 2004, 73, 451–456. [Google Scholar] [CrossRef] [PubMed]
  31. Zhang, H.; Zhang, H.; Jin, W.; Ding, L. Determination of dissociation constants of cyclodextrin—Ligand inclusion complexes by electrospray ionization mass spectrometry. Eur. J. Mass Spectrom. 2006, 12, 291–299. [Google Scholar] [CrossRef] [PubMed]
  32. Sortino, S.; Giuffrida, S.; Fazio, S.; Monti, S. Spectroscopic characterization and photochemical behavior of host–guest complexes between β-cyclodextrin and drugs containing a biphenyl-like chromophore. New J. Chem. 2001, 25, 707–713. [Google Scholar] [CrossRef]
  33. Uccello-Barretta, G.; Chiavacci, C.; Bertucci, C.; Salvadori, P. Stereochemistry and dynamics of the inclusion complex of (S)-(+)-fenoprofen with cyclomaltoheptaose (beta-cyclodextrin). Carbohydr. Res. 1993, 243, 1–10. [Google Scholar] [CrossRef]
  34. Kavetsou, E.; Tsoukalas-Koulas, C.; Katopodi, A.; Kalospyros, A.; Alexandratou, E.; Detsi, A. Ιnclusion complexes of magnesium phthalocyanine with cyclodextrins as potential photosensitizing agents. Bioengineering 2023, 10, 244. [Google Scholar] [CrossRef] [PubMed]
  35. Ueda, K.; Higashi, K.; Moribe, K. Amorphous drug solubility and maximum free drug concentrations in cyclodextrin solutions: A quantitative study using NMR diffusometry. Mol. Pharm. 2021, 18, 2764–2776. [Google Scholar] [CrossRef] [PubMed]
  36. Udo, H.; Brinker, D.K. Circular dichroism of cyclodextrin complexes. In Cyclodextrins and Their Complexes; WILEY-VCH Verlag GmbH & Co.: Weinheim, Germany, 2006. [Google Scholar]
  37. Gawroński, J.; Grajewski, J. The significance of induced circular dichroism. Org. Lett. 2003, 5, 3301–3303. [Google Scholar] [CrossRef]
  38. Harata, K.; Uedaira, H. The Circular Dichroism spectra of the β-cyclodextrin complex with naphthalene derivatives. Bull. Chem. Soc. Jpn. 1975, 48, 375–378. [Google Scholar] [CrossRef]
  39. Kajtár, M.; Horváth-Toró, C.; Kuthi, É.; Szejtli, J. A simple rule for predicting circular dichroism induced in aromatic guests by cyclodextrin hosts in inclusion complexes. In Proceedings of the First International Symposium on Cyclodextrins; Szejtli, J., Ed.; Springer Netherlands: Dordrecht, The Netherlands, 1982; pp. 181–193. [Google Scholar]
  40. Kodaka, M. Application of a general rule to induced circular dichroism of naphthalene derivatives complexed with cyclodextrins. J. Phys. Chem. A 1998, 102, 8101–8103. [Google Scholar] [CrossRef]
  41. Sohajda, T.; Hu, W.H.; Zeng, L.L.; Li, H.; Szente, L.; Noszál, B.; Béni, S. Evaluation of the interaction between sitagliptin and cyclodextrin derivatives by capillary electrophoresis and nuclear magnetic resonance spectroscopy. Electrophoresis 2011, 32, 2648–2654. [Google Scholar] [CrossRef] [PubMed]
  42. Tűz, B.; Noszál, B.; Hosztafi, S.; Mazák, K. β-cyclodextrin complex formation and protonation equilibria of morphine and other opioid compounds of therapeutic interest. Eur. J. Pharm. Sci. 2022, 171, 106120. [Google Scholar] [CrossRef] [PubMed]
  43. Bland, J.M.; Altman, D.G. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet 1986, 327, 307–310. [Google Scholar] [CrossRef]
  44. R Core Team. R: A Language and Environment for Statistical Computing. 2022. Available online: https://www.R-project.org/ (accessed on 22 April 2022).
Figure 1. The structural formulae of the studied guest compounds (fenbufen and fenoprofen racemate) and the BCyD host compound, together with the numbering of H atoms of interest for 1H NMR.
Figure 1. The structural formulae of the studied guest compounds (fenbufen and fenoprofen racemate) and the BCyD host compound, together with the numbering of H atoms of interest for 1H NMR.
Ijms 24 07544 g001
Figure 2. The 1H NMR assignation of (a) fenbufen, (b) fenoprofen, and (c) BCyD.
Figure 2. The 1H NMR assignation of (a) fenbufen, (b) fenoprofen, and (c) BCyD.
Ijms 24 07544 g002
Figure 3. (a) Induced CD spectra of HPBCyD(4.5)-fenbufen complexes with increasing host concentration (highest host concentration corresponds to the tallest peak) in 90% H2O/10% D2O. (b) The regression model of the ellipticity data obtained from the CD measurements at λmax. (c) 1H NMR spectra series (aromatic region) of HPBCyD(4.5)-fenbufen complexes at different host–guest ratios in 90% H2O/10% D2O. (d) The regression model of the H2 NMR signal (obtained from global fit).
Figure 3. (a) Induced CD spectra of HPBCyD(4.5)-fenbufen complexes with increasing host concentration (highest host concentration corresponds to the tallest peak) in 90% H2O/10% D2O. (b) The regression model of the ellipticity data obtained from the CD measurements at λmax. (c) 1H NMR spectra series (aromatic region) of HPBCyD(4.5)-fenbufen complexes at different host–guest ratios in 90% H2O/10% D2O. (d) The regression model of the H2 NMR signal (obtained from global fit).
Ijms 24 07544 g003
Figure 4. (a) The induced CD spectra of CyD-fenoprofen 1:10 in 90% H2O/10% methanol-d4. (b) The regression model of the ellipticity data for BCyD-fenoprofen obtained from the CD measurements at λmax. (c) The 1H NMR spectra of CyD-fenoprofen 1:10 complexes in 90% H2O/10% methanol-d4. (d) The regression model of the H2 NMR signal of BCyD-fenoprofen (obtained from global fit).
Figure 4. (a) The induced CD spectra of CyD-fenoprofen 1:10 in 90% H2O/10% methanol-d4. (b) The regression model of the ellipticity data for BCyD-fenoprofen obtained from the CD measurements at λmax. (c) The 1H NMR spectra of CyD-fenoprofen 1:10 complexes in 90% H2O/10% methanol-d4. (d) The regression model of the H2 NMR signal of BCyD-fenoprofen (obtained from global fit).
Ijms 24 07544 g004
Figure 5. The ROESY NMR spectrum of fenoprofen-BCyD 1:5 complex in 90% H2O/10% methanol-d4. Assignation of the signals in blue is for fenoprofen, in black is for BCyD.
Figure 5. The ROESY NMR spectrum of fenoprofen-BCyD 1:5 complex in 90% H2O/10% methanol-d4. Assignation of the signals in blue is for fenoprofen, in black is for BCyD.
Ijms 24 07544 g005
Figure 6. Left: The deviance of each logK obtained from single NMR signals from its concomitant globally fit logK as a function of absolute maximal chemical shift response. Right: The Bland–Altman plot of the complex stability constants for CyD-fenbufen and CyD-fenoprofen complexes obtained with two methods: CD and 1H NMR spectroscopy.
Figure 6. Left: The deviance of each logK obtained from single NMR signals from its concomitant globally fit logK as a function of absolute maximal chemical shift response. Right: The Bland–Altman plot of the complex stability constants for CyD-fenbufen and CyD-fenoprofen complexes obtained with two methods: CD and 1H NMR spectroscopy.
Ijms 24 07544 g006
Figure 7. Specific absorptions of fenbufen (black) and fenoprofen calcium (blue) in methanolic solutions with extrema.
Figure 7. Specific absorptions of fenbufen (black) and fenoprofen calcium (blue) in methanolic solutions with extrema.
Ijms 24 07544 g007
Table 1. Comparison of stability constants. Abbreviations: CE—capillary electrophoresis; PS—phase solubility measurement, NMR—1H nuclear magnetic resonance spectroscopy.
Table 1. Comparison of stability constants. Abbreviations: CE—capillary electrophoresis; PS—phase solubility measurement, NMR—1H nuclear magnetic resonance spectroscopy.
CyD/Guest CompoundlogK Determined by CD [14]logK from Literature
BCyD/
Bifonazole3.403.40–3.52 (CE) [23]
Clotrimazole2.652.66 (PS) [24]
DIMEB/
Tioconazole3.863.84 (CE) [25]
HPBCyD/
Bifonazole4.463.66 (NMR) [26]
Clotrimazole2.092.65 (PS) [27]
Miconazole2.712.56 (PS) [28], 2.41 (PS) [29]
Tioconazole3.302.86 (CE) [25]
SBEBCyD/
Bifonazole4.723.94 (NMR) [26]
Table 2. The complex stability constants (in log10 units ± standard error) of fenbufen with the various cyclodextrins determined by CD and NMR spectroscopy in 90% H2O/10% D2O.
Table 2. The complex stability constants (in log10 units ± standard error) of fenbufen with the various cyclodextrins determined by CD and NMR spectroscopy in 90% H2O/10% D2O.
CyclodextrinlogKCDlogKNMRDiff.
BCyD (H2O) 13.35 ± 0.03NANA
BCyD3.34 ± 0.023.24 ± 0.060.10
CRYSMEB3.39 ± 0.083.33 ± 0.050.06
DIMEB503.56 ± 0.123.62 ± 0.07−0.06
DIMEB804.12 ± 0.284.00 ± 0.140.12
DIMEB954.06 ± 0.073.95 ± 0.100.10
RAMEB3.59 ± 0.083.47 ± 0.020.12
HPBCyD(4.5)3.39 ± 0.043.53 ± 0.06−0.14
HPBCyD(6.3)3.53 ± 0.043.61 ± 0.06−0.08
1 Initial trial experiment performed in H2O without any deuterated solvent.
Table 3. The complex stability constants of CyD-fenbufen complexes (every odd line) followed by the maximal chemical shift response (every even line) obtained from regression fits performed on the data of each individual NMR peak for CyD-fenbufen complexes. The uncertainty of chemical shift response values (on ppm scale) is on the order of 0.001.
Table 3. The complex stability constants of CyD-fenbufen complexes (every odd line) followed by the maximal chemical shift response (every even line) obtained from regression fits performed on the data of each individual NMR peak for CyD-fenbufen complexes. The uncertainty of chemical shift response values (on ppm scale) is on the order of 0.001.
CyclodextrinH1,1′H2,2′H3,3′H4,4′H5
BCyDNA3.413.302.952.70
−0.179−0.273−0.128−0.095
CRYSMEBNA3.313.373.193.23
−0.185−0.1850.0290.070
DIMEB50NA3.833.944.283.42
−0.201−0.1890.0340.109
DIMEB80NA4.094.063.813.55
−0.207−0.2000.0260.100
DIMEB95NA4.084.174.783.58
−0.220−0.2080.0140.090
RAMEBNA3.683.703.092.97
−0.176−0.1830.0410.095
HPBCyD(4.5)2.933.623.59NA3.09
0.036−0.134−0.1680.065
HPBCyD(6.3)3.103.713.67NA3.19
0.044−0.132−0.1670.069
Table 4. The complex stability constants (in log10 units ± standard error) of fenoprofen with the various cyclodextrins determined by CD and 1H NMR spectroscopy in 90% H2O/10% methanol-d4.
Table 4. The complex stability constants (in log10 units ± standard error) of fenoprofen with the various cyclodextrins determined by CD and 1H NMR spectroscopy in 90% H2O/10% methanol-d4.
CyclodextrinlogKCDlogKNMRDiff.
BCyD (H2O) 13.12 ± 0.02NANA
BCyD3.06 ± 0.052.98 ± 0.040.08
CRYSMEB3.10 ± 0.023.05 ± 0.010.05
DIMEB503.10 ± 0.083.13 ± 0.02−0.03
DIMEB803.33 ± 0.093.36 ± 0.01−0.03
DIMEB953.32 ± 0.093.26 ± 0.010.06
RAMEB3.34 ± 0.073.15 ± 0.0040.19
HPBCyD(4.5)3.18 ± 0.083.09 ± 0.010.09
HPBCyD(6.3)3.19 ± 0.053.14 ± 0.0040.05
1 Initial trial experiment performed in H2O without any deuterated solvent.
Table 5. The complex stability constants of CyD-fenoprofen complexes (every odd line) followed by the maximal chemical shift response (every even line) obtained from regression fits performed on the data of each individual NMR peak for CyD-fenoprofen complexes. The uncertainty of chemical shift response values (on ppm scale) is on the order of 0.001.
Table 5. The complex stability constants of CyD-fenoprofen complexes (every odd line) followed by the maximal chemical shift response (every even line) obtained from regression fits performed on the data of each individual NMR peak for CyD-fenoprofen complexes. The uncertainty of chemical shift response values (on ppm scale) is on the order of 0.001.
CyclodextrinH1H2H3H4H5,5′H6,6′H7
BCyD3.172.992.992.982.993.022.96
0.0260.068−0.040−0.390−0.179−0.0080.043
CRYSMEB3.103.093.023.053.053.383.11
−0.0180.033−0.090−0.415−0.1400.0100.051
DIMEB503.163.113.163.123.163.203.16
−0.0760.088−0.126−0.556−0.1230.0320.094
DIMEB803.463.253.433.353.413.603.45
−0.0760.071−0.112−0.519−0.1390.0180.082
DIMEB953.393.193.313.293.393.533.40
−0.0790.074−0.118−0.532−0.1380.0190.083
RAMEB3.183.133.143.153.133.223.13
−0.0400.064−0.136−0.510−0.0720.0680.109
HPBCyD(4.5)NA3.073.123.093.113.003.04
NA0.059−0.071−0.395−0.1180.0270.063
HPBCyD(6.3)2.943.153.133.143.123.343.16
0.0160.069−0.066−0.395−0.1180.0310.067
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

Kraszni, M.; Ágh, F.; Horváth, D.; Mirzahosseini, A.; Horváth, P. Effect of Substitution Degree and Homogeneity on Cyclodextrin-Ligand Complex Stability: Comparison of Fenbufen and Fenoprofen Using CD and NMR Spectroscopy. Int. J. Mol. Sci. 2023, 24, 7544. https://doi.org/10.3390/ijms24087544

AMA Style

Kraszni M, Ágh F, Horváth D, Mirzahosseini A, Horváth P. Effect of Substitution Degree and Homogeneity on Cyclodextrin-Ligand Complex Stability: Comparison of Fenbufen and Fenoprofen Using CD and NMR Spectroscopy. International Journal of Molecular Sciences. 2023; 24(8):7544. https://doi.org/10.3390/ijms24087544

Chicago/Turabian Style

Kraszni, Márta, Ferenc Ágh, Dániel Horváth, Arash Mirzahosseini, and Péter Horváth. 2023. "Effect of Substitution Degree and Homogeneity on Cyclodextrin-Ligand Complex Stability: Comparison of Fenbufen and Fenoprofen Using CD and NMR Spectroscopy" International Journal of Molecular Sciences 24, no. 8: 7544. https://doi.org/10.3390/ijms24087544

APA Style

Kraszni, M., Ágh, F., Horváth, D., Mirzahosseini, A., & Horváth, P. (2023). Effect of Substitution Degree and Homogeneity on Cyclodextrin-Ligand Complex Stability: Comparison of Fenbufen and Fenoprofen Using CD and NMR Spectroscopy. International Journal of Molecular Sciences, 24(8), 7544. https://doi.org/10.3390/ijms24087544

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