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.

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.

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.

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.

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.

Fentanyl Family at the Mu-Opioid Receptor: Uniform Assessment of Binding and Computational Analysis.

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).

The Vibrational Circular Dichroism Pattern of the ν(C=O) Bands in Isoindolinones.

The IR and vibrational circular dichroism (VCD) spectra of both enantiomers of Me-, iPr-, nBu-, Ph-, and CH2 Ph-substituted isoindolinones in solution and KBr pellets were measured and interpreted by DFT calculations. The spectra in solution revealed no important differences in the C=O stretching vibration region while the interpretation of very distinct spectra taken in pellets required determining the crystal structures. The studied compounds crystallized in the P21 21 21 (Me, iPr, CH2 Ph), P31 (nBu), and P21 (Ph) space groups. We found that the quality of simulated spectra strongly depends on the substituent, the structure of the molecular cluster assumed, basis set, and use of the dispersion correction. The IR spectra can be reproduced well based on the simplest linear arrangement of hydrogen-bonded chains mimicking the molecular arrangement in the crystals. We found no common approach to reproduce all the registered VCD spectra in the crystal phase. For the Me and nBu isoindolinones, the VCD pattern was the best reproduced by full optimization of the selected large molecular clusters. For iPr, Ph and CH2 Ph derivatives optimizing only the position of H-atoms in a fragment frozen as in the crystal provides the best results. Such an approach can reduce the computation time from months to one week.

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

sEDA, pEDA, and cSAR descriptors of the substituent effect were determined for >30 monosubstituted benzenes in the first excited singlet S1 state at the LC-ωB97XD/aug-cc-pVTZ level. It was found that in the S1 state, the σ- and π-valence electrons are a bit less and a bit more affected, respectively, than in the S0 state, but basically, the effect in both states remains the same. In the S0 and S1 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(S1), LUMO(S1), and HOMO(S1)-LUMO(S1) gap correlate also with the pEDA(S1) and cSAR(S1) 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(S0) and HOMO(S1) 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.

Structural Insights into σ1 Receptor Interactions with Opioid Ligands by Molecular Dynamics Simulations.

Despite considerable advances over the past years in understanding the mechanisms of action and the role of the σ1 receptor, several questions regarding this receptor remain unanswered. This receptor has been identified as a useful target for the treatment of a diverse range of diseases, from various central nervous system disorders to cancer. The recently solved issue of the crystal structure of the σ1 receptor has made elucidating the structure-activity relationship feasible. The interaction of seven representative opioid ligands with the crystal structure of the σ1 receptor (PDB ID: 5HK1) was simulated for the first time using molecular dynamics (MD). Analysis of the MD trajectories has provided the receptor-ligand interaction fingerprints, combining information on the crucial receptor residues and frequency of the residue-ligand contacts. The contact frequencies and the contact maps suggest that for all studied ligands, the hydrophilic (hydrogen bonding) interactions with Glu172 are an important factor for the ligands' affinities toward the σ1 receptor. However, the hydrophobic interactions with Tyr120, Val162, Leu105, and Ile124 also significantly contribute to the ligand-receptor interplay and, in particular, differentiate the action of the agonistic morphine from the antagonistic haloperidol.

Chiral Thiophene Sulfonamide-A Challenge for VOA Calculations.

Two enantiomers of 2-methyl-N-(1-thien-2-ylethyl)propane-2-sulfonamide (TSA) were synthesized, and their VCD, ROA, IR, and Raman spectra were registered. The solved (S)-TSA X-ray structure shows a disorder connected to the presence of two TSA conformers differing by a slight rotation of the thiophene ring. Two molecules in the unit cell of the monoclinic P21 crystal form a net of NH···OS and C*H···OS hydrogen bonds. Out of a series of computational levels tested to interpret the spectra, the B3LYP functional combined with the def2TZVP basis set satisfactorily reproduces the experimental VCD and ROA spectra. To simulate the VCD spectra of TSA enantiomers in KBr pellets, dimers and tetramers, with two different positions of the thiophene ring, were considered. The VCD spectra measured in CDCl3 are completely different from those taken in KBr due to the conformational freedom of TSA in chloroform. Seven TSA conformers fall into two groups of opposite configurations at the pyramidal N atom forming the additional stereogenic center. However, the barriers between conformers in each group are lower than the energy of thermal motions at 300 K. Thus, all conformers, but the most stable in each group, are likely to be metastable states. The calculated IR, VCD, Raman, and ROA spectra of the conformers depend not only on the type of stereogenic N atom but also on the thiophene ring rotation. Yet, they are likely to coexist because of low barriers between them. Three approaches were tested to reproduce the chiroptical spectra in solution using PCM and hybrid solvation models. As a consequence, it was found that a model in which all conformers contribute to the spectra with equal population factors seems to best reproduce the experimental data. Such a result suggests that in a dissolved state in 300 K TSA occurs in a very shallow potential well and all of its conformers coexist.

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.

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.

The sEDA(=) and pEDA(=) descriptors of the double bonded substituent effect.

New descriptors of the double bonded substituent effect, sEDA(=) and pEDA(=), were constructed based on quantum chemical calculations and NBO methodology. They show to what extent the σ and π electrons are donated to or withdrawn from the substituted system by a double bonded substituent. The new descriptors differ from descriptors of the classical substituent effect for which the pz orbital of the ipso carbon atom is engaged in the π-electron system of the two neighboring atoms in the ring. For double bonded substituents, the pz orbital participates in double bond formation with only one external atom. Moreover, the external double bond forces localization of the double bond system of the ring, significantly changing the core molecule. We demonstrated good agreement between our descriptors and the Weinhold and Landis' "natural σ and π-electronegativities": so far only descriptors allowing for evaluation of the substitution effect by a double bonded atom. The equivalency between descriptors constructed for 5- and 6-membered model structures as well as linear dependence/independence of the constructed parameters was discussed. Some interrelations between sEDA(=) and pEDA(=) and the other descriptors of (hetero)cyclic systems such as aromaticity and electron density in the ring and bond critical points were also examined.

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.

Heteroatom incorporation effect in σ- and π-electron systems: the sEDA(II) and pEDA(II) descriptors.

The effect of heteroatom or heteroatomic group incorporation into unsaturated five- and six-membered cyclic systems was studied by means of DFT/B3LYP/aug-cc-pVDZ calculations. Two descriptors of the incorporation effect, sEDA(II) and pEDA(II), reflecting the influence of the incorporated atom or group on the population of the σ and π valence electrons, were constructed on the basis of natural bond orbital analysis. The sEDA(II) and pEDA(II) descriptors were shown to be linearly independent; the former correlated very well with electronegativity scales, whereas the latter correlated with NICS(1)(ZZ) and HOMA(CC) aromaticity indices. The two descriptors seem to be universal tools for analyzing different chemical and physicochemical effects occurring in unsaturated heterocyclic systems.

Chirality measures of α-amino acids.

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.

On vibrational circular dichroism chirality transfer in electron donor-acceptor complexes: a prediction for the quinine···BF3 system.

Vibrational circular dichroism (VCD) chirality transfer occurs when an achiral molecule interacts with a chiral one and becomes VCD-active. Unlike for H-bonds, for organic electron donor-acceptor (EDA) complexes this phenomenon remains almost unknown. Here, the VCD chirality transfer from chiral quinine to achiral BF3 is studied at the B3LYP/aug-cc-pVDZ level. Accessibility of four quinine electron donor sites changes with conformation. Therefore, the quinine conformational landscape was explored and a considerable agreement between X-ray and the most stable conformer geometries was achieved. The BF3 complex through the aliphatic quinuclidine N atom is definitely dominating and is predicted to be easily recognizable in the VCD spectrum. Out of several VCD chirality transfer modes, the ν(s)(BF3) mode, the most intense in the entire VCD spectrum, satisfies the VCD mode robustness criterion and can be used for monitoring the chirality transfer phenomenon in quinine···BF3 system.

On stability, chirality measures, and theoretical VCD spectra of the chiral C58X2 fullerenes (X = N, B).

The stability of all 23 C(58)N(2) and C(58)B(2) heterofullerenes in the singlet and triplet states was determined at the B3LYP/6-31G** level. In equilibrium mixture the achiral (1,4) C(58)N(2) isomer would be populated in ca. 95.8%, the chiral (1,16) one in ca. 3.3%, and the achiral (1,4) C(58)B(2) in 100%, whereas all triplet state isomers are less stable. Fourteen out of 23 C(58)X(2) are chiral. Four different chirality measures were calculated by our own CHIMEA program: pure geometrical, labeled, mass, and charge. Intercorrelations between the measures for all chiral compounds indicate that the pure geometrical chirality measure is unstable and should not be used in QSAR predictions of the other molecular properties, while the labeled and mass-weighted ones are promising QSAR descriptors. For each chiral C(58)N(2) molecule, some very strong VCD bands, of intensity comparable with that in the IR spectra, can serve in identification and characterization of the isomers.