Design Strategies for Oxytocin Analogues. OT, a neurohormone conserved across species for 700 million years,11,26 plays a major role in many fundamental physiological processes in both the central and peripheral nervous system.11,18,26,27 Its primary role in behavior is “to facilitate stability in changing environments, rooted in its evolutionary origins as a system to promote survival.”28 This biological significance led to sustained interest in its pharmaceutical application to a range of disorders,27 including pain. Its extensive conservation suggests that at least some of its native structure is advantageous. We made minimal structural changes to address the pharmacological limitations of oxytocin: lack of receptor specificity and poor bioavailability. We aimed to optimize its druggability and elicit antinociception in mouse models.
Specifically, we modified the native amino acid sequence and glycosylated the C-terminus of the peptide chain. Replacing the labile disulfide bond with a more robust amide bond produced highly potent glycopeptides, which are exquisitely selective for the OT receptor, brain penetrant, and antinociceptive after intravenous administration in mice. The amide bond was also used to vary the size of the resultant macrocyclic ring. Our synthetic modifications drew on established design strategies in peptide drug design, such as macrocyclization and glycosylation. Macrocyclization is used for increasing stability towards proteolysis29,30 and providing an entropic advantage in receptor binding (as compared to more flexible linear peptides).31 Glycosylation slows the metabolism of peptide drugs, increasing their bioavailability, and promotes BBB penetration.32–38 It is hypothesized that glycosylation increases receptor selectivity because of “membrane hopping” due to “biousian behavior,” and that the carbohydrate acts as a shield to protect the peptide from enzymatic cleavage.32,33,38,39
Receptor Potencies in CHO Cells. All novel oxytocin derivatives were tested for their potency and selectivity at the Oxytocin Receptor (OTR) and closely related Arginine-Vasopressin 1a/1b Receptors (AVP1a/b). They were evaluated using an In Cell Western assay for ERK MAPK activation in CHO cells expressing human OTR/AVP. We found that oxytocin (OT) and vasopressin (VP) control peptides activated the cells with high potency, as expected (Table 1). The 8 OT derivatives showed a range of receptor potencies at the OTR, from the very high (8) to low (1–2) with ligands performing in between. Several compounds improved on the potency of OT, demonstrating that our modifications can enhance both pharmacodynamic and pharmacokinetic performance (3, 5, 6, and 8). Notably, all derivatives showed no activity at the AVP1a/1b receptors, suggesting that they may be selective in vivo and avoid blood pressure-related side effects caused by AVP activation (Table 1).
Table 1
In Vitro Molecular Pharmacology of Novel Oxytocin Derivatives. Bold amino acids indicate cyclization site. Dap = L-1,3-diaminopropionic acid, Aib = α-amino-isobutyric acid. NA = no activity. Potency values at the 3 receptors reported as the mean ± SEM of N = 3 independent experiments.
|
Structure
|
MW / RTHPLC
|
OTR EC50 (nM)
|
AVP1a EC50 (nM)
|
AVP1b EC50 (nM)
|
OT
|
C[Cys-Tyr-Ile-Gln-Asn-Cys]-Pro-Leu-Gly-NH2
|
1007.19 / —
|
2.3 ± 1.7
|
–
|
–
|
VP
|
C[Cys-Tyr-Phe-Gln-Asn-Cys]-Pro-Arg-Gly-NH2
|
1008.2
|
–
|
32 ± 23
|
11 ± 8
|
1
|
C[Asp-Tyr-Ile-Gln-Asn-Dap]-Pro-Leu-Gly-NH2
|
986.61 / 9.82
|
3800 ± 3100
|
NA
|
NA
|
2
|
C[Glu-Tyr-Ile-Gln-Asn-Dap]-Pro-Leu-Gly-NH2
|
1178.71 / 9.29
|
7000 ± 3000
|
NA
|
NA
|
3
|
C[Asp-Tyr-Ile-Gln-Asn-Dap]-Pro-Leu-Ser(Lact)-NH2
|
1340.76 / 8.42
|
0.54 ± 0.48
|
NA
|
NA
|
4
|
C[Asp-Tyr-Ile-Gln-Asn-Dap]-Pro-Leu-Gly-Ser(β-Lact)-NH2
|
1398.69 / 9.50
|
160 ± 160
|
NA
|
NA
|
5
|
C[Asp-Tyr-Ile-Gln-Asn-Dap]-Pro-Aib-Ser(β-Glc)-NH2
|
1150.64 /8.84
|
0.28 ± 0.17
|
NA
|
NA
|
6
|
C[Glu-Tyr-Ile-Gln-Asn-Dap]-Pro-Leu-Ser(β-Glc)-NH2
|
1197.72 / 9.66
|
0.73 ± 0.49
|
NA
|
NA
|
7
|
C[Glu-Tyr-Ile-Gln-Asn-Dap]-Pro-Aib-Ser(β-Glc)-NH2
|
1164.67 / 8.50
|
220 ± 110
|
NA
|
NA
|
8
|
C[Glu-Tyr-Ile-Gln-Asn-Dap]-Pro-Aib-Gly-Ser(β-Glc)-NH2
|
1221.67 / 8.74
|
0.06 ± 0.01
|
NA
|
NA
|
Based on their activity in cells, we next selected oxytocin derivatives to evaluate their antinociceptive ability in vivo. We chose 1 as a comparison since it is not glycosylated and only contains the stabilized amide bond, and selected 8 as the most potent glycoside (38-fold more potent than OT). All drugs were delivered by the intravenous route to mice. When we tested antinociceptive activity by the tail flick test, we found robust and efficacious pain relief for both 1 and 8, with 8 slightly higher than 1. Vehicle and oxytocin control produced no effect, while morphine produced robust antinociception, validating the assay. These results strongly suggest that our stabilized oxytocin derivatives can penetrate the brain to produce pain relief.
Conformational Studies of Oxytocin Analogues by NMR. Following the cell-based efficacy data from in vitro assays, we selected compound 8 for extensive NMR characterization. Importantly, researchers in the past have found the 20-membered ring of OT to be essential for its binding and activity.40,41 To our knowledge, only one other 21-membered analogue has produced OTR agonism, and this alkene-containing derivative required much higher dosage than our analogue does to elicit effects.42 The results from analogue 8, suggest the 20-membered ring motif is not essential, since 8 has greater specificity for and activity at OTR despite its 21-membered ring. We obtained NMR spectra in solution corresponding to data in the prior aqueous NMR studies of oxytocin.43–45 High-field (1H and 13C 800 MHz and 200 MHz, respectively) 1D and 2D (1H-13C-HSQC, 1H-13C-HMBC, 1H-1H-COSY, 1H-1H NOESY, 1H-1H TOCSY, and 1H-15N HSQC) NMR studies were used to unambiguously confirm the structure of the compound and obtain sequence specific resonance assignments. Water suppression was achieved with excitation sculpting in the pulse sequence. Sequential NOEs (nuclear Overhauser effects) between the HG/HG’ of Glu1 and HG in Dap6 confirm the cyclization between those amino acids (Fig. 2). Sequential HN-HA NOEs were not significantly different in strength as compared to the sequential HN-HN NOEs. Chemical shifts for Ha, Ca, and backbone carbonyl carbons were compared to Wishart’s Chemical-Shift Index (CSI) to predict secondary structure46,47 Variable temperature 1H NMR (VT NMR) data collected at five temperatures was used to quantify chemical shift changes of the amide protons (temperature coefficient, TC, DdNH/DT), which are used to predict the presence of intramolecular hydrogen bonding (IMHB).48,49 TC calculations showed evidence of intramolecular hydrogen bonding between residues in the cyclic moiety, although that evidence is weakened by the CSI and NOESY analysis. Preliminary TC data in ionic (SDS) micelles provides much greater evidence for these interactions, although the extent of the contribution from the charge has yet to be determined (manuscript in preparation). Additional spectra, assignments, TC calculations, and CSI analysis can be found in the supplemental information.
Molecular Modeling Studies. Published studies based on molecular modeling, photoaffinity labeling, and site-directed mutagenesis indicate that the oxytocin macrocycle Cys1-Tyr2-Ile3-Gln4-Asn5-Cys6 binds into the upper third of the oxytocin receptor pocket at the transmembrane helices TM3, TM4, and TM6. The linear tail Pro-Leu-Gly-amide is bound closer to the surface and interacts with TM2 and TM3.50,51 Oxytocin-like biological activity (e.g. selectivity for OT vs AVP) is marked by the presence of two conserved intramolecular hydrogen bonds in the ligand backbone. This bonding system produces a Type I β-turn within the macrocyclic ring (see Supplemental Info for details).52,53 In addition, the activity of oxytocin is also associated with more compact structure as indicated with radius of gyration (Rg) centered at 5Å and below.53 In our studies all analogues were modeled using MOE®54 by building the structures of the compounds guided by either crystal structure of deamino-oxytocin,52 oxytocin,20 or by MOE® based energy of each amino acid sequence. Each structure was subjected to charge correction, energy minimization and conformational search using the LowModeMD method available in MOE®. Analysis of each ensemble for the presence of β-turn conformers, Rg and dE (i.e. strain energy of the conformation relative to the lowest energy conformation) across all the conformers showed that they all retained an accessible Type I β-turn. Descriptors of solubility, lipophilicity and total polar surface area were calculated for all analogues under evaluations. See supplementary methods for details of the calculations.
Stability Assessments. The stability of these compounds was assessed in vitro in artificial cerebrospinal fluid (aCSF) to verify the utility of our structural modifications (see Supplemental Information for data). Substitution of the disulfide cyclization with the amide linker shows a significant increase in tested compounds as compared to oxytocin. Interestingly, unlike most (glyco)peptides assessed in our lab, the half-lives of the glycosylated analogues was not statistically different than their un-glycosylated “partners,” which indicates the amide cyclization substitution is a greater contributing factor to stability in vitro. The candidates were also studied in vivo to assess their serum stability and BBB penetration. The in vivo results are also different than most comparisons of glycosylated/un-glycosylated peptides in our lab. In the case of the studied OT analogues, the glycosylated peptide was slightly less stable in vivo, degrading initially to the unglycosylated peptide. However, in every case the glycosylated “partner” showed greater CSF concentration than the corresponding un-glycosylated peptide, but this improvement was not always statistically significant. BBB penetration was probed by measuring CSF concentration using microdialysis and MS2. The tested OT compounds showed a peak in concentration around 10-20 minutes, and the CSF levels were sufficiently high to suggest BBB penetration, consistent with the in vivo antinociception studies.
Discussion and Outlook. OT remains a poor drug candidate, however, 700 million years of evolutionary conservation across species and almost 100 years of clinical applications in humans support the merit of much of its structure. Following this logic and guidance from Sir Francis Bacon who taught, “Nature to be commanded must be obeyed,”55 we made minimal synthetic modifications. OT is known for its short in vivo half-life, primarily due to its unstable disulfide bond. Replacement of the labile S–S bond with the more robust lactam (HN–C=O) increased the in vitro serum stability from 1.4 ± 0.6 min to 25.8 ± 0.9 min, and glycosylation increased the in vitro stability to 34.5 ± 0.9 min. The increased BBB penetration observed with the OT glycosides seems largely due to the increased stability of the lactam rings. (See SI for details.)
One of the most significant pharmacodynamic shortcomings of OT is its lack of receptor selectivity. It is well known that promiscuous binding can lead to unwanted and potentially dangerous side effects.56,57 Since OT binds to a wide range of receptors, including those involved in physiological regulatory mechanisms we have no interest in modulating (e.g. V1a receptor binding has been linked to bradycardia and body temperature changes after peripheral OT administration)41, we aim to design analogues that are selective for OTR. Additional support for this aim comes from Brackley et al., who found that V1aR cross-activation from high doses of their OT agonist compounds compromised OIRD reversal in their study of OTR agonism for fentanyl-induced OIRD.6
NMR data of 8 in aqueous solution do not indicate a preference for any specific secondary structure. Increasing the intensity of the spectra produces signals that indicate conformational isomers, which is expected for peptides in solution and especially for those containing proline residues.58 Glycopeptide 8 showed changes in TC (ΔδNH/ΔT) that indicate the potential for some structural preference, but this evidence is challenged by the CSI analysis for both proton and carbon spectra (see Supplemental Information). Conformational data of OT in aqueous solution has been published, and the receptor-bound conformation of OT•OTR was recently determined via cryo-EM.19 The crystal structure of OT bound to OTR has been recently determined via cryo-EM.20,21 Further NMR studies in micelles (both charged and neutral) will be used and compared with previous literature examining OT in nonpolar solvents (e.g., CDCl3, DMSO)41 and in charged micelles.59
Several of our OT analogues showed increased potency as compared to OT, demonstrating that our modifications can improve both pharmacodynamic and pharmacokinetic performance. All of the compounds produced no activity at the AVP1a/1b receptors, suggesting that they may be selective in vivo and avoid side effects caused by AVP activation. Antinociception was evaluated using the 52°C tail flick test in analogue 1 and its glycosylated “partner” 8, which showed robust, efficacious pain relief. Intravenous delivery in mice required BBB penetration, and the antinociceptive effects suggest that penetration was successful. Serum stability studies in vivo support this belief from increased dialysate levels 10-20 minutes after administration. These results strongly suggest that our stabilized oxytocin derivatives can penetrate the brain to produce antinociceptive effects. The half-lives of several glycopeptides were measured both in vitro and in vivo and showed no statistically significant differences between the glycosylated and un-glycosylated peptides of the same sequence, indicating that the amide replacement of the labile disulfide bond is much more significant in stabilizing the compounds than the glycosylation state.
An important consideration of novel pain therapeutics is avoiding the negative side effects of opioids, including reinforcing behavior. OT has been investigated as a potential OUD treatment in animals. Exogenous OT administration has been shown to mitigate classical addiction metrics, including self-administration, withdrawal and tolerance.14-17,60 More recently, clinical trials have lacked consistency and yielded mixed results,14 though strong evidence indicates that intranasal administration of OT leads to increased compliance with therapeutic interventions, addressing a well-established barrier to effective OUD/SUD treatment. OT was found to be a safe, effective attenuation for nonopioid respiratory depression in patients at increased risk for OIRD.6,61 Cardiorespiratory depression caused by fentanyl in rats was dose-dependently reversed by OTR activation, both with OT and with a nonpeptide partial agonist in a recent study by Brackley et al.6
Preliminary data will guide our future work generating additional efficacious analogues by consideration of conformational analysis, cell-based potencies, and extent of BBB penetration to guide our design efforts. This data suggests that the 20-membered ring in OT is not required for its activity at OTR, contributing to the growing insights into binding and conformation, and will inform our future drug design efforts. We will modify the structure of 8 and then study its receptor binding affinity, antinociceptive activity, and serum stability. NMR studies (both solution and micellar conditions) will be used to provide additional insight into in silico modeling efforts. Ultimately, in vivo studies will assess the addictive potential of these compounds to confirm that their administration yields no reward that might lead to overuse and/or misuse analogous to that which is associated with opioids.
Despite decades of efforts, chronic pain treatment remains challenging to researchers and clinicians alike.62 Growing insight into pain perception mechanisms are aiding the search for effective treatments. However, prescription opioids remain the primary approach to pain treatment, and their use is irrefutably dangerous. It is imperative for researchers to discover alternative non-opioid analgesics to address the persistent problem of pain without exacerbating the opioid epidemic.