Ferran Feixas

On the performance of some aromaticity indices: a critical assessment using a test set.

Aromaticity is a central chemical concept widely used in modern chemistry for the interpretation of molecular structure, stability, reactivity, and magnetic properties of many compounds. As such, its reliable prediction is an important task of computational chemistry. In recent years, many methods to quantify aromaticity based on different physicochemical properties of molecules have been proposed. However, the nonobservable nature of aromaticity makes difficult to assess the performance of the numerous existing indices. In the present work, we introduce a series of fifteen aromaticity tests that can be used to analyze the advantages and drawbacks of a group of aromaticity descriptors. On the basis of the results obtained for a set of ten indicators of aromaticity, we conclude that indices based on the study of electron delocalization in aromatic species are the most accurate among those examined in this work.

A Critical Assessment of the Performance of Magnetic and Electronic Indices of Aromaticity

Abstract: The lack of reference aromatic systems in the realm of inorganic aromatic compounds makes the evaluation of aromaticity in all-metal and semimetal clusters a difficult task. To date, calculation of nucleus-independent chemical shifts (NICS) has been the most widely used method to discuss aromaticity in these systems. In the first part of this work, we briefly review our previous studies, showing some pitfalls of the NICS indicator of aromaticity in organic molecules. Then, we refer to our study on the performance of some aromaticity indices in a series of 15 aromaticity tests, which can be used to analyze the advantages and drawbacks of aromaticity descriptors. It is shown that indices based on the study of electron delocalization are the most accurate among those analyzed in the series of proposed tests, while NICS(1) zz and NICS(0) πzz present the best behavior among NICS indices. In the second part, we discuss the use of NICS and electronic multicenter indices (MCI) in inorganic clusters. In particular, we evaluate the aromaticity of two series of all-metal and semimetal clusters with predictable aromaticity trends by means of NICS and MCI. Results show that the expected trends are generally better reproduced by MCI than NICS. It is concluded that NICS(0)

Electron Localization Function at the Correlated Level: A Natural Orbital Formulation.

In this work we present a 2-fold approximation for the calculation of the electron localization function (ELF) which avoids the use of the two-particle density (2-PD). The first approximation is used for the calculation of the ELF itself and the second one is used to approximate pair populations integrated in the ELF basins. Both approximations only need the natural orbitals and their occupancies, which are available for most methods used in electronic structure calculations. In this way, methods such as CCSD and MP2 can be used for the calculation of the ELF despite the lack of a pertinent definition of the 2-PD. By avoiding the calculation of the 2-PD, the present formulation provides the means for routine calculations of the ELF in medium-size molecules with correlated methods. The performance of this approximation is shown in a number of examples.

Exploring the Role of Receptor Flexibility in Structure-Based Drug Discovery

The proper understanding of biomolecular recognition mechanisms that take place in a drug target is of paramount importance to improve the efficiency of drug discovery and development. The intrinsic dynamic character of proteins has a strong influence on biomolecular recognition mechanisms and models such as conformational selection have been widely used to account for this dynamic association process. However, conformational changes occurring in the receptor prior and upon association with other molecules are diverse and not obvious to predict when only a few structures of the receptor are available. In view of the prominent role of protein flexibility in ligand binding and its implications for drug discovery, it is of great interest to identify receptor conformations that play a major role in biomolecular recognition before starting rational drug design efforts. In this review, we discuss a number of recent advances in computer-aided drug discovery techniques that have been proposed to incorporate receptor flexibility into structure-based drug design. The allowance for receptor flexibility provided by computational techniques such as molecular dynamics simulations or enhanced sampling techniques helps to improve the accuracy of methods used to estimate binding affinities and, thus, such methods can contribute to the discovery of novel drug leads.

Electron delocalization and aromaticity in low-lying excited states of archetypal organic compounds

Aromaticity is a property usually linked to the ground state of stable molecules. Although it is well-known that certain excited states are unquestionably aromatic, the aromaticity of excited states remains rather unexplored. To move one step forward in the comprehension of aromaticity in excited states, in this work we analyze the electron delocalization and aromaticity of a series of low-lying excited states of cyclobutadiene, benzene, and cyclooctatetraene with different multiplicities at the CASSCF level by means of electron delocalization measures. While our results are in agreement with Bairds rule for the aromaticity of the lowest-lying triplet excited state in annulenes having 4nπ-electrons, they do not support Soncini and Fowlers generalization of Bairds rule pointing out that the lowest-lying quintet state of benzene and septet state of cyclooctatetraene are not aromatic.

Accelerated molecular dynamics simulations of protein folding

Folding of four fast‐folding proteins, including chignolin, Trp‐cage, villin headpiece and WW domain, was simulated via accelerated molecular dynamics (aMD). In comparison with hundred‐of‐microsecond timescale conventional molecular dynamics (cMD) simulations performed on the Anton supercomputer, aMD captured complete folding of the four proteins in significantly shorter simulation time. The folded protein conformations were found within 0.2–2.1 Å of the native NMR or X‐ray crystal structures. Free energy profiles calculated through improved reweighting of the aMD simulations using cumulant expansion to the second‐order are in good agreement with those obtained from cMD simulations. This allows us to identify distinct conformational states (e.g., unfolded and intermediate) other than the native structure and the protein folding energy barriers. Detailed analysis of protein secondary structures and local key residue interactions provided important insights into the protein folding pathways. Furthermore, the selections of force fields and aMD simulation parameters are discussed in detail. Our work shows usefulness and accuracy of aMD in studying protein folding, providing basic references in using aMD in future protein‐folding studies.

A Dissected Ring Current Model for Assessing Magnetic Aromaticity: A General Approach for both Organic and Inorganic Rings

A model based on classical electrodynamics is used to measure the strength of ring currents of different molecular orbitals, i.e., σ‐ and π‐orbitals, and characteristics of ring current loops, i.e., ring current radii and height of current loops above/below the ring planes, among a number of organic as well as inorganic molecules. For the π‐current, the present model represents an improvement of previous approaches to determine ring current intensity. It is proven that the present model is more precise than previous models as they could not explain presence of the minimum in the plot of NICSπzz versus distance close to the ring plane. Variations in the charge of molecules and the types of constituent atoms of each species affect the ring current radii of both σ‐ and π‐current loops as well as the height of π‐current loops above/below the ring plane. It is suggested that variation in the distribution of the one‐electron density in different systems is the main source of differences of the ring current characteristics.

Electronic States of o-Nitrobenzaldehyde: A Combined Experimental and Theoretical Study

The experimental UV/vis absorption spectrum of ortho-nitrobenzaldehyde (o-NBA) has been assigned by means of MS-CASPT2/CASSCF, TD-DFT, and RI-CC2 theoretical computations. Additional information on the nature of the absorbing bands was obtained by comparing the o-NBA spectrum with that of related compounds, as, e.g., nitrobenzene and benzaldehyde. For wavelengths larger than approximately 280 nm, the absorption spectrum of o-NBA is dominated by a series of weak n pi* absorptions from the NO2 and CHO groups. These weak transitions are followed in energy by a more intense band, peaking at 250 nm and arising from charge transfer pi pi* excitations involving mainly benzene and nitro orbitals. Finally, the most intense band centered at 220 nm has its origin in the overlap of two different absorptions: the first one localized in the NO2 substituent and the second one arising from a charge transfer excitation involving the NO2 and the CHO fragments, respectively.

Peculiarities of multiple Cr-Cr bonding. Insights from the analysis of domain-averaged Fermi holes.

The recently proposed methodology known as the analysis of domain-averaged Fermi holes was applied to reveal the nature and peculiarity of metal-metal bonding interactions in two recently reported complexes with an ultrashort, presumably quintuple Cr-Cr bond. The results of the analysis straightforwardly confirm the considerable reduction of the Cr-Cr bond order resulting from depletion of the electron density from one of the electron pairs involved in the metal-metal bonding. Because of this depletion, the Cr-Cr bond can be best classified as the effective quadruple bond, with the contribution of another weak component corresponding to the antiferromagnetic coupling of electrons in one of the available delta-electron pairs.

Chemical bonding and aromaticity in metalloporphyrins

We report here the chemical bonding and aromaticity patterns in metalloporphyrins, which were obtained with density functional theory (DFT) calculations at the OPBE/TZP level. This level of theory was previously shown to be very accurate for determining spin-state splittings [J. Chem. Theory Comput. 2008, 4, 2057] of transition-metal complexes. We considered metalloporphyrins along the first-row transition metals (Sc–Zn) extended with alkaline-earth metals (Mg, Ca) and several second-row transition metals (Ru, Pd, Ag, Cd). An energy decomposition analysis was performed to study the metal–ligand interactions, which showed that almost all complexes are significantly stabilized through (covalent) orbital interactions. The only exception is with calcium as the central metal, which interacts with the porphyrin mainly through electrostatic interactions. Furthermore, we studied aromaticity patterns for these complexes by looking at a number of (structural and electronic) aromaticity descriptors, for both the inn...