Multinuclear 1H/13C/15N chemical shift assignment of therapeutic octreotide acetate performed at natural abundance.

Octreotide acetate, the active pharmaceutical ingredient in the long-acting release (LAR) drug product Sandostatin®, is a cyclic octapeptide that mimics the naturally occurring somatostatin peptide hormone. Modern NMR can be a robust analytical method to identify and quantify octreotide molecules. Previous 1H chemical shift assignments were mostly performed in organic solvents, and no assignments for heteronuclear 13C, 15N, and aromatic 1H nuclei are available. Here, using state-of-the-art 1D and 2D homo- and heteronuclear NMR experiments, octreotide was fully assigned, including water exchangeable amide protons, in aqueous buffer except for 13CO and 15NH of F1, 15NH of C2, and 15NζHζ of K5 that were not observed because of water exchange or conformational exchange. The solution NMR spectra were then directly compared with 1D 1H/13C/15N solid-state NMR (SSNMR) spectra showing the potential applicability of 13C/15N SSNMR for octreotide drug product characterization.

Proton Affinity and Conformational Integrity of a 24-Atom Triazine Macrocycle across Physiologically Relevant pH.

For 24-atom triazine macrocycles, protonation of the heterocycle leads to a rigid, folded structure presenting a network of hydrogen bonds. These molecules derive from dynamic covalent chemistry wherein triazine monomers bearing a protected hydrazine group and acetal tethered by the amino acid dimerize quantitatively in an acidic solution. Here, lysine is used, and the product is a tetracation. The primary amines of the lysine side chains do not interfere with quantitative yields of the desired bis(hydrazone) at concentrations of 5-125 mg/mL. Mathematical modeling of data derived from titration experiments of the macrocycle reveals that the pKa values of the protonated triazines are 5.6 and 6.7. Changes in chemical shifts of resonances in the 1H NMR spectra corroborate these values and further support assignment of the protonation sites. The pKa values of the lysine side chains are consistent with expectation. Upon deprotonation, the macrocycle enjoys greater conformational freedom as evident from the broadening of resonances in the 1H and 13C NMR spectra indicative of dynamic motion on the NMR time scale and the appearance of additional conformations at room temperature. While well-tempered metadynamics suggests only a modest difference in accessible conformational footprints of the protonated and deprotonated macrocycles, the shift in conformation(s) supports the stabilizing role that the protons adopt in the hydrogen-bonded network.

Adaptation of Empirical Methods to Predict the LogD of Triazine Macrocycles.

Octanol/water partition coefficients guide drug design, but algorithms do not always accurately predict these values. For cationic triazine macrocycles that adopt a conserved folded shape in solution, common algorithms fall short. Here, the logD values for 12 macrocycles differing in amino acid choice were predicted and then measured experimentally. On average, AlogP, XlogP, and ChemAxon predictions deviate by 0.9, 2.8, and 3.9 log units, with XlogP overestimating lipophilicity and AlogP and ChemAxon underestimating lipophilicity. Importantly, however, a linear relationship (R2 > 0.98) exists between the values predicted by AlogP and the experimentally determined logD values, thus enabling more accurate predictions.

Controlling Swing Rates in Macrocyclic Molecular Mortise Hinges.

Hinge motion is observed in macrocyclic, mortise-type molecular hinges using variable temperature NMR spectroscopy. The data is consistent with dynamic hinging from a folded-to-extended-to-folded enantiomeric state. Crystallographic and solution structures of the folded states are reported. Chemical shift predictions derived from crystallographic data corroborate fully revolute hinge motion. The rate of hinging is affected by steric congestion at the hinge axis. A macrocycle containing glycine, 1, hinges faster than one comprising aminoisobutyric acid, 2. The free energies of activation, ΔG≠ , for 1 and 2 were determined to be 13.3±0.3 kcal/mol and 16.3±0.3 kcal/mol, respectively. This barrier is largely independent of solvent across those surveyed (CD3 OD, CD3 CN, DMSO-d6 , pyridine-d5 , D2 O). Experiment and computation predict energy barriers that are consistent with disruption of an intramolecular network of hydrogen bonds. DFT calculations reveal a pathway for hinge motion.

Computational and Experimental Evidence for Templated Macrocyclization: The Role of a Hydrogen Bond Network in the Quantitative Dimerization of 24-Atom Macrocycles.

In the absence of preorganization, macrocyclization reactions are often plagued by oligomeric and polymeric side products. Here, a network of hydrogen bonds was identified as the basis for quantitative yields of macrocycles derived from the dimerization of monomers. Oligomers and polymers were not observed. Macrocyclization, the result of the formation of two hydrazones, was hypothesized to proceed in two steps. After condensation to yield the monohydrazone, a network of hydrogen bonds formed to preorganize the terminal acetal and hydrazine groups for cyclization. Experimental evidence for preorganization derived from macrocycles and acyclic models. Solution NMR spectroscopy and single-crystal X-ray diffraction revealed that the macrocycles isolated from the cyclization reaction were protonated twice. These protons contributed to an intramolecular network of hydrogen bonds that engaged distant carbonyl groups to realize a long-range order. DFT calculations showed that this network of hydrogen bonds contributed 8.7 kcal/mol to stability. Acyclic models recapitulated this network in solution. Condensation of an acetal and a triazinyl hydrazine, which adopted a number of conformational isomers, yielded a hydrazone that adopted a favored rotamer conformation in solution. The critical hydrogen-bonded proton was also evident. DFT calculations of acyclic models showed that the rotamers were isoenergetic when deprotonated. Upon protonation, however, energies diverged with one low-energy rotamer adopting the conformation observed in the macrocycle. This conformation anchored the network of hydrogen bonds of the intermediate. Computation revealed that the hydrogen-bonded network in the acyclic intermediate contributed up to 14 kcal/mol of stability and preorganized the acetal and hydrazine for cyclization.

A Model for the Rapid Assessment of Solution Structures for 24-Atom Macrocycles: The Impact of ß-Branched Amino Acids on Conformation.

Experiment and computation are used to develop a model to rapidly predict solution structures of macrocycles sharing the same Murcko framework. These 24-atom triazine macrocycles result from the quantitative dimerization of identical monomers presenting a hydrazine group and an acetal tethered to an amino acid linker. Monomers comprising glycine and the ß-branched amino acids threonine, valine, and isoleucine yield macrocycles G-G, T-T, V-V, and I-I, respectively. Elements common to all members of the framework include the efficiency of macrocyclization (quantitative), the solution- and solid-state structures (folded), the site of protonation (opposite the auxiliary dimethylamine group), the geometry of the hydrazone (E), the C2 symmetry of the subunits (conserved), and the rotamer state adopted. In aggregate, the data reveal metrics predictive of the three-dimensional solution structure that derive from the fingerprint region of the 1D 1H spectrum and a network of rOes from a single resonance. The metrics also afford delineation of more nuanced structural features that allow subpopulations to be identified among the members of the framework. Well-tempered metadynamics provides free energy surfaces and population distributions of these macrocycles. The areas of the free energy surface decrease with increasing steric bulk (G-G > V-V ∼ T-T > I-I). In addition, the surfaces are increasingly isoenergetic with decreasing steric bulk (G-G > V-V ∼ T-T > I-I).

Well-Tempered Metadynamics Simulations Predict the Structural and Dynamic Properties of a Chiral 24-Atom Macrocycle in Solution.

Inspired by therapeutic potential, the molecular engineering of macrocycles is garnering increased interest. Exercising control with design, however, is challenging due to the dynamic behavior that these molecules must demonstrate in order to be bioactive. Herein, the value of metadynamics simulations is demonstrated: the free-energy surfaces calculated reveal folded and flattened accessible conformations of a 24-atom macrocycle separated by barriers of ∼6 kT under experimentally relevant conditions. Simulations reveal that the dominant conformer is folded-an observation consistent with a solid-state structure determined by X-ray crystallography and a network of rOes established by 1H NMR. Simulations suggest that the macrocycle exists as a rapidly interconverting pair of enantiomeric, folded structures. Experimentally, 1H NMR shows a single species at room temperature. However, at lower temperature, the interconversion rate between these enantiomers becomes markedly slower, resulting in the decoalescence of enantiotopic methylene protons into diastereotopic, distinguishable resonances due to the persistence of conformational chirality. The emergence of conformational chirality provides critical experimental support for the simulations, revealing the dynamic nature of the scaffold-a trait deemed critical for oral bioactivity.