Quantification of thermal ring flexibilities of aromatic and heteroaromatic compounds

The consequences of thermal fluctuations occurring at room temperatures on the aromatic character of a broad group of compounds were analyzed in three distinct ways. First of all, the ring deformations were modeled along normal coordinates coming from quantum thermo-chemistry computations. The amplitudes of vibrations were estimated according to absorbed energies at room temperature. Alternatively, in-plane and out-of-plane ring deformations were modeled via scanning procedure with partial relaxation of the molecular geometry. The influence of ring deformations on π–electron delocalization was expressed in terms of HOMA values. Besides, the ring deformability was defined as the averaged change of bond angles or dihedral angles constituting the ring that was associated with 1.5 kcal mol^-1 increase of the system energy. The molecules structures adopted during vibrations at room temperature can lead to significant heterogeneity of structural index of aromaticity. The broad span of HOMA values was obtained for analyzed five- or six-membered aromatic and heteroaromatic rings. However, the averaged values obtained for such fluctuations almost perfectly match HOMA values of molecule in the ground state. It has been demonstrated that the ring deformability imposed by bond angle changes is much smaller than for dihedral angles with the same rise of system energy. Interestingly in the case of out-of-plane vibrations modeled by scanning procedure there is observed linear correlation between ring deformability and HOMA values. Proposed method for inclusion of thermal vibrations in the framework of π–electron delocalization provides natural shift of the way of thinking about aromaticity from a static quantity to a dynamic and heterogeneous one due to inclusion of a more realistic object of analysis – thermally deformed structures. From this perspective the thermal fluctuations are supposed to be non-negligible contributions to aromaticity phenomenon.     

Homodesmotic reactions in helicene and phenacene molecules

Helicene molecules were analysed on the basis of their intrinsic aromaticity, and compared with the corresponding phenacene. The analysis of aromaticity was performed by homodesmotic reactions, and by the nucleus independent chemical shifts (NICS) index in the centre of the aromatic rings as well as in the centre of the helix. The results are very interesting, because the NICS shows that there is a strong antiaromatic zone in the centre of the molecule, and a strong aromatic zone in the strip. Therefore, these molecules represent a medley of electronic effects.     

A study of the aromaticity of heteroannelated cyclooctatetraene derivatives

The aromaticity of the rings of thiophene, pyrrole, furan, and benzene annelated cyclooctatetraene (COT) derivatives and of their double charged ions was studied using the graph-theoretical theory of aromaticity. On the basis of topological resonance energy, it was found that the global aromaticity is dependent upon on the arrangement of heteroatoms in the given molecule. Relative stability of these molecules when in different charged states can been explained in terms of the topological charge stabilization rule. We expect that fusing the COT ring with an increasing number of aromatic rings will lead to an increase in the aromaticity of the molecule. According to the bond resonance energy (BRE) and circuit resonance energy (CRE) indices, local antiaromaticity of the COT ring is weakened as the number of fused rings increases, and these changes play a significant role in the global aromaticity of the molecule. For some compounds, our BRE and CRE indices do not predict the same order of magnitude of the local aromatic character of certain rings that the nucleus independent chemical shift (NICS(0) and (NICS(1)) methods predict. Finally, for the available compounds, correlations between the diatropic and paratropic chemical shifts of the protons and our ring current results were analyzed and good agreement was found.     

MacMolecular: A program for visualization of molecular structures on the Macintosh

MacMolecular displays small- to medium-sized biomolecules, with particular emphasis on peptides. It has been developed to run on color Macintosh computers. The display can be stick, ball and stick, depth cued by thickness stick, or several types of space-filling representations. The program takes input from standard PDB files, simple Cartesian coordinate files, and, in addition, from Kinemage files in which atom information has been included. The program allows color changes of various types as well as the normal functions of translation, rotation, and zooming. In addition, animation files may be produced for subsequent display. Bonding of atoms is done by a distance algorithm (standard) or sequentially to properly display Cα traces and traces of peptides containing simplified representations of amino acids. Stereo viewing is available, and manipulated structures which were drawn from PDB files can be written out to new PDB files. In addition, PICT files of the drawing window can be generated.     

Solvent impact on the aromaticity of benzene analogues: implicit versus explicit solvent approach

Solvent impact on the structural index of aromaticity was modelled by polarised continuum field approximation (IEFPCM) and hybrid quantum chemistry (QM/MM) method. Significant solvent related relaxation of the solutes geometries were noticed especially for highly polar species. The significant reduction of the aromaticity was observed for some aromatic compounds in water solution compared to gas phase. The rationale of this fact was provided based on dipole moment changes, energy penalty for polarisation of solute and the distribution of frontier orbital densities. The incoherent predictions of explicit and implicit solvation models are noticed since in some cases the PCM approach artificially exaggerate the geometry relaxation in solution which is not observed if explicit solvent molecules are taken into account.     

A comparative study of the aromaticity of pyrrole, furan, thiophene, and their aza-derivatives

The relative aromaticity of pyrrole, furan, thiophene, and their aza-derivatives has been examined using TRE (topological resonance energy), MRE (magnetic resonance energy), ring current (RC), and ring current diamagnetic susceptibility (χ_G) methods. The results obtained were compared with results obtained by others who used the energetic method ASE (aromatic stabilization energy), the geometric method HOMA (harmonic oscillator model of aromaticity), and the magnetic method NICS(1) (nucleus-independent chemical shift). The impact of nitrogen atoms on the aromaticity of the aza-derivatives of pyrrole, furan, and thiophene is discussed. An excellent correlation was found between the energetic (TRE, MRE) and magnetic (RC and χ_G) criteria of aromaticity for all compounds. It was expected that inclusion of a heteroatom would decrease the aromaticity relative to the cyclopentadienyl anion. Our results show that the type of the first heteroatom, which donates two electrons to the system, as well as the number of nitrogen atoms and their positions in the molecule have a strong effect on aromaticity. In general, aromaticity is enhanced when the nitrogen atom is adjacent to the first heteroatom. The magnitude of aromaticity is related closely with the uniformity of distribution of π-electrons in the molecule.     

Aromatic character of heptafulvene and its complexes with halogen atoms

Geometry optimization of heptafulvene-halogen complexes (halogens: F, Cl, Br, I, and At) carried out at the B3LYP/6-311+G(d,p) level of theory allowed us to estimate the geometry-based aromaticity index HOMA, the magnetism-based indices NICS, NICS(1), and NICS(1)(ZZ), as well as the energy of complex formation. Application of the NBO method allowed us to estimate the pEDA characteristics of the π-electron distribution in complexes (i.e., the electron excess/deficiency of the π-electron system in the ring). All of the characteristics of the complexes were found to be mutually interrelated, exhibiting good or at least acceptable correlation coefficients. It was also noted that halogen atoms with greater radii yielded weaker complexes and lower aromaticities, as shown by the HOMA, NICS, and pEDA indices. The energy of complex formation was observed to be linearly correlated with the degree of aromatization of the heptafulvene ring, as expressed by these indices.     

Location of the extrinsic subunit PsbP in photosystem II studied by pulsed electron-electron double resonance

The binding site of the extrinsic protein PsbP in plant photosystem II was mapped by pulsed electron-electron double resonance, using mutant spinach PsbP (Pro20Cys, Ser82Cys, Ala111Cys, and Ala186Cys) labeled with 4-maleimido-TEMPO (MSL) spin label. The distances between the spin label and the Tyr160 neutral radical (Y_D) in PsbD, the D2 subunit of plant photosystem II, were 50.8?±?3.5?Å, 54.9?±?4.0?Å, 57.8?±?4.9?Å, and 58.4?±?14.1?Å, respectively. The geometry inferred from these distances was fitted to the PsbP crystal structure (PDB: 4RTI) to obtain the coordinates of Y_D relative to PsbP. These coordinates were then fitted under boundary conditions to the structure of cyanobacterial photosystem II (PDB: 4UB6), by rotating on Euler angles centered at fixed Y_D coordinates. The result proposed two models which show possible acidic amino acid residues in CP43, CP47 and D2 that can bind the basic amino acids Arg48, Lys143, and Lys160 in PsbP.     

Specificities of boron disubstituted sumanenes

In this article we focused on computational research of sumanenes disubstituted by boron where the two carbon atoms are substituted by two boron atoms. Disubstitution of rim carbon atoms with boron atoms significantly affected the geometry of the bowl. The main stability factors were used to determine the stability of isomers. The most stable, the shallowest and the deepest isomers were subjected to further study of NMR parameters, chemical shielding and NICS, aromaticity, bowl to bowl inversion barrier and NBO/NPA analysis. The introduction of boron atoms significantly affected the above parameters, changing the aromatic nature of rings, reducing bowl to bowl inversion barrier and produced charge transfer. The NICS are correlated with bowl depth having the result that the function of the fourth degree of bowl depth does not only correlate well to the bowl to bowl inversion barrier with bowl depth, but also finely correlates the change of the NICS and NICSzz with bowl depth.     

Using quantum mechanics to improve estimates of amino acid side chain rotamer energies

Amino acid side chains adopt a discrete set of favorable conformations typically referred to as rotamers. The relative energies of rotamers partially determine which side chain conformations are more often observed in protein structures and accurate estimates of these energies are important for predicting protein structure and designing new proteins. Protein modelers typically calculate side chain rotamer energies by using molecular mechanics (MM) potentials or by converting rotamer probabilities from the protein database (PDB) into relative free energies. One limitation of the knowledge‐based energies is that rotamer preferences observed in the PDB can reflect internal side chain energies as well as longer‐range interactions with the rest of the protein. Here, we test an alternative approach for calculating rotamer energies. We use three different quantum mechanics (QM) methods (second order Møller‐Plesset (MP2), density functional theory (DFT) energy calculation using the B3LYP functional, and Hartree‐Fock) to calculate the energy of amino acid rotamers in a dipeptide model system, and then use these pre‐calculated values in side chain placement simulations. Energies were calculated for over 36,000 different conformations of leucine, isoleucine, and valine dipeptides with backbone torsion angles from the helical and strand regions of the Ramachandran plot. In a subset of cases these energies differ significantly from those calculated with standard molecular mechanics potentials or those derived from PDB statistics. We find that in these cases the energies from the QM methods result in more accurate placement of amino acid side chains in structure prediction tests. Proteins 2008. © 2007 Wiley‐Liss, Inc.     

Ab initio and DFT study of the aromaticity of some Fulvalenes derived from Methylidenecyclopropabenzene

Methylidencyclopropabenzene (MCPB) 1 and Fulvalenes 2–4 are molecules of special interest due to the relation between structure and aromaticity. The aim of this work was to analyze this relation and to quantify the aromaticity in 1–4 using different methods. Magnetic properties are directly related with aromaticity; here we studied the magnetic susceptibility and the anisotropy of the magnetic susceptibility. Nucleus indepedent chemical shift (NICS) and the anisotropy of the induced current density (ACID) were also employed. Tools of very different nature, geometric indexes HOMA and Bird, were determinated too for 1–4. All of these measures were found to be in agreement. Figure Both spatial NICS and ACID plot allow to show the aromaticity/antiaromaticity of a ring     

Aromaticity of graphene nanoflakes in a new way: fragment analysis by combination of the nucleus-independent chemical shifts and the anisotropy of current induced density

A graphene nanoflake (GNF) is a polycyclic aromatic hydrocarbon (PAH) with a huge two-dimensional π-conjugated carbon material in which a central benzene ring is surrounded by identical benzene-type rings through infinite alternant method. In this paper, we explore the structure-aromaticity relationship of the GNFs and the GNFs with hollow sites (GNFHs) by combining the nucleus-independent chemical shifts (NICS) with the anisotropy of the current induced density (ACID). Firstly, the benzene is a typical aromatic molecule (NICS = ?9.671 ppm), GNFs 1-6 is darned with benzene and the corresponding GNFHs 1′-6′. Secondly, the NICS values of GNFs 1-6 alternately vary: ?1.214 (1) > ?13.847 (2) < ?2.662 (3) > ?14.530 (4) < ?3.932 (5) > ?13.978 (6) ppm, the GNFs (2, 4, 6) with even fragments of annulene have larger aromaticity than that of GNFs (1, 3, 5) with odd fragments of annulene. Significantly, the NICS values of GNFs 1-6 can also be fragment analyzed by the NICS values and ACID of benzene and corresponding GNFHs 1′-6′. The NICS values for GNFs (2, 4, 6) can be roughly estimated by the NICS value of benzene minus the NICS value of the GNFHs (2′, 4′, 6′), respectively. The NICS values for GNFs (1, 3, 5) can be roughly estimated by the NICS value of the GNFHs (1′, 3′, 5′) minus the NICS value of benzene, respectively. We hope that the present work can provide a simple and reliable method for the rational design of the GNF with aromaticity, which may be used to understand the origin of the graphene nanoflake aromatic properties.     

Disulfide by Design: a computational method for the rational design of disulfide bonds in proteins

SUMMARY: Disulfide by Design is a program for the design of novel disulfide bonds in proteins. Protein structure files in PDB format are analyzed to identify residue pairs that are likely to form a disulfide bond if the respective amino acids are mutated to cysteines. The output displays residue pairs having the appropriate geometric characteristics for disulfide formation and provides automated generation of modified PDB files including modeled disulfides. Validation demonstrates a high level of accuracy for the algorithm. AVAILABILITY: http://www.ehscenter.org/dbd/ Supplementary information: http://www.ehscenter.org/dbd/     

MutaProt: a web interface for structural analysis of point mutations

A web-based tool, termed 'MutaProt', is described which analyses pairs of PDB files whose members differ in one, or two, amino acids. MutaProt examines the micro environment surrounding the exchanged residue(s) and can be searched by specifying a PDB ID, keywords, or any pair of amino acids. Detailed information about accessibility of the exchanged residue(s) and its atomic contacts are provided based on CSU software (Sobolev et al., Bioinformatics, 15, 327-332, 1999). An interactive 3D presentation of the superimposed regions around the mutation(s) is included. MutaProt is updated weekly.     

Substrate diffusion and oxidation in GMC oxidoreductases: an experimental and computational study on fungal aryl-alcohol oxidase

AAO (aryl-alcohol oxidase) provides H?O? in fungal degradation of lignin, a process of high biotechnological interest. The crystal structure of AAO does not show open access to the active site, where different aromatic alcohols are oxidized. In the present study we investigated substrate diffusion and oxidation in AAO compared with the structurally related CHO (choline oxidase). Cavity finder and ligand diffusion simulations indicate the substrate-entrance channel, requiring side-chain displacements and involving a stacking interaction with Tyr?2. Mixed QM (quantum mechanics)/MM (molecular mechanics) studies combined with site-directed mutagenesis showed two active-site catalytic histidine residues, whose substitution strongly decreased both catalytic and transient-state reduction constants for p-anisyl alcohol in the H502A (over 1800-fold) and H546A (over 35-fold) variants. Combination of QM/MM energy profiles, protonation predictors, molecular dynamics, mutagenesis and pH profiles provide a robust answer regarding the nature of the catalytic base. The histidine residue in front of the FAD ring, AAO His??2 (and CHO His???), acts as a base. For the two substrates assayed, it was shown that proton transfer preceded hydride transfer, although both processes are highly coupled. No stable intermediate was observed in the energy profiles, in contrast with that observed for CHO. QM/MM, together with solvent KIE (kinetic isotope effect) results, suggest a non-synchronous concerted mechanism for alcohol oxidation by AAO.     

Theoretical study of the novel sandwich compound [Au<Subscript>3</Subscript>Cl<Subscript>3</Subscript>Tr<Subscript>2</Subscript>]<Superscript>2+</Superscript>

A theoretical study of a sandwich compound with a metal monolayer sheet between two aromatic ligands is presented. A full geometry optimization of the [Au₃Cl₃Tr₂]²⁺ (1) compound, which is a triangular gold(I) monolayer sheet capped by chlorines and bounded to two cycloheptatrienyl (Tr) ligands was carried out using perturbation theory at the MP2 computational level and DFT. Compound (1) is in agreement with the 18–electron rule, the bonding nature in the complex may be interpreted from the donation interaction coming from the Tr rings to the Au array, and from the back-donation from the latter to the former. NICS calculations show a strong aromatic character in the gold monolayer sheet and Tr ligands; calculations done with HOMA, also report the same aromatic behavior on the cycloheptatrienyl fragments giving us an insight on the stability of (1). The Au –Au bond lengths indicate that an intramolecular aurophilic interaction among the Au(I) cations plays an important role in the bonding of the central metal sheet. Figure (a) Ground state geometry of complex 1; (b) Top view of compound 1 and Wiberg bond orders computed with the MP2/B1 computational method; (c) Lateral view of compound 1 and NICS values calculated with the MP2/B1 method; the values in parenthesis were obtained at the VWN/TZP level