Synthesis of Fluorescent Pyrrolo[1,2-a]pyrimidines from Fischer Carbene Complexes as Building Blocks.

A novel synthetic methodology is reported for the synthesis of fluorescent pyrrolo[1,2-a]pyrimidines. Fischer carbene complexes served as the synthetic platform for (3+3) cyclization to form the heterocyclic moiety. The reaction process furnished two products, their ratio being modulated by the metal, base, and solvent. The selectivity exhibited was studied by analyzing the potential energy surface with density functional theory tools. The photophysical properties of absorption and emission were also evaluated. The dyes absorbed at wavelengths of 240-440 nm, depending on the substituents. The maximum emission wavelength was in the range of 470-513 nm, with quantum yields of 0.36-1.0 and a high Stokes shift range of 75-226 nm.

Modeling Adsorption of CO2 in Rutile Metallic Oxide Surfaces: Implications in CO2 Catalysis.

CO2 is the most abundant greenhouse gas, and for this reason, it is the main target for finding solutions to climatic change. A strategy of environmental remediation is the transformation of CO2 to an aggregated value product to generate a carbon-neutral cycle. CO2 reduction is a great challenge because of the large C=O dissociation energy, ~179 kcal/mol. Heterogeneous photocatalysis is a strategy to address this issue, where the adsorption process is the fundamental step. The focus of this work is the role of adsorption in CO2 reduction by means of modeling the CO2 adsorption in rutile metallic oxides (TiO2, GeO2, SnO2, IrO2 and PbO2) using Density Functional Theory (DFT) and periodic DFT methods. The comparison of adsorption on different metal oxides forming the same type of crystal structure allowed us to observe the influence of the metal in the adsorption process. In the same way, we performed a comparison of the adsorption capability between two different surface planes, (001) and (110). Two CO2 configurations were observed, linear and folded: the folded conformations were observed in TiO2, GeO2 and SnO2, while the linear conformations were present in IrO2 and PbO2. The largest adsorption efficiency was displayed by the (001) surface planes. The CO2 linear and folded configurations were related to the interaction of the oxygen on the metallic surface with the adsorbate carbon, and the linear conformations were associated with the physisorption and folded configurations with chemisorption. TiO2 was the material with the best performance for CO2 interactions during the adsorption.

Structure and Conformation of Novel BODIPY Ugi Adducts.

Two novel BODIPY-Ugi (boron dipyrromethene) adducts exhibit peculiar room temperature (T=20 °C) H-1 NMR spectra in that several protons located at the aromatic aniline-type ring are lost in the baseline. This observation revealed the existence of a dynamic conformational process where rotation around the C-N bond is hindered. Variable-temperature H-1 and C-13 NMR spectroscopic analysis confirmed this conclusion; that is, low-temperature spectra show distinct signals for all four aromatic protons below coalescence, whereas average signals are recorded above coalescence (T=+120 °C). Particularly interesting was the rather large difference in chemical shifts for the ortho protons below coalescence, Δδ=1.45 ppm, which was explained based on DFT computational analysis. Indeed, the calculated lowest-energy gas-phase conformation of the BODIPY Ugi adducts locates one half of the aniline-type ring in the shielding anisotropic cone of the bridge phenyl ring in the BODIPY segment. This is in contrast to the solid-state conformation established by X-ray diffraction analysis that shows a nearly parallel arrangement of the aromatic rings, probably induced by crystal packing forces.

Grammatical evolution-based design of SARS-CoV-2 main protease inhibitors.

A series of SARS-CoV-2 main protease (SARS-CoV-2-Mpro) inhibitors were modeled using evolutive grammar algorithms. We have generated an automated program that finds the best candidate to inhibit the main protease, Mpro, of SARS-CoV-2. The candidates were constructed based on a pharmacophore model of the above-mentioned target; relevant moieties of such molecules were modified using data-basis sets with similar chemical behavior to the reference moieties. Additionally, we used the SMILES language to translate 3D chemical structures to 1D words; then, an evolutive grammar algorithm was used to explore the chemical space and obtain new candidates, which were evaluated via the binding energy of molecular coupling assays as an evaluation function. Finally, sixteen molecules were obtained in 3 runs of our program, three of which show promising binding properties as SARS-CoV-2-Mpro inhibitors. One of them, TTO, maintained its relevant binding properties during 100 ns molecular dynamics experiments. For this reason, TTO is the best candidate to inhibit SARS-CoV-2-Mpro. The software we developed for this contribution is available at the following URL: https://github.com/masotelof/GEMolecularDesign.

Effect of the substituents of new coumarin-imidazo[1,2-a]heterocyclic-3-acrylate derivatives on nonlinear optical properties: a combined experimental-theoretical approach.

A series of new coumarin-imidazo[1,2-a]heterocyclic-3-acrylate derivatives 7a-h were synthesized by the Heck reaction between the corresponding 3-(imidazo[1,2-a]pyrimidines)-(2-yl)-2H-chromen-2-ones 4a-e and methyl acrylate in 45-87% yields. The effect of the distinct substituents on third-order nonlinear optical properties was examined, experimentally measuring their nonlinear refractive indexes by the Z-scan technique. Density functional theory and time-dependent density functional theory were utilized with the B3LYP, CAM-B3LYP, PBE (PBEPBE), and M062X functionals on Gaussian09 software to calculate the vertical excitation, relaxation of the brightest excited states, conformation, HOMO-LUMO gaps, oscillator strength, polarizability, and hyperpolarizabilities of all derivatives. Although all acrylates showed a nonlinear response at a certain level of power, the compounds bearing a diethylamino electron-donating group exhibited higher nonlinear refractive index values (∼10-9 cm2 W-1), which is in agreement with the trend in the computational calculations of the first and second hyperpolarization. According to the structural analysis, the electron-withdrawing group (acrylate) is mainly responsible for the loss of coplanarity because of increasing the dihedral angle between the coumarin and imidazo[1,2-a]heterocyclic moieties (to 39.1°). On the other hand, the unsubstituted compound 4a presented the greatest nonlinearity due to its almost coplanar structure (n2 ∼ 10-8 cm2 W-1), highlighting the importance of this feature.

Exploring novel capping framework: high substituent pyridine-hydroxamic acid derivatives as potential antiproliferative agents.

PURPOSE: Histone deacetylases (HDACs) play a vital role in the epigenetic regulation of gene expression due to their overexpression in several cancer forms. Therefore, these enzymes are considered as a potential anticancer drug target. Different synthetic and natural structures have been studied as HDACs inhibitors; based on available structural design information, the capping group is important for the biological activity due to the different interactions in the active site entrance. The present study aimed to analyze high substituted pyridine as a capping group, which included carrying out the synthesis, antiproliferative activity analysis, and docking studies of these novel compounds. METHODS: To achieve the synthesis of these derivatives, four reaction steps were performed, generating desired products 15a-k. Their effects on cell proliferation and gene expression of p21, cyclin D1, and p53 were determined using the sulphorhodamine B (SRB) method and quantitative real-time polymerase chain reaction. The HDAC1, HDAC6, and HDAC8 isoforms were used for performing docking experiments with our 15a-k products. RESULT: The products 15a-k were obtained in overall yields of 40-71%. Compounds 15j and 15k showed the highest antiproliferative activity in the breast (BT-474 and MDA-MB-231) and prostate (PC3) cancer cell lines at a concentration of 10 µM. These compounds increased p21 mRNA levels and decreased cyclin D1 and p53 gene expression. The docking study showed an increment in the strength, and in the number of interactions performed by the capping moiety of the tested molecules compared with SAHA; interactions displayed are mainly van der Waals, π-stacking, and hydrogen bond. CONCLUSION: The synthesized compounds 2-thiophene (15j) and 2-furan (15k) pyridine displayed cell growth inhibition, regulation of genes related to cell cycle progression in highly metastatic cancer cell lines. The molecular coupling analysis performed with HDAC1, HDAC6 and HDAC8 showed an increment in the number of interactions performed by the capping moiety and consequently in the strength of the capping group interaction.

Modeling Adsorption and Optical Properties for the Design of CO2 Photocatalytic Metal-Organic Frameworks.

Four Metal-Organic Frameworks (MOFs) were modeled (IRMOF-C-BF2, IRMOF-C-(2)-BF2, IRMOF-C'-BF2, and IRMOF-C-CH2BF2) based on IRMOF-1. A series of linkers, based on Frustrated Lewis Pairs and coumarin moieties, were attached to IRMOF-1 to obtain MOFs with photocatalytic properties. Four different linkers were used: (a) a BF2 attached to a coumarin moiety at position 3, (b) two BF2 attached to a coumarin moiety in positions 3 and 7, (c) a BF2 attached in the coumarin moiety at position 7, and (d) a CH2BF2 attached at position 3. An analysis of the adsorption properties of H2, CO2, H2O and possible CO2 photocatalytic capabilities was performed by means of computational modeling using Density Functional Theory (DFT), Time-Dependent Density Functional (TD-DFT) methods, and periodic quantum chemical wave function approach. The results show that the proposed linkers are good enough to improve the CO2 adsorption, to hold better bulk properties, and obtain satisfactory optical properties in comparison with IRMOF-1 by itself.

Predictive Global Models of Cruzain Inhibitors with Large Chemical Coverage.

Chagas disease affects 8-11 million people worldwide, most of them living in Latin America. Moreover, migratory phenomena have spread the infection beyond endemic areas. Efforts for the development of new pharmacological therapies are paramount as the pharmacological profile of the two marketed drugs currently available, nifurtimox and benznidazole, needs to be improved. Cruzain, a parasitic cysteine protease, is one of the most attractive biological targets due to its roles in parasite survival and immune evasion. In this work, we compiled and curated a database of diverse cruzain inhibitors previously reported in the literature. From this data set, quantitative structure-activity relationship (QSAR) models for the prediction of their pIC50 values were generated using k-nearest neighbors and random forest algorithms. Local and global models were calculated and compared. The statistical parameters for internal and external validation indicate a significant predictability, with q loo 2 values around 0.66 and 0.61 and external R 2 coefficients of 0.725 and 0.766. The applicability domain is quantitatively defined, according to QSAR good practices, using the leverage and similarity methods. The models described in this work are readily available in a Python script for the discovery of novel cruzain inhibitors.

Design of Fluorescent Coumarin-Hydroxamic Acid Derivatives as Inhibitors of HDACs: Synthesis, Anti-Proliferative Evaluation and Docking Studies.

Coumarin-hydroxamic acid derivatives 7a-k were herein designed with a dual purpose: as antiproliferative agents and fluorescent probes. The compounds were synthesized in moderate yields (30-87%) through a simple methodology, biological evaluation was carried out on prostate (PC3) and breast cancer (BT-474 and MDA-MB-231) cell lines to determine the effects on cell proliferation and gene expression. For compounds 7c, 7e, 7f, 7i and 7j the inhibition of cancer cell proliferation was similar to that found with the reference compound at a comparable concentration (10 µM), in addition, their molecular docking studies performed on histone deacetylases 1, 6 and 8 showed strong binding to the respective active sites. In most cases, antiproliferative activity was accompanied by greater levels of cyclin-dependent kinase inhibitor p21, downregulation of the p53 tumor suppressor gene, and regulation of cyclin D1 gene expression. We conclude that compounds 7c, 7e, 7f, 7i and 7j may be considered as potential anticancer agents, considering their antiproliferative properties, their effect on the regulation of the genes, as well as their capacity to dock to the active sites. The fluorescent properties of compound 7j and 7k suggest that they can provide further insight into the mechanism of action.

Interaction and Reactivity of Cisplatin Physisorbed on Graphene Oxide Nano-Prototypes.

The physical adsorption of cisplatin (CP) on graphene oxide (GO) and reduced graphene oxide (rGO) is investigated at the DFT level of theory by exploiting suitable molecular prototypes representing the most probable adsorbing regions of GO and rGO nano-structures. The results show that the CP binding energy is enhanced with respect to that for the interaction with pristine graphene. This is due to the preferential adsorption of the drug in correspondence of the epoxy and hydroxy groups located on GO basal plane: an energy decomposition analysis of the corresponding binding energy reveals that the most attractive contribution comes from the electrostatic attraction between the -NH 3 ends of CP and the oxygen groups on (r)GO, which can be associated with hydrogen bonding effects. Moreover, it is found that the reactivity of the physically adsorbed CP is practically unaltered being the free energy variation of the first hydrolysis reaction almost matching that of its free (unadsorbed drug) counterpart. The reported results suggest that the CP physical adsorption on GO and rGO carriers is overall feasible being an exergonic process in aqueous solution. The CP adsorption could facilitate its solubility and transport in water solutions, exploiting the high hydrophilicity of the peripheral carboxylic groups located on the edge of the GO and rGO nano-structures. Moreover, the the higher affinity of CP with respect to the oxidized sites suggests a possible dependence of drug loading and release on pH conditions, which would highly facilitate its specific delivery.

Binding of Pb-Melatonin and Pb-(Melatonin-metabolites) complexes with DMT1 and ZIP8: implications for lead detoxification.

We have applied the docking methodology to characterize the binding modes of the divalent metal transporter 1 (DMT1) and the zinc transporter 8 (ZIP8) protein channels with: melatonin, some melatonin metabolites, and a few lead complexes of melatonin and its metabolites, in three different coordination modes (mono-coordinated, bi-coordinated and tri-coordinated). Our results show that bi-coordinated and tri-coordinated lead complexes prefer to bind inside the central region of ZIP8. Moreover, the interaction strength is larger compared with that of the free melatonin and melatonin metabolites. On the other hand, the binding modes with DMT1 of such complexes display lower binding energies, compared with the free melatonin and melatonin metabolites. Our results suggest that ZIP8 plays a major role in the translocation of Pb, bi or tri coordinated, when melatonin metabolites are present. Finally, we have characterized the binding modes responsible for the ZIP8 large affinities, found in bi-coordinated and tri-coordinated lead complexes. Our results show that such interactions are greater, because of an increase of the number of hydrogen bonds, the number and intensity of electrostatic interactions, and the interaction overlay degree in each binding mode. Our results give insight into the importance of the ZIP8 channel on lead transport and a possible elimination mechanism in lead detoxification processes. Graphical abstract .

Computational study of metal complexes formed with EDTA, melatonin, and its main metabolites: implications in lead intoxication and clues to a plausible alternative treatment.

Melatonin has been proposed as an alternative treatment to the usage of EDTA for lead intoxication. In this computational paper, since previous work has not systematically studied the complexes that may be formed in the existing and proposed treatments, we study 45 possible complexes that we suggest may be formed between Pb and some essential metals with melatonin, melatonin metabolites, and EDTA, analyzing the stability and viability of these through the Gibbs free energy of complexation (ΔΔG), molecular orbitals, and energy decomposition analysis at the DFT level of theory PBE/TZ2P. Our findings show that most complexes present exergonic energies of reaction, and thus spontaneous complex formation. In addition, we show that the AMK and 3OHM melatonin metabolites possess electronic and thermodynamic properties adequate to act as lead trapping molecules due to the lower Pauli repulsion energies involved in the complexes they form and their large negative values of ΔΔG. Therefore, it is shown that both melatonin and some of its metabolites may be employed in a viable treatment for lead intoxication through formation of stable Pb-complexes. Graphical abstract Metal complexes formed with EDTA, melatonin, and its main metabolites.

Effects of lead and lead-melatonin exposure on protein and gene expression of metal transporters, proteins and the copper/zinc ratio in rats.

Human lead (Pb) exposure induces many adverse health effects, including some related to lead accumulation in organs. Although lead bio-distribution in the body has been described, the molecular mechanism underlying distribution and excretion is not well understood. The transport of essential and toxic metals is principally mediated by proteins. How lead affects the expression of metal transporter proteins in the principal metal excretory organs, i.e., the liver and kidney, is unknown. Considering that co-administration of melatonin and lead reduces the toxic effects of lead and lead levels in the blood in vivo, we examined how lead and co-administration of lead and melatonin affect the gene and protein expression of metal transporter proteins (ZIP8, ZIP14, CTR1 and DMT1) in these organs. Rats were exposed intraperitoneally to lead or lead-melatonin. Our results show that Pb exposure induces changes in the protein and gene expression of ZIP8, ZIP14 and CTR1. Alterations in the copper/zinc ratio found in the blood, liver and kidney were likely related to these changes. With DMT1 expression (gene and protein), a positive correlation was found with lead levels in the kidney. Co-administration of melatonin and lead reduced lead-induced DMT1 expression through an unknown mechanism. This effect of melatonin relates to reduced lead levels in the blood and kidney. The metal transport protein function and our results suggest that DMT1 likely contributes to lead accumulation in organs. These data further elucidate the effects of lead on Cu and Zn and the molecular mechanism underlying lead bio-distribution in animals.

Noncovalent interactions between cisplatin and graphene prototypes.

Cisplatin (CP) has been widely used as an anticancer drug for more than 30 years despite severe side effects due to its low bioavailability and poor specificity. For this reason, it is paramount to study and design novel nanomaterials to be used as vectors capable to effectively deliver the drug to the biological target. The CP square-planar geometry, together with its low water solubility, suggests that it could be possibly easily adsorbed on 2D graphene nanostructures through the interaction with the related highly conjugated π-electron system. In this work, pyrene has been first selected as the minimum approximation to the graphene plane, which allows to properly study the noncovalent interactions determining the CP adsorption. In particular, electronic structure calculations at the MP2C and DFT-SAPT levels of theory have allowed to obtain benchmark interaction energies for some limiting configurations of the CP-pyrene complex, as well as to assess the role of the different contributions to the total interaction: it has been found that the parallel configurations of the aggregate are mainly stabilized around the minimum region by dispersion, in a similar way as for complexes bonded through π-π interactions. Then, the benchmark interaction energies have been used to test corresponding estimations obtained within the less expensive DFT to validate an optimal exchange-correlation functional which includes corrections to take properly into account for the dispersion contribution. Reliable DFT interaction energies have been therefore obtained for CP adsorbed on graphene prototypes of increasing size, ranging from coronene, ovalene, and up to C150 H30 . Finally, DFT geometry optimizations and frequency calculations have also allowed a reliable estimation of the adsorption enthalpy of CP on graphene, which is found particularly favorable (about -20 kcal/mol at 298 K and 1 bar) being twice that estimated for the corresponding benzene adsorption. © 2017 Wiley Periodicals, Inc.

Hydrogen bond cooperativity and anticooperativity within the water hexamer.

The hydrogen bond (HB), arguably the most important non-covalent interaction in chemistry, is getting renewed attention particularly in materials engineering. We address herein HB non-additive features by examining different structures of the water hexamer (cage, prism, book, bag and ring). To that end, we rely on the interacting quantum atoms (IQA) topological energy partition, an approach that has been successfully used to study similar effects in smaller water clusters (see Chem. - Eur. J., 19, 14304). Our IQA interaction energies, , are used to classify the strength of HBs in terms of the single/double character of the donor and acceptor H2O molecules involved in the interaction. The strongest hydrogen bonds on this new scale entail double donors and acceptors that show larger values of than those observed in homodromic cycles, paradigms of cooperative effects. Importantly, this means that besides the traditional HB anticooperativity ascribed to double acceptors and donors, the occurrence of these species is also related to HB strengthening. Overall, we hope that the results of this research will lead to a further understanding of the HB non-additivity in intramolecular and intermolecular interactions.

Electron correlation in the interacting quantum atoms partition via coupled-cluster lagrangian densities.

The electronic energy partition established by the Interacting Quantum Atoms (IQA) approach is an important method of wavefunction analyses which has yielded valuable insights about different phenomena in physical chemistry. Most of the IQA applications have relied upon approximations, which do not include either dynamical correlation (DC) such as Hartree-Fock (HF) or external DC like CASSCF theory. Recently, DC was included in the IQA method by means of HF/Coupled-Cluster (CC) transition densities (Chávez-Calvillo et al., Comput. Theory Chem. 2015, 1053, 90). Despite the potential utility of this approach, it has a few drawbacks, for example, it is not consistent with the calculation of CC properties different from the total electronic energy. To improve this situation, we have implemented the IQA energy partition based on CC Lagrangian one- and two-electron orbital density matrices. The development presented in this article is tested and illustrated with the H2 , LiH, H2 O, H2 S, N2 , and CO molecules for which the IQA results obtained under the consideration of (i) the CC Lagrangian, (ii) HF/CC transition densities, and (iii) HF are critically analyzed and compared. Additionally, the effect of the DC in the different components of the electronic energy in the formation of the T-shaped (H2 )2 van der Waals cluster and the bimolecular nucleophilic substitution between F(-) and CH3 F is examined. We anticipate that the approach put forward in this article will provide new understandings on subjects in physical chemistry wherein DC plays a crucial role like molecular interactions along with chemical bonding and reactivity. © 2016 Wiley Periodicals, Inc.

Hydrogen-bond cooperative effects in small cyclic water clusters as revealed by the interacting quantum atoms approach.

The cooperative effects of hydrogen bonding in small water clusters (H2 O)n (n=3-6) have been studied by using the partition of the electronic energy in accordance with the interacting quantum atoms (IQA) approach. The IQA energy splitting is complemented by a topological analysis of the electron density (ρ(r)) compliant with the quantum theory of atoms-in-molecules (QTAIM) and the calculation of electrostatic interactions by using one- and two-electron integrals, thereby avoiding convergence issues inherent to a multipolar expansion. The results show that the cooperative effects of hydrogen bonding in small water clusters arise from a compromise between: 1) the deformation energy (i.e., the energy necessary to modify the electron density and the configuration of the nuclei of the isolated water molecules to those within the water clusters), and 2) the interaction energy (Eint ) of these contorted molecules in (H2 O)n . Whereas the magnitude of both deformation and interaction energies is enhanced as water molecules are added to the system, the augmentation of the latter becomes dominant when the size of the cluster is increased. In addition, the electrostatic, classic, and exchange components of Eint for a pair of water molecules in the cluster (H2 O)n-1 become more attractive when a new H2 O unit is incorporated to generate the system (H2 O)n with the last-mentioned contribution being consistently the most important part of Eint throughout the hydrogen bonds under consideration. This is opposed to the traditional view, which regards hydrogen bonding in water as an electrostatically driven interaction. Overall, the trends of the delocalization indices, δ(Ω,Ω'), the QTAIM atomic charges, the topology of ρ(r), and the IQA results altogether show how polarization, charge transfer, electrostatics, and covalency contribute to the cooperative effects of hydrogen bonding in small water clusters. It is our hope that the analysis presented in this paper could offer insight into the different intra- and intermolecular interactions present in hydrogen-bonded systems.

Performance of the density matrix functional theory in the quantum theory of atoms in molecules.

The generalization to arbitrary molecular geometries of the energetic partitioning provided by the atomic virial theorem of the quantum theory of atoms in molecules (QTAIM) leads to an exact and chemically intuitive energy partitioning scheme, the interacting quantum atoms (IQA) approach, that depends on the availability of second-order reduced density matrices (2-RDMs). This work explores the performance of this approach in particular and of the QTAIM in general with approximate 2-RDMs obtained from the density matrix functional theory (DMFT), which rests on the natural expansion (natural orbitals and their corresponding occupation numbers) of the first-order reduced density matrix (1-RDM). A number of these functionals have been implemented in the promolden code and used to perform QTAIM and IQA analyses on several representative molecules and model chemical reactions. Total energies, covalent intra- and interbasin exchange-correlation interactions, as well as localization and delocalization indices have been determined with these functionals from 1-RDMs obtained at different levels of theory. Results are compared to the values computed from the exact 2-RDMs, whenever possible.

Nature of Chemical Interactions from the Profiles of Electron Delocalization Indices.

We analyze the behavior of the profiles of delocalization indices (DIs) between relevant pairs of atoms along reaction coordinates for a set of model chemical processes. A relationship between the topology of the DI and the nature of the underlying chemical change is reported. As shown, exponential shapes correspond to the traditional category of repulsive/nonbonded interactions, while sigmoidal profiles signal the formation/breaking of chemical links.

Energetic and electron density analysis of hydrogen dissociation of protonated benzene.

The electronic structure of the benzenium cation, [C6H7]+, the simplest intermediate of electrophilic aromatic substitution reactions, was analyzed in terms of the properties of electron densities obtained from multiconfigurational quantum theoretical methods. The indirect C-H coupling constants and the physical contributions to their values were calculated and rationalized in terms of the electron delocalization between the quantum topological atoms in the molecule. The evolution of the electronic structure for the intramolecular proton migration and for the dissociation into [C6H6]+ + H, or [C6H5]+ + H2 was also studied. The potential energy surface for intramolecular H migration has six equivalent transition states and two equivalent two-fold saddles, whereas each dissociation process occurs without the presence of any transition state. The calculated energy barriers of 9.5, 80.3 and 72.4 kcal mol(-1) for the intramolecular proton migration, H and H2 eliminations, respectively, agree with experimental reports. The quantitative chemical descriptors based on the electron density of the benzenium cation provide insight on the nature of the chemical bond, including electron delocalization, and its evolution during chemical transformations of the molecule.