Tricarbonyl rhenium(i) complexes with 8-hydroxyquinolines: structural, chemical, antibacterial, and anticancer characteristics.

Twelve tricarbonyl rhenium(i) complexes in the '2 + 1' system with the anionic bidentate N,O-donor ligand (deprotonated 8-hydroxyquinoline (HQ) or its 2-methyl (MeHQ) or 5-chloro (ClHQ) derivative) and neutral N-donor diazoles (imidazole (Him), 2-methylimidazole (MeHim), 3,5-dimethylpyrazole (Hdmpz), and 3-phenylpyrazole (HPhpz)) were synthesized: [Re(CO)3(LN,O)LN] (LN,O = Q-, MeQ-, ClQ-; LN = Him, MeHim, Hdmpz, HPhpz). Their crystal structures were determined by the scXRD method, compared with the DFT-calculated ones, and characterized by analytical (EA) and spectroscopic techniques (FT-IR, NMR, and UV-Vis) interpreted with DFT and TD-DFT calculations. Most of the Re(i) complexes did not show relevant antibacterial activity against Gram-negative and Gram-positive bacterial strains. Only [Re(CO)3(MeQ)Him] demonstrated significant action 4-fold better against Gram-negative Pseudomonas aeruginosa than the free MeHQ ligand. The cytotoxicity of the compounds was estimated using human acute promyelocytic leukemia (HL-60), ovarian (SKOV-3), prostate (PC-3), and breast (MCF-7) cancer, and breast non-cancerous (MCF-10A) cell lines. Only HQ and ClHQ ligands and [Re(CO)3(Q)Hdmpz] complex had good selectivity toward MCF-7 cell line. HL-60 cells were sensitive to all complexes (IC50 = 1.5-14 µM). Still, pure HQ and ClHQ ligands were slightly more active than the complexes.

VCD spectra of chiral naphthalene-1-carboxamides in the solid-state.

The VCD spectra of chiral 2,3-dihydro-1H-benzo[de]isoquinolin-1-one (8-substituted naphthalene-1-carboxamide, BIQ) were studied in KBr pellets. The X-ray diffractometry revealed that the Me, Ph, and pClPh BIQs crystalize in the monoclinic P21, while nBu, pMePh, and oMeOPh BIQs in the orthorhombic P212121 space group. Only the Me-BIQ crystal exhibits the presence of cyclic amide dimers, while the others contain chains of the amid group hydrogen bonds. For all BIQs, except pMePh, the most intense IR band in the 1750-1550 cm-1 region is located at ca. 1680 cm-1 and is accompanied by two weak ones at ca. 1618 and 1590 cm-1. For the pMePh derivative, four almost equally intense IR bands at 1662, 1639, 1614, and 1588 cm-1 are observed. This region of the IR spectra of BIQs, but pMePh, is well reproduced by calculations based on BIQ monomers. On the other hand, the complex IR pattern of pMePh is computationally reproduced when larger crystal fragments, like octamers, are considered. Registration of the VCD spectra enabled recognizing the complexity of IR contours at ca. 1680 cm-1 by the corresponding VCD motives. For (i) Me, Ph and pClPh (R)-enantiomers, two (+)(-) bands were distinguished and for (ii) nBu and pMePh ones, one VCD band with right-side asymmetry was found. For (iii) oMeOPh the VCD pattern cannot be unambiguously assigned. Thus, the VCD spectra in the ν(C=O) range diverse the studied compounds. Among the set of molecules, pMePh has exceptional crystal geometry. Therefore, its most intense ν(C=O) band position and shape can be connected with the geometry of the hydrogen bonds, interactions, and crystal packing. Interpretation of the VCD spectra is based on linear and packed BIQ octamers. This cluster model can reproduce the main features of the solid-state VCD of BIQs.

HOMA Index Establishes Similarity to a Reference Molecule.

The article shows that the definition of the HOMA index of geometrical aromaticity satisfies the axioms of a similarity function between the examined and benzene ring. Consequently, for purely mathematical reasons, the index works exceptionally well as an index of aromaticity: it expresses a geometric similarity to the archetypal aromatic benzene. Thus, if the molecule is geometrically similar to benzene, then it is also chemically similar, and therefore, it is aromatic. However, the similarity property legitimizes using the HOMA-like indices to express similarity to molecules other than benzene, whether cyclic or linear and existing or hypothetical. The paper demonstrates an example of HOMA-similarity to cyclohexane, which expresses a (relaxed)-saturicity property not accompanied by strong structural strains or steric hindrances. Further, it is also shown that the HOMA index can evaluate the properties of whole molecules, such as 25 unbranched catacondensed isomers of hexacene. The index exhibits a significant quadratic correlation with the total energy differences of planar isomers from which the nonplanar ones deviate. Moreover, the HOMA index of hexacene isomers significantly correlates with the Kekulé count connected to the resonance energy in the Hückel approximation. As a result, the study shows that the HOMA index can be used not only for aromaticity analyses but also as a general chemical descriptor applicable to rings, chains, composed molecular moieties, or even whole molecules.

On Integral INICS Aromaticity of Pyridodiazepine Constitutional Isomers and Tautomers.

The structure, energetics, and aromaticity of c.a. 100 constitutional isomers and tautomers of pyrido[m,n]diazepines (m = 1, 2; n = 2, 3, 4, 5; m ≠ n) were studied at the B3LYP/cc-pVTZ level. The pyrido[1,3]diazepines appear the most, while pyrido[2,4]diazepines are the least stable (ca. 26 kcal/mol). In the pyrido[1,n]diazepine group (n = 2-5), the [1,5] isomers are higher in energy by ca. 4.5 kcal/mol and the [1,4] ones by ca. 7 kcal/mol, and the pyrido[1,2]diazepines are the least stable (ca. 20 kcal/mol). All the most stable pyrido[1,n]diazepines have N-atoms near the ring's junction bond but on opposite sites. The most stable [2,n]-forms are also those with the pyridine ring N6-atom near the junction bond. Surprisingly, for the [1,2]-, [1,3]-, and [1,4]-isomer condensation types of pyridine and diazepine rings, the same N9 > N7 > N6 > N8 stability pattern obeys. The stability remains similar in a water medium simulated with the Polarizable Continuum Model of the solvent and is conserved when calculated using the CAM-B3LYP or BHandHlyp functionals. The ring's aromaticity in the pyridine[m,n]diazepines was established based on the integral INICS index resulting from the NICSzz-scan curves' integration. The integral INICS index is physically justified through its relation to the ringcurrent as demonstrated by Berger, R.J.F., et al. Phys. Chem. Chem. Phys. 2022, 24, 624. The six-membered pyrido rings have negative INICSZZ indices and can be aromatic only if they are not protonated at the N-atom. All protonated pyrido and seven-membered rings exhibit meaningful positive INICSZZ values and can be assigned as antiaromatic. However, some non-protonated pyrido rings also have substantial positive INICSZZ indices and are antiaromatic. A weak linear correlation (R2 = 0.72) between the INICSZZ values of the pyridine I(6) and diazepine I(7) rings exists and is a consequence of the communication between the π-electron systems of the two rings. The juxtaposition of the INICS descriptor of the six- and seven-membered rings and diverse electron density parameters at the Ring Critical Points (RCP) revealed good correlations only with the Electrostatic Potentials from the electrons and nuclei (ESPe and ESPn). The relationships with other RCP parameters like electron density and its Laplacian, total energy, and the Hamiltonian form of kinetic energy density were split into two parts: one nearly constant for the six-membered rings and one linearly correlating for the seven-membered rings. Thus, most of the electron density parameters at the RCP of the six-membered rings of pyridodiazepines practically do not change with the diazepine type and the labile proton position. In contrast, those of the seven-membered rings display aromaticity changes in the antiaromatic diazepine with its ring structural modifications.

UV-vis and ECD spectroelectrochemistry of atropisomeric naphthalenediimide derivative.

The UV-vis and ECD spectroelectrochemistry (SEC) of a chiral binaphthalenylamine derivative of the N-butyl naphthalenediimide (NDIB-NH2) enantiomers were applied to measure UV-vis and ECD spectra of NDIB-NH2 radicals and dianion formed in the reduction and oxidation processes observed in cyclic voltammetry (CV). The CV curves and EPR spectroelectrochemistry enabled us to establish conditions at which a radical-anion [NDIB-NH2]̇.-, a dianion [NDIB-NH2]2-, and a radical-cation [NDIB-NH2]̇.+ are formed. The DFT restricted open-shell CAM-B3LYP-D3/def2TZVP/PCM calculations demonstrated that in the radical-anion [NDIB-NH2]̇.-, spin is spread over the NDI system while in the radical-cation [NDIB-NH2]̇+ it is spread over the aminonaphthalene moiety. The UV-vis spectra of radical-anion and dianion show the most significant changes in the 400-800 nm range. In that range, the ECD spectra varied with the change of electrode potential more than the UV-vis did and enabled the identification of a new ECD band of [NDIB-NH2]̇.- at ca. 400 nm hidden in the background in the UV spectra at -1000 mV. A broad structured ECD pattern with a maximum at ca. 530 nm was observed for [NDIB-NH2]̇.- (-1000 mV), while a single smooth ECD band of [NDIB-NH2]2- was located at 520 nm (-1750 mV). For the first time, an isosbestic point (455 nm) was found in ECD spectroelectrochemical measurements for the radical-cation [NDIB-NH2]̇.+ in equilibrium with the NDIB-NH2 neutral form. The TD-DFT CAM-B3LYP-D3/6-31G** calculations combined with the hybrid (explicit combined with implicit) solvation model fairly well reproduced the UV-vis and ECD SEC of neutral and redox forms of NDIB-NH2 but the ECD spectrum of [NDIB-NH2]̇.+ above 390 nm.

On Aromaticity of the Aromatic α-Amino Acids and Tuning of the NICS Indices to Find the Aromaticity Order.

The NICS aromaticity indices of the rings in flexible phenylalanine (Phe), tryptophan (Trp), tyrosine (Tyr), and histidine (His) chiral molecules were analyzed. These molecules have several dozens of conformers, and their rings are slightly non-planar. Therefore, the population-averaged NICSpav index was defined, and the NICS scans had to be performed with respect to planes found by the least-squares routine. A rule differentiating an obverse and a reverse ring face in aromatic amino acids was formulated. The NICS scan minima corresponding to the obverse and reverse face were unequal, which prompted us to use the term ring face aromaticity/ring face tropicity. It appeared that for Phe, Trp, Tyr, and His, the reverse face has always had higher ring face aromaticity/ring face tropicity than the obverse one. Despite the NICS modifications, uncertainty about the amino acid aromaticity order remained. This motivated us to use the integral INICS index newly proposed by Stanger as well. Then, the following sequence was obtained: Trp(phenyl) > Phe > Trp(pyrrole) > His > Tyr. The juxtaposition of the INICS indices of amino acids with that of some model rings revealed a fair transferability of the values. Finally, analysis of the substituent effect on INICS demonstrated that the aromaticity of Tyr is the lowest due to the strength of the OH group π-electron-donating effect able to perturb enough the ring charge distribution and its magnetic aromaticity. The NICS calculations were executed using the ARONICS program written within the project.

New chiral ECD-Raman spectroscopy of atropisomeric naphthalenediimides.

In this study, we found that a recently discovered ECD-Raman effect dominated over the natural Raman optical activity in a series of atropisomeric naphthalenediimides, and we investigated the kind of information about the molecular structure that could be obtained from the spectra. The ECD-Raman effect is polarised Raman scattering modulated by electronic circular dichroism. We showed that the spectra significantly depended on the substitution of the solute and/or the change of the solvent. Moreover, the spectra could be well-predicted by the theory, thus providing an interesting tool to monitor the chirality of the binaphthyl compounds.

The solid state VCD of a novel N-acylhydrazone trifluoroacetate.

A novel N-acylhydrazone with pharmaceutical importance was subject of structural and IR/VCD investigations in the solid state. In the crystal structure, dimers of anion-cation pairs are stabilized by H-bonding and ionic interactions. Some less common interaction types, like C=N···C-NH3+ (σ-hole) interactions, hydrazone-aromatic interactions and dispersive contacts of the CF3 groups are also present in the crystal. Satisfactory reproduction of the solid state IR and VCD spectra required that quantum-chemical calculations be done on a tetramer (four cation-anion pairs) cut out from the crystal structure, exhibiting key intermolecular interactions. Ten DFT functionals were assessed as to the agreement between the calculated and experimental spectra. Various approaches to scaling of the calculated frequencies were applied. The best results were yielded with individual (optimized) frequency scaling factors (FSFs) and band half-widths at half maximum-(HWHM) for four separate spectral subregions. The best matching between the experimental and theoretical spectra (according to SimIR, SimVCD and SimVDF indices) was found for the B3PW91 functional, however, a few other functionals follow closely in the ranking. Based on the quantum chemical calculations, spectral assignments have been made.

Substituent Effect in the Cation Radicals of Monosubstituted Benzenes.

In 30 monosubstituted benzene cation radicals, studied at the ωB97XD/aug-cc-pVTZ level, the phenyl rings usually adopt a compressed form, but a differently compressed form-equivalent to an elongated one-may coexist. The computational and literature ionization potentials are well correlated. The geometrical and magnetic aromaticity, estimated using HOMA and NICS indices, show the systems to be structurally aromatic but magnetically antiaromatic or only weakly aromatic. The partial charge is split between the substituent and ring and varies the most at C(ipso). In the ring, the spin is 70%, concentrated equally at the C(ipso) and C(p) atoms. The sEDA(D) and pEDA(D) descriptors of the substituent effect in cation radicals, respectively, were determined. In cation radicals, the substituent effect on the σ-electron system is like that in the ground state. The effect on the π-electron systems is long-range, and its propagation in the radical quinone-like ring is unlike that in the neutral molecules. The pEDA(D) descriptor correlates well with the partial spin at C(ipso) and C(p) and weakly with the HOMA(D) index. The correlation of the spin at the ring π-electron system and the pEDA(D) descriptor shows that the electron charge supplied to the ring π-electron system and the spin flow oppositely.

ECD spectroelectrochemistry: A review.

The electronic circular dichroism (ECD) spectroscopy is probably the most important chiraloptical method, and the role of chirality in contemporary chemistry, pharmacy, and material science constantly increases. On the other hand, the electrochemical methods are also very sensitive tools for studying multivarious redox processes. Nevertheless, the first ECD spectroelectrochemical (SEC) study was only published by Daub, Salbeck and Aurbach in 1988, and since then, the ECD SEC method has been mentioned in only thirty papers. By the summer of 2020, the ECD SEC studies were mainly focused around molecular systems for organic, and marginally, inorganic chiroptical switching studies of biochemical redox reactions. The review provides more details about the ECD SEC studies carried out so far. At the end, we suggest some future applications for the ECD spectroelectrochemistry.

Spectral Probe for Electron Transfer and Addition Reactions of Azide Radicals with Substituted Quinoxalin-2-Ones in Aqueous Solutions.

The azide radical (N3●) is one of the most important one-electron oxidants used extensively in radiation chemistry studies involving molecules of biological significance. Generally, it was assumed that N3● reacts in aqueous solutions only by electron transfer. However, there were several reports indicating the possibility of N3● addition in aqueous solutions to organic compounds containing double bonds. The main purpose of this study was to find an experimental approach that allows a clear assignment of the nature of obtained products either to its one-electron oxidation or its addition products. Radiolysis of water provides a convenient source of one-electron oxidizing radicals characterized by a very broad range of reduction potentials. Two inorganic radicals (SO4●-, CO3●-) and Tl2+ ions with the reduction potentials higher, and one radical (SCN)2●- with the reduction potential slightly lower than the reduction potential of N3● were selected as dominant electron-acceptors. Transient absorption spectra formed in their reactions with a series of quinoxalin-2-one derivatives were confronted with absorption spectra formed from reactions of N3● with the same series of compounds. Cases, in which the absorption spectra formed in reactions involving N3● differ from the absorption spectra formed in the reactions involving other one-electron oxidants, strongly indicate that N3● is involved in the other reaction channel such as addition to double bonds. Moreover, it was shown that high-rate constants of reactions of N3● with quinoxalin-2-ones do not ultimately prove that they are electron transfer reactions. The optimized structures of the radical cations (7-R-3-MeQ)●+, radicals (7-R-3-MeQ)● and N3● adducts at the C2 carbon atom in pyrazine moiety and their absorption spectra are reasonably well reproduced by density functional theory quantum mechanics calculations employing the ωB97XD functional combined with the Dunning's aug-cc-pVTZ correlation-consistent polarized basis sets augmented with diffuse functions.

Thallium(I) Tropolonates: Synthesis, Structure, Spectral Characteristics, and Antimicrobial Activity Compared to Lead(II) and Bismuth(III) Analogues.

Synthesis, single-crystal X-ray determination diffraction and FT-IR, NMR (1H, 13C, 19F and 205Tl), UV-vis, and luminescence spectra characteristics were described for series of thallium(I) compounds: thallium(I) triflate (Tl(OTf)), 1:1 co-crystals of thallium(I) triflate and tropolone (Htrop), Tl(OTf)·Htrop, as well as simple thallium(I) chelates: Tl(trop) (1), Tl(5-metrop) (2), Tl(hino) (3), with Htrop, 5-methyltropolone (5-meHtrop), 4-isopropyltropolone (hinokitiol, Hhino), respectively, and additionally more complex {Tl@[Tl(hino)]6}(OTf) (4) compound. Comparison of their antimicrobial activity with selected lead(II) and bismuth(III) analogs and free ligands showed that only bismuth(III) complexes demonstrated significant antimicrobial activity, from two- to fivefold larger than the free ligands.

Enantioselective sensing of (S)-Thalidomide in blood plasma with a chiral naphthalene diimide derivative.

Fast, simple in use and highly effective voltammetric enantiosensor dedicated for determination of thalidomide (TD) enantiomers (especially towards the toxic (S)-enantiomer) in blood plasma is still desirable. Here we have proven that newly synthesized chiral naphthalene diimide (NDI) derivatives are excellent electroactive materials for TD enantiosensors. The recognition process relies on the specific interaction between the chiral NDI receptor and the thalidomide enantiomer of the opposite configuration. This unique specific interaction between (S)-thalidomide and (R)-NDI derivative counterparts, evident in the DPV voltammograms, was confirmed by molecular modeling. The demonstrated voltammetric enantiosensors are characterized by the low detection limit at the level of µg·L-1, wide analytical range from 5·10-4 - 10 mg·L-1, high selectivity and long lifetime. The results of the recovery rates showed a very good degree of accuracy towards the determination of (S)-thalidomide in the blood samples, so it can be successfully used in the analysis of clinical samples.

Chiral Lanthanide Complexes with l- and d-Alanine: An X-ray and Vibrational Circular Dichroism Study.

A whole series of [Ln(H2O)4(Ala)2]26+ dimeric cationic lanthanide complexes with both L- and D-alanine enantiomers was synthesized. The single-crystal X-ray diffraction at 100 and 292 K shows the formation of two types of dimers (I and II) in crystals. Between the dimer centers, the alanine molecules behave as bridging (µ2-O,O'-) and chelating bridging (µ2-O,O,O'-) ligands. The first type of bridge is present in dimers I, while both bridge forms can be observed in dimers II. The IR and vibrational circular dichroism (VCD) spectra of all L- and D-alanine complexes were registered in the 1750-1250 cm-1 range as KBr pellets. Despite all the studied complexes are exhibiting similar crystal structures, the spectra reveal correlations or trends with the Ln-O1 distances which exemplify the lanthanide contraction effect in the IR spectra. This is especially true for the positions and intensities of some IR bands. Unexpectedly, the ν(C=O) VCD bands are quite intense and their composed shapes reveal the inequivalence of the C=O vibrators in the unit cell which vary with the lanthanide. Unlike in the IR spectra, the ν(C=O) VCD band positions are only weakly correlated with the change of Ln and the VCD intensities at most show some trends. Nevertheless, this is the first observation of the lanthanide contraction effect in the VCD spectra. Generally, for the heavier lanthanides (Ln: Dy-Lu), the VCD band maxima are very close to each other and the mirror reflection of the band of two enantiomers is usually better than that of the lighter Lns. DFT calculations show that the higher the multiplicity the higher the stability of the system. Actually, the molecular geometry in crystals (at 100 K) is well predicted based on the highest-spin structures. Also, the simulated IR and VCD spectra strongly depend on the Ln electron configuration but the best overall agreement was reached for the Lu complex, which is the only system with a fully filled f-shell.

Substituent Effect in the First Excited Triplet State of Monosubstituted Benzenes.

The structure of 30 monosubstituted benzenes in the first excited triplet T1 state was optimized with both unrestricted (U) and restricted open shell (RO) approximations combined with the ωB97XD/aug-cc-pVTZ basis method. The substituents exhibited diverse σ- and π-electron-donating and/or -withdrawing groups. Two different positions of the substituents are observed in the studied compounds in the T1 state: one distorted from the plane and the other coplanar with a quinoidal ring. The majority of the substituents are π-electron donating in the first group while π-electron withdrawing in the second one. Basically, U- and RO-ωB97XD approximations yield concordant results except for the B-substituents and a few of the planar groups. In the T1 state, the studied molecules are not aromatic, yet aromaticity estimated using the HOMA (harmonic oscillator model of aromaticity) index increases from ca. -0.2 to ca. 0.4 with substituent distortion, while in the S1 state, they are only slightly less aromatic than in the ground state (HOMA ≈0.8 vs ≈1.0, respectively). Unexpectedly, the sEDA(T1) and pEDA(T1) substituent effect descriptors do not correlate with analogous parameters for the ground and first excited singlet states. This is because in the T1 state, the geometry of the ring changes dramatically and the sEDA(T1) and pEDA(T1) descriptors do not characterize only the functional group but the entire molecule. Thus, they cannot provide useful scales for the substituents in the T1 states. We found that the spin density in the T1 states is accumulated at the Cipso and Cp atoms, and with the substituent deformation angle, it nonlinearly increases at the former while decreases at the latter. It appeared that the gap between singly unoccupied molecular orbital and singly occupied molecular orbital (SUMO-SOMO) is determined by the change of the SOMO energy because the former is essentially constant. For the nonplanar structures, SOMO correlates with the torsion angle of the substituent and the ground-state pEDA(S0) descriptor of the π-electron-donating substituents ranging from 0.02 to 0.2 e. Finally, shapes of the SOMO-1 instead of SOMO frontier orbitals in the T1 state somehow resemble the highest occupied molecular orbital ones of the S0 and S1 states. For several planar systems, the shape of the U- and RO-density functional theory-calculated SOMO-1 orbitals differs substantially.

Endohedral isomerism in model achiral and chiral La@C58N2 systems.

Endohedral structures with La0 or La3+ encapsulated in chiral (1,16)C58N2 or achiral (1,4)C58N2 diazafullerenes were studied at the B3LYP/G-31G*/SDD level. Two stable locations of La0 and La3+ are possible in each cage but only with La0@(1,16)C58N2 can the two isomers coexist. We found that an AIM determined hapticity of the endohedral species selectively differentiates the systems. We predict that there will always exist IR and Raman bands which allow for them to be identified in the presence of the parent cage. For the La0@(1,16) C58N2 molecules and the parent diazafullerene, the Raman spectra are likely to reveal a pre-resonance effect even at 785 nm and it seems possible to selectively excite only one isomer. The calculated electronic spectra suggested a chance to determine the less populated diazafullerene in the presence of the more populated one, be it chiral or achiral. For the chiral endohedral isomers, the calculated VCD spectra are quite dissimilar and the two endohedral isomers and the parent heterofullerene seem to be easily detected. Eventually, we defined the endohedral isomerism as follows: The endohedral isomerism is the phenomenon whereby an internal individuum captured in a cage can occupy more than one stable position without changing the cage connectivity.

On Raman optical activity sign-switching between the ground and excited states leading to an unusual resonance ROA induced chirality.

Raman optical activity (ROA) spectra recorded for a chiral naphthalene diimide derivative (nBu-NDI-BINAM) dissolved in a series of solvents exhibit strong solute to solvent induced chirality with: (1) dominating bands of solvents, (2) nBu-NDI-BINAM resonance ROA (RROA) bands which are barely visible, (3) monosignate RROA Solvent spectra with an unexpected sign concordant with that of the ECD band of the resonant electronic state, (4) bisignate RROA bands for a few solvents, and (5) superposition of non-resonant and resonant ROA bands of the chiral solvents. The unusual ROA enhancement was explained in terms of resonance energy transfer with resonant Raman emission. The surprising RROA sign-switching was found to be due to specific conformational equilibria where one solute conformer dominates in the ground and the other in the first excited singlet state, and, the signs of the related ECD bands of these two conformers are opposite.

Three Queries about the HOMA Index.

HOMA (Harmonic Oscillator Model of Aromaticity) is a simple, successful, and widely used geometrical aromaticity index. However, HOMA can also be used as a general molecular descriptor appropriate for any type of molecule. It reaches the global maximum for benzene, whereas the potent magnetic aromaticity NICS index has no lower or upper limits. Hence, questions arise and go beyond mere differences between the geometric and magnetic aspects of aromaticity: (1) Does a molecule of aromaticity greater than that of benzene, but undisclosed by the HOMA definition, exist? (2) Can the Kekuléne cyclohexatriene moiety with HOMA = 0 exist as a part of a larger system? (3) Can the geometrical aromaticity index be defined better? Our answer to the first query is "It is not likely enough", to the second, "Why not define HOMA using a less mysterious molecule than cyclohexatriene?", and to the third, "It is possible to construct another fair geometrical index, but is it better for evaluating aromaticity?" To find these answers, we have studied: (1) the HOMA and NICS indices of over 50 hexahomosubstituted benzenes, (2) the HOMA, as well as EN and GEO, indices of over 100 triply fused hexasubstituted benzenes, and (3) the HOMA and new Geometrical Auxiliary Index (GAI) , of different unsaturated and saturated, aromatic and aliphatic hydrocarbons including all alkane constitutional isomers composed of up to nine carbon atoms.

Affinity of fentanyl and its derivatives for the σ1-receptor.

Fentanyl and its 11 commercially available derivatives were investigated as to their affinity for the σ1 receptor. The parent compound is a rather poor binder (IC50 = 4973 nM), but its close derivatives (benzylfentanyl or p-fluorofentanyl) have submicromolar affinities. Modelling provides a structural basis for the observed trends in activity.

Molecular dynamics of fentanyl bound to µ-opioid receptor.

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 Trp2936.48 than observed for morphine or unliganded receptor. Another difference is tighter arrangement of the interaction between Asp1473.32 and Tyr3267.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.