Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 2680 (2023 ) Cite this article Plastic Column Formwork
Cyclodextrins (CDs) have been previously shown to display modest equilibrium binding affinities (Ka ~ 100–200 M-1) for the synthetic opioid analgesic fentanyl. In this work, we describe the synthesis of new CDs possessing extended thioalkylcarboxyl or thioalkylhydroxyl moieties and assess their binding affinity towards fentanyl hydrochloride. The optimal CD studied displays a remarkable affinity for the opioid of Ka = 66,500 M−1, the largest value reported for such an inclusion complex to date. One dimensional 1H Nuclear Magnetic Resonance (NMR) as well as Rotational Frame Overhauser Spectroscopy (2D-ROESY) experiments supported by molecular dynamics (MD) simulations suggest an unexpected binding behavior, with fentanyl able to bind the CD interior in one of two distinct orientations. Binding energies derived from the MD simulations work correlate strongly with NMR-derived affinities highlighting its utility as a predictive tool for CD candidate optimization. The performance of these host molecules portends their utility as platforms for medical countermeasures for opioid exposure, as biosensors, and in other forensic science applications.
Fentanyl-based synthetic opioids have found numerous applications in the medical field as effective anesthetics and are considered safer than morphine during perioperative procedures1,2,3. Their potent activity stems from their ability to bind various nociceptive receptors in the nervous system, with one of their main targets being the μ-opioid receptor4,5. Like morphine and other alkaloids designed to combat pain, fentanyls unfortunately display high risks of dependency and addiction due to their euphoric effect on the central nervous system (CNS). Unfortunately, their illicit use has led to a high rate of overdose-related death in what is largely recognized as a fentanyl epidemic6,7,8,9.
The flagship compound in this opioid class is fentanyl, which has a potency of ~ 100 times that of morphine (Fig. 1a)10,11. Though both fentanyl and morphine display toxic potential if not carefully administered to an individual, fentanyl is often of higher concern due to its much simpler production. Furthermore, given the opioid’s ease of synthesis, the production of novel analogs characterized by increased analgesic profiles is not difficult12,13,14. In addition to their medical uses, some of these opioids such as acetylfentanyl have been linked to fatal overdose cases from recreational use, while carfentanil and remifentanil have been linked to use as chemical warfare agents (Fig. 1b)15,16. As fentanyl-based opioids have increased in both potency and use over the last several decades17, significant research has focused on the development of efficient analytical methods for their detection9,16,18 and on effective antidotes to counteract their physiological effects19,20. Antidotes such as naloxone and naltrexone (Fig. 1c) have become the principal methods of overdose treatment, however, they still lack the ability to provide prolonged protection against these synthetic opioids. Such reduced effectiveness is due to the opioid’s long circulation half-life in the blood (t1/2 ~ 7–9 h) whereas the circulation half-lives for naloxone and naltrexone are significantly shorter at t1/2 ~ 2 h and t1/2 ~ 4 h, respectively19. To address this issue, research groups have launched programs aimed at developing medical countermeasure (MCMs) candidates that provide longer-term protection against fentanyls, particularly when the affected individual is removed from emergency medical intervention. Some of these efforts have revolved around the notion of employing host molecules that could in principle capture synthetic opioids in their interior and depending on their binding affinity towards it, these could be effective scaffolds to be used as MCMs. One class of host systems that has been studied in detail by the Isaacs group at the University of Maryland for this purpose and also for the reversal of anesthesia has been the cucurbit[n]urils21,22. Their work that started with these cyclic-based hosts led to the eventual discovery of a powerful host for methamphetamine23 and fentanyl24 was curiously acyclic in nature, therefore providing it with enough flexibility to modify its effective trapping of the opioid. Another set of similar hosts but based on a carbohydrate backbone are the cyclodextrins, which unfortunately have been marginally evaluated for these applications25.
(a) Chemical structures of fentanyl and morphine, the potency of the opioid relative to morphine (~ 100x) is provided along with its systemic circulation half-life (t1/2); (b) structures of two of the most common antidotes used to treat fentanyl poisoning, naloxone and naltrexone along with their systemic circulation half-lives; (c) structures of fentanyl analogs encountered in several overdose cases in the US (acetylfentanyl) as well as others used as incapacitating agents (carfentanil and remifentanil), their potencies relative to morphine are provided; (d) structure of β-cyclodextrin showing the seven glucose units linked in an α1,4-fashion giving rise to a conical structure (frustrum) open at both ends with a hydrophobic interior and a hydrophilic exterior; (e) representation of a hypothetical inclusion complex formed between β-CD and fentanyl.
Cyclodextrins (CDs), cyclic oligosaccharide motifs composed of glucosyl units linked together via α1,4-glycosidic bonds (Fig. 1d), have recently garnered attention as potential medical countermeasure candidates. The cyclic arrangement of these glucose units gives rise to a macromolecular scaffold that resembles a frustrum with a bottom rim larger than its top (Fig. 1e). This arrangement provides CDs with a hydrophilic exterior and a hydrophobic interior able to accommodate organic molecules, a key physical property of these compounds26,27. This ability to serve as a molecular host has made CDs versatile platforms for use in industry28, medicine29, and materials science30. Naturally, the larger the number of glucose units, the larger the frustrum’s cavity is and thus its ability to host larger guests in its interior. A crucial and useful experimental handle researchers have is the host:guest complexes, also called inclusion complexes, can be studied by various spectroscopic methods that include UV-Vis31, infra-red32, and nuclear magnetic resonance spectroscopies33 as well as spectrometric methods like liquid chromatography mass spectrometry (LC–MS)34. Our previous work has relied principally on NMR to study complexes between native CDs and fentanyl and other toxins (Fig. 1e)25,35. Making use of extensive NMR titration studies in conjunction with two-dimensional Rotational Frame Overhauser Spectroscopy (2D-ROESY), one can determine CD:fentanyl binding affinities as well as the structure of the complex. Even though it was previously demonstrated inclusion complexes form between native CDs and fentanyl, the binding affinities were low and not significantly strong to make a case for their use as a potential fentanyl binder. In this report, we present the synthesis of more complex CD scaffolds observed to form strong inclusion complexes with fentanyl. Our efforts have leveraged NMR titrations to determine the binding strength between the current lead candidate, called “subetadex” (SBX), and fentanyl. In addition, important insight into the structure of the SBX:fentanyl complex was accomplished using 2D-ROESY NMR experiments. These experimental data were found to be in good agreement with computational modeling calculations that were also used to explore additional CD scaffold candidates with enhanced fentanyl affinities.
This work focuses on the exploration of more complex CD scaffolds encouraged by our initial studies on the capture of fentanyl with native and simply substituted CDs25 as well as the knowledge of other CD scaffolds that have successfully found major biomedical applications for their ability to bind and effectively neutralize specific targets. One such CD is Sugammadex (SGX), which is a FDA-approved drug for the treatment of patients who experience prolonged post-operative anesthetic effects of rocuronium bromide (RocBr)36,37,38,39,40. SGX was discovered in 2001 at Organon Laboratories36 and displayed a strong affinity for RocBr while minimizing potentially toxic side effects. After years of evaluation, SGX was approved by the FDA in 2015 and has now become a common treatment to alleviate the anesthetic symptoms associated with the use of RocBr in the surgical setting due to its non-toxic nature and the strong complex association constant for the anesthetic (Ka ~ 750,000 M-1) determined by nuclear magnetic resonance (NMR) titration experiments41. Like SGX, the present work sought to identify a CD candidate with a non-toxic profile while displaying a high affinity for fentanyl and related synthetic opioids.
Structurally, SGX is based on a γ-CD core having its upper rim chemically modified with 3-mercaptopropanoic acid chains. These additional moieties serve to elongate the native CD hydrophobic interior. Based on our previous studies25, the γ-CD cavity was found to be too large to accommodate fentanyl in its interior, and MD calculations showed the opioid spent the majority of the simulation time unbound. For this reason, initial screening efforts were focused on β-cyclodextrin-(β-CD)-based scaffolds, as their interior better accommodated the opioid. Our previous NMR data and MD simulations both showed fentanyls bind native β-CD at the fentanyl amide phenyl ring, though observed binding affinities were quite low (Ka < 250 M-1). The structures of these complexes determined from these data were found to leave the majority of the fentanyl molecule outside of the CD, with its phenethyl group largely solvated by water. These enthalpically unfavorable interactions in part explain in the observed weak complexes. To improve complexation, CDs were envisioned that enclose a far greater fraction of the opioid thereby providing more favorable interactions, principally van der Waals and Coulombic34. The library screened in this current work is based on the β-CD host with slight variations on a common modification, namely the replacement of all 6-position hydroxyl groups with straight-chain thioalkylcarboxylates or thioalkylalcohols (Fig. 2).
Basic structural scheme of a βCD with its upper rim, C6-hydroxyl moiety replaced with a generalized R group. There are two types of moieties investigated in this work, the anionic ones comprising CDs featuring the 2-mercaptoacetyl (SBX-1, 1), 3-mercaptopropionoyl (SBX, 2) and 4-mercaptobutanoyl (SBX + 1, 3) groups. The second type of moieties involve ones bearing neutral moieties in the upper rim featuring the 2-mercaptoethanol (SBN-1, 4), 3-mercaptopropanol (SBN, 5), and 4-mercaptobutanol (SBN + 1, 6). In addition to these βCD-based targets, SAX (Sualphadex, 7) and SGX (Sugammadex, 8) were also synthesized and evaluated separately, however, both did not show any significant binding to fentanyl.
Equilibrium binding affinities between the CD panel and fentanyl·HCl were determined using 1H NMR titration experiments. In order to derive the most accurate binding constant, experiments were designed such that the CD concentration [H]T (in units of M) satisfied the inequality [H]T K < 1 whenever possible. Initial studies required best-guess estimates of K since they were a priori unknown. The 1H NMR chemical shifts of reporter protons on CD (H3 and H5) were monitored while the fentanyl titrant solution was added incrementally. Representative data from these titrations are found in Fig. 2 of Supporting Information. Non-linear least squares fitting was performed using mathematical expressions appropriate for the binding data obtained. This fitting was performed on all reporter proton titration curves simultaneously yielding a global binding affinity for each CD:fentanyl pair25. Results from this fitting procedure are given in the first data column of Table 1.
In most cases, the NMR data could be fit to a model that describes the binding as a simple 1:1 system. NMR Job plot data, however, which is used to justify the use of such a model showed distinct deviation from a 1:1 binding model for both SBX and SBX + 1, with maxima slightly off from the expected fraction, r = 0.5 (see Fig. 1 in Supporting Information). These data urged the consideration of a two-state 1:1 model which was able to describe both the observed Job plot behavior and that from the titrations themselves. Additionally, the titration data showed distinctly non-monotonic behavior for multiple reporter protons (see Fig. 2 in the Supporting Information).
Figure 3 compares NMR-determined binding affinities for 1:1 CD:fentanyl complexes to MD simulation binding energies (values given in right-most columns of Table 1). Data for α-CD and β-CD from our previous work are included for completeness. For the anionic SBX-1, SBX and SBX + 1 (2–4), NMR results show strong correlation with behavior predicted by the MD simulations. That is, binding affinity increases with binding energy likely due to enhanced van der Waals interactions afforded by the hydrophobic nature of the elongated backbone of the anionic thioalkylcarboxylates. This assertion is supported by additional computational work showing an plateau in affinity for CDs displaying further-extended anionic arms (e.g., SBX + 2 and SBX + 3 (not synthesized)). This observation indicates there is little advantage to thiocarboxylate arms longer than four carbons (i.e., SBX + 1), which are sufficient to encase completely the fentanyl molecule (i.e., optimal size complementarity). We refrain, however, from overinterpreting any direct NMR/MD comparisons, as MD-derived quantities neglect the effect of conformational entropy in the overall thermodynamics. In other words, a free energy is not strictly calculated. Instead, an average energy is calculated, referred to as < U + W > as previously described41.
NMR binding affinities (simple 1:1 model) versus enthalpic energies determined from simulation for the dominant conformer of fentanyl bound to hosts 1–6 and to α- and β-cyclodextrin31. Note the affinity reported here for SBX + 1 is the dominant conformation from the more complex two-state 1:1 model (67 mM-1) Neutral CD scaffolds given as circles; all others as filled squares. The solid line is a linear fit to filled data points considering measurement uncertainties. The dashed line has slope RT (0.602 kcal mol-1 at 303 K).
For the neutral CDs (i.e., SBN-1, SBN and SBN + 1) initial MD simulations predicted an enhanced binding affinity with increasing arm length in similar fashion to their anionic counterparts. NMR data, however, show a dramatic decrease in binding affinity, one that actually worsens with increasing arm length (see Table 1, Fig. 1). NMR structural data from ROESY experiments show clear correlations between the CD’s interior H3 and H5 protons and those along its thioalkylhydroxyl arms (see Supporting Information). Additional MD simulations of these neutral CDs alone in solution supported this behavior. Ultimately, the hydrophobic driving force that exists for the hydrophobic thioalkylhydroxyl arms to orient themselves inside the CD cavity seems to preclude fentanyl from binding via a competitive, self-binding effect with longer arms binding better within the cavity. We then infer the thiocarboxylate arms of the charged CDs not only provide the necessary aqueous solubility of the CD but, arguably more importantly, likely hold the anionic SBX open wide enough for fentanyl to bind via electrostatic repulsion. Note MD suggests there is little interaction between the carboxylate terminuses and the fentanyl molecule, as the charged groups orient into their aqueous surrounding.
The results presented thus far represent a simplified view of structural details of the SBX:fentanyl complexes. As mentioned earlier, NMR data in the form of Job plots (see Fig. 1 in Supporting Information), which are used to determine host:guest stoichiometries42,43,44, revealed more complex binding behavior initially thought to originate from a 1:2 or 2:1 host:guest complex. MD simulations, however, showed CD could not accommodate a second fentanyl molecule. Additional simulations, however, showed the charged CDs could accommodate fentanyl in two distinct conformations (Fig. 4). In the case of SBX-1, the orientation with the amide half of fentanyl pointed “down” toward the anionic chain ends is slightly energetically favored over fentanyl in the opposite direction (see Fig. 4a). For SBX + 1, the reversed configuration is favored with the amide half pointed “up” toward the unmodified secondary rim of CD (see Fig. 4b). Both of these conformations were observed to be favorable, with relatively small differences in energies between them (0.8–5.8 kcal mol-1 depending on host molecule, Fig. 4c). For SBX, the two conformations are equally favorable.
MD results for the host:guest complex with two orientations: (a) fentanyl (carbon atoms in cyan) aligned “down” with amide half near primary rim of SBX (carbon atoms in green), (b) fentanyl in alternate, “up” position within SBX, c) binding energies for all three SBX complexes.
MD findings are consistent with rotating-frame nuclear Overhauser effect spectroscopy (ROESY) NMR experiments. Figure 5 gives the two-dimensional ROESY spectrum of a SBX + 1:fentanyl mixture (see also Fig. 3 in Supporting Information). Of particular importance is the complete absence of correlations between methylene protons on the anionic CD arm and aromatic protons on the amide nitrogen-bound phenyl ring. Rather, only correlations between these methylenes and the phenethyl aromatic protons of fentanyl were found. This finding indicates the majority of fentanyl binds in the “up” orientation, consistent with MD results. ROESY data for a SBX-1:fentanyl mixture shows a similar agreement to MD data, namely, amide phenyl protons correlate to interior CD protons near its unmodified, secondary rim. These peaks provide clear evidence for a preference of the complex to adopt the “down” orientation. For all three charged CDs, Job plot data (see select data in Fig. 1 of Supporting Information) display evidence of a two-state binding system regardless of a lack of clear correlations in the ROESY spectra. Job plots are much more sensitive at detecting these types of phenomena, as they cover a substantially wider relative concentration range (i.e. [H]0/[G]0 ratio) than titration experiments.
Partial ROESY spectrum of SBX + 1:fentanyl HCl complex. Correlations between fentanyl aromatic protons from the phenethyl ring and SBX + 1 methylenes are highlighted by the dashed box. Full spectrum is given in Fig. 3 of Supporting Information.
Evidence of an equilibrium coexistence of the two conformations implies the 1:1 binding model invoked previous is not strictly appropriate for the extraction of equilibrium binding affinities. We therefore developed a fitting model that explicitly considers both orientations on the chemical shift titration data. Details of this “self-competitive” model are given in Supporting Information. The results of fitting to such a model yielded two equilibrium binding constants. For SBX:fentanyl, the larger of the two Ka values was observed to be statistically similar to that derived from the 1:1 model. Secondary binding constants, indicative of the weaker complex, are also extracted, but it was possible with titration data alone to unambiguously assign affinities to particular conformations due to its associated relative population being so low. The apparent agreement of NMR ROESY data with the MD results, however, seems to suggest the dominant conformations established from NMR data are those with the highest computed energies (i.e., most enthalpically favorable) as given in Fig. 3. However, the fitting parameters from this more physically relevant fit are much less robust than those obtained from the 1:1 model due to the number of and high correlation among fitting parameters. For this reason, we plotted binding constants obtained for the more robust 1:1 model in Fig. 1 (with SBX + 1 being the exception since the data could not be fit reasonably well with the simple model).
Reagents and solvents were purchased and used as received. Methyl 3-mercaptopropionate, cesium carbonate, and 2-mercaptoethan-1-ol were purchased from Alfa Aesar (Ward Hill, MA). Methyl thioglycolate was purchased from Pfaltz and Bauer, Inc. (Waterbury, CT). Methyl 4-sulfanylbutanoate was purchased from Enamine Ltd. (Kyiv, Ukraine). N-methyl-2-pyrrolidone was purchased from Applied Biosystems, Thermo Fisher (Grand Island, NY). 4-Mercaptobutan-1-ol and 3-mercaptopropan-1-ol were synthesized as previously described45,46 purified by flash column chromatography in silica gel cartridges with a Biotage Isolera purification system, and their purity assessed by 1H NMR and GC-MS47,48,49. Deuterated water (D2O) was purchased from Cambridge Isotope Laboratories, Inc. (Tewksbury, MA). Hexakis-, heptakis- and octakis-C6-brominated cyclodextrins were purchased from Arachem/Cyclodextrin Shop (Tilberg, Netherlands). Thiol-mediated displacement reactions on the C6-per-brominated cyclodextrins followed by methyl ester hydrolysis (where applicable) were accomplished with a modified version of the protocol originally published by Adam et al36. For the synthesis of mercaptoalkylcarboxylic acid methyl ester intermediates, thin layer chromatography (TLC) was used to monitor their production using Merck kieselgel 60-F254 sheets and detection accomplished with UV light (λ = 254 nm) in conjunction with development of color with Ce(SO4)2/(NH4)6Mo7O24·4H2O/H2SO4/H2O (CAM)50,51,52,53 and iodine vapor54,55,56. Centrifugation was performed in a Eppendorf centrifuge model 5810R at 3220 rpm for 5 min, using Falcon tubes (50 mL) purchased from VWR (Aurora, CO.). Solvents were removed using an IKA RV8 model rotary evaporator coupled to a KNF-Lab vacuum filtration pump and a VWR RS-232 cooling/circulating system (9 °C water/ethylene glycol 1:1).
1H NMR (600 MHz), 13C NMR (150 MHz) and 13C-DEPT NMR (150 MHz) spectra were recorded in D2O unless otherwise specified. Spectra were obtained using a Bruker Avance III 600 MHz instrument equipped with a Bruker TCI 5 mm cryoprobe (Bruker Biospin, Billerica, MA) at 30.0 ± 0.1 °C. NMR data for the synthesized compounds detailed below are reported as follows: chemical shift (δ) (parts per million, ppm); multiplicity: s (singlet), d (doublet), t (triplet), dd (doublet of doublets), td (triplet of doublets), m (multiplet), app (apparent) and br (broad); coupling constants (J) are given in Hertz (Hz). For the experimental data used in characterization of host:guest complexes, One-dimensional 1H NMR data were collected with water suppression by excitation sculpting with gradients57. For each 1-D experiment, 16 to 128 transients (depending on analyte concentration) were collected using 65,536 data points, a 4.0 s acquisition time, and a 1.0 s relaxation delay. Prior to the Fourier transformation the free induction decays (FIDs) were apodized with an exponential decay equivalent to 0.25 Hz line broadening. Two-dimensional ROESY spectra were acquired using 16,384 data points with 1024 increments, 8—32 scans for each increment, and a continuous wave spin lock with a 200 ms mixing time and a frequency of 3.57 kHz. Further details can be found in our previous publication25. Exact masses of the synthesized compounds below were obtained on a Bruker micrOTOF-Q III (Bruker Daltonics, Billerica, MA) equipped with an electrospray ionization (ESI) source. Time-of-flight mass spectroscopic detection was performed in negative ion mode. Mass calibration was performed in the final minutes of each run using Agilent’s ESI-L Low Concentration Tuning Mix (G1969-85,000). An Agilent 1260 Liquid Chromatograph (LC) equipped with an Atlantis T3 Reverse Phase Column (C18, 150 mm × 2.1 mm, 3 μm particle size, Waters, Milford, MA) was used for the separation of all compounds. The LC mobile phase consisted of water with 0.1% formic acid (A) and acetonitrile with 0.1% formic acid (B). The gradient profile started with 95% A for 2 min, ramped to 5% A at 18 min, held for 13.5 min, ramped quickly back to 95% over 0.5 min, and held for 10 min for column regeneration. This method was used for all samples except those for the vinyl tile surface study, for which the final 5% A was held for 18.5 min. Ten microliters of the liquid sample were introduced via an autosampler (Agilent B1329B) to the injection port.
Equilibrium binding constants between the CD family and fentanyl HCl were determined using 1H NMR titration experiments. In order to derive the most accurate binding constant, the cyclodextrin concentration [H]T (in units of M) satisfied the inequality [H]TK < 1 whenever possible. Initial studies required best-guess estimates of K since they were a priori unknown. Ultimately, this condition was satisfied for all compounds save for SBX + 1, where the required concentration yielded an unacceptable signal-to-noise ratio. Use of this non-ideal concentration leads to a larger uncertainty in the derived binding constant (vide infra). Due to the limited aqueous solubilities of the neutral compounds 5–7 (< 1 mg·mL-1 in all cases), all titrations were conducted at [H]T < 0.1 mM. Each of the modified CD was dissolved into deuterated water (D2O) to a total volume of 500 μL with 1 μL of acetonitrile added as a chemical shift reference point. Fentanyl HCl was then titrated into the solution. The total (bound to cyclodextrin + free in solution) fentanyl HCl concentration [G]T (in units of M) was varied to cover the range of concentrations between \(\frac{\left(0.2{\left[H\right]}_{T}+0.25\right)}{K}\) and \(\frac{\left(0.8{\left[H\right]}_{T}+4\right)}{K}\) in order to allow accurate extraction of binding constants. At each titration point, the proton chemical shifts (referenced to acetonitrile at 2.014 ppm) of the H1, H3, H5, and H6 protons on the CD were recorded; subtraction of the chemical shift obtained from the first data point, consisting of the cyclodextrin alone, yielded a difference between bound and free chemical shift \(\Delta {\delta }_{i}\) for each of reporter proton i. The resultant curves of \(\Delta {\delta }_{i}\) versus [G]T were then initially fit (Mathematica 10, Wolfram) using an non-linear least squares algorithm discussed previously (see additional mathematical details in Supporting Information)31. As in our previous work, the binding constant K and a chemical shift sensitivity factor \(\Delta {\updelta }_{c,i}\) for each reporter proton are the fitting parameters. Restricting K to a single value for each NMR titration curve yields a global binding constant incorporating data from all reporter protons58. All experiments were conducted in triplicate. Reported values of K for each compound are the average and standard deviation across the three experiments.
Molecular dynamics (MD) were performed with AMBER59 (version 12) employing the recent charges and parameters applied for cyclodextrin scaffolds60,61 and the GAFF force field for the various fentanyls62 in their protonated states63. Fentanyl charges were generated by AM1-BCC calculations55 in the program ANTECHAMBER. The program CHIMERA was employed for modeling the CD:fentanyl complexation processes64. The CD and CD:fentanyl complexes were solvated in a box of TIP3P water sufficient in size to have at least 15 Å of water between the solute and the solvent interface (~ 51 × 51 × 51 Å3 initial box size)65. To neutralize the systems, sufficient sodium (Na+) ions (typically one) were added to the system. The systems consisted of about 12,500 atoms (~ 4,100 water molecules). Each system was energy minimized using 250 steps of steepest descent and 1500 steps of conjugate gradient. Constant temperature and pressure dynamics (NPT) were performed on these minimized systems66,67. Periodic boundary conditions were used, and electrostatic interactions treated by particle mesh Ewald methods using a 9 Å cut-off in direct space with a 1 Å grid68. Bonds containing hydrogen were constrained using SHAKE and a time step of 2 fs was used for each simulation69. The systems were initially coupled to a heat bath at 100 K for the first 100 ps, then increased to 200 K for the next 100 ps, and finally raised to 300 K for the remainder of the simulation. Each simulation was performed for a total of 10 ns or 30 ns depending on the CD guest. Initial simulations involved only single trajectories, but it was universally observed that the βCD structure did not change significantly during the simulation, and that the root mean squared deviation of the non-hydrogen atoms relative to the average structure quickly plateaued by 1 ns. The first 2.5 ns of the 300 K dynamics were used for equilibration. To obtain more comprehensive data and statistical information, ten replicate simulations (10 ns each) were performed for complexes between all four charged fentanyls and βCD where the guest molecule was parallel to the wider CD rim. These replicate simulations were characterized by changing the seed value for generating the initial velocities for the system.
The free energy of binding between the cyclodextrins and fentanyls were estimated using the Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method from snapshots of the solvated trajectories70,71,72,73. MM-GBSA energy calculations were performed on replicate simulations and then averaged to obtain the average binding energy for a particular fentanyl:CD system. The binding free energy was estimated by the equation:
were each term, G is estimated as the sum of gas-phase molecular mechanics energy Egas and the solvation energy Gsol:
The contribution of entropy was neglected in these free energy calculations. The solvation free energy (Gsol) is the sum of the polar and nonpolar solvation energies of the molecules determined by solving the Generalized Born (GB) equation. The binding free energies for the complexes were calculated using the MMPBSA.py script in AMBER12 on snapshots from each 7.5 ns trajectory. The modified Onufriev-Bashford-Case-I GB (ib = 2) model was used for the calculation with a fentanyl salt concentration of 0.24 mM74. The surface tension used to calculate the nonpolar contribution to the free energy of solvation was 0.0072 kcal mol-1 Å-2.
In this work, we have studied the inclusion complexes between a number of chemically modified CDs with the synthetic opioid fentanyl. The evaluated panel of CDs exhibited either anionic thioalkylcarboxylic acid side chains or neutral thioalkyl alcohols side chains on their top rim. The charged CDs were found to have extended hydrophobic interiors kept open by mutual electrostatic side chain repulsion, thus providing a sufficient binding pocket for the opioid. MD simulations of the SBX:fentanyl inclusion complexes were found to be consistent with binding affinity trends extracted from NMR experiments. That is, binding strength increases as the length of the modified arm increases for the three charged CDs studied. For CDs with neutral thioalkyl side chains, binding affinities were observed to decrease with increasing arm length, likely a result of competitive binding of the more hydrophobic neutral side chains also seeking the hydrophobic SBN interiors. MD simulations also supported this competitive binding phenomenon revealing that longer arms tended to interact more strongly with the hydrophobic interior. For the panel of promising charged SBXs, the maximum fentanyl binding affinity was found to be Ka = 66,500 M-1 for SBX + 1. To our knowledge this is the largest equilibrium binding affinity reported for CD:fentanyl complexes. We believe that these CD scaffolds could be developed into promising, novel medical countermeasures candidates against these potential deadly opioids. Another important avenue for future investigation but is not covered presently: the potential power of simultaneously fitting Job plot and titration data to extract more confident values of highly correlated fitting parameters. Our group is also focusing on exploiting these compounds for environmental scavenging of fentanyls as well as for detection of these opioids in complex biological forensic samples.
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Clotz , MA & Nahata , MC Clinical uses of fentanyl , sufentanil , and alfentanil .Clin.Pharm.10, 581–593 (1991).
Miró, Ò., Gil, V. & Peacock, W. F. Morphine in acute heart failure: Good in relieving symptoms, bad in improving outcomes. J. Thorac. Dis. 9, 1831–1834 (2017).
Chen, A. & Ashburn, M. A. Cardiac effects of opioid therapy. Pain Med. 16, S27–S31 (2015).
Lipinski, P. F. J. et al. Fentanyl family at the mu-opioid receptor: Uniform assessment of binding and computational analysis. Molecules 24, 740 (2019).
Article CAS PubMed PubMed Central Google Scholar
Lipinski, P. F. J., Jaronczyk, M., Dobrowolski, J. C. & Sadlej, J. Molecular dynamics of fentanyl bound to μ-opioid receptor. J. Mol. Model 25, 144 (2019).
Gladden, R. M., Martinez, P. & Seth, P. Fentanyl law enforcement submissions and increases in synthetic opioid-involved overdose deaths - 27 states, 2013–2014. MMWR Morb. Mortal. Wkly. Rep. 65, 837–843 (2016).
Marinetti, L. J. & Ehlers, B. J. A series of forensic toxicology and drug seizure cases involving illicit fentanyl alone and in combination with heroin, cocaine or heroin and cocaine. J. Anal. Toxicol. 38, 592–598 (2014).
Article CAS PubMed Google Scholar
Han, Y. et al. The rising crisis of illicit fentanyl use, overdose, and potential therapeutic strategies. Transl. Psychiat. 9, 282 (2019).
Valdez, C. A. Gas chromatography-mass spectrometry analysis of synthetic opioids belonging to the fentanyl class: A review. Crit. Rev. Anal. Chem. https://doi.org/10.1080/10408347.2021.1927668 (2021).
Stanley, T. H. Fentanyl. J. Pain Symptom Manag. 29, S67–S71 (2005).
Stanley, T. H. The history and development of the fentanyl series. J. Pain Symptom Manag. 7, S3–S7 (1992).
Janssen, P. A. J. & Eddy, N. B. Compounds related to pethidine-IV. New general chemical methods of increasing the analgesic activity of pethidine. J. Med. Chem. 2, 31–45 (1960).
Valdez, C. A., Leif, R. N. & Mayer, B. P. An efficient, optimized synthesis of fentanyl and related analogs. PLoS ONE 9, e108250 (2014).
Article ADS PubMed PubMed Central Google Scholar
Walz, A. J. & Hsu, F.-L. An operationally simple synthesis of fentanyl citrate. Org. Prep. Proc. Int. 49, 467–470 (2017).
Wax, P. M., Becker, C. E. & Curry, S. C. Unexpected, “Gas” casualties in Moscow: A medical toxicology perspective. Ann. Emerg. Med. 41, 700–705 (2003).
Riches, J. R., Read, R. W., Black, R. M., Cooper, N. J. & Timperley, C. M. Analysis of clothing and urine from Moscow theatre siege casual-ties reveals carfentanil and remifentanil use. J. Anal. Toxicol. 36, 647–656 (2012).
Article CAS PubMed Google Scholar
Pardo, B. et al. The Future of Fentanyl and Other Synthetic Opioids (RAND Corporation, 2019).
Valdez, C. A., Leif, R. N., Hok, S. & Hart, B. R. Analysis of chemical warfare agents by gas chromatography-mass spectrometry: Methods for their direct detection and derivatization approaches for the analysis of their degradation products. Rev. Anal. Chem. 37, 1–26. https://doi.org/10.1515/revac-2017-0007 (2018).
Yeung, D. T., Bough, K. J., Harper, J. R. & Platoff, G. E. Jr. National Institutes of Health (NIH) executive meeting summary: Developing medical countermeasures to rescue opioid-induced respiratory depression (a Trans-Agency Scientific Meeting)-August 6/7, 2019. J. Med. Toxicol. 16, 87–105 (2020).
France, C. P. et al. Countermeasures for preventing and treating opioid overdose. Clin. Pharmacol. Therap. 109, 578–590 (2021).
Deng, C. L., Murkli, S. L. & Isaacs, L. D. Supramolecular hosts as in vivo sequestration agents for pharmaceuticals and toxins. Chem. Soc. Rev. 49, 7516–7532 (2020).
Article CAS PubMed PubMed Central Google Scholar
Diaz-Gil, D. et al. A novel strategy to reverse general anesthesia by scavenging with the acyclic cucurbit[n] uril-type molecular container calabadion 2. Anesthesiology 125, 333–345 (2016).
Article CAS PubMed Google Scholar
Brockett, A. T. et al. In vitro and in vivo sequestration of methamphetamine by a Sulfated Acyclic CB[n]-type receptor. Chem. Eur. J. 27, 17476–17486 (2021).
Article CAS PubMed Google Scholar
Thevathasan, T. et al. Calabadion 1 selectively reverses respiratory and central nervous system effects of fentanyl in rat model. Br. J. Anaesth. 125, e140–e147 (2020).
Article CAS PubMed PubMed Central Google Scholar
Mayer, B. P., Kennedy, D. J., Lau, E. Y. & Valdez, C. A. Solution-state structure and affinities of cyclodextrin: Fentanyl complexes by nuclear magnetic resonance spectroscopy and molecular dynamics simulation. J. Phys. Chem. B 120, 2423–2433 (2016).
Article CAS PubMed Google Scholar
Del Valle, E. M. M. Cyclodextrins and their uses: A review. Process Biochem. 39, 1033–1046 (2004).
Crini, D. Review: A history of cyclodextrins. Chem. Rev. 114, 10940–10975 (2014).
Article CAS PubMed Google Scholar
Easton, C.J., Lincoln, S.F. Modified cyclodextrins: Scaffolds and templates for supramolecular chemistry:1–283 (Imperial College Press, 1999).
Braga, S. S. Cyclodextrins: Emerging medicines of the new millennium. Biomolecules 9, 801 (2019).
Article CAS PubMed PubMed Central Google Scholar
Liu, Y. et al. Research progress on synthesis and applications of cyclodextrin polymers. Molecules 26, 1090 (2021).
Article CAS PubMed PubMed Central Google Scholar
Jiang, H., Xu, Y., Na, L., Jin, R. & Zhang, S. UV-vis spectral analysis of inclusion complexes between beta-cyclodextrin and aromatic/aliphatic guest molecules. Curr. Drug Discov. Technol. 5, 173–176 (2008).
Article CAS PubMed Google Scholar
Garnero, C., Aiassa, V. & Longhi, M. Sulfamethoxazole: Hydroxypropyl-b-cyclodextrin complex: Preparation and characterization. J. Pharm. Biomed. Anal. 63, 74–79 (2012).
Article CAS PubMed Google Scholar
Mathiron, D. et al. A combined approach of NMR and mass spectrometry techniques applied to the a-Cyclodextrin/Moringin complex for a novel bioactive formulation. Molecules 23, 1714 (2018).
Article PubMed PubMed Central Google Scholar
Rekharsky, M. V. & Inoue, Y. Complexation thermodynamics of cyclodextrins. Chem. Rev. 98, 1875–1918 (1998).
Article CAS PubMed Google Scholar
Mayer, B. P., Albo, R. L. F., Hok, S. & Valdez, C. A. NMR spectroscopic investigation of inclusion complexes between cyclodextrins and the neurotoxin tetramethylenedisulfotetramine. Magn. Res. Chem. 50, 229–235 (2012).
Adam, J. M. et al. Cyclodextrin-derived host molecules as reversal agents for the neuromuscular blocker rocuronium bromide: Synthesis and structure-activity relationships. J. Med. Chem. 45, 1806–1816 (2002).
Article CAS PubMed Google Scholar
Cameron, K. S. et al. Modified γ-Cyclodextrins and their rocuronium complexes. Org. Lett. 4, 3403–3406 (2002).
Article CAS PubMed Google Scholar
Nag, K. et al. Sugammadex: A revolutionary drug in neuromuscular pharmacology. Anesth. Essays Res. 7, 302–306 (2013).
Article PubMed PubMed Central Google Scholar
Mirakhur, R. K. Sugammadex in clinical practice. Anesthesia 64, 45–54 (2009).
Bom, A. et al. A novel concept of reversing neuromuscular block: Chemical encapsulation of rocuronium bromide by a cyclodextrin-based synthetic host. Angew. Chem. Int. Ed. Engl. 41, 266–270 (2002).
Fielding, L. Determination of association constants (Ka) from solution NMR data. Tetrahedron 56, 6151–6170 (2000).
Chen, W., Chang, C.-E. & Gilson, M. K. Calculation of cyclodextrin binding affinities: energy, entropy, and implications for drug design. Biophys. J. 87, 3035–3039 (2004).
Article ADS CAS PubMed PubMed Central Google Scholar
Saha, S., Roy, A., Roy, K. & Roy, M. N. Study to explore the mechanism to form inclusion complexes of β-cyclodextrin with vitamin molecules. Sci. Rep. 6, 35764 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Rajbanshi, B. et al. Study to probe subsistence of host-guest inclusion complexes of α and β-Cyclodextrins with biologically potent drugs for safety regulatory dischargement. Sci. Rep. 8, 13031 (2018).
Article ADS PubMed PubMed Central Google Scholar
Carneiro, S. B. et al. Cyclodextrin-drug inclusion complexes: In vivo and in vitro approaches. Int. J. Mol. Sci. 20, 642 (2019).
Article CAS PubMed PubMed Central Google Scholar
Baker, B. R., Querry, M. V., Bernstein, S., Safir, S. R. & Subbarow, Y. Biotin .IV. A second synthesis of 2-(δ-carboxybutyl)thiophane-3,4-dicarboxylic acid. J. Org. Chem. 12, 167–173 (1947).
Article CAS PubMed Google Scholar
Ekeroth, J., Borgh, A., Konradsson, P. & Liedberg, B. Synthesis and monolayer characterization of phosphorylated amino acid analogs. J. Colloid Interface Sci. 254, 322–330 (2002).
Article ADS CAS PubMed Google Scholar
Valdez, C. A., Leif, R. N. & Alcaraz, A. Effective methylation of phosphonic acids related to chemical warfare agents mediated by trimethyloxonium tetrafluoroborate for their qualitative detection and identification by gas chromatography-mass spectrometry. Anal. Chim. Acta 933, 134–143 (2016).
Article CAS PubMed Google Scholar
Valdez, C. A. et al. Structural modification of fentanyls for their retrospective identification by gas chromatographic analysis using chloroformate chemistry. Sci. Rep. 11, 22489 (2021).
Article ADS CAS PubMed PubMed Central Google Scholar
Valdez, C. A. et al. Methylation protocol for the retrospective detection of isopropyl-, pinacolyl- and cyclohexylmethylphosphonic acids, indicative markers for the nerve agents sarin, soman and cyclosarin, at low levels in soils using EI-GC-MS. Sci. Total Environ. 683, 175–184 (2019).
Article ADS CAS PubMed Google Scholar
Wnuk, S. F. et al. Doubly homologated dihalovinyl and acetylene analogues of adenosine: Synthesis, interaction with S-adenosyl-L-homocysteine hydrolase, and antiviral and cytostatic effects. J. Med. Chem. 43, 1180–1186 (2000).
Article CAS PubMed Google Scholar
Wnuk, S. F. et al. Sugar-modified conjugated diene analogues of adenosine and uridine: Synthesis, interaction with S-adenosyl-L-homocysteine hydrolase, and antiviral and cytostatic effects. J. Med. Chem. 45, 2651–2658 (2002).
Article CAS PubMed Google Scholar
Showalter, B. M. et al. Diazeniumdiolate ions as leaving groups in anomeric displacement reactions: A protection-deprotection strategy for ionic diazeniumdiolates. J. Am. Chem. Soc. 127, 14188–14189 (2005).
Article CAS PubMed Google Scholar
Valdez, C. A. et al. Hydrolytic reactivity trends among potential prodrugs of the O2-glycosylated diazeniumdiolate family. Targeting nitric oxide to macrophages for antileishmanial activity. J. Med. Chem. 51, 3961–3970 (2008).
Article CAS PubMed PubMed Central Google Scholar
Koziol, L. et al. Toward a small molecule, biomimetic carbonic anhydrase model: Theoretical and experimental investigations of a panel of zinc(II) aza-macrocyclic catalysts. Inorg. Chem. 51, 6803–6812 (2012).
Article CAS PubMed Google Scholar
Floyd, W. C. et al. Evaluation of a carbonic anhydrase mimic for industrial carbon capture. Environ. Sci. Technol. 47, 10049–10055 (2013).
Article ADS CAS PubMed Google Scholar
Hwang, T. L. & Shaka, A. J. Water suppression that works. Excitation sculpting using arbitrary wave-forms and pulsed-field gradients. J. Magn. Reson. Ser. A 112, 275–279 (1995).
Article ADS CAS Google Scholar
Tordarsson, P. Determining association constants from titration experiments in supramolecular chemistry. Chem. Soc. Rev. 40, 1305–1323 (2011).
Case, D. A. et al. The amber biomolecular simulation programs. J. Comput. Chem. 26, 1668–1688 (2005).
Article CAS PubMed PubMed Central Google Scholar
Kalydi, E. et al. Fully symmetric cyclodextrin polycarboxylates: How to determine reliable protonation constants from NMR titration data. Int. J. Mol. Sci. 23, 14448 (2022).
Article CAS PubMed PubMed Central Google Scholar
Cézard, C., Trivelli, S., Aubry, F., Djedaïni-Pilard, F. & Dupreadeu, F. Y. Molecular dynamics studies of native and substituted cyclodextrins in different media: 1. Charge derivation and force field performances. Phys. Chem. Chem. Phys. 13, 15103–15121 (2011).
Thurlkill, R. L., Cross, D. A., Scholtz, J. M. & Pace, N. pKa of fentanyl varies with temperatures: Implications for acid-base management during extremes of body temperature. J. Cardiothor. Vasc. Anesth. 19, 759–762 (2005).
Wang, J., Wolf, R. M., Caldwell, J. W., Kollman, P. A. & Case, D. A. Development and testing of a general amber force field. J. Comput. Chem. 25, 1157–1174 (2004).
Article CAS PubMed Google Scholar
Jakalian, A., Bush, B. L., Jack, D. B. & Bayly, C. I. Fast, efficient generation of high-quality atomic charges. AM1-BCC model: I. Method. J. Comput. Chem. 21, 132–146 (2000).
Pettersen, E. F. et al. UCSF chimera: A visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612 (2004).
Article CAS PubMed Google Scholar
Jorgensen, W. L., Chandrasekhar, J., Madura, J. D., Impey, R. & Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 79, 926 (1983).
Article ADS CAS Google Scholar
Miller, B. R. III. et al. MMPBSA.py: An efficient program for end-state free energy calculations. J. Chem. Theory Comput. 8, 3314–3321 (2012).
Article CAS PubMed Google Scholar
Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A. & Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 81, 3684–3690 (1984).
Article ADS CAS Google Scholar
Darden, T., York, D. & Pedersen, L. Particle mesh Ewald: An N log (N) method for ewald sums in large systems. J. Chem. Phys. 98, 10089–10092 (1993).
Article ADS CAS Google Scholar
Ryckaert, J.-P., Cicotti, G. & Berendsen, H. J. C. Numerical integration of the cartesian equations of motion of a system with constraints: Molecular dynamics of n-alkanes. J. Comput. Phys. 23, 327–341 (1977).
Article ADS CAS Google Scholar
Kollman, P. A. et al. Calculating structures and free energies if complex molecules: Combining molecular mechanics and continuum models. Acc. Chem. Res. 33, 889–897 (2000).
Article CAS PubMed Google Scholar
Satcher, J. H. Jr. et al. Modeling, synthesis and characterization of zinc containing carbonic anhydrase active site mimics. Energy Procedia 4, 2090–2095 (2011).
Wong, S. E. et al. Designing small-molecule catalysts for CO2 capture. Energy Procedia 4, 817–823 (2011).
Onufriev, A., Bashford, D. & Case, D. A. Exploring protein native states and large scale conformational changes with a modified generalized born model. Proteins: Struc Func. Bioinf. 55, 383–394 (2004).
This document (LLNL-JRNL-844442) was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes.
This work was performed under the auspices of the U. S. Department of Energy by Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. The work was funded fully by a Laboratory Directed Research and Development Grant (14-ERD-048) to C. A. V.
Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
Brian P. Mayer, Daniel J. Kennedy, Edmond Y. Lau & Carlos A. Valdez
Nuclear and Chemical Sciences Division, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, CA, L-090 94550, USA
Brian P. Mayer, Daniel J. Kennedy, Edmond Y. Lau & Carlos A. Valdez
Forensic Science Center, Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
Brian P. Mayer, Daniel J. Kennedy & Carlos A. Valdez
Biosciences and Biotechnology Division, Lawrence Livermore National Laboratory, Livermore, CA, 94550, USA
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
C. A. V. synthesized and purified all compounds used in this work. B. P. M. and D. J. K. conducted all the 1D and 2D NMR experiments. E. Y. L. conducted all the molecular dynamic simulations and docking experiments. C. A. V and B. P. M. wrote the manuscript. All authors reviewed the manuscript.
Correspondence to Carlos A. Valdez.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Mayer, B.P., Kennedy, D.J., Lau, E.Y. et al. Evaluation of polyanionic cyclodextrins as high affinity binding scaffolds for fentanyl. Sci Rep 13, 2680 (2023). https://doi.org/10.1038/s41598-023-29662-1
DOI: https://doi.org/10.1038/s41598-023-29662-1
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Formwork accessories Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.