Structure-Kinetics Relationships of Opioids from Metadynamics and Machine Learning Analysis.

The nation's opioid overdose deaths reached an all-time high in 2021. The majority of deaths are due to synthetic opioids represented by fentanyl. Naloxone, which is a FDA-approved reversal agent, antagonizes opioids through competitive binding at the µ-opioid receptor (mOR). Thus, knowledge of the opioid's residence time is important for assessing the effectiveness of naloxone. Here, we estimated the residence times (τ) of 15 fentanyl and 4 morphine analogs using metadynamics and compared them with the most recent measurement of the opioid kinetic, dissociation, and naloxone inhibitory constants (Mann et al. Clin. Pharmacol. Therapeut. 2022, 120, 1020-1232). Importantly, the microscopic simulations offered a glimpse at the common binding mechanism and molecular determinants of dissociation kinetics for fentanyl analogs. The insights inspired us to develop a machine learning approach to analyze the kinetic impact of fentanyl's substituents based on the interactions with mOR residues. This proof-of-concept approach is general; for example, it may be used to tune ligand residence times in computer-aided drug discovery.

Structure-Kinetics Relationships of Opioids from Metadynamics and Machine Learning.

The nation's opioid overdose deaths reached an all-time high in 2021. The majority of deaths are due to synthetic opioids represented by fentanyl. Naloxone, which is an FDA-approved reversal agent, antagonizes opioids through competitive binding at the mu-opioid receptor (mOR). Thus, knowledge of opioid's residence time is important for assessing the effectiveness of naloxone. Here we estimated the residence times of 15 fentanyl and 4 morphine analogs using metadynamics, and compared them with the most recent measurement of the opioid kinetic, dissociation, and naloxone inhibitory constants (Mann, Li et al, Clin. Pharmacol. Therapeut. 2022). Importantly, the microscopic simulations offered a glimpse at the common binding mechanism and molecular determinants of dissociation kinetics for fentanyl analogs. The insights inspired us to develop a machine learning (ML) approach to analyze the kinetic impact of fentanyl's substituents based on the interactions with mOR residues. This proof-of-concept approach is general; for example, it may be used to tune ligand residence times in computer-aided drug discovery.

How µ-opioid receptor recognizes fentanyl.

Roughly half of the drug overdose-related deaths in the United States are related to synthetic opioids represented by fentanyl which is a potent agonist of mu-opioid receptor (mOR). In recent years, X-ray crystal structures of mOR in complex with morphine derivatives have been determined; however, structural basis of mOR activation by fentanyl-like opioids remains lacking. Exploiting the X-ray structure of BU72-bound mOR and several molecular simulation techniques, we elucidated the detailed binding mechanism of fentanyl. Surprisingly, in addition to the salt-bridge binding mode common to morphinan opiates, fentanyl can move deeper and form a stable hydrogen bond with the conserved His2976.52, which has been suggested to modulate mOR's ligand affinity and pH dependence by previous mutagenesis experiments. Intriguingly, this secondary binding mode is only accessible when His2976.52 adopts a neutral HID tautomer. Alternative binding modes may represent a general mechanism in G protein-coupled receptor-ligand recognition.

Kinetics and Mechanism of Fentanyl Dissociation from the µ-Opioid Receptor.

Driven by illicit fentanyl, opioid related deaths have reached the highest level in 2020. Currently, an opioid overdose is resuscitated by the use of naloxone, which competitively binds and antagonizes the µ-opioid receptor (mOR). Thus, knowledge of the residence times of opioids at mOR and the unbinding mechanisms is valuable for assessing the effectiveness of naloxone. In the present study, we calculate the fentanyl-mOR dissociation time and elucidate the mechanism by applying an enhanced sampling molecular dynamics (MD) technique. Two sets of metadynamics simulations with different initial structures were performed while accounting for the protonation state of the conserved H2976.52, which has been suggested to modulate the ligand-mOR affinity and binding mode. Surprisingly, with the Nδ-protonated H2976.52, fentanyl can descend as much as 10 Å below the level of the conserved D1473.32 before escaping the receptor and has a calculated residence time τ of 38 s. In contrast, with the Nϵ- and doubly protonated H2976.52, the calculated τ are 2.6 and 0.9 s, respectively. Analysis suggests that formation of the piperidine-Hid297 hydrogen bond strengthens the hydrophobic contacts with the transmembrane helix (TM) 6, allowing fentanyl to explore a deep pocket. Considering the experimental τ of ∼4 min for fentanyl and the role of TM6 in mOR activation, the deep insertion mechanism may be biologically relevant. The work paves the way for large-scale computational predictions of opioid dissociation rates to inform evaluation of strategies for opioid overdose reversal. The profound role of the histidine protonation state found here may shift the paradigm in computational studies of ligand-receptor kinetics.

How mu-Opioid Receptor Recognizes Fentanyl.

The opioid crisis has escalated during the COVID-19 pandemic. More than half of the overdose-related deaths are related to synthetic opioids represented by fentanyl which is a potent agonist of mu-opioid receptor (mOR). In recent years, crystal structures of mOR complexed with morphine derivatives have been determined; however, structural basis of mOR activation by fentanyl-like synthetic opioids remains lacking. Exploiting the X-ray structure of mOR bound to a morphinan ligand and several state-of-the-art simulation techniques, including weighted ensemble and continuous constant pH molecular dynamics, we elucidated the detailed binding mechanism of fentanyl with mOR. Surprisingly, in addition to the orthosteric site common to morphinan opiates, fentanyl can move deeper and bind mOR through hydrogen bonding with a conserved histidine H297, which has been shown to modulate mOR's ligand affinity and pH dependence in mutagenesis experiments, but its precise role remains unclear. Intriguingly, the secondary binding mode is only accessible when H297 adopts a neutral HID tautomer. Alternative binding modes and involvement of tautomer states may represent general mechanisms in G protein-coupled receptor (GPCR)-ligand recognition. Our work provides a starting point for understanding mOR activation by fentanyl analogs that are emerging at a rapid pace and assisting the design of safer analgesics to combat the opioid crisis. Current protein simulation studies employ standard protonation and tautomer states; our work demonstrates the need to move beyond the practice to advance our understanding of protein-ligand recognition.

How µ-Opioid Receptor Recognizes Fentanyl.

In 2019, drug overdose has claimed over 70,000 lives in the United States. More than half of the deaths are related to synthetic opioids represented by fentanyl which is a potent agonist of mu-opioid receptor (mOR). In recent years, the crystal structures of mOR in complex with morphine derivatives have been determined; however, structural basis of mOR activation by fentanyl-like synthetic opioids remains lacking. Exploiting the X-ray structure of mOR bound to a morphinan ligand and several state-of-the-art simulation techniques, including weighted ensemble and continuous constant pH molecular dynamics, we elucidated the detailed binding mechanism of fentanyl with mOR. Surprisingly, in addition to forming a salt-bridge with Asp1473.32 in the orthosteric site common to morphinan opiates, fentanyl can move deeper and bind mOR through hydrogen bonding with a conserved histidine His2976.52, which has been shown to modulate mOR's ligand affinity and pH dependence in mutagenesis experiments, but its precise role remains unclear. Intriguingly, the secondary binding mode is only accessible when His297 adopts a neutral HID tautomer. Alternative binding modes and involvement of tautomer states may represent general mechanisms in G protein-coupled receptor (GPCR)-ligand recognition. Our work provides a starting point for understanding the molecular basis of mOR activation by fentanyl which has many analogs emerging at a rapid pace. The knowledge may also inform the design of safer analgesics to combat the opioid crisis. Current protein simulation studies employ standard protonation and tautomer states; our work demonstrates the need to move beyond the practice to advance our understanding of protein-ligand recognition.

Microscopic Behaviors of Tri- n-Butyl Phosphate, n-Dodecane, and Their Mixtures at Air/Liquid and Liquid/Liquid Interfaces: An AMBER Polarizable Force Field Study.

In solvent extraction processes for recovering metal ions from used nuclear fuel, as well as other industrial applications, a better understanding of the metal complex phase transfer phenomenon would greatly aid ligand design and process optimization. We have approached this challenge by utilizing the classical molecular dynamics simulations technique to gain visual appreciation of the vapor/liquid and liquid/liquid interface between tri- n-butyl phosphate (TBP) and n-dodecane with air and water. In this study, we successfully reparameterized polarizable force fields for TBP and n-dodecane that accurately reproduced several of their thermophysical properties such as density, heat of vaporization, and dipole moment. Our models were able to predict the surface and interfacial tension of different systems when compared to experimental results that were also performed by us. Through this study, we gained atomistic understanding of the behaviors of TBP and n-dodecane at the interface against air and water, useful in further computational studies of such systems. Finally, our studies indicate that the initial configuration of a simulation may have a large effect on the final result.

Quantifying Dimer and Trimer Formation by Tri-n-butyl Phosphates in n-Dodecane: Molecular Dynamics Simulations.

Tri-n-butyl phosphate (TBP), a representative of neutral organophosphorous ligands, is an important extractant used in the solvent extraction process for the recovery of uranium and plutonium from spent nuclear fuel. Microscopic pictures of TBP isomerism and its behavior in n-dodecane diluent were investigated utilizing MD simulations with previously optimized force field parameters for TBP and n-dodecane. Potential mean force (PMF) calculations on a single TBP molecule show seven probable TBP isomers. Radial distribution functions (RDFs) of TBP suggest the existence of TBP trimers at high TBP concentrations in addition to dimers. 2D PMF calculations were performed to determine the angle and distance criteria for TBP trimers. The dimerization and trimerization constants of TBP in n-dodecane were obtained and match our own experimental values using the FTIR technique. The new insights into the conformational behaviors of the TBP molecule as a monomer and as part of an aggregate could greatly aid in the understanding of the complexation between TBP and metal ions in a solvent extraction system.

Quantifying Dimer and Trimer Formation of Tri-n-butyl Phosphates in Different Alkane Diluents: FTIR Study.

Tri-n-butyl phosphate (TBP), a representative of neutral organophosphorous metal-ion-extracting reagents, is an important ligand used in solvent extraction processes for the recovery of uranium and plutonium from spent nuclear fuel, as well as other non-nuclear applications. Ligand-ligand and organic solvent-ligand interactions play an important role in these processes. The self-association behavior of TBP in various alkane diluents of different chain lengths (8, 12, and 16 carbons) and a branched alkane (iso-octane) was investigated by Fourier transform infrared spectroscopic measurements. By careful deconvolution of the spectra into multiple peaks, our results indicate that TBP self-associates to form not only dimers, as previous studies showed, but also trimers in the practical concentration range. Using a mathematical fitting procedure, the dimerization and trimerization constants were determined. As expected, these equilibrium constants are dependent on the solvent used. As the alkane chain for linear hydrocarbon solvents becomes longer, dimerization decreases whereas trimerization increases. For the more branched hydrocarbon, we observe a significantly higher dimerization constant. These effects are most likely due to the intermolecular van der Waals interactions between the butyl tails of each TBP molecule and the diluent hydrocarbon chain as all solvents in this study are relatively nonpolar.

Computational study of molecular structure and self-association of tri-n-butyl phosphates in n-dodecane.

Tri-n-butyl phosphate (TBP) is an important extractant used in the solvent extraction process for recovering uranium and plutonium from used nuclear fuel. An atomistic molecular dynamics study was used to understand the fundamental molecular-level behavior of extracting agents in solution. Atomistic parametrization was carried out using the AMBER force field to model the TBP molecule and n-dodecane molecule, a commonly used organic solvent. Validation of the optimized force field was accomplished through various thermophysical properties of pure TBP and pure n-dodecane in the bulk liquid phase. The mass density, dipole moment, self-diffusion coefficient, and heat of vaporization were calculated from our simulations and compared favorably with experimental values. The molecular structure of TBPs in n-dodecane at a dilute TBP concentration was examined based on radial distribution functions. 1D and 2D potential mean force studies were carried out to establish the criteria for identifying TBP aggregates. The dimerization constant of TBP in the TBP/n-dodecane mixture was also obtained and matched the experimental value.

The ATM/p53 pathway is commonly targeted for inactivation in squamous cell carcinoma of the head and neck (SCCHN) by multiple molecular mechanisms.

The ATM/p53 pathway plays a critical role in maintenance of genome integrity and can be targeted for inactivation by a number of characterized mechanisms including somatic genetic/epigenetic alterations and expression of oncogenic viral proteins. Here, we examine a panel of 24 SCCHN tumors using various molecular approaches for the presence of human papillomavirus (HPV), mutations in the p53 gene and methylation of the ATM promoter. We observed that 30% of our SCCHN samples displayed the presence of HPV and all but one was HPV type 16. All HPV E6 gene-positive tumors exhibited E6 transcript expression. We observed 21% of the tumors harbored p53 mutations and 42% of tumors displayed ATM promoter methylation. The majority of tumors (71%) were positive for at least one of these events. These findings indicate that molecular events resulting in inactivation of the ATM/p53 pathway are common in SCCHN and can arise by a number of distinct mechanisms.

The ATM gene is a target for epigenetic silencing in locally advanced breast cancer.

Several epidemiological studies on ataxia-telangiectasia families indicate that obligate ATM heterozygotes display an elevated risk for developing breast cancer. However, a molecular basis for a potential link between diminished ATM function and sporadic breast malignancy remains elusive. Here, we show that 78% (18 out of a panel of 23) of surgically removed breast tumors (stage II or greater) displayed aberrant methylation of the ATM proximal promoter region as judged by methylation-specific PCR. Aberrant methylation of the ATM promoter was independently confirmed in several tumors by bisulfite sequencing. Moreover, bisulfite sequencing indicated that this region of the genome is subject to dense methylation. Further, we found a highly significant correlation (P = 0.0006) between reduced ATM mRNA abundance, as measured by real-time RT-PCR, and aberrant methylation of the ATM gene promoter. These findings indicate that epigenetic silencing of ATM expression occurs in locally advanced breast tumors, and establish a link at the molecular level between reduced ATM function and sporadic breast malignancy.

Ataxia-telangiectasia-mutated (ATM) gene in head and neck squamous cell carcinoma: promoter hypermethylation with clinical correlation in 100 cases.

The Ataxia-telangiectasia-mutated (ATM) gene product is a well-characterized tumor suppressor that plays a key role in maintenance of genomic stability. We have recently documented that the ATM promoter is a target for epigenetic silencing in cultured tumor cells. Here we show that aberrant methylation of the ATM promoter occurs in a significant percentage (25%) of head and neck squamous cell carcinomas. The presence of methylated ATM promoter shows a statistically significant correlation with an earlier age of initial diagnosis and decreased overall survival, particularly in early-stage tumors. These findings indicate that ATM promoter hypermethylation occurs in head and neck squamous cell carcinoma, and this feature is a potentially useful prognostic marker in this tumor type.

Aberrant methylation of the ATM promoter correlates with increased radiosensitivity in a human colorectal tumor cell line.

Recent findings suggest that DNA alkylating agents trigger cellular responses that overlap those activated after ionizing radiation. Moreover, activation of these responses is dependent upon a functional mismatch repair (MMR) system. These developments led us to test if MMR-deficient cells may be compromised in their ability to activate appropriate cellular signaling pathways after ionizing radiation. An initial experiment to address this notion was to determine the level of radiosensitivity of several MMR-deficient cell lines derived from patients with Hereditary Non-Polyposis Colorectal Cancer (HNPCC). While two of the three HNPCC lines investigated show levels of radiosensitivity consistent with that displayed by normal human fibroblasts, HCT-116 cells display moderate radiosensitivity compared to the other MMR-deficient lines. This increased sensitivity to ionizing radiation correlates with lowered levels of ATM expression in HCT-116. Analysis of genomic DNA from HCT-116 cells determined that these cells possess aberrant methylation of multiple CpG dinucleotides within the proximal promoter region of the ATM gene. The significance of this finding is underscored by our observations that co-culturing HCT-116 cells with the DNA demethylating agent 5-azacytidine reverses promoter methylation, promotes normal levels of ATM expression, and restores normal radiosensitivity. The proximal ATM promoter is a approximately 520 bp region shared with the NPAT gene, and current evidence suggests that this region functions as a bi-directional promoter. We found that, unlike ATM, the methylation status of this intergenic region does not effect the expression of the NPAT gene. In sum, these observations indicate that the ATM gene is a novel target for epigentic silencing through inappropriate methylation of its proximal promoter region.