The molecular dynamics simulations of fentanyl complexed with the μ-opioid receptor (μOR) were studied using both inactive 4DKL and active 5C1M opioid receptor crystal structures. Analogous simulations in morphine with or without a ligand were done for comparison. Simulations of the inactive states were carried out in the absence and presence of the Na + ion. The obtained fentanyl's binding mode agrees with some of the mutagenesis data, and it overlaps with that of morphine only to a minor extent. Notably, fentanyl stabilizes different rotameric states of Trp293 6.48 than observed for morphine or unliganded receptor. Another difference is tighter arrangement of the interaction between Asp147 3.32 and Tyr326 7.43 (a link between helices TM3 and TM7) in the presence of fentanyl. Principal component analysis reveals differences in the trajectories dependent on the ligand bound. The differences found could be linked to ligand-dependent efficacy with respect to receptor intracellular signaling events.

The NICS(1) magnetic aromaticity index is split into NICS(1) and NICS(−1) indices when the points 1 Å above and below the ring center are inequivalent by symmetry. The two indices characterize the aromaticity of the two ring faces rather than the ring itself.

Interactions of 21 fentanyl derivatives with μ-opioid receptor (μOR) were studied using experimental and theoretical methods. Their binding to μOR was assessed with radioligand competitive binding assay. A uniform set of binding affinity data contains values for two novel and one previously uncharacterized derivative. The data confirms trends known so far and thanks to their uniformity, they facilitate further comparisons. In order to provide structural hypotheses explaining the experimental affinities, the complexes of the studied derivatives with μOR were modeled and subject to molecular dynamics simulations. Five common General Features (GFs) of fentanyls’ binding modes stemmed from these simulations. They include: GF1) the ionic interaction between D147 and the ligands’ piperidine NH+ moiety; GF2) the N-chain orientation towards the μOR interior; GF3) the other pole of ligands is directed towards the receptor outlet; GF4) the aromatic anilide ring penetrates the subpocket formed by TM3, TM4, ECL1 and ECL2; GF5) the 4-axial substituent (if present) is directed towards W318. Except for the ionic interaction with D147, the majority of fentanyl-μOR contacts is hydrophobic. Interestingly, it was possible to find nonlinear relationships between the binding affinity and the volume of the N-chain and/or anilide’s aromatic ring. This kind of relationships is consistent with the apolar character of interactions involved in ligand–receptor binding. The affinity reaches the optimum for medium size while it decreases for both large and small substituents. Additionally, a linear correlation between the volumes and the average dihedral angles of W293 and W133 was revealed by the molecular dynamics study. This seems particularly important, as the W293 residue is involved in the activation processes. Further, the Y326 (OH) and D147 (Cγ) distance found in the simulations also depends on the ligands’ size. In contrast, neither RMSF measures nor D114/Y336 hydrations show significant structure-based correlations. They also do not differentiate studied fentanyl derivatives. Eventually, none of 14 popular scoring functions yielded a significant correlation between the predicted and observed affinity data (R < 0.30, n = 28).

sEDA, pEDA, and cSAR descriptors of the substituent effect were determined for >30 monosubstituted benzenes in the first excited singlet S state at the LC-ωB97XD/aug-cc-pVTZ level. It was found that in the S state, the σ- and π-valence electrons are a bit less and a bit more affected, respectively, than in the S state, but basically, the effect in both states remains the same. In the S and S states, the d(C-X) distances to the substituent's first atom and the ring perimeter correlate with the sEDA and pEDA in the appropriate states, respectively. The energies and the gap of the frontier orbitals in the two states are linearly correlated and for the HOMO(S), LUMO(S), and HOMO(S)-LUMO(S) gap correlate also with the pEDA(S) and cSAR(S) descriptors. In all studied correlations, three similar groups of substituents can be distinguished, for which correlations (i) are very good, (ii) deviate slightly, and (iii) deviate significantly. Comparison of the shape of the HOMO(S) and HOMO(S) orbitals shows that for case (i) HOMO orbitals exhibit almost perfect antisymmetry against the benzene plane, for case (ii) the antisymmetry of HOMO in one of the states is either perturbed or changed, and for case (iii) one HOMO state has σ-character.

No abstract

To measure molecular chirality, the molecule is treated as a finite set of points in the Euclidean R(3) space supplemented by k properties, p(1)((i)), p(2)((i)), ..., p(k)((i)) assigned to the ith atom, which constitute a point in the Property P(k) space. Chirality measures are described as the distance between a molecule and its mirror image minimized over all its arbitrary orientation-preserving isometries in the R(3) × P(k) Cartesian product space. Following this formalism, different chirality measures can be estimated by taking into consideration different sets of atomic properties. Here, for α-amino acid zwitterionic structures taken from the Cambridge Structural Database and for all 1684 neutral conformers of 19 biogenic α-amino acid molecules, except glycine and cystine, found at the B3LYP/6-31G** level, chirality measures have been calculated by a CHIMEA program written in this project. It is demonstrated that there is a significant correlation between the measures determined for the α-amino acid zwitterions in crystals and the neutral forms in the gas phase. Performance of the studied chirality measures with changes of the basis set and computation method was also checked. An exemplary quantitative structure–activity relationship (QSAR) application of the chirality measures was presented by an introductory model for the benchmark Cramer data set of steroidal ligands of the sex-hormone binding globulin.

The area of multitarget compounds, joining two pharmacophores within one molecule, is a vivid field of research in medicinal chemistry. Not only pharmacophoric elements are essential for the design and activity of such compounds, but the type and length of linkers used to connect them are also crucial. In the present contribution, we describe compound 1 in which a typical opioid peptide sequence is combined with a fragment characteristic for neurokinin-1 receptor (NK1R) antagonists through a hydrazone bridge. The compound has a high affinity for μ-and δ-opioid receptors (IC 50 = 12.7 and 74.0 nM, respectively) and a weak affinity for the NK1R. Molecular modeling and structural considerations explain the observed activities. In in vivo test, intrathecal and intravenous administrations of 1 exhibited a strong analgesic effect, which indicates potential BBB penetration. This letter brings an exemplary application of the hydrazone linker for fast, facile, and successful preparation of chimeric compounds.

No abstract